第4章 塔里木河流域洪水过程分析

4.1　洪水类型与特征

塔里木河流域地域广阔,气候和地形变化复杂,从而形成类型多样、不同特征的突发性洪水,按洪水成因可分为暴雨洪水、融雪洪水、混合型洪水和冰川湖溃决洪水。

4.1.1　暴雨洪水

暴雨洪水是干旱区的主要洪水类型,由区域性暴雨和局地性暴雨形成,与水汽来源、影响暴雨的天气系统、暴雨中心的活动路径、暴雨强度以及下垫面流域特性等因素有关。按形成暴雨洪水的原因和暴雨洪水的特性来分,可分为大尺度天气系统形成的区域性暴雨洪水和中、小尺度天气系统形成的局地性暴雨洪水。

塔里木河的暴雨洪水是针对上游源流而言的,据资料,塔里木河源流区暴雨以阿克苏与喀什出现次数相对较多,由于其发生范围小,即使强度高,总降水量小,产流量不大,除对位于暴雨中心的小河流及山前地段会造成局部危害外,对塔里木河干流河道影响不大。由于源流区有许多机械风化的碎屑物质和冰碛物,在发生暴雨洪水时也经常伴随着泥石流的发生。从垂直地带分布规律可以看出,高山带由于地势较高,液态降水范围有限,从而制约了暴雨洪水的规模,暴雨洪水多发生在中低山带,与气象条件和下垫面有着密切关系,如植被覆被情况及受水区的地形条件及其组合,与地面组成物质的渗水性能也有很大关系。

塔里木河源流区暴雨洪水的特点是来势凶猛、陡涨陡落、年际变化大、挟沙量大、破坏力强,通常暴雨后24h内洪峰即至,历时短暂,由几小时到几天不等,在流量过程线上呈明显的尖峰状,流速很大,但洪量不大。如天山中段南坡的托什干河在1969年7月15日前后发生一场局地暴雨洪水(图4-1),洪峰流量达1070m3/s,洪峰流量是起涨流量的4.4倍,洪水历时56h,场次洪水总量仅1.0×108 m3,不足年总量的2.7%;昆马力克河流域1969年7月发生的暴雨洪水过程见图4-2。

图4-1　托什干河1969年7月暴雨洪水过程

图4-2　昆马力克河1969年7月暴雨洪水过程

塔里木河源流区暴雨的变率远远超过年降水量,与温度条件密切相关的春季融雪洪水和夏季消融洪水相比,暴雨洪水年际变化极大。在低山带的一些山洪沟,平时干涸,只有暴雨后才有水流,这种情况可能几年甚至数十年才能出现一次。暴雨洪水发生时挟沙量大,有些暴雨洪水伴随有泥石流,同时也溶解了大量的盐分,矿化度可达1～3g/L。

4.1.2　融雪洪水

融雪型洪水发生在夏季持续高温天气的情况下。在新疆,主要是由西风带槽脊的变化和西藏高压北挺共同形成的高压脊,其中心位于西藏高原和南疆上空,呈南北或东北—西南向的暖高压脊。这种高压脊造成新疆持续高温,促成高山冰雪强烈形成冰川融水型洪水。冰雪的积累和消融,取决于太阳辐射、气温、冰川面积和积雪量、冰川的类型和分布、流域内地形地貌等气候特征。这类洪水的特点是由于气温和辐射平衡都有明显的日变化过程,洪水总流量大但峰值不高,洪峰流量过程线肥胖,持续时间长。按其形成的原因可分为季节积雪融化洪水和高山冰川消融洪水两种情况。

4.1.2.1 季节积雪融化洪水

季节积雪融化洪水是指上半年冬季和当年春季积雪,在春季和夏初气温回升时,零度等温线上升,位于中山和部分高山区的积雪开始消融产生融雪径流,汇入河道形成融雪洪水。由于洪水是山区的季节积雪融化形成的,所以融雪洪水的大小与积雪覆盖面积和积雪深以及影响积雪消融的热量因素有关,融水洪水出现的时间与积雪消融区的高程有关。这类洪水一般峰低量小,造成灾害的大小主要取决于积雪量和升温幅度,若前期积雪量大,且升温幅度也大,则易于造成大的灾害。

4.1.2.2 高山冰川消融洪水

塔里木河源流山区分布着众多的现代冰川和永久性积雪,高山冰川消融洪水是指夏季随气温升高及零度等温线不断上升,高山区的冰川和永久积雪开始强烈消融补给河流,形成冰雪融雪洪水。因此,高山冰雪融雪洪水大小不仅取决于冰雪储量、面积及其分布情况,还取决于太阳辐射和气温等热量因素。高山区5—9月一般为冰川、永久积雪融水期,同时是高山气温由负转正的时段,特别是当零度等温线上升至4000～5000m并持续相当长的一段时间,造成大量冰川积雪融水补给河流,在持续高温天气下,冰川强烈消融期形成冰川消融型洪水。

