1957-泾河干流径流量变化趋势分析.docx

? ? 1957 泾河干流径流量变化趋势分析 ? ? 陈沛源，俞 巧，李金文，马金珠 （兰州大学 资源环境学院，甘肃 兰州 730000） 随着社会经济的发展，人民生活水平逐渐提高，水资源利用形势严峻，水资源供需矛盾加剧。 我国不断优化水资源分布、提升水资源利用效率，就是为了不让水资源利用问题成为掣肘社会经济和谐发展的因素［1］。 另一方面，气候环境的变化在一定程度上改变了降雨、蒸发、地表径流等水循环过程［2］。 近些年，世界各地河川径流量出现不同程度的减少，尤其是在干旱和半干旱地区，水资源问题更为突出［3］。 近些年来，有大量的学者对河流径流量进行了一系列研究，且取得了重要的成果［4–11］。 有学者采用小波分析法对沣河和黑河干流径流量年际周期变化和变化原因进行了分析，发现径流量均呈周期性递减趋势，径流量减少受降雨量变化和人类活动的影响，且降雨量变化为主要影响因素［12－13］。 确定径流量突变点、对比流域水文事件发生点，可以有效分析径流量变化原因。 M－K 突变检验、Pettit 检验［14］、有序聚类法、滑动t 检验［15］等方法近些年被广泛应用于径流量序列的分析研究。 王娇等［16］采用R／S 方法对泾河张家山水文站的未来径流量、汛期与非汛期径流量进行预测，发现该水文站的径流量未来呈持续减少趋势。 王随继等［17］利用累积量斜率变化率比较方法分析皇甫川流域径流量变化原因，得出人类活动对径流量减少的贡献率高于降雨量变化的。 泾河流域甘肃段是陇东黄土高原地区的主要水源地，对当地的经济发展尤其重要。 有学者研究发现，受人类活动影响，20 世纪80—90 年代泾河流域水资源量出现逐年递减趋势，进而导致水资源问题更加突出［9］。 随着我国“西部大开发”政策和“一带一路”倡议的实施，当地经济快速发展，对水资源的需求量越来越大，合理分配水资源已成为当地迫切需要解决的问题［18］，对泾河干流径流量变化的研究显得尤为迫切和重要。 因此，笔者以泾河干流甘肃段1957—2016 年的实测径流量资料为基础，利用Mann－Kendall 法、小波变化法和累积距平法对泾河径流量的变化趋势、突变年份和变化周期等进行分析，并初步探讨其产生的原因，以期为当地水资源的合理分配和防洪提供科学的理论依据。 1 研究区概况 泾河是渭河最大支流，是黄河流域十大水系之一，位于黄土高原中部（东经106°14′—108°42′、北纬34°46′—37°19′），发源于宁夏泾源县，流经甘肃省和陕西省，在高陵县注入渭河。 泾河干流全长455.1 km，流域面积45 421 km2，流域内年降水量350 ～600 mm，属于典型的温带大陆性气候区［19－21］。 泾河流域生态环境脆弱、水资源短缺和水资源分配利用的区域差异大等问题，一定程度上制约了区域社会经济的发展［22］。 泾河干流甘肃段流经平凉市崆峒区、崇信县、泾川县，主要支流有颉河、芮河、洪河、蒲河、马莲河、黑河等，其中马莲河为最大的一级支流。 2 资料与方法 2.1 数据来源 数据来源于黄河水利委员会在泾河干流设置的水文站。 各站实测资料序列长短不一，其中崆峒峡站1973—1976 年因修建水库而无实测资料、泾川（三）站1962—1970 年数据缺失。 考虑到径流序列的完整性和一致性，对上述序列进行插补延长后，选取崆峒峡、泾川（三）、杨家坪（二）、平凉4 站1957—2016 年（其中平凉站设站时间为1974 年）的实测径流量和天然径流量序列进行分析。 各站控制面积分别为：崆峒峡站597 km2、平凉站1 305 km2、泾川（三）站3 145 km2、杨家坪（二）站14 156 km2。 2.2 方 法 2.2.1 Mann－Kendall 法 Mann－Kendall 法（M－K 法）为非参数方法，不需要样本遵从一定的分布，不受异常值的干扰［23］，适用于水文气象序列分析［24］。 目前对长时间序列数据的趋势分析常用M－K 法，该方法首先确定原假设H0为趋势不显著，备择假设H1是双边检验，对于所有的k≤n、j≤n且k≠j，Xk和Xj的分布是不同的，检验统计量S计算公式如下： S为正态分布，均值为0，方差Var（S）＝n（n－1）（2n＋5）／18。 当n＞10 时，标准正态变量Z通过下式计算： 若｜Z｜≥Z1－α／2，则原假设不可接受，即在α置信水平上，时间序列数据存在明显的上升或下降趋势。对于统计量Z，大于0 时是上升趋势，小于0 时是下降趋势。 当Z的绝对值大于等于1.28、1.64和2.32 时，代表的分别是通过了置信度90%、95%和99%的显著性检验。 M－K 法可用来检验数据序列的突变点。 首先用逆序排列原始样本序列去计算新序列M－K 检验值，进而得到UBi＝－UFi（i＝n，n－1，…，2，