1. 引言
山洪泥石流是指在山谷、陡峭地区由于大雨,暴雪或其他自然灾害引起的,并携带大量沉积物和石块的现象。山洪泥石流具有突发性、不确定性的特点。进入21世纪,全球平均每年有近千人死于山洪灾害,造成的经济损失高达每年463亿美元 [1] 。预警作为应对山洪灾害的防御措施,引起了世界各国的重视。中国从20世纪80年代起开始山洪灾害预警相关研究,并在21世纪初明确了以山洪预警等非工程措施为主的山洪灾害防御和综合治理思路 [2] 。
从上世纪五十年代,计算机和遥感测量技术的结合以及地理信息系统的快速发展,形成了多种风险评估方法,遥感和GIS技术在大范围的洪水淹没灾情分析的应用上也更为成熟可靠 [3] 。赵士鹏 [4] 通过山洪灾害的强频分析,以灾变系数为指标针对闽江上游地区的山洪灾害进行了风险评估。唐川、师玉娥 [5] 在其基础上,提出了一套基于GIS的从数据采集到空间属性数据库建立、评价指标体系选择、预测评价分析及山洪灾害危险性评价与风险区划的技术路线和方法体系。管珉、陈兴旺 [6] 利用该思路,以江西省小流域地理地图为基础,对暴雨气候、地形坡度、河网分布等因子进行分析,得到江西省山洪灾害危险性评价图,并将其与易损性评价图进行叠加,完成江西省山洪灾害风险区划研究。唐川 [7] 、朱静 [8] 在对云南省文山县城进行城市山洪灾害易损性分析、风险评价时引入了山洪淹没范围、淹没水深两个指标,提高了山洪易损性、预期损失评估的准确性,进一步完善了山洪风险评价的指标体系。林孝松等 [9] 构建了由微地貌类型、地形坡位、汇流累积量和植被覆盖度等组成的山区镇域山洪危险分区指标体系,并用于重庆市跳石镇山洪灾害危险性分区。许小华、何雯 [10] 将土地利用作为一个风险评估的新因子,得到了更加合理可靠的江西省山洪灾害风险等级划分。周毅 [11] 利用HEC-RAS、ArcView GIS耦合HEC-GeoRAS,对疏勒河流域的昌马水库溃坝洪水演进进行了模拟。霍艾迪 [12] 采用单元分类抽样法利用GIS技术做出地质灾害敏感性评价。
本文通过对河南省灵宝市大湖峪山洪泥石流沟进行淹没分析,在水文模拟基础上,用HEC-RAS、HEC-GeoRAS耦合ArcGIS对山洪泥石流进行分析,实现了山洪泥石流淹没范围的三维显示,获得泥石流淹没水深及淹没面积。进而对大湖峪主峪河道及其支沟进行风险评估,所得的泥石流风险评估结果为泥石流的防治和预警提供了科学依据。
2. 研究区概况及数据准备
研究区概况
大湖峪位于陕西省和河南省两省交界处的小秦岭金矿带东段的北麓,海拔680~2120 m,相对高差1440 m。是河南省灵宝市阳平河的山区河流,为常年性河流。主干水系为黄河一级支流大湖峪河,呈南北走向,发源于小秦岭山地并向北汇入黄河。研究区大湖峪河流总长约10.61 km,主沟平均纵比降135.7%。大湖峪属于暖温带大陆性季风气候,研究区年平均气温13.8℃。多年平均降水量为645.8 mm,最大年降水量984.7 mm (1958年),最小年降水量318.7 mm (1997年)。降水多集中在7月、8月、9月,这3个月的降水量可占全年降水量的50.8%,且多为暴雨。由于河流均处于山区上游且多数为季节性河流,所以旱季流量很小甚至干枯,雨季水量变化很大。多年平均蒸发量1616.4 mm,年最大蒸发量1972.2 mm,年最小蒸发量1221 mm。山势陡峭,沟床比降大,平均都在10%~15.2%。地形有利于坡面径流迅速汇集,具有强大的侵蚀和搬运能力,每逢暴雨,容易形成山洪泥石流。
3. 研究方法和数据来源
3.1. 淹没分析方法
本研究由DEM数据利用SWAT模型把研究区域划分成一系列子流域,根据实测气象、径流数据及流域土地利用图、土壤类型图构建研究区域降水-径流模型,并对该模型进行校准和敏感性分析,以日为时间步长模拟该流域径流量,根据气象部门的降水预报,对该流域未来的径流量进行模拟预报研究。
根据HEC-RAS软件,由降雨预测的径流推求各断面的水位,得到各流域不同控制断面的水位高程。对于淹没范围研究,本文根据有源淹没,建立基于GIS的适合山区丘陵地区的淹没范围模型。利用ArcGIS中的Spatial Analysis/Tope to raster对河道控制断面上的水位进行空间插值计算,生成符合实际的洪水高程面。运用GIS中的栅格计算器,计算实际淹没水深及淹没范围,淹没水深 = 空间插值洪水高程面 − 底层DEM高程,并对计算水深结果进行分类,分为大于0和小于0,大于0的部分即为此洪水位高程时的洪水淹没范围水深分布。把水深图层与研究区的卫星图像相叠加,用GIS的可视化功能,可以直观的显示淹没区范围,计算淹没区面积,评估淹没区损失。
利用HEC-RAS模型分析河道山洪泥石流需要研究区一系列的几何资料,包括河网结构及走向、断面分布等。河网几何数据通过HEC-GeoRAS耦合ArcView GIS来方便地建立。由于研究区域河流纵比降较大,根据实际情况共布设了24个断面。因洪峰流量往往决定是否会发生山洪泥石流淹没及淹没的最大范围。本研究根据大湖峪泥石流流域的山洪泥石流资料,分别模拟了三十年一遇、六十年一遇、一百二十年一遇的山洪泥石流演进过程。本研究采用河网中河段起始断面的正常水深作为模型边界条件,流态采用急流来进行山洪演进分析,得到山洪的水位水面线图。当HEC-RAS来模拟山洪泥石流计算完成后,根据计算的洪峰流量,输入到HEC-RAS洪水演进模型中进行山洪(泥石流)淹没演进分析,将计算结果导入到GIS中进行山洪(泥石流)淹没范围三维显示,从三维山洪(泥石流)淹没范围图中可以直观的读取淹没范围。该淹没范围的计算可为山洪泥石流风险分析提供基础依据。技术流程图如图1所示:
3.2. 山洪泥石流风险评估方法
泥石流灾害风险评估的总的思路是灾害识别 → 灾害调查 → 风险分析 → 风险评估 [13] 。对泥石流灾害的风险评估首先要确定可能的泥石流的沟道及泥石流可能威胁的范围;在确定评估范围的前提下,应该详细研究泥石流通道,主要包括泥石流基本特征调查及承灾体调查两方面。以详实的调查为依托,从危险性和危害性两方面定性分析,随后进行定量计算泥石流风险;最后根据风险许可标准确定泥石流风险。
灾害风险等级评估方法采用国际通用的灾害风险的定义计算泥石流的风险,计算公式如下:
泥石流灾害风险 = 泥石流发生概率 * 泥石流可能造成的后果 (1)
泥石流是山区独特的一种自然地质现象,其发生必须具备3个基本条件:有陡峭便于集水集物的地形;上游堆积有丰富的松散固体物质;短期内有突然性的大量流水。小秦岭地区已经具备泥石流的地形条件及物源条件。由于小秦岭地区山势陡峻、岩石裸露,对雨水的入渗、截流、调节作用小,地表径流强度大,因此当有暴雨天气时,极易引发研究区山洪、泥石流等地质灾害。通过淹没分析得到的淹没范围结合研究区的降雨条件,即可得到研究区的山洪泥石流发生概率图。基于HEC的山洪泥石流淹没分析,模拟研究区内五个居民点在三种不同频率下泥石流发生时的淹没面积及水深,即可得到研究区山洪泥石流发生的后果图。
泥石流发生概率按表1进行评定,泥石流可能造成的后果按表2进行评定,然后利用上述泥石流灾害风险计算公式计算泥石流灾害的风险,最后通过表3评定泥石流灾害的风险等级。
![](Images/Table_Tmp.jpg)
Table 1. Risk probability of mountain torrents
表1. 泥石流发生概率表
![](Images/Table_Tmp.jpg)
Table 2. Consequences caused by mountain torrents
表2. 泥石流可能造成的后果表
![](Images/Table_Tmp.jpg)
Table 3. Disaster risk classification of mountain torrents
表3. 泥石流灾害风险分级
对灵宝市山洪泥石流的发生进行风险评估。
在对大湖峪地区做出初步评估分析后,通过实地考察灵宝市的自然动力因素及人类活动因素,结合室内淹没分析,采用和大湖峪相同的评估方法和标准对灵宝市山洪泥石流的发生进行较准确的风险评估,为今后该区域防灾减灾提供理论支撑。
4. 结果及分析
4.1. 淹没分析结果
根据计算的洪峰流量,输入到HEC-RAS洪水演进模型中进行山洪(泥石流)淹没演进分析,将计算结果导入到GIS中进行山洪(泥石流)淹没范围三维显示,从三维山洪(泥石流)淹没范围图中可以直观的读取淹没范围。得到三种不同频率下的山洪泥石流淹没水深的三维显示图,淹没面积如表4所示。
![](Images/Table_Tmp.jpg)
Table 4. Inundation is of different frequency of mountain torrents
表4. 研究区不同山洪泥石流频率下淹没面积
将三种频率下山洪泥石流淹没情况对比可得:
1) 六十年一遇山洪泥石流淹没范围较三十年一遇山洪泥石流有所增大,其中蜜蜂湾淹没面积增大142.44 m2;黄家瑶淹没面积增大2734.69 m2;南河村淹没面积增大715.19 m2;大湖村淹没面积增大10565.89 m2;大湖北桥淹没面积增大3129.99 m2。
2) 当发生一百二十年一遇山洪泥石流时,蜜蜂湾淹没面积较六十年一遇泥石流增大682.91 m2;黄家瑶淹没面积增大2101.83 m2;南河村淹没面积增大461.13 m2;大湖村淹没面积增大6472.88 m2;大湖北桥淹没面积增大5222.12 m2。
4.2. 风险分析结果
将泥石流发生概率赋值标准及可能造成的结果赋值标准,分别结合大湖峪实际情况,可得大湖峪山洪泥石流发生概率图(图2)及山洪泥石流可能造成的后果图(图3)。
大湖峪里的矿山大约从上世纪80年代开始进行大规模无序的矿山开采,截止调查日期,采矿坑口数千个,形成矿渣方量约5000 × 104 m3。根据河南省泥石流灾害记载,大湖峪过去曾多次发生泥石流灾害,给当地居民带来了巨大的人身安全威胁和经济损失。
通过山洪泥石流发生概率图(图2)可以看出大湖峪地区存在泥石流灾害形成的重大隐患,且主要发生在大湖峪主沟和大湖峪西峪。一旦大湖峪地区达到泥石流发生的条件,该地区大多数情况下就会发生泥石流。根据图1 b可分析出,当山洪泥石流发生时,将大湖峪主沟和大湖峪西峪造成极严重的后果:通常死亡人数大于100,受威胁人数大于1000,直接经济损失大于1000万元,潜在经济损失大于10,000万元。且该地区不同地方均会受到不同程度的影响。
通过图2、图3,结合公式(1),可得到了大湖峪山洪泥石流风险计算图(如图4),可看出大湖峪地区山洪泥石流发生概率大,而且一旦发生,造成的后果也十分严重。大湖峪主峪道发生山洪泥石流风险度极高,可达25,需立即采取行动以减轻可能的灾害影响。
因为大湖峪地区多山地,地形十分陡峭,地质条件复杂,沟谷发育,地形起伏大。岩石较破碎,夏季多暴雨,地表植被覆盖率较低。所以处于山洪泥石流灾害高发区且泥石流规模较大,对研究区内国民经济建设和人民生命财产造成了极为重大的损失。由图5可得,大湖峪地区北部山区山洪泥石流灾害分级可达V级,有极高风险。
![](//html.hanspub.org/file/18-2950497x11_hanspub.png)
Figure 2. Risk probability of mountain torrents in Dahu valley, Xiaoqinling
图2. 大湖峪山洪泥石流发生概率图
![](//html.hanspub.org/file/18-2950497x12_hanspub.png)
Figure 3. Consequences caused by mountain torrents in Dahu valley, Xiaoqinling
图3. 大湖峪山洪泥石流可能造成的后果图
![](//html.hanspub.org/file/18-2950497x13_hanspub.png)
Figure 4. Risk calculation of mountain torrents in Dahu valley, Xiaoqinling
图4. 大湖峪山洪泥石流风险计算图
![](//html.hanspub.org/file/18-2950497x14_hanspub.png)
Figure 5. Disaster risk classification of mountain torrents in Dahu valley, Xiaoqinling
图5. 大湖峪山洪泥石流灾害风险分级图
根据上述公式、表格、结论,及大湖峪的具体情况来看,大湖峪矿渣的随意堆放且与地表物质没有任何粘结关系,在暴雨情况下势必会发生泥石流灾害。可以初步判断大湖峪矿区泥石流隐患相当严重,其发生泥石流的物源条件、地形条件都已具备。一旦水源情况达到,泥石流就会危及到流域下游居民区的生命及财产安全。这一风险评估结果与大湖峪历史发生的泥石流灾害和灾情的实际情况基本相符。
在对大湖峪地区做出初步评估分析后,采用相同的评估方法和标准对灵宝市山洪泥石流的发生进行风险评估,得到灵宝市山洪泥石流风险分级图(图6)。
![](//html.hanspub.org/file/18-2950497x15_hanspub.png)
Figure 6. Disaster risk classification of mountain torrents in Linbao city
图6. 灵宝市山洪泥石流风险分级图
灵宝市泥石流灾害高易发区分布在西峪–大湖峪–武家山以南小秦岭中山,区内共有泥石流沟22条 [14] 。山洪泥石流灾害分级可达V级,危害程度一般为重大级、特大级,破坏性强,一旦发生将给人民生命财产造成极大的损失。
通过实地调查知,灵宝市除了有形成泥石流必备的水源条件、地形条件等自然动力因素外,采矿及交通设施建设对地质环境产生了严重的影响。采矿活动不仅使山(岩)体受到破坏,增加了不稳定因素,而且产生的大量废渣堆积在坑口的沟河之内,严重堵塞了沟道,成为泥石流灾害形成的主要物源。当有暴雨情况时,需要立即行动以减轻可能灾害的影响。
5. 结论
针对大湖峪典型山洪泥石流区的特点、地形状况来研究山洪泥石流演进状况,实现了三维淹没分析与灾害风险评估的结合,为山洪泥石流预报研究和风险评估提供了新的思路:通过对山洪泥石流地区进行淹没分析,在水文模拟基础上,用HEC-RAS、HEC-GeoRAS耦合ArcGIS模型对山洪泥石流进行分析,实现了山洪泥石流淹没范围的三维显示,获得泥石流淹没水深及淹没面积。通过分析并结合实际进而对泥石流地区进行风险评估,所得的泥石流风险评估结果为泥石流的防治和预警提供了科学依据。
通过野外实地调查、室内数据分析,对研究区大湖峪矿渣型泥石流进行了深入的研究探讨,得到以下结论:
1) 本研究对不同频率的山洪泥石流进行演进模拟,模拟结果显示五处居民点均有风险。其中黄家瑶、蜜蜂湾在60年一遇以上山洪泥石流时有了明显淹没,流域下游部分的大湖村和大湖北桥淹没水深和淹没范围都比较严重,对整个村庄造成的灾害损失也会比较大。
2) 大湖峪矿区泥石流隐患已经相当严重,属于具有高风险的泥石流沟。其发生泥石流的物源条件、地形条件都已具备。一旦水源情况达到,泥石流就会危及到流域下游居民区的生命及财产安全。这一风险评估结果与大湖峪历史发生的泥石流灾害和灾情的实际情况基本相符。
灵宝市也处于山洪泥石流灾害高发区,具有极高风险,一旦发生泥石流,会对研究区内国民经济建设和人民生命财产造成极为重大的损失,需要高优先级进行风险防范。
NOTES
*第一作者。
#通讯作者。