高山冰川消融洪水的特点是日变化规律显著,持续时间长,变化缓慢,洪量大而洪峰不高,对下游造成的灾害相对较大。冰川消融型洪水的洪水过程线呈锯齿状,是与日气温变化呈同步关系的峰、谷交替变化过程,场次洪水过程表现为基流渐次增加再减少的台阶式复峰形态(图4-3)。常与暴雨洪水形成混合型洪水,在冰川消融洪水发生期内,若出现强降雨过程,消融型洪水与暴雨洪水遭遇形成混合型洪水,特别是在消融洪水消退初期发生遭遇,极易形成较高洪峰,如1999年7月末由于持续升温引发了托什干河消融型洪水过程(图4-4),8月1日局地暴雨洪水与消融型洪水相遇,洪峰流量高达1250m3/s,洪峰流量是起涨流量的8.2倍,洪水历时380h,成为当年最大洪峰。

图4-3　昆马力克河1973年8月高山冰川消融型洪水过程

图4-4　托什干河1999年7月高山冰川消融型洪水过程

4.1.3　混合型洪水

塔里木河流域河川径流补给形式多样,从高山至山前平原区,要经过几种不同的垂直带,不同高度带补给水源形成的不同类型洪水,遭遇混合在一起所形成的混合型洪水。主要的混合型洪水是由暴雨洪水和融雪洪水相互叠加在一起形成的,雨雪混合洪水形成的条件是:流域积雪、较大的降雨强度和较高的地表气温。按其成因分为季节积雪融水与暴雨混合洪水、高山冰雪融水与暴雨混合洪水两种情况。

(1)季节积雪融水与暴雨混合洪水。由于在春季或初夏气温升高,中、低山带季节性积雪大量消融使河流流量增加,若适遇暴雨则加速积雪消融,同时洪峰流量也相应增加。这类混合型洪水历时较长,洪量较大,通常在融雪洪水的基础上,有一个尖瘦的暴雨洪峰,日变化过程不明显,如1963年托什干河,由于上年冬季降雪量多,春季气温偏低,入夏后气温升高很快,积雪大量融化,5月19、20日又有一次大的降水过程,产生春季较大的混合型洪水(图4-5)。

图4-5　托什干河1963年5月混合型洪水过程

(2)高山冰雪融水与暴雨混合洪水大多数发生在夏季,在前期持续高温天气下,高山冰雪消融形成洪水下泄至中低山地带,如适遇暴雨洪水汇合形成混合型洪水,如托什干河1963年7月、8月初出现持续高温,8月为历年当月降雨量最多的年份,为122.7mm,在8月22日、23日两天降雨量分别为10.7mm、23mm,8日23时洪峰流量达到459m3/s,洪峰流量是起涨流量的2.6倍,洪峰历时28小时(图4-6)。

4.1.4　冰川湖突发洪水

由于冰川活动阻塞河谷和高山冰川的消融形成冰川阻塞湖,冰川阻塞湖的突然排水或崩溃而引发所形成突发性洪水。在国际上称为:Glacier Dammed Lake Outburst Flood,或Jokulhlaup(冰岛语,意为溃决、暴发、突发性洪水);国内命名为:冰川阻塞湖引起的突发性洪水,简称冰川湖突发洪水。在我国主要发生在西藏和新疆的冰川作用中心区,但两地的冰川湖类型有所不同,在西藏主要为冰碛阻塞湖的漫溢或溃决造成,新疆主要为冰川阻塞湖冰下、冰内的突然排水形成。

图4-6　托什干河1963年8月混合型洪水过程

冰川湖突发洪水的排水机制是冰川坝在湖水巨大静水压力作用下浮起而使湖水排出。根据水力学原理,此种情况一般在水深超过冰坝高的9/10时发生,同时在静水压力及冰川流动所产生的剪切应力作用和热力作用下,湖水使冰面产生裂隙和断裂向外扩张使冰内、冰下水道扩大而突然排水,继后因冰川塑性变形,导致排水道的收缩率增大,直至大于湖水热侵蚀的扩大率时为止。因为冰川排水道的收缩率是以时间的9次幂速度闭合的,故排水量减少很快,这就是这类洪水落水过程非常迅速的原因。

