CN112036092A - River flow prediction method based on relation between river intermediate surface speed and river width - Google Patents
River flow prediction method based on relation between river intermediate surface speed and river width Download PDFInfo
- Publication number
- CN112036092A CN112036092A CN202010662956.3A CN202010662956A CN112036092A CN 112036092 A CN112036092 A CN 112036092A CN 202010662956 A CN202010662956 A CN 202010662956A CN 112036092 A CN112036092 A CN 112036092A
- Authority
- CN
- China
- Prior art keywords
- river
- flow
- intermediate surface
- relation
- average
- 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 43
- 230000014509 gene expression Effects 0.000 claims description 9
- 230000005888 antibody-dependent cellular phagocytosis Effects 0.000 abstract description 3
- 239000000463 material Substances 0.000 abstract description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 11
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 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
- 239000007788 liquid Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention discloses a river flow prediction method based on the relation between river intermediate surface speed and river width, which comprises the steps of firstly fitting the relation between the intermediate surface speed and the average speed according to the existing hydrological data (average depth, average speed and intermediate surface speed), measuring the river width, and then calculating the flow by adopting a flow velocity area method; because the invention only needs to measure the speed and the width of one point, compared with a Venturi method, the method has the advantages that the required hydrological parameter precision is low, and the obtained flow precision is higher; compared with ADCP, the method is simple to operate, saves a large amount of manpower, material resources and financial resources, and is suitable for some emergency situations; compared with the Manning formula for measuring the flow, the method does not need to measure the roughness and the hydraulic gradient.
Description
Technical Field
The invention belongs to the field of wide river flow prediction, and particularly relates to a river flow prediction method based on the relation between river intermediate surface speed and river width.
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.
Disclosure of Invention
The invention provides a river flow prediction method based on the relation between river intermediate surface speed and river width, aiming at solving the problems of measuring river flow by adopting the traditional manual operation flow velocity area method.
In order to solve the problems, the technical scheme is as follows:
a river flow prediction method based on the relation between river intermediate surface speed and river width comprises the following steps:
s1: acquiring river hydrological data, wherein the hydrological data comprises average depth, average flow speed and intermediate surface speed;
s2: fitting the relation between the intermediate surface speed and the average flow speed and the relation between the average flow speed and the average depth according to the hydrological data in S1;
s3: measuring the width of the river;
s4: the flow rate was measured by the flow velocity area method.
Preferably, the intermediate superficial velocity and the average flow velocity are functionally related, and the functional expression is:
v=f(vmiddle watch)
In the formula: v is the average flow velocity, vMiddle watchIs the river intermediate surface velocity.
Preferably, the best function expression for correlation is fitted according to Matlab as follows:
preferably, the average flow velocity and the average depth are as follows:
in the formula: h is the average depth, and p and q are the generation coefficients.
Preferably, the parameters obtained by the functional relationship between the average flow velocity and the average depth need to be classified into flow grades, and the denser the flow grade classification is, the higher the accuracy is.
Preferably, the flow rate measured by the flow velocity area method is in a functional relationship of: q ═ B × v × h, in which: q is the flow and B is the width.
Has the advantages that:
the method only needs to measure the speed and the width of one point, and compared with a Venturi method, the method has the advantages that the required hydrological parameter precision is low, and the obtained flow precision is high. Compared with ADCP, the method is simple to operate, saves a large amount of manpower, material resources and financial resources, and is suitable for some emergency situations. Compared with the Manning formula for measuring the flow, the method does not need to measure the roughness and the hydraulic gradient.
The river intermediate surface speed, the average river flow speed and the average river depth are combined, and a novel river flow prediction method is provided, so that the river flow prediction method is simpler and faster.
Drawings
FIG. 1 is a schematic process flow diagram.
FIG. 2 is a graph comparing data analysis.
FIG. 3 is a relational table diagram.
Fig. 4 is a functional table diagram.
Detailed Description
The present invention is further illustrated by the following specific examples, which are intended to be illustrative only and not limiting.
Fig. 1 shows a river discharge prediction method based on the relationship between river intermediate surface speed and river width, which comprises the following steps:
s1: acquiring river hydrological data, wherein the hydrological data comprises average depth, average flow speed and intermediate surface speed;
s2: fitting the relation between the intermediate surface speed and the average flow speed and the relation between the average flow speed and the average depth according to the hydrological data in S1;
s3: measuring the width of the river;
s4: the flow rate was measured by the flow velocity area method.
The intermediate superficial velocity and the average flow velocity are functionally related, and the functional expression is as follows:
v=f(vmiddle watch)
In the formula: v is the average flow velocity, vMiddle watchIs the river intermediate surface velocity.
Fitting the best function expression for correlation according to Matlab as follows:
the average flow rate and average depth are as follows:
in the formula: h is the average depth, and p and q are the generation coefficients.
The parameters calculated by the functional relation between the average flow speed and the average depth need to be divided into flow grades, and the denser the flow grade division is, the higher the precision is.
The flow rate measured by the flow velocity area method has the functional relationship that: q ═ B × v × h, in which: q is the flow and B is the width.
The principle of the method is as follows:
(1) and fitting a function expression between the two according to the existing river intermediate surface speed and average flow speed, wherein different rivers have different function expressions.
v=f(vMiddle watch)
In the formula: v is the average flow velocity, vMiddle watchIs the river intermediate surface velocity.
(2) The absolute energy is the sum of potential energy and kinetic energy of unit mass of liquid, and the absolute energy calculation formula isThe same flow rate will have different absolute energies at different slope drops. Fig. 3 shows that, for example, in the same case of h, the larger v, the larger absolute energy is, which is caused by the river slope. Thus, it is possible to provideThere is a certain relation to the slope. According to the formula of ManningWhen the river is wide, the hydraulic radius can be regarded as the average water depth and thusAndthere is a certain relation, and according to the experimental result, it can be obtained that:
and substituting the existing average flow velocity and average depth into the following formula to obtain related parameters (the parameters obtained by the method need to be classified into flow grades, the more densely the flow grades are, the higher the precision is, namely, when the flow is between 0 and 500, the related parameters are obtained from 500 to 1000, and the like).
In the formula: h is the average depth, and p and q are the generation coefficients.
(3) The width B of the river is measured by other methods such as related instruments.
(4) The flow rate was measured by the flow velocity area method.
Q=B*v*h
In the formula: q is the flow rate.
Example (b): (1) from the existing intermediate superficial velocity and average flow velocity, the function expression with the best correlation was selected using the function fitted by Matlab (shown in fig. 4).
(2) Substituting the average depth and average flow velocity intoThe selected data is data of the flow rates 2000 to 3000.
(3) Substituting the river width, average flow velocity and average depth into Q ═ B × v × h
As shown in fig. 2, according to the above embodiment, 6 sets of station data are given, specifically 31, 32, 33, 34, 35, and 36 stations, and surface midpoint velocity, average velocity (average flow velocity), calculated velocity, average depth, calculated depth, river width, ADCP flow, calculated flow, and relative error are given, respectively, and through comparison of data analysis in fig. 2, it can be seen that the error ratio obtained by the prediction method of the present invention is small, and the demand for river flow prediction can be satisfied.
The above description is only an embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention.
Claims (6)
1. A river flow prediction method based on the relation between river intermediate surface speed and river width is characterized by comprising the following steps:
s1: acquiring river hydrological data, wherein the hydrological data comprises average depth, average flow speed and intermediate surface speed;
s2: fitting the relation between the intermediate surface speed and the average flow speed and the relation between the average flow speed and the average depth according to the hydrological data in S1;
s3: measuring the width of the river;
s4: the flow rate was measured by the flow velocity area method.
2. The river discharge prediction method based on the relation between the river intermediate surface speed and the river width as claimed in claim 1, wherein the relation between the intermediate surface speed and the average flow velocity has a functional relation, and the functional expression is as follows:
v=f(vmiddle watch)
In the formula: v is the average flow velocity, vMiddle watchIs the river intermediate surface velocity.
4. a river discharge prediction method based on the relationship between river intermediate surface speed and river width as claimed in claim 3, wherein the function relationship between the average flow velocity and the average depth is as follows:
in the formula: h is the average depth, and p and q are the generation coefficients.
5. The method as claimed in claim 4, wherein the parameters derived from the functional relationship between the average flow velocity and the average depth are classified into classes, and the more dense the classes, the higher the accuracy.
6. A river discharge prediction method based on the relationship between river intermediate surface speed and river width as claimed in claim 5, wherein the functional relationship of the discharge measured by the flow area method is: q ═ B × v × h, in which: q is the flow and B is the width.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662956.3A CN112036092B (en) | 2020-07-10 | 2020-07-10 | River flow prediction method based on river intermediate surface speed and river width relation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662956.3A CN112036092B (en) | 2020-07-10 | 2020-07-10 | River flow prediction method based on river intermediate surface speed and river width relation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112036092A true CN112036092A (en) | 2020-12-04 |
CN112036092B CN112036092B (en) | 2023-11-28 |
Family
ID=73579024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010662956.3A Active CN112036092B (en) | 2020-07-10 | 2020-07-10 | River flow prediction method based on river intermediate surface speed and river width relation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112036092B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113093531A (en) * | 2021-03-01 | 2021-07-09 | 武汉大学 | Large pipe-channel combined system emergency dispatching control method based on model predictive control |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008216010A (en) * | 2007-03-02 | 2008-09-18 | Tokyo Electric Power Co Inc:The | River flow rate calculating device and method, and computer program |
CN106595777A (en) * | 2016-12-01 | 2017-04-26 | 广西师范大学 | Calculation method for detecting flow of section of river in non-contact manner |
CN107044875A (en) * | 2017-03-07 | 2017-08-15 | 上海航征测控***有限公司 | A kind of flow-measuring method |
US20180010936A1 (en) * | 2016-07-07 | 2018-01-11 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | River discharge and depth estimation |
CN109308375A (en) * | 2018-08-20 | 2019-02-05 | 河海大学 | A kind of measuring method of the basin optimal flow rate based on landforms parameter |
CN109635435A (en) * | 2018-12-12 | 2019-04-16 | 中山大学 | A kind of natural river course stage discharge relation based on bayesian theory determines method |
-
2020
- 2020-07-10 CN CN202010662956.3A patent/CN112036092B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008216010A (en) * | 2007-03-02 | 2008-09-18 | Tokyo Electric Power Co Inc:The | River flow rate calculating device and method, and computer program |
US20180010936A1 (en) * | 2016-07-07 | 2018-01-11 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | River discharge and depth estimation |
CN106595777A (en) * | 2016-12-01 | 2017-04-26 | 广西师范大学 | Calculation method for detecting flow of section of river in non-contact manner |
CN107044875A (en) * | 2017-03-07 | 2017-08-15 | 上海航征测控***有限公司 | A kind of flow-measuring method |
CN109308375A (en) * | 2018-08-20 | 2019-02-05 | 河海大学 | A kind of measuring method of the basin optimal flow rate based on landforms parameter |
CN109635435A (en) * | 2018-12-12 | 2019-04-16 | 中山大学 | A kind of natural river course stage discharge relation based on bayesian theory determines method |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113093531A (en) * | 2021-03-01 | 2021-07-09 | 武汉大学 | Large pipe-channel combined system emergency dispatching control method based on model predictive control |
CN113093531B (en) * | 2021-03-01 | 2022-03-15 | 武汉大学 | Large pipe-channel combined system emergency dispatching control method based on model predictive control |
Also Published As
Publication number | Publication date |
---|---|
CN112036092B (en) | 2023-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Sustainability analysis for Yellow River water resources using the system dynamics approach | |
CN108254032B (en) | River ultrasonic time difference method flow calculation method | |
Jalón-Rojas et al. | Turbidity in the fluvial Gironde Estuary (southwest France) based on 10-year continuous monitoring: sensitivity to hydrological conditions | |
CN112215461B (en) | Long-distance water diversion project water delivery scheduling state analysis and early warning system | |
CN104727264B (en) | A kind of flood at dam site method of reducing based on Copula function | |
CN107766985B (en) | Estuary salinity prediction method | |
CN207741811U (en) | Discharge of river real time on-line monitoring device | |
CN110186440B (en) | Dynamic monitoring method for tidal river reach channel renovation project | |
CN103235883A (en) | Non-contact river surface velocity measurement based regression calculation method for instantaneous flow rate | |
CN104652347A (en) | Method for evaluating relation between non-static water level and population affected by submerging in mountain region | |
CN109783934A (en) | A kind of mean velocity in section fitting rating method based on H-ADCP | |
Shi et al. | Fluctuations in the tidal limit of the Yangtze River estuary in the last decade | |
CN115422850B (en) | Method for judging pollution contribution rate of river inflow river drain outlet based on EFDC and edge calculation | |
CN108612037A (en) | A kind of method and its system determining river bed reference elevation based on big cross section measurement data | |
CN112036092A (en) | River flow prediction method based on relation between river intermediate surface speed and river width | |
CN107423544A (en) | A kind of determination method and device of the non-stationary influence index of River Basin Hydrology extreme value | |
Bao et al. | Dynamic correction of roughness in the hydrodynamic model | |
CN114819322B (en) | Forecasting method for flow of lake entering lake | |
CN111428915A (en) | Method for predicting silt returning amount of deep-water channel at river estuary based on big data | |
CN114001710B (en) | High-precision channel flow cross-section area measuring device and measuring method | |
CN106320257A (en) | Lake and reservoir channel storage curve determining method based on hydrometry | |
CN106320255A (en) | Overbank flow floodplain and channel excess flow computing method | |
CN102810183A (en) | Water resource energy efficiency extension ecological evaluation method based on improved distance | |
Li et al. | On the probabilistic distribution of wind speeds: theoretical development and comparison with data | |
CN114595596B (en) | Simulation method of dense sedimentation rate of reservoir sediment |
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 |