CN115326026A - Method and device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation - Google Patents

Method and device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation Download PDF

Info

Publication number
CN115326026A
CN115326026A CN202210824755.8A CN202210824755A CN115326026A CN 115326026 A CN115326026 A CN 115326026A CN 202210824755 A CN202210824755 A CN 202210824755A CN 115326026 A CN115326026 A CN 115326026A
Authority
CN
China
Prior art keywords
hydrodynamic
assimilation
calculation
contact measurement
model
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.)
Pending
Application number
CN202210824755.8A
Other languages
Chinese (zh)
Inventor
黄凯霖
陈华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Publication of CN115326026A publication Critical patent/CN115326026A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • 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

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Hydrology & Water Resources (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Fluid Mechanics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measuring Volume Flow (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

The invention provides a method and a device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation, which can acquire the hydraulic characteristics of a target monitoring station in a high-efficiency and stable manner in real time. The method for acquiring the hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation is characterized by comprising the following steps of: step I, hydrologic monitoring is carried out on the research area based on a non-contact measurement mode, hydrologic monitoring data are obtained, and a surface flow velocity measurement value of the river cross section of the research area is obtained based on the hydrologic monitoring data; step II, constructing a hydrodynamic model based on hydrologic monitoring data of the research area; step III, hydraulic characteristic assimilation calculation: and calculating to obtain the optimal dynamic parameters of the river hydrodynamic system and the optimal assimilation flow velocity and flow corresponding to the flow velocity measurement value based on the hydrodynamic model according to the surface flow velocity measurement value obtained in a non-contact measurement mode.

Description

Method and device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation
Technical Field
The invention belongs to the technical field of hydrological monitoring, and particularly relates to a method and a device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation.
Background
Hydrological monitoring is an important national foundation. The hydrologic monitoring and collecting object comprises basic data such as water level, flow velocity and flow. Real-time hydrologic monitoring is an important technical support for realizing hydrologic modernization construction. The traditional hydrology test mainly takes contact measurement as a main means, and hydrology test facilities suitable for national standard specifications comprise a rotor current meter, an acoustic Doppler current profiler and the like, so that the problems of high labor and equipment cost, potential safety hazards in the flood period and the like exist.
With the development of modern information technology, radar and video gradually receive more and more attention in the field of hydrological non-contact measurement. For the radar, the civil engineering cost for installing the radar is higher, and different processing modes are required for monitoring the flow speed with different measuring ranges. For videos, the speed of acquiring monitoring data is high, the advantages of low installation cost are achieved under the capability of providing historical monitoring scenes which can be duplicated, and the flow velocity of a target area can be obtained only after the target area is installed and fixed on the bank side of a river channel.
The current main flow calculation method of the non-contact equipment is a flow velocity-area method, and the specific method is as shown in formula (1), namely, the surface flow velocity v of each calculation area is calculated by dividing the section into n calculation areas i And by converting the coefficient alpha i Obtaining the depth average flow velocity, and finally, the depth average flow velocity and the section area A i And multiplying and summing to obtain the flow of the section.
Figure BDA0003743584690000011
Although the non-contact method can acquire the surface flow velocity of the river channel at a lower cost and a higher speed, stable full-section flow velocity monitoring is still more difficult due to the change of an external light source and the existence of external environmental factors such as abnormal noise, and a challenge is also formed for continuous observation of flow.
Disclosure of Invention
The present invention has been made to solve the above-described problems, and an object of the present invention is to provide a method and an apparatus for acquiring a hydraulic characteristic based on non-contact measurement and hydrodynamic fusion assimilation, which can acquire a hydraulic characteristic of a target monitoring station in real time in an efficient and stable manner.
In order to achieve the purpose, the invention adopts the following scheme:
< method >
The invention provides a method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation, which is characterized by comprising the following steps of:
step I, hydrologic monitoring is carried out on the research area based on a non-contact measurement mode, hydrologic monitoring data are obtained, and a surface flow velocity measurement value of the river cross section of the research area is obtained based on the hydrologic monitoring data;
step II, constructing a hydrodynamic model based on hydrologic monitoring data of the research area;
step III, hydraulic characteristic assimilation calculation:
and calculating to obtain the optimal dynamic parameters of the river hydrodynamic system and the optimal assimilation flow speed and flow corresponding to the flow speed measurement value based on the hydrodynamic model according to the surface flow speed measurement value obtained in a non-contact measurement mode. The optimal assimilation flow rate means a corrected optimal value of the flow rate (excluding adverse influence factors such as external environment) that is the value closest to the true value. The same principle is applied to optimal assimilation flow.
Preferably, in the method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention, step I may further include the following steps:
step 1, flow measurement field pretreatment: installing and arranging equipment for acquiring monitoring data in a non-contact mode, and positioning coordinates;
step 2, analyzing the cross section flow velocity: and analyzing and calculating the acquired video to obtain a surface flow velocity measurement value of the river channel section based on the selected measurement scheme.
Preferably, the method for acquiring the hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: in step I, for the on-site space positioning, firstly, a corresponding marking plate is arranged at the near bank of the river channel, and the European space coordinates (x, y, z) of the marker are set by using the related positioning equipment T Positioning and recording; for equipment installation, after processing infrastructures such as power supply, moisture protection and network, a camera is started to record images of a monitored target, and image pixel coordinates [ i, j ] corresponding to a marker in the previous step are recorded] T 。
Preferably, the method for acquiring the hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: in step II, grids are divided according to the collected hydrological data of the research area, a model is preliminarily built through interpolation, boundary conditions are determined, model calibration is carried out, and a hydrodynamic model after calibration is obtained.
Preferably, the method for acquiring the hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: in step II, the model construction unit performs model establishment, boundary condition determination, and model calibration; for the part of model establishment, meshing and interpolating are carried out by relying on the terrain data collected in the previous step; in the boundary condition determining section, the upstream and downstream boundary positions and the types of boundary conditions of the entrance and the exit are determined based on the established computational mesh. Finally, the model is calibrated according to the collected observation (measurement) data.
Preferably, the method for acquiring the hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: in step II, the boundary of the modeled area is determined based on the topographic data collected in step I, and the calculation area Σ is decomposed into finite grid sets Σ that do not overlap each other i ＝{(V i ,E i ) In which V is i Set of all spatially discrete points of the ith mesh, E i The set of edges formed for all the discrete points of the ith mesh satisfies
Figure BDA0003743584690000031
And carrying out interpolation on the area grids:
determination of boundary conditions:
Figure BDA0003743584690000032
using the above control equation, where ρ represents the density of the fluid, t represents the time variable,
Figure BDA0003743584690000033
is indicative of the velocity vector of the fluid,
Figure BDA0003743584690000034
representing stress tensor, R representing Reynolds stress, and g representing gravitational acceleration;
calculating the model by respectively specifying corresponding boundary conditions (including but not limited to the combination of a water level sequence and a flow rate sequence) at upstream and downstream;
and (3) model calibration: and calibrating the model according to the collected historical hydrological data.
Preferably, the method for acquiring the hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: the non-contact measurement mode is a video measurement mode or a radar measurement mode.
Preferably, the method for acquiring the hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: in the step III, the coupling relation between the monitoring state and the physical numerical model is established by taking the information of the fusion monitoring perception source as a target, and the coupling detection method which is most suitable for the measured variable and the physical numerical model is searched. Taking a fusion method based on bayesian theory as an example, similar approaches can also adopt a variation Assimilation method (Variational analysis) such as Ensemble Kalman (Ensemble Kalman filter), and a hydrodynamic equation of motion as a basis for calculation according to the surface flow velocity of the river channel observed (measured) by video or radar or other non-contact methods. And finally, analyzing the posterior distribution of the parameters to obtain an uncertain interval of the monitored flow, thereby obtaining the distribution condition of the flow and the flow velocity on the cross section.
Preferably, the method for acquiring the hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: when the non-contact measurement mode is a video measurement mode and a fusion method based on Bayesian theory is adopted, the step III comprises the following steps:
step III-1, sampling: setting the parameter of the control equation established in the step II as theta, and setting the corresponding prior distribution as p (theta) opt ) Show that first from p (θ) opt ) Middle sample, the sampled sample being denoted by theta * ；
Step III-2, calculation: according to the sampled parameter theta * Calculating model, and recording characteristic variable of the calculating model as S m ；
Step III-3, judging, calculating a posterior parameter distribution p (theta | d (S)) based on the following calculation index function according to the calculation result of the step III-2 m ,S obs ) ); after the posterior parameter distribution is obtained by calculation, the following formula can be further adopted to judge whether the updating is needed or not:
d(S m ,S obs )≤ε；
wherein ε is an acceptable threshold;
III-4, analyzing and evaluating uncertainty of flow calculation results, determining posterior distribution of dynamic parameters of the river channel hydrodynamic system, and obtaining optimal hydrodynamic system parameters theta according to maximum posterior opt And determining an optimal assimilation flow rate and flux measurements corresponding to the flow rate measurements; at the same time, the confidence interval of the flow can be analyzed according to the calculated posterior distribution.
< means >
Further, the invention also provides a device for acquiring river hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation, which is characterized by comprising:
the acquisition part is used for hydrologically monitoring the research area based on a non-contact measurement mode, acquiring hydrologically monitoring data and obtaining a surface flow velocity measurement value of the river section of the research area based on the hydrologically monitoring data;
a model construction unit that constructs a hydrodynamic model based on hydrologic monitoring data of a study area;
a calculation unit for hydraulic characteristic assimilation calculation: calculating to obtain an optimal power parameter of the river channel hydrodynamic system and an optimal assimilation flow velocity and flow measurement value corresponding to flow velocity measurement data based on a hydrodynamic model according to a surface flow velocity measurement value obtained in a non-contact measurement mode; and
and the control part is in communication connection with the acquisition part, the model construction part and the calculation part and controls the operation of the acquisition part, the model construction part and the calculation part.
Preferably, the device for acquiring the river hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention further comprises: and the input display part is in communication connection with the acquisition part, the model construction part, the calculation part and the control part and is used for allowing a user to input an operation instruction and performing corresponding display.
Preferably, the device for acquiring the river hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: the acquisition part comprises a pretreatment unit and a flow velocity analysis unit; the pretreatment unit is used for recording relevant information of non-contact flow measurement on-site pretreatment; the flow velocity analysis unit obtains hydrological monitoring data of hydrological monitoring of a research area based on a non-contact measurement mode, and the hydrological monitoring data are analyzed and calculated to obtain a surface flow velocity measurement value of a river section of the research area.
Preferably, the device for acquiring the river hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: in the model construction part, the boundary of the modeling area is determined based on the topographic data collected by the acquisition part, and the calculation area sigma is decomposed into non-overlapping finite grid sets sigma i ＝{(V i ,E i ) In which V is i All spatially discrete point sets of the ith grid, E i The set of edges formed for all the discrete points of the ith mesh satisfy
Figure BDA0003743584690000051
And interpolating the area grid:
determination of boundary conditions:
Figure BDA0003743584690000052
using the above control equation, where ρ represents the density of the fluid, t represents the time variable,
Figure BDA0003743584690000053
which is indicative of the velocity vector of the fluid,
Figure BDA0003743584690000054
representing stress tensor, R representing Reynolds stress, and g representing gravitational acceleration;
respectively appointing corresponding boundary conditions (a water level sequence and a flow sequence) at upstream and downstream to calculate the model;
and (3) model calibration: and calibrating the model according to the collected historical hydrological data.
Preferably, the device for acquiring the river hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention can also have the following characteristics: when the non-contact measurement mode is a video measurement mode and a fusion method based on Bayesian theory is adopted, the calculation part adopts the following steps to carry out hydraulic characteristic assimilation calculation:
step III-1, sampling: setting the parameter of the control equation established in the step II as q, and setting the corresponding prior distribution as p (theta) opt ) Show that first from p (θ) opt ) Middle sample, the sampled sample being denoted by theta * ；
Step III-2, calculation: according to the sampled parameter theta * Calculating model, and recording characteristic variable of the calculating model as S m ；
Step III-3, judging: from the calculation result of step III-2, the posterior parameter distribution p (θ | d (S) is calculated based on the following calculation index function m ,S obs ) ); after the posterior parameter distribution is obtained by calculation, the following formula can be further adopted to judge whether the updating is needed or not:
d(S m ,S obs )≤ε；
wherein ε is an acceptable threshold;
step III-4, analyzing and evaluating uncertainty of flow calculation results, determining posterior distribution of dynamic parameters of the river channel hydrodynamic system, and obtaining optimal hydrodynamic system parameters theta according to maximum posterior opt And determining an optimal assimilation flow rate and flux corresponding to the flow rate measurements; at the same time, the confidence interval of the flow can be analyzed according to the calculated posterior distribution.
The method and the device for acquiring the hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation provided by the invention are used for acquiring surface flow velocity data based on a non-contact measurement means aiming at the field of hydrological testing of a river channel, fusing the surface flow velocity with a hydrodynamic model, taking the fused hydrological testing result as a final measuring result, namely obtaining the measuring result which is closer to a real measuring reference value.
Drawings
Fig. 1 is a flowchart of a method for acquiring a hydraulic characteristic based on non-contact measurement-hydrodynamic fusion assimilation according to an embodiment of the present invention.
FIG. 2 is a schematic diagram of on-site video stream processing according to an embodiment of the present invention;
fig. 3 is a flowchart of a hydraulic characteristic assimilation calculation method according to an embodiment of the present invention;
fig. 4 is a comparison diagram related to the embodiment of the present invention, in which curve 1 represents the monitored flow rate corresponding to the fused optimal parameter, and curve 2 represents the video-measured river flow rate.
Detailed Description
The following describes in detail specific embodiments of the method and apparatus for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation according to the present invention with reference to the accompanying drawings.
As shown in fig. 1, the method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation provided by this embodiment includes 1) video flow measurement field preprocessing; 2) Monitoring video flow velocity analysis; 3) Constructing a field hydrodynamic model; 4) Bayes flow assimilation calculation; in four respects, each of which is described below.
Step 1, as shown in fig. 2, performing video stream measurement on-site preprocessing:
step 1-1, space positioning and equipment installation: the selected monitoring area should have good lighting conditions while avoiding environments that are prone to damage to the equipment. And configuring infrastructures such as power supply and transmission required by video recording and transmission. The infrastructure comprises a high-definition camera, a power supply, a fixed carrier and the like; a bridging device for transmitting data; the remote PC equipment and be used for the sign and the measuring instrument of demarcation. Before surveying, a measuring instrument such as a total station is used for carrying out relevant calibration positioning processing, and the setting of a calibration object should meet the requirements of high contrast background and easy searching. The specific positioning method is to note l n Determining the world coordinates [ x, y, z ] of the marker for the parameter to be determined] T To pixel coordinate u, v] T Projection relation proj ([ u, v) ]] T )→[x,y,z] T As shown in equation (1). The determination method is that at least 6 calibration plates are arranged, and the corresponding coefficient l is obtained by solving through a simultaneous linear equation set n 。
Figure BDA0003743584690000071
Step 1-2, topographic exploration: if topographic data of a monitored site does not exist in advance, topographic data of the upstream and the downstream of the monitored site needs to be measured, and if the topographic data refers to local reorganization basic data.
Step 2, analyzing the flow rate of the section video:
and 2-1, selecting a measurement scheme, and selecting a proper measurement scheme aiming at the measured scene. The invention is not limited to algorithms and devices for estimating flow velocity from video, as long as the method gives a measure of the surface flow velocity along a section. Take a flow field estimation method based on a variational principle and a space-time image method as an example. And (3) recording the video on the hardware basis completed in the step (1) after the selection is completed, and preparing the video before calculation.
Step 2-2, calculating the flow rate, and dividing the flow rate into two parts:
and 2-2a, aiming at a river channel surface flow velocity calculation method based on a variational principle, estimating the flow velocity by minimizing an energy functional shown in (2). In the formula, I represents a pixel set of an ith row in an image, j represents a pixel set of a jth column in the image, t is an image shooting time variable, w (I, j) represents a pixel flow field to be solved, I = I (I, j, t) represents a shot image sequence, and D is a diffusion coefficient;
Figure BDA0003743584690000072
is a gradient operator, delta is a Laplace operator, and alpha is a weight coefficient; phi, phi are functions for flow and gradient, respectively. And projecting the calculated flow velocity field to a world coordinate system to obtain the corresponding flow velocity.
Figure BDA0003743584690000073
Step 2-2b, aiming at the river surface flow velocity calculation method of the space-time image, the technical scheme adopted by the partial method is that through a flow velocity line drawn along the flow direction in an obliquely shot video, the change of the brightness of pixel points on the flow velocity line along with time is drawn on an image, the horizontal axis represents the length of the flow velocity line, the vertical axis represents a time axis, the generated image is called the space-time image (hereinafter referred to as STIV for short), and texture angles in the STIV image are marked as the space-time image (hereinafter referred to as STIV for short)
Figure BDA0003743584690000074
The relationship between the flow velocity and the texture angle is satisfied
Figure BDA0003743584690000081
Calculated maximum sequence density estimated value w m And = u is the calculated flow velocity, and assuming that in a three-dimensional space, the moving distance of the characteristic quantity moving along the speed measurement line in the time T is S, and meanwhile, the pixel point moves by i pixels in the k frame, the flow velocity of the speed measurement line is as shown in (3). Wherein fps is the frame rate of the camera and the unit is frame/s,
Figure BDA0003743584690000082
is the tangent of the slope, P, of the STIV spatio-temporal image x Is the spatial distance represented by the pixel.
Figure BDA0003743584690000083
And 3, constructing a hydrodynamic model of the research area, wherein the invention content of the part is divided into three contents of model establishment, boundary condition determination and model calibration.
Step 3-1, establishing a model, determining a region boundary for modeling aiming at the topographic data collected in the step 1, wherein the grid division meets the following requirement, and decomposing a calculation region into a non-overlapping limited grid set sigma for the calculation region sigma i ＝{(V i ,E i ) In which V is i Set of all spatially discrete points of the ith mesh, E i The set of edges formed for all the discrete points of the ith mesh. Satisfy the requirements of
Figure BDA0003743584690000084
And interpolates the area grid.
Step 3-2, determining boundary conditions: a turbulence model is used as a control equation, wherein rho represents the density of the fluid, t represents a time variable,
Figure BDA0003743584690000085
representing the velocity vector of the fluid.
Figure BDA0003743584690000086
Representing the stress tensor. R represents Reynolds stress.
Figure BDA0003743584690000087
After the terrain is processed in the step 3-1, the boundary can be specified, and the model is calculated by respectively specifying corresponding boundary conditions (water level sequence and flow rate sequence) at the upstream and the downstream.
And 3-3, carrying out model calibration, and optimizing the parameters of the model according to the collected historical data.
And 4, bayes flow velocity-flow assimilation calculation, wherein a visual fusion flow velocity-flow assimilation model is established by adopting surface flow field information obtained through video calculation and adopting a hydrodynamic motion equation and a Bayes theory. The specific technical scheme is that the optimal parameter theta of the river hydrodynamic model is determined on the basis of the three-dimensional closed turbulence dynamic model established on the basis of the steps opt . The specific implementation is shown in fig. 3:
step 4-1, sampling, wherein the parameter of the turbulence model established in the step is set as theta, and the corresponding prior distribution is p (theta) opt ) Show that first from p (θ) opt ) Middle sample, the sampled sample being denoted as θ * 。
Step 4-2, calculating according to the sampled parameter theta * A calculation model, wherein the characteristic variables of the calculation model are denoted as S m 。
Step 4-3, judging, namely adopting (6) to calculate an index function to judge whether to update the posterior parameter distribution p (theta | d (S) according to the calculation result of the step 4-2 m ,S obs ))。
d(S m ,S obs )≤ε (5)
Finally, the uncertainty of the flow calculation result is analyzed and evaluated. From the above, the posterior distribution of the dynamic parameters of the river hydrodynamic system and the optimal hydrodynamic system parameters found from the maximum posterior, as well as the optimal assimilation flow observations corresponding to the flow velocity observations, can be determined. At the same time, the confidence interval of the flow can be analyzed according to the calculated posterior distribution.
As shown in fig. 4, in this embodiment, a result of on-site measurement of a certain river reach in Fujian is given, and the measured flow rate can efficiently and stably obtain hydraulic characteristics such as the river flow of a target monitoring site and the distribution of the flow rate on a section in real time through the algorithm flow provided by the present invention, and obtain a measurement result closer to a true measurement reference value.
Further, the embodiment also provides a device capable of automatically realizing the method, and the device comprises an acquisition part, a model construction part, a calculation part, an input display part and a control part.
The acquisition part carries out hydrological monitoring on the research area based on a non-contact measurement mode, acquires hydrological monitoring data and obtains a surface flow velocity measurement value of the river section of the research area based on the hydrological monitoring data. In this embodiment, the acquisition unit includes a preprocessing unit and a flow rate analysis unit. The preprocessing unit executes the content described in the step 1, performs non-contact flow measurement on-site preprocessing, and records relevant information of the non-contact flow measurement on-site preprocessing. And the flow velocity analysis unit executes the content described in the step 2, acquires hydrologic monitoring data for hydrologic monitoring of the research area based on a non-contact measurement mode, and analyzes and calculates the hydrologic monitoring data to obtain a surface flow velocity measurement value of the river cross section of the research area.
The model building unit performs the process described in step 3 above to build a hydrodynamic model based on the hydrographic monitoring data of the region of interest.
The calculation unit performs the flow velocity-flow assimilation calculation to obtain the hydraulic characteristics of the river in the study area, as described in step 4 above. The method comprises the following specific steps: and calculating to obtain an optimal dynamic parameter of the river channel hydrodynamic system and an optimal assimilation flow observation value corresponding to the flow observation data based on the hydrodynamic motion model according to the surface flow velocity measurement value obtained in the non-contact measurement mode.
And the input display part is in communication connection with the acquisition part, the model construction part and the calculation part and is used for allowing a user to input an operation instruction and performing corresponding display.
The control part is in communication connection with the acquisition part, the model construction part and the calculation part and controls the operation of the acquisition part, the model construction part and the calculation part.
The above embodiments are merely illustrative of the technical solutions of the present invention. The method and device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation are not limited to what is described in the above embodiments, but are subject to the scope defined by the claims. Any modification or supplement or equivalent replacement made by a person skilled in the art on the basis of this embodiment is within the scope of the invention as claimed in the claims.

