Page 125 - 2023年第54卷第10期
P. 125
由式(22)(24)(26)得到抛物线型断面的马斯京根法分河段参数K、X与水动力模型河宽 B、河
i i p
道比降i、最大水深 H 之间的转换关系。
0p p
根据K、X与等效河道要素的转换关系,通过马斯京根参数K、X、总河段长 L、洪峰流量 Q、糙
i
i
i
i
率 n、平均流速转换波速系数 α及分河段数 N,可以计算得到不同概化断面下等效河道的河宽 B、河
道比降i、最大水深 H等河道信息。本文只提出了抛物线型和矩形这两种概化断面下的转换关系,根
0
据相同的原理,其他形状的概化断面(如三角形、梯形等)也可推导求出。
4 基于 K、X等效河道的具体应用及合理性分析
4.1 基于 K、X等效河道的具体应用 何惠等 [22] 提出了用最小二乘法对马斯京根法的 K、X参数
进行最优估计, 并 以 1961年 海 河 流 域 南 运 河 称 沟 湾 至 临 清 段 的 一 次 洪 水 过 程 为 例, 其 中 K=
13.05、X =- 0 .2716。现以该方法率定 的 K、 X值 运 用 于 本 文 提 出 的 等 效 河 道 方 法, 对 该 场 次 洪
水进行演算。
南运河位于河北省南部,南运河称沟湾至临清段全长 83.8km,如图 4(a)所示。本文提出的等
效河道方法对该河段的概化如图 4(b)(c)所示。整个河道概化为一条长 83.8km的长河段,取 N为
3
100,将河段等距离划分为 100个子河段,取糙率 n为 0.02,洪峰流量 Q为 550m ?s。经计算,当
河道断面概化为矩形时,河道全程比降为 0.0045%,河宽为 36.04m;当河道断面概化为抛物线型
时,河道全程比降为 0.0042%,河宽为 51.62m。K、X值对 应 的 等 效 河 道 信 息, 转 换 后 可 用 于 水
动力模型的边界条件。由临清断面形状信息,结合曼宁公式及流速流量关系可推出临清断面水 位
流量关系。取南运河称沟湾实测流量为上边界条件,临清断面的水位流量关系为下边界条件进 行
模拟。
图 4 实际河道及下边界水位流量关系
图 5和表 1为出口临清断面实测流量、参考文献 [22] 方法演算流量及本文方法演算流量对比。
2
从表 1统计结果来看,矩形断面下,均方根误差 RMSE为 10.71,比文献方法大 4.64;可决系数 R为
0.99,与文献方法相近,均在 0.99以上。抛物线型断面下,均方根误差 RMSE为 10.76,比文献方法大
2
4.69;可决系数 R为 0.99,与文献方法相近,均在 0.99以上。本算例下,两种概化方式演算结果相差
不大,矩形的概化方式表现更好。从图 5流量过程线来看,本文提出的两种概化方法演算的临清流量
2
— 1 6 1 —