CN110426088B - Flow monitoring method, device and equipment and readable storage medium - Google Patents
Flow monitoring method, device and equipment and readable storage medium Download PDFInfo
- Publication number
- CN110426088B CN110426088B CN201910584822.1A CN201910584822A CN110426088B CN 110426088 B CN110426088 B CN 110426088B CN 201910584822 A CN201910584822 A CN 201910584822A CN 110426088 B CN110426088 B CN 110426088B
- Authority
- CN
- China
- Prior art keywords
- flow
- speed
- speed sensor
- monitoring device
- formula
- 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.)
- Active
Links
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/66—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
- G01F1/667—Arrangements of transducers for ultrasonic flowmeters; Circuits for operating ultrasonic flowmeters
- G01F1/668—Compensating or correcting for variations in velocity of sound
Landscapes
- Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
- Electromagnetism (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
Abstract
The invention discloses a flow monitoring method, a device, equipment and a readable storage medium, wherein the flow monitoring method comprises the following steps: acquiring the water level in the monitoring device; judging whether the water level reaches the height of the monitoring device or not; when the water level reaches the height of the monitoring device, acquiring a preset first flow formula and flow rates measured by each speed sensor; and inputting the flow speed measured by each speed sensor into a first flow formula to obtain the overflow flow of the water cross section in the monitoring device. According to the invention, a proper flow formula is obtained by judging the water level in the monitoring device, so that the calculated overflow flow of the water passing section is more accurate.
Description
Technical Field
The invention relates to the technical field of hydrological monitoring, in particular to a flow monitoring method, a flow monitoring device, flow monitoring equipment and a readable storage medium.
Background
The water measuring equipment on the market at present comprises a positive displacement flowmeter, a differential pressure flowmeter, a float flowmeter, a mass flowmeter, a vortex flowmeter, a rotary wheel type water meter, a shunt rotary wing type flowmeter, an electromagnetic flowmeter, an ultrasonic flowmeter, a turbine flowmeter, a rectangular box culvert water meter and a micro-power consumption electronic water level flowmeter. The passive type mainly refers to a listening method, calculates the fluid speed by means of noise generated by fluid flow, and mainly comprises a propagation speed time difference method, a Doppler method and a beam offset method. The ultrasonic flow measurement probes are arranged in various ways, the installation and arrangement positions of the probes of different water measuring equipment are different, and the probes are divided into single channels, double channels and multiple channels according to different ultrasonic sound channels. The method can be divided into a parallel arrangement and a symmetrical cross arrangement of double sound channels according to different sound channel installation modes, such as a Z-type method, a V-type method, an X-type method and an N-type method. However, due to the reasons of technical and theoretical research, the existing water measuring equipment has poor measurement accuracy.
Disclosure of Invention
In view of this, embodiments of the present invention provide a method, an apparatus, a device and a readable storage medium for flow monitoring, so as to solve the problem in the prior art that the flow monitoring measurement accuracy is poor.
According to a first aspect, an embodiment of the present invention provides a flow rate monitoring method applied to a monitoring device provided with a plurality of speed sensors, including the following steps:
acquiring the water level in the monitoring device;
judging whether the water level reaches the height of the monitoring device or not;
when the water level reaches the height of the monitoring device, acquiring a preset first flow formula and flow rates measured by all the speed sensors;
and inputting the flow speed measured by each speed sensor into the first flow formula to obtain the overflow flow of the water cross section in the monitoring device.
The flow monitoring method provided by the embodiment of the invention comprises the steps of obtaining the water level in a monitoring device, judging whether the water level reaches the height of the monitoring device, and obtaining a preset first flow formula and the flow velocity measured by each velocity sensor when the water level reaches the height of the monitoring device; the flow velocity measured by each velocity sensor is input into the first flow formula to obtain the overflow flow of the water passing section in the monitoring device.
With reference to the first aspect, in a first embodiment of the first aspect, the flow monitoring method further includes the steps of:
when the water level does not reach the height of the monitoring device, determining the submerged speed sensor in the monitoring device by using the water level and the installation position of each speed sensor in the monitoring device;
when the submerged speed sensor exists, acquiring the flow rate measured by the submerged speed sensor and a preset second flow formula;
matching according to the number of the submerged speed sensors to obtain a correction coefficient and a comprehensive correction coefficient in the second flow formula;
and inputting the flow speed measured by the submerged speed sensor, the correction coefficient and the comprehensive correction coefficient into the second flow formula to obtain the overflow flow of the water cross section in the monitoring device.
With reference to the first embodiment of the first aspect, in a second embodiment of the first aspect, the flow monitoring method further includes the steps of:
when the submerged speed sensor does not exist, acquiring a preset third flow formula;
and calculating the overflow flow of the water passing section in the monitoring device by using the third flow formula.
With reference to the first aspect, in a third implementation manner of the first aspect, the first flow formula is:
in the first flow formula, B is the width of the water passing section; h is the height of the water passing section; n is the number of channels in the measuring device; h is1The distance between the 1 st speed sensor and the bottom of the measuring device is defined, wherein the 1 st speed sensor is closest to the bottom of the measuring device; h isn+1The distance between the nth speed sensor and the top of the measuring device is set, wherein the nth speed sensor is closest to the top of the measuring device; v. ofiThe flow speed measured by the ith speed sensor is the flow speed measured by the ith speed sensor, wherein i is more than or equal to 1 and less than or equal to n, and n is a positive integer; a is v1The correction coefficient of (2); c is vnThe correction coefficient of (2);
is a comprehensive correction coefficient of the flow speed.
With reference to the first embodiment of the first aspect, in a fourth embodiment of the first aspect, the second flow rate formula is:
in the second flow formula, B is the width of the water cross section; h isIThe height of the water passing section; n is the number of channels in the measuring device; h is1The distance between the 1 st speed sensor and the bottom of the measuring device is defined, wherein the 1 st speed sensor is closest to the bottom of the measuring device; delta L is the distance between the mth speed sensor and the (m +1) th speed sensor, wherein m is more than or equal to 1 and less than or equal to n-1; a is v1The correction coefficient of (2); c is vnThe correction coefficient of (2);
is a comprehensive correction coefficient of the flow speed.
With reference to the second embodiment of the first aspect, in a fifth embodiment of the first aspect, the third flow rate formula is:
in the third flow formula, C is 0.602+0.083h/P, h is the height of the water passing section, P is the distance from the bottom of the channel to the inner wall surface of the bottom of the box body, g is the weight acceleration, b is the width of the water passing section, hεThe distance from the bottom of the channel to the water surface.
With reference to the third and fourth embodiments of the first aspect, in the sixth embodiment of the first aspect, the correction coefficients a and c and the comprehensive correction coefficient in the first and second flow rate formulas
The determination method comprises the following steps:
acquiring the flow velocity of each velocity sensor at different sound channel angles under the same working condition, and obtaining the functional relation between the flow velocity and the installation position of the velocity sensor through fitting;
bringing the installation position of the 1 st speed sensor into the functional relation to obtain the flow speed of the 1 st speed sensor, and bringing the installation position of the nth speed sensor into the functional relation to obtain the flow speed of the nth speed sensor;
obtaining a corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the 1 st speed sensor to the average inlet flow speed; obtaining c corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the nth speed sensor to the average inlet flow speed;
inputting a and c corresponding to different sound channel angles under the working condition into a preset flow velocity formula to obtain the sound channel angles corresponding to the different sound channel angles under the working condition
Corresponding angles a, c and c of different sound channels under different working conditions,
Averaging to obtain correction coefficients a and c and comprehensive correction coefficient in the first flow formula and the second flow formula
According to a second aspect, an embodiment of the present invention provides a flow monitoring device, including:
the acquisition module is used for acquiring the water level in the monitoring device;
the judging module is used for judging whether the water level reaches the height of the monitoring device or not;
the first processing module is used for acquiring a preset first flow formula and flow rates measured by the speed sensors when the water level reaches the height of the monitoring device;
and the second processing module is used for inputting the flow speed measured by each speed sensor into the first flow formula to obtain the overflow flow of the water cross section in the monitoring device.
According to a third aspect, an embodiment of the present invention provides a flow monitoring device, including: the flow monitoring method comprises a water level collector, a speed collector, a memory and a processor, wherein the water level collector, the speed collector, the memory and the processor are in communication connection with each other, a computer instruction is stored in the memory, and the processor executes the computer instruction so as to execute the flow monitoring method in the first aspect or any implementation manner of the first aspect.
According to a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to execute the flow monitoring method according to the first aspect or any one of the implementation manners of the first aspect.
Drawings
The features and advantages of the present invention will be more clearly understood by reference to the accompanying drawings, which are illustrative and not to be construed as limiting the invention in any way, and in which:
FIG. 1 is a schematic flow chart of a flow monitoring method in example 1 of the present invention;
FIG. 2 is a schematic flow chart of a flow monitoring method in embodiment 2 of the present invention;
fig. 3 is a schematic structural view of a flow rate monitoring device in embodiment 3 of the present invention;
FIG. 4 is a schematic view showing the arrangement of ultrasonic sound channels in embodiment 4 of the present invention;
FIG. 5 is a schematic diagram of a prior art computational model;
fig. 6 is a schematic view of the installation position of the speed sensor in embodiment 4 of the present invention;
fig. 7 is a schematic structural diagram of a flow rate monitoring device in embodiment 5 of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, 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 some, but not all, embodiments of the present invention. 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.
Example 1
s101: and acquiring the water level in the monitoring device. As a specific embodiment, the monitoring device may have a structure similar to a rectangular parallelepiped, and the cross section of the water passing hole is rectangular, and has a width B and a height H. For example, an electronic water gauge may be provided in the monitoring device, and the water level in the monitoring device may be obtained through the electronic water gauge.
S102: and judging whether the water level reaches the height of the monitoring device.
S103: and when the water level reaches the height of the monitoring device, acquiring a preset first flow formula and the flow rate measured by each speed sensor.
Specifically, the first flow formula is as follows:
in the first flow formula, B is the width of the water passing section; h is the height of the water passing section; n is the number of sound channels in the measuring device and is equal to the number of speed sensors installed in the measuring device; h is1The distance between the 1 st speed sensor and the bottom of the measuring device is defined, wherein the 1 st speed sensor is closest to the bottom of the measuring device; h isn+1The distance between the nth speed sensor and the top of the measuring device is set, wherein the nth speed sensor is closest to the top of the measuring device; v. ofiThe flow speed measured by the ith speed sensor is the flow speed measured by the ith speed sensor, wherein i is more than or equal to 1 and less than or equal to n, and n is a positive integer; a is v1The correction coefficient of (2); c is vnThe correction coefficient of (2);
is a comprehensive correction coefficient of the flow speed. The overflow flow of the water cross section can be accurately calculated through the first flow formula in the embodiment 1 of the invention.
As a specific implementation mode, in the embodiment 1 of the present invention, the correction coefficients a and c and the comprehensive correction coefficient in the first flow rate formula
The determination method can adopt the following technical scheme:
acquiring the flow velocity of each velocity sensor at different sound channel angles under the same working condition, and obtaining the functional relation between the flow velocity and the installation position of the velocity sensor through fitting;
bringing the installation position of the 1 st speed sensor into the functional relation to obtain the flow speed of the 1 st speed sensor, and bringing the installation position of the nth speed sensor into the functional relation to obtain the flow speed of the nth speed sensor;
obtaining a corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the 1 st speed sensor to the average inlet flow speed; obtaining c corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the nth speed sensor to the average inlet flow speed;
inputting a and c corresponding to different sound channel angles under the working condition into a preset flow velocity formula to obtain the sound channel angles corresponding to the different sound channel angles under the working condition
Corresponding angles a, c and c of different sound channels under different working conditions,
Averaging to obtain correction coefficients a and c and comprehensive correction coefficient in the first flow formula and the second flow formula
Specifically, the flow rate formula is as follows:
in the flow rate formula, in the above-mentioned flow rate formula,
is the average inlet flow rate, in known amounts.
As a specific implementation, the n sound channel lines divide the water cross section into (n +1) incompletely equal regions, wherein the heights of the 2 nd region to the nth region are equal, and the height of the 1 st region and the height of the n +1 th region are less than the height of the 2 nd region. Because the inside gate delivery port that is close to of box and the distribution of speed of being close to box top are comparatively complicated under the circumstances that becomes the aperture, utilize inhomogeneous scheme of arranging, encrypted both ends to a certain extent, capture velocity of flow information that can be better, assurance measuring stability that can be better. The two speed sensors respectively positioned on the two inner walls of the monitoring device form a sound channel line, and the formed sound channel line is parallel to the bottom surface of the monitoring device.
S104: and inputting the flow speed measured by each speed sensor into the first flow formula to obtain the overflow flow of the water cross section in the monitoring device.
In the flow monitoring method provided in embodiment 1 of the present invention, the water level in the monitoring device is obtained, whether the water level reaches the height of the monitoring device is determined, and when the water level reaches the height of the monitoring device, a preset first flow formula and the flow rates measured by the speed sensors are obtained; the flow velocity measured by each velocity sensor is input into the first flow formula to obtain the overflow flow of the water passing section in the monitoring device.
Example 2
An embodiment 2 of the present invention provides a flow monitoring method applied to a monitoring device provided with a plurality of speed sensors, where fig. 2 is a schematic flow diagram of the flow monitoring method in the embodiment 2 of the present invention, and as shown in fig. 2, the flow monitoring method in the embodiment 2 of the present invention includes the following steps:
s201: and acquiring the water level in the monitoring device.
S202: and judging whether the water level reaches the height of the monitoring device.
S203: and when the water level does not reach the height of the monitoring device, determining the submerged speed sensor in the monitoring device according to the water level and the installation position of each speed sensor in the monitoring device. When there is a submerged speed sensor, go to S204; when there is no submerged speed sensor, go to S207.
S204: when there is a submerged speed sensor, the flow rate measured by the submerged speed sensor and the preset second flow formula are obtained, and the process goes to S205.
Specifically, the second flow formula is as follows:
in the second flow formula, B is the width of the water cross section; h isIThe height of the water passing section is equal to the water level in the measuring device; n is the number of sound channels in the measuring device and is equal to the number of speed sensors installed in the measuring device; h is1The distance between the 1 st speed sensor and the bottom of the measuring device is defined, wherein the 1 st speed sensor is closest to the bottom of the measuring device; delta L is the distance between the mth speed sensor and the (m +1) th speed sensor, wherein m is more than or equal to 1 and less than or equal to n-1; a is v1The correction coefficient of (2); c is vnThe correction coefficient of (2);
is a comprehensive correction coefficient of the flow speed. The overflow flow of the water cross section can be accurately calculated through the first flow formula in the embodiment 2 of the invention.
Under the condition of non-full pipe, air exists in the box body, and the space between the top surface and the water surface in the box body is filled with the air, so that the flow state in the box body is influenced to a certain extent, and the accurate flow can be obtained by utilizing the calculation of the embodiment 2 of the invention.
As a specific implementation mode, in the embodiment 2 of the present invention, the correction coefficients a and c and the comprehensive correction coefficient in the first flow rate formula
The determination method can adopt the following technical scheme:
acquiring the flow velocity of each velocity sensor at different sound channel angles under the same working condition, and obtaining the functional relation between the flow velocity and the installation position of the velocity sensor through fitting;
bringing the installation position of the 1 st speed sensor into the functional relation to obtain the flow speed of the 1 st speed sensor, and bringing the installation position of the nth speed sensor into the functional relation to obtain the flow speed of the nth speed sensor;
obtaining a corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the 1 st speed sensor to the average inlet flow speed; obtaining c corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the nth speed sensor to the average inlet flow speed;
inputting a and c corresponding to different sound channel angles under the working condition into a preset flow velocity formula to obtain the sound channel angles corresponding to the different sound channel angles under the working condition
Corresponding angles a, c and c of different sound channels under different working conditions,
Averaging to obtain correction coefficients a and c and comprehensive correction coefficient in the first flow formula and the second flow formula
Specifically, the flow rate formula is as follows:
in the flow rate formula, in the above-mentioned flow rate formula,
is the average inlet flow rate, in known amounts.
S205: and matching the correction coefficient and the comprehensive correction coefficient in the second flow formula according to the number of the submerged speed sensors, and turning to S206.
Wherein the correction coefficient in the second flow formula is denoted by v1The correction coefficient a of (a); v. ofnThe correction coefficient c of (a); the integrated correction factor in the second flow equation refers to the integrated correction factor for the flow rate
S206: and inputting the flow speed measured by the submerged speed sensor, the correction coefficient and the comprehensive correction coefficient into the second flow formula to obtain the overflow flow of the water cross section in the monitoring device.
S207: when there is no submerged speed sensor, a preset third flow formula is obtained, and the process goes to S208.
Specifically, the third flow formula is as follows:
in the third flow formula, C is 0.602+0.083h/P, h is the height of the water passing section, P is the pressure above the water surface, g is the weight acceleration, b is the width of the water passing section, and h isεThe distance from the bottom of the channel to the water surface.
S208: and calculating the overflow flow of the water passing section in the monitoring device by using the third flow formula.
According to the flow monitoring method provided by the embodiment 2 of the invention, a proper flow formula can be obtained by judging the water level in the monitoring device and the number of the submerged speed sensors, so that the calculated overflow flow of the water-passing section is more accurate.
Example 3
Embodiment 3 of the present invention provides a flow rate monitoring device, and fig. 3 is a schematic structural diagram of the flow rate monitoring device in embodiment 3 of the present invention, and as shown in fig. 3, the flow rate monitoring device in embodiment 3 of the present invention includes: the device comprises an acquisition module 30, a judgment module 32, a first processing module 34 and a second processing module 36.
The obtaining module 30 is configured to obtain a water level in the monitoring device.
And the judging module 32 is used for judging whether the water level reaches the height of the monitoring device.
And the first processing module 34 is configured to obtain a preset first flow formula and flow rates measured by the speed sensors when the water level reaches the height of the monitoring device.
And the second processing module 36 is configured to input the flow rates measured by the speed sensors into the first flow formula to obtain the overflow flow of the water cross section in the monitoring device.
Further, in the flow rate monitoring device according to embodiment 3 of the present invention, the first processing module 34 is further configured to determine, when the water level does not reach the height of the monitoring device, a submerged speed sensor in the monitoring device by using the water level and the installation position of each speed sensor in the monitoring device; when the submerged speed sensor exists, acquiring the flow rate measured by the submerged speed sensor and a preset second flow formula; and matching the correction coefficient and the comprehensive correction coefficient in the second flow formula according to the number of the submerged speed sensors.
The second processing module 36 is further configured to obtain a correction coefficient and a comprehensive correction coefficient in the second flow formula according to the number of the submerged speed sensors; and inputting the flow speed measured by the submerged speed sensor, the correction coefficient and the comprehensive correction coefficient into the second flow formula to obtain the overflow flow of the water cross section in the monitoring device.
Example 4
To illustrate the flow monitoring method of the present invention in more detail, a more detailed embodiment is given.
1. Speed calculation
Fig. 4 is a schematic layout diagram of the ultrasonic sound channels in embodiment 4 of the present invention. In which fig. 4(a) is a top view, fig. 4(b) is a middle cross-sectional view, and fig. 4(c) is a velocity profile. The sound wave is propagated in the fluid, the propagation speed of the sound wave is increased due to the superposition of components in the water flow direction in the downstream direction, the water flow speed is reduced due to the influence of the components of the water flow speed in the upstream direction, different propagation time exists due to the change of the speed in the same propagation distance, the flow speed is obtained by using the relation between the propagation speed difference and the flow speed of the measured fluid, the propagation speed difference method is called, and the average flow speed V of the water flow line of each channel can be obtained by the following calculation model.
Suppose A-The time of B sound wave propagating in the downstream direction is t1And the time of the reverse propagation of the B-A sound wave is t2Then, there are:
the time of counter-current propagation was:
the downstream propagation time is:
further assume that:
ΔT＝t2-t1(1-3)
the following compounds are obtained by (1-1), (1-2), (1-3) and (1-4):
the 5 quantities related to the speed in the formula are changed into 4, and the sound speed does not need to be considered in the calculation, so that the errors of the temperature and the sound speed on the measurement results are reduced as far as possible. The water flow speed of different water layers of the channel is different and changes along with the change of the liquid level, so the expression of the equation is written as follows:
in the formula (1-6), LiThe length of the channel through which the sound wave propagates, θ, the acoustic path angle, i, is the number of channels. Because the flow velocity actually measured by the ultrasonic probe is the linear average flow velocity of the line on the channel, the average flow velocity of the water cross section needs to be measured in general engineering, so the average flow velocity needs to be corrected, and if the correction coefficient is K:
Vi＝KVS(1-7)
Vion different sound channelsLinear average flow velocity of VS-cross-sectional average flow velocity.
2. First flow formula derivation process
I. Assuming that the cross section of the water flow is rectangular, the width is B and the height is H, n sound channels are utilized to divide the cross section of the water flow into n +1 equal parts, and the average flow velocity of the cross section is
The prior art generally proceeds as follows, and fig. 5 is a schematic diagram of a prior art computational model.
Defining:
γ＝Y/H (1-8)
order to
The flow calculation formula can be written as:
assume that n sound channel lines divide the water cross section into (n +1) incompletely equal regions, the speed sensors are arranged at equal intervals in the central section, and the installation distance between the center and the end is different from the installation distance between the center and the end, specifically, h2 ═ h3 … … ═ hn, h1 < h2, and hn +1 < h2, and fig. 6 is a schematic view of the installation positions of the speed sensors in embodiment 4 of the present invention.
In example 4 of the present invention, the calculation formula of the flow rate and the cross-sectional average flow velocity obtained by deformation according to 1-8 to 1-13 is:
deforming 1-14 to obtain 1-15.
Combined by 1-15
To obtain 1-16.
Defining:
1-20, 1-21, 1-22 and 1-24 can be obtained by carrying out derivation according to 1-16 combined with 1-17, 1-18 and 1-19.
B, H, h therein1、hn+1For a known parameter, v1、viMeasured by a speed sensor, a and c are numerical experiment calculated values,
is the integrated velocity correction factor.
Specifically, a, c, and C in the above formulas 1 to 23,
For the requested quantity, the correction factors a, c and the overall correction factor in the first flow equation are determined
The following method is adopted:
and acquiring the flow velocity of each velocity sensor at different sound channel angles under the same working condition, and obtaining the functional relation between the flow velocity and the installation position of the velocity sensor through fitting. For example, setting conditions, namely boundary conditions, of two working conditions, namely 5# and 6# are given in table 1, setting modes of different channel angles under the same working condition are given in table 2, and a functional relation between the flow velocity and the installation position of the speed sensor obtained by fitting is given in table 2, wherein the functional relation is a cubic function or a logarithmic function, e in table 2 represents a constant, and f, g and h represent coefficients.
TABLE 1 calculation parameter Table
| Working conditions | Opening degree of gate | Upstream depth of water (m) | Downstream depth of water (m) | Track angle (°) | Average inlet flow velocity (m/s) | Outlet (Pa) |
| 5# | 0.2 | 0.85 | 0.11 | 48 | 0.264 | Free outflow |
| 6# | 0.4 | 0.92 | 0.54 | 48 | 0.539 | Free outflow |
TABLE 2 a, c, k value calculation Table
And bringing the installation position of the 1 st speed sensor into the functional relation to obtain the flow speed of the 1 st speed sensor, and bringing the installation position of the nth speed sensor into the functional relation to obtain the flow speed of the nth speed sensor.
Obtaining a corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the 1 st speed sensor to the average inlet flow speed; and obtaining c corresponding to different sound channel angles under the working condition by utilizing the ratio of the flow speed of the nth speed sensor to the average inlet flow speed. Where the average inlet flow rate is a known quantity.
Inputting a and c corresponding to different sound channel angles under the working condition into a preset flow velocity formula to obtain the sound channel angles corresponding to the different sound channel angles under the working condition
Corresponding angles a, c and c of different sound channels under different working conditions,
Averaging to obtain correction coefficients a and c and comprehensive correction coefficient in the first flow formula and the second flow formula
Specifically, the following are defined:
aja values, c, for different channel angles in the j-th operating modejC values, k, for different channel angles in the j-th operating modejFor k values corresponding to different channel angles for the j-th condition, the correlation coefficients are all found to be within an acceptable range by SPSS curve regression.
3. Flow monitoring during full pipe
In embodiment 4 of the present invention, when the water level in the monitoring device reaches the height of the monitoring device, it is called full pipe. The following can be obtained through the calculation: a is 1.1504, c is 1.116,
e in the table represents a constant, the linear speed on the sound channel shows logarithmic change along with the installation position of the sound channel when the opening degree of the regression function model is 0.2, 0.4 and 0.6, and shows quadratic function change when the opening degree is 0.8 and 1.0,at this time, the values of a-1.1504, c-1.116,
and all the box geometrical structure parameters are substituted into (1-24) and simplified to obtain:
further, the cross-sectional average flow velocity is calculated by the following formula:
4. flow monitoring during non-full pipe
In embodiment 4 of the present invention, when the water level in the monitoring device does not reach the height of the monitoring device, it is called a non-full pipe. Firstly, the number of the submerged sound channels (namely the number of the submerged speed sensors) is judged, and the electronic water gauge is arranged under the condition that the pipe is not full, and the water level measured by the electronic water gauge to the inside of the box body is hIBy comparison of hIAnd (h)1+ n Δ L) the flow calculation formula for the case of different channels is selected, where Δ L is 90 mm.
At non-full pipe, the second flow calculation formula is:
in the formula, B is the net width of the box body, h1Initial position of speed sensor installation, hIThe water depth, DeltaL-speed sensor installation center distance, other letters and full pipe mean the same. When the gate is fully opened, assuming that the water level inside the equipment square box body is 0.11m, 0.20m, 0.29m, 0.38m, 0.47m, 0.56m, 0.65m and 0.74m, the number of sound channels submerged in each water depth depends onThe times are 1, 2, 3, 4, 5, 6, 7 and 8. The upper layer is air and the lower layer is water. When the pipe is not full, a, c,
The determination method of (2) is the same as that in the case of the non-full pipe, except that the setting of the boundary condition is not used. When the pipe is not full, the boundary conditions are set to be that the boundary condition of the upper layer is a pressure inlet, the incoming flow direction is a pressure inlet, the boundary condition of the outlet is a pressure outlet, the bottom and two sides of the box are fixed wall surface boundaries, the calculation is carried out by CFD software, and the calculation of water-gas two-phase flow is carried out by combining VOF and a standard K-epsilon model to obtain a, c, B, C, B, C,
The specific results are as follows:
when h is generated1≤hI<h1+ Δ L, submerging 1 channel: 1.116, 0.941 for a and c,
when h is generated1+ΔL≤hI<h1+2 Δ L, submerging 2 channels: a is 1.068, c is 0.721,
when h is generated1+2ΔL≤hI<h1+3 Δ L, submerging 3 channels: 1.051, c 0.804,
when h is generated1+3ΔL≤hI<h1+4 Δ L, submerging 4 channels: a is 1.080, c is 0.924,
when h is generated1+4ΔL≤hI<h1+5 Δ L, submerging 5 channels: a is 1.143, c is 0.851,
when, h1+5ΔL≤hI<h1+6 Δ L, submerging 6 channels: a is 1.215, c is 0.905,
when h is generated1+6ΔL≤hI<h1+7 Δ L, submerging 7 channels: 1.1804, c 1.016,
when h is generated1+7ΔL≤hI<H, submerge 8 channels: a is 1.1504, c is 1.116,
in addition, when 0 is less than or equal to hI<h1I.e., when 0 channels are drowned out,
wherein G is the weight acceleration, C is 0.602+0.083h/p, b is the width of the water passing section (in this embodiment, b is 0.762m), and h is the height of the water passing section, in this embodiment, h is G0+0.04, wherein G0The distance from the bottom of the tank to the water surface of the tank, and the distance from the bottom of the channel to the inner wall surface of the bottom of the tank, P being 0.04m in this example of the facility.
Namely, the flow calculation formula is as follows:
5. flow comparative analysis with rotary propeller flow measuring instrument
To verify the accuracy of the first and second flow equations of the present invention, experimental values calculated according to the first/second flow equations were compared with simulated values obtained by a rotometer, the comparison results are shown in table 3.
It can be seen from table 3 that the flow calculated by using the calculation model is compared and calibrated with the calibrated propeller flow meter, and the error between the flow and the calibrated propeller flow meter is within 10%, which indicates that the calculation method is more accurate.
Computational model calculation flow verification as used in Table 3
Example 5
An embodiment of the present invention further provides a flow monitoring device, as shown in fig. 7, the flow monitoring device may include a water level collector 70, a speed collector 71, a processor 72, and a memory 73, where the processor 72 and the memory 73 may be connected by a bus or in another manner, and fig. 7 takes the connection by the bus as an example.
The processor 72 may be a Central Processing Unit (CPU). The Processor 72 may also be other general purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components, or combinations thereof.
The memory 73 is a non-transitory computer readable storage medium, and can be used to store non-transitory software programs, non-transitory computer executable programs, and modules, such as program instructions/modules corresponding to the flow monitoring method in the embodiment of the present invention (for example, the first obtaining module 30, the determining module 32, the second obtaining module 34, and the first calculating module 3 shown in fig. 3). The processor 72 executes various functional applications and data processing of the processor by executing non-transitory software programs, instructions and modules stored in the memory 73, namely, implements the flow monitoring method in the above method embodiment.
The memory 73 may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created by the processor 72, and the like. Further, the memory 73 may include high speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory 73 may optionally include memory located remotely from the processor 72, and such remote memory may be connected to the processor 72 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The one or more modules are stored in the memory 73 and, when executed by the processor 72, perform the flow monitoring method in the embodiment shown in fig. 1-2.
The details of the flow monitoring device may be understood with reference to the corresponding descriptions and effects in the embodiments shown in fig. 1 to fig. 2, and are not described herein again.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, abbreviated as HDD) or a Solid State Drive (SSD), etc.; the storage medium may also comprise a combination of memories of the kind described above.
Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations fall within the scope defined by the appended claims.
Claims (8)
1. A flow monitoring method applied to a monitoring device provided with a plurality of speed sensors is characterized by comprising the following steps:
acquiring the water level in the monitoring device;
judging whether the water level reaches the height of the monitoring device or not;
when the water level reaches the height of the monitoring device, acquiring a preset first flow formula and flow rates measured by all the speed sensors;
inputting the flow speed measured by each speed sensor into the first flow formula to obtain the overflow flow of the water cross section in the monitoring device;
wherein the first flow equation is:
in the first flow formula, B is the width of the water passing section; h is the height of the water passing section; n is the number of channels in the measuring device; h is1The distance between the 1 st speed sensor and the bottom of the measuring device is defined, wherein the 1 st speed sensor is closest to the bottom of the measuring device; h isn+1The distance between the nth speed sensor and the top of the measuring device is set, wherein the nth speed sensor is closest to the top of the measuring device;v i the flow speed measured by the ith speed sensor is the flow speed measured by the ith speed sensor, wherein i is more than or equal to 1 and less than or equal to n, and n is a positive integer;ais composed ofv 1 The correction coefficient of (2);cis composed ofv n The correction coefficient of (2);
is a comprehensive correction coefficient of the flow speed.
2. The flow monitoring method according to claim 1, further comprising:
when the water level does not reach the height of the monitoring device, determining submerged speed sensors in the monitoring device according to the water level and the installation positions of the speed sensors in the monitoring device;
when the submerged speed sensor exists, acquiring the flow rate measured by the submerged speed sensor and a preset second flow formula;
matching according to the number of the submerged speed sensors to obtain a correction coefficient and a comprehensive correction coefficient in the second flow formula;
inputting the flow speed measured by the submerged speed sensor, the correction coefficient and the comprehensive correction coefficient into the second flow formula to obtain the overflow flow of the water cross section in the monitoring device;
wherein the second flow equation is:
in the second flow formula, B is the width of the water cross section; h isIThe height of the water passing section; n is the number of channels in the measuring device; h is1The distance between the 1 st speed sensor and the bottom of the measuring device is defined, wherein the 1 st speed sensor is closest to the bottom of the measuring device; h isn+1The distance between the nth speed sensor and the top of the measuring device is set, wherein the nth speed sensor is closest to the top of the measuring device; delta L is the distance between the mth speed sensor and the (m +1) th speed sensor, wherein m is more than or equal to 1 and less than or equal to n-1;ais composed ofv 1 The correction coefficient of (2);cis composed ofv n The correction coefficient of (2);
is a comprehensive correction coefficient of the flow speed.
3. The flow monitoring method according to claim 2, further comprising:
when the submerged speed sensor does not exist, acquiring a preset third flow formula;
calculating the overflow flow of the water passing section in the monitoring device by using the third flow formula;
the third flow equation is:
in the third flow formula, h is the height of the water passing section, P is the distance from the bottom of the channel to the inner wall surface of the bottom of the box body, g is the weight acceleration, b is the width of the water passing section,
the distance from the bottom of the channel to the water surface.
4. The method of claim 1, wherein the correction factors a, c and the integrated correction factor in the first flow equation
The determination method comprises the following steps:
acquiring the flow velocity of each velocity sensor at different sound channel angles under the same working condition, and obtaining the functional relation between the flow velocity and the installation position of the velocity sensor through fitting;
bringing the installation position of the 1 st speed sensor into the functional relation to obtain the flow speed of the 1 st speed sensor, and bringing the installation position of the nth speed sensor into the functional relation to obtain the flow speed of the nth speed sensor;
obtaining the ratio of the flow speed of the 1 st speed sensor to the average flow speed of the inlet under the working conditiona(ii) a Obtaining the corresponding sound channel angles under the working condition by utilizing the ratio of the flow speed of the nth speed sensor to the average inlet flow speedc；
Inputting a and c corresponding to different sound channel angles under the working condition into a preset flow velocity formula to obtain the sound channel angles corresponding to the different sound channel angles under the working condition
；
For different workersWith angular correspondence of different vocal tracta、c、 
Averaging to obtain correction coefficients a and c and comprehensive correction coefficient in the first flow formula
。
5. The flow monitoring method of claim 2, wherein the correction factors a, c and the integrated correction factor in the second flow equation
The determination method comprises the following steps:
acquiring the flow velocity of each velocity sensor at different sound channel angles under the same working condition, and obtaining the functional relation between the flow velocity and the installation position of the velocity sensor through fitting;
bringing the installation position of the 1 st speed sensor into the functional relation to obtain the flow speed of the 1 st speed sensor, and bringing the installation position of the nth speed sensor into the functional relation to obtain the flow speed of the nth speed sensor;
obtaining the ratio of the flow speed of the 1 st speed sensor to the average flow speed of the inlet under the working conditiona(ii) a Obtaining the corresponding sound channel angles under the working condition by utilizing the ratio of the flow speed of the nth speed sensor to the average inlet flow speedc；
Inputting a and c corresponding to different sound channel angles under the working condition into a preset flow velocity formula to obtain the sound channel angles corresponding to the different sound channel angles under the working condition
；
For angles of different sound channels under different conditionsa、c、 
Averaging to obtain correction coefficients a and c and comprehensive correction coefficient in the second flow formula
。
6. A flow monitoring device, comprising:
the acquisition module is used for acquiring the water level in the monitoring device;
the judging module is used for judging whether the water level reaches the height of the monitoring device or not;
the first processing module is used for acquiring a preset first flow formula and flow rates measured by the speed sensors when the water level reaches the height of the monitoring device;
the second processing module is used for inputting the flow speed measured by each speed sensor into the first flow formula to obtain the overflow flow of the water cross section in the monitoring device;
wherein the first flow equation is:
in the first flow formula, B is the width of the water passing section; h is the height of the water passing section; n is the number of channels in the measuring device; h is1The distance between the 1 st speed sensor and the bottom of the measuring device is defined, wherein the 1 st speed sensor is closest to the bottom of the measuring device; h isn+1The distance between the nth speed sensor and the top of the measuring device is set, wherein the nth speed sensor is closest to the top of the measuring device;v i the flow speed measured by the ith speed sensor is the flow speed measured by the ith speed sensor, wherein i is more than or equal to 1 and less than or equal to n, and n is a positive integer;ais composed ofv 1 The correction coefficient of (2);cis composed ofv n The correction coefficient of (2);
is a comprehensive correction coefficient of the flow speed.
7. A flow monitoring device, comprising:
the flow monitoring device comprises a water level collector, a speed collector, a memory and a processor, wherein the water level collector, the speed collector, the memory and the processor are in communication connection with each other, computer instructions are stored in the memory, and the processor executes the computer instructions so as to execute the flow monitoring method according to any one of claims 1 to 5.
8. A computer-readable storage medium having stored thereon computer instructions for causing a computer to perform the flow monitoring method of any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910584822.1A CN110426088B (en) | 2019-07-01 | 2019-07-01 | Flow monitoring method, device and equipment and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910584822.1A CN110426088B (en) | 2019-07-01 | 2019-07-01 | Flow monitoring method, device and equipment and readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110426088A CN110426088A (en) | 2019-11-08 |
| CN110426088B true CN110426088B (en) | 2020-05-05 |
Family
ID=68409918
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910584822.1A Active CN110426088B (en) | 2019-07-01 | 2019-07-01 | Flow monitoring method, device and equipment and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110426088B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114280326B (en) * | 2021-12-23 | 2023-08-25 | 广西大学 | Karst depression waterlogging water falling hole flow velocity measuring device and overflow capacity monitoring method |
| CN118310594B (en) * | 2024-06-11 | 2024-08-30 | 安徽汉威电子有限公司 | Strong filtering high-precision metering intelligent water meter |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5271675A (en) * | 1992-10-22 | 1993-12-21 | Gas Research Institute | System for characterizing pressure, movement, temperature and flow pattern of fluids |
| JP2009264906A (en) * | 2008-04-24 | 2009-11-12 | Ricoh Elemex Corp | Flow meter |
| US9816848B2 (en) * | 2014-01-23 | 2017-11-14 | Ultimo Measurement Llc | Method and apparatus for non-invasively measuring physical properties of materials in a conduit |
| CN108021762B (en) * | 2017-12-14 | 2021-07-16 | 黑龙江省水利水电勘测设计研究院 | Method and device for determining width of channel bottom |
| CN108519129B (en) * | 2018-04-25 | 2024-01-19 | 中国灌溉排水发展中心 | Channel flow monitoring method and system |
-
2019
- 2019-07-01 CN CN201910584822.1A patent/CN110426088B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110426088A (en) | 2019-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ma et al. | Shock-capturing non-hydrostatic model for fully dispersive surface wave processes | |
| CN108254032B (en) | River ultrasonic time difference method flow calculation method | |
| Cassan et al. | Experimental and numerical investigation of flow under sluice gates | |
| CN110426088B (en) | Flow monitoring method, device and equipment and readable storage medium | |
| CN105051503B (en) | flow measurement | |
| CN104390680B (en) | Transit time ultrasonic ripple heat quantity flow error correcting method | |
| CN118464133A (en) | Flow measurement | |
| CN108775936B (en) | Flow metering device, metering method and measurement and control integrated gate system | |
| CN110135069B (en) | Method and device for acquiring silt characteristics during water delivery of water delivery tunnel and computer equipment | |
| Sharifipour et al. | Numerical investigation of flow field and flowmeter accuracy in open-channel junctions | |
| CN113553737B (en) | Valve flow prediction method based on valve pressure difference | |
| Barjastehmaleki et al. | Spillway stilling basins lining design via Taylor hypothesis | |
| Schulte | Water flow through a 20-pore perforation plate in vessels of Liquidambar styraciflua | |
| CN111122116A (en) | Experimental device for simulating coarse sublayer flow velocity distribution and method for determining coarse sublayer flow velocity distribution | |
| Thanh et al. | Physical and numerical model of flow through the spillways with a breast wall | |
| Goel et al. | Open channel flow measurement of water by using width contraction | |
| Mulahasan et al. | Simulation of flow around a permeable dike using physical and 3D-CFD models | |
| Demirel | Measured and simulated flow downstream of the submerged sluice gate | |
| Yang et al. | Discharge estimation of the Shin-Yuan Canal using indirect method | |
| CN114965199A (en) | KaF-RBF collaborative fusion-based suspended load sand content measurement method | |
| CN114485863A (en) | Flow error correction method, system, computer and medium for ultrasonic water meter | |
| Helmi et al. | Experimental and numerical investigations of flow through free double baffled gates | |
| Abbaspuor et al. | Experimental and numerical evaluation of discharge capacity of sharp-crested triangular plan form weirs | |
| Petrut et al. | Flow rate computation of highly turbulent pipe flows using wide-band signals and matched filter-based approach in ultrasonic multi-element configuration | |
| Cassan et al. | Experimental and numerical studies of the flow structure generated by a submerged sluice gate |
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 |