邹嘉琪, 权赫春
( 延边大学 工学院, 吉林 延吉 133002 )
随着地球气候的变化和人类活动的影响,洪水灾害日益频发,给人类社会带来了严重威胁.延边朝鲜族自治州因区内河网密布、地形复杂、降雨集中,近年来也频遭洪水侵害,造成了重大经济损失.例如2016年延边州因台风“狮子山”遭受重大洪灾,经济损失高达20亿元,受灾人口约10万人.
近年来,3S(RS、GIS、GPS的统称)技术因其具有强大的空间分析功能和众多的数据管理工具,以及直观的表现手段,逐渐被利用在洪水灾害的研究中,并取得了较好的研究成果[1].例如:Shanker等将导致洪灾的降水和气象因子以及人类活动等归纳为洪水风险影响因子,利用GIS和HEC - 6水力模型,对Pari河流域的洪水灾害进行了评价[2];傅春梅以重庆市沙坪坝区陈家桥为例,利用GIS空间分析和遥感技术制作出该地区洪水灾害危险性和易损性分布图,并将两图加权叠加得到研究区的洪水灾害风险性成果图[3];何报寅等以湖北省为例,通过对构成洪水灾害危险性的4个因子(降水、地形、河网以和历史上洪灾发生的频次)进行分析,得到了湖北省洪涝灾害危险性评价图[4].图们江流域作为延边州最大的水系,目前学者利用3S技术对其洪水灾害研究的还相对较少,如冯恒栋等通过模拟洪水水位条件,得出了不同水位下的洪水淹没范围[5].基于此,本文利用层次分析法和GIS空间分析方法,考虑降雨、地形、水系、水库分布和土壤类型5个因素,对延边州图们江流域洪水灾害危险性进行分析,以此为图们江流域的防灾减灾提供参考.
延边朝鲜族自治州位于吉林省东部的长白山区,北纬41°59′~44°30′,东经127°27′~131°18′.图们江流域主要位于延边朝鲜族自治州内,地处中、俄、朝3国交界,东临日本海[6].本文选取延边州境内图们江流域所在的7个县市(安图县、汪清县、和龙市、龙井市、延吉市、图们市、珲春市)作为研究区.图们江流域的主要支流包括嘎呀河、布尔哈通河、海兰河、珲春河、二道白河、三道白河等.
层次分析法(AHP)是由美国学者T.L.Saaty于20世纪70年代提出的一种多准则决策方法[7],因该方法计算简单、可靠,所以被广泛应用于地质灾害预测中.在运用层次分析法解决问题时,首先需根据问题设定总目标,然后构造层次判断矩阵,最后通过计算判断矩阵的特征值和相应特征向量来获得各因子的权重系数.获得权重系数后,还需利用随机一致性指标和随机一致性比率进行一致性检验,以此判断权重是否合理[8].
洪水的形成主要取决于天气因素和下垫面因素,其中下垫面因素包括地形、植被、土壤、地质、土地利用类型、湖泊及水库的分布等[9].考虑到与洪水灾害危险性的相关性大小和数据的可获取性,选取地形、降雨、河流、土壤类型和水库分布5个因素对洪水灾害危险性进行评价.
根据图们江流域等高线图生成TIN(不规则三角网的简称),然后按30 m×30 m分辨率将TIN转为栅格图,并进行重分类,得到高程重分类后的图层.地形变化利用坡度来表示,用高程相对标准差表示坡度.高程相对标准差用Arcmap中Spatial Analyst模块中的邻域分析计算获得[10],并将地形标准差分为3级: 1级(小于12.5 m), 2级(12.5~37.6 m), 3级(大于37.6 m).由表1可以看出,绝对高程h(m)和高程标准差越小,其因子影响度越大.按照该原则对不同的地形条件进行影响度赋值,得到研究区综合地形影响因子图(图1).
降雨的多少直接影响洪灾的严重程度.根据2005—2015年图们江流域境内各县市的平均降雨量进行分级.分级的原则是:降雨量越大,对洪水灾害的影响度越高,见表2.根据该原则对不同地区赋予不同的影响度值,得到研究区的降雨影响因子图(图2).
表1 综合地形因子影响度
表2 降雨因子影响度
图们江流域内水系众多,因此利用Arcmap中的筛选功能选取研究区内一级支流进行分析.经综合考虑河流级别及其所处的地形,对一级河流建立两级缓冲区,并确定不同缓冲区对应的宽度(表3).河流所处的地形越平坦,洪灾危险性越大.根据文献[4],对非缓冲区赋0.5, 一级缓冲区赋0.9, 二级缓冲区赋0.8, 从而得到研究区综合河流影响因子图(图3).
当土壤中黏粒含量多时,透水性差,易产生地面径流;沙粒含量多时,透水性好,不易产生径流[7].将图们江流域的土壤要素图转为栅格图后,对土壤类型进行划分.划分原则是土壤黏粒含量越多,洪水灾害发生的概率越大.依据该原则对不同土壤类型赋予不同的影响度值,如图4所示.
表3 河流缓冲区等级与宽度值划分标准 m
图1 综合地形影响因子图
图3 综合河流影响因子图
图2 降雨影响因子图
图4 土壤影响因子图
水库属于防护工程.但是,当水库蓄水至一定程度需要开闸泄洪时,离水库越近的地区越危险.因处在不同高程的水库其影响范围不同,故将图们江流域中水库的分布点划分为5个高程等级,并分别建立二级缓冲区.缓冲区的高程划分和缓冲区宽度的选择与河流因子一致.根据水库所在地形越平坦,洪灾危险性越大的原则,对一级、二级与非缓冲区赋0.9、0.8和0.5,如图5所示.
图5 水库影响因子图
H=0.373 6I1+0.248 6I2+0.175 5I3+
0.122 5I4+0.079 8I5.
(1)
表4 构造判断矩阵
根据公式(1)和地形、降雨、河流、土壤类型和水库分布影响度的栅格图,利用ACRGIS中栅格计算器对各图层进行代数叠加,得到图们江流域的洪水灾害的危险性分布图.采用自然断点法,将危险性分为极低、较低、中等、较高和极高5类,由此得到图们江流域洪水灾害危险性的等级分布图,如图6所示.
图6 洪水灾害危险等级分布图
图6中,数值越高的区域越危险,即受洪水灾害的可能性越大.从整体上看,受灾可能性最大的区域位于安图县的二道白河水系和图们市的图们江水系.
1)危险性极高的区域面积约为3 950.69 km2,占研究区总面积的12.29%,主要分布在:安图县的两江镇、二道白河镇和永庆乡(其中以白河村、永红村、西江村等尤为危险),珲春市的敬信镇、板石镇,龙井市的龙江村、老头沟镇,图们市的安山村、农安村、向阳村等.这些地区离水系很近,水库分布集中,而且地势平坦,土壤中黏粒含量多,所以最容易发生洪水灾害.
2)危险性较高的区域面积约为5 066.14 km2,占研究区总面积的15.76%,主要分布在:安图县中北部,如清山村、江源村、金化村、江北村等村镇;龙井市的东平村、铜南村;和龙市的八家子镇、东城镇;图们市的河口村、曲水村;汪清县的东明村、吉祥村;珲春市的马川子乡等.这些地区离水系比较近,绝对高程较低,地势相对平坦,所以发生洪水灾害的可能性也较高.
3)危险性中等的区域面积约为6 297.32 km2,占研究区总面积的19.59%,主要分布在:安图县的环山村、东安村、长白村,延吉市的山根村、三峰村,龙井市的兴南村,汪清县的百草沟镇,珲春市的杨泡乡等地区.这些地区离水系较远,高程标准差相对较大,土壤中沙粒和黏粒含量相当,因此发生洪水灾害的可能性相对较低.
4)危险性较低的区域面积约为8 126.40 km2,占研究区总面积的25.28%,主要分布在安图县南部的奶头村和汪清县西南部的鸡冠乡等.这些地区因远离主干河流,地势相对起伏较大,水库较少,因此洪水灾害性较低.
5)危险性极低的区域面积约为8 705.02 km2,占研究区总面积的27.08%,主要分布在:和龙市的南部和东部,汪清县的中部和东部,延吉市的北部.这些地区因地形起伏大,远离水系,水库较少,土壤中沙粒含量多,因此洪水灾害危险性极低.
本文选取地形、降雨、河流、水库分布和土壤类型5个因子,根据层次分析法计算各因子的加权后,利用GIS空间分析方法得到了图们江流域洪水灾害危险性的等级分布图.研究结果表明,图们江流域洪水灾害的高发区位于安图县的二道白河水系和图们市的图们江水系.该结论与近年来延边州洪水灾害的实际情况基本一致,该结果可为政府相关部门制定防洪减灾方案提供一定参考.洪水灾害危险性分析所涉及的因素众多,本文在研究中仅选取了5个影响因子,因此在今后的研究中,我们将选取更多的影响因子,如防洪工程建设、防洪政策等,以得出更加准确的分析结构.