经过夏季冰川消融期的积累,冰川阻塞湖的蓄水量不断增加,直至静水压力足够大造成冰川坝溃决,引起洪水暴发。这类洪水涨落迅猛,洪峰高,突发性强,难以预测预防,破坏性极大。在新疆境内主要发生在喀喇昆仑山乔戈里峰冰川作用区、天山托木尔峰冰川作用区和帕米尔公格尔峰冰川作用区,在塔里木河流域冰川阻塞湖形成于喀喇昆仑山乔戈里峰冰川作用中心区和托木尔峰冰川作用中心区,而发源于喀喇昆仑山的叶尔羌河和位于天山西部山区的昆马力克河是冰川湖突发性洪水多发河流。这些地区海拔较高,终年积雪,冬季受西伯利亚冷空气影响,积雪大量囤积;夏秋季受副热带高压的影响,冰雪融化,雪线升高;且山脉与盆地之间的海拔落差比较大,使洪水有较大的势能、暴发性地冲灌到低海拔区,形成洪水灾害。

冰川阻塞湖一般有两种情况:一是支冰川阻塞主谷型,是由支流冰川前进阻塞干流河谷而形成;另一是主冰川阻塞支谷型,是由支流冰川与干流冰川交汇,由于支流冰川后退形成冰川湖。

乔戈里峰冰川区的冰川阻塞湖属于支冰川阻塞主谷型,叶尔羌河河源的克勒青河发育有克亚吉尔、特拉木坎力、斯坦格尔、乌尔多克和迦雪布鲁姆等横向冰川,它们均有阻塞干流河谷形成冰川阻塞湖的可能。克亚吉尔冰川阻塞克勒青河谷而形成的克亚吉尔特索湖,其溃决是形成叶尔羌河特大洪水的主要根源;特拉木坎力冰川虽然也堵塞克勒青河河谷,但由于其冰坝下泄水通道发育,使湖泊时有时无。冰川阻塞湖以可塑性冰坝为堤坝,平时没有固定的排水道或排水孔,因而湖水水位随着来水量的增加而升高。当水位达到临界值时,湖水和冰坝之间的静水压力增大,使冰坝下或冰坝内的排水道逐渐扩大而发生溃决从而形成特大洪水。以冰川冰为堤坝的冰川阻塞湖溃决洪水的规模和暴发频率与冰川进退密切相关。叶尔羌河河源克亚吉尔和特拉木坎力等冰川目前处于相对稳定、变薄阶段。预计21世纪30—50年代,气温继续升高,而降水量增加所引起的有利于冰川扩张作用又抵消不了因气温升高所引起的加速冰川退缩作用时,冰川将继续后退和变薄,冰川阻塞溃决洪水的规模将相应减小,暴发频率减低,因而洪水灾害也将日益减轻。

主冰川阻塞支谷型的冰川湖形成于汗腾格里峰—托木尔峰西坡的冰川密集分布区。这里发育有南、北伊诺勒切克两条巨大的树枝状山谷冰川。由于北伊诺勒切克冰川后退与主冰川南伊诺勒切克冰川相分离,南伊诺勒切克冰川阻塞北伊诺勒切克冰川后退空出的谷地,并接纳这两条冰川的融水蓄积成麦兹巴赫冰川湖(现位于吉尔吉斯斯坦共和国境内,距我国境内协合拉水文站约200km)。麦兹巴赫冰川湖突发性洪水暴发频率大大高于叶尔羌河上游的克亚吉尔特索湖,昆马力克河上游冰川湖溃决洪水的发生时间多在冰雪强烈消融时期,常叠加较大的高山冰川消融洪水,一般从发生地点到塔里木河仅需2～3天,巨大的洪量和高频率的暴发,常给下游阿克苏地区带来巨大的经济损失。据记载,受洪水威胁的农田近2×104 hm2,阿瓦提、温宿及阿克苏市均受到洪水的危害,洪水所输送的大量泥沙和碎石,常使水库和灌渠严重淤积,利用系数降低。

冰川湖溃决往往引起较大的洪水,叶尔羌河卡群站和昆马力克河协合拉站均观测到多次冰川湖溃决洪水,从塔河实测最大洪峰流量(1999年8月24日,2280m3/s)看,对塔里木河影响最大的突发性洪水是来自于昆马力河上游麦兹巴赫冰川,现吉尔吉斯斯坦境内的冰川湖溃决洪水。20世纪50年代以来,受全球变暖的影响,与冰川湖相连的北伊诺勒切克冰川不断变薄后退,麦兹巴赫冰川湖的蓄水量逐渐增大,冰川湖突发洪水的频率不断增加。由于南、北伊诺勒切克冰川规模都很大,而且南伊诺勒切克主冰川近几十年一直处于稳定状态,因而冰川阻塞湖还会长期存在,冰川湖溃決洪水的洪峰流量和洪水总量呈增加趋势,造成的危害将更严重。

4.2　塔里木河流域洪水

