翻斗式雨量计翻斗旋转时间特性分析

2020-09-07廖敏涵刘九夫廖爱民黄一昕李薛刚

水利水运工程学报 2020年4期

廖敏涵，刘九夫，廖爱民，黄一昕，蔡钊，闵星，李薛刚

(1.南京水利科学研究院水文水资源与水利工程科学国家重点实验室，江苏南京 210029；2.南京水利科学研究院滁州综合水文实验基地，安徽滁州 239080；3.河海大学水文水资源学院，江苏南京 210098)

降水是一种最基本的水文现象，降水量是水资源量中最重要的组成部分，因此，精确观测降水量在水资源估算和水文水利计算中有着重要的地位[1-2]。据文献记载，降水量的观测距今已有2 000 多年的历史[3]。人类早期的降水测量是将自然降水通过漏斗引入校准容器中进行降水量的测量[4]，随后各种种类和用途的雨量器、雨量计应运而生。翻斗式雨量计是一种现代化的测量仪器，具有自动化程度高、获取降水量及时性强、降水量资料易于保存和传输的特点，现已被国内外水文、气象和科研部门广泛应用于降雨观测和科学研究中[5-8]。然而，由于翻斗式雨量计仪器本身不够精密以及人为操作不当的原因，不可避免地会产生仪器误差和测量误差。为尽可能减小这种误差，许多科研机构做了大量研究工作，例如世界气象组织（WMO），不仅对翻斗式雨量计进行了实验室比对试验，还进行了场地比测试验[9-13]。这是WMO 首次在实验室条件下组织的雨量计比对试验，主要用于测试不同类型雨量计的工作性能。研究结果表明未经校正的雨量计测量误差较大，而采用正确的校正软件校正过的雨量计则可以准确地测量降水量[6,9]。因此，降水量观测前，需先对雨量计进行率定。

目前，翻斗式雨量计的率定仍是国内外雨量计研究中的热点问题[10-15]，从最初的静态率定，到动态率定，再到现在的智能化率定，逐步趋于成熟和完善。由于翻斗式雨量计翻斗的旋转速度较快，人眼难以捕捉这一过程，常规的测量手段难以达到时间精度要求，所以翻斗的旋转时间测量一直以来都是一个难点且研究严重不足。翻斗旋转时间的研究，有助于进一步了解翻斗在旋转过程中的“水量损失”，以弥补该项研究的不足。Duchon 等[16]曾利用高速摄像装置（帧率为500 fps）对翻斗旋转时间进行试验研究，测得翻斗旋转时间约为0.52 s，但试验的次数较少（仅8 翻斗），且选取的翻斗式雨量计单一（Met One 雨量计），试验结果代表性不足；此外，试验中采用的高速摄像机设备过于昂贵，普通科研工作者难以利用。

因此，本文拟采用数字图像技术（包括视频录制设备和图像处理中的帧间差分法）对国内常见的几种翻斗式雨量计的翻斗旋转时间特性进行系统研究，为设计和生产高精度测量的翻斗式雨量计提供建议与技术支撑。

1 试验设备及方法

1.1 试验设备

本文选择国内最常见的4 种单层翻斗式雨量计进行试验研究，分别是JDZ02，JDZ05，CQS·JD02 和CQS·JD05（见图1），基本参数见表1。4 种雨量计在不同雨强条件下的翻斗旋转时间特性研究内容包括：（1）研究翻斗启动到水平位置的时间（T1）；（2）研究翻斗从启动到另一侧的时间T2（不考虑翻斗的反弹过程）；（3）研究不同雨强下，上述两个时间的变化规律。

图1 4 种型号的翻斗式雨量计Fig.1 Photos of four types of tipping bucket rain gauges

1.2 试验方法

实验室目前具备一种智能化的翻斗式雨量计率定系统（图2，具备供水、自动计数、防风等功能），可模拟不同雨强对翻斗式雨量计进行试验研究。按GB/T 21978.2—2014 《降水量观测仪器第二部分：翻斗式雨量传感器》规定，在0～4 mm/min 雨强范围内，采用5 种雨强（0.4，1.0，2.0，3.0，4.0 mm/min）进行对比试验。在模拟降雨的同时，摄像录制设备（帧率为100 fps）记录单个翻斗式雨量计翻斗启动→翻斗达到水平位置→旋转停止的过程。记录完成后，可用软件进行逐帧读取图像画面，可获得翻斗从启动到翻斗水平的时间和翻斗从启动到停止的时间。每种雨强下分别进行20 翻斗试验。

