Page 115 - 2025年第56卷第2期
P. 115

顶向下游的倾斜面与堰顶的夹角为 45°,堰宽与测控闸门孔口的宽度保持一致,流量计算方法采用
              Kindsvater - Carter公式、Bijankhan公式和 Bijankhan - Mahdavi公式    [19] 。通过现场观察确认了闸门流态
              状态后,分别进行了孔流、堰流和临界流态(堰、孔流转换)的试验。当闸门流态为孔流时,在不同开
              度 e下(在 0.02~0.30m范围内调节),根据每种开度 e的设定值来调整水泵变频器,以获取不同工况
              下的闸前水位和过闸流量真 实值。在孔 流条件 下,闸 门 开度 和 闸 前水 位的 调节 幅度 分 别 为 0.02和
              0.01m。而当闸门转为堰流状态时,则借助调节主渠上游水位来确定过闸流量。每组工况下都对闸门
              开度 e、过闸真实流量 Q 、闸前水位 Y、闸上水头 H和闸后水位 h进行了测量。
                                                                           t
                                    m
                  ( 3)数据采集。本试验使用变频器调节水位至待测水位值,待系统完全稳定后,按顺序分别测量
              闸前和闸后的水位值,每个测点重复 3次。若每次测量水位的组间差异小于 0.1mm,则取测量结果的
              平均值作为该点水位值的实测值,否则需要重新测量。闸前水流为亚临界流,将闸前水位测点选在距
              离闸门 1.0m处;闸后流态受水跃影响明显,因此,将闸后水位测点设在距离闸门 2.5m处。水位采
              用精度为±0.1mm的钢尺测量。结合我国末级渠系的灌溉条件                          [17,20] ,试验采用的流量范围为 5~120L?s,
              流量梯度控制在 3~6L?s以内。
                  ( 4)数据分析。本研究流程如图 3所示,分闸门流态识别、流量模型构建以及精度评价三个阶段
              探究测控闸门的水力性能。首先通过试验确定闸门流态判定阈值 γ和大、小孔口分界点 η值;其次采
              用布金汉定理建立闸门流量系数计算公式,并在考虑不同函数类型(包括线性、指数、对数、幂函数、
                                                                                                      2
              多项式等)的基础上,借助最小二乘辨识估计模型参数,综合评估各拟合式的残差、修正后的 R值以
              及方差分析结果;最终选取拟合优度最高且形式简洁的公式作为闸门流量计算模型。




















                                                  图 3 测控闸门流量计算流程

                  为验证流量预测的准确性,本研究采用 75%的数据作为拟合样本,其余 25%的数据用于验证。选
              用平均绝对误差(MAPE)、最大相对误差绝对值(MRE)、误差绝对值小于 5%的数据比例(CP )、拟
                                                                                                    5%
                                 2
              合优度或决定系数( R)和均方根误差(RMSE)5个统计评价指标对流量模型进行评价,各指标的计算
              公式如下:

                                                        100  n  Q -Q  mi
                                                                 ci
                                                MAPE=      ∑                                            (1)
                                                         n i =1   Q mi
                                                            Q- Q
                                                           ( )     × 100
                                                             c
                                                                 m
                                                  MRE = max                                             (2)
                                                              Q
                                                               m
                                                            n 5
                                                      CP = × 100                                        (3)
                                                         5%
                                                             n
                                                           n
                                                         ∑  (Q -Q )    2
                                                               ci
                                                                    mi
                                                  2
                                                 R =1-    i =1                                          (4)
                                                           n
                                                         ∑  (Q -珚 Q )  2
                                                               mi
                                                                     m
                                                          i =1
                                                                                                —  2 5 5 —
   110   111   112   113   114   115   116   117   118   119   120