塔里木河流域洪水由源流山区多种形式的洪水共同叠加而成,主要以暴雨洪水、季节性积雪融化洪水、高山冰川消融洪水,以及叶尔羌河和昆马力克河冰川湖溃决形成的突发性洪水为主,这些洪水特性决定于产流区冰川、积雪规模(面积、储量)、热力条件(高空气温的升幅、积温和持续时间)和各源流来洪的时空分布情况等。暴雨洪水在天山南坡相对较多,昆仑山中低山带亦有出现,暴雨洪水一般表现为峰高、量小、历时短,但因具有突发性,对塔里木河源流区的危害性大,对塔里木河干流一般不会造成大的威胁。

塔里木河流域夏季常处于高压天气系统的控制之下,天气晴朗,光热充足,能提供冰雪融水的热量条件,如遇气温升幅大、高温持续时间长的气候条件,会发生融雪洪水。各源流冰川融水在时间上常具有同步性,由于源流流域范围广、沿程传播距离(河道长度)不同和产汇流条件的差异,洪水汇入塔里木河时一般不相遇,这是造成阿拉尔断面年内发生多次洪峰的主要原因。

经实测资料分析,阿拉尔站的洪水过程呈单峰型或连续多峰型特征。单峰型洪水过程是由某一源流或三源流洪峰遭遇形成,三源流洪水一旦遭遇,将形成塔里木河的较大洪水或大洪水,其特点是洪峰高、洪量大,对塔里木河威胁严重。如1999年8月洪水是阿拉尔站实测历年最大洪水(图4-7),洪峰流量为2280m3/s,主要是由阿克苏河的洪水所致,洪量则是由阿克苏河与和田河洪水共同遭遇形成,该场洪水历时约30d,最大30d洪量达27.42×108 m3。

图4-7　塔里木河1999年单峰型洪水过程

连续多峰型洪水过程是由三源流洪水交错叠加形成,这类洪水的过程线宽,洪峰一般不高,但洪量较大,历时较长,洪水沿程削减相对较少,对塔里木河威胁也较严重。如1978年7月、8月洪水(图4-8),阿拉尔站超过1000m3/s的洪峰发生8次,整个洪水过程历时60d,最大洪峰流量1574m3/s,最大30d洪量26.41×108 m3。

图4-8　塔里木河1978年多峰型洪水过程

阿拉尔站多年平均年径流量为46×108 m3,6—9月是塔里木河的汛期,水量占全年总水量的70%～80%,而春季只占到8%左右。据阿拉尔站1958—2002年45年实测资料统计(表4-1),塔里木河的年最大洪水主要发生在7月下旬—8月上旬,其中7月发生11次,占24.4%,8月发生30次,占66.7%,9月发生3次,占6.7%,其余1次发生在5月,由此可见8月是年最大洪水多发期。阿克苏河年最大洪峰流量主要发生在8月,洪水发生频次为50%,其中86%的年最大洪峰流量集中在7月和8月,集中程度较高(表4-1)。

表4-1　 塔里木河与阿克苏河年最大洪水(1958—2002年)频次统计表

塔里木河自身不产流,主要由发源于冰川区的阿克苏河、叶尔羌河、和田河3条源流补给,经实测资料统计分析(1958—2002年),阿克苏河、和田河和叶尔羌河的洪水平均占阿拉尔洪峰流量比重分别为62.2%、36.3%和1.5%;占3d洪量比重分别为54.3%、44.7%和1%;占7d洪量比重分别为55.3%、43.7%和1%;占15d洪量比重分别为56.6%、42.1%和1.3%;占30d洪量比重分别为59.7%、37.7%和2.6%(表4-2),这表明塔里木河的洪水主要来源于阿克苏河,其次是和田河,叶尔羌河所占比重最小,但叶尔羌河的突发性洪水仍对塔里木河造成影响;随统计时段的增加,阿克苏河洪水比重逐渐增大,和田河洪水比重逐渐减小,这也说明阿克苏河已成为塔里木河水量的主要来源区。

表4-2　塔里木河洪水来源及组成统计表%

另外,从阿拉尔断面的较大洪水与阿克苏河、和田河、叶尔羌河相同年份的洪水过程可以看出,当阿拉尔发生较大洪水时,阿克苏河一般也发生较大洪水,和田河、叶尔羌河不一定发生较大洪水,这也说明阿克苏河是塔里木河洪水的主要源流区(图4-9)。

图4-9　塔里木河和阿克苏河1958—2002年年最大洪峰流量变化图

4.3　阿克苏河流域洪水

4.3.1　托什干河洪水

托什干河流域平均海拔低,集水区面积大,干流沿途汇集多条发源于中低山区的河流,液态降水区范围大,是新疆天山南坡的主要雨区之一,局地降雨易引发暴雨洪水。其主要发生春季融雪洪水、夏秋季冰川消融洪水、暴雨洪水及混合型洪水四种类型的洪水。暴雨洪水发生的频次高,常与融雪型洪水和消融型洪水相遭遇,形成混合型洪水。托什干河的年最大洪峰流量发生在4—8月,其中5月、7月、8月年最大洪峰再现频率高,8月是年最大洪峰流量出现频率最高的月份,出现频率为37%(表4-3)。

