CN111666667A - Method for determining flow of riverbed making of swimming river - Google Patents
Method for determining flow of riverbed making of swimming river Download PDFInfo
- Publication number
- CN111666667A CN111666667A CN202010466315.0A CN202010466315A CN111666667A CN 111666667 A CN111666667 A CN 111666667A CN 202010466315 A CN202010466315 A CN 202010466315A CN 111666667 A CN111666667 A CN 111666667A
- Authority
- CN
- China
- Prior art keywords
- flow
- bed
- river
- making
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention relates to a method for determining the bed-building flow of an oscillating river channel, which belongs to the field of hydraulic engineering, and the IPC classification number of the method is E02B 1/00. The invention absorbs the basic idea of the classical makavjer method for analyzing and calculating the bed-building flow, and two important parameters are defined: flow level length parameter Δ Q1And flow level movement distance parameter Δ Q2The flow magnitude moving accumulation bed making capacity is obtained sequentially, the flow process data are fully and effectively utilized, and errors caused by manual dividing and distributing of the flow magnitude are avoided. Flow level length parameter Δ Q1Flow level movement distance parameter Δ Q2Both the parameters are easy to select objectively, the change of the two parameters in a large range has no influence on the determination of the first and second bed making flow, and the calculated bed making flow is more objective and accurate.
Description
Technical Field
The invention relates to a method for determining the bed-building flow of an oscillating river channel, which belongs to the field of hydraulic engineering, and the IPC classification number of the method is E02B 1/00.
Background
Most of the wandering rivers develop in plain areas and are characterized in that the mainstream swings indefinitely and the river situation changes complicatedly. In plain areas, the population is dense, and in addition, the flow process of the swimming riverways is changed violently, so that the flood control pressure of the riverways is huge, and particularly the downstream river sections of the yellow rivers are most famous. When river flood control and regulation are carried out, a key flow parameter equivalent to the comprehensive bed-building effect in the flow process for many years, namely the bed-building flow, needs to be determined.
The current commonly used methods for determining the flow rate of a building bed comprise a flat flow method (such as Chinese patent application CN110555248A), a sand conveying rate method (such as Chinese patent CN106759063B and CN107401140B) and a sand conveying capacity method. The most classical calculation method proposed by makavejev, a presorter of soviet union, namely the makavejev method, is the most applicable to engineering practice. This method considers that the bedding effect of a certain flow is related to both the sand transport capacity of that flow and the duration of that flow. The former is proportional to the product of the m-th power of the flow Q and the drop J, and the latter is expressed in terms of the frequency P at which the flow occurs. The expression of the bed making capacity of the method is QmJP, where m is the sand transport rate G of the measured section on a log-log coordinatesThe slope of a line fitted to the corresponding flow rate Q; dividing the flow process of the river reach section changing along with time into a plurality of flow levels with equal intervals, and counting the frequency P of the flow occurring in each flow level. Finally, the process will typically result in two build flows, a first build flow and a second build flow.
A plurality of scholars adopt a Makawei Yefu method to widely analyze the change law of the bed-building flow of a downstream river channel of the yellow river, a middle-downstream river channel of a Songhua river, a river channel of the Yangtze river, a river channel of a Huaihe river and a snow-melting river channel in Xinjiang. However, the flow process of the wandering river varies drastically, and the flow extremum range is large. When the makavejeff method is adopted, the frequency P value corresponding to each flow is obviously influenced by the artificial dividing flow magnitude, and further great errors are brought to the calculation result of the bed-building flow: if the flow magnitude is small, the real value of a certain small flow magnitude is covered due to overhigh occurrence frequency, and the flow magnitude value is more and not obvious; if the flow magnitude is too large, the calculated bed forming flow error is large. At present, no unified theorem exists on how to select the flow magnitude, and most scholars select the flow magnitude according to experience during calculation.
Disclosure of Invention
The invention aims to solve the defects in the prior art, provides a water storage and sand discharge separation scheduling method in flood season, which is convenient to implement, low in cost and high in safety, aims to solve the problems, and provides a method for determining the bed forming flow of the oscillating river channel.
In order to achieve the purpose, the invention adopts the following scheme:
a method for determining the flow of a swimming riverway bed-building bed comprises the following steps:
in the formulaUpper and lower bounds of flow magnitude, respectively, Δ Q1For each set of flow level length, defined asΔQ2Is the distance traveled by the flow level, defined asniIs the number of flow levels.
[Qmin+(i-1)ΔQ2,Qmin+(i-1)ΔQ2+ΔQ1)(i=1,2,…ni) (2)
selecting Δ Q1Is QmaxInteger in hundred around 15, Δ Q2Is selected as the unit flow value.
Step 4, drawingIn the figure, the area enclosed by the figure is obtained by integrationTo pairPerforming a normalization process, i.e.Drawing Qi~ΦiGraph, finding phi in the graphiThe flow rate corresponding to the maximum value of (1) is the first bed-making flow rate, find out phiiThe flow rate corresponding to the second largest value of (a) is the second bed forming flow rate.
Compared with the prior art, the invention has the beneficial effects that:
the invention absorbs the basic idea of the classical makavjer method for the analysis and calculation of the bed-building flow, and fully and effectively utilizes the flow process data by solving the flow magnitude moving accumulation bed-building capacity in sequence, thereby avoiding the error caused by dividing the flow magnitude manually. In the method for calculating the bed-making flow, each group of flow level length parameters delta Q1Flow level movement distance parameter Δ Q2Both the parameters are easy to select objectively, the change of the two parameters in a large range has no influence on the determination of the first and second bed making flow, and the calculated bed making flow is more objective and accurate.
Drawings
FIG. 1 shows the sand transport rate G of the measured cross section in the 1986-2015 dual logarithmic coordinates of the garden opening at the downstream of the yellow riversA graph of Q versus corresponding flow;
FIG. 2 shows the garden of the yellow river in 1986About 2015 years of flow level of 100m3A Makavjeff method bed-making flow analysis chart under the condition of/s;
FIG. 3 shows that the flow rate of the garden opening downstream of the yellow river is 800m between 1986 and 20153A Makavjeff method bed-making flow analysis chart under the condition of/s;
FIG. 4 is a flow level length parameter Δ Q involved in an embodiment of the present invention1And flow level movement distance parameter Δ Q2A schematic diagram;
FIG. 5 is a graph of delta Q according to an embodiment of the present invention1Is 500m3S and Δ Q2Is 1m3A bed-making flow analysis chart from 1986 to 2015 of a garden opening at the downstream of the yellow river in/s;
FIG. 6 is a graph of delta Q according to an embodiment of the present invention1Is 500m3S and Δ Q2Is 1m3A bed-making flow analysis diagram after normalization treatment in 1986-2015 of garden mouths at the downstream of the yellow river in/s;
FIG. 7 is a graph of delta Q according to an embodiment of the present invention1Is 500m3/s、ΔQ2Are respectively 1, 50, 100, 300 and 600m3A bed-making flow analysis chart from 1986 to 2015 of a garden opening at the downstream of the yellow river in/s;
FIG. 8 is a graph of delta Q according to an embodiment of the present invention2Is 1m3/s、ΔQ1Respectively 200, 300, 400, 500, 600m3The bed-making flow analysis chart of garden mouths at the downstream of the yellow river in the time of/s from 1986 to 2015.
Detailed Description
Firstly, the defects of the prior method are indicated by practical calculation examples with the attached drawings.
Taking the garden port river channel downstream of the yellow river as an example, based on the measured hydrological data of a garden port hydrological station (the south bank of the yellow river at 17 km of the suburb in the Heizhou province, Zheng city, Henan province) in 1986 to 2015, the conventional makavejeff method is adopted to determine the bed-building flow rate of the garden port river channel downstream of the yellow river: (1) the sand conveying rate G of the section on the double logarithmic coordinates from 1986 to 2015 of the garden opening is countedsThe slope of the line was fitted to the flow rate Q to give m of 2.1, as shown in FIG. 1. (2) The minimum value of the daily average flow of the garden opening from 1986 to 2015 is 7.58m3At a maximum of 7270m3And s. The measured gradient was averaged to obtain J0.0002. When the artificially divided flow level is 100m3When the flow is/s, respectively counting the garden opening to measure the daily average flow falling into each flow magnitude [0m3/s,100m3/s)、[100m3/s,200m3/s)······[7200m3/s,7300m3/s) frequency P value. (3) Drawing Q to QmJP, as shown in FIG. 2, calculates the first and second bed forming flow rates (circle mark) and (square mark) to 858m, respectively3S and 1260m3/s。
When the artificially divided flow level is 800m3When the flow is/s, the step (2) is changed to respectively count the daily average flow falling into each flow magnitude [0m3/s,800m3/s)、[800m3/s,1600m3/s)······[7200m3/s,8000m3P value of frequency,/s); the steps (1) and (3) are unchanged. FIG. 3 shows Q to Q obtainedmJP relational graph shows that the calculated first and second bed forming flow rates were 1060m, respectively3S and 5053m3/s。
It can be seen that, with the prior art, the artificial division flow magnitude has a great influence on the result of the ma kavejeff-method bed-building flow calculation: when the flow level is taken to be small, a real value of a certain small flow level is covered due to overhigh occurrence frequency, the flow level values are more, the inflection point of the bed making capacity curve is not obvious, and the subjective factor is larger when the bed making flow is determined; if the flow magnitude is large, part of inflection points cannot be calculated, and the calculation result precision is not high.
Hereinafter, a specific embodiment of the method for determining the bed-building flow rate of the oscillating river according to the present invention will be described in detail.
Similarly, the river reach analyzed by the embodiment is the river reach of the garden downstream of the yellow river, the adopted hydrological data are the data of daily average flow, daily average sand transportation rate, average ratio drop and the like of the garden mouth hydrological station from 1986 to 2015, and the calculation result can be compared with the calculation result of the makavejeff method.
The invention relates to a method for determining the flow of a swimming riverway bed, which comprises the following steps:
step 4, drawingGraph, as shown in FIG. 5, the area enclosed by the graph is obtained by integrationIn order to facilitate the comparison of the bed-making capacities of different flow magnitudes, the pairIs subjected to normalization treatmentDrawing Qi~ΦiGraph, as shown in FIG. 6, in which Φ is foundiMaximum value of 4.45 × 10-4(marked by circles in the figure) and the corresponding flow rate is the first bed-making flow rate, namely 910m3/s,Φi2.62 × 10 of the second largest peak-4(marked at block in the figure) corresponding to a second bed forming flow, 2670m3And s. Compared with the flat flow obtained by statistics from 1986 to 2015 at the garden opening, the calculated bed-making flow is consistent with the average value of the flat flow for many years, which shows that the calculation result is objective and accurate.
The invention relates to a method for determining the flow of an oscillating riverway bed, which relates to two set parameters, namely the length delta Q of each flow level1And distance of flow stage movement Δ Q2. Still based on the actual hydrological data of the garden opening section of the downstream of the yellow river from 1986 to 2015, the delta Q is fixed1=500m3Respectively calculating when Δ Q2Is 1, 50, 100, 300, 600m3Q at si~ΦiA relationship curve, as shown in FIG. 7; fixed Δ Q2=1m3Respectively calculating when Δ Q 1200, 300, 400, 500, 600m3Q at si~ΦiThe relationship is shown in FIG. 8. As can be seen, even when Δ Q is measured1And Δ Q2At a set value (Delta Q)1Is QmaxInteger in hundred around 15, Δ Q2Unit flow rate value), the first and second bed forming flow rates are not determined by Δ Q1、ΔQ2The effect of a change in value. Thus, the process of the invention Δ Q1And Δ Q2The parameter value is taken in a large range, the determination of the bed-making flow is not influenced, the calculated bed-making flow is objective and accurate, the influence of the manual dividing flow magnitude in the Makavejeff method on the calculation result is avoided, and the method has good practicability and universality.
The method of the invention fully and effectively utilizes hydrological data of the riverway by obtaining the flow level moving accumulation bed-building capability in sequence, the two related setting parameters are easy to select, the determination of the first and second bed-building flows is not influenced in a larger range, meanwhile, the flow level in the Makavejev method is prevented from being difficult to select objectively, and the calculation error caused by dividing the flow level manually is avoided, the calculation result is objective and accurate, the precision is high, and the method has good practicability and universality.
Claims (1)
1. A method for determining the flow of a riverbed making bed on a swinging river course is characterized by comprising the following steps:
step 1, defining a series of flow levels as follows:
in the formulaUpper and lower bounds of flow magnitude, respectively (according to fig. 4), Δ Q1For each set of flow level length, defined asΔQ2Is the distance traveled by the flow level, defined asniIs the number of flow levels.
Step 2, arranging the actually measured flow data of the analyzed river according to the numerical value from small to large, and finding out the minimum value QminAnd maximum value Qmax(ii) a 1 st flow stageLower boundary of (1)Assigned a value of QminA series of flow levels defined by equation (1)The expression is as follows:
[Qmin+(i-1)ΔQ2,Qmin+(i-1)ΔQ2+ΔQ1)(i=1,2,…ni) (2)
selecting Δ Q1Is QmaxInteger in hundred around 15, Δ Q2Is selected as the unit flow value.
Step 3, for the ith traffic levelDefine its characterized flow asCounting all the flow rates falling into the ith flow rate class in the analyzed river channel measured flow dataAdding the bed making capacity to obtain the characteristic bed making capacity of the ith flow magnitudeWherein n isjTo fall within the ith flow levelThe number of (2); m is the sand conveying rate G of the actually measured section on the double logarithmic coordinatessThe slope of a line fitted to the corresponding flow rate Q; j is the average slope of the river; upper bound of current magnitudeChanging to be more than or equal to the maximum value Q in the analyzed river channel actual measurement flow datamaxAnd then, the last counting is finished.
Step 4, drawingIn the figure, the area enclosed by the figure is obtained by integrationTo pairPerforming a normalization process, i.e.Drawing Qi~ΦiGraph, finding phi in the graphiThe flow rate corresponding to the maximum value of (1) is the first bed-making flow rate, find out phiiThe flow rate corresponding to the second largest value of (a) is the second bed forming flow rate.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010466315.0A CN111666667B (en) | 2020-05-29 | 2020-05-29 | Method for determining bed-making flow of wandering river |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010466315.0A CN111666667B (en) | 2020-05-29 | 2020-05-29 | Method for determining bed-making flow of wandering river |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111666667A true CN111666667A (en) | 2020-09-15 |
CN111666667B CN111666667B (en) | 2023-06-13 |
Family
ID=72384902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010466315.0A Active CN111666667B (en) | 2020-05-29 | 2020-05-29 | Method for determining bed-making flow of wandering river |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111666667B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112989565A (en) * | 2021-02-03 | 2021-06-18 | 长江水利委员会水文局 | Method for calculating bed-building flow based on water flow sand-carrying capacity |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008145356A (en) * | 2006-12-13 | 2008-06-26 | Sato Kogyo Co Ltd | Method for measuring tension of embedded rod member |
CN101713178A (en) * | 2009-09-29 | 2010-05-26 | 黄河水利委员会黄河水利科学研究院 | Method for creating three-stage flow passage of meandering riverway in lower Yellow River |
CN106759063A (en) * | 2016-12-05 | 2017-05-31 | 河海大学 | A kind of computational methods of tidal reach channel forming disahcge |
CN107401140A (en) * | 2017-08-29 | 2017-11-28 | 武汉大学 | The determination method of tributary lower reaches section channel forming disahcge |
-
2020
- 2020-05-29 CN CN202010466315.0A patent/CN111666667B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008145356A (en) * | 2006-12-13 | 2008-06-26 | Sato Kogyo Co Ltd | Method for measuring tension of embedded rod member |
CN101713178A (en) * | 2009-09-29 | 2010-05-26 | 黄河水利委员会黄河水利科学研究院 | Method for creating three-stage flow passage of meandering riverway in lower Yellow River |
CN106759063A (en) * | 2016-12-05 | 2017-05-31 | 河海大学 | A kind of computational methods of tidal reach channel forming disahcge |
CN107401140A (en) * | 2017-08-29 | 2017-11-28 | 武汉大学 | The determination method of tributary lower reaches section channel forming disahcge |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112989565A (en) * | 2021-02-03 | 2021-06-18 | 长江水利委员会水文局 | Method for calculating bed-building flow based on water flow sand-carrying capacity |
CN112989565B (en) * | 2021-02-03 | 2022-03-11 | 长江水利委员会水文局 | Method for calculating bed-building flow based on water flow sand-carrying capacity |
Also Published As
Publication number | Publication date |
---|---|
CN111666667B (en) | 2023-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105005707B (en) | A kind of stage discharge relation of Heavenly Stems and Earthly Branches stream confluence determines method | |
CN107401140B (en) | The determination method of tributary lower reaches section channel forming disahcge | |
CN110458359A (en) | It is a kind of based on the plain river network water quantity and quality optimization regulating method not cut completely under the conditions of dirt | |
CN107490410A (en) | It is a kind of based on represent vertical line point flow velocity measurement Gao Hong contactless flow measurement method | |
CN112257352A (en) | Coupling method and system of one-dimensional hydrodynamic model and two-dimensional hydrodynamic model | |
CN107045568B (en) | River course roughness inversion method based on dynamic programming successive approximation method | |
CN112595489B (en) | Method and system for calculating dynamic sand transporting capacity of river channel | |
CN114580793A (en) | Water environment optimization method and system based on multi-source complementation and local regulation | |
CN108625337B (en) | Method for determining regulation water level of sandy riverbed section below tidal current boundary | |
CN111104746B (en) | River flood beach elevation determination method based on wavelet analysis | |
CN106598918A (en) | Non-uniform designed flood calculation method based on quantile regression | |
CN113642259B (en) | River flow field reconstruction method based on multilayer flow velocity measurement | |
CN115587546A (en) | Method for identifying flood level variation trend and cause of multiple river branches under dam | |
CN111666667A (en) | Method for determining flow of riverbed making of swimming river | |
CN109145396B (en) | River course roughness zoning calibration method based on vegetation distribution | |
CN106320255B (en) | A kind of computational methods of Flood Plain Flow swale inflow-rate of water turbine | |
CN114528624A (en) | Water flow acceleration method and system for water delivery open channel | |
CN113505339B (en) | Index model-based two-dimensional flow field prediction method for riverway with non-submerged vegetation | |
CN107167188A (en) | A kind of method applied to rectangular open channel flow measurement | |
CN106320257A (en) | Lake and reservoir channel storage curve determining method based on hydrometry | |
CN113585161A (en) | Construction method of alluvial river deep body lateral migration prediction model | |
CN110847112B (en) | River flood discharge early warning method based on hydraulics simulation | |
CN109992868B (en) | River channel flood forecasting method based on heterogeneous-parameter discrete generalized Nash confluence model | |
CN117033888A (en) | Watershed confluence unit line pushing method based on segmentation base flow | |
CN115879816A (en) | Winter river ice-melting water amount analysis model and method based on alpine region |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |