新疆阿克苏河流域洪水演变趋势研究

蒋军新，蔡 明，徐永军，方功焕，李 稚，陈永金

（1.新疆塔里木河流域阿克苏管理局，新疆阿克苏843000；2.中国科学院新疆生态与地理研究所荒漠与绿洲生态国家重点实验室，新疆乌鲁木齐830011；3.聊城大学环境与规划学院，山东聊城252059）

0 引言

洪水灾害是一种突发性强、发生频率高、危害严重的气候灾害［1］。在过去几十年，世界各地因洪水造成的损失多达数百亿美元，造成几千人死亡［1］。因此，洪水变化研究成为水文学中的热点，如洪水频率与演变分析、洪水工程设计、洪水预估与预警等［2-5］。

在气候变化背景下，天山山区水循环与水文过程发生了改变，洪水的发生特征与演变趋势也相应发生了变化［4，6］。阿克苏河发源于天山南坡，由高山区冰川融水、积雪融水和降水共同补给，对气候变化响应敏感［7］，是北半球中纬度高山区最具代表性的河流。理解阿克苏河流域洪水变化对于理解气候变化背景下中纬度高寒山区的洪水演变过程具有重要意义。

阿克苏河担负着850万亩农田的灌溉，养育了阿克苏地区维吾尔族、汉族、回族、蒙古族、哈萨克族、柯尔克孜族等150万人口［8］。近300年来旱洪灾害呈增加趋势，且洪水事件较干旱事件明显［9］。洪水的发生频率和强度增强，有记录的3次最强洪水均发生在1990年以来［10-11］。毛炜峄等［12］指出阿克苏河的洪水与0℃层高度密切相关；沈永平等［13］分析了近百年来阿克苏河流域麦茨巴赫冰川湖溃决突发洪水的演变，发现冰湖溃决突发洪水的总量由20世纪70年代的1×108m3增加到1990年代以来的3×108～4×108m3。随着气候变化，阿克苏河流域的洪水变的更加复杂，暴雨洪水和融冰融雪洪水发生了改变。同时，随着冰川加速融化，末端退缩或平均冰川厚度减薄，导致冰川作用区的冰湖面积增大，冰湖溃决灾害事件增多［14］，严重威胁下游公路、桥梁、水电站、水库等水利工程及工农业发展［15-16］。在气候变化背景下，研究阿克苏河流域的洪水变化趋势及形成机制对下游绿洲环境及财产安全保护有重要意义。

本文基于长时间序列的阿克苏河流域逐日出山口流量数据，利用块最大值抽样方法和超定量阈值（POT）抽样方法提取了有器测资料以来阿克苏河流域两条支流的洪水事件，分析了洪水强度和发生时间的变化趋势，并解析了年和季节最大洪水与气候变化的关系。研究成果不仅可服务于地方防洪与社会经济发展规划，而且对理解天山洪水演变趋势及灾害评估可以提供科学依据。

1 研究区概况与研究方法

1.1 阿克苏河流域

阿克苏河属于典型的北半球中纬度高山河流，流域总面积5×104km2。阿克苏河有两条支流，其中，西支为托什干河，北支为库玛拉克河，两条支流的多年平均径流量分别为26.63×108m3和47.88×108m3（图1）。阿克苏河是唯一一条常年向塔里木河输水的河流，多年平均下泄塔里木河水量为35.61×108m3，占塔里木河干流补给量的70%～80%，是塔里木河干流最大的补给来源［17］。

阿克苏河主要由降水、冰川与积雪融水和基岩裂隙水补给［18］。由于补给类型的多样性，阿克苏河流域的洪水类型也复杂多样，不仅包括暴雨洪水，还有冰川积雪融水产生的洪水和冰川湖溃决突发洪水。其中，冰湖溃决突发洪水主要发生在库玛拉克河上的麦茨巴赫冰川湖。麦兹巴赫湖位于中国—吉尔吉斯斯坦边界的天山托木尔—汗腾格里山区，在北伊利尔切克冰川表面，受阻于南伊利尔切克冰川并接受来自两支冰川的融水蓄积而成［15，19-20］（图1）。麦兹巴赫湖是天山山区典型的冰川堵塞湖，长3.4 km，宽1.2 km，面积4 km2，平均水深44 m［19］。麦兹巴赫湖的排水机制主要为冰内和冰下水道扩大排水［13，21-23］。只有温度上升到0℃以上时才能产生消融。冰湖溃决突发洪水受水位和入湖水温度的双重影响，当达到一定条件时便发生溃决。

图1 阿克苏河流域水系分布与麦兹巴赫冰川湖位置Fig.1 The hydrological system of the Aksu River and the location of the Metzbach Lake

本文用到的水文数据是库玛拉克河协合拉水文站和托什干河沙里桂兰克水文站的1958—2011年的日流量数据（图1），用来提取阿克苏河两条支流的洪水事件。由于阿克苏河的西大桥水文站受人类活动影响较大，人为调控因素较强，本文仅分析两条支流的洪水变化特征，不考虑干流的洪水变化。本文用到的气象数据是有长时间序列的阿克苏气象站和阿合奇气象站的1958—2011年的日气温和日降水数据（图1），分别代表库玛拉克河和托什干河的气候情况。

1.2 阿克苏河流域洪水提取方法

块最大值采样方法就是取一个时间步长（如年、季节、月份等）内的最大日流量作为洪水事件（公式1），但是该方法容易忽略该步长内发生的低于最大洪水事件的其他较大洪水信息，而且对于没有发生洪水的时间步长，块最大值法同样采集了样本数据［24］。为了克服块最大值抽样所具有的缺陷，本文同时利用POT采样方法对块最大值采样进行补充，例如洪水发生次数、相应量级以及峰现时间等信息。本方法在干旱区内陆河流域得到了较好的应用［11，25］。

式中：Qi为第i个时间步长（如第i年，第i个春季等）内提取的洪水事件，q1，q2，…，qj，…，qn为第i时间步长内的日流量序列。

POT方法通过选取超过某一阈值的流量，组成洪水极值序列。阈值的选取非常重要，但是目前还缺乏公认的客观方法［26］。选择洪峰的主要思路是：在连续的峰值流量过程中，只能选择一个最大的峰值。由于阿克苏河属于典型的季节性河流，夏季流量高，冬季低，本文分别考虑不同季节的流量特征，通过测试不同的阈值，以逐渐增加到作为洪水发生阈值的备选值［27］。另外，两个洪峰还必须满足洪峰独立性标准。结合年平均洪水发生次数，选择满足平均发生次数在1.60～3.00范围内的较大流量作为阈值［28］，最终选取作为不同季节的阈值，进而选取出阿克苏河流域洪水极值事件的序列［24，29-30］。在形成序列前，采用美国水资源协会提出的判别标准对洪峰独立性进行判别［式（2）］。

式中：D为连续两个洪峰的间隔时间；A为两条支流的汇流面积，单位为平方英里，合2.59 km2；Q1和Q2分别为连续两个洪峰的量级，单位为m3·s-1。考虑到不同流域的汇流时间，两条支流的两次洪水事件的时间应该不小于13 d。不满足上述条件的连续洪峰中，只取其中最大一次洪峰。

1.3 洪水变化趋势检验方法

本研究利用非参数Mann-Kendall秩次统计方法来分析阿克苏河两条支流的年和季节最大洪水的变化趋势。非参数Mann-Kendall趋势检验方法广泛的应用于水质、径流、温度、降水等水文气象时间序列的趋势变化检验中，其特点是不必先假定样本的统计分布，在趋势分析中具有广泛的应用［31-32］。Mann-Kendall秩 次 统 计 方 法 可 用 下 式表示：

其中：

式中：xk、xi为连续的数据值；n为数据资料的长度，本文为年数。在Mann-Kendall检验中，当|Zc|＞Z1-α/2时，拒绝零假设。其中，Z1-α/2为标准正态方差；α为显著性检验水平；Zc为正值表示增加趋势。当Zc的绝对值大于1.28、1.64、2.32时分别表示通过了信度90%、95%、99%的显著性检验。

2 结果与分析

2.1 基于POT方法的洪水事件选取

根据阿克苏河两条支流的年内分布特征（图2），11月—翌年4月为库玛拉克河的干季，5—10月为湿季。对于托什干河，由于融雪水占比较高，4月份流量较高，春季洪峰明显，因此，本文将11月—翌年3月定义为托什干河的干季，4—10月为湿季。利用POT采样方法，确定协合拉干季和湿季的阈值均为即77 m3·s-1和730 m3·s-1。在1958—2011年，基于此阈值选取的洪水事件在库玛拉克河共出现106次，其中干季17次，湿季89次。托什干河沙里桂兰克站干季和湿季的阈值分别为40 m3·s-1和360 m3·s-1，在1958—2011年共出现112次洪水事件，其中干季14次，湿季98次，平均每年发生洪水次数为2.07次。

