基于单值化分析方法的下巴沟水文站水位流量关系研究

李晓胜

(甘肃省临洮水文水资源勘测局，兰州 730000)

1 背景资料

1.1 自然地理

洮河是黄河流域较大的支流之一，发源于青海省海南藏族自治州西倾山东麓，于甘肃省永靖县汇入黄河干流。洮河全程长673.1 km，比降0.27%～0.36%，流域面积2.55×104km2。洮河河源与河口落差达2 630 m以上，蕴藏有丰富的水能资源[1]。洮河多年平均径流量为53×108m3，在黄河众多支流中，其年均来水量仅次于渭河，居于第二位。

洮河流域降水量空间分布不均，受海拔、气象等因素的影响，洮河流域年降水量在300～600 mm范围内波动变化。洮河流域年降水量在年内分布上的差异也较为明显，突出表现为5-10月份(夏、秋)气温高、降水多，且以7月份为最大[2]，而以1和12月份的降水量为最小。洮河径流过程季节变化线与降水呈相似规律，但峰值相对降雨稍显滞后，以7-10月份为丰水期，11月份开始锐减，12-2月份为枯水期，3-6月份则水量逐渐增多。流域丰水期径流量占全年径流总量约为56%。

1.2 水文测站概况

下巴沟水文站是国家重要水文测站，建站于1960年1月，位于E103°01′、N34°42′。该水文站是洮河重要控制站，对收集洮河流域水雨情信息和流域防灾减灾发挥着不可替代的作用。

1.3 水文测站断面情况

下巴沟水文站测验河段控制较好，上下游河道呈S型状，断面较稳定。基本水尺断面上游1 km里处建有安果水电站，最大引水流量96.0 m3/s。下巴沟水文站控制断面于每年4和9月份分别测有断面详细高程数据。以每年两次测得的数据均值作为当年断面形状，得到2010-2017年断面套图(图1)，并据此可知下巴沟水文站测验段河道断面相对稳定。平均河底高程与历年变化线(图2)同时显示，下巴沟水文站断面高程受冲淤等影响不大，冲淤交替进行，保证了控制断面良好的一致性。稳定的河床断面有利于该水文站控制断面水位-流量单值化分析研究的开展。

图1 下巴沟水文站2010-2017年断面套图

图2 下巴沟水文站2010-2017年断面平均河底高程与逐年变化线关系

1.4 单值化分析的必要性与可行性

传统水文测验方式必须通过大量测次的水文测验掌握水位流量中的复式绳套变化，以满足现行人工整编或电算整编的要求，以至于测验次数既难以精简次数，又不能停测或者改流量站为水位站。传统水文测验方式已与现代水文工作的高效要求相违背[3]。

水位-流量单值化分析工作，将大幅减轻水文测验工作的频次和强度，显著提高工作效率。这既可优化外业测量的实施，又可以提升测量成果质量，实现以人为本、科学发展的管理模式[4-5]，这将极大推进下巴沟水文站的现代化进程和信息化程度。

下巴沟水文站建立于1960年1月，已采取驻测测验方式不间断收集58年以来的逐日水文资料，并在连续监测的基础上进行数据分析，完成了本流域的基本水文特征数据收集工作。该站控制断面经多年观测和大洪水考验，断面基本稳定；其水位流量关系亦符合《水文巡测规范》[6]等有关单值化分析资料的相关要求。因此，洮河下巴沟水文站水位-流量关系单值化分析是可行的。

2 水位流量关系单值化分析

2.1 资料选取

此次水位流量关系单值化分析采用下巴沟水文站2010-2017年共8年的逐日水位流量资料。样本资料包含了自2010-2017年逐日水位流量资料，数据的测验方法均严格按照《河流流量测验规范》(GB 50179-1993)[7]的相关要求进行收集，资料的可靠性、一致性、代表性审查结果均为良好。选取下巴沟水文站上下游相邻测站作为辅助站，两者距下巴沟水文站分别为13和22 km。

2.2 分析方法选取

