CN107084765A - A kind of canal capacity computational methods and device - Google Patents
A kind of canal capacity computational methods and device Download PDFInfo
- Publication number
- CN107084765A CN107084765A CN201710203894.8A CN201710203894A CN107084765A CN 107084765 A CN107084765 A CN 107084765A CN 201710203894 A CN201710203894 A CN 201710203894A CN 107084765 A CN107084765 A CN 107084765A
- Authority
- CN
- China
- Prior art keywords
- section
- flow
- channel
- canal
- cross
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
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/76—Devices for measuring mass flow of a fluid or a fluent solid material
- G01F1/86—Indirect mass flowmeters, e.g. measuring volume flow and density, temperature or pressure
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The present invention provides a kind of canal capacity computational methods and device, and methods described includes：The flow velocity of default each flow velocity measurement point, determines the velocity flow profile factor of the canal cross section in distance and the canal cross section based on channel water level to channel bottom；The velocity flow profile factor based on the canal cross section, determines the flow speed value of each position point on the canal cross section；Based on the cross-sectional area of the flow speed value of each position point and the channel on the canal cross section, the flow of the channel is calculated.A kind of canal capacity computational methods proposed by the present invention and device, by determining the velocity flow profile factor of canal cross section in channel, to calculate flow velocity in channel, and then calculate canal capacity, calculating process lifts computational accuracy independent of channel environmental factor.
Description
Technical field
The present invention relates to hydraulic engineering field, more particularly, to a kind of canal capacity computational methods and device.
Background technology
At present, it is to realize the accurate water distribution in irrigated area, scientific dispatch and information system management to be measured based on final stage Flow Regulation in Irrigation Canal Systems
Key.
Final canal system is generally trapezoid canal with arc-bottom, and velocity flow profile will be different from traditional rectangular channel or trapezoidal open channel, and then
So that flow, which is inquired into, is different from traditional rectangular or trapezoidal open channel.The change of flow state of current is various in the open channel of irrigated area, therefore inquires into logical
Flow relocity calculation method is very necessary.In the prior art, it is general to useCalculate canal
Road flow, wherein Q are that flow, A are that channel cross-sectional area, V are that channel flow velocity, C are gool arc bottom radius, i to thank ability coefficient, R
It is the disconnected roughness in channel side for channel base slope and n, wherein the ability coefficient of thanking breaks coarse with section configuration, size and channel side
Rate is relevant.
Existing canal capacity computational methods need the influence of strict control environmental factor when calculating flow, and irrigated area water
The live change of flow state of stream is various, can not obtain the data of environmental factor in real time, this strong influence canal capacity
Computational accuracy.
The content of the invention
The present invention provides a kind of canal capacity calculating side for overcoming above mentioned problem or solving the above problems at least in part
Method and device.
According to an aspect of the present invention there is provided a kind of canal capacity computational methods, including：
The flow velocity of default each flow velocity measurement point in distance and canal cross section based on channel water level to channel bottom, really
The velocity flow profile factor of the fixed canal cross section；
The velocity flow profile factor based on the canal cross section, determines the flow speed value of each position point on the canal cross section；
Based on the cross-sectional area of the flow speed value of each position point and the channel on the canal cross section, the channel is calculated
Flow.
Wherein, default each flow velocity measurement point in the distance and canal cross section based on channel water level to channel bottom
Flow velocity, determine the velocity flow profile factor of the canal cross section, including：
First flow velocity distribution factor is determined by following formula：
α1=(gSh-u* 2)/u* 2
Wherein, α1It is that acceleration of gravity, S are that the gradient of channel, h are the channel for the first flow velocity distribution factor, g
Distance of the water level to channel bottom, u*For the local resistance flow velocity on canal cross section；
Second flow speed distribution factor is determined by following formula：
Wherein, α2For longitudinal and vertical flow velocity that the second flow speed distribution factor, u and v are each flow velocity measurement point；
The velocity flow profile factor is determined by following formula：
α=α1+α2
Wherein, α is the velocity flow profile factor.
Wherein, the velocity flow profile factor based on the canal cross section, determines each position point on the canal cross section
Velocity flow profile, including：
The flow speed value of each position point on the canal cross section is determined by following formula：
Wherein, κ is that karman constant, y are each position point position, y on the canal cross section0It is cross section correct constant, u with y
The flow speed value for being distributed as each position point on the canal cross section.
Wherein, the cross-sectional area based on the flow speed value of each position point and the channel on the canal cross section, meter
The flow of the channel is calculated, including：
Calculate flow speed value of each position point along canal cross section normal direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section normal direction it is disconnected
Face area, calculates the flow of the channel.
Wherein, the cross-sectional area based on the flow speed value of each position point and the channel on the canal cross section, meter
The flow of the channel is calculated, including：
Calculate flow speed value of each position point along canal cross section vertical line direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section vertical line direction it is disconnected
Face area, calculates the flow of the channel.
On the other hand, the present invention provides a kind of canal capacity computing device, including：
First determining module, it is default in the distance based on channel water level to channel bottom and the canal cross section
The flow velocity of each flow velocity measurement point, determines the velocity flow profile factor of the canal cross section；
Second determining module, for the velocity flow profile factor based on the canal cross section, is determined each on the canal cross section
The flow speed value of location point；
Computing module, for the section face based on the flow speed value of each position point and the channel on the canal cross section
Product, calculates the flow of the channel.
Wherein, the first determining module specifically for：
First flow velocity distribution factor is determined by following formula：
α1=(gSh-u* 2)/u* 2
Wherein, α1It is that acceleration of gravity, S are that the gradient of channel, h are the channel for the first flow velocity distribution factor, g
Distance of the water level to channel bottom, u*For the local resistance flow velocity on canal cross section；
Second flow speed distribution factor is determined by following formula：
Wherein, α2For longitudinal and vertical flow velocity that the second flow speed distribution factor, u and v are each flow velocity measurement point；
The velocity flow profile factor is determined by following formula：
α=α1+α2
Wherein, α is the velocity flow profile factor.
Wherein, the second determining module specifically for：
The flow speed value of each position point on the canal cross section is determined by following formula：
Wherein, κ is that karman constant, y are each position point position, y on the canal cross section0It is cross section correct constant, u with y
The flow speed value for being distributed as each position point on the canal cross section.
Wherein, the computing module specifically for：
Calculate flow speed value of each position point along canal cross section normal direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section normal direction it is disconnected
Face area, calculates the flow of the channel.
Wherein, the computing module specifically for：
Calculate flow speed value of each position point along canal cross section vertical line direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section vertical line direction it is disconnected
Face area, calculates the flow of the channel.
A kind of canal capacity computational methods proposed by the present invention and device, by determine the velocity flow profile of canal cross section because
Son, to calculate flow velocity in channel, and then calculates canal capacity, calculating process is independent of channel environmental factor, lifting
Computational accuracy.
Brief description of the drawings
Fig. 1 is a kind of canal capacity computational methods flow chart provided in an embodiment of the present invention；
Fig. 2 is a kind of canal capacity computing device structure figure provided in an embodiment of the present invention；
Embodiment
With reference to the accompanying drawings and examples, the embodiment to the present invention is described in further detail.Implement below
Example is used to illustrate the present invention, but is not limited to the scope of the present invention.
Fig. 1 is a kind of canal capacity computational methods flow chart provided in an embodiment of the present invention, including S1 to S3：
Default each flow velocity measurement point on S1, the distance based on channel water level to channel bottom and the canal cross section
Flow velocity, determines the velocity flow profile factor of the canal cross section；
S2, the velocity flow profile factor based on channel face absolutely, determine the flow velocity of each position point on the canal cross section
Value；
S3, the cross-sectional area based on the flow speed value of each position point and the channel on the canal cross section, calculate described
The flow of channel.
In S1, the channel is generally trapezoid canal with arc-bottom, it is to be understood that trapezoid canal with arc-bottom is final canal system
Typical Representative, the present invention is not specifically limited to channel type.
In S1, the layout rules embodiment of the present invention of default each flow velocity measurement point is specifically limited not to this, but
Arrange that quantity can not typically arrange more than 30 measurement points very little.
In S1, the velocity flow profile factor is the parameter of influence velocity flow profile, and the determination of the velocity flow profile factor is based on section
The determination of drag velocity and the bursting phenomenon of vertical flow velocity.
In S2, the flow speed value of each position point is the flow speed value of any location point on canal cross section on the canal cross section, can
With understanding, each position point described herein includes the default flow velocity measurement point described in S1.
In S3, the flow for calculating channel is that the mean flow rate in face multiplies absolutely according to the cross-sectional area and channel of channel
Product try to achieve, the channel absolutely face mean flow rate be according on channel absolutely face in a direction all location point flow velocitys be averaged
Gained.It should be noted that a direction refers to each measurement point along canal cross section vertical line or normal side on the canal cross section
To flow velocity.
Specifically, the velocity distribution experiment of research channel, builds the tank of a Semi round bottom, the upstream water leg that breaks is symmetrical
With the center line of channel so that the current injection condition at feeder connection is consistent.Electronic flowmeter is fixed on pipeline,
For monitoring canal capacity.
64 flow velocity measurement points are set in canal cross section, and 64 flow velocity measurement points are measured under conditions of h/D=0.333
The time average velocity at place, the error that measured value calculates formula calculated value with velocity flow profile is analyzed along normal direction, and vertically simultaneously
Orientation measured value calculates the error of formula calculated value with velocity flow profile, as a result as shown in table 1.
Velocity flow profile error table of comparisons during 1 h/D=0.333 of table
Wherein, h is longitudinally height, arc bottom diameter that D be channel bottom of the channel water level to channel bottom.As shown in table 1,
When being calculated along canal cross section normal direction, the mean error of measured value and actual value is that error is small in 2.6%, 64 measurement points
90.8% is accounted in 5%, and overall data is not above 8% error dot and meets computational accuracy, along canal cross section vertical line
During direction calculating, the mean error of measured value and actual value is 4.8%, and the point that error is less than 5% accounts for 46.9%, meets and calculates essence
Degree, thus the accuracy rate of the provable flow rate calculation of the embodiment of the present invention of experimental data is improved.
Thus the accuracy rate of the provable flow relocity calculation of the embodiment of the present invention of experimental data is improved.
The embodiment of the present invention is by determining the velocity flow profile factor of canal cross section, to calculate flow velocity in channel, and
Comparing calculation is carried out from along canal cross section vertical line and normal direction both direction, canal capacity can be accurately calculated, accuracy rate is obtained
Lifting is arrived.
On the basis of embodiment described in Fig. 1, the distance and channel based on channel water level to channel bottom are absolutely
The flow velocity of default each flow velocity measurement point on face, determines the velocity flow profile factor of the canal cross section, including：
First flow velocity distribution factor is determined by following formula：
α1=(gSh-u* 2)/u* 2
Wherein, α1It is that acceleration of gravity, S are that the gradient of channel, h are the channel for the first flow velocity distribution factor, g
Water level is to the local resistance flow velocity that the distance of channel bottom, u* are on the canal cross section；
Second flow speed distribution factor is determined by following formula：
Wherein, α2For longitudinal and vertical flow velocity that the second flow speed distribution factor, u and v are each flow velocity measurement point；
The velocity flow profile factor is determined by following formula：
α=α1+α2
Wherein, α is the velocity flow profile factor.
Specifically, the embodiment of the present invention is entered by proposing velocity flow profile factor-alpha to flow rate of water flow distribution scenario in channel
Row explanation, but not limited to this title.Velocity flow profile factor-alpha is influences the parameter of velocity flow profile, by the bursting phenomenon of vertical flow velocity
With section drag velocity obtained by, specifically, the velocity flow profile factor can be divided into first-class fast factor-alpha1With second flow speed distribution because
Sub- α2, when current flow downward, because there is edge down phenomenon, α=α when current flow downward in flow velocity1+α2＞ 0, α1
＞ 0,Similarly understand, when current are flowed up, α=α1+α2＜ 0, α1＞ 0 ,-α2＞ α1＞ 0.
The embodiment of the present invention is calculated by calculating the velocity flow profile factor come the overall velocity flow profile to current in channel,
It need not obtain in the calculation on environmental factors such as canal cross section materials, computational methods are general, convenient and accuracy rate is higher.
On the basis of above-described embodiment, the velocity flow profile factor based on the canal cross section determines the channel
The velocity flow profile of each position point on section, including：
The flow speed value of each position point on the canal cross section is determined by following formula：
Wherein, κ is that karman constant, y are channel each position point position, y on face absolutely0For cross section correct constant, u with
The y flow speed value for being distributed as each position point on the canal cross section.
Specifically, main inquire into according to Reynolds average equation group progress flow velocity under the conditions of smooth water flow, the Reynolds average
Equation group such as formula (1) and formula (2)：
Wherein u, v, w represent x, y respectively, and the mean flow rate in z directions, g is acceleration of gravity, and S is the gradient, and μ cuts for Reynolds
Tangent tensor part.For full-blown uniform flow, (2) formula can be written as：
Formula (3) can be further written as：
The viscosity that wherein ρ is fluid density, μ is fluid；U', v', w' represent x, and y, z directions Turbulence Flow is fluctuated,
In two dimension flow, formula (4) can be reduced to：
Formula (7) integration is obtained：
In formula (8),For drag velocity, h be channel water level to the height of channel bottom, R is hydraulic radius,
The equation of momentum can be written as：
Formula (9) is substituted into Boussinesq to assume to obtain：
According to eddy viscosity formula：
Wherein, νtRepresent eddy viscosity, κ and represent karman constant.
Formula (10) and formula (11) are substituted into formula (9) to obtain：
Formula (12) is integrated and can obtained：
In formula (13), yo=ν/(cu*(z)), becauseFormula (13) can be reduced to：
In formula (14), α is that the velocity flow profile factor, κ are that karman constant, y are each position point position on the canal cross section
Put, y0For cross section correct constant, u with y the flow speed value for being distributed as each position point on the canal cross section.
It is disconnected that the embodiment of the present invention calculates that the point flow velocity of the optional position in current is determined by substituting into the velocity flow profile factor
It is concise in face mean flow rate, calculating process, it is possible to increase computational efficiency.
It is described based on the flow speed value of each position point and institute on the canal cross section on the basis of embodiment described in Fig. 1
The cross-sectional area of channel is stated, the flow of the channel is calculated, including：
Calculate flow speed value of each position point along canal cross section normal direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section normal direction it is disconnected
Face area, calculates the flow of the channel.
It should be noted that the embodiment of the present invention is flowed by determining in channel it can be seen from the table 1 in Fig. 1 embodiments
The velocity flow profile factor of fast measurement point, to calculate flow velocity in channel, and selection is counted along canal cross section normal direction
Calculate, can accurately calculate channel flow velocity, accuracy rate is improved.In order to prove the conclusion, the embodiment of the present invention changes experiment bar
Part, is compared to channel flow velocity further experiment.
Specifically, the embodiment of the present invention is from 112 flow velocity measurement points, and along normal side under conditions of h/D=0.500
The error of formula calculated value is calculated with velocity flow profile to analysis measured value, and vertically Orientation measured value divides with flow velocity simultaneously
The error of cloth calculation formula calculated value, experimental result is as shown in table 2.
Velocity flow profile error table of comparisons during 2 h/D=0.500 of table
As shown in table 2, it is 2.9% along the mean error of canal cross section method line computation in h/D=0.500, and along channel
The mean error of section vertical line direction calculating then reaches 11.8%.So preferred, the embodiment of the present invention is used along canal cross section
It is relatively reasonable that normal orientation calculates flow velocity.
The embodiment of the present invention is used along canal cross section normal orientation and calculates canal capacity, and the calculating side can be obtained through experimental data
Method accuracy rate is high.
On the basis of embodiment described in Fig. 1, it is described based on the velocity flow profile of each position point on the canal cross section and
The cross-sectional area of the channel, calculates the flow of the channel, including：
Calculate flow speed value of each position point along canal cross section vertical line direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section vertical line direction it is disconnected
Face area, calculates the flow of the channel.
It should be noted that the result calculated along canal cross section normal direction is better than along along canal cross section vertical line direction
Result of calculation, there has also been larger along the calculating in canal cross section vertical line direction relative to existing method of calculating flux in accuracy rate
Lifting, is the experimental result for the h/D=0.4 for employing 100 data measurement points as shown in table 3.
Velocity flow profile error table of comparisons during 3 h/D=0.4 of table
As shown in table 3, it is 6.8% along the mean error of canal cross section vertical line direction calculating, this error is relative to existing
Technology is also obviously improved, although not as the calculating accuracy rate height along canal cross section normal direction, but also meet computational accuracy.
Fig. 2 is a kind of canal capacity computing device structure figure provided in an embodiment of the present invention, including：
First determining module 1, the second determining module 2 and computing module 3, wherein：
First determining module, it is default in the distance based on channel water level to channel bottom and the canal cross section
The flow velocity of each flow velocity measurement point, determines the velocity flow profile factor of the canal cross section；
Second determining module, for the velocity flow profile factor based on the canal cross section, is determined each on the canal cross section
The flow speed value of location point；
Computing module, for the section face based on the flow speed value of each position point and the channel on the canal cross section
Product, calculates the flow of the channel.
Specifically, the bursting phenomenon of determination and vertical flow velocity of first determining module 1 based on section drag velocity.
Second determining module 2 determines that the flow speed value of each position point on the canal cross section is any location point on canal cross section
Velocity flow profile situation, it is to be understood that each position point described herein includes the default flow velocity measurement point described in S1.
The flow that computing module 3 calculates the calculating channel is the mean flow of cross-sectional area according to channel and canal cross section
The product of speed is tried to achieve, and the mean flow rate of the canal cross section is according to all location point flow velocitys are put down in a direction on canal cross section
Obtained by.It should be noted that a direction refers to each measurement point along canal cross section vertical line or normal on the canal cross section
The flow velocity in direction.
The embodiment of the present invention determines the velocity flow profile factor of canal cross section by the first determining module, and computing module comes to canal
Flow velocity is calculated in road, and carries out comparing calculation from along canal cross section vertical line and normal direction both direction, can accurately be counted
Canal capacity is calculated, accuracy rate is improved.
On the basis of embodiment described in Fig. 2, the first determining module specifically for：
First flow velocity distribution factor is determined by following formula：
α1=(gSh-u* 2)/u* 2
Wherein, α1It is that acceleration of gravity, S are that the gradient of channel, h are the channel for the first flow velocity distribution factor, g
Distance of the water level to channel bottom, u*For the local resistance flow velocity on the canal cross section；
Second flow speed distribution factor is determined by following formula：
Wherein, α2For longitudinal and vertical flow velocity that the second flow speed distribution factor, u and v are each flow velocity measurement point；
The velocity flow profile factor is determined by following formula：
α=α1+α2
Wherein, α is the velocity flow profile factor.
How specific first determining module determines that the velocity flow profile factor can be found in above-described embodiment, and here is omitted.
On the basis of above-described embodiment, the second determining module specifically for：
The flow speed value of each position point on the canal cross section is determined by following formula：
Wherein, κ is that karman constant, y are each position point position, y on the canal cross section0It is cross section correct constant, u with y
The flow speed value for being distributed as each position point on the canal cross section.
How specific second determining module determines that velocity flow profile can be found in above-described embodiment, and here is omitted.
On the basis of embodiment described in Fig. 2, the computing module specifically for：
Calculate flow speed value of each position point along canal cross section normal direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section normal direction it is disconnected
Face area, calculates the flow of the channel.
It is preferred that, computing module is used to be calculated along canal cross section normal direction, table 1, the and of table 2 provided by above-described embodiment
Data must can be demonstrate,proved in table 3.
The embodiment of the present invention calculates canal capacity by computing module along canal cross section normal orientation, improves flow rate calculation
Accuracy rate.
On the basis of embodiment described in Fig. 2, the computing module specifically for：
Calculate flow speed value of each position point along canal cross section vertical line direction on the canal cross section；
Flow speed value and the channel based on each position point on the canal cross section along canal cross section vertical line direction it is disconnected
Face area, calculates the flow of the channel.
Although computing module uses the result along canal cross section vertical line direction calculating not as along canal cross section normal direction
Result of calculation is accurate, but calculates accuracy rate and also improved a lot than prior art.
It is preferred that, flow velocity and the channel based on each position point on the canal cross section along canal cross section normal direction
Cross-sectional area, the flow for calculating the channel is contrasted with electromagnetic flowmeter flow, under three kinds of operating modes, h/D=0.333,
When 0.4 and 0.5, mean error is respectively 5.5%, 6.5% and 4.7%.The results show embodiment of the present invention is used along canal
The computational methods accuracy rate of road section normal direction direction calculating canal capacity is high.
A kind of canal capacity computational methods proposed by the present invention and device, by the flow velocity point for determining canal cross section in channel
The cloth factor, to calculate flow velocity in channel, and then calculates mean velocity in section and canal capacity, calculating process independent of
Channel environmental factor, lifts computational accuracy.
Finally, the present processes are only preferably embodiment, are not intended to limit the scope of the present invention.It is all
Within the spirit and principles in the present invention, any modification, equivalent substitution and improvements made etc. should be included in the protection of the present invention
Within the scope of.
Claims (10)
1. a kind of canal capacity computational methods, it is characterised in that including：
The flow velocity of default each flow velocity measurement point in distance and the canal cross section based on channel water level to channel bottom, really
The velocity flow profile factor of the fixed canal cross section；
The velocity flow profile factor based on the canal cross section, determines the flow speed value of each position point on the canal cross section；
Based on the cross-sectional area of the flow speed value of each position point and the channel on the canal cross section, the stream of the channel is calculated
Amount.
2. according to the method described in claim 1, it is characterised in that the distance based on channel water level to channel bottom and
The flow velocity of default each flow velocity measurement point on canal cross section, determines the velocity flow profile factor of the canal cross section, including：
First flow velocity distribution factor is determined by following formula：
α1=(gSh-u* 2)/u* 2
Wherein, α1It is that acceleration of gravity, S are that the gradient of channel, h are the channel water level for the first flow velocity distribution factor, g
Distance, u to channel bottom*For the local resistance flow velocity on the canal cross section；
Second flow speed distribution factor is determined by following formula：
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>u</mi>
<mi>v</mi>
</mrow>
<mrow>
<msup>
<msub>
<mi>u</mi>
<mo>*</mo>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
Wherein, α2For longitudinal and vertical flow velocity that the second flow speed distribution factor, u and v are each flow velocity measurement point；
The velocity flow profile factor is determined by following formula：
α=α1+α2
Wherein, α is the velocity flow profile factor.
3. method according to claim 2, it is characterised in that the velocity flow profile factor based on the canal cross section,
The flow speed value of each position point on the canal cross section is determined, including：
The flow speed value of each position point on the canal cross section is determined by following formula：
<mrow>
<mfrac>
<mi>u</mi>
<msub>
<mi>u</mi>
<mo>*</mo>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>&kappa;</mi>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>y</mi>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>&alpha;</mi>
<mi>&kappa;</mi>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>y</mi>
<mi>h</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
Wherein, κ is that karman constant, y are each position point position, y on the canal cross section0For cross section correct constant, u with y distribution
For the flow speed value of each position point on the canal cross section.
4. according to the method described in claim 1, it is characterised in that the flow velocity based on each position point on the canal cross section
The cross-sectional area of value and the channel, calculates the flow of the channel, including：
Calculate flow speed value of each position point along canal cross section normal direction on the canal cross section；
The section face of flow speed value and the channel based on each position point on the canal cross section along canal cross section normal direction
Product, calculates the flow of the channel.
5. according to the method described in claim 1, it is characterised in that the flow velocity based on each position point on the canal cross section
The cross-sectional area of value and the channel, calculates the flow of the channel, including：
Calculate flow speed value of each position point along canal cross section vertical line direction on the canal cross section；
The section face of flow speed value and the channel based on each position point on the canal cross section along canal cross section vertical line direction
Product, calculates the flow of the channel.
6. a kind of canal capacity computing device, it is characterised in that including：
First determining module, default each stream in the distance based on channel water level to channel bottom and the canal cross section
The flow velocity of fast measurement point, determines the velocity flow profile factor of the canal cross section；
Second determining module, for the velocity flow profile factor based on the canal cross section, determines each position on the canal cross section
The flow speed value of point；
Computing module, for the cross-sectional area based on the flow speed value of each position point and the channel on the canal cross section, meter
Calculate the flow of the channel.
7. device according to claim 6, it is characterised in that the first determining module specifically for：
First flow velocity distribution factor is determined by following formula：
α1=(gSh-u* 2)/u* 2
Wherein, α1It is that acceleration of gravity, S are that the gradient of channel, h are the channel water level for the first flow velocity distribution factor, g
Distance, u to channel bottom*For the local resistance flow velocity on canal cross section；
Second flow speed distribution factor is determined by following formula：
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>u</mi>
<mi>v</mi>
</mrow>
<mrow>
<msup>
<msub>
<mi>u</mi>
<mo>*</mo>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
Wherein, α2For longitudinal and vertical flow velocity that the second flow speed distribution factor, u and v are each flow velocity measurement point；
The velocity flow profile factor is determined by following formula：
α=α1+α2
Wherein, α is the velocity flow profile factor.
8. device according to claim 7, it is characterised in that the second determining module specifically for：
The flow speed value of each position point on the canal cross section is determined by following formula：
<mrow>
<mfrac>
<mi>u</mi>
<msub>
<mi>u</mi>
<mo>*</mo>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>&kappa;</mi>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>y</mi>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>&alpha;</mi>
<mi>&kappa;</mi>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>y</mi>
<mi>h</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
Wherein, κ is that karman constant, y are each position point position, y on the canal cross section0For cross section correct constant, u with y distribution
For the flow speed value of each position point on the canal cross section.
9. device according to claim 6, it is characterised in that the computing module specifically for：
Calculate velocity flow profile of each position point along canal cross section normal direction on the canal cross section；
The section face of flow speed value and the channel based on each position point on the canal cross section along canal cross section normal direction
Product, calculates the flow of the channel.
10. device according to claim 6, it is characterised in that the computing module specifically for：
Calculate flow speed value of each position point along canal cross section vertical line direction on the canal cross section；
The section face of flow speed value and the channel based on each position point on the canal cross section along canal cross section vertical line direction
Product, calculates the flow of the channel.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710203894.8A CN107084765A (en) | 2017-03-30 | 2017-03-30 | A kind of canal capacity computational methods and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710203894.8A CN107084765A (en) | 2017-03-30 | 2017-03-30 | A kind of canal capacity computational methods and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107084765A true CN107084765A (en) | 2017-08-22 |
Family
ID=59614468
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710203894.8A Pending CN107084765A (en) | 2017-03-30 | 2017-03-30 | A kind of canal capacity computational methods and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107084765A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108319745A (en) * | 2017-12-18 | 2018-07-24 | 中国水利水电科学研究院 | Channel unsteady flow computational methods and device |
| CN108491561A (en) * | 2018-01-26 | 2018-09-04 | 中国水利水电科学研究院 | A kind of gates of segmental shape based on flow regime identification crosses flow calculation methodologies |
| CN108519129A (en) * | 2018-04-25 | 2018-09-11 | 中国灌溉排水发展中心 | Canal capacity monitoring method and system |
| CN110568213A (en) * | 2019-07-10 | 2019-12-13 | 湖南省计量检测研究院 | Canal water flow velocity calculation method and canal |
| CN111524035A (en) * | 2020-05-13 | 2020-08-11 | 水利部交通运输部国家能源局南京水利科学研究院 | Hydrodynamic regulation and control threshold determination method for improvement of river network water environment in plain city |
| CN113074786A (en) * | 2021-03-12 | 2021-07-06 | 中国农业大学 | Curve channel flow measuring method and device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5467650A (en) * | 1990-02-12 | 1995-11-21 | Cushing; Vincent J. | Open channel flowmeter |
| CN104535123A (en) * | 2014-12-29 | 2015-04-22 | 太原理工大学 | Online automatic channel flow measuring method |
| CN105157759A (en) * | 2015-05-13 | 2015-12-16 | 柴元方 | Method for measuring U-shaped open channel flow based on Thiessen polygon multi-point method |
| CN105841753A (en) * | 2016-05-26 | 2016-08-10 | 中国水利水电科学研究院 | Channel flow acquiring method |
-
2017
- 2017-03-30 CN CN201710203894.8A patent/CN107084765A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5467650A (en) * | 1990-02-12 | 1995-11-21 | Cushing; Vincent J. | Open channel flowmeter |
| CN104535123A (en) * | 2014-12-29 | 2015-04-22 | 太原理工大学 | Online automatic channel flow measuring method |
| CN105157759A (en) * | 2015-05-13 | 2015-12-16 | 柴元方 | Method for measuring U-shaped open channel flow based on Thiessen polygon multi-point method |
| CN105841753A (en) * | 2016-05-26 | 2016-08-10 | 中国水利水电科学研究院 | Channel flow acquiring method |
Non-Patent Citations (1)
| Title |
|---|
| YANG SQ等: "Velocity Distribution and Dip-Phenomenon in Smooth Uniform Open Channel Flows", 《JOURNAL OF HYDRAULIC ENGINEERING》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108319745A (en) * | 2017-12-18 | 2018-07-24 | 中国水利水电科学研究院 | Channel unsteady flow computational methods and device |
| CN108319745B (en) * | 2017-12-18 | 2020-12-08 | 中国水利水电科学研究院 | Channel unsteady flow calculation method and device |
| CN108491561A (en) * | 2018-01-26 | 2018-09-04 | 中国水利水电科学研究院 | A kind of gates of segmental shape based on flow regime identification crosses flow calculation methodologies |
| CN108491561B (en) * | 2018-01-26 | 2021-05-11 | 中国水利水电科学研究院 | Radial gate over-flow calculation method based on flow state identification |
| CN108519129A (en) * | 2018-04-25 | 2018-09-11 | 中国灌溉排水发展中心 | Canal capacity monitoring method and system |
| CN108519129B (en) * | 2018-04-25 | 2024-01-19 | 中国灌溉排水发展中心 | Channel flow monitoring method and system |
| CN110568213A (en) * | 2019-07-10 | 2019-12-13 | 湖南省计量检测研究院 | Canal water flow velocity calculation method and canal |
| CN110568213B (en) * | 2019-07-10 | 2021-09-03 | 湖南省计量检测研究院 | Canal water flow velocity calculation method and canal |
| CN111524035A (en) * | 2020-05-13 | 2020-08-11 | 水利部交通运输部国家能源局南京水利科学研究院 | Hydrodynamic regulation and control threshold determination method for improvement of river network water environment in plain city |
| CN113074786A (en) * | 2021-03-12 | 2021-07-06 | 中国农业大学 | Curve channel flow measuring method and device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107084765A (en) | A kind of canal capacity computational methods and device | |
| Raudkivi | Study of sediment ripple formation | |
| CN108254032A (en) | River ultrasonic wave time difference method method of calculating flux | |
| Straub | Effect of channel‐contraction works upon regimen of movable bed‐streams | |
| CN111021304B (en) | Dam breaking test composite simulation device and method based on real-time reservoir capacity adjustment | |
| Liu et al. | Impact characteristics of free over-fall in pool zone with upstream bed slope | |
| CN108104052B (en) | Tidal reach tributary inlet door ejectment water hinge entirety fluid flowing model test method | |
| CN113705127A (en) | Dam break risk prediction and emergency response method and device for tailing pond and electronic equipment | |
| CN111044410A (en) | Device and method for detecting rheological property of coal slurry based on safety ring pipe | |
| CN101699264B (en) | Testing device and testing method for emission reduction of maintenance liquid in pipeline | |
| Salazar-Campoy et al. | A physical model to study the effects of nozzle design on dense two-phase flows in a slab mold casting Ultra-Low carbon steels | |
| Tullis | Size-scale effects of labyrinth weir hydraulics | |
| CN107167188A (en) | A kind of method applied to rectangular open channel flow measurement | |
| CN209296060U (en) | A kind of rectangle of measuring automatically flow is without larynx flume | |
| CN107090798A (en) | Restoration ＆ saturation coefficient computational methods under a kind of reservoir over-saturation sediment transport state | |
| DE102008035349A1 (en) | Plant and method for dispensing liquid from a multi-chamber tank car essentially by gravity | |
| CN103852113B (en) | A kind of flume and its metering method suitable for U-shaped channel | |
| CN110186649A (en) | Ice water two-phase flow water delivery experimental rig and its test method | |
| CN105444828A (en) | Design method for rectifier of ultrasonic flowmeter | |
| Bailey | Development and testing of experimental equipment to measure pore pressure and dynamic pressure at points outside a pipe leak | |
| CN202372236U (en) | Constant water head water supply system of flow measuring device | |
| Arora et al. | Resistance to Flow and Velocity Distribution in Rigid Boundary Channels Carrying Sediment‐Laden Flow | |
| Ramamurthy et al. | End depth–discharge relation at free overfall of trapezoidal channels | |
| Roshan et al. | The effect of reservoir geometry on the critical submergence depth in hydroelectric power plants intake | |
| CN111121854A (en) | Device for measuring energy dissipation rate of trajectory jet energy dissipater and using method |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170822 |
|
| RJ01 | Rejection of invention patent application after publication |