表4-3　托什干河1958—2002年最大洪水频次统计表

春季融雪洪水(简称春洪)是托什干河较为典型的洪水类型,其主要特征是发生时间早,一般出现在4月至5月上、中旬,与出现在7月、8月的冰川消融洪水从流量过程线很容易区分开来。托什干河在4月下旬至5月上、中旬都有春洪发生,洪水峰高量大,历时长,洪水规模取决于冬季山区降雪强度(积雪厚度)及翌年春季的气温变化,1957—2002年有26%左右春洪为最大洪峰流量,其中1995年5月出现的特大春洪成为历史最大洪峰流量,春洪多年平均洪峰流量为356m3/s。春洪几乎每年都会发生(图4-10),个别年份春季出现2～3次洪峰,未出现春洪的有9年,即1957年、1958年、1961年、1969年、1971年、1974年、1987年、1998年和2001年,变差系数CV为0.37。春洪形成取决于温度条件和冬季积雪厚度。从阿合奇气象站气温可以看出,出现最大洪峰流量当日气温都较同期温度要高,1995年5月11日历史最大春洪,当日气温高达20.5℃(表4-4)。

图4-10　托什干河历年春季融雪洪水最大洪峰流量

表4-4　托什干河春季融雪洪水频次统计表

春季融雪洪水具有洪量大而洪峰高的特点,与夏洪截然不同,从洪峰过程线形状上看,春洪起涨点低,而夏洪起涨点高;春洪回落点比起涨点高,而夏洪起涨点和回落点基本一致。这是由于在托什干河流域,春洪出现之前即4月上、中旬夜间气温常达0℃以下,日均气温不到10℃,流量一般不超过100m3/s。春洪过后,气温升高,基流相应增大,回落点降不到起涨点的水平上。而夏洪的形成主要依赖于一次大的天气过程或气温的升高,洪水过后流量迅速回落到起涨点。春洪是季节性积雪融水形成的,而夏洪主要由高山冰川消融洪水和暴雨洪水而形成,两者补给来源有差别,所以造成不同的峰值和流量过程线。

4.3.2　昆马力克河洪水

昆马力克河发源于天山最广阔的冰川群——汗腾格里峰冰川区,昆马力克河源流区冰川资源十分丰富,因夏季气温升高而产生的冰川消融型洪水是昆马力克河最常见的洪水类型,其发生区域广、持续时间长、缓涨缓落、洪峰不高但洪量较大,最大洪峰流量常在1500m3/s以下,最大洪峰流量与起始流量的平均比率为3.54,洪水日变化过程显著。昆马力克河上游麦兹巴赫冰川湖突然排水造成冰川湖溃决洪水,这类洪水多叠加于消融型洪水之上,洪峰流量高,洪水流量大,洪水过程退水迅速。另外,局地暴雨形成了昆马力克河流域中小型洪水,由于流域平均海拔高,高山区降水多为固态,因而暴雨洪水多由中低山区降雨形成,雨区面积不大,不易形成较大的洪水。

分析昆马力克河45年内(1957—2002年)最大洪水的洪峰流量及洪水特征(表4-5),在主汛期7月、8月出现的最大洪水频率高达83%。年最大洪水是由高山冰川消融洪水和冰川湖溃决洪水两类洪水形成,其中冰川湖溃决洪水频率高达69%,而第三类冰川湖溃决洪水(叠加于消融洪水之上的冰川湖溃决洪水)频率为52%,是昆马力克河最主要的洪水类型。

表4-5　昆马力克河1957—2002年最大洪水频次统计表

4.3.2.1 冰川湖溃决洪水

麦兹巴赫湖是昆马力克河源流区最大的高山冰川湖,位于托木尔峰西北30km吉尔吉斯斯坦共和国境内,地理位置在东经79°52′,北纬42°13′,位于昆马力克河主要支流依诺勒切克河上游的北依诺勒切克冰川。1902年德国人G.麦兹巴赫对北伊诺勒切克冰川进行考察时发现并命名。麦兹巴赫冰川湖湖面海拔3600m,在协合拉水文观测站上游约208km处,目前最大蓄水量4.3×108 m3,高水位时湖泊长4.5km,平均宽度为1.5km,冰坝出口附近最大水深130m。

麦兹巴赫冰川湖及其周围发育的冰川群成为昆马力克河源流的重要补给源,这里发育有南、北伊诺勒切克两条巨大的树枝状山谷冰川,南伊诺勒切克冰川是天山山区最大的冰川,该冰川长6.1×103 km,面积5.7×104 km2。北伊诺勒切克冰川长4.1×103 km,面积2.4×104 km2。高山冰川湖位于距冰川末端的消融区中部,由于北伊诺勒切克冰川退缩与南伊诺勒切克主冰川相分离,在空出的冰川谷地中形成的主冰川阻塞支谷型冰川湖,补给水量主要来自北伊诺勒切克和南伊诺勒切克冰川。南、北伊诺勒切克冰川属大陆型冰川,冰面裂隙相当发育,冰面融水通过冰裂隙、冰井进入冰内。由于热力作用使冰内温度升高,促进冰内融化。同时地热作用也可使冰下消融,形成冰内融水排泄的通道。随着温度的升高,冰面消融量增加,冰内水量随冰面来水量的增大而加大,冰内水压超过通道所能承受的最大压力时,通道溃决,迅速排水,正常的水文过程被破坏,产生突发性、历时短、洪量大的冰川湖溃決洪水。

麦兹巴赫湖突发洪水观测历史可以追溯到1902年G.麦茨巴赫发现冰川湖时;新疆维吾尔自治区水利厅阿克苏水文水资源勘测大队调查历史洪水后证实,1906年该地区发生了洪峰流量为2520m3/s的大洪水,在1932—1956年麦茨巴赫湖突发洪水14次;据协合拉站水文资料,1957—2002年发生各类冰川湖溃决洪水共计43次(表4-6),频率高达91%。

冰川湖溃决洪水的特点是洪峰流量高、洪水总量大、爆发频率高,从5月至12月都有可能发生,是昆马力克河、阿克苏河乃至塔里木河最主要的灾害性洪水。冰川湖溃决洪水形成机制相同,但因其发生时间不同而在洪峰流量、洪水历时、涨洪历时、落洪历时、起始流量、洪峰流量与起始流量的比率等方面表现出不同的特征,因而可以将冰川湖溃决洪水分为3种类型(表4-6)。

表4-6　昆马力克河冰川湖溃决洪水统计特征参数(值区/平均值)表

第一类为冰川湖溃决洪水发生在非汛期的5—6月或10—12月,特点是洪峰流量不大,此时河道基流一般为30～100m3/s,洪峰流量通常是基流的10倍以上;洪水历时均在10d以上,涨洪历时为8～15d,落洪历时为2～5d。

第二类冰川湖溃决洪水主要发生在主汛期后的8月、9月,特点是洪水历时短,总历时不超过10d,平均为6.4d;落水极为迅速,最长时间为2.2d,最短仅为0.8d,平均为1.3d;涨洪段通常为3～8d,平均为5.5d。平均洪峰流量为980.6m3/s,起始流量为150～450m3/s,洪峰流量通常是基流的2～6倍。

第三类冰川湖溃决洪水发生在汛期7月、8月,叠加于高山冰川消融洪水之上,洪峰流量大幅度提高、洪量剧增,洪峰流量均为年最大洪峰流量,实测最大洪峰流量达到2700m3/s(1994年7月23日);起始流量比较大,常大于300m3/s;洪水历时最小值仅为3d,最长历时超过30d,取决于冰川湖溃决洪水遭遇消融洪水的3种不同方式:一是在消融洪水涨升段发生冰川湖溃决洪水[图4-11(a)],平均洪水历时为4.8d;二是在主汛期河流发生持续性消融洪水时,在消融洪水的高峰期[图4-11(b)]或其后保持高基流情况下发生冰川湖溃决洪水[图4-11(c)],平均洪水历时为14.1d;三是在消融洪水消退过程的初期发生冰川湖溃决洪水,平均洪水历时为8.0d[图4-11(d)]。

图4-11　昆马力克河第三类冰川湖溃决洪水

4.3.2.2 气温影响分析

20世纪50年代以来,全球气候变暖,冰川处于强烈消融并产生大量融水,冰川湖蓄水规模迅速增加,并引发冰川湖突发洪水频率和洪水总量均有所增加。协合拉站实测日平均气温数据可以看出,1957—2002年,年平均气温的递增率为0.0175℃/年,其中冬季(12月—次年2月)气温递增率最大,为0.0414℃/a,夏、秋、春季的平均气温递增率依次为0.0132℃/a、0.0105℃/a、0.0095℃/a,气温升高是冰川湖溃决洪水频繁发生的主导因素。图4-12中表示了冰川湖溃决洪水1957—2002年的多年变化特征,可以看出冰川湖突发洪水频率有所增加,特别是在1980年以后几乎年年发生,且洪峰流量和洪水总量都有所增加。

图4-12　昆马力克河1957—2002年冰川湖溃决洪水洪峰流量

麦兹巴赫湖冰内冰下排水与冰川湖的水量、水位,阻塞冰坝的特征(高度、宽度、结构及物质组成)和垮坝机制,以及下游地貌、沉积特征等多种因素有关。由于冰川湖的突发性排水常发生在冰川强烈消融期及其后的高水位期,气温升高引起的冰川融水迅猛增加,是冰川湖溃决洪水暴发的主要根源。气候变暖对冰川冰运动有着重大影响,随气温升高,冰床软化变形加大,冰底部滑动加速,促进冰川流动与消融;冰川流速增大,阻塞河谷,促进冰川湖形成;冰内温度升高,冰湖水排水通道更易打开;冰川消融加剧,冰川湖库容增大,洪水发生频率、洪峰流量及洪水总量相应增加。

第一类冰川湖溃决洪水共出现7次,规律性不明显,除了1963年6月7日和1980年5月27日发生的洪水洪峰流量成为当年最大洪峰流量,其余洪峰流量一般不大,大致以6～9年为发生周期。

第二类冰川湖溃决洪水共发生12次,均在1990年以前发生,1990年后还没有发生过,其中有6场洪水为年最大洪峰流量,大致以2～4年为发生周期。

第三类冰川湖溃决洪水出现24次,每场洪水洪峰流量都是当年最大洪峰流量,且在1983年后洪峰流量和暴发频率有明显的增加趋势,平均洪峰流量在1983年前为1281m3/s,而1983年后为1678m3/s,增加了25%;在1983年前以2～5年为发生周期,而1983年后以1～3年为发生周期。

4.3.2.3 干燥度指数影响分析

干燥度指数以水分收支与热量平衡的比值来表征某个地区干湿程度的指标,其倒数称为湿润度指数。干燥度指数近年来成为全球变化研究中经常涉及的气候指标之一。通常用降水量与可能蒸发量的对比关系表征某个地区湿润或干燥程度。Selianinov在1937年提出的干燥度经验公式,利用温度与降水量计算干燥度,原公式中的经验系数为0.10,我国科学家根据国内的实际情况经过大量推算将0.10改为0.16。

式中:K为干燥度;0.16∑t表示可能蒸发量;∑t为日平均温度≥10℃稳定期的积温,℃;r为同一时期的降水量,mm。由式(4-1)可知,当K≥4.0即为干旱。

按修正后的Selianinov公式,近年来流域干燥度指数趋于增加,如1994年7月23日最大冰川湖溃决洪水的年干燥度高达13.64;而最干旱的1997年K值为16.4,当年最大洪水洪峰流量为1900m3/s。从各种洪水类型发生时的年平均燥度指数与年平均气温(表4-7)可以看出:第三类冰川湖溃决洪水出现时干燥程度最大,特别是1983年后,干燥度随流域气温升高而增大;第一类冰川湖溃决洪水发生时对应干燥程度最小,但年平均气温较高,促使在非汛期的5—6月或10—12月内因冰川冰雪消融而引发突发性洪水。

表4-7　昆马力克年最大洪水的洪峰流量与干燥度及气温关系表

昆马力克河流域的最小K值为2.4,最大K值为16.4,多年平均值为7.5,属于干旱地区。以K值为4.0和11.0为临界值划为湿润和干旱年,在干旱年份(K＞11.0)和湿润年份(K＜4.0),冰川湖溃决洪水的洪峰流量(分别为1751m3/s和1470m3/s)和年平均洪峰流量(分别为1601m3/s和1355m3/s)都比较大(图4-13),都远高于其他年份洪水的洪峰流量,这表明在干旱少雨和湿润多雨年份出现大洪水的可能性较大。由于在干旱少雨年份气温偏高,高山冰川冰雪消融量增加,发生大洪水的可能性和洪水洪峰流量增大;而在湿润年份降水偏多,山区降水(雨)频繁发生,冰川湖蓄水量较大,积雪和土壤同样有一定的持水能力,雨强较大的降雨能提供融雪所必要的热量,并可使积雪变松,加速融雪出水过程;同时雨水的冲刷作用促使冰雪尽快融化,雨雪融水一泄而下极易形成较大的洪峰过程,但降水时高山冰川区消融量随着温度的降低而减少,因而出现洪水的洪峰流量低于干旱年份洪水。

图4-13　干燥度与冰川湖溃决洪水洪峰流量和年平均洪峰流量对比图

4.3.3　洪水传播特征

