CN113074786A - Curve channel flow measuring method and device - Google Patents
Curve channel flow measuring method and device Download PDFInfo
- Publication number
- CN113074786A CN113074786A CN202110272485.XA CN202110272485A CN113074786A CN 113074786 A CN113074786 A CN 113074786A CN 202110272485 A CN202110272485 A CN 202110272485A CN 113074786 A CN113074786 A CN 113074786A
- Authority
- CN
- China
- Prior art keywords
- channel
- flow velocity
- curve
- section
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/002—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow wherein the flow is in an open channel
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention provides a curve channel flow measuring method and device. Wherein, the method comprises the following steps: determining a flow velocity distribution model of a curve channel; determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the parting line characteristics; and dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method, and obtaining the flow information of the curve channel based on a speed area method. The flow measuring process can be simplified under different hydraulic conditions by adopting the method for measuring the flow of the curved channel, and the accuracy and efficiency of flow measurement of the curved channel are effectively improved.
Description
Technical Field
The invention relates to the technical field of computer application, in particular to a method and a device for measuring flow of a curve channel. In addition, an electronic device and a non-transitory computer readable storage medium are also related.
Background
Measuring open channel flow is an important issue in hydraulic engineering. The most common open channel section type in the irrigated area is a trapezoidal section, and the outstanding advantages of the channel are regular section, stable side slope, small difficulty in construction technology, strong applicability and the like. However, the trapezoidal section is generally weak in expansion and freezing resistance, and the curved channel section has excellent hydraulic conditions of large flow, strong sand transportation capability and remarkable seepage prevention effect. In particular, as various curved channel sections (such as circular, U-shaped, parabolic, etc.) have been popularized and applied to irrigation drainage engineering and sewage treatment plants, the water flow motion law and hydraulic characteristics of open channels have become hot spots for current technicians to study. In the conventional specification, the channel section is generally divided by a vertical line, and hydraulic information such as flow velocity of a fixed point on the vertical line is measured.
Therefore, how to quickly and accurately measure the flow rate of the curved channel on the basis of determining the flow velocity distribution of the cross section of the curved channel in the flow rate calculation process of the curved channel becomes a technical problem to be solved by the technical personnel in the field.
Disclosure of Invention
Therefore, the invention provides a method and a device for measuring the flow of a curved channel, which aim to solve the problem that the current actual use requirement cannot be met due to lower accuracy and efficiency of measuring the flow of the section of the curved channel in the prior art.
The invention provides a curve channel flow measuring method, which comprises the following steps:
determining a flow velocity distribution model of a curve channel;
determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the parting line characteristics;
and dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method, and obtaining the flow information of the curve channel based on a speed area method.
Further, the expression corresponding to the curve-type channel flow velocity distribution model is as follows:
wherein: A. b, C is the correlation coefficient; u is local friction flow rate; nu is the motion viscosity coefficient of water; l is the distance of the curve channel normal from each point on the curve channel normal to the side wall; klIs the vertical distance from the side wall to the dividing line.
Further, the expression of the point position correspondence of the normal average flow velocity on the normal of the curvilinear channel is as follows:
wherein l1、l2Representing the point location for the normal average flow velocity,/1、l2The speed is the average speed of the speed measurement normal; klIs the distance from the curved channel bottom to the water surface line along the normal direction, and l is the normal average flow velocity representative point along KlTo the bottom of the curved trench.
Further, the method for dividing the cross-sectional area of the curved channel by using a preset three-line-polygon division method according to the position of the normal average flow velocity representative point and obtaining the flow information of the curved channel based on a velocity area method specifically comprises the following steps:
dividing the area of the section of the curved channel according to the position of the representative point of the normal average flow velocity by using a preset three-line-polygon division method, and determining a plurality of partition polygons corresponding to the section of the curved channel;
calculating the flow of the curved channel according to the plurality of partition polygons by using a preset channel flow calculation formula to obtain the flow information of the curved channel;
the channel flow calculation formula is as follows:
Qgeneral assembly＝2×[α·u1·(S1+S2)]+β·u2·(S3+S4)
Wherein u is1Is the flow velocity, u, of the target point on the first and third normal lines2Is the flow velocity of the target point on the second normal; alpha and beta are flow correction coefficients; s1、S2、S3And S4The areas of the plurality of partition polygons are respectively.
Further, the determining the position of the representative point of the normal average flow velocity of the curved channel section specifically includes: determining three corresponding normal lines on the section of the curve-type channel; respectively determining corresponding normal average flow velocity representative point positions on the three normals;
the method for determining the plurality of partition polygons corresponding to the curved channel section by utilizing a preset three-line-polygon partition method and dividing the area of the curved channel section according to the position of the normal average flow velocity representative point specifically comprises the following steps: connecting the normal average flow velocity representative points on the three normals according to a preset mode to form a plurality of corresponding triangles; determining an intersection point of the outer center position of the triangle and the water surface and an intersection point of the outer center position of the triangle and the canal bottom according to the perpendicular bisectors corresponding to the three sides of the triangle; taking the outer center position as a partition boundary point, and connecting the partition boundary point to obtain a corresponding partition polygon; wherein, each partition polygon includes a normal average flow velocity representative point.
The invention also provides a curve type channel flow measuring device, comprising:
the algorithm model determining unit is used for determining a curve type channel flow velocity distribution model;
the representative point position determining unit is used for determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the partition line characteristics;
and the channel flow calculating unit is used for dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method and obtaining the flow information of the curve channel based on a speed area method.
Further, the expression corresponding to the curve-type channel flow velocity distribution model is as follows:
wherein: A. b, C is the correlation coefficient; u is local friction flow rate; nu is the motion viscosity coefficient of water; l is the distance of the curve channel normal from each point on the curve channel normal to the side wall; klIs the vertical distance from the side wall to the dividing line.
Further, the expression of the point position correspondence of the normal average flow velocity on the normal of the curvilinear channel is as follows:
wherein l1、l2Representing the point location for the normal average flow velocity,/1、l2The speed is the average speed of the speed measurement normal; klIs the distance from the curved channel bottom to the water surface line along the normal direction, and l is the normal average flow velocity representative point along KlTo the bottom of the curved trench.
Further, the method for dividing the cross-sectional area of the curved channel by using a preset three-line-polygon division method according to the position of the normal average flow velocity representative point and obtaining the flow information of the curved channel based on a velocity area method specifically comprises the following steps:
dividing the area of the section of the curved channel according to the position of the representative point of the normal average flow velocity by using a preset three-line-polygon division method, and determining a plurality of partition polygons corresponding to the section of the curved channel;
calculating the flow of the curved channel according to the plurality of partition polygons by using a preset channel flow calculation formula to obtain the flow information of the curved channel;
the channel flow calculation formula is as follows:
Qgeneral assembly＝2×[α·u1·(S1+S2)]+β·u2·(S3+S4)
Wherein u is1Is the flow velocity, u, of the target point on the first and third normal lines2Is the flow velocity of the target point on the second normal; alpha and beta are flow correction coefficients; s1、S2、S3And S4The areas of the plurality of partition polygons are respectively.
Further, the representative point position determining unit is specifically configured to: determining three corresponding normal lines on the section of the curve-type channel; respectively determining corresponding normal average flow velocity representative point positions on the three normals;
the channel flow calculating unit is specifically configured to: connecting the normal average flow velocity representative points on the three normals according to a preset mode to form a plurality of corresponding triangles; determining an intersection point of the outer center position of the triangle and the water surface and an intersection point of the outer center position of the triangle and the canal bottom according to the perpendicular bisectors corresponding to the three sides of the triangle; taking the outer center position as a partition boundary point, and connecting the partition boundary point to obtain a corresponding partition polygon; wherein, each partition polygon includes a normal average flow velocity representative point.
The invention also provides a method for acquiring the flow velocity distribution model of the curve channel, which comprises the following steps:
determining an initial curve type channel flow velocity distribution model according to the curve type channel;
on the basis of the initial curve channel flow velocity distribution model, a curve channel full-section flow velocity distribution model is obtained based on a three-point parameter quantification method of a section division theory; the three-point parameter quantification method based on the section segmentation theory obtains a full-section flow velocity distribution model of a curve channel, and specifically comprises the following steps: determining three normal lines corresponding to the section of the curve-type channel, and respectively acquiring the flow velocities of target points on the three normal lines; inputting the flow velocity of the target point into the initial curve type channel flow velocity distribution model to obtain corresponding quantitative undetermined parameters; and determining the full-section flow velocity distribution model of the curve type channel based on the quantitative undetermined parameter and the initial curve type channel flow velocity distribution model.
Correspondingly, the invention also provides a device for obtaining the flow velocity distribution model of the curve channel, which comprises:
the initial flow velocity distribution model determining unit is used for determining an initial curve type channel flow velocity distribution model according to the curve type channel;
a full-section flow velocity distribution model determining unit, configured to obtain a full-section flow velocity distribution model of the curvilinear channel based on a three-point parameter quantification method of a section segmentation theory on the basis of the initial curvilinear channel flow velocity distribution model; the three-point parameter quantification method based on the section segmentation theory obtains a full-section flow velocity distribution model of a curve channel, and specifically comprises the following steps: determining three normal lines corresponding to the section of the curve-type channel, and respectively acquiring the flow velocities of target points on the three normal lines; inputting the flow velocity of the target point into the initial curve type channel flow velocity distribution model to obtain corresponding quantitative undetermined parameters; and determining the full-section flow velocity distribution model of the curve type channel based on the quantitative undetermined parameter and the initial curve type channel flow velocity distribution model.
The invention also provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the curve type channel flow measuring method.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the curvilinear channel flow determination method as described in any one of the above.
The method for measuring the flow of the curved channel can measure the flow velocity of the section of the curved channel under different hydraulic conditions according to the real flow velocity distribution condition of the section of the curved channel, simplifies the flow measurement process, has small error and effectively improves the accuracy and efficiency of flow measurement of the curved channel.
By adopting the method for acquiring the flow velocity distribution model of the curved channel, the flow velocity distribution model of the curved channel can be rapidly determined, so that the accuracy and the efficiency of flow measurement of the curved channel are improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a schematic flow chart of a method for measuring a flow rate of a curved channel according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a point location represented by the average flow velocity of the normal of a curved channel according to an embodiment of the present invention;
fig. 3 shows a curve channel flow Q of 17m according to an embodiment of the present invention3Flow distribution/h;
fig. 4 is a schematic position diagram of three normal lines in the method for obtaining a flow velocity distribution model of a curved channel according to the embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating experimental values of average flow velocity under different conditions in comparison with theory according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of an application of the three-line polygon partition method according to an embodiment of the present invention;
FIG. 7 is a schematic sectional view of a curved channel according to an embodiment of the present invention;
FIG. 8 is a schematic view of a complete flow chart of a curvilinear channel flow measurement provided by an embodiment of the present invention;
FIG. 9 is a schematic structural diagram of a curved channel flow measurement device according to an embodiment of the present invention;
fig. 10 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention;
FIG. 11 is a schematic flow chart illustrating a method for obtaining a flow velocity distribution model of a curved channel according to an embodiment of the present invention;
fig. 12 is a schematic structural diagram of an apparatus for obtaining a flow velocity distribution model of a curved channel according to an embodiment 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.
The following describes an embodiment of the method for measuring a flow rate through a curved channel according to the present invention in detail. As shown in fig. 1, which is a schematic flow chart of a method for measuring a flow rate of a curved channel according to an embodiment of the present invention, the specific process includes the following steps:
step 101: and determining a curve type channel flow velocity distribution model.
In the step, after the test data are analyzed, firstly, target data on a perpendicular bisector of the curve channel are extracted to analyze the flow velocity distribution rule, a logarithmic formula which is corrected by adding a wake function is selected to fit an empirical formula of the flow velocity distribution rule on the perpendicular bisector of the curve channel, a target function formula suitable for the flow velocity distribution rule on the perpendicular bisector is preliminarily established, and corresponding coefficients are determined through the test data. And determining an initial curve type channel flow velocity distribution model corresponding to the flow velocity distribution rule on the vertical line in the curve type channel according to the objective function formula. Further, the initial curve type channel flow velocity distribution model is verified, and if the verification is successful, the curve type channel flow velocity distribution model is obtained.
In a specific implementation process, the initial curvilinear channel flow velocity distribution model is verified, and if the verification is successful, the curvilinear channel flow velocity distribution model is obtained, wherein the specific implementation process comprises the following steps: quantifying a correlation coefficient in the initial curve type channel flow velocity distribution model through preset test data and simulation data; utilizing the initial curve type channel flow velocity distribution model to test the flow velocity distribution characteristics of the full section of the curve type channel and outputting a verification result; and if the fitted initial curve channel flow velocity distribution model is determined to describe the flow velocity distribution rule in the whole section of the curve channel to meet the preset condition according to the verification result, obtaining the curve channel flow velocity distribution model. The target function formula is popularized and analyzed, and the measured flow speed data of the full section of the inspection curve channel is brought into the target function formula to judge whether the flow speed data accords with the target function formula. The fitted objective function formula is substituted into other normals for comparison, and the result proves that the fitted flow velocity distribution empirical formula can better describe the flow velocity distribution rule in the whole section, so that the whole section can be described by using one empirical formula, and the description of the flow velocity distribution of the open channel water flow becomes more concise. The expression corresponding to the curve type channel flow velocity distribution model is as follows:
wherein: A. b, C is a correlation coefficient, and is obtained by fitting experimental data; u is local friction flow rate; nu is the motion viscosity coefficient of water; l is the distance of the curve channel normal from each point on the curve channel normal to the side wall; klIs the vertical distance from the side wall to the dividing line.
Step 102: and determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the partition line characteristics.
In the embodiment of the invention, the wake flow part of the curve-type channel flow velocity distribution model obtained by fitting can be changed as follows:
according to the trigonometric function relationship:
thus, equation (2) can be written as:
as shown in FIG. 2, K is the length and width of the film at any normal positionlDw, since dw tends towards 0, the flow through the infinitesimal can be expressed as:
Q＝u·Kl·dw (5)
in the infinitesimal body, an infinitesimal rectangle (shown in fig. 2 as a shaded portion) having a length and a width dl and dw, respectively, is taken, and the flow rate passing through the infinitesimal rectangle is:
dQ＝u·dw·dl (6)
integrating dl in the y' direction, substituting formula (4) to obtain formula I, substituting formula I to formula (6) to obtain formula II, and simultaneously solving the two formulas:
wherein t is
From a number of experimental fits, the values of the correlation coefficients a and B can be determined:
A＝2.5
therefore, the expression of the point position corresponding to the normal average flow speed on the normal of the curvilinear channel is as follows:
wherein l1、l2Representing the point location for the normal average flow velocity,/1、l2The speed is the average speed of the speed measurement normal; klIs the distance from the curved channel bottom to the water surface line along the normal direction, and l is the normal average flow velocity representative point along KlTo the bottom of the curved trench.
Step 103: and dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method, and obtaining the flow information of the curve channel based on a speed area method.
In the embodiment of the present invention, the method for obtaining the flow information of the curved channel based on the velocity area method includes the following specific steps: dividing the area of the section of the curved channel according to the position of the representative point of the normal average flow velocity by using a preset three-line-polygon division method, and determining a plurality of partition polygons corresponding to the section of the curved channel; and calculating the flow of the curve channel according to the plurality of partition polygons by using a preset channel flow calculation formula to obtain the flow information of the curve channel.
By the flow velocity area method, the channel flow calculation formula is as follows:
Qgeneral assembly＝2×[α·u1·(S1+S2)]+β·u2·(S3+S4) (10)
Wherein u is1The flow velocities at points on the first normal line (i.e., normal 1 in FIG. 6 or 7) and the third normal line (i.e., normal 3 in FIG. 6 or 7), u2Is the second normal (i.e. the method in fig. 6 or 7)Flow rates at various points on line 2); l1、l2Is determined by equation (9); alpha and beta are flow correction coefficients, for example, 0.9 is taken respectively; s1、S2、S3And S4The areas of the plurality of partition polygons are respectively.
The channel flow calculation formula can be expressed as:
Qgeneral assembly＝2×[α·u1·(S1+S2)]+β·u2·(S3+S4) (11)
The channel flow is calculated by the above formula, and the relative error e is as follows:
the actual channel flow refers to the flow displayed by the electromagnetic flowmeter in the test, and the theoretical calculation flow is the result obtained by calculation by using an expression.
Further, the determining of the normal average flow velocity representative point position of the curved channel section specifically includes: determining three corresponding normal lines on the section of the curve-type channel; and respectively determining the corresponding normal average flow velocity representative point positions on the three normals. Correspondingly, the area of the section of the curved channel is divided by using a preset three-line-polygon division method according to the position of the representative point of the normal average flow velocity, and a plurality of partition polygons corresponding to the section of the curved channel are determined, wherein the specific implementation process comprises the following steps: connecting the normal average flow velocity representative points on the three normals according to a preset mode to form a plurality of corresponding triangles; determining an intersection point of the outer center position of the triangle and the water surface and an intersection point of the outer center position of the triangle and the canal bottom according to the perpendicular bisectors corresponding to the three sides of the triangle; taking the outer center position as a partition boundary point, and connecting the partition boundary point to obtain a corresponding partition polygon; wherein, each partition polygon includes a normal average flow velocity representative point.
The device of the whole hydraulic circulating system can comprise a water supply pump and a water storage tankThe device comprises a water supply pipeline, an electric valve, a manual valve, an electromagnetic flowmeter, a water stabilizing tank, a turbulence grid, a curved organic glass channel, a water return pipeline, a PIV three-dimensional particle velocimeter, an ADV Doppler current meter, a mobile measuring needle and the like. In the practical implementation process, the model test aspect is as follows: a curve channel flow velocity distribution model is constructed, hydraulic factors such as flow, water depth and bottom slope are used as variables, the water depth in a curve channel is controlled by controlling a valve system on a water return pipeline, so that water flow in the whole channel section can easily reach an open channel uniform flow state, a selectable flow velocity measurement position is located in the middle and lower reaches (two thirds of the whole channel section) of the whole water channel section, and the flow velocity of the cross section of the channel under different hydraulic conditions is measured by using an ultrasonic current meter (ADV), a three-dimensional particle imaging system (PIV) and a high-speed photography dynamic imaging technology. And (3) numerical simulation aspect: the three-dimensional flow velocity distribution of the water flow of the curve channel can be simulated based on CFD calculation software Fluent numerical value, and the model is verified based on the test result. The test conditions are as follows: the water inlet of the curved channel is provided with the steady flow honeycomb, so that the effect of stabilizing water is achieved, the water outlet of the curved channel is provided with the tail gate, the required ideal water depth can be obtained by adjusting the tail gate, and the water flow of the curved channel is kept in a uniform flow state. The flow of the whole test operation is controlled by an electromagnetic flowmeter, and the required flow can be obtained by effectively combining the opening degree of a valve and the operation frequency. For example, the actually set test flow rates can be respectively 15m3/h-27m3/h、30m3/h-65m3H, and simulating the flow of 100m by using ANSYS Fluent fluid simulation software3/h、110m3/h、120m3The flow velocity distribution data under three different flow conditions are shown in fig. 3.
In addition, since the flow velocity distribution of the cross section of the curved channel should be symmetrically distributed about the perpendicular bisector due to the symmetry of the cross section, the flow velocity distribution of the channel on the half side may be measured only with the perpendicular bisector as a boundary when measuring the flow velocity. After the measurement density and the survey line distribution of the cross-section flow velocity are determined, the instantaneous flow velocity data of each survey point is recorded through a computer.
In the specific implementation process, due to the close proximity of waterRepresentative points of the surface are susceptible to water surface turbulence and in the following analysis, representative points near the bottom of the trench may be selected. Theoretical value (u) of point represented by normal average flow velocityt) As abscissa, test value (u)m) On the ordinate, as can be seen from FIG. 5, the flow rate data under each condition is distributed at ut/umAround the straight line of 1. Through comparative analysis, the test value is closer to the calculated value of the flow velocity distribution model of the curve channel, and the relative errors are within 10 percent, which shows that the flow velocity distribution model of the curve channel can accurately describe the real flow velocity distribution of the curve section.
According to an objective function formula of the curve type channel flow velocity distribution model, average flow velocity representative point positions can be obtained, and measuring points can be arranged at the average flow velocity representative points in a test. Three normals of a perpendicular bisector and a boundary line are taken in a curve type channel: respectively determining the positions of two normal average flow velocity representative points on the three normals; connecting the representative points to form a plurality of triangles, and making perpendicular bisectors of three sides of the triangles to determine the positions of outscenters of the triangles and the intersection points of the three sides of the triangles with the water surface and the channel bottom (as shown in fig. 6). The position of the outer center is the boundary point of the subarea, and the boundary points are connected to obtain subarea polygons, and each polygon contains a normal flow velocity representative point (as shown in fig. 7).
As shown in fig. 8, a flow velocity distribution rule of the cross section of the curved channel is proposed by introducing a partition theory of the cross section of the curved channel through a physical model test and simulation. And deducing to obtain a theoretical calculation expression of the representative point position of the average flow velocity of the section of the curve channel based on the flow velocity partition theory of the curve channel and in combination with a flow velocity distribution model of the section of the curve channel. The flow measuring method of the curved channel is obtained by reasonably dividing the section area of the curved channel by a three-line-polygon division method.
By adopting the method for measuring the flow of the curved channel, the flow of the section of the curved channel under different hydraulic conditions can be measured according to the real flow distribution condition of the section of the curved channel, the flow measuring process is simplified, the error is small, and the accuracy and the efficiency of the flow measurement of the curved channel are effectively improved.
Corresponding to the curve channel flow measuring method, the invention also provides a curve channel flow measuring device. Since the embodiment of the device is similar to the method embodiment, the description is simple, and the related points can be referred to the description of the method embodiment, and the embodiment of the curve type channel flow measuring device described below is only schematic. Fig. 9 is a schematic structural diagram of a curved channel flow measuring device according to an embodiment of the present invention. The curve type channel flow measuring device specifically comprises the following parts:
an algorithm model determining unit 901, configured to determine a curve-type channel flow velocity distribution model;
a representative point position determining unit 902, configured to determine a representative point position of a normal average flow velocity of a cross section of the curved channel according to the flow velocity distribution model of the curved channel, the flow velocity partition of the curved channel, and the characteristics of the partition line;
and the channel flow calculating unit 903 is configured to segment the area of the cross section of the curved channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon segmentation method, and obtain flow information of the curved channel based on a velocity area method.
By adopting the curve channel flow measuring device provided by the embodiment of the invention, the flow velocity of the section of the curve channel under different hydraulic conditions can be measured according to the real flow velocity distribution condition of the section of the curve channel, the flow measuring process is simplified, the error is small, and the flow measuring precision and efficiency of the curve channel are effectively improved.
Corresponding to the curve channel flow measuring method, the invention also provides electronic equipment. Since the embodiment of the electronic device is similar to the above method embodiment, the description is simple, and please refer to the description of the above method embodiment, and the electronic device described below is only schematic. Fig. 10 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention. The electronic device may include: a processor (processor)1001, a memory (memory)1002, and a communication bus 1003, wherein the processor 1001 and the memory 1002 communicate with each other via the communication bus 1003. The processor 1001 may invoke logic instructions in the memory 1002 to perform a curvilinear channel flow determination method comprising: determining a flow velocity distribution model of a curve channel; determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the parting line characteristics; and dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method, and obtaining the flow information of the curve channel based on a speed area method.
In addition, the logic instructions in the memory 1002 may be implemented in software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, an embodiment of the present invention further provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium, the computer program including program instructions, when the program instructions are executed by a computer, the computer being capable of executing the curvilinear channel flow determination method provided by the above-mentioned method embodiments, the method including: determining a flow velocity distribution model of a curve channel; determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the parting line characteristics; and dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method, and obtaining the flow information of the curve channel based on a speed area method. It should be noted that the electronic device refers to a computer device used in cooperation with a PIV three-dimensional particle velocimeter. The PIV three-dimensional particle velocimeter can acquire related flow velocity data and transmit the related flow velocity data to computer equipment to analyze the data by using the method. A PIV three-dimensional particle velocimeter is an instrument for realizing velocity measurement by using a transient, multi-point and non-contact fluid mechanics velocity measurement method. The PIV technology is characterized in that the velocity distribution information of a large number of spatial points can be recorded in the same transient state, the limitation of single-point flow measurement is broken through, and the PIV technology is a comprehensive result of the development of high and new technologies such as a laser technology, a digital signal processing technology, a chip technology, a computer technology and an image processing technology. The three-dimensional PIV system mainly comprises test tracer particles, a flow field illumination system, an image acquisition and processing system and the like.
In yet another aspect, an embodiment of the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, the computer program being implemented by a processor to perform the method for determining a channel flow according to the curve type provided in the foregoing embodiments, the method including: determining a flow velocity distribution model of a curve channel; determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the parting line characteristics; and dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method, and obtaining the flow information of the curve channel based on a speed area method.
In addition, the application also provides a method and a device for acquiring the flow velocity distribution model of the curve channel. As shown in fig. 11, the method for obtaining the curve-type channel flow velocity distribution model specifically includes the following steps:
step 1101: determining an initial curve type channel flow velocity distribution model according to the curve type channel;
step 1102: and on the basis of the initial curve channel flow velocity distribution model, obtaining a curve channel full-section flow velocity distribution model based on a three-point parameter quantification method of a section segmentation theory.
The three-point parameter quantification method based on the section segmentation theory obtains a full-section flow velocity distribution model of a curve channel, and specifically comprises the following steps: determining three normal lines corresponding to the cross section of the curved channel, as shown in fig. 4, where the three normal lines are respectively: perpendicular bisector 1, normal 2 at left side bottom (AD) curve length 1/2, and normal 3 at right side bottom (AE) curve length 1/2. Respectively acquiring the flow velocities of target points on the three normals, and inputting the flow velocities of the target points into the initial curve type channel flow velocity distribution model to obtain corresponding quantitative undetermined parameters; and determining the full-section flow velocity distribution model of the curve type channel based on the quantitative undetermined parameter and the initial curve type channel flow velocity distribution model.
As shown in fig. 12, the apparatus for obtaining the flow velocity distribution model of the curvilinear channel specifically includes the following components:
an initial flow velocity distribution model determining unit 1201, configured to determine an initial curved channel flow velocity distribution model according to a curved channel;
a full-section flow velocity distribution model determining unit 1202, configured to obtain a full-section flow velocity distribution model of the curvilinear channel based on a three-point parameter quantification method of a section segmentation theory on the basis of the initial curvilinear channel flow velocity distribution model; the three-point parameter quantification method based on the section segmentation theory obtains a full-section flow velocity distribution model of a curve channel, and specifically comprises the following steps: determining three normal lines corresponding to the section of the curve-type channel, and respectively acquiring the flow velocities of target points on the three normal lines; inputting the flow velocity of the target point into the initial curve type channel flow velocity distribution model to obtain corresponding quantitative undetermined parameters; and determining the full-section flow velocity distribution model of the curve type channel based on the quantitative undetermined parameter and the initial curve type channel flow velocity distribution model.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A curve type channel flow measuring method is characterized by comprising the following steps:
determining a flow velocity distribution model of a curve channel;
determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the parting line characteristics;
and dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method, and obtaining the flow information of the curve channel based on a speed area method.
2. The curvilinear channel flow measurement method of any one of claim 1, wherein the curvilinear channel flow velocity distribution model corresponds to an expression:
wherein: A. b, C is the correlation coefficient; u. of*The local friction flow rate; nu is the motion viscosity coefficient of water; l is the distance of the curve channel normal from each point on the curve channel normal to the side wall; klIs the vertical distance from the side wall to the dividing line.
3. The curvilinear channel flow measurement method of claim 1, wherein the normal average flow velocity on the curvilinear channel normal represents a point location correspondence expression:
wherein l1、l2Representing the point location for the normal average flow velocity,/1、l2The speed is the average speed of the speed measurement normal; klIs the distance from the curved channel bottom to the water surface line along the normal direction, and l is the normal average flow velocity representative point along KlTo the bottom of the curved trench.
4. The method for measuring the flow rate of the curvilinear channel according to claim 1, wherein the method for obtaining the flow rate information of the curvilinear channel based on the velocity area method by dividing the cross-sectional area of the curvilinear channel according to the position of the representative point of the normal average flow velocity by using a preset three-line-polygon division method specifically comprises:
dividing the area of the section of the curved channel according to the position of the representative point of the normal average flow velocity by using a preset three-line-polygon division method, and determining a plurality of partition polygons corresponding to the section of the curved channel;
calculating the flow of the curved channel according to the plurality of partition polygons by using a preset channel flow calculation formula to obtain the flow information of the curved channel;
the channel flow calculation formula is as follows:
Qgeneral assembly＝2×[α·u1·(S1+S2)]+β·u2·(S3+S4)
Wherein u is1Is the flow velocity, u, of the target point on the first and third normal lines2Is the flow velocity of the target point on the second normal; alpha and beta are flow correction coefficients; s1、S2、S3And S4The areas of the plurality of partition polygons are respectively.
5. The curvilinear channel flow measurement method of claim 4, wherein determining the normal average flow velocity representative point position of the curvilinear channel cross section comprises: determining three corresponding normal lines on the section of the curve-type channel; respectively determining corresponding normal average flow velocity representative point positions on the three normals;
the method for determining the plurality of partition polygons corresponding to the curved channel section by utilizing a preset three-line-polygon partition method and dividing the area of the curved channel section according to the position of the normal average flow velocity representative point specifically comprises the following steps: connecting the normal average flow velocity representative points on the three normals according to a preset mode to form a plurality of corresponding triangles; determining an intersection point of the outer center position of the triangle and the water surface and an intersection point of the outer center position of the triangle and the canal bottom according to the perpendicular bisectors corresponding to the three sides of the triangle; taking the outer center position as a partition boundary point, and connecting the partition boundary point to obtain a corresponding partition polygon; wherein, each partition polygon includes a normal average flow velocity representative point.
6. A method for obtaining a flow velocity distribution model of a curvilinear channel, comprising:
determining an initial curve type channel flow velocity distribution model according to the curve type channel;
on the basis of the initial curve channel flow velocity distribution model, a curve channel full-section flow velocity distribution model is obtained based on a three-point parameter quantification method of a section division theory;
the three-point parameter quantification method based on the section segmentation theory obtains a full-section flow velocity distribution model of a curve channel, and specifically comprises the following steps: determining three normal lines corresponding to the section of the curve-type channel, and respectively acquiring the flow velocities of target points on the three normal lines; inputting the flow velocity of the target point into the initial curve type channel flow velocity distribution model to obtain corresponding quantitative undetermined parameters; and determining the full-section flow velocity distribution model of the curve type channel based on the quantitative undetermined parameter and the initial curve type channel flow velocity distribution model.
7. A curved channel flow measurement device, comprising:
the algorithm model determining unit is used for determining a curve type channel flow velocity distribution model;
the representative point position determining unit is used for determining the normal average flow velocity representative point position of the section of the curve channel according to the curve channel flow velocity distribution model, the curve channel flow velocity partition and the partition line characteristics;
and the channel flow calculating unit is used for dividing the section area of the curve channel according to the position of the normal average flow velocity representative point by using a preset three-line-polygon division method and obtaining the flow information of the curve channel based on a speed area method.
8. An apparatus for obtaining a flow velocity distribution model of a curvilinear channel, comprising:
the initial flow velocity distribution model determining unit is used for determining an initial curve type channel flow velocity distribution model according to the curve type channel;
a full-section flow velocity distribution model determining unit, configured to obtain a full-section flow velocity distribution model of the curvilinear channel based on a three-point parameter quantification method of a section segmentation theory on the basis of the initial curvilinear channel flow velocity distribution model; the three-point parameter quantification method based on the section segmentation theory obtains a full-section flow velocity distribution model of a curve channel, and specifically comprises the following steps: determining three normal lines corresponding to the section of the curve-type channel, and respectively acquiring the flow velocities of target points on the three normal lines; inputting the flow velocity of the target point into the initial curve type channel flow velocity distribution model to obtain corresponding quantitative undetermined parameters; and determining the full-section flow velocity distribution model of the curve type channel based on the quantitative undetermined parameter and the initial curve type channel flow velocity distribution model.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of the curvilinear channel flow determination method of any of claims 1-5.
10. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the curvilinear channel flow determination method of any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110272485.XA CN113074786B (en) | 2021-03-12 | 2021-03-12 | Curve channel flow measuring method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110272485.XA CN113074786B (en) | 2021-03-12 | 2021-03-12 | Curve channel flow measuring method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113074786A true CN113074786A (en) | 2021-07-06 |
| CN113074786B CN113074786B (en) | 2022-05-06 |
Family
ID=76612724
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110272485.XA Active CN113074786B (en) | 2021-03-12 | 2021-03-12 | Curve channel flow measuring method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113074786B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114543903A (en) * | 2022-01-28 | 2022-05-27 | 水利部水土保持监测中心 | Portable torque type flow measuring device and flow measuring method |
| CN116755477A (en) * | 2023-08-16 | 2023-09-15 | 西安倍得新数据科技有限公司 | Automatic flow control and regulation method and system for fluid channel |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1559999A1 (en) * | 2004-01-30 | 2005-08-03 | Fuji Electric Systems Co., Ltd. | Wedge and wedge unit for use in ultrasonic doppler flow meter |
| JP2011122831A (en) * | 2009-12-08 | 2011-06-23 | Tokyo Electric Power Co Inc:The | Ultrasonic flow rate measurement method and ultrasonic flow rate measurement device |
| JP2017096773A (en) * | 2015-11-24 | 2017-06-01 | 株式会社Ｃｏｒｅ技術研究所 | Nondestructive inspection method and nondestructive inspection device and information specifying method in elastic wave tomography and information specifying device |
| CN107084765A (en) * | 2017-03-30 | 2017-08-22 | 中国农业大学 | A kind of canal capacity computational methods and device |
| CN108534845A (en) * | 2018-03-12 | 2018-09-14 | 江苏大渠自动化科技有限公司 | Open channel flow rate measurement device and assay method |
| CN109255185A (en) * | 2018-09-11 | 2019-01-22 | 中国水利水电科学研究院 | An a kind of two-dimentional hydrodynamics coupling analytical method based on urban surface underground pipe network |
| CN109916458A (en) * | 2019-04-12 | 2019-06-21 | 江苏亚楠电子科技有限公司 | A kind of decomposition cross-sectional flow method |
-
2021
- 2021-03-12 CN CN202110272485.XA patent/CN113074786B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1559999A1 (en) * | 2004-01-30 | 2005-08-03 | Fuji Electric Systems Co., Ltd. | Wedge and wedge unit for use in ultrasonic doppler flow meter |
| JP2011122831A (en) * | 2009-12-08 | 2011-06-23 | Tokyo Electric Power Co Inc:The | Ultrasonic flow rate measurement method and ultrasonic flow rate measurement device |
| JP2017096773A (en) * | 2015-11-24 | 2017-06-01 | 株式会社Ｃｏｒｅ技術研究所 | Nondestructive inspection method and nondestructive inspection device and information specifying method in elastic wave tomography and information specifying device |
| CN107084765A (en) * | 2017-03-30 | 2017-08-22 | 中国农业大学 | A kind of canal capacity computational methods and device |
| CN108534845A (en) * | 2018-03-12 | 2018-09-14 | 江苏大渠自动化科技有限公司 | Open channel flow rate measurement device and assay method |
| CN109255185A (en) * | 2018-09-11 | 2019-01-22 | 中国水利水电科学研究院 | An a kind of two-dimentional hydrodynamics coupling analytical method based on urban surface underground pipe network |
| CN109916458A (en) * | 2019-04-12 | 2019-06-21 | 江苏亚楠电子科技有限公司 | A kind of decomposition cross-sectional flow method |
Non-Patent Citations (2)
| Title |
|---|
| YU HAN，TONGSHU LI，SHIYU WANG，JIAN CHEN: "A Flow-Measuring Algorithm of Arc-Bottomed Open Channels through Multiple Characteristic Sensing Points of the Flow-Velocity Sensor in Argricultural Irrigation Areas", 《SENSORS》 * |
| 韩宇，曹文轩，李彤姝，吴章生，陈建: "圆弧底渠道测点布控方法研究", 《中国农村水利水电》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114543903A (en) * | 2022-01-28 | 2022-05-27 | 水利部水土保持监测中心 | Portable torque type flow measuring device and flow measuring method |
| CN114543903B (en) * | 2022-01-28 | 2022-11-04 | 水利部水土保持监测中心 | Portable torque type flow measuring device and flow measuring method |
| CN116755477A (en) * | 2023-08-16 | 2023-09-15 | 西安倍得新数据科技有限公司 | Automatic flow control and regulation method and system for fluid channel |
| CN116755477B (en) * | 2023-08-16 | 2023-11-03 | 西安倍得新数据科技有限公司 | Automatic flow control and regulation method and system for fluid channel |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113074786B (en) | 2022-05-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113074786B (en) | Curve channel flow measuring method and device | |
| Parsaie et al. | CFD modeling of flow pattern in spillway’s approach channel | |
| CN106503396A (en) | The multidimensional Hydraulic Power System transition analogy method coupled with finite volume method based on finite difference calculus | |
| Allen et al. | Characterizing the impact of geometric simplification on large woody debris using CFD | |
| CN105205841B (en) | The ground drawing generating method and system of GIS-Geographic Information System | |
| Andersson et al. | A study of the location of the entrance of a fishway in a regulated river with CFD and ADCP | |
| RAN et al. | Three-dimensional numerical simulation of flow in trapezoidal cutthroat flumes based on FLOW-3D | |
| KR101106548B1 (en) | Method for analyzing shallow water flow using the two-dimensional river flow model with tensor-type eddy viscosity | |
| Shan et al. | A simple method for estimating bed shear stress in smooth and vegetated compound channels | |
| Hodskinson | Computational fluid dynamics as a tool for investigating separated flow in river bends | |
| CN109190187A (en) | A kind of wind power plant wake flow calculation method and system based on more physical models | |
| Curulli et al. | Improved formulation for the geometric characteristics of the scour hole and deposition mound caused by a rotating propeller jet on a mobile bed | |
| Zhou et al. | Numerical modeling of 3D flow field among a compound stilling basin | |
| CN113742988B (en) | Centrifugal pump numerical simulation method and system based on gridless particle method | |
| Maddahi et al. | Determining minimum numbers of transects for accurate flow measurements using moving-vessel ADCPs | |
| JP5213212B2 (en) | Flow rate calculation system, flow rate calculation program, and flow rate calculation method | |
| CN114417473A (en) | Engineering arrangement method for optimizing fairway confluence area based on particle tracking | |
| Aricò et al. | Anisotropic potential of velocity fields in real fluids: Application to the MAST solution of shallow water equations | |
| Ferraro et al. | Unconfined propeller scour in waterways: The role of flow intensity | |
| Adamkowski et al. | The comparative analysis of the current-meter method and the pressure-time method used for discharge measurements in the Kaplan turbine penstocks | |
| CN109141542B (en) | Rectangular channel non-groove water metering method based on critical water depth groove working principle | |
| Jin et al. | Hydraulic characteristics and flow measurement performance analysis of portable fish-shaped flumes in U-shaped channels | |
| Yu et al. | Experimental investigation on interaction modes and geometric features of two ventilated cavities | |
| Krzemianowski et al. | Use of non-uniform rational b-splines for discharge calculation in the velocity area method | |
| Bizjan et al. | Flow image velocimetry method based on advection-diffusion equation |
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 |