现行水位-流量关系单值化分析方法较为常用的有组合因子法[8]、落差指数法[9]等。考虑到目前下巴沟水文站测验断面相对稳定，水位-流量关系主要受洪水涨落和下游回水顶托的影响，利用下巴沟河段现有的水文站网资料，最终采用综合落差指数法[10]开展水位-流量关系单值化研究。

综合落差指数法计算水位-流量关系的公式为：

(1)

式中：q为校正流量，m3/s；Qm为实测流量，m3/s；ΔZm为综合落差，m；α为落差指数。

ΔZm为综合落差，是考虑了上游洪水涨落与下游回水顶托的综合影响，可将其分解为：

ΔZm=km1·ΔZm1+km2·ΔZm2

(2)

式中：ΔZm1为上游辅助水尺至本断面的落差，m；ΔZm2为本断面到下游辅助水尺的落差，m；km1为上游辅助水尺至本断面的落差系数；km2为本断面到下游辅助水尺的落差系数。

整理得出综合落差指数法计算水位-流量关系的公式为：

q=Qm/(km1·ΔZm1+km2·ΔZm2)α

(3)

2.3 落差系数及其确定

落差指数初选值按照取α=0.5，km1、km2的初选值取测站上游辅助水尺至下巴沟水文站基本水尺之间的距离占上下游辅助水尺之间总距离的比值，即km1=13/35=0.371，km2=22/35=0.629。采用试错法，通过MATLAB编程计算，确定最终km1、km2取值。方法为首先固定α=0.5不变，依据方差或标准差最小的原则，然后按合适的计算步长不断调整km1、km2的值，并将其代入式(3)中计算实测流量对应的校正流量，分别绘制水位-校正流量关系曲线。根据适线情况反复调整参数，直到各年水位-校正流量关系曲线有较好的线型分布，且呈单一曲线关系。各年曲线通过符号检验、适线检验和偏离数值检验，最终取值为km1=0.354，km2=0.646，α=0.63。

由此计算所得的下巴沟水文站水位-校正流量关系公式为：

q=Qm/(0.354ΔZm1+0.646ΔZm2)0.63

(4)

3 成果检验

3.1 历年水位流量演变

选取2010-2017年下巴沟水文站水位-流量数据，并点绘历年水位-流量关系图。由图3可知，2012-2014年水位-流量曲线呈现逐年左移的趋势，原因在于水文站上游附近违法占用河道现象严重，使河道收束。经2014年底整治后，恢复河道原貌。但由于气候原因，2010-2013年的水位-流量关系曲线均呈现单一关系曲线，而2014-2017年均出现2条临时低水曲线。

图3 下巴沟水文站历年水位-流量关系图

3.2 水位-流量关系曲线检验

根据《水文巡测规范》(SL 195-1997)的要求，采用 2010-2017年实测流量中分别选取2011、2013和2017年作为丰、平、枯3个代表年份进行单值化分析，利用综合线进行流量推定，并推算年径流总量、汛期总量和一次洪水总量，与原有整编成果进行统计分析，统计结果见表1。结果显示，按照综合落差指数法推算所得的结果，各时段径流量在丰水年结果较原整编方法稍大，而在枯水年则稍小于原整编成果。但通过该方法法得出的下巴沟水文站水位-流量单值化分析成果的精度，符合《水文巡测规范》(SL 195-1997)的相关要求。

表1 下巴沟水文站各时段推算径流量及其误差分析表

4 结 论

利用下巴沟河段现有的站网水位流量资料，采用综合落差指数法开展水位流量关系单值化研究，综合各方面因素得出下巴沟水文站水位-流量单值关系公式。并根据《水文巡测规范》(SL 195-1997)对各年的年总量、汛期总量、一次洪水总量分析显示，该单值关系公式的误差符合相关规定，是可靠的。

将该成果应用于下巴沟水文站的实际工作中，加快本站水文现代化进程[5]，可将大量的人力资源从传统的测验方式中解放出来，同时可以将有限的人员调余济缺最大限度优化人力资源，亦是将有限的技术资源集中起来，做好服务型水文的工作，为社会提供更加优质服务的重要途径[11-12]。