Claims (10)

1. The method for acquiring the hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation is characterized by comprising the following steps of:
step I, hydrologic monitoring is carried out on the research area based on a non-contact measurement mode, hydrologic monitoring data are obtained, and a surface flow velocity measurement value of the river cross section of the research area is obtained based on the hydrologic monitoring data;
step II, constructing a hydrodynamic model based on hydrologic monitoring data of the research area;
step III, hydraulic characteristic assimilation calculation:
and calculating to obtain the optimal dynamic parameters of the river hydrodynamic system and the optimal assimilation flow speed and flow corresponding to the flow speed measurement value based on the hydrodynamic model according to the surface flow speed measurement value obtained in a non-contact measurement mode.
2. The method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 1, wherein the method comprises the following steps:
wherein, step I includes the following steps:
step 1, flow measurement on-site pretreatment: installing and arranging equipment for acquiring monitoring data in a non-contact mode, and positioning coordinates;
step 2, analyzing the cross section flow velocity: and analyzing and calculating the acquired video to obtain a surface flow velocity measurement value of the river channel section based on the selected measurement scheme.
3. The method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 1, wherein the method comprises the following steps:
in the step II, grids are divided according to the collected hydrological data of the research area, a model is preliminarily built through interpolation, boundary conditions are determined, model calibration is carried out, and the calibrated hydrodynamic model is obtained.
4. The method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 1, wherein the method comprises the following steps:
in step II, the boundary of the modeling area is determined based on the terrain data collected in step I, and for the calculation area sigma, the calculation area is decomposed into non-overlapping finite grid set sigma i ＝{(V i ,E i ) In which V is i All spatially discrete point sets of the ith grid, E i The set of edges formed for all the discrete points of the ith mesh satisfy
Figure FDA0003743584680000011
And interpolating the area grid:
determination of boundary conditions:
Figure FDA0003743584680000021
using the above control equation, where ρ represents the density of the fluid, t represents the time variable,
Figure FDA0003743584680000022
which is indicative of the velocity vector of the fluid,
Figure FDA0003743584680000023
showing stress tensionAmount, R represents Reynolds stress, g represents gravitational acceleration;
respectively appointing corresponding boundary conditions at upstream and downstream to calculate the model;
and (3) model calibration: and calibrating the model according to the collected historical hydrological data.
5. The method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 1, wherein the method comprises the following steps:
the non-contact measurement mode is a video measurement mode or a radar measurement mode.
6. The method for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 1, wherein the method comprises the following steps:
when the non-contact measurement mode is a video measurement mode and a fusion method based on Bayesian theory is adopted, the step III comprises the following steps:
step III-1, sampling: setting the parameter of the control equation established in the step II as theta, and setting the corresponding prior distribution as p (theta) opt ) Denotes, first from p (θ) opt ) Middle sample, the sampled sample being denoted by theta * ；
Step III-2, calculation: according to the sampled parameter theta * A calculation model, wherein the characteristic variables of the calculation model are denoted as S m ；
Step III-3, judging, calculating a posterior parameter distribution p (theta | d (S)) based on the following calculation index function according to the calculation result of the step III-2 m ,S obs ))；
III-4, analyzing and evaluating uncertainty of flow calculation results, determining posterior distribution of dynamic parameters of the river channel hydrodynamic system, and obtaining optimal hydrodynamic system parameters theta according to maximum posterior opt And determining an optimal assimilation flow rate and flux corresponding to the flow rate measurement; at the same time, the confidence interval of the flow can be analyzed according to the calculated posterior distribution.
7. Device based on non-contact measurement-hydrodynamic force fuses assimilation and acquires hydraulic power characteristic, its characterized in that includes:
the acquisition part is used for hydrologically monitoring the research area based on a non-contact measurement mode, acquiring hydrologic monitoring data and obtaining a surface flow velocity measurement value of the river cross section of the research area based on the hydrologic monitoring data;
a model construction unit for constructing a hydrodynamic model based on hydrologic monitoring data of a study area;
a calculation unit for hydraulic characteristic assimilation calculation: calculating to obtain optimal dynamic parameters of the river hydrodynamic system and optimal assimilation flow velocity and flow corresponding to the flow velocity measurement value based on a hydrodynamic motion model according to a surface flow velocity measurement value obtained in a non-contact measurement mode; and
and the control part is in communication connection with the acquisition part, the model construction part and the calculation part and controls the operation of the acquisition part, the model construction part and the calculation part.
8. The device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 7, further comprising:
and the input display part is in communication connection with the acquisition part, the model construction part, the calculation part and the control part and is used for allowing a user to input an operation instruction and performing corresponding display.
9. The device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 7, wherein:
wherein the acquisition section includes: a pretreatment unit and a flow velocity analysis unit;
the pretreatment unit is used for recording relevant information of non-contact flow measurement on-site pretreatment; the flow velocity analysis unit acquires hydrological monitoring data of hydrological monitoring of a research area based on a non-contact measurement mode, and analyzes and calculates the hydrological monitoring data to obtain a surface flow velocity measurement value of a river cross section of the research area.
10. The device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation as claimed in claim 7, wherein:
when the non-contact measurement mode is a video measurement mode and a fusion method based on the Bayesian theory is adopted, the calculation part performs hydraulic characteristic assimilation calculation by adopting the following steps:
step III-1, sampling: setting the parameter of the control equation established in the step II as theta, and setting the corresponding prior distribution as p (theta) opt ) Denotes, first from p (θ) opt ) Middle sample, the sampled sample being denoted by theta * ；
Step III-2, calculation: according to the sampled parameter theta * Calculating model, and recording characteristic variable of the calculating model as S m ；
Step III-3, judging: from the calculation result of step III-2, the posterior parameter distribution p (θ | d (S) is calculated based on the following calculation index function m ,S obs ))；
Step III-4, analyzing and evaluating uncertainty of flow calculation results, determining posterior distribution of dynamic parameters of the river channel hydrodynamic system, and obtaining optimal hydrodynamic system parameters theta according to maximum posterior opt And determining an optimal assimilation flow rate and flux corresponding to the flow rate measurements; at the same time, the confidence interval of the flow can be analyzed according to the calculated posterior distribution.
CN202210824755.8A 2022-07-05 2022-07-13 Method and device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation Pending CN115326026A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2022107918752 2022-07-05
CN202210791875 2022-07-05