表1 4 种翻斗式雨量计的基本参数Tab.1 Basic characteristics of four types of tipping bucket rain gauges

图2 智能化的翻斗式雨量计率定系统Fig.2 Intelligent tipping bucket rain gauge calibration system

翻斗旋转是一个变加速的过程，启动时非常缓慢，然后慢慢加速到水平位置，最后再极快地旋转到另外一侧，完成1 次翻斗过程。为了解决这一问题，本文拟采用图像处理中的帧间差分法来判断物体是否运动。帧间差分法常用在交通领域中，判断车辆是否移动，或用在智能监控领域中[17-19]。该方法在本文不做重复叙述，图3 为该方法的结果。通过连续增大的SUM 值来判断图像中存在运动，即翻斗旋转；无连续增大的区域，说明图像是静止的，翻斗也未旋转。基于帧间差分法，并使用Matlab 软件编程来进行图像的逐帧处理，自动保留视频序列中运动的部分，静止的部分则剔除。保留下的运动部分则用来分析翻斗旋转的时间特性。图4 为CQS·JD05 在4 mm/min 雨强下一次翻斗旋转的过程，选取了8 个画面进行展示。图4（a）为算法识别的起始状态，记为T=0 s；图4（b）至图4（e）都是翻斗从启动到翻斗到达水平位置的状态；图4（f）为翻斗到达水平位置的时刻；图4（g）为翻斗从水平位置开始向右侧旋转；图4（h）为翻斗旋转到达右侧并触底，完成1 次旋转过程。

图3 帧间差分法的结果Fig.3 The result of the inter-frame difference algorithm

图4 翻斗旋转过程的一个例子Fig.4 An example of the bucket rotation

2 试验结果与讨论

2.1 视频帧数对试验结果的影响

在试验条件相同的情况下，以JDZ05 翻斗式雨量计为例，模拟雨强使翻斗进行20 次旋转，用50 fps与100 fps 的试验结果进行对比。表2 为50 fps 视频试验结果与100 fps 结果的对比。可见：50 fps 的最小值都偏大，但最大值、极差、标准差和平均值都偏大，说明50 fps 的视频测量精度不够，100 fps 的视频测量精度要更高；另一方面，由于翻斗旋转的速度较快，50 fps 视频有的时候不能获取到翻斗到达水平位置的图像，而100 fps 的视频基本不会出现这种情况。然而，需要用更高帧率的视频进行拍摄吗？如达到国外学者使用500 fps 的高速摄像装置进行视频拍摄[16]。笔者认为答案是否定的，不一定要使用更高的帧率进行拍摄，原因如下：(1)在进行50 fps 与100 fps 视频对比分析时发现，100 fps 的数据文件要比50 fps 的大1 倍，即处理时间上要多出1 倍的时间；(2)在阅读国外学者使用500 fps 高速摄像装置所做研究后，发现该学者只分析了8 次翻斗旋转的时间。Duchon 使用的方法非常新颖，但是工作量却非常少，原因可能在于数据量太大，难以做更多的斗数来分析。综上所述，帧率越高，效率可能并不会提高，只要精度满足了要求即可。故本试验采用100 fps 的视频进行拍摄，获取试验数据进行分析。其中，由于该项分析仅对帧率进行对比，且0.4 mm/min 雨强进行试验较为耗时，故未设置该雨强进行试验。

表2 视频帧数对比试验结果Tab.2 Comparison of video frame rates

2.2 不同类型翻斗式雨量计翻斗旋转时间特性分析

