CN111985164A - River flow prediction method based on fluid power and piecewise Manning formula - Google Patents
River flow prediction method based on fluid power and piecewise Manning formula Download PDFInfo
- Publication number
- CN111985164A CN111985164A CN202010662953.XA CN202010662953A CN111985164A CN 111985164 A CN111985164 A CN 111985164A CN 202010662953 A CN202010662953 A CN 202010662953A CN 111985164 A CN111985164 A CN 111985164A
- Authority
- CN
- China
- Prior art keywords
- river
- formula
- section
- flow
- velocity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 239000012530 fluid Substances 0.000 title claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 29
- 238000004080 punching Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 4
- 238000012271 agricultural production Methods 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention discloses a river flow prediction method based on hydrodynamic force and a segmented Manning formula, which can be used for solving river flow under some special conditions by only needing one stress body and according to river depth and stress body included angle.
Description
Technical Field
The invention belongs to the field of river flow prediction, and particularly relates to a river flow prediction method based on fluid power and a piecewise Manning formula.
Background
With the increase of water consumption in industrial and agricultural production and people's life, "the contradiction between water supply and demand is increasingly sharp" the planned water use and water allocation work is more and more important. The proposal and the attention of water resource problems are that the flow of an open channel can be rapidly and accurately measured in many occasions, and the method is prominently reflected in the aspects of water quantity calculation, pollutant total quantity control and water resource scheduling and distribution of irrigation areas and water diversion projects. The method consistently used in China mainly adopts a flow velocity area method, and manual operation is required, and the requirement of automation cannot be met although the measurement precision is high. Although the water level flow relation method used in a small amount can automatically measure the flow, the accuracy is not high, and the requirement of accurate water quantity measurement cannot be met. In addition, there are few methods for automatically measuring flow by means of imported instruments (such as acoustic doppler profilers), which are expensive and high in maintenance cost, and thus cannot meet the requirements for automatic measurement of flow of open channels with a large number of wide areas.
Therefore, the river flow prediction method based on the fluid power and the piecewise Manning formula is provided, and the method needs a stress body, places the stress body in water, measures an included angle, a river depth and a distance from an offshore edge, and then respectively substitutes the included angle, the river depth and the distance into the formula of the invention to calculate the river flow.
Disclosure of Invention
The invention aims to provide a river flow prediction method based on fluid power and a piecewise Manning formula, which can solve the problem of dependence on an import instrument by only needing one stress body under some special conditions and solving the river flow according to the river depth and the stress body included angle.
In order to solve the problems, the technical scheme is as follows:
a river flow prediction method based on fluid power and a piecewise Manning formula comprises the following steps:
s1: putting a stress body on the surface of a river in a crossing mode, and measuring an included angle alpha between the stress body and the vertical direction and the length d of the stress body exposed out of the water surface;
s2: substituting the length, width, weight, thickness and included angle of the stress body into a formula to obtain the average flow velocity of the position of the stress body; the formula is as follows:
s3: dividing the river section into N small sections, and substituting the distance between the small sections and the shore and the water depth into the following formula Vi=VFlat plateCalculating the average flow velocity of each small section;
s4: the cross-section river flow is obtained according to a flow velocity area method, and the formula is as follows: qj=Sj*Vj;
S5: and summing the river flows of the N sections to predict the river flow.
Wherein: qjFlow rate of j-position section, SjIs the area of the cross-section at the j position, VjThe average velocity of the j position section.
Preferably, in S2, before determining the average flow velocity of the force-bearing body, the flow velocity at a certain point is determined, and then the flow velocity functional relationship at a certain point is obtained by:
Wherein: the length of the force-bearing body is L, the length d above the water surface is W, the width is G, the thickness is G, the weight of the water is gamma, and the mass of the force-bearing body is rho.
the impulse force of the stress body is as follows: l isPunching machine=(L-d)cosα;
With the upper end of the object as the center, the following can be obtained according to the law of conservation of moment:
Lheavy loadFHeavy loadsinα-LPunching machineFPunching machinecosα-LFloating bodyFFloating bodysinα=0
Combining the above formula to obtain:
preferably, the surface velocity has a linear relationship with the average flow velocity based on a river whose river flow velocity is not affected by two walls as follows
VFlat plate=0.8871VWatch (A)
Preferably, the hydraulic radius at any one pointFor any riverLeft wet cycle is L1Right wet cycle is L2Left side water passing area is S1Right side water passing area is S2(ii) a Hydraulic radius of any placeThe average flow velocity at any position can be obtained by substituting the slope and the roughness into a formula, wherein the specific formula is as follows:
wherein n is roughness, RRIs the hydraulic radius and J is the slope.
Preferably, the above formula needs to measure the gradient and the roughness, but in the practical process, the roughness is inversely calculated according to the Manning formula, the speed of each position can be compared with the speed of a certain position, two unknown parameters of the slope (J) and the roughness (n) can be eliminated,
has the advantages that:
according to the moment conservation law, in the river flow velocity measurement, the average of the position can be calculated only by measuring the included angle of the stress body, and two parameters of roughness and slope drop in a Manning formula are effectively avoided when the average flow velocity of each small section is calculated.
Drawings
Fig. 1 is a schematic view of a stress body suspended in a river in a crossing manner.
Fig. 2 is a schematic view of a river.
Fig. 3 is a schematic sectional view.
Detailed Description
The present invention is further illustrated by the following specific examples, which are intended to be illustrative only and not limiting.
Fig. 2 and 3 show examples of river widths of 70m, with the following other parameters:
distance from left bank | 10 | 20 | 30 | 40 | 50 | 60 |
Depth of field | 3 | 3.3 | 3.5 | 3.3 | 3.1 | 3 |
Left side water passing area S1 | 30 | 63 | 98 | 131 | 162 | 192 |
Left Wet week L1 | 9.5 | 18.9 | 28.3 | 37.8 | 47.6 | 57.1 |
Right side water passing area S2 | 192 | 162 | 129 | 94 | 61 | 30 |
Right humidity week L2 | 57.1 | 47.6 | 38.2 | 28.8 | 19.3 | 9.5 |
A force-bearing body (aluminum, length 1m, width 0.1m, mass 27kg) was taken, the force-bearing body was located 10m from the left bank, the length exposed on the water surface was 0.3m, and α was 30 ° (shown in fig. 1).
Substituting the parameters of the force-bearing body into the following formula:
Vflat plate=0.99
The average flow velocity determined here is the average flow velocity of the cross section where the force-receiving body is located.
The above-obtained result and river-related parameters are substituted into the following formula:
substituting the above speed resultsInto Qj=Sj*Vj;
Total flow Q ═ Sigma Q of cross section of riverj;Q=399.7。
Claims (6)
1. A river flow prediction method based on fluid power and a piecewise Manning formula is characterized by comprising the following steps:
s1: putting a stress body on the surface of a river in a crossing mode, and measuring an included angle alpha between the stress body and the vertical direction and the length d of the stress body exposed out of the water surface;
s2: substituting the length L of the force-bearing body, the length d exposed out of the water surface, the width W, the mass rho and the thickness G of the force-bearing body and the included angle into a formula to obtain the average flow velocity of the position of the force-bearing body; the formula is as follows:
s3: dividing the river section into N small sections, and substituting the distance between the small sections and the shore and the water depth into the following formula Vi=VFlat plateCalculating the average flow velocity of each small section;
s4: the cross-section river flow is obtained according to a flow velocity area method, and the formula is as follows: qj=Sj*Vj;
S5: and summing the river flows of the N sections to predict the river flow.
Wherein: qjFlow rate of j-position section, SjIs the area of the cross-section at the j position, VjThe average velocity of the j position section.
2. The method for predicting river discharge according to claim 1, wherein in step S2, before the average flow velocity of the stressed body is determined, the flow velocity at a certain point is determined, and the flow velocity functional relationship at a certain point is obtained by the following steps:
Wherein: the length of the stress body is L, the length d above the water surface is W, the width is G, the thickness is G, the weight of water is gamma, and the mass of the stress body is rho;
the impulse force of the stress body is as follows: l isPunching machine=(L-d)cosα;
With the upper end of the object as the center, the following can be obtained according to the law of conservation of moment:
Lheavy loadFHeavy loadsinα-LPunching machineFPunching machinecosα-LFloating bodyFFloating bodysinα=0
Combining the above formula to obtain:
3. the river discharge prediction method based on fluid dynamics and piecewise Manning's equation according to claim 2, wherein based on a river whose river discharge velocity is not affected by two walls, the surface velocity has a linear relationship with the average discharge velocity as follows
VFlat plate=0.8871VWatch (A)
4. The river discharge prediction method based on fluid power and piecewise Manning's equation according to claim 3,
hydraulic radius at any pointFor any riverLeft wet cycle is L1Right wet cycle is L2Left side water passing area is S1Right side water passing area is S2(ii) a The hydraulic radius, slope drop and roughness at any position are substituted into a formula to obtain the average flow speed at any position, and the specific formula is as follows:
wherein: qjFlow rate of j-position section, SjIs the area of the cross-section at the j position, VjThe average velocity of the j position section.
Wherein n is roughness, RRIs the hydraulic radius and J is the slope.
5. The river flow prediction method based on hydrodynamic force and piecewise Manning formula in claim 4, wherein the above formula is used to measure gradient and roughness, but in practice roughness is calculated according to Manning formula, and the speed of each position can be compared with the speed of a position to eliminate two unknown parameters of n and J,
6. the method of claim 1The river flow prediction method based on the fluid power and the piecewise Manning formula is characterized in that the total flow Q of the river section is sigma Qj。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662953.XA CN111985164B (en) | 2020-07-10 | 2020-07-10 | River flow prediction method based on hydrodynamic force and sectional Manning formula |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662953.XA CN111985164B (en) | 2020-07-10 | 2020-07-10 | River flow prediction method based on hydrodynamic force and sectional Manning formula |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111985164A true CN111985164A (en) | 2020-11-24 |
CN111985164B CN111985164B (en) | 2024-06-18 |
Family
ID=73439014
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010662953.XA Active CN111985164B (en) | 2020-07-10 | 2020-07-10 | River flow prediction method based on hydrodynamic force and sectional Manning formula |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111985164B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117571505A (en) * | 2024-01-12 | 2024-02-20 | 水利部交通运输部国家能源局南京水利科学研究院 | Device and method for measuring critical shear force of fine groove erosion |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030000980A (en) * | 2001-06-27 | 2003-01-06 | 한국수자원공사 | Method of Stream Liquid Velocity Determination and Apparatus for thereof |
KR101029282B1 (en) * | 2010-03-30 | 2011-04-18 | 한국건설기술연구원 | Method for estimation of discharge using the measured stage data at the cross-sections with lateral change of roughness |
CN104132710A (en) * | 2014-06-27 | 2014-11-05 | 华北水利水电大学 | Water level flow sensor, water level flow monitoring system and water level flow monitoring method |
CN106021854A (en) * | 2016-05-06 | 2016-10-12 | 西安理工大学 | A natural river course cross section stage-discharge relation determining method |
CN106777775A (en) * | 2017-01-10 | 2017-05-31 | 清华大学 | A kind of neural net method based on many section water level forecast river discharges |
CN109840368A (en) * | 2019-01-17 | 2019-06-04 | 广东省水利水电科学研究院 | A kind of irrigated area channel stage discharge relation rating method |
CN110672163A (en) * | 2019-08-21 | 2020-01-10 | 中国水利水电科学研究院 | Method for testing ice period flow of canal in cold region |
-
2020
- 2020-07-10 CN CN202010662953.XA patent/CN111985164B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030000980A (en) * | 2001-06-27 | 2003-01-06 | 한국수자원공사 | Method of Stream Liquid Velocity Determination and Apparatus for thereof |
KR101029282B1 (en) * | 2010-03-30 | 2011-04-18 | 한국건설기술연구원 | Method for estimation of discharge using the measured stage data at the cross-sections with lateral change of roughness |
CN104132710A (en) * | 2014-06-27 | 2014-11-05 | 华北水利水电大学 | Water level flow sensor, water level flow monitoring system and water level flow monitoring method |
CN106021854A (en) * | 2016-05-06 | 2016-10-12 | 西安理工大学 | A natural river course cross section stage-discharge relation determining method |
CN106777775A (en) * | 2017-01-10 | 2017-05-31 | 清华大学 | A kind of neural net method based on many section water level forecast river discharges |
CN109840368A (en) * | 2019-01-17 | 2019-06-04 | 广东省水利水电科学研究院 | A kind of irrigated area channel stage discharge relation rating method |
CN110672163A (en) * | 2019-08-21 | 2020-01-10 | 中国水利水电科学研究院 | Method for testing ice period flow of canal in cold region |
Non-Patent Citations (2)
Title |
---|
MOHAMMAD SHARIFI: "《Determination of river design discharge (Tar River case study)》", 《JOURNAL OF WATER AND CLIMATE CHANGE》 * |
曹浩: "《大样本数据模型方法在中小河流流量测验中的应用》", 《长江科学院院报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117571505A (en) * | 2024-01-12 | 2024-02-20 | 水利部交通运输部国家能源局南京水利科学研究院 | Device and method for measuring critical shear force of fine groove erosion |
CN117571505B (en) * | 2024-01-12 | 2024-03-22 | 水利部交通运输部国家能源局南京水利科学研究院 | Device and method for measuring critical shear force of fine groove erosion |
Also Published As
Publication number | Publication date |
---|---|
CN111985164B (en) | 2024-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108254032B (en) | River ultrasonic time difference method flow calculation method | |
CN103353908B (en) | A kind of pipe resistance coefficient Method for Accurate Calculation based on numerical computations | |
CN106446411B (en) | Flow velocity, impact force spatial distribution measurement method inside viscous mud-flow | |
CN107490410A (en) | It is a kind of based on represent vertical line point flow velocity measurement Gao Hong contactless flow measurement method | |
Bagheri et al. | Application of free vortex theory to estimating discharge coefficient for sharp-crested weirs | |
CN106886614B (en) | A kind of measuring method of the river manning roughness based on steady nonuniform flow | |
DE602005024121D1 (en) | METHOD AND DEVICE FOR MEASURING PARAMETERS OF A COATED RIVER | |
CN102023034A (en) | Novel open channel water metering method based on multipoint water levels | |
Afzalimehr et al. | Non-uniform flow over cobble bed with submerged vegetation strip | |
CN113737710B (en) | Estimation method of sediment transport rate of natural river bed load | |
CN111985164A (en) | River flow prediction method based on fluid power and piecewise Manning formula | |
CN106869077A (en) | The method that cold region river freezeup period flow is estimated using Channel Group's phase stage discharge relation | |
CN110006511B (en) | Correction method for actually measured water level of river pressure type water level meter | |
Li et al. | Influence of auxiliary geometric parameters on discharge capacity of piano key weirs | |
CN116337163A (en) | Automatic calibration method for flow measurement system by time difference method of arcsine vertical flow velocity distribution | |
CN113505339B (en) | Index model-based two-dimensional flow field prediction method for riverway with non-submerged vegetation | |
CN112036092B (en) | River flow prediction method based on river intermediate surface speed and river width relation | |
CN106320257A (en) | Lake and reservoir channel storage curve determining method based on hydrometry | |
CN105547638A (en) | Thin layer water flow rolling wave measurement system based on pressure sensor and method thereof | |
Lee et al. | Experiments on hydraulic relations for flow over a compound sharp-crested weir | |
CN109371897B (en) | Fishway inlet water supplementing and energy dissipating model test device and test method | |
CN206756872U (en) | Horizontal ADCP flow measurement platforms | |
CN116735150A (en) | Device for measuring single-width flow transverse distribution of vegetation water flow of water tank | |
Mohanty et al. | Energy and momentum coefficients for wide compound channels | |
CN105466527A (en) | Thin sheet flow roll wave measurement system and method based on electromagnetic sensors |
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 |