《濯河坝水文站水平ADCP测验资料整编方法研究.》由会员分享，可在线阅读，更多相关《濯河坝水文站水平ADCP测验资料整编方法研究.（15页珍藏版）》请在金锄头文库上搜索。
1、濯河坝水文站水平ADCP测验资料整编方法研究邓 彬1 董泽纯1 牛睿平2 王艳红3 于宝荣1（1.重庆市水文水资源勘测局 重庆 401147）（2.南京水利水文自动化研究所 江苏南京 210000）（3.重庆市江北区农林水利局 重庆市 401147）摘要：为有效解决濯河坝水文站受上游水电站调蓄和下游电站回水影响下的流量水文测验难题，2010年濯河坝水文站引进安装了Sontek淘金者SL500水平ADCP，通过近几年的运行，所收集的水位、流速资料系列较为完整，水平ADCP流量与转子流速仪流量的相关性较好。为了深入分析该水平ADCP的特性以及其使用条件，本文对濯河坝水文站2011年收集的水文资料进
2、行了分析和整理，总结水平ADCP使用经验，以期更好地用于重庆水文测验工作之中。关键字：水文站；ADCP；测验；资料整编1 概述流量测验是水文测站工作的主要任务之一。传统的流量测验手段由于测流历时长、自动化程度低、操作复杂等缺点，越来越不能适应水文现代化事业的发展要求。随着社会经济的发展和河流流域的逐步开发，受水利工程影响的水文测站将逐渐增多，使用传统的水文测验方法使得流量测验变得更加困难。 水平式ADCP利用声学多普勒原理通过自动测定河流代表流层的指标流速的方式来进行流量测验，具有不扰动流场、历时短、自动化程度高、实时性好等优点，可有效地解决受水利工程影响的河道流量测验问题，极大地减轻了水文测
3、验劳动强度、提高测验工作效率，是流量测验现代化发展的趋势之一。为有效解决濯河坝水文站受上游水电站调蓄和下游电站回水等影响下的水文测验难题，2010年濯河坝水文站引进和安装了Sontek淘金者SL500水平ADCP。通过近几年的运行，所收集的水位、流速资料系列较为完整，初步分析得出ADCP流量与转子流速仪流量的相关性较好。为了深入分析该水平ADCP的特性以及其使用条件，本文对濯河坝水文站2011年收集的水文资料进行了分析和整理，总结水平ADCP使用经验，以期更好地用于重庆水文测验工作之中。2 水平ADCP流量测验原理ADCP的名称叫声学多普勒流速剖面仪，其利用压电陶瓷作为声学换能器，向水流中发射
4、脉冲声波，然后接收被水中颗粒物散射回来的声波，依据声波脉冲频移，计算出水流的流速和流向。水平ADCP通常安装在河流岸边，从岸边向河流中央水平发射有一定夹角的两条波束，实时测量某一水层的流速分布，作为指标流速。水平ADCP通常要与走航ADCP或转子流速仪等进行比测，率定水平ADCP测量的指标流速与断面平均流速的相关关系。Sontek淘金者SL500水平ADCP还内置水温、纵横摇倾角、超声波水位、压力水位等传感器，用于实时水位采集、温度补偿和参数调整等。根据指标流速与断面平均流速的相关关系，可以将ADCP实测的指标流速反求出断面平均流速，再将断面平均流速乘以过水断面面积求得实时流量。其计算基本公式
5、如下：1、，断面平均流速是指标流速的函数；2、，过水断面面积为水位的函数；3、，流量计算公式；其中： Vc-断面平均流速（m/s），Vi-指标流速或ADCP流速（m/s），H-水位（m），A-过水断面面积（m2），Q-流量（m3/s）。3 水平ADCP测验资料的整编方法3.1 水平ADCP测验资料的预处理 由于自然灾害、仪器故障、水中漂浮物缠绕等原因，水平式ADCP测验数据可能缺测，在资料整理时需要对缺测资料进行插补，通常可以采用上下游水位过程线法、建立数学模型插补等方法。水平ADCP所测得的水文资料时间系列是非平稳的，呈现各种各样的特性，在资料整编前需要对其进行平滑处理。目前小波分析方法和滑
6、动平均法均可以用于水文时间系列的去噪处理。对于受水利工程影响的水文资料系列存在骤变等特性，在平滑处理时采用小波分析方法更为适合一些，利用小波分析方法排除水文数据中的噪音，充分利用小波变换的多分辨率多尺度分析的思想和时频局部化特点，提高了ADCP测验数据的精度和质量。根据近年来小波分析方法在水文时间系列处理结果来看，采用Daubechies小波基函数对水文时间序列进行多层分解，采用硬阀值或软阀值去处水文时间序列中的高频部分，可以有效去除水文时间系列中的随机噪声和突变点1,2。3.2 水位与断面过水面积推算水平ADCP测得的水位数据是由水位传感器测得的水深加上水位传感器高程，由于水平ADCP安装支
7、架的沉降位移、传感器零点漂移等影响，仪器水位资料经过预处理后，按照水位观测规范要求，需用人工水尺水位对水平ADCP水位数据进行订正，得到连续、准确的水位监测数据。由于水文站测验断面存在一定的冲淤变化，在汛期要加强水道断面的测量。在大断面借用时，为保证推算面积的准确性，需要根据推算流量的具体时间选用适合的大断面成果。3.3 指标流速与断面平均流速关系率定指标流速法是美国地质调查局(USGS)于1982年提出并应用，其原理是利用相关分析法或回归法建立局部流速与断面平均流速之间的相关关系，由局部流速推算断面平均流速，进而实现断面流量的实时监测。水平ADCP指标流速与断面平均流速回归模型的选择受仪器安
8、装方式、测验参数设置、断面情况、测站条件、资料质量要求等多方面因素影响。目前，常用的指标流速回归模型5有：一元线性、二元线性、复合线性、一元二次多项式、幂函数、径向神经网络等。在具体应用中可以选用多种回归模型进行比较分析，确定适合的回归模型和比测方式。在水文站测验过程中，水平ADCP可能无法测得测站的整个水位变幅和流量过程，此时可采用其他测流方式作为补充，如枯水和低水可用转子流速仪涉水测验、高水采用水文缆道、浮标、电波流速仪等进行监测。3.4水平ADCP测验资料的精简分析水平ADCP的自动化程度很高，一年采集的数量量很大，以每5分钟自动采集一个数据计算，一年就有10万多个数据。数据量的增大在提
9、高水文资料质量的同时也增加了资料整理工作的难度。可采用精简分析法，从大量的水文测验数据中选出较少的测次来表达水文要素的变化过程，达到减少数据处理量，提高工作效率的目的。水文资料精简方法有手工摘录、计算机自动摘录等。手工摘录的方式存在工作效率低、易错等问题。计算机自动摘录中的折线逼近方法适宜于对遥测水文资料进行的精简处理4。4 应用实例濯河坝水文站位于重庆市黔江区冯家镇阿蓬江流域，东经10846、北纬2921，于2010年1月1日由原濯河坝水位站上迁6km而来，现行集水面积3180平方千米，距河口距离84.5千米。断面形状呈U型，测验河段顺直长约400米，无区间入流，河床由卵石夹砂组成，两岸均为
10、自然护坡。断面控制较好，基下800m为浅滩控制，为低水临时测流断面，各级水位主泓基本稳定，全年最大流速在3.0m/s左右，平均流速在0.102.50m/s，水面宽75140m，最大水深变幅在7.015.0m之间。断面冲淤变化小，汛期水位上涨时漂浮物较多，含沙量较大，落水及平水时漂浮物较少。枯季江水较清。水位测验方式采用人工水尺观测，流量测验方式采用电动水文缆道测验。流量资料采用临时曲线法整编。随着阿蓬江流域的梯级开发，濯河坝水文站测验断面上游约6km建有渔滩电站，下游14km建有箱子岩电站。由于受上游电站调蓄发电和下游电站水库回水的影响，水文测验工作量较大，如图1。为解决水文测验受人为影响的问
11、题，2010年6月濯河坝水文站安装了Sontek公司的淘金者SL500kHz水平ADCP，仪器安装于测验断面左岸，具体位置如图2。图1 2011年6月水位流量过程线图Fig.1 The stage and discharge hydrographdiagram of ZhuoHebahydrological stationin 2011June图2 濯河坝水文站水平ADCP安装位置图Fig.2The horizontal ADCP installationlocation map of ZhuoHebahydrological station首先，将2011年濯河坝水平ADCP数据进行汇总分析
12、，对因设备故障等原因导致的水平ADCP缺测数据进行插补，得到完整的水文时间序列。采用“DB5”小波基函数对2011年水文时间序列3层分解，对高频部分进行软阀值去噪。成果数据显示，小波变换前后2011年水文时间序列的均值不变，方差降低不明显。其次，濯河坝水文站2011年实测流量58份，除去在临时断面测流、无ADCP流量资料对应等无法参与率定的19份，选取39份转子流速仪实测流量资料摘录与之对应的水平ADCP流量资料进行率定分析。分析率定采用两种模型：一是二元线性回归模型，二是径向神经网络模型。二元线性回归模型采用关系式，其中，a、b、c为系数，Vc为断面平均流速（m/s），Vi为指标流速或ADC
13、P流速（m/s），Z为水位（m）。率定过程中为防止系数出现大值，对水位、断面流速数据进行归一化处理，得到基本关系式。径向神经网络模型中采用关系式，其中，Vi为指标流速或ADCP流速（m/s），Z为水位（m）。率定前对水位、流速数据进行归一化处理，率定过程中调整神经网络隐形节点格式使网络零点误差集中。网络误差评定见图3-5。图3 神经网络结构图Fig.3 Neural networkstructure diagram图4 神经网络误差分布图Fig.4 Neural network errordistribution map图5 神经网络相关关系图Fig.5 Neural networkcorre
14、lationdiagram结合濯河坝水文站2011年水文资料整编成果分析，在水位低于406.70m时，水平ADCP流量测验误差大于整编规范3允许范围。由于濯河坝水文站低水水位流量关系呈稳定单一线（见图6），因此，对于水平ADCP低水流量资料需用该站低水水位流量关系线进行改正。图6 濯河坝水文站低水水位流量关系曲线Fig.6 The low stage discharge relationcurve of Zhuodamhydrological station濯河坝水文站2011年水平ADCP流量测验资料有34674个数据，冗余数据量较大，采用折线逼近法4分别对上述模型转换的流量数据，精简后分别
15、得到5477个数据和5488个数据，数据精简到原始数据的15%左右，有效降低了整编工作量。最后，将精简的2011年濯河坝水文站水平ADCP流量时间序列导入“南方资料整编软件2.0”中进行整编计算，得出流量测验成果。将ADCP测验成果与转子流速仪测验成果进行比较（见表1、2、3），可以得出：1、濯河坝水文站水平ADCP指标流速与断面平均流速关系式采用二元线性回归和径向神经网络拟合得出的流量成果基本一致。2、由于濯河坝水文站受电站调蓄影响频繁，人工观测水位很难控制到整个水位变化过程，因此转子流速仪法与水平ADCP的部分月平均流量、月最小流量相差较大，建议安装遥测水位计弥补人工观测水位的不足。表1 流速仪与水平ADCP测法月平均流量对照表Table1 Every monthaverage discharge comparison table 时间 流量测法一月二月三月四月五月六月七月八月九月十月十一月十二月全年转子流速仪法流量(m3/s)20.819.221.857.129.511740.516645.912211749.867.4水平ADCP(线性回归模型)流量(m3/s)19.423.92565.730.112038.516046.4125124