Publications (1)

Publication Number Publication Date
CN115326026A true CN115326026A (en) 2022-11-11

Family

ID=83917808

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210824755.8A Pending CN115326026A (en) 2022-07-05 2022-07-13 Method and device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation

Country Status (1)

Country Link
CN (1) CN115326026A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115773797A (en) * 2023-02-10 2023-03-10 成都秦川物联网科技股份有限公司 Intelligent gas flow correction method, internet of things system, device and medium
CN116222676A (en) * 2023-05-08 2023-06-06 成都赐华科技有限公司 Millimeter wave water flow monitoring method and system with accurate positioning

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115773797A (en) * 2023-02-10 2023-03-10 成都秦川物联网科技股份有限公司 Intelligent gas flow correction method, internet of things system, device and medium
US11867548B2 (en) 2023-02-10 2024-01-09 Chengdu Qinchuan Iot Technology Co., Ltd. Methods, Internet of Things systems, and mediums for correcting smart gas flow
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

Similar Documents

Publication Publication Date Title
KR101978351B1 (en) System and Method for Measuring Real Time Water Discharge based on CCTV Image
CN111798386B (en) River channel flow velocity measurement method based on edge identification and maximum sequence density estimation
CN115326026A (en) Method and device for acquiring hydraulic characteristics based on non-contact measurement-hydrodynamic fusion assimilation
Le Coz et al. Performance of image-based velocimetry (LSPIV) applied to flash-flood discharge measurements in Mediterranean rivers
CN105045950B (en) A kind of bridge security assessment system based on 3 D laser scanning
Fujita et al. Application of aerial LSPIV to the 2002 flood of the Yodo River using a helicopter mounted high density video camera
Muste et al. Considerations on direct stream flow measurements using video imagery: Outlook and research needs
Ran et al. Application of an automated LSPIV system in a mountainous stream for continuous flood flow measurements
Tauro et al. Flow monitoring with a camera: a case study on a flood event in the Tiber river
CN111781113B (en) Dust grid positioning method and dust grid monitoring method
ITRM20070134A1 (en) METHOD OF MEASURING THE FLOW RATE IN A WATER COURSE AND ITS SYSTEM
CN113689473B (en) River channel surface flow velocity calculation method and device based on variational principle
CN117408173B (en) Hydrologic flow recompilation intelligent model construction method based on machine learning
CN114117913A (en) Dam post-dam open-air seepage monitoring system and measured value attribution analysis method
CN115422766A (en) Debris flow monitoring method and system based on digital twinning technology
CN116663762A (en) Urban planning underground space investigation and mapping method and system
CN115585731A (en) Space-air-ground integrated hydropower station space state intelligent monitoring management system and method thereof
CN116776444A (en) Digital twinning-based immersed tube tunnel construction method and system
CN106303412B (en) Refuse dump based on monitoring image is displaced remote real time monitoring method
CN104597286A (en) Yellow River ice run ice floating speed measuring method
CN114998425B (en) Target object geographic coordinate positioning method and device based on artificial intelligence
Hou et al. Adaptive large-scale particle image velocimetry method for physical model experiments of flood propagation with complex flow patterns
CN116858098B (en) Automatic acquisition method and system for multi-element information of tunnel in construction period
CN116203277B (en) Sea surface small-scale flow field measurement method based on PTV and PIV technologies
Tauro et al. Surface water velocity observations from a camera: a case study on the Tiber River

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