图5 为4 种类型的翻斗式雨量计旋转时间平均值（T1与T2）与雨强之间的关系。JDZ02 和CQS·JD02的平均值随雨强的增大而减小，R2分别为0.99 和0.86，有较好的线性关系；JDZ05 和CQS·JD05 线性关系不明显，R2均小于0.6，相关性较差。Duchon 使用Met One 翻斗式雨量计进行试验研究，认为翻斗旋转的时间与雨强无关[16]。WMO 曾在2006 年发表的《WMO laboratory intercomparison of rainfall intensity gauges》报告中指出：“翻斗旋转时间ΔT 是雨强的函数”，即雨强与ΔT 存在关系，但ΔT 很难测量[9]。Marsalek 认为翻斗式雨量计的旋转时间在0.3～0.6 s[11]，这与本次试验结果接近。

图5 翻斗旋转时间与雨强的关系Fig.5 Relationship of rotation time of tipping bucket with rainfall intensity

表3 是4 种类型的翻斗式雨量计在不同雨强下的试验结果，试验斗数均为20 斗。由试验数据可得：

（1）从极差上看，JDZ02 和JDZ05 在4.0 mm/min 雨强下极差最小，在0.4 mm/min 雨强下极差最大，极差随雨强的减小而增大。CQS·JD02 和CQS·JD05 的极差与雨强变化无直接关系，极差随雨强变化不大。

表3 4 种类型的翻斗式雨量计在不同雨强下的试验结果Tab.3 Test results of four types of tipping bucket rain gauges under different rainfall intensities

（2）从平均值上看，JDZ02 和CQS·JD02 的平均值在4.0 mm/min 雨强下最小，翻斗旋转最快，在0.4 mm/min雨强下最大，旋转最慢，平均值随雨强的增大而减小。但是对于JDZ05 和CQS·JD05 两款分辨率为0.5 mm的雨量计来说线性关系不明显，这两款雨量计翻斗旋转时间的平均值在4.0 mm/min 雨强下最小，但在0.4 mm/min 雨强下并不是最大。

（3）从标准差上看，4 款雨量计在4.0 mm/min 雨强下标准差最小，在0.4 和1.0 mm/min 雨强下标准差最大。方差的意义是说明翻斗旋转时间的波动情况，结果表明在4.0 mm/min 雨强下，翻斗旋转时间波动最小，在0.4 和1.0 mm/min 雨强下波动最大。

综上所述，雨量计的翻斗在4.0 mm/min 雨强下的旋转时间较为稳定，波动很小，但在0.4 mm/min 雨强下的旋转时间却不稳定，波动较大。翻斗式雨量计由于其自身结构的特点，在大雨强时会造成翻斗旋转过程中的“水量损失”，使得记录的水量偏小，造成测量误差。研究翻斗旋转时间特征可为进一步认识翻斗式雨量计的计量误差特性提供理论依据。水量损失的计算见式（1），误差见式（2）。

式中：Δm 为翻斗旋转过程中的水量损失；T1为翻斗从启动到水平位置时间；I 为雨强；ρ 为去离子水的密度（1 g/cm3），本试验采用去离子水，减少密度对试验结果的影响。

式中：Er为翻斗式雨量计的计量误差；m1为理论上1 翻斗的质量。

针对翻斗式雨量计翻斗旋转时间在小雨强下不稳定，波动较大的问题，建议使用虹吸部件。虹吸部件安装在承雨器与翻斗漏斗之间，起调节雨强作用，使得不均匀的降雨能够以相对均匀的方式流入到翻斗中（经过虹吸部件调节后，流入翻斗的雨强会变大），此时调整翻斗对应时刻流入雨强的计量误差使之趋近于0，以减少雨强对翻斗计量的误差[14,20]。

2.3 翻斗旋转角度对旋转时间影响分析

翻斗式雨量计翻斗的旋转受诸多因素影响，如雨强、翻斗质量、翻斗翻转角度等。本研究主要分析翻斗旋转角度对旋转时间的影响。4 种翻斗式雨量计的翻斗旋转角度已在表1 中给出，由图5 可知，JDZ02和JDZ05 两种相同厂家生产的雨量计，翻斗旋转角度越大，T1与T2也越大；对于CQS·JD02 和CQS·JD05两种相同厂家生产的雨量计情况也如此。但对于4 种翻斗式雨量计来说，并不是翻斗旋转角度越大，旋转时间也就越长，主要的区别在于翻斗的质量不同。