Yellow river downstream sediment transport water demand calculation method, system and equipment based on riverbed evolution equilibrium theory and medium
Through the calculation method based on the riverbed evolution equilibrium theory, the problem of calculating the sand transport water demand in the lower reaches of the Yellow River is solved, and the accurate estimate of the sand transport water demand is achieved, ensuring the flood transportation capacity of the river and the reasonable allocation of water resources.
Patent Information
- Application Number
- CN202510327007.2
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-03-19
- Publication Date
- 2025-06-13
AI Technical Summary
It is difficult to accurately calculate the amount of sand transport in the lower reaches of the Yellow River, especially under the influence of climate change and human activities, the amount of sand entering the yellow water decreases, the sand quality of the riverbed becomes coarse, and the sand transport efficiency decreases.
Based on the riverbed evolution equilibrium theory, a calculation method is provided, including obtaining river channel information and sediment parameters, calculating the comprehensive stability index of the riverbed, connecting the water flow continuity equation, the water flow motion equation and river relationship, calculating the equilibrium river parameters, and calculating the saturated sand transport capacity through the sediment transport equation, and then calculating the sand transport water demand under different incoming sand volumes.
This method can accurately reflect the dynamic response relationship between water and sand situation - riverbed evolution - sand transport water demand, ensure that the river channel in the lower reaches of the Yellow River maintains a certain flood transportation capacity, and is suitable for estimating sand transport water demand in the future period. It is of great significance to the rational allocation of water resources and river sediment management.
Smart Images
Figure FT_1 
Figure FT_2 
Figure FT_3
Abstract
Description
Technical Field
[0001] The invention belongs to the technical field of water conservancy engineering, and relates to a method, system, equipment and medium for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory. Background Art
[0002] The Yellow River is the river with the highest sediment content on earth. Due to "less water and more sediment, and the uncoordinated relationship between water and sediment", the lower reaches of the Yellow River have been seriously silted up and shrunk. The water demand for sediment transport is crucial to maintaining the scale of the main channel of the lower reaches of the Yellow River and ensuring the water and sediment transport channel. At present, many studies have been carried out on the water demand for sediment transport in the lower reaches of the Yellow River, but most of the results are based on the water and sediment conditions and riverbed boundary conditions before 2000. Since the beginning of the 21st century, on the one hand, due to climate change and strong human activities, the amount of water and sediment entering the Yellow River has been greatly reduced; on the other hand, the operation of the Xiaolangdi Reservoir has led to the continuous coarsening of the bed sand in the lower reaches of the Yellow River, and the sediment transport efficiency has been significantly reduced. As a result, the use of historical empirical relationships to estimate the water demand for sediment transport in the lower reaches of the Yellow River is no longer in line with current objective conditions. For the water demand for sediment transport in the lower reaches of the Yellow River in the future, Zhang Hongwu et al. (2021, China Water Resources) used the river dynamics method to carry out relevant calculations, but this method regards the river width and water depth as fixed values without considering the dynamic adjustment of the river channel morphology. Summary of the invention
[0003] In order to solve the problems of the above-mentioned prior art, the present invention provides a method, system, equipment and medium for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory, so as to accurately calculate the water demand for sediment transport under different sediment inflows and riverbed boundary conditions in the future.
[0004] The present invention is achieved through the following technical solutions: In a first aspect, the present invention provides a method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory, comprising: S1, obtain the river channel information and sediment parameters of the target river section in the lower reaches of the Yellow River; calculate the comprehensive riverbed stability index based on the river channel information and sediment parameters; S2, based on river channel information, sediment parameters and comprehensive riverbed stability index, the equilibrium river channel parameters are calculated by simultaneously solving the water flow continuity equation, water flow motion equation and river phase relationship; S3, substitute the equilibrium river channel parameters and sediment parameters into the sediment transport equation to obtain the saturated sediment transport capacity of the equilibrium river channel; calculate the sediment transport water demand under different sediment inflow scenarios based on the saturated sediment transport capacity of the equilibrium river channel.
[0005] Preferably, in S1, the river channel information includes flat-shoal flow, riverbed gradient and hydraulic gradient; and the sediment parameters include bed sediment median particle size and suspended sediment median particle size.
[0006] Further, in S1, the comprehensive stability index of the riverbed is calculated based on the river channel information and sediment parameters, specifically as follows:
[0007] In the formula, is the comprehensive stability index of the riverbed, is the riverbed slope, is the median diameter of the bed material, is the acceleration due to gravity, is the kinematic viscosity of the water flow, are all coefficients, calibrated through the measured hydrological data and sediment data of the target river reach by formula (1).
[0008] Further, in S2, the water flow continuity equation, water flow motion equation, and river regime relationship are respectively: Water flow continuity equation:
[0009] Water flow motion equation:
[0010] River regime relationship:
[0011] Among them,
[0012]
[0013] In the formula, is the bankfull discharge, is the river width under bankfull conditions, is the water depth under bankfull conditions, is the hydraulic gradient, is the riverbed slope, is the flow velocity under bankfull conditions, radius, n is the roughness coefficient, γ is the unit weight of water flow, γ s is the unit weight of sediment, is the Froude number, is the friction parameter corresponding to the bed material, is the comprehensive stability index of the riverbed; The equilibrium river channel parameters are the river width , water depth and flow velocity .
[0014] Further, in S3, the sediment transport equation is:
[0015] Among them,
[0016]
[0017]
[0018] Substitute into the following formula to obtain the saturated sediment transport capacity of the equilibrium river channel:
[0019] In the formula, is the sediment carrying capacity of the water flow, is the sediment volume concentration, is the von Karman constant of the sediment in clear water ( = 0.4), is the von Karman constant of the sediment in the muddy water, is the average settling velocity of the non-uniform sediment in clear water, is the group settling velocity of the sediment in the muddy water, is the unit weight of the water flow, is the unit weight of the sediment, is the unit weight of the muddy water, is the saturated sediment transport capacity of the equilibrium river channel, is the bankfull discharge, is the water depth under the bankfull condition, is the flow velocity under the bankfull condition, is the median diameter of the bed sediment, is the median diameter of the suspended sediment.
[0020] Furthermore, in S3, calculating the water demand for sediment transport under different incoming sediment scenarios according to the saturated sediment transport capacity of the equilibrium river channel specifically includes:
[0021]
[0022] Wherein, is the duration of saturated sediment transport, is the incoming sediment volume, is the water demand for sediment transport.
[0023] Preferably, the target river reach is the bottleneck reach of the flow capacity in the lower Yellow River.
[0024] In a second aspect, the present invention provides a system for calculating the water demand for sediment transport in the lower Yellow River based on the theory of riverbed evolution equilibrium, including: An information acquisition module, configured to acquire the river channel information and sediment parameters of the target river reach in the lower Yellow River; calculate the comprehensive riverbed stability index according to the river channel information and sediment parameters; The balanced river channel parameter calculation module is used to calculate the balanced river channel parameters based on the river channel information, sediment parameters and comprehensive riverbed stability index by simultaneously solving the water flow continuity equation, water flow motion equation and river phase relationship; The sediment transport water demand calculation module is used to substitute the equilibrium river channel parameters and sediment parameters into the sediment transport equation to obtain the saturated sediment transport capacity of the equilibrium river channel; and calculate the sediment transport water demand under different sediment delivery scenarios based on the saturated sediment transport capacity of the equilibrium river channel.
[0025] In a third aspect, the present invention provides a computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory when executing the computer program.
[0026] In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory as described above.
[0027] Compared with the prior art, the present invention has the following beneficial effects: The lower reaches of the Yellow River are alluvial rivers. When the water and sediment or boundary conditions in the basin change, it will adjust its own shape to reach a state of equilibrium that adapts to the changes. Based on this understanding, the present invention starts from the theory of riverbed evolution equilibrium and establishes a method for calculating the water demand for sediment transport in the lower reaches of the Yellow River with appropriate adjustment of water, sediment and riverbed, which serves the estimation of water demand for sediment transport in the future. This method objectively and accurately reflects the dynamic response relationship between water and sediment situation, riverbed evolution and water demand for sediment transport. It can not only use the flat beach flow to smoothly and efficiently transport sand into the sea in the future, but also ensure that the lower reaches of the Yellow River have a certain flood discharge capacity. This is of great significance for the rational allocation of water resources in the Yellow River basin and the scientific management of sediment in the downstream river in the future. BRIEF DESCRIPTION OF THE DRAWINGS
[0028] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.
[0029] Figure 1 Sediment inflow obtained by fitting based on the water and sediment data of the Gaocun-Sunkou river section with a sedimentation ratio of ≤10% during flood season from 1950 to 1985 Water demand for sediment transport The relationship curve between Figure 2 is the median grain size of bed sediment in the Gaosun River section for different years; Figure 3 is the average settling velocity data corresponding to different median grain sizes of suspended sediment. Specific implementation manners
[0030] The following uses specific specific examples to illustrate the implementation manners of the present invention. Those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. The present invention can also be implemented or applied through other different specific implementation manners. Various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention.
[0031] It should be noted that the process equipment or devices not specifically noted in the following embodiments all adopt conventional equipment or devices in the art.
[0032] It should be noted that the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products, or devices. Moreover, unless otherwise specified, the numbers of the method steps are only convenient tools for identifying the method steps, rather than limiting the arrangement order of the method steps or the scope in which the present invention can be implemented. The change or adjustment of their relative relationships, without substantially changing the technical content, should also be regarded as the scope in which the present invention can be implemented.
[0033] The method for calculating the water requirement for sediment transport in the lower Yellow River based on the theory of riverbed evolution equilibrium of the present invention has the following specific operation steps: Step 1: Determine the target river section, and determine the boundary conditions of the equilibrium river channel in the future period by sorting and analyzing historical measured data, involving the screening of information such as the "hump" river section (i.e., the flow bottleneck river section), floodplain discharge, riverbed slope, hydraulic slope, median grain size of bed sediment, median grain size of suspended sediment, and comprehensive stability index of the riverbed.
[0034] Step 1.1 Determine the target river section The lower reaches of the Yellow River are a complex riverbed of the "river channel + beach" type. The river channel, as the core channel for flood discharge and sand transport, plays a decisive role in the downstream flood discharge capacity. The flow scale of the main channel of the lower reaches of the Yellow River is "large at both ends and small in the middle", and Sunkou is the location of the flow bottleneck, which means that maintaining the scale of the main channel of the Sunkou river section is the key to ensuring the flood discharge capacity of the downstream river. In addition, measured data show that the unit water demand for sand transport in the Gaocun-Aishan river section during the flood season is greater than that of the upstream and downstream rivers. The water demand for sand transport determined by the present invention must not only ensure that all the sand from the upstream is transported to the sea, but also ensure that a certain scale of the main channel is maintained in the lower reaches of the Yellow River. Therefore, the Gaocun-Sunkou river section (referred to as the Gaosun river section) is selected as the target river section of the embodiment of the present invention.
[0035] Step 1.2 Determine boundary conditions and obtain river channel information and sediment parameters Considering the future flood discharge capacity demand of the lower Yellow River and the impact of the combined use of reservoirs on the flow capacity of the river channel, a representative flat-shoal flow was selected. The riverbed gradient and hydraulic gradient were determined based on the measured hydrological data of the target river section. The representative median particle size of the bed sand was selected based on the measured sediment data of the target river section. The representative median particle size of suspended sediment was selected based on the Yellow River Sediment Bulletin. The riverbed gradient and the median particle size of the bed sand were substituted into the empirical formula (1) to calculate the comprehensive riverbed stability index of the target river section. The comprehensive riverbed stability index can approximately characterize the river phase relationship in the lower reaches of the Yellow River.
[0036] (1) In the formula, is the riverbed gradient, is the median particle size of bed sand, is the acceleration due to gravity, is the kinematic viscosity of water, are coefficients, which are obtained by calibration based on the measured hydrological data and sediment data of the target river section and formula (1). =-1.2754, =0.609, =0.9585.
[0037] Step 2: Based on the river channel information, sediment parameters and comprehensive riverbed stability index obtained in step 1, the equilibrium river channel parameters such as river width, water depth and flow velocity are calculated by combining the water flow continuity equation, water flow motion equation and river phase relationship.
[0038] Substitute the representative flat-shoal flow, riverbed gradient, median particle size of bed sand, and comprehensive riverbed stability index obtained in step 1 into the equation group consisting of "water flow continuity equation, water flow motion equation, and river phase relationship" to calculate the equilibrium river channel parameters under the flat-shoal condition, including river width, water depth, and flow rate.
[0039] Continuity equation for water flow: (2) Flow motion equation: (3) River regime relationship: (4) Among them, (5) (6) In the formula, the subscript " bf " represents the bankfull situation, is the bankfull discharge, is the river width under bankfull conditions, is the water depth under bankfull conditions, is the hydraulic gradient, is the flow velocity under bankfull conditions, is the hydraulic radius, n is the roughness coefficient, γ is the unit weight of water flow, γ s is the unit weight of sediment, is the Froude number, is the friction parameter corresponding to the bed sediment (for the lower Yellow River channel, take the value of 19).
[0040] Step 3: Substitute the equilibrium river channel parameters obtained in Step 2 and the sediment parameters obtained in Step 1 into the sediment transport equation to obtain the saturated sediment transport capacity of the equilibrium river channel, and then calculate the water demand for sediment transport under different sediment inflow scenarios.
[0041] Step 3.1 Calculate the saturated sediment transport capacity ① Substitute the equilibrium river channel parameters obtained in Step 2 (i.e., water depth , flow velocity ) and the sediment parameters obtained in Step 1 (i.e., median diameter of bed sediment 、 median diameter of suspended sediment ) into the sediment transport formula (7) to obtain the sediment carrying capacity of the equilibrium river channel (corresponding to the sediment volume concentration ).
[0042] (7) Among them, (8) (9) (10) In the formula, is the Karman constant of sediment in clear water ( = 0.4), is the Karman constant of sediment in the muddy water, is the average settling velocity of non-uniform sediment in clear water (obtained according to "Hydrological Data of the Yellow River Basin"), is the collective settling velocity of sediment in the muddy water, is the unit weight of the water flow, is the unit weight of the sediment, is the unit weight of the muddy water, is the median diameter of the suspended sediment.
[0043] ② Substitute the sediment-carrying capacity of the water flow into the sediment transport rate formula , and obtain the saturated sediment transport capacity of the equilibrium river channel .
[0044] Step 3.2 Calculate the water demand for sediment transport ① For a given sediment inflow , substitute the saturated sediment transport capacity obtained in Step 3.1 into the formula , and the saturated sediment transport duration can be obtained.
[0045] ② Substitute the saturated sediment transport duration into the formula , and the water demand for sediment transport under different sediment inflow situations can be obtained.
[0046] Embodiment To verify the reliability of the present invention, a comparison was made between the historical measured values of the water demand for sediment transport and the calculated values of the present invention: ① Based on the water and sediment data of the flood period from 1950 to 1985 in the Gaocun-Sunkou reach with a sedimentation ratio ≤ 10% (a sedimentation ratio ≤ 10% can be approximately regarded as scouring and silting balance), a linear relationship between the sediment inflow and the water demand for sediment transport was fitted as follows: , as shown in Figure 1 . ② The multi-year average of the median diameter of the bed sediment at the Gaocun and Sunkou stations from 1950 to 1985 was about 0.091 mm, which was used as the representative value of the bed sediment particle size in the Gaocun-Sunkou reach during this period; the flood discharge of the Gaocun-Sunkou reach from 1950 to 1985 was roughly between 3500 - 8000 m³ / s, and 4000 m³ / s was selected as the representative value of the flood discharge during this period. ③ According to the calculation process of the present invention, under the conditions of = 0.091 mm and = 4000 m³ / s, the water demand for sediment transport when the sediment inflow = 200 - 300 million t was obtained, and it can be seen that the calculated data points fell within the measured range, reflecting the reliability of the present invention.
[0047] To verify the feasibility of the present invention, based on the theory of riverbed evolution equilibrium, the water demand for sediment transport in the Gaocun-Sunkou reach (referred to as the Gaosun reach) of the lower Yellow River in the future was estimated, and the calculation process is described as follows.
[0048] Step 1: Determine the calculation reach and boundary conditions Step 1.1 Determine the calculation reach Select the Gaosun reach in the lower Yellow River with limited flow capacity as the calculation reach.
[0049] Step 1.2 Determine the boundary conditions ① Based on the understanding of the future water and sediment conditions entering the lower Yellow River, relevant studies believe that the flood-discharge capacity of the main river channel in the lower Yellow River needs to be maintained at about 4000 m³ / s for a long time to ensure the basic water and sediment transport channel. Therefore, 4000 m³ / s is selected as the representative value of the flood-discharge capacity in the future. ② Measured data show that the change in the riverbed slope in the lower Yellow River is limited, and the current riverbed slope of the Gaosun reach = 0.13‰ can be selected as the riverbed slope in the future. ③ Due to the violent change in the water surface slope caused by the transmission of flood waves during the flood period, the average water surface slope of 0.169‰ during the flood period is selected as the hydraulic slope of the Gaosun reach in the future. ④ Select the median grain size of the bed sediment commonly found in the lower Yellow River = 0.10 mm as the median grain size of the bed sediment in the Gaosun reach in the future, which is in line with the continuous refinement trend of the median grain size of the bed sediment in the Gaosun reach in recent years (as Figure 2 shown). ⑤ According to the Yellow River sediment bulletin, the multi-year average value of the median grain size of the suspended sediment at each station in the lower Yellow River is stable at 0.02 mm. Therefore, = 0.02 mm is selected as the median grain size of the suspended sediment in the Gaosun reach in the future; according to the "Hydrological Data of the Yellow River Basin", the average settling velocity of the non-uniform sediment corresponding to = 0.02 mm in clear water is = 0.193 cm / s (as Figure 3 shown). ⑥ Substitute the riverbed slope = 0.13‰ and the median grain size of the bed sediment = 0.10 mm into the empirical formula (1) to obtain the riverbed comprehensive stability index = 8.4.
[0050] Step 2: Calculate the equilibrium river channel parameters Substitute the flood-discharge capacity = 4000 m 3 / s, riverbed slope = 0.13‰, hydraulic slope = 0.169‰, median grain size of the bed sediment = 0.10 mm, and riverbed comprehensive stability index Substitute = 8.4 into the system of equations consisting of the "water flow continuity equation, water flow motion equation, and river regime relationship" to calculate the river width under the flat floodplain condition. = 599.7 m, water depth = 2.85 m and flow velocity = 2.34 m / s.
[0051] Step 3: Calculate the saturated sediment transport capacity and the water demand for sediment transport Step 3.1 Calculate the saturated sediment transport capacity ① Substitute the equilibrium river channel parameters obtained in Step 2 (i.e., water depth = 2.85 m, flow velocity = 2.34 m / s) and the sediment parameters obtained in Step 1 (i.e., median bed sediment diameter = 0.10 mm 、 median suspended sediment diameter = 0.02 mm, average settling velocity = 0.193 cm / s) into the sediment transport formula (7) to obtain the sediment carrying capacity of the equilibrium river channel = 33.0 kg / m 3 .
[0052] ② Substitute the sediment carrying capacity = 33.0 kg / m 3 into the sediment transport rate formula to obtain the saturated sediment transport capacity of the equilibrium river channel = 132000 kg / s = 132 t / s.
[0053] Step 3.2 Calculate the water consumption for sediment transport ① Relevant research shows that the average annual sediment inflow in the lower reaches of the Yellow River in the future period is 245 million tons, so the present invention sets = 250 million tons as the long-term sediment inflow. Substitute the saturated sediment transport capacity = 132 t / s obtained in Step 3.1 into the formula to obtain the saturated sediment transport duration ≈ 22 d.
[0054] ② Substitute the saturated sediment transport duration ≈ 22 d into the formula to obtain the water demand for sediment transport when the sediment inflow = 250 million tons 3 .
[0055] Considering that the above case only calculated the sediment transport water demand under typical working conditions, the sediment transport water demand under different sediment inflows and riverbed boundary conditions in the future period is further analyzed and compared, as shown in Table 1. It can be seen that with the increase of flat beach flow and the coarsening of bed sand particle size, the sediment transport water demand required to transport the same weight of sediment shows an increasing trend.
[0056] Table 1
[0057] The following are device embodiments of the present invention, which can be used to perform method embodiments of the present invention. For details not disclosed in the device embodiments, please refer to the method embodiments of the present invention.
[0058] In yet another embodiment of the present invention, a system for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory is provided, comprising: The information acquisition module is used to obtain the river channel information and sediment parameters of the target river section in the lower reaches of the Yellow River; the comprehensive riverbed stability index is calculated based on the river channel information and sediment parameters; The balanced river channel parameter calculation module is used to calculate the balanced river channel parameters based on the river channel information, sediment parameters and comprehensive riverbed stability index by simultaneously solving the water flow continuity equation, water flow motion equation and river phase relationship; The sediment transport water demand calculation module is used to substitute the equilibrium river channel parameters and sediment parameters into the sediment transport equation to obtain the saturated sediment transport capacity of the equilibrium river channel; and calculate the sediment transport water demand under different sediment delivery scenarios based on the saturated sediment transport capacity of the equilibrium river channel.
[0059] In another embodiment of the present invention, a computer device is provided. The computer device includes a processor and a memory. The memory is used to store a computer program, and the computer program includes program instructions. The processor is used to execute the program instructions stored in the computer storage medium. The processor may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing core and control core of the terminal, and is suitable for implementing one or more instructions. Specifically, it is suitable for loading and executing one or more instructions in the computer storage medium to implement the corresponding method flow or corresponding function. The processor described in the embodiment of the present invention can be used for the operation of the method for calculating the water requirement for sediment transport in the lower reaches of the Yellow River based on the theory of riverbed evolution equilibrium.
[0060] In another embodiment of the present invention, a storage medium is also provided, specifically a computer-readable storage medium (Memory). The computer-readable storage medium is a memory device in the computer device and is used to store programs and data. It can be understood that the computer-readable storage medium here can include both the built-in storage medium in the computer device and, of course, the extended storage medium supported by the computer device. The computer-readable storage medium provides a storage space, and the operating system of the terminal is stored in this storage space. And, one or more instructions suitable for being loaded and executed by the processor are also stored in this storage space. These instructions can be one or more computer programs (including program codes). It should be noted that the computer-readable storage medium here can be a random access memory (RAM), or a non-volatile memory, such as at least one disk memory. One or more instructions stored in the computer-readable storage medium can be loaded and executed by the processor to implement the corresponding steps of the method for calculating the water requirement for sediment transport in the lower reaches of the Yellow River based on the theory of riverbed evolution equilibrium in the above embodiments.
[0061] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system or a computer program product. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk memory, optical memory, etc.) that contain computer-usable program code.
[0062] The present invention is described with reference to the flowcharts and / or block diagrams of methods, systems (devices) and computer program products according to the embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram, as well as the combination of flows and / or blocks in the flowchart and / or block diagram, can be realized by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for realizing the functions specified in one Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.
[0063] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing devices to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured article including an instruction device, and the instruction device realizes the functions specified in one Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.
[0064] These computer program instructions can also be loaded onto a computer or other programmable data processing devices, so that a series of operation steps are executed on the computer or other programmable devices to generate a computer-implemented process. Therefore, the instructions executed on the computer or other programmable devices provide steps for realizing the functions specified in one Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.
[0065] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: still can modify the specific implementation manners of the present invention or make equivalent replacements. Any modification or equivalent replacement without departing from the spirit and scope of the present invention should be covered by the protection scope of the claims of the present invention.
Claims
1. The method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory is characterized by: include: S1, obtain the river channel information and sediment parameters of the target river section in the lower reaches of the Yellow River; The comprehensive riverbed stability index is calculated based on river channel information and sediment parameters; S2, based on river channel information, sediment parameters and comprehensive riverbed stability index, the equilibrium river channel parameters are calculated by simultaneously solving the water flow continuity equation, water flow motion equation and river phase relationship; S3, substitute the equilibrium river channel parameters and sediment parameters into the sediment transport equation to obtain the saturated sediment transport capacity of the equilibrium river channel; calculate the sediment transport water demand under different sediment inflow scenarios based on the saturated sediment transport capacity of the equilibrium river channel.
2. The method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory according to claim 1 is characterized in that: In S1, the river channel information includes flat-shoal flow, riverbed gradient and hydraulic gradient; the sediment parameters include bed sediment median particle size and suspended sediment median particle size.
3. The method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory according to claim 2 is characterized in that: In S1, the comprehensive riverbed stability index is calculated based on the river channel information and sediment parameters, specifically: In the formula, It is an indicator of comprehensive riverbed stability. is the riverbed gradient, is the median particle size of bed sand, is the acceleration due to gravity, is the kinematic viscosity of water, are coefficients, which are calibrated by formula (1) based on the measured hydrological data and sediment data of the target river section.
4. The method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory according to claim 2 is characterized in that: In S2, the water flow continuity equation, water flow motion equation and river phase relationship are respectively: Continuity equation for water flow: The equation of motion for water flow is: River phase relationship: in, In the formula, It is flat-beach flow. is the width of the river in a flat situation. It is the water depth under flat conditions. is the hydraulic gradient, is the riverbed gradient, is the flow velocity under flat conditions, radius, n is the roughness coefficient, γ is the water flow density, γ s is the bulk density of sediment, is the Froude number, is the friction parameter corresponding to the bed sand, It is an indicator of comprehensive riverbed stability; The balanced river channel parameter is the river width , water depth and flow rate .
5. The method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory according to claim 2 is characterized in that: In S3, the sediment transport equation is: in, Substitute the following formula to obtain the saturated sediment transport capacity of the equilibrium river: In the formula, It is the force of water carrying sand. is the sediment volume concentration, is the Karman constant of sediment in clear water ( =0.4), is the Karman constant for sediment in turbid water, is the average sinking velocity of non-uniform sand in clear water, is the mass sinking speed of sediment in muddy water. is the water flow density, is the bulk density of sediment, It is the density of muddy water. It is the saturated sediment transport capacity of the balanced river channel. It is flat beach flow. It is the water depth under flat conditions. is the flow velocity under flat conditions, is the median particle size of bed sand, is the median particle size of suspended sediment.
6. The method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory according to claim 5 is characterized in that: In S3, the water demand for sediment transport under different sediment inflow scenarios is calculated based on the saturated sediment transport capacity of the balanced river channel, specifically: in, is the saturated sediment transport duration, It's the amount of sand. It is the water requirement for sediment transport.
7. The method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory according to claim 1 is characterized in that: The target river section is the bottleneck section of the flow capacity of the lower reaches of the Yellow River.
8. The calculation system of water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory is characterized by: include: Information acquisition module, used to obtain river channel information and sediment parameters of the target river section in the lower reaches of the Yellow River; The comprehensive riverbed stability index is calculated based on river channel information and sediment parameters; The balanced river channel parameter calculation module is used to calculate the balanced river channel parameters based on the river channel information, sediment parameters and comprehensive riverbed stability index by simultaneously solving the water flow continuity equation, water flow motion equation and river phase relationship; The sediment transport water demand calculation module is used to substitute the equilibrium river channel parameters and sediment parameters into the sediment transport equation to obtain the saturated sediment transport capacity of the equilibrium river channel; and calculate the sediment transport water demand under different sediment delivery scenarios based on the saturated sediment transport capacity of the equilibrium river channel.
9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory as described in any one of claims 1 to 7 is implemented.
10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by the processor, the method for calculating the water demand for sediment transport in the lower reaches of the Yellow River based on the riverbed evolution equilibrium theory as described in any one of claims 1 to 7 is implemented.