Method for controlling ecological flow of river channel by coupling riverbed erosion and deposition evolution and wetted perimeter response
Patent Information
- Application Number
- CN202611041449.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2026-07-14
- Publication Date
- 2026-10-27
AI Technical Summary
[0005]本申请的主要目的在于提供一种耦合河床冲淤演变与湿周响应的河道生态流量控制方法,旨在解决对于多沙河流及冲淤活跃河段,如何提高对生态流量控制的精确度与可靠性的技术问题
获取目标河流的基础数据集,基于基础数据集构建水动力泥沙耦合模型,并通过水动力泥沙耦合模型分别模拟目标河流在多个工况下的河床冲淤过程,得到各工况下的河床冲淤厚度;基于河床冲淤厚度,更新目标河流的河道断面几何形态,得到更新后的河道断面几何形态;基于更新后的河道断面几何形态,重构湿周与流量之间的动态响应关系,并基于动态响应关系,确定目标河流的生态流量阈值;根据生态流量阈值生成并输出控制指令,控制指令用于控制目标河流上游水利工程泄放流量。获取目标河流基础数据并构建水动力泥沙耦合模型以模拟多工况下河床冲淤厚度,解决了传统方法默认断面形态不变导致无法反映真实河床演变的技术问题,实现了对冲淤演变过程的精细化量化表征,确保了后续计算基础的真实性。采用了基于冲淤厚度动态更新河道断面几何形态的技术手段,解决了静态断面与动态河流实际形态脱节的技术问题,实现了对河道过水能力的实时校准,消除了因断面变形带来的计算偏差。采用了基于更新后断面重构湿周与流量动态响应关系并确定生态流量阈值的技术手段,解决了传统湿周法在冲淤河段阈值识别失真及不稳定的技术问题,实现了生态流量阈值的自适应动态修正,显著提升了计算结果的精确度与可靠性。根据生态流量阈值直接生成并输出控制指令以调控上游水利工程泄放流量,实现了从理论计算到工程执行的闭环控制,满足了多沙河流生态流量精细化管理的实际应用需求。
Smart Images