图2 库玛拉克河协合拉水文站和托什干河沙里桂兰克水文站年内流量分布（a）、（c）和年洪水发生次数-POT采样阈值关系图（b）、（d）Fig.2 The intra-annual distribution of streamflow（a），（c）and the relationship between the POT thresholds and flood times（b），（d）of the Xiehela station of the Kumarak River and the Shaliguilank station of the Toshgan River

根据利用POT方法提取的库玛拉克河的106次洪水过程和托什干河的112次洪水流量过程，两条河流的洪水平均持续时间（从开始上涨到完全消退）分别为10.4 d和8.8 d，上涨时间分别为5.8 d和4.1 d，下降时间分别为4.6 d和4.7 d。其中，库玛拉克河的最高洪水持续时间可达26 d，从1988年11月26日起涨，在12月12日达到洪峰，12月22日结束；托什干河的最高洪水持续时间为21 d，从1967年4月11日起涨，在4月26日达 到 峰 值，5月2日结束。

不同补给类型河流的洪水的年内分布不同。以融冰融雪水补给为主的库玛拉克河，洪水主要发生在7—8月，其中7月和8月的洪水发生次数分别为37和42次，占洪水总发生次数的74.5%。对于托什干河，洪水主要发生在6—8月，分别发生了27、32和20次洪水事件，占洪水总发生次数的70.5%（图3）。

图3 库玛拉克河协合拉水文站和托什干河沙里桂兰克水文站洪水发生年内分布Fig.3 Intra-annual distribution of flood occurrence at the Xiehela station of the Kumarak River and the Shaliguilank station of the Toshgan River

2.2 基于POT方法的阿克苏河洪水变化特征分析

阿克苏河流域两条支流库玛拉克河和托什干河的洪水对降水的响应过程可见图4。由于库玛拉克河和托什干河的洪水主要发生在夏季，图4仅给出第180～240 d的降水和洪水关系。协合拉水文站洪水发生前15 d降水相对较小，洪水受暴雨影响较小，而对于托什干河，洪水发生前的降水频率和强度均较高。由此表明，托什干河更容易形成暴雨洪水，这与库玛拉克河和托什干河的径流补给特点相关。库玛拉克河主要是冰川融水补给，而托什干河主要是降雨和融雪水补给，库玛拉克河的洪水对气温变化更为敏感［33-35］。

图4 库玛拉克河和托什干河的洪水与降水关系Fig.4 Relationship between the timing of floods and precipitation of the Kumarak River and the Toshgan River

2.3 基于块最大值法方法的阿克苏河洪水变化趋势

对于库玛拉克河，年最大洪水表现出总体增加趋势，增加速率为8.48 m3·s-1，Mann-Kendall趋势检验的Zc值可达3.425，显著性水平为P＜0.01。对于春季最大洪水，同样表现为增加趋势，Zc值为2.104，显著性水平为P＜0.05。夏季最大洪水的变化特征与年最大洪水变化特征类似。对于秋季最大洪水，Zc值为-1.388，表现为不显著下降趋势（表1和图5）。

对于托什干河，年最大洪水表现出总体不显著增加趋势，增加速率为3.40 m3·s-1，Mann-Kendall趋势检验的Zc值为1.567。虽然四个季节的洪水均表现为增加趋势，但是只有春季、秋季和冬季的洪水表现为显著增加趋势，Zc值分别为2.328、2.209和4.619。夏季最大洪水虽然也表现为增大趋势，但是没有通过显著性检验。

分析库玛拉克河和托什干河的洪水与气候因子的关系发现，库玛拉克河的春季洪水与最高连续7 d气温具有相同的变化特征，二者的相关性较强，Spearman相关系数为0.44（图5和表1），推测库玛拉克河春季洪水主要是由融雪水补给为主。夏季、秋季和冬季洪水大小也与7 d最高气温呈正相关关系，并且夏季的洪水大小与气温的相关性通过了显著性检验（P＜0.05）。这主要是由于库玛拉克河夏季径流由冰川融水补给为主，冰川融水补给比例高达60%［36］。值得注意的是，冬季最大洪水在1988年和1996年分别达到了528 m3·s-1和720 m3·s-1，远大于冬季最大洪水的中位数值38.7 m3·s-1，这是由于在1988年和1996年的12月份发生了冰湖溃决突发洪水，导致河道流量骤增。由于库玛拉克河的洪水大部分与麦兹巴赫冰川湖的溃决相关，而本冰川湖的溃决与气温和降水没有直接关系［22，37］，而与冰坝的稳定性、湖面温度、补给速度、冰川内部储存水的快速释放有关［23］，具有非常强的随机性，因此，洪水的变化规律与最高连续7 d气温和最高连续3 d降水的变化规律不甚一致。托什干河的冰川面积相对较小，仅占流域面积的3.43%，冰川融水径流不足出山口径流量的20%［33］。总体来看，托什干河的洪水增加与气温关系微弱，而受最高连续3 d降水量影响较大。秋季洪水的增加主要是由于最高连续3 d降水增加所致。

