CN113390471A - River flow estimation method based on GNSS reflected signals - Google Patents
River flow estimation method based on GNSS reflected signals Download PDFInfo
- Publication number
- CN113390471A CN113390471A CN202110650289.1A CN202110650289A CN113390471A CN 113390471 A CN113390471 A CN 113390471A CN 202110650289 A CN202110650289 A CN 202110650289A CN 113390471 A CN113390471 A CN 113390471A
- Authority
- CN
- China
- Prior art keywords
- gridding
- river
- flow velocity
- section
- water level
- 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 30
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 69
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 13
- 238000001514 detection method Methods 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/002—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow wherein the flow is in an open channel
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention relates to a river flow estimation method based on GNSS reflected signals, which comprises the following steps: acquiring a water level sequence of a river section; performing Lagrange interpolation according to the water level sequence to determine a gridding water level; acquiring riverbed elevation information; performing Lagrange interpolation according to the riverbed elevation information to determine the gridded riverbed elevation; determining gridding water depth according to the gridding water level and the gridding riverbed elevation; acquiring a flow velocity value of a river section; performing Lagrange interpolation according to the flow velocity value to determine a gridding surface flow velocity; determining the average flow velocity of the gridding section according to the gridding surface flow velocity; and determining the section flow of the river according to the average flow velocity of the gridding section and the gridding water depth. The river flow is determined through the cross-section flow velocity and the water level information, and the accuracy of river flow estimation is improved.
Description
Technical Field
The invention relates to the field of river flow estimation, in particular to a river flow estimation method based on GNSS reflected signals.
Background
The river flow measuring method mainly comprises a direct measuring method and a remote sensing method. The direct measurement method adopts a measurement method of direct contact with rivers, mainly comprises a hydrological station based on a flow meter and an acoustic Doppler flow profiler, but is influenced by economic factors and complex environments, the deployment difficulty is large, and the detection cost is high. The traditional remote sensing method is mainly characterized in that an empirical model is established by combining parameters such as river width, water level or flow velocity and the like detected by a microwave remote sensor with actual river flow. According to the method, river cross-section flow velocity distribution and river bed elevation information are not needed, the detection difficulty is reduced, the ubiquitous capability is poor, a specific empirical model needs to be established according to a specific observation point, the space-time resolution is low, the size and the power consumption are large, and river flow information cannot be obtained through a single sensor.
Disclosure of Invention
The invention aims to provide a river flow estimation method based on GNSS reflected signals, which determines river flow through cross-sectional flow velocity and water level information and improves river flow estimation accuracy.
In order to achieve the purpose, the invention provides the following scheme:
a river flow estimation method based on GNSS reflected signals comprises the following steps:
acquiring a water level sequence of a river section;
performing Lagrange interpolation according to the water level sequence to determine a gridding water level;
acquiring riverbed elevation information;
performing Lagrange interpolation according to the riverbed elevation information to determine the gridded riverbed elevation;
determining gridding water depth according to the gridding water level and the gridding riverbed elevation;
acquiring a flow velocity value of a river section;
performing Lagrange interpolation according to the flow velocity value to determine a gridding surface flow velocity;
determining the average flow velocity of the gridding section according to the gridding surface flow velocity;
and determining the section flow of the river according to the average flow velocity of the gridding section and the gridding water depth.
Optionally, before the acquiring the water level sequence of the river section, the method further includes:
and carrying out meshing treatment on the river section.
Optionally, the gridding water level calculation formula is as follows:
wherein HriFor gridding the water level, BiIs the grid position, xkIs the position of the kth measurement point, xjIs the position of the jth measurement point, yjIs jth water level information.
Optionally, the calculation formula of the gridding riverbed elevation is as follows:
wherein HbiFor gridding the bed elevation, zjThe bed height of the jth measurement point location.
Optionally, the calculation formula of the gridding water depth is as follows:
hri=Hri-Hbi
wherein h isriTo grid the water depth.
Optionally, the calculation formula of the gridding surface flow velocity is as follows:
wherein v isiFor gridding the surface flow velocity, vjThe flow rate corresponding to the jth measurement point.
Optionally, the calculation formula of the average flow velocity of the gridding section is as follows:
wherein,is the average flow velocity of the gridded cross section, m is the first riverbed roughness related parameter, ksIs the second and riverbed roughness related parameter.
Optionally, the calculation formula of the cross-sectional flow of the river is as follows:
wherein Q isrThe cross-sectional flow of the river, BupAn upper bound of river width, BdownThe lower bound of the river width.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides a river flow estimation method based on GNSS reflected signals, which comprises the steps of carrying out Lagrange interpolation according to a water level sequence and determining a gridding water level; performing Lagrange interpolation according to the riverbed elevation information to determine the gridded riverbed elevation; determining the gridding water depth according to the gridding water level and the gridding riverbed elevation; performing Lagrange interpolation according to the flow velocity value to determine the gridding surface flow velocity; determining the average flow velocity of the gridding section according to the gridding surface flow velocity; and determining the section flow of the river according to the average flow velocity of the gridding section and the gridding water depth. The river flow is determined through the cross-section flow velocity and the water level information, and the accuracy of river flow estimation is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a method for estimating river discharge based on GNSS reflected signals according to the present invention;
FIG. 2 is a top view of a river;
FIG. 3 is a cross-sectional view of a river;
fig. 4 is a specific working flow chart of the method for estimating river discharge based on GNSS reflected signals according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a river flow estimation method based on GNSS reflected signals, which determines river flow through cross-sectional flow velocity and water level information and improves river flow estimation accuracy.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, the method for estimating river discharge based on GNSS reflected signals according to the present invention includes:
step 101: and acquiring a water level sequence of the river section.
Step 102: and performing Lagrange interpolation according to the water level sequence to determine the gridding water level. Wherein, the gridding water level calculation formula is as follows:
wherein HriIs the gridded water level, i.e. the distance from the river surface to the bottom of the section, BiIs the grid position, xkIs the position of the kth measurement point, xjIs the position of the jth measurement point, yjIs jth water level information.
Step 103: and acquiring river bed elevation information.
Step 104: and performing Lagrange interpolation according to the riverbed elevation information to determine the gridded riverbed elevation. Wherein, the calculation formula of the gridding riverbed elevation is as follows:
wherein HbiFor gridding the bed elevation, i.e. the bed elevation corresponding to grid i, zjThe bed height of the jth measurement point location.
Step 105: and determining the gridding water depth according to the gridding water level and the gridding riverbed elevation.
Wherein, the calculation formula of the gridding water depth is as follows:
hri=Hri-Hbi
wherein h isriTo grid the water depth.
Step 106: and acquiring the flow velocity value of the river section.
Step 107: and performing Lagrange interpolation according to the flow velocity value to determine the gridding surface flow velocity. Wherein, the calculation formula of the gridding surface flow velocity is as follows:
wherein v isiFor gridding the surface flow velocity, vjThe flow rate corresponding to the jth measurement point.
Step 108: and determining the average flow velocity of the gridding section according to the gridding surface flow velocity. Wherein, the calculation formula of the average flow velocity of the gridding section is as follows:
wherein,for gridded cross-sectional average flow velocity, m isFirst bed roughness related parameter, ksIs the second and riverbed roughness related parameter.
Step 109: and determining the section flow of the river according to the average flow velocity of the gridding section and the gridding water depth. Wherein, the calculation formula of the cross-section flow of the river is as follows:
wherein Q isrThe cross-sectional flow of the river, BupAn upper bound of river width, BdownThe lower bound of the river width.
In practical applications, step 101 further includes: and carrying out meshing treatment on the river section.
As shown in fig. 4, the present invention further provides a specific workflow of a method for estimating river discharge based on GNSS reflected signals, the steps of which are as follows:
in order to calculate the river flow, the river section is gridded. As shown in FIG. 2, assuming the grid width is Δ B, the river is divided into n grids, denoted as { B }i|i=1,2,…,n}。
Wherein, BupAnd BdownRespectively representing the upper and lower bounds of the river width,the upper rounding operator is represented. River flow was calculated every 2 hours. During the period, the mirror reflection point tracks of M satellites are distributed in the land/river or river, and N flow velocity values and K water level values are obtained in total after time-space domain mapping.
The first step is as follows: and calculating the depth of the section water.
1) Reading a water level sequence H with the length of K acquired within two hoursr(Bs) In which B issIndicating the water level HrCorresponding position. Water level HrAbout the position of the measuring point BsIs calculated as | Bs-BiI, sorting according to size, selecting 4 absolute values with minimum value, namely distance grid BiThe nearest 4 points, the position of the corresponding measuring point can be recorded as x1,x2,x3,x4,x1,x2,x3,x4At the measuring point position BsIn (1), its corresponding water level is recorded as y1,y2,y3,y4Calculating gridding water level { H using Lagrange interpolation ri1,2, …, n }, wherein:
Hrithe distance from the surface of the river to the bottom of the cross section is shown.
2) Reading riverbed elevation information Hb(Bt) (local topographic mapping or low-frequency geological radar is obtained once), the riverbed elevation information is the position of a riverbed elevation measuring point, the riverbed elevation is also a sequence related to the position of the measuring point, as shown in figure 3, the read riverbed elevation information is divided into grids, similarly to a method for calculating the water level, for each grid, 4 measuring points closest to the center of the grid are taken, and the positions of the measuring points are marked as x1,x2,x3,x4And its corresponding riverbed height is recorded as z1,z2,z3,z4And performing Lagrange interpolation on the n-th order to obtain n gridded { HbiI |, 1,2, …, n }. Wherein the elevation of the riverbed corresponding to the grid i is Hbi. Calculating gridding riverbed elevation { Hbi|i=1,2,…,n}:
3) Real water depth h corresponding to grid iriFor the depth H of the corresponding sectionriMinus the height H of the bedbi. From this, the gridded water depth h is calculatedri:
hri=Hri-Hbi
The second step is that: calculating the average flow velocity of the cross section
1) Reading N flow velocity values v (B) obtained within two hourss). Selecting distance grid B by calculating gridding water leveliThe position of the nearest 4 flow velocity values corresponding to the measuring points is recorded as x1,x2,x3,x4The corresponding flow velocity is denoted v1,v2,v3,v4Performing Lagrange interpolation to calculate gridding surface flow velocity { vi|i=1,2,…,n}:
Wherein m and ksIs a parameter related to the roughness of the riverbed, and m is 1/6, ks0.2. (m and k)sThe values of (A) are common to all people, belong to the common general knowledge in the field and need not give any explanation to the places)
The third step: calculating the cross-sectional flow
The cross-sectional flow of river water is the integral of the average flow velocity of the gridding cross-section on the cross-section:
wherein Q isrThe estimated river section flow is obtained. The first part of the formula is a theoretical expression. The theoretical expressions cannot be calculated, so discretization is performed. The second step is to adjust the water levelDiscretization is carried out, and the third step is to discretize the river width.
The method combines the river width, water level and flow rate obtained by GNSS reflected signal measurement with the actual flow, and establishes an empirical model for river flow estimation through a calculation method of river section average flow rate and water depth distribution through spatial interpolation and full parameter integration based on section flow rate and water level information, thereby realizing accurate estimation of river flow.
The invention has the advantages that:
the invention provides a full-parameter integral river flow estimation model by a method for calculating the average flow velocity and the water depth distribution of a river section through spatial interpolation and based on the section flow velocity and water level information. The GNSS reflected signal technology is used for river flow detection, algorithm complexity is low, and detection precision is high.
The GNSS reflected signal technology used by the invention has the advantages of abundant signal resources, global full-time and all-weather coverage, wide detection range and low cost.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (8)
1. A river flow estimation method based on GNSS reflected signals is characterized by comprising the following steps:
acquiring a water level sequence of a river section;
performing Lagrange interpolation according to the water level sequence to determine a gridding water level;
acquiring riverbed elevation information;
performing Lagrange interpolation according to the riverbed elevation information to determine the gridded riverbed elevation;
determining gridding water depth according to the gridding water level and the gridding riverbed elevation;
acquiring a flow velocity value of a river section;
performing Lagrange interpolation according to the flow velocity value to determine a gridding surface flow velocity;
determining the average flow velocity of the gridding section according to the gridding surface flow velocity;
and determining the section flow of the river according to the average flow velocity of the gridding section and the gridding water depth.
2. The method as claimed in claim 1, wherein the step of obtaining the water level sequence of the river section further comprises:
and carrying out meshing treatment on the river section.
5. The GNSS reflected signal based river discharge estimation method of claim 4, wherein the gridded water depth is calculated according to the following formula:
hri=Hri-Hbi
wherein h isriTo grid the water depth.
7. The method as claimed in claim 6, wherein the calculation formula of the average flow velocity of the gridding section is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110650289.1A CN113390471B (en) | 2021-06-10 | 2021-06-10 | River flow estimation method based on GNSS reflected signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110650289.1A CN113390471B (en) | 2021-06-10 | 2021-06-10 | River flow estimation method based on GNSS reflected signals |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113390471A true CN113390471A (en) | 2021-09-14 |
CN113390471B CN113390471B (en) | 2022-03-01 |
Family
ID=77620366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110650289.1A Active CN113390471B (en) | 2021-06-10 | 2021-06-10 | River flow estimation method based on GNSS reflected signals |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113390471B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116222676A (en) * | 2023-05-08 | 2023-06-06 | 成都赐华科技有限公司 | Millimeter wave water flow monitoring method and system with accurate positioning |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569988A (en) * | 2015-01-30 | 2015-04-29 | 长江水利委员会水文局 | Echo sounding-based correction method for great sounding |
CN104865581A (en) * | 2015-03-19 | 2015-08-26 | 中国科学院空间科学与应用研究中心 | GNSS-R technology-based river flow velocity measurement method and system |
CN109541162A (en) * | 2018-11-21 | 2019-03-29 | 大连融利达新能源汽车有限公司 | A kind of hydrological monitoring device based on camera and GNSS |
CN110779498A (en) * | 2019-09-19 | 2020-02-11 | 中国科学院测量与地球物理研究所 | Shallow river water depth mapping method and system based on unmanned aerial vehicle multi-viewpoint photography |
CN112415219A (en) * | 2020-11-09 | 2021-02-26 | 长江水利委员会水文局 | Positioning and correcting method for fixed speed measurement vertical line for river non-contact side-scan flow velocity measurement |
CN212808332U (en) * | 2020-09-04 | 2021-03-26 | 重庆长毛象科技有限公司 | GNSS intelligence velocity of flow measurement system |
CN112765909A (en) * | 2021-01-22 | 2021-05-07 | 武汉大学 | Method for unifying grid unit numbers of elevation data and water depth data |
WO2021088479A1 (en) * | 2019-11-05 | 2021-05-14 | 诺领科技(南京)有限公司 | A-gnss auxiliary data request method in internet of things |
-
2021
- 2021-06-10 CN CN202110650289.1A patent/CN113390471B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569988A (en) * | 2015-01-30 | 2015-04-29 | 长江水利委员会水文局 | Echo sounding-based correction method for great sounding |
CN104865581A (en) * | 2015-03-19 | 2015-08-26 | 中国科学院空间科学与应用研究中心 | GNSS-R technology-based river flow velocity measurement method and system |
CN109541162A (en) * | 2018-11-21 | 2019-03-29 | 大连融利达新能源汽车有限公司 | A kind of hydrological monitoring device based on camera and GNSS |
CN110779498A (en) * | 2019-09-19 | 2020-02-11 | 中国科学院测量与地球物理研究所 | Shallow river water depth mapping method and system based on unmanned aerial vehicle multi-viewpoint photography |
WO2021088479A1 (en) * | 2019-11-05 | 2021-05-14 | 诺领科技(南京)有限公司 | A-gnss auxiliary data request method in internet of things |
CN212808332U (en) * | 2020-09-04 | 2021-03-26 | 重庆长毛象科技有限公司 | GNSS intelligence velocity of flow measurement system |
CN112415219A (en) * | 2020-11-09 | 2021-02-26 | 长江水利委员会水文局 | Positioning and correcting method for fixed speed measurement vertical line for river non-contact side-scan flow velocity measurement |
CN112765909A (en) * | 2021-01-22 | 2021-05-07 | 武汉大学 | Method for unifying grid unit numbers of elevation data and water depth data |
Non-Patent Citations (1)
Title |
---|
朱新元等: "天然复式河道的流量计算", 《四川大学学报(工程科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116222676A (en) * | 2023-05-08 | 2023-06-06 | 成都赐华科技有限公司 | Millimeter wave water flow monitoring method and system with accurate positioning |
CN116222676B (en) * | 2023-05-08 | 2023-07-28 | 成都赐华科技有限公司 | Millimeter wave water flow monitoring method and system with accurate positioning |
Also Published As
Publication number | Publication date |
---|---|
CN113390471B (en) | 2022-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mason et al. | Construction of an inter‐tidal digital elevation model by the ‘Water‐Line’Method | |
CN109060056A (en) | A kind of river cross-section method of calculating flux of contactless radar flow measurement | |
CN108469620B (en) | Underwater terrain measurement method suitable for shallow water sea area of radiation sand ridge group | |
CN103389077B (en) | Seabed sand wave geomorphology movement detection method based on MBES (multi-beam echo sounding) | |
Dinehart | Bedform movement recorded by sequential single-beam surveys in tidal rivers | |
CN109031301A (en) | Alpine terrain deformation extracting method based on PSInSAR technology | |
Rutten et al. | Accuracy of nearshore bathymetry inverted from ${x} $-band radar and optical video data | |
CN109543356A (en) | Consider the ocean interior temperature-salinity structure remote sensing inversion method of Space atmosphere | |
CN104652347A (en) | Method for evaluating relation between non-static water level and population affected by submerging in mountain region | |
CN103235883A (en) | Non-contact river surface velocity measurement based regression calculation method for instantaneous flow rate | |
CN110208807A (en) | A kind of raininess grade inversion method based on marine radar image detection region otherness parameter | |
Moulton et al. | Improving the time resolution of surfzone bathymetry using in situ altimeters | |
CN113390471B (en) | River flow estimation method based on GNSS reflected signals | |
CN111220146B (en) | Underwater terrain matching and positioning method based on Gaussian process regression learning | |
Yutsis et al. | Bottom topography, recent sedimentation and water volume of the Cerro Prieto Dam, NE Mexico | |
Matte et al. | Quantifying lateral and intratidal variability in water level and velocity in a tide‐dominated river using combined RTK GPS and ADCP measurements | |
CN117516636B (en) | Coastal dyke safety monitoring and early warning method and system | |
Jahanmard et al. | Towards realistic dynamic topography from coast to offshore by incorporating hydrodynamic and geoid models | |
CN115808686A (en) | Vertical line deviation gridding method and system based on waveform retracing quality weighting | |
Swales | Geostatistical estimation of short-term changes in beach morphology and sand budget | |
Schmitt et al. | Use of swath bathymetry in the investigation of sand dune geometry and migration around a near shore ‘banner’tidal sandbank | |
Varbla et al. | Utilizing airborne laser scanning and geoid model for near-coast improvements in sea surface height and marine dynamics | |
Maxwell et al. | Generating river bottom profiles with a dual-frequency identification sonar (DIDSON) | |
Byrnes et al. | Quantifying potential measurement errors and uncertainties associated with bathymetric change analysis | |
Travers | Quantifying sea-ice volume flux using moored instrumentation in the Bering Strait |
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 |