〇、闸门出流流态判别流程
不同时段内闸门因闸上水位、闸下水位和闸门开度的相互关系不同而会出现闸孔自由出流、闸孔淹没出流、堰自由出流、堰淹没出流等不同的流态。
对于同一种流态,根据闸下水位是否能够淹没跃后水深可分为两种不同流态。对于平面闸门,通常使用闸门的相对开度(即开度e和上游水位hu的比值e/hu)判断出流形式为堰流或着闸孔出流,闸孔出流流态下自由出流和淹没出流的判别标准为闸下水位hd和跃后水深hz’的大小关系,堰流流态下自由出流和淹没出流的判别标准为闸上水位hu和闸下水位hd的比值,因此,闸门出流流态的判别可按下图所示流程图进行判别1:
具体可描述为:若相对开度小于0.1,则认为闸门处于关闭状态;若相对开度处于0.1-0.65之间,则通过闸门的流态为闸孔出流,根据闸下水位与跃后水深的大小关系,如果闸下水位小于跃后水深,则当前闸门流态为闸孔自由出流,否则当前闸门流态为闸孔淹没出流;当相对开度大于0.65时,则通过闸门的流态为宽顶堰出流;根据闸下水位与闸上水位的比值,若该比值处于0-0.8之间,则通过闸门的流态为宽顶堰自由出流,若比值大于0.8,则以0.9作为宽顶堰低淹没出流和高淹没出流的分界值。
作者简介
很高兴认识您!
我叫卢家波,河海大学水文学及水资源博士研究生,研究兴趣为高效洪水淹没预测、洪水灾害预警、机器学习、替代模型和降阶模型。
变化环境下,极端洪水事件多发,我希望能通过研究为水灾害防御做出贡献,为人民服务。
欢迎交流讨论和研究合作,vx Jiabo_Lu。
主页 https://lujiabo98.github.io
简历 https://lujiabo98.github.io/file/CV_JiaboLu_zh.pdf
博客 https://blog.csdn.net/weixin_43012724?type=blog
来信请说明博客标题及链接,谢谢。
一、传统过闸流量公式
1.1 武水公式
1.1.1 闸孔出流公式
Q
=
σ
m
b
e
2
g
H
0
Q=\sigma m be \sqrt{2gH_0}
Q=σmbe2gH0
式中:Q为过闸流量;
σ
\sigma
σ为淹没系数;m为自由孔流综合流量系数;b为过水断面宽度;e为闸门开启高度;g为重力加速度;
H
0
H_0
H0为计入行近流速水头的堰闸前水深。
1.1.2 自由孔流综合流量系数
(1)弧形闸门闸孔自由出流流量系数2:
m
=
(
0.97
−
0.81
θ
180
°
)
−
(
0.56
−
0.81
θ
180
°
)
e
H
m=(0.97-0.81\frac{\theta}{180\degree})-(0.56-0.81\frac{\theta}{180\degree})\frac{e}{H}
m=(0.97−0.81180°θ)−(0.56−0.81180°θ)He
应用范围:
25
°
<
θ
≤
90
°
25\degree<\theta\le90\degree
25°<θ≤90°,
0.1
<
e
H
<
0.65
0.1<\frac{e}{H}<0.65
0.1<He<0.65
(2)平板闸门闸孔自由出流流量系数:
m
=
0.60
−
0.18
e
H
m=0.60-0.18\frac{e}{H}
m=0.60−0.18He
应用范围:
0.1
<
e
H
<
0.65
0.1<\frac{e}{H}<0.65
0.1<He<0.65
1.1.3 闸孔出流的淹没系数
闸底坎为宽顶堰(包括平底)的淹没出流,淹没系数
σ
\sigma
σ可查下图。
▲
σ
=
f
(
e
H
,
H
−
t
H
)
\sigma=f(\frac{e}{H},\frac{H-t}{H})
σ=f(He,HH−t)关系曲线
1.2 南科院公式
1.2.1 闸孔出流公式
(1)自由孔流
Q
=
m
b
e
2
g
H
0
Q=mbe\sqrt{2gH_0}
Q=mbe2gH0
式中:Q为过闸流量;m为自由孔流综合流量系数;b为过水断面宽度;e为闸门开启高度;g为重力加速度;
H
0
H_0
H0为计入行近流速水头的堰闸前水深。
(2)淹没孔流3
Q
=
μ
b
e
2
g
(
H
0
−
H
2
)
Q=\mu be\sqrt{2g(H_0-H_2)}
Q=μbe2g(H0−H2)
式中:Q为过闸流量;
μ
\mu
μ为流量系数;b为过水断面宽度;e为闸门开启高度;g为重力加速度;
H
0
H_0
H0为计入行近流速水头的堰闸前水深,
H
2
H_2
H2为下游控制断面的水深。
1.2.2 自由孔流流量系数
同1.1.2。
1.2.3 淹没孔流流量系数
(1)弧形闸门3
μ
=
0.88
−
0.32
(
H
2
H
0
−
0.45
)
\mu=0.88-0.32(\frac{H_2}{H_0}-0.45)
μ=0.88−0.32(H0H2−0.45)
(2)平板闸门
μ
=
0.76
−
0.15
(
H
2
H
0
−
0.45
)
\mu=0.76-0.15(\frac{H_2}{H_0}-0.45)
μ=0.76−0.15(H0H2−0.45)
1.3 Henry公式
1.3.1 孔流公式
Q
=
C
d
b
e
2
g
H
0
Q=C_dbe\sqrt{2gH_0}
Q=Cdbe2gH0
式中:Q为过闸流量;
C
d
C_d
Cd为流量系数;b为过水断面宽度;e为闸门开启高度;g为重力加速度;
H
0
H_0
H0为计入行近流速水头的堰闸前水深。
1.3.2 自由孔流流量系数
判别条件4:
H
0
≥
0.81
h
t
(
h
t
e
)
0.72
H_0\ge 0.81h_t(\frac{h_t}{e})^{0.72}
H0≥0.81ht(eht)0.72
流量系数:
C
d
=
0.611
(
H
0
−
e
H
0
+
15
e
)
0.072
C_d=0.611(\frac{H_0-e}{H_0+15e})^{0.072}
Cd=0.611(H0+15eH0−e)0.072
1.3.3 淹没孔流流量系数
判别条件:
h
t
<
H
0
<
0.81
h
t
(
h
t
e
)
0.72
h_t <H_0 < 0.81h_t(\frac{h_t}{e})^{0.72}
ht<H0<0.81ht(eht)0.72
流量系数:
C
d
=
0.611
(
H
0
−
e
H
0
+
15
e
)
0.072
(
H
0
−
h
t
)
0.7
0.32
[
0.81
h
t
(
h
t
e
)
0.72
−
H
0
]
0.7
+
(
H
0
−
h
t
)
0.7
C_d=\frac{0.611(\frac{H_0-e}{H_0+15e})^{0.072}(H_0-h_t)^{0.7}}{0.32[0.81h_t(\frac{h_t}{e})^{0.72}-H_0]^{0.7}+(H_0-h_t)^{0.7}}
Cd=0.32[0.81ht(eht)0.72−H0]0.7+(H0−ht)0.70.611(H0+15eH0−e)0.072(H0−ht)0.7
1.4 美国陆军工程兵团公式
1.4.1 弧形闸门孔流公式
(1)自由孔流5
Q
=
C
1
C
2
b
e
2
g
h
Q=C_1C_2be\sqrt{2gh}
Q=C1C2be2gh
式中:Q为过闸流量;
C
1
C_1
C1为无底坎时的流量系数;
C
2
C_2
C2为底坎校正系数;b为过水断面宽度;e为闸门开启高度;g为重力加速度;h为上游水深。
(2)淹没孔流
Q
=
C
s
b
h
s
2
g
h
Q=C_sbh_s\sqrt{2gh}
Q=Csbhs2gh
式中:Q为过闸流量;
C
s
C_s
Cs为淹没流量系数;b为过水断面宽度;
h
s
h_s
hs为槛上的尾水深度;g为重力加速度;h为库水位与尾水位的总水头差(包括行近流速水头)。
1.4.2 自由出流流量系数
1.4.3 淹没出流流量系数
二、创新过闸流量公式
2.1 闸门参数
2.1.1 基本公式
Q
=
m
b
e
2
g
H
0
Q=mbe\sqrt{2gH_0}
Q=mbe2gH0
式中:Q为过闸流量;m为流量系数;b为过水断面宽度;e为闸门开启高度;g为重力加速度;
H
0
H_0
H0为计入行近流速水头的堰闸前水深。
2.1.2 自由出流流量系数
m
=
a
l
n
(
e
H
u
−
H
d
)
+
b
m=aln(\frac{e}{H_u-H_d})+b
m=aln(Hu−Hde)+b
式中:m为流量系数;
H
u
H_u
Hu、
H
d
H_d
Hd分别为上下游水头;e为闸门开启高度;a,b为闸门参数(根据实测数据拟合)。
2.1.3 淹没出流流量系数
m
=
a
(
e
H
u
−
H
d
)
2
+
b
(
e
H
u
−
H
d
)
+
c
m=a(\frac{e}{H_u-H_d})^2+b(\frac{e}{H_u-H_d})+c
m=a(Hu−Hde)2+b(Hu−Hde)+c
式中:m为流量系数;
H
u
H_u
Hu、
H
d
H_d
Hd分别为上下游水头;e为闸门开启高度;a,b,c为闸门参数(根据实测数据拟合)。
2.1.4 闸门参数初始值
根据已有过闸流量实测数据,通过最小二乘法,计算出在自由出流和淹没出流情形下的各闸门参数。
下图为曲线拟合过程。
2.2 闸门参数率定
2.2.1 算法简介
差分进化算法(DE)是进化算法的一种,与遗传算法(GA)、粒子群算法(PSO)、蚁群算法、狼群算法等群体智能算法一样,常用来求解复杂优化问题的全局最优解。
2.2.2 目标函数
以模型模拟值与实测值差值的平方和最小为优化目标6,目标函数公式为:
m
i
n
f
=
∑
i
=
1
n
(
Z
(
m
o
d
,
i
)
−
Z
(
o
b
s
,
i
)
)
2
minf= \sum\limits_{i=1}^n(Z_{(mod,i)}-Z_{~(obs,i)})^2
minf=i=1∑n(Z(mod,i)−Z (obs,i))2
式中:f为目标函数;i为时段序号;n为总时段数;Z(mod,i)表示第i时段的模拟值;Z(obs,i)表示第i时段的实测值。
2.3 闸门参数不确定性分析
模型参数经过率定校准后,运行产生的输出结果并不是准确的,只是方差最小的结果。因此存在基于预测结果的概率分布。这就是模型参数的不确定性。可采用零空间蒙特卡洛法、线性分析法(FOSM一阶二次矩法)对闸门参数做不确定性分析。
三、内边界条件线性化
在输水渠道的数值模拟中,考虑到计算的连续性,将节制闸、倒虹吸、涵洞等水工建筑物作为内边界,利用泰勒公式将过闸流量公式线性化,联立连续方程,求解带形矩阵线性方程组得到下一时刻水工建筑物上下断面即内边界的水位及流量。计算堰流流量的经典传统计算公式为:
Q
=
σ
m
b
H
2
g
H
Q=\sigma mb H \sqrt{2gH}
Q=σmbH2gH
将
H
H
H看作自变量,
Q
Q
Q看作因变量,可得到如下函数关系式
Q
(
H
)
=
σ
m
b
2
g
H
1.5
Q(H)=\sigma mb \sqrt{2g} H^{1.5}
Q(H)=σmb2gH1.5
泰勒定理7
设 n n n 是一个正整数。如果定义在一个包含 a a a 的区间上的函数 f f f 在 a a a 点处 n + 1 n+1 n+1 次可导,那么对于这个区间上的任意 x x x,都有:
f ( x ) = f ( a ) + f ′ ( a ) 1 ! ( x − a ) + f ( 2 ) ( a ) 2 ! ( x − a ) 2 + . . . + f ( n ) ( a ) n ! ( x − a ) n + R n ( x ) f(x) = f(a)+ \frac{f'(a)}{1!}(x-a) + \frac{f^{(2)}(a)}{2!}(x-a)^2 +...+ \frac{f^{(n)}(a)}{n!}(x-a)^n + R_n(x) f(x)=f(a)+1!f′(a)(x−a)+2!f(2)(a)(x−a)2+...+n!f(n)(a)(x−a)n+Rn(x)
其中的多项式称为函数在 a a a 处的泰勒展开式,剩余的 R n ( x ) R_n(x) Rn(x)是泰勒公式的余项,是 ( x − a ) n (x-a)^n (x−a)n 的高阶无穷小。
按照泰勒定理,将上述堰流公式函数关系式
Q
(
H
)
Q(H)
Q(H)在
H
i
j
H_i^j
Hij处泰勒展开
Q
(
H
)
=
Q
(
H
i
j
)
+
Q
′
(
H
i
j
)
1
!
(
H
−
H
i
j
)
+
Q
2
(
H
i
j
)
2
!
(
H
−
H
i
j
)
2
+
.
.
.
+
Q
n
(
H
i
j
)
n
!
(
H
−
H
i
j
)
n
+
R
n
(
x
)
Q(H) = Q(H_i^j) + \frac{Q'(H_i^j)}{1!} (H - H_i^j) + \frac{Q^{2}(H_i^j)}{2!} (H - H_i^j)^2 + ... + \frac{Q^{n}(H_i^j)}{n!} (H - H_i^j)^n + R_n(x)
Q(H)=Q(Hij)+1!Q′(Hij)(H−Hij)+2!Q2(Hij)(H−Hij)2+...+n!Qn(Hij)(H−Hij)n+Rn(x)
考虑到从本时刻到下一时刻的堰上水头不会剧烈变化,泰勒展开式可仅保留一阶项,则堰流公式线性化为
Q
(
H
)
≈
Q
(
H
i
j
)
+
Q
′
(
H
i
j
)
1
!
(
H
−
H
i
j
)
=
Q
(
H
i
j
)
+
1.5
σ
m
b
2
g
H
i
j
(
H
−
H
i
j
)
Q(H) \approx Q(H_i^j) + \frac{Q'(H_i^j)}{1!} (H - H_i^j) = Q(H_i^j) + 1.5 \sigma mb \sqrt{2g H_i^j} (H - H_i^j)
Q(H)≈Q(Hij)+1!Q′(Hij)(H−Hij)=Q(Hij)+1.5σmb2gHij(H−Hij)
因此
j
+
1
j + 1
j+1时刻的堰流流量为
Q
(
H
i
j
+
1
)
≈
Q
(
H
i
j
)
+
1.5
σ
m
b
2
g
H
i
j
(
H
i
j
+
1
−
H
i
j
)
Q(H_i^{j+1}) \approx Q(H_i^j) + 1.5 \sigma mb \sqrt{2g H_i^j}(H_i^{j+1} - H_i^j)
Q(Hij+1)≈Q(Hij)+1.5σmb2gHij(Hij+1−Hij)
水头差与水位差相等
H
i
j
+
1
−
H
i
j
=
Z
i
j
+
1
−
Z
i
j
H_i^{j+1} - H_i^j = Z_i^{j+1} - Z_i^j
Hij+1−Hij=Zij+1−Zij
因此
Q
(
H
i
j
+
1
)
≈
Q
(
H
i
j
)
+
1.5
σ
m
b
2
g
H
i
j
(
Z
i
j
+
1
−
Z
i
j
)
Q(H_i^{j+1}) \approx Q(H_i^j) + 1.5 \sigma mb \sqrt{2g H_i^j}(Z_i^{j+1} - Z_i^j)
Q(Hij+1)≈Q(Hij)+1.5σmb2gHij(Zij+1−Zij)
简记为
Q
i
j
+
1
≈
k
Z
i
j
+
1
+
Q
i
j
−
k
Z
i
j
Q_i^{j+1} \approx k Z_i^{j+1} + Q_i^j - k Z_i^j
Qij+1≈kZij+1+Qij−kZij
其中
k
=
1.5
σ
m
b
2
g
H
i
j
k = 1.5 \sigma m b \sqrt {2g H_i^j}
k=1.5σmb2gHij
式中:
σ
\sigma
σ为淹没系数,
m
m
m为流量系数,
b
b
b为闸门总净宽(
m
m
m),
g
g
g为重力加速度(
m
/
s
2
m/s^2
m/s2),
H
H
H为包括行近流速水头的堰上水头(
m
m
m) ,
Z
Z
Z表示水位(
m
m
m),
Q
Q
Q表示流量(
m
3
/
s
m^3/s
m3/s),变量上标表示计算时刻,下标表示断面号。
在内边界所处上下断面间的相容公式为
f
(
x
)
=
{
Q
i
j
+
1
=
Q
i
+
1
j
+
1
Q
i
j
+
1
=
k
Z
i
j
+
1
+
Q
i
j
−
k
Z
i
j
f(x)=\left\{ \begin{aligned} Q_i^{j+1} & = & Q_{i+1}^{j+1} \\ Q_i^{j+1} & = & k Z_i^{j+1} + Q_i^j - k Z_i^j \end{aligned} \right.
f(x)={Qij+1Qij+1==Qi+1j+1kZij+1+Qij−kZij
代入五对角矩阵即可参与水力学计算。
四、参考文献
梁益闻.城市河湖闸泵群防洪排涝优化调度模型研究[D].武汉:华中科技大学,2018. ↩︎
李炜.水力计算手册[K].北京:中国水利水电出版社,2006. ↩︎
毛昶熙.堰闸隧洞的泄流能力计算公式商榷[J].水利学报,1999,(10):38-44. DOI:10.3321/j.issn:0559-9350.1999.10.007. ↩︎ ↩︎
刘孟凯,王长德,闫奕博, 等.弧形闸门过闸流量公式比较分析[J].南水北调与水利科技,2009,7(3):18-19,26. DOI:10.3969/j.issn.1672-1683.2009.03.005. ↩︎
美国陆军工程兵团.水力设计准则[M].北京:水利出版社,1982. ↩︎
Jeremy T. White,Michael N. Fienen,Paul M. Barlow,Dave E. Welter. A tool for efficient, model-independent management optimization under uncertainty[J]. Environmental Modelling and Software,2018,100. ↩︎