图5 库玛拉克河和托什干河年际和各季节洪水、最高连续7 d气温（T7）和最高连续3 d降水（P3）的变化趋势Fig.5 The annual and seasonal variations of flood magnitude，seven-day maximum temperature，and three-day maximum precipitation in the Kumarak River and the Toshgan River

表1 库玛拉克河和托什干河最大洪水变化的Mann-Kendall检验结果，以及洪水与最高连续7 d气温（T7）和最高连续3 d降水（P3）之间的相关系数Table 1 Trend analysis of the annual and seasonal maximum floods of the Kumarak River and the Toshgan River based on the Mann-Kendall test.The correlation coefficients between flood magnitude and seven-day maximum temperature，three-day maximum precipitation were shown.The superscript star（*）and（**）represent the p-value of the trend analysis or the correlation are less than 0.01 or 0.05，respectively

2.4 基于POT方法的阿克苏河流域洪水发生时间变化特征

库玛拉克河和托什干河的洪水发生次数有增加的趋势，并且有略提前的趋势（图6）。对于库玛拉克河，1958—2011年共发生7次10年一遇的洪水，发生日期分别为1983年8月21日、1987年8月19日、1994年7月23日、1996年8月14日、1997年7月31日、1999年7月31日和2005年7月15日，对应的最高洪峰流量分别为1 700 m3·s-1、1 600 m3·s-1、2 180 m3·s-1、1 590 m3·s-1、1 710 m3·s-1、1 850 m3·s-1和1 620 m3·s-1。对于托什干河，1958—2011年共发生5次10年一遇的洪水，分别发生在1969年8月3日、1995年5月10日、1995年6月8日、1999年7月31日和2001年8月17日，对应的最高流量分别为989 m3·s-1、1 020 m3·s-1、922 m3·s-1、1 270 m3·s-1和1 110 m3·s-1，对应的重现期分别为20年、23年、15年、58年和33年。可以看出，库玛拉克河和托什干河高强度洪水事件在20世纪90年代以来频繁发生，分别占有器测资料以来总发生次数的71%和80%。

图6 阿克苏流域洪水发生时间变化特征Fig.6 Variations of flood magnitude and timing of the Kumarak River（a）and the Toshgan River（b）from 1958 to 2011.The flood timings of the Kumarak River and the Toshgan River of specific period were shown in subplot（c）and（d）

在年内分布上，以冰川融水补给为主的库玛拉克河的洪水事件趋向于提前发生，而以积雪融水洪水和暴雨洪水为主的托什干河的洪水发生时间变得更加分散，表现为春季洪水发生提前而秋季洪水发生推迟。这主要是由于在全球变暖背景下，阿克苏河流域的温度不断上升，库玛拉克河的冰川融化速度加快，融水增多，导致夏季洪水提前，湿季洪水事件的平均发生时间约提前了6 d［图6（c）］。而对于托什干河，主要是4—5月份的洪水有提前趋势，主要是由于春季融雪水增加所致，平均洪水发生时间提前了18 d［图6（d）］。

3 结论

阿克苏河流域的产汇流机制复杂，多数洪水属于混合型洪水。通过分析1958年以来的洪水变化趋势，可以得出以下结论：

（1）阿克苏河流域的库玛拉克河和托什干河季节性强，当用POT方法提取洪水事件时，分枯水期和丰水期采用不同的洪水阈值可以更加实际的提取出洪水事件。在1958—2011年，库玛拉克河和托什干河分别提取出洪水事件106次和112次，其中干季分别为17和14次，湿季分别为89和98次。

（2）1958—2011年库玛拉克河和托什干河的洪水均表现出增加趋势。除了库玛拉克河的秋季最大洪水表现出降低趋势外，两河的季节洪水均表现为不同程度的增加态势。10年一遇洪水事件绝大多数发生在20世纪90年代以后。

（3）在洪水发生时间上，库玛拉克河的湿季洪水发生时间有提前趋势，而对于以暴雨产流和融雪产流为主的托什干河，春季洪水的发生时间有提前趋势而秋季洪水发生时间有推后现象，洪水发生时间更加分散，为水资源管理者带来挑战。对于以冰川融水补给为主的库玛拉克河，春季、夏季洪水发生与最高连续7 d气温相关，而托什干河年洪水以及春季和秋季洪水的发生与最高连续3 d降水量相关。