Figure CN122987731A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of water conservancy and environmental engineering technology, and in particular to a method for controlling river ecological flow by coupling riverbed scouring and deposition evolution with wetted perimeter response. Background Technology
[0002] Ecological flow is the hydrological basis for maintaining the ecological structure and function of rivers. It is directly related to the effectiveness of river ecological protection, water resource allocation and river ecological restoration. Especially for rivers with high sediment content and active scouring and deposition, it is crucial whether the ecological base flow can truly reflect the dynamic changes in river hydraulic conditions.
[0003] Currently, traditional methods for controlling river ecological flow are based on calculating ecological flow thresholds at fixed river cross-sections. However, for river sections with significant scour and deposition evolution, the riverbed topography continuously changes with the process of water and sediment transport. Traditional methods do not consider the dynamic impact of scour and deposition evolution on cross-sectional morphology and wetted perimeter response, which can easily lead to problems such as distorted ecological flow threshold identification, insufficient calculation stability, and poor adaptability to rivers with high sediment content, thus affecting the accuracy and reliability of ecological flow control.
[0004] Therefore, improving the accuracy and reliability of ecological flow control is an urgent problem to be solved for rivers with high sediment content and active erosion and deposition. Summary of the Invention
[0005] The main objective of this application is to provide a method for controlling ecological flow in rivers that couples riverbed scour and deposition evolution with wetted period response, aiming to solve the technical problem of how to improve the accuracy and reliability of ecological flow control for sediment-laden rivers and river sections with active scour and deposition.
[0006] To achieve the above objectives, this application proposes a method for controlling river ecological flow that couples riverbed scouring and deposition evolution with wetted perimeter response. The method includes: Obtain the basic dataset of the target river, construct a hydrodynamic-sediment coupling model based on the basic dataset, and simulate the riverbed scouring and deposition process of the target river under multiple working conditions using the hydrodynamic-sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition. Based on the scouring and sedimentation thickness of the riverbed, the cross-sectional geometry of the target river is updated to obtain the updated cross-sectional geometry. Based on the updated river cross-sectional geometry, the dynamic response relationship between wetted perimeter and flow is reconstructed, and based on the dynamic response relationship, the ecological flow threshold of the target river is determined. Control commands are generated and output based on the ecological flow threshold, and the control commands are used to control the discharge flow of the upstream water conservancy project of the target river.
[0007] In one embodiment, the steps of obtaining a basic dataset of the target river, constructing a hydrodynamic-sediment coupling model based on the basic dataset, and simulating the riverbed scouring and deposition process of the target river under multiple working conditions using the hydrodynamic-sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition include: Obtain the basic dataset of the target river, and vectorize the target river based on the basic dataset to create a river channel geometry file. The basic dataset includes surveying data, hydrological data, and meteorological data. Hydrodynamic control equations are constructed based on the conservation of mass and momentum, and sediment transport control equations are constructed based on the sediment continuity equation. These sediment transport control equations are used to calculate riverbed elevation changes through sediment transport imbalance. The initial and boundary conditions of the hydrodynamic-sediment coupling model are determined based on the aforementioned river channel geometry file. Based on the initial and boundary conditions, the hydrodynamic control equations, and the sediment transport control equations, a hydrodynamic-sediment coupling model is determined. The riverbed roughness, sediment porosity, and hydraulic gradient are input into the hydrodynamic sediment coupling model to simulate the riverbed scouring and deposition process of the target river under various working conditions. The changes in riverbed elevation at each elevation point on the characteristic cross section of the target river under each working condition are obtained, and the changes in riverbed elevation are used as the riverbed scouring and deposition thickness. The working conditions include the high water period, the normal water period, and the low water period.
[0008] In one embodiment, the step of updating the cross-sectional geometry of the target river based on the riverbed erosion and deposition thickness to obtain the updated cross-sectional geometry includes: Extract the riverbed elevation change at each elevation point on the characteristic cross-section of the target river; The changes in riverbed elevation at each elevation point are superimposed onto the original elevation of the corresponding elevation point to obtain the updated elevation value of each elevation point; Based on the updated elevation values of all elevation points, the updated river channel cross-sectional geometry is generated.
[0009] In one embodiment, after the step of obtaining the riverbed elevation change at each elevation point on the characteristic cross-section of the target river under various working conditions and using the riverbed elevation change as the riverbed scouring and deposition thickness, the method further includes: Based on the changes in riverbed elevation at each elevation point on the characteristic cross-section and the corresponding water depth at each elevation point, the riverbed evolution rate is calculated. The riverbed evolution rate is used to indicate the degree of drastic change in the geometric shape of the river channel cross-section. Determine whether the riverbed evolution rate during the current calculation period is greater than a preset threshold; If the riverbed evolution rate during the current calculation period is greater than a preset threshold, then the step of reconstructing the dynamic response relationship between wetted perimeter and flow rate based on the updated river cross-sectional geometry is executed. If the riverbed evolution rate in the current calculation period is less than or equal to a preset threshold, then the dynamic response relationship or the ecological flow threshold of the previous calculation period will be used.
[0010] In one embodiment, the step of reconstructing the dynamic response relationship between wetted perimeter and flow based on the updated river cross-sectional geometry, and determining the ecological flow threshold of the target river based on the dynamic response relationship, includes: Based on the updated river cross-sectional geometry, the wetted perimeter value corresponding to each flow value under the preset flow sequence is calculated using the derivation relationship between the Manning formula and the Chezy formula, and a wetted perimeter sequence is generated. The wetted period sequence and the preset flow rate sequence are fitted using a continuously differentiable function to construct a fitting function that characterizes the correspondence between wetted period and flow rate, and the dynamic response relationship is determined based on the fitting function. Based on the dynamic response relationship, the curvature of the fitting function at each flow value is calculated according to the curvature calculation formula. The ecological flow threshold is determined based on the change in curvature, where the flow value corresponding to the peak curvature is the ecological flow threshold of the target river.
[0011] In one embodiment, the control command includes an initial control command, and the step of generating and outputting the control command based on the ecological flow threshold includes: The ecological flow threshold is compared with the actual operational constraints of the water conservancy project to generate an initial control command, which is used to adjust the opening of the water conservancy project's discharge facilities; wherein, the actual operational constraints include at least one of the following: generator unit operating status, flood control limit water level during the flood season, and flood season replenishment and dry season scheduling rules; The initial control command is sent to the actuator of the water conservancy project, which includes at least one of a gate controller, a valve controller, or a unit speed governor.
[0012] In one embodiment, the control command further includes a target control command, and after the step of sending the initial control command to the actuator of the water conservancy project, the method further includes: After the water conservancy project executes the initial control command, it acquires the actual discharge flow monitoring data; Based on the actual discharge flow monitoring data, the simulation calculation of the hydrodynamic sediment coupling model is re-executed to obtain the updated riverbed scouring and deposition thickness. Based on the updated riverbed scour and sedimentation thickness, the process returns to the step of updating the cross-sectional geometry of the target river based on the updated riverbed scour and sedimentation thickness, until the difference between the ecological flow thresholds calculated in two adjacent iterations is less than the preset convergence threshold. Then, the iteration stops and the current control command is output as the target control command.
[0013] Furthermore, to achieve the above objectives, this application also proposes a river ecological flow control device that couples riverbed scour and deposition evolution with wetted perimeter response. The river ecological flow control device that couples riverbed scour and deposition evolution with wetted perimeter response includes: The data acquisition module is used to acquire the basic dataset of the target river, construct a hydrodynamic-sediment coupling model based on the basic dataset, and simulate the riverbed scouring and deposition process of the target river under multiple working conditions through the hydrodynamic-sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition. The morphology update module is used to update the cross-sectional geometry of the target river based on the riverbed erosion and deposition thickness, so as to obtain the updated cross-sectional geometry. The flow determination module is used to reconstruct the dynamic response relationship between wetted perimeter and flow based on the updated river cross-sectional geometry, and to determine the ecological flow threshold of the target river based on the dynamic response relationship. The instruction output module is used to generate and output control instructions based on the ecological flow threshold, and the control instructions are used to control the discharge flow of the upstream water conservancy project of the target river.
[0014] In addition, to achieve the above objectives, this application also proposes an electronic device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, the computer program being configured to implement the steps of the river ecological flow control method coupled with riverbed scouring and deposition evolution and wetted perimeter response as described above.
[0015] In addition, to achieve the above objectives, this application also proposes a non-transitory computer storage medium storing a computer program, which, when executed by a processor, implements the steps of the river ecological flow control method that couples riverbed scouring and deposition evolution with wetted period response as described above.
[0016] One or more technical solutions proposed in this application have at least the following technical effects: This study acquires a basic dataset of the target river, constructs a hydrodynamic-sediment coupling model based on this dataset, and simulates the riverbed scouring and deposition processes under multiple operating conditions using this model, obtaining the scouring and deposition thickness under each condition. Based on the riverbed scouring and deposition thickness, the cross-sectional geometry of the target river is updated, resulting in the updated cross-sectional geometry. Based on the updated cross-sectional geometry, the dynamic response relationship between wetted perimeter and flow is reconstructed, and based on this dynamic response relationship, the ecological flow threshold of the target river is determined. Control commands are generated and output based on the ecological flow threshold, and these commands are used to control the discharge flow of upstream water conservancy projects on the target river. Acquiring basic data of the target river and constructing a hydrodynamic-sediment coupling model to simulate riverbed scouring and deposition thickness under multiple operating conditions solves the technical problem of traditional methods that default to unchanged cross-sectional morphology, thus failing to reflect the actual riverbed evolution. This achieves a refined and quantitative characterization of the scouring and deposition evolution process, ensuring the authenticity of the subsequent calculation basis. A technique based on dynamically updating the river cross-sectional geometry based on scour and sediment thickness was adopted, solving the technical problem of the disconnect between static cross-sections and the actual dynamic river morphology. This enabled real-time calibration of the river's water-carrying capacity and eliminated calculation errors caused by cross-sectional deformation. A technique based on reconstructing the dynamic response relationship between wetted perimeter and flow based on the updated cross-section and determining the ecological flow threshold was also employed. This solved the technical problems of threshold identification distortion and instability in scour and sediment-affected river sections using the traditional wetted perimeter method, achieving adaptive dynamic correction of the ecological flow threshold and significantly improving the accuracy and reliability of the calculation results. Control commands were directly generated and output based on the ecological flow threshold to regulate the discharge flow of upstream water conservancy projects, realizing closed-loop control from theoretical calculation to engineering execution and meeting the practical application needs of refined ecological flow management in sediment-laden rivers. Attached Figure Description
[0017] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application.
[0018] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, for those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0019] Figure 1 This is a flowchart illustrating the first embodiment of the river ecological flow control method that couples riverbed scouring and deposition evolution with wetted perimeter response in this application. Figure 2 This is a schematic diagram showing the morphological changes of section 5 under different hydrological conditions in this application. Figure 3This is a schematic diagram of the calculation process for the river ecological flow control method that couples riverbed scouring and deposition evolution with wetted perimeter response, as described in this application. Figure 4 This is a schematic diagram of the module structure of the river ecological flow control device that couples riverbed scouring and deposition evolution with wetted period response, as described in this application.
[0020] The purpose, features, and advantages of this application will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation
[0021] The technical solutions in the embodiments of this application will be clearly and completely described below with reference to the embodiments of this application. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of this application.
[0022] This application provides a method for controlling river ecological flow that couples riverbed scouring and deposition evolution with wetted perimeter response. Specifically, refer to... Figure 1 , Figure 1 This is a flowchart illustrating the first embodiment of the river ecological flow control method coupling riverbed scour and deposition evolution with wetted perimeter response, as described in this application. In this embodiment, the river ecological flow control method coupling riverbed scour and deposition evolution with wetted perimeter response includes steps S10-S40:
[0023] Step S10: Obtain the basic dataset of the target river, construct a hydrodynamic-sediment coupling model based on the basic dataset, and simulate the riverbed scouring and deposition process of the target river under multiple working conditions through the hydrodynamic-sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition.
[0024] It should be noted that the target river can be understood as the river studied in this application's method for river ecological flow control coupled with riverbed scour and deposition evolution and wetted perimeter response, specifically referring to river sections with high sediment content or significant riverbed scour and deposition evolution characteristics. The basic dataset is the set of original observational data used to support model construction, including topographic and geomorphological data, hydrological survey data, sediment transport data, and boundary condition data of the target river. The hydrodynamic-sediment coupling model is a computational model that simulates the interaction between water flow and sediment transport through numerical calculations. In this model, hydrodynamic conditions drive sediment transport, and the scour and deposition changes in sediment, in turn, alter the riverbed boundary and affect the flow structure. The operating condition refers to the river's operational state under specific time periods or specific water and sediment inflow conditions, corresponding to various hydrological situations, such as the high-water season, normal-water season, low-water season, or different typical years. The riverbed scour and deposition process refers to the dynamic changes in the riverbed caused by water flow, either scouring or deposition. Riverbed scouring and deposition thickness is an indicator used to quantify the riverbed scouring and deposition process. It can be understood as the change in the riverbed surface elevation relative to its initial state per unit time or under specific working conditions.
[0025] Step S20: Based on the thickness of riverbed scouring and sedimentation, update the cross-sectional geometry of the target river to obtain the updated cross-sectional geometry.
[0026] It should be noted that the cross-sectional geometry of a river channel refers to the shape of the cross section perpendicular to the direction of river flow. It can be composed of a series of elevation point coordinates and describes geometric characteristics such as river width, water depth, flow area, and wetted perimeter. The updated cross-sectional geometry refers to the cross-sectional shape of the current riverbed under the actual state, which is obtained by recalculating after considering the scouring and deposition thickness under specific working conditions.
[0027] Step S30: Based on the updated river cross-sectional geometry, reconstruct the dynamic response relationship between wetted perimeter and flow, and determine the ecological flow threshold of the target river based on the dynamic response relationship.
[0028] It should be noted that wetted perimeter refers to the circumference of the water body in contact with the riverbed at the cross-section of the river. The dynamic response relationship can be understood as follows: as the riverbed cross-sectional shape changes due to scouring and deposition, the relationship between wetted perimeter and flow rate is no longer a fixed correspondence, but rather a mapping relationship that dynamically adjusts with time or operating conditions. The ecological flow threshold refers to the minimum or suitable flow rate value required to meet the basic needs of the river ecosystem.
[0029] Step S40: Generate and output control commands based on the ecological flow threshold.
[0030] It should be noted that the control command is used to control the discharge flow of upstream water conservancy projects on the target river. For example, upstream water conservancy projects refer to hydraulic structures located upstream of the target river that have flow regulation functions, such as reservoirs, hydropower stations, and dams. Discharge flow refers to the amount of water released by the water conservancy project into the downstream river channel through the opening and closing of gates or other control facilities.
[0031] In this embodiment, basic data of the target river is acquired and a hydrodynamic-sediment coupling model is constructed to simulate the riverbed scour and deposition thickness under multiple working conditions. This solves the technical problem of traditional methods, which default to unchanged cross-sectional morphology and cannot reflect the actual riverbed evolution. It achieves a refined quantitative characterization of the scour and deposition evolution process, ensuring the authenticity of the subsequent calculation basis. A technique for dynamically updating the river cross-sectional geometry based on scour and deposition thickness is adopted, solving the technical problem of the disconnect between static cross-sections and the actual dynamic river morphology. This enables real-time calibration of the river's water-carrying capacity and eliminates calculation deviations caused by cross-sectional deformation. A technique for reconstructing the dynamic response relationship between wetted perimeter and flow based on the updated cross-section and determining the ecological flow threshold is adopted, solving the technical problems of threshold identification distortion and instability in scour and deposition river sections in the traditional wetted perimeter method. This achieves adaptive dynamic correction of the ecological flow threshold, significantly improving the accuracy and reliability of the calculation results. Control commands are directly generated and output based on the ecological flow threshold to regulate the discharge flow of upstream water conservancy projects, realizing closed-loop control from theoretical calculation to engineering execution, and meeting the practical application needs of refined ecological flow management in sediment-rich rivers.
[0032] In one implementation, step S10 includes: Step S101: Obtain the basic dataset of the target river, and vectorize the target river based on the basic dataset to create a river channel geometry file.
[0033] For example, the basic dataset for a target river may include surveying data, hydrological data, and meteorological data. Specifically, the surveying data includes measured cross-sectional morphology, cross-sectional spacing, and river vector data; the hydrological data includes upstream and downstream boundary channel water levels, flow rates, flow velocities, and suspended solids concentrations; and the meteorological data includes rainfall, evaporation, wind direction, and wind speed in the area where the target river is located. Non-digital geographic information of the target river, such as remote sensing imagery, can be converted into computer-readable vector data, such as .shp format, using GIS tools for accurate depiction of the river channel's spatial morphology. The river channel geometry file can be understood as a file storing vectorized geometric information of the target river channel, such as its boundaries, cross-sectional locations, and cross-sectional spacing, serving as the input basis for the hydrodynamic-sediment coupling model.
[0034] For example, a representative monitoring section #5 of a typical alluvial river in the middle reaches of the Yangtze River can be selected as the research object. The total length of the river section is 38.6 km. The river channel is frequently alternating between scouring and silting. During the flood season, the high sediment load of the water flow causes the riverbed to be silted up and lifted. During the dry season, the discharge of clear water causes the main channel to be scoured and cut down. The morphology of the river channel section and the wetted perimeter-discharge relationship have significant time-varying characteristics.
[0035] Step S102: Based on the conservation of mass and momentum, construct the hydrodynamic control equations and the sediment transport control equations based on the sediment continuity equation. The sediment transport control equations are used to calculate the changes in riverbed elevation through sediment transport imbalance.
[0036] It should be noted that the hydrodynamic governing equations are mathematical equations derived based on the conservation of mass and momentum, used to simulate the changes in hydrodynamic elements such as velocity, water level, and flow rate within a river channel. The sediment continuity equation is a mathematical equation describing the conservation of total sediment load during sediment transport, used to construct the sediment transport governing equations. These equations are derived from the sediment continuity equation and are used to calculate the amount of sediment transported under the influence of water flow and changes in riverbed elevation.
[0037] For example, the hydrodynamic governing equations can be constructed based on the one-dimensional unsteady flow hydrodynamic equations of mass conservation and momentum conservation, and the specific formula is expressed as follows: in, Q For traffic; Z Water level; A The cross-sectional area of the main channel for water flow; x t represents the spatial coordinates along the direction of water flow; t represents the time coordinates along the direction of water flow. Bw To adjust the width; B Mainly the width of the cross-section; R The hydraulic radius; n For riverbed roughness; u The average velocity across the cross section; g It is the acceleration due to gravity; q This refers to the side inflow flow rate.
[0038] The sediment transport control equation uses the sediment continuity equation to characterize the riverbed erosion and deposition process, and calculates the riverbed elevation change through sediment transport imbalance. The specific formula is as follows: When considering source and sink terms, such as tributary inflows or artificial sand mining, it can be expressed as: in, Riverbed elevation; For time; This represents the porosity of riverbed sediments, typically ranging from 0.3 to 0.4. The river is wide; Volumetric sediment transport rate ( It is often processed as mass transport rate or volume transport rate, depending on the density of the input sediment; This represents the distance along the river channel. This refers to the volume of sediment per unit length of river channel.
[0039] Step S103: Determine the initial and boundary conditions of the hydrodynamic-sediment coupling model based on the river channel geometry file.
[0040] It should be noted that initial conditions refer to the initial state parameters when the model starts, such as the initial riverbed elevation, initial water level, and initial flow rate. Boundary conditions are the boundary constraint parameters of the model's computational domain, such as upstream inflow, downstream water level, and lateral inflow, which limit the model's operating range.
[0041] Step S104: Based on the initial and boundary conditions, the hydrodynamic control equations, and the sediment transport control equations, determine the hydrodynamic-sediment coupling model.
[0042] It should be noted that the hydrodynamic-sediment coupling model can be understood as a numerical model that solves the hydrodynamic control equations and the sediment transport control equations simultaneously, realizing bidirectional feedback calculation of water flow and sediment transport.
[0043] Step S105: Input the riverbed roughness, sediment porosity and hydraulic gradient into the hydrodynamic sediment coupling model to simulate the riverbed scouring and deposition process of the target river under various working conditions. Obtain the riverbed elevation change at each elevation point on the characteristic section of the target river under each working condition. Use the riverbed elevation change as the riverbed scouring and deposition thickness. The working conditions include the high water period, the normal water period and the low water period.
[0044] It should be noted that riverbed roughness is a parameter reflecting the impact of riverbed surface roughness on water flow resistance; a larger value indicates a rougher riverbed and greater water flow resistance. Sediment porosity refers to the proportion of pore volume to total volume in riverbed sediments, typically ranging from 0.3 to 0.4, but here a value of 0.35 can be used to calculate the effective volume for sediment transport. Hydraulic gradient is the water level difference per unit river length, reflecting the dynamic conditions of water flow and affecting sediment transport capacity. Riverbed scouring and deposition processes represent the dynamic changes in the riverbed caused by the interaction of water flow and sediment. A characteristic cross-section can be understood as a representative monitoring section in the target river, used for focused analysis of riverbed morphological changes. Elevation points refer to discrete sampling points on the characteristic cross-section used to describe the riverbed morphology; each point corresponds to a specific horizontal position and vertical elevation. Riverbed elevation change is the difference between an elevation point on the characteristic cross-section and its initial elevation during the calculation period, directly reflecting the magnitude of scouring and deposition at that point. Riverbed scour and deposition thickness can be understood as the total amount of scour and deposition in the riverbed obtained by summing up the changes in riverbed elevation at various elevation points, and is used to quantify the degree of riverbed evolution.
[0045] For example, after inputting core parameters such as riverbed roughness, sediment porosity, and hydraulic gradient into the hydrodynamic-sediment coupling model, the model is calibrated and verified. The absolute error of the water level can be set to no more than 10 cm, and the relative errors of flow rate and sediment deposition to no more than 10%. Furthermore, using the calibrated and verified hydrodynamic-sediment coupling model, the riverbed sediment deposition process is simulated under different hydrological conditions and sediment concentrations during the flood season, normal water season, and dry season. The calculation results show that during the high-water season, the sedimentation amplitude of the main channel at section 5 is 0.60–1.10 m, and the sedimentation amplitude of the side banks is 0.30–0.50 m; during the normal water season, the sedimentation amplitude at each measuring point shows no uniform pattern, with significant local differences; during the dry season, the maximum downcut of the main channel is 0.60–0.80 m, and the side banks only experience slight scouring of 0.10–0.20 m. The cross-sectional morphology is significantly distorted in the three periods, lacking parallel rise and fall characteristics, requiring dynamic updates of the cross-section for subsequent calculations.
[0046] In this implementation, a technique was employed to acquire and vectorize the target river's basic dataset to create a river channel geometry file. This addresses the lack of precise spatial morphological input in traditional models, enabling accurate digital representation of river channel geometric features and providing a reliable foundation for subsequent model construction. Hydrodynamic control equations were constructed based on mass and momentum conservation, and sediment transport control equations were built based on the sediment continuity equation. This solved the problem of traditional models not coupling water-sediment interactions, achieving bidirectional feedback simulation of water flow and sediment transport, and improving the simulation accuracy of riverbed scouring and deposition processes. The technique of determining initial and boundary conditions based on the river channel geometry file standardized the model's operating conditions, ensuring consistency in simulation results. By integrating initial conditions, boundary conditions, and control equations to determine the hydrodynamic-sediment coupling model, the technique solved the problem of traditional models being unable to simultaneously simulate water-sediment processes, enabling dynamic quantitative simulation of riverbed scouring and deposition processes and providing data support for subsequent cross-section updates. By employing the technical means of inputting riverbed roughness, sediment porosity, and hydraulic gradient, and simulating the riverbed scouring and deposition process under multiple working conditions, the evolution characteristics of the riverbed in different hydrological periods were fully captured, ensuring the scenario adaptability of subsequent ecological flow calculations.
[0047] In one implementation, step S20 includes: Step S201: Extract the riverbed elevation change at each elevation point on the characteristic cross section of the target river.
[0048] It should be noted that the elevation points are discrete sampling points distributed laterally along the characteristic cross-section, each point corresponding to a specific horizontal position and vertical height. The change in riverbed elevation is the difference between a certain elevation point on the characteristic cross-section and the initial reference elevation during the calculation period. Positive values indicate riverbed siltation and uplift, while negative values indicate riverbed scouring and downcutting.
[0049] Step S202: The changes in riverbed elevation at each elevation point are superimposed onto the original elevation of the corresponding elevation point to obtain the updated elevation value of each elevation point.
[0050] It should be noted that the original elevation can be understood as the baseline height value of each elevation point on the characteristic cross-section in its initial state, i.e., before the scouring and deposition simulation. The updated elevation value is the result of algebraically adding the change in riverbed elevation to the original elevation of the corresponding elevation point. The updated elevation value reflects the actual elevation state of the riverbed at the current time.
[0051] Step S203: Based on the updated elevation values of all elevation points, generate the updated river channel cross-sectional geometry.
[0052] For example, the updated river channel cross-sectional geometry can be understood as a new cross-sectional shape formed by connecting the updated elevation values of all elevation points, including key geometric parameters such as river width, flow area, and wetted perimeter, reflecting the actual morphological characteristics of the current riverbed. The riverbed elevation change ΔZi at each elevation point in section #5 can be extracted, and the time-series simulation results of the riverbed scour and deposition spatial distribution are shown below. Figure 2 As shown.
[0053] In this embodiment, by extracting the changes in riverbed elevation at each elevation point on the feature cross-section, the technical problem of traditional methods being unable to quantify local riverbed morphological changes is solved, achieving precise capture of the spatial distribution of riverbed scouring and deposition, and avoiding morphological distortion caused by overall averaging. By employing the technique of superimposing the changes in riverbed elevation onto the corresponding original elevation to obtain updated elevation values, the technical problem of traditional static cross-sections being unable to reflect dynamic evolution is solved, achieving precise correction of the height of individual elevation points and ensuring the precision of cross-section morphological updates. Furthermore, by employing the technique of generating new river channel cross-section geometry based on all updated elevation values, the technical problem of traditional cross-section morphology lagging behind actual riverbed changes is solved.
[0054] In one implementation, after step S10, the method further includes: calculating the riverbed evolution rate based on the riverbed elevation change at each elevation point on the characteristic cross-section and the water depth corresponding to each elevation point; the riverbed evolution rate is used to indicate the degree of drastic change in the geometric shape of the river cross-section; determining whether the riverbed evolution rate in the current calculation period is greater than a preset threshold; if the riverbed evolution rate in the current calculation period is greater than the preset threshold, then step S30 is executed; if the riverbed evolution rate in the current calculation period is less than or equal to the preset threshold, then the dynamic response relationship or ecological flow threshold of the previous calculation period is used.
[0055] It should be noted that the riverbed evolution rate is used to quantitatively characterize the degree of drastic change in the river channel cross-section morphology. It can be obtained by calculating the relative ratio of the change in riverbed elevation at each elevation point of the cross-section to the corresponding water depth. A larger value indicates a more drastic change in riverbed morphology, while a smaller value indicates a more stable morphology. Water depth refers to the water depth at a specific elevation point in the river channel cross-section. The preset threshold can be a critical value for the riverbed evolution rate pre-set based on engineering experience or management precision requirements. When the calculated value exceeds this limit, it is considered that the cross-section morphology has changed significantly, and the ecological flow needs to be recalculated; otherwise, the original result can be maintained. The current calculation period refers to the time step or specific time period of the hydrodynamic-sediment coupling model simulation, such as ten days, months, or seasons. The previous calculation period refers to the previous adjacent time unit of the current calculation moment. When the riverbed evolution rate does not exceed the preset threshold, the dynamic response relationship or ecological flow threshold calculated in the previous period is directly adopted, and the cross-section reconstruction and threshold identification calculation process is not repeated.
[0056] For example, the riverbed evolution rate is calculated under different hydrological conditions based on the formula for calculating the riverbed evolution rate: in, The degree of riverbed change at the i-th point of the characteristic cross-section. Let |R(t)| be the water depth corresponding to that point, and n be the total number of elevation points. A larger |R(t)| indicates a more drastic evolution, and vice versa. In this embodiment, |R(t)| is set to... A value greater than 0.01 indicates significant cross-sectional evolution, necessitating reconstruction of the wetted perimeter-discharge relationship. Calculations show that the riverbed evolution rate |R(t)| during the wet, normal, and dry seasons are 0.0294, 0.0072, and 0.0126, respectively. The wet and dry seasons exhibit more dramatic cross-sectional morphological evolution, requiring reconstruction of the wetted perimeter-discharge relationship and identification of ecological flow thresholds. The riverbed evolution rate during the normal season is moderate, allowing for selective reconstruction based on engineering accuracy requirements. If the evolution rate is low during stable periods, the reconstruction step can be skipped, and the ecological flow results from the previous stage can be used, saving model computation costs and improving operational efficiency.
[0057] In this embodiment, the technique of calculating the riverbed evolution rate based on the changes in riverbed elevation and corresponding water depth at various elevation points on the characteristic cross-section solves the technical problem of being unable to quantitatively judge the severity of cross-sectional morphological changes. This achieves accurate perception and quantitative assessment of dynamic changes in riverbed scouring and deposition, providing a scientific basis for subsequent calculation decisions. By employing a technique to determine whether the riverbed evolution rate in the current calculation period is greater than a preset threshold, necessary and unnecessary calculation scenarios are effectively distinguished. By executing step S30 to reconstruct the relationship and determine the threshold when the riverbed evolution rate is greater than the preset threshold, real-time correction of the ecological flow threshold is achieved during periods of drastic riverbed changes, ensuring the accuracy of ecological flow control. By employing a technique to reuse the results of the previous calculation period when the riverbed evolution rate is less than or equal to the preset threshold, the technical problem of low efficiency and slow response caused by redundant and complex calculations during stable cross-section periods is solved. This achieves optimized allocation of computing resources and a significant improvement in computational efficiency, while enhancing the stability and determinism of ecological flow management.
[0058] In one implementation, step S30 includes: Step S301: Based on the updated river cross-sectional geometry, calculate the wetted perimeter value corresponding to each flow value under the preset flow sequence using the derivation relationship between the Manning formula and the Chezy formula, and generate the wetted perimeter sequence.
[0059] For example, the cross-sectional flow and corresponding wetted perimeter sequence after scouring and sedimentation simulation are calculated based on the derivation relationship between Manning's formula and Chezy's formula. The corresponding formula is: in, QFor river flow, m 3 / s; A Main channel cross-sectional area, m 2 ; Let m be the wetted perimeter of the river cross-section; j Hydraulic gradient; n For example, the overall roughness of a river channel. n =0.03, hydraulic gradient j =0.0008, calculate the flow range Q Control it within 10~80m3 / s.
[0060] Step S302: Fit the wetted period sequence and the preset flow rate sequence according to the continuously differentiable function to construct a fitting function to characterize the correspondence between wetted period and flow rate, and determine the dynamic response relationship according to the fitting function.
[0061] For example, Table 1 can be used as a reference for the wet cycle-corresponding flow sequence of the current cross-section and the cross-sections during the wet, normal and dry seasons.
[0062] Table 1. Relationship between wetted perimeter and corresponding discharge sequence after scour and deposition evolution at section 5 under different hydrological conditions. Various methods can be used to construct fitting functions for the cross-sectional wetted cycle sequence and the corresponding flow rate sequence, such as linear functions, quadratic functions, logarithmic functions, and exponential functions. For example, here we use a cubic polynomial function to fit the current, wet and dry season cross-sectional wetted cycle-corresponding flow rate sequences. The fitting function expression and R value are shown below.
[0063] The formula corresponding to the current situation is: The formula corresponding to the high-water season is: The formula corresponding to the dry season is: Step S303: Based on the dynamic response relationship, calculate the curvature of the fitted function at each flow rate value according to the curvature calculation formula.
[0064] For example, the curvature of the fitting function P(Q) is calculated for the wet season and the dry season respectively to determine the ecological flow. The formula for solving the curvature is expressed as:
[0065] Step S304: Determine the ecological flow threshold based on the change in curvature, wherein the flow value corresponding to the peak curvature is the ecological flow threshold of the target river.
[0066] For example, specific methods for dynamically identifying ecological flow thresholds include the method of identifying regions with significant curvature changes or the curvature peak value. The solution method uses the sensitive flow range or curvature peak as the ecological flow threshold. The current ecological flow thresholds for the wet season and dry season are 31.65 m³ / s, 28.42 m³ / s, and 35.76 m³ / s, respectively.
[0067] In this embodiment, based on the updated river cross-sectional geometry, the technique of calculating the wetted perimeter value corresponding to each flow rate under a preset flow sequence and generating a wetted perimeter sequence using the derivation relationship between the Manning formula and the Chezy formula solves the technical problem of the traditional method relying on static cross-sections, which leads to a disconnect between the wetted perimeter calculation and the actual riverbed state. This achieves accurate quantification of the wetted perimeter parameter as it dynamically evolves with the riverbed, providing a reliable data foundation for subsequent relationship construction. The wetted perimeter sequence is fitted to the preset flow sequence using a continuously differentiable function to construct a fitting function representing the correspondence between wetted perimeter and flow rate, and to determine the dynamic response relationship. This solves the technical problem that discrete data points cannot reflect continuous change patterns, providing a continuously differentiable function carrier for curvature calculation. By adopting a technique based on the dynamic response relationship and calculating the curvature of the fitting function at each flow rate value according to the curvature calculation formula, the technique of the traditional method being unable to accurately locate the sensitive area of wetted perimeter change solves the technical problem. This achieves a quantitative characterization of the bending characteristics of the wetted perimeter during flow rate changes, providing a mathematical basis for identifying critical turning points. By employing a technique that determines the ecological flow threshold based on curvature changes and uses the flow value corresponding to the peak curvature as the target river's ecological flow threshold, the technical problem of threshold identification distortion caused by cross-sectional distortion in traditional inflection point methods in sediment-laden rivers has been solved. This achieves adaptive dynamic correction of the ecological flow threshold, significantly improving the accuracy and reliability of the calculation results in sediment-laden rivers and actively eroded and deposited river sections.
[0068] In one embodiment, the control command is used to adjust the opening degree of the spillway of the water conservancy project, and the control command includes an initial control command. Based on this, step S40 includes: comparing the ecological flow threshold with the actual operating constraints of the water conservancy project to generate an initial control command; wherein, the actual operating constraints include at least one of the generator unit operating status, the flood control limit water level during the flood season, and the flood season replenishment and dry season scheduling rules; and sending the initial control command to the actuator of the water conservancy project, the actuator including at least one of the gate controller, valve controller, or generator speed governor.
[0069] Among them, the generator unit operating status refers to the different operating conditions of the hydropower station units, such as power generation, standby, maintenance, or shutdown, which directly affects the allocation of flow available for ecological release. The flood control limit water level during the flood season refers to the highest permissible water level strictly controlled during the flood season to reserve flood control storage capacity. The flood season replenishment and dry season scheduling rules can be understood as a scheduling strategy that utilizes surplus water during the flood season to replenish or store water downstream to maintain ecological flow during the dry season.
[0070] It should be noted that after step S40, the process also includes: after the water conservancy project executes the initial control command, acquiring actual discharge flow monitoring data; based on the actual discharge flow monitoring data, re-executing the simulation calculation of the hydrodynamic sediment coupling model to obtain the updated riverbed scour and deposition thickness; based on the updated riverbed scour and deposition thickness, returning to step S20, until the difference between the ecological flow thresholds obtained from two adjacent iterations is less than the preset convergence threshold, stopping the iteration and outputting the current control command as the target control command.
[0071] The actual discharge flow monitoring data refers to the real-time data of the actual discharge volume collected by sensors such as flow meters and water level gauges installed downstream of the water conservancy project. After inputting the actual discharge flow as a new boundary condition into the sediment model, the riverbed scouring and deposition changes are re-simulated, resulting in the updated riverbed scouring and deposition thickness. A preset convergence threshold is used to determine the allowable error range for determining whether the calculation results are stable. When the difference between two adjacent ecological flow thresholds is less than the preset convergence threshold, the calculation is considered to have reached a stable state. For example, 1% to 3% of the current ecological flow threshold can be used as the preset convergence threshold. The calculation flow diagram of the above embodiments and implementation methods of this application is shown below. Figure 3 As shown.
[0072] In this embodiment, the technical means of comparing the ecological flow threshold with the actual operational constraints of the water conservancy project to generate initial control commands solves the technical problem of the disconnect between theoretical ecological flow and actual engineering operation conditions, achieving a balance between ecological needs and engineering safety, and ensuring the executability of control commands. By acquiring several monitoring segments of actual downstream flow after sending commands, real-time perception of the actual hydrological situation downstream is achieved, providing real data for model feedback. The technical means of re-executing the hydrodynamic-sediment coupling model simulation based on monitoring data to obtain the updated riverbed scour and deposition thickness solves the technical problem that a single calculation cannot reflect the long-term effects of real water flow, achieving bidirectional dynamic coupling between water and sediment processes and riverbed evolution. The technical means of returning to the updated cross-section based on the updated scour and deposition thickness and iteratively calculating until the difference between two adjacent thresholds is less than the preset convergence threshold solves the technical problems of result drift and instability in traditional static calculations in sediment-rich rivers, achieving adaptive convergence and dynamic optimization of the ecological flow threshold. The final output target control command has extremely high accuracy and environmental adaptability.
[0073] Based on the same inventive concept, such as Figure 4 As shown, this application also provides a river ecological flow control device that couples riverbed scour and deposition evolution with wetted perimeter response. The river ecological flow control device that couples riverbed scour and deposition evolution with wetted perimeter response includes: Data acquisition module 10 is used to acquire the basic dataset of the target river, construct a hydrodynamic sediment coupling model based on the basic dataset, and simulate the riverbed scouring and deposition process of the target river under multiple working conditions through the hydrodynamic sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition. The morphology update module 20 is used to update the cross-sectional geometry of the target river based on the riverbed scour and sedimentation thickness, so as to obtain the updated cross-sectional geometry. The flow determination module 30 is used to reconstruct the dynamic response relationship between wetted perimeter and flow based on the updated river cross-sectional geometry, and to determine the ecological flow threshold of the target river based on the dynamic response relationship. The instruction output module 40 is used to generate and output control instructions based on the ecological flow threshold, and the control instructions are used to control the discharge flow of the upstream water conservancy project of the target river.
[0074] It should be noted that the river ecological flow control device coupled with riverbed scour and sedimentation evolution and wetted period response provided in this application embodiment and the river ecological flow control method coupled with riverbed scour and sedimentation evolution and wetted period response provided in this application embodiment are based on the same application concept. Therefore, the specific implementation of this embodiment can refer to the implementation of the aforementioned river ecological flow control method coupled with riverbed scour and sedimentation evolution and wetted period response, and the repeated parts will not be described again.
[0075] In some embodiments, an electronic device provided in this application includes: a memory, a processor, and a computer program stored in the memory and executable on the processor. The computer program is configured to implement the above-described method for controlling river ecological flow that couples riverbed scouring and deposition evolution with wetted perimeter response.
[0076] Specifically, the processor may include, for example, a general-purpose microprocessor, an instruction set processor and / or an associated chipset and / or a special-purpose microprocessor (e.g., an application-specific integrated circuit (ASIC)), etc. The processor may also include onboard memory for caching purposes. The processor may be a single processing unit or multiple processing units for performing different actions of the method flow according to embodiments of this application.
[0077] Memory can be any medium capable of containing, storing, transmitting, propagating, or transmitting instructions. For example, memory can include, but is not limited to, electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices, instruments, or propagation media. Specific examples of memory include: magnetic storage devices such as magnetic tape or hard disk drives (HDDs); optical storage devices such as optical discs (CD-ROMs); and also random access memory (RAM) or flash memory; and / or wired / wireless communication links.
[0078] This application also provides a non-transitory computer storage medium storing a computer program. When executed by a processor, this program implements the aforementioned method for river ecological flow control that couples riverbed scouring and deposition evolution with wetted perimeter response. This storage medium may be included in the device / apparatus / system described in the above embodiments; or it may exist independently and not be assembled into that device / apparatus / system. The aforementioned non-transitory computer storage medium carries one or more programs, which, when executed, implement the method as described in the embodiments or implementations of this application.
[0079] According to embodiments of this application, a non-transitory computer storage medium can be a computer-readable signal medium or a computer-readable storage medium, or any combination thereof. A computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this application, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in connection with an instruction execution system, apparatus, or device. A computer-readable signal medium can also be any storage medium other than a computer-readable storage medium that can transmit, propagate, or transfer a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the storage medium can be transmitted using any suitable medium, including but not limited to: wireless, wired, optical fiber, radio frequency signals, etc., or any suitable combination thereof.
[0080] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.
[0081] Those skilled in the art will understand that the features described in the various embodiments of this application can be combined and / or combined in various ways, even if such combinations or combinations are not explicitly described in this application. In particular, the features described in the various embodiments of this application can be combined and / or combined in various ways without departing from the spirit and teachings of this application. All such combinations and / or combinations fall within the scope of this application. Therefore, the scope of this application should not be limited to the above embodiments. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.
Claims
1. A method for controlling river ecological flow by coupling riverbed scouring and deposition evolution with wetted perimeter response, characterized in that, The method includes: Obtain the basic dataset of the target river, construct a hydrodynamic-sediment coupling model based on the basic dataset, and simulate the riverbed scouring and deposition process of the target river under multiple working conditions using the hydrodynamic-sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition. Based on the scouring and sedimentation thickness of the riverbed, the cross-sectional geometry of the target river is updated to obtain the updated cross-sectional geometry. Based on the updated river cross-sectional geometry, the dynamic response relationship between wetted perimeter and flow is reconstructed, and based on the dynamic response relationship, the ecological flow threshold of the target river is determined. Control commands are generated and output based on the ecological flow threshold, and the control commands are used to control the discharge flow of the upstream water conservancy project of the target river.
2. The method as described in claim 1, characterized in that, The steps of obtaining the basic dataset of the target river, constructing a hydrodynamic-sediment coupling model based on the basic dataset, and simulating the riverbed scouring and deposition process of the target river under multiple working conditions using the hydrodynamic-sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition include: Obtain the basic dataset of the target river, and vectorize the target river based on the basic dataset to create a river channel geometry file. The basic dataset includes surveying data, hydrological data, and meteorological data. Hydrodynamic control equations are constructed based on the conservation of mass and momentum, and sediment transport control equations are constructed based on the sediment continuity equation. These sediment transport control equations are used to calculate riverbed elevation changes through sediment transport imbalance. The initial and boundary conditions of the hydrodynamic-sediment coupling model are determined based on the aforementioned river channel geometry file. Based on the initial and boundary conditions, the hydrodynamic control equations, and the sediment transport control equations, a hydrodynamic-sediment coupling model is determined. The riverbed roughness, sediment porosity, and hydraulic gradient are input into the hydrodynamic sediment coupling model to simulate the riverbed scouring and deposition process of the target river under various working conditions. The changes in riverbed elevation at each elevation point on the characteristic cross section of the target river under each working condition are obtained, and the changes in riverbed elevation are used as the riverbed scouring and deposition thickness. The working conditions include the high water period, the normal water period, and the low water period.
3. The method as described in claim 1, characterized in that, The step of updating the cross-sectional geometry of the target river based on the riverbed erosion and deposition thickness to obtain the updated cross-sectional geometry includes: Extract the riverbed elevation change at each elevation point on the characteristic cross-section of the target river; The changes in riverbed elevation at each elevation point are superimposed onto the original elevation of the corresponding elevation point to obtain the updated elevation value of each elevation point; Based on the updated elevation values of all elevation points, the updated river channel cross-sectional geometry is generated.
4. The method as described in claim 2, characterized in that, After the step of obtaining the riverbed elevation change at each elevation point on the characteristic cross section of the target river under various working conditions, and using the riverbed elevation change as the riverbed scouring and deposition thickness, the method further includes: Based on the changes in riverbed elevation at each elevation point on the characteristic cross-section and the corresponding water depth at each elevation point, the riverbed evolution rate is calculated. The riverbed evolution rate is used to indicate the degree of drastic change in the geometric shape of the river channel cross-section. Determine whether the riverbed evolution rate during the current calculation period is greater than a preset threshold; If the riverbed evolution rate during the current calculation period is greater than a preset threshold, then the step of reconstructing the dynamic response relationship between wetted perimeter and flow rate based on the updated river cross-sectional geometry is executed. If the riverbed evolution rate in the current calculation period is less than or equal to a preset threshold, then the dynamic response relationship or the ecological flow threshold of the previous calculation period will be used.
5. The method as described in claim 1, characterized in that, The step of reconstructing the dynamic response relationship between wetted perimeter and flow based on the updated river cross-sectional geometry, and determining the ecological flow threshold of the target river based on the dynamic response relationship, includes: Based on the updated river cross-sectional geometry, the wetted perimeter value corresponding to each flow value under the preset flow sequence is calculated using the derivation relationship between the Manning formula and the Chezy formula, and a wetted perimeter sequence is generated. The wetted period sequence and the preset flow rate sequence are fitted using a continuously differentiable function to construct a fitting function that characterizes the correspondence between wetted period and flow rate, and the dynamic response relationship is determined based on the fitting function. Based on the dynamic response relationship, the curvature of the fitting function at each flow value is calculated according to the curvature calculation formula. The ecological flow threshold is determined based on the change in curvature, where the flow value corresponding to the peak curvature is the ecological flow threshold of the target river.
6. The method as described in claim 1, characterized in that, The control commands include initial control commands, and the step of generating and outputting control commands based on the ecological flow threshold includes: The ecological flow threshold is compared with the actual operational constraints of the water conservancy project to generate an initial control command, which is used to adjust the opening of the water conservancy project's discharge facilities; wherein, the actual operational constraints include at least one of the following: generator unit operating status, flood control limit water level during the flood season, and flood season replenishment and dry season scheduling rules; The initial control command is sent to the actuator of the water conservancy project, which includes at least one of a gate controller, a valve controller, or a unit speed governor.
7. The method as described in claim 6, characterized in that, The control command further includes a target control command, and after the step of sending the initial control command to the actuator of the water conservancy project, it further includes: After the water conservancy project executes the initial control command, it acquires the actual discharge flow monitoring data; Based on the actual discharge flow monitoring data, the simulation calculation of the hydrodynamic sediment coupling model is re-executed to obtain the updated riverbed scouring and deposition thickness. Based on the updated riverbed scour and sedimentation thickness, the process returns to the step of updating the cross-sectional geometry of the target river based on the updated riverbed scour and sedimentation thickness, until the difference between the ecological flow thresholds calculated in two adjacent iterations is less than the preset convergence threshold. Then, the iteration stops and the current control command is output as the target control command.
8. A river ecological flow control device that couples riverbed scouring and deposition evolution with wetted period response, characterized in that, The river ecological flow control device that couples riverbed scouring and deposition evolution with wetted perimeter response includes: The data acquisition module is used to acquire the basic dataset of the target river, construct a hydrodynamic-sediment coupling model based on the basic dataset, and simulate the riverbed scouring and deposition process of the target river under multiple working conditions through the hydrodynamic-sediment coupling model to obtain the riverbed scouring and deposition thickness under each working condition. The morphology update module is used to update the cross-sectional geometry of the target river based on the riverbed scour and sedimentation thickness, so as to obtain the updated cross-sectional geometry. The flow determination module is used to reconstruct the dynamic response relationship between wetted perimeter and flow based on the updated river cross-sectional geometry, and to determine the ecological flow threshold of the target river based on the dynamic response relationship. The instruction output module is used to generate and output control instructions based on the ecological flow threshold, and the control instructions are used to control the discharge flow of the upstream water conservancy project of the target river.
9. An electronic device, characterized in that, The electronic device includes: a memory, a processor, and a computer program stored in the memory and executable on the processor, the computer program being configured to implement the steps of the river ecological flow control method coupled with riverbed scouring and deposition evolution and wetted perimeter response as described in any one of claims 1 to 7.
10. A non-transitory computer storage medium, characterized in that, The non-transitory computer storage medium stores a computer program, which, when executed by a processor, implements the steps of the river ecological flow control method as described in any one of claims 1 to 7, which couples riverbed scouring and deposition evolution with wetted perimeter response.