阿克苏河由托什干河和昆马力克河两大支流汇合而成,由于阿克苏河流域覆盖面广,两大支流地形、地貌特征、水流条件、流域形状等形成河流洪水的基本要素之间存在着较大的差异。托什干河主河道与山脉平行,流域形状成长条状,且河道纵坡较缓,流域平均高程相对较低,所以洪水汇流集中条件不及昆马力克河,洪峰流量相对小些。昆马力克河流域形状呈扇形,主河道垂直于山脉,由北而南,河道比较短,而且纵坡陡峻,其左边支流发源于汗腾格里峰与托木尔峰冰川集中作用区,其附近有20多座高于6000m的高峰,所以流域平均高程较托什干河流域高。同时,由于昆马力克河上游存在冰川阻塞湖,洪水主要为冰川湖突发洪水和高山冰川消融洪水,其洪量高于托什干河的夏秋季消融型洪水。

选取托什干河洪水109场和昆马力克河洪水67场的传播过程(到阿克苏河西大桥站)进行分析(表4-8),两大支流洪水沿程传播时间(速度)随洪水类型及洪峰流量的差异而变化。总的趋势为,昆马力克河洪峰流量高于托什干河,传播速度大于托什干河洪水;昆马力克河各种类型洪水,传播速度由大到小依次为第三类冰川湖溃决洪水(Ⅲ)、第二类冰川湖溃決洪水(Ⅱ)、第一类冰川湖溃决洪水(Ⅰ)、消融洪水;托什干河洪水传播速度最快的为暴雨洪水,其次为夏季消融型洪水,春季融雪型洪水洪峰流量不大,传播速度也最小;就夏季高山冰雪消融型洪水的传播速度而言,托什干河洪水大于昆马力克河。

托什干河和昆马力克河同时发生洪水,若两河形成洪峰时刻相差通常不超过24h,两场洪峰可能汇合并且同时到达西大桥(新大河)站。分析沙里桂兰克站、协合拉站、西大桥(新大河)站1957—2002年汛期洪峰过程线,两大支流托什干河和昆马力克河的洪水过程遭遇共39次(表4-9),且洪水遭遇情况较复杂,最多是来自两大支流的消融型洪水,共计25次,其他混合型洪水也有发生。

表4-8　阿克苏河流域洪峰传播特征值统计表

表4-9　阿克苏河流域洪水遭遇类型统计分析表　单位:次

4.3.4　洪水系列一致性分析

在阿克苏河流域区间用水(包括引水、退水和排水量)是影响参证站(沙里桂兰克站、协合拉站和西大桥站)洪水系列一致性的主要因素。由于协合拉站、沙里桂兰克站位于两大支流山口区,断面以上多为游牧区,无大规模灌溉引水,两站径流及洪水过程受人类活动影响极小,因而其洪水系列不存在不一致性问题。而协合拉站及沙里桂兰克站以下冲积平原区是阿克苏河流域的水利开发区,河流沿岸建有多处较大引水枢纽,从理论上讲,这些引水枢纽的运行在一定程度上改变了径流的产生条件,使下游测站存在着一定的洪水序列不一致性问题。但由于20世纪50—70年代多数引水设施简陋,引水量观测资料缺乏或已散失,加之年代久远,难以搜集到可靠的引水资料来还原参证站洪水系列。

从阿克苏地区水利局提供的流域内主要引水枢纽多年平均汛期各月平均引水量分析,相对于场次大洪水而言,区间引水量所占比例不大。另外从西大桥站(阿克苏河)与协合拉站洪峰流量及洪量相关图(图4-14)可知,除了特大洪水年份(1994年、1997年、1999年、2001年),数据点没有明显的时段性系统偏离,这表明自1952年以来人类活动尚未对阿克苏河洪水系列一致性产生显著的影响。

图4-14(一)　西大桥站与沙里桂兰克站、协合拉站洪峰流量和洪量相关图

图4-14(二)　西大桥站与沙里桂兰克站、协合拉站洪峰流量和洪量相关图

分析1957—2002年沙里桂兰克站、协合拉站与西大桥站(新大河)年最大洪峰流量(m3/s),年最大1日、3日、7日、15日和30日洪量(m3)相关性表明,阿克苏河洪水变化特征与昆马力克河有较大的相关性(图4-14),而与托什干河洪水序列相关关系不明显,表明阿克苏河洪水在很大程度上承袭了昆马力克河洪水动态特征。在西大桥站(新大河)与协合拉站洪峰流量及洪量相关图[图4-14(a)]中,洪峰流量或洪水总量都有极大值时发生偏离的现象,如洪峰流量出现的特大洪水年份1994年、1999年和2001年出现偏离;2001年洪峰流量虽高,但洪水总量不大,而1997年洪峰流量不甚高,但洪水总量较大,在最大1日、3日、7日和15日西大桥站(新大河)与协合拉站洪量相关图[图4-14(c)(e)(g)(i)]上,1997年、1999年和1994年出现偏离;在最大30日洪量相关图上无显著的偏离年份[图4-14(k)],并且年最大30日西大桥站(新大河)与协合拉站和沙里桂兰克站洪量相关性[图4-14(l)]最好,分别为92%和72%。