江西省吉安县地质灾害气象预警研究

2021-12-31张国华何学文

华东地质 2021年4期

张国华,何学文

(江西省地质局第二地质大队,江西九江 332000)

暴雨型地质灾害具有突发性、群发性、区域性等特点,大气降雨或与气象相关的因素成为约三分之二的突发性地质灾害的直接诱因[1],因此,开展地质灾害气象预警工作是抓住地质灾害预警的关键。近年来,众多学者对地质灾害气象预警进行了分析研究[1-13],预警模型由初期仅考虑临界降雨量的隐式统计模型发展到现今综合考虑临界降雨量与地质环境条件等多因素的显式统计模型[10-11,14],提高了预警的精度。2003年以来的实践表明,地质灾害气象预警在防灾减灾中具有较为明显的成效[2-5,9-12,15-16],地质灾害气象预警在当前经济条件下不失为一种宏观层面上简单、经济、有效的防灾手段[11]。地质灾害气象预警是以地质环境条件为基础,根据前期实际降雨量和未来1～3天预报降雨量,对降雨可能诱发的突发性地质灾害在空间、时间和危险性等级上进行预测,并通过电视、电台、互联网、电话、短信等方式向社会公众发出预警。

目前,吉安市气象局已建立了自动气象站雨量实时监测报警平台,实现了吉安市域内国家气象自动观测站、区域自动气象站和山洪自动雨量站降雨的实时监控、资料统计查询、阈值设置和报警等功能[17],为开展暴雨型地质灾害气象预警提供了便利条件,为提高地质灾害气象预警精度提供了保证。本文依托“1∶5万吉安县地质灾害调查”项目,基于吉安县地形地貌、岩土体类型、人类工程活动等地质环境条件,以及江西省气象局提供的降雨量数据,采用信息量法进行吉安县地质灾害易发程度分区,根据降雨量与崩塌、滑坡、泥石流地质灾害的相关性确定地质灾害气象预警临界值,建立吉安县崩塌、滑坡、泥石流地质灾害区域气象预警模型,从而实现地质灾害气象预警。

1 区域环境背景

吉安县位于江西省中部(图1)。地理坐标: 114°24′18″～115°03′25″E, 26°49′50″～27°38′04″N。南北长89.43 km,东西宽64.03 km,总面积2 125.43 km2。铁路、高速公路、国道、省道、县道及乡村道路网状分布,连接县内各乡镇及行政村,交通便利。铁路、高速公路边坡多已采取护坡措施(昌吉赣高铁在吉安境内多为架桥通过),国道、省道护坡路段相对较少,县、乡道路则较少有工程护坡措施。

1.1 气象特征与水力水电建设

据江西省气象局提供资料,1952～2017年吉安县多年平均降雨量为1 524.7 mm。最大年降雨量2 208.7 mm,最小年降雨量897.6 mm;最大月降雨量574 mm(1982年6月);最大日降雨量249.3 mm(2013年7月14日);最大小时降雨量64.8 mm(1998年7月7日)。降雨量年内分配不均。3～8月为雨季,降雨量164.0～238.5 mm;2月、9月为平水期,月降雨量80.1～95.3 mm;10月至次年1月为枯水季,月降雨量49.3～70.7 mm[18]。各月平均降雨量见表1。

表1 吉安县多年月平均降雨量统计结果(1952—2017年)[18]

吉安县降雨量也具有空间分布的不均匀性。表现为山区多年平均降雨量大于河谷平原,县境南东部凤凰—横江一带年均降雨量较小,县境北部—南西部的油田、天河、敖城等地年均降雨量较大(表2)。

表2 吉安县各气象站年平均降雨量统计结果(2008—2017年)[18]

根据2008—2017年降雨量资料[18]统计:全县17个站点年均降雨日148～179天,天河、指阳两站点最多(179天),凤凰站最少(148天);年均暴雨日2～5日,指阳站2日,登龙、小湖村两站5日;连续暴雨最长17～22天,官田、流芳村、鹤州村、小湖村17天,永阳、桐坪22天;连续暴雨最长天数累计降雨量152.1～279.4 mm,凤凰站最小,固江站最大;最大过程降雨量218.2～376.8 mm,澧田站最小,官田站最大;最大24 h暴雨降雨量81.9～238.9 mm,指阳站最小,小湖村最大。

县内水力资源丰富,水利水电设施较完善。据吉安县水利局资料,全县有中型水库5座,小(一)型水库30座,小(二)型水库176座,山塘3 128座。全县有大小引水工程共212处,引水工程以陂坝为主。因水利建设,产生了众多高陡人工边坡,有的位于坝肩,有的在库岸修公路切坡,一些灌溉、小水电引水渠的渗漏也可能构成滑坡诱因。

1.2 地理及地质特征

吉安县处于武功山南翼、罗霄山脉中段与雩山隆褶带交接地带,县境北端、西部和西南侧为山地,由低山、丘陵地形组成,中部和东部为盆地,由低丘、岗地地形和河谷冲积平原构成。全县总体北部、西部和西南侧隆起,山岭起伏,地势高,北部属武功山南翼,主峰海拔580.3 m,西南为罗霄山脉中段,主峰龙山海拔728.7 m,为境内最高点;中部和东部为盆地,地势低洼,海拔一般56～60 m,形成窄长的河谷平原。

吉安县地貌类型有侵蚀构造低山地貌、构造剥蚀高丘陵地貌、构造剥蚀低丘陵地貌、侵蚀堆积河谷平原(阶地)地貌。低山地貌占全县面积的17.0%,高丘陵地貌占全县面积的23.1%,低丘陵约占42.5%,河谷平原地貌约占17.4%。

吉安县岩土体类型有第四系松散岩类、变质岩类、红色碎屑岩类、坚硬—较坚硬一般碎屑岩类、软弱相间的一般碎屑岩类、碳酸盐岩类。境内岩浆岩不发育,仅分布0.668 km2,且位于沟谷坡脚地带,岩浆岩类分布面积太小,几乎可以忽略。

全县已开发利用矿种数11种,主要有砖瓦用页岩、砖瓦用黏土、建筑用石料、煤、铁、水泥用灰岩、高岭土等。全县现有各类矿山企业94个。矿山开采方式主要为露天开采,露采矿山对地质环境的破坏主要是造成地形地貌景观破坏,尤以油田镇为甚。位于天河、安塘、官田等乡镇的煤矿则采取地下开采方式。随着经济发展,城镇建设如火如荼,在城镇建设中,存在较多的依山而建的建筑,工程建设中削坡现象较为普遍,山区农村建房切坡现象较多。

2 地质灾害(隐患)发育特征

2.1 地质灾害(隐患)概况

1∶5万吉安县地质灾害调查项目查明,自20世纪70年代至2019年4月底,全县共发生295处地质灾害,其中滑坡124处、崩塌104处、地面塌陷65处(含岩溶地面塌陷53处、采空地面塌陷12处)、泥石流2处。崩塌隐患172处、滑坡隐患19处、泥石流隐患3处等403处地质灾害隐患点。

1976—2019年，吉安县共发生崩塌、滑坡、泥石流灾害230处,其中只有213处灾害调查到了发生年份,2008—2016年发生了190处,占总数的89.2%。

2.2 地质灾害发育特征

吉安县地质灾害具有如下特征:灾害类型以崩塌、滑坡、地面塌陷为主,泥石流较少;灾害规模以小型为主。崩塌、滑坡、泥石流地质灾害多出现在易于风化破碎的岩土分布区,滑坡以土质滑坡为主,崩塌以岩质崩塌为主;低山区、丘陵区崩塌、滑坡、泥石流地质灾害比较发育;1979—2019年有发生月份记载的地质灾害中,年内发生时间主要集中在4～6月汛期,与降雨量有较好的一致性(图2)。地质灾害与人类工程活动关系密切,81%的地质灾害与人类工程活动有关。

图2 吉安县地质灾害与月均降雨量关系图(1979—2019年)Fig. 2 Relation between geological disasters and average monthly rainfall in Ji’an County(1979—2019)

3 气象预警技术方法

采用显式统计模型法,该方法具体有图层叠加法和指数法(栅格法)[11]两种,本文采用图层叠加法,即以崩塌、滑坡、泥石流地质灾害易发性分区为基础,以降雨作为气象预警分区的主导因素,结合江西省气象局提供的多年的降雨量数据,选取有具体发生时间的地质灾害与气象数据,进行地质灾害与降雨量相关性分析,建立地质灾害时空分布与降雨过程的统计关系,形成吉安县地质灾害气象预警模型。

3.1 地质灾害易发程度区划

以往区划采用的是定性—半定量的袭扰系数法并结合地质分析类比法进行地质灾害易发程度分区[19],为了提高易发程度分区的精度,本次采用基于GIS技术的信息量法定量评价方法进行全县崩塌、滑坡、泥石流地质灾害易发程度分区。信息量模型是基于地质灾害发生的区域规律,采用信息熵的概念来衡量各种影响因素对预测事物不确定性减少的程度。崩塌、滑坡、泥石流灾害预测成功与否是与预测过程中所获取的信息数量和质量有关,用信息量来衡量。信息量越大,地质灾害发生的可能性越大。

首先确定控制和影响地质灾害发生的主要条件。对于地质灾害发育的影响因子选取,不同学者有不同选择[8-9,13,20-21],但对于地层岩性、地质构造、人类工程活动、地形坡度等的选择有较好的一致性。本文通过对吉安县地质灾害及其隐患特征分析认为,地质灾害发育与地形地貌(高程、地形坡度)、地层岩性、交通建设切坡、建房切坡密切相关,地形地貌与地层岩性是地质灾害发育的基础条件,不同的地貌条件和地层岩性条件下,地质灾害的发育程度不同,人类工程活动破坏了原有地质应力平衡,加剧地质灾害发育,是影响地质灾害发育的主要因数。水系与地质构造虽与现有地质灾害发育关系密切度较小,但他们在地质灾害发育中也有不可忽视作用。因此,选取高程、地形坡度、水系、地层岩性、地质构造、交通建设、建房切坡等七大因素作为地质灾害易发程度区划中的地质灾害影响因子。然后按25 m×25 m方格将全县进行栅格化处理,对每个单元格内单因素信息进行量化处理,根据信息量模型先求得各因子信息量值,再将每个单元格内7个因子的信息量叠加得到每单元格内总信息量值,再根据信息量大小,利用统计学自然断点法、易发性信息量分布曲线等方法确定易发区分界信息量值。通过对栅格单元图合并同类项,并考虑类似的地质、自然环境具有类似的地质灾害问题原则,适当进行归并与单元边界线的修改,得到吉安县地质灾害易发程度分区结果。

利用信息量法将吉安县崩塌、滑坡、泥石流地质灾害易发程度按高易发区、中易发区、低易发区、不易发区4个等级进行划分,划分为3个高易发区、11个中易发区、14个低易发区和8个不易发区(图3)。

1.地质灾害高易发区; 2.地质灾害中易发区; 3.地质灾害低易发区; 4.地质灾害不易发区; 5.地质灾害易发程度分区界线; 6.断层; 7.水系; 8.铁路; 9.高速公路; 10.国道; 11.县级行政单位所在地; 12.乡镇所在地; 13.县级境界线图3 吉安县崩塌、滑坡、泥石流地质灾害易发程度分区图Fig. 3 Susceptibility zoning of collapse, landslide, debris flow in Ji’an County

以往利用袭扰系数法并结合地质分析类比法[19]将全县地质灾害划分为易发区、较易发区、少发区3个等级,属于定性的方法,与之相比,本次采用定量法将地质灾害易发区分为4个等级,结果更加细化、精度更高。该方法的地形、地质等基础资料精度有所提高,综合考虑地质灾害本身和形成地质灾害的各种因素,减小了人为对区划结果的干扰,结合此次野外调查结果,对区划结果进行适当的修正,明显提高了区划结果的准确性。

3.2 地质灾害气象预警临界值的确定

根据江西省气象局提供的吉安县气象数据,选取2008—2016年时间段内58处有效数据点进行统计,对地质灾害与降雨量进行相关性分析,建立地质灾害时空分布与降雨过程的统计关系,确定降雨量特征值，建立有效过程降雨量累计值模型,计算临界降雨量特征值(有效降雨量阈值)，最后依据有效降雨量阈值,进行有效降雨量阈值和地质灾害背景的耦合分析,完成吉安县地质灾害气象预警区划。

与地质灾害密切相关的降雨强度特征一般有连续多天过程降雨量和日降雨量。根据吉安县2008—2016年3—8月仅有的30处地质灾害与多天过程降雨量关系数据分析(图4),地质灾害的发生数量有因过程降雨量增大而增多的现象。地质灾害的发生数量与日降雨量关系比较密切,这与吉安县地质灾害基本属于小型规模有关。有关研究表明,江西省小型地质灾害在降雨当日或前期连续降雨的发生率92.7%[22],这与斜坡短直有关,在极端降雨条件下时,短时间内就可使短直斜坡的岩土体饱含水,诱发地质灾害[23]。因此,本文主要选用日降雨量临界值这一参数作为地质灾害气象预警临界值,根据吉安县日降雨量资料,参考江西省地质灾害气象预警日降雨量临界阈值,确定吉安县地质灾害气象预警临界阈值。

图4 吉安县多天过程降雨量与发生灾害数量关系图(2008—2016年)Fig. 4 Relation between rainfall of multi-day and disaster number in Ji’an County(2008—2016)

2008—2016年发生的58处崩塌、滑波、泥石流灾害点发生时日降雨量均已达到50 mm以上,降雨量均达到了暴雨级别。对灾害点按不同降雨量阈值统计情况见表3,按不同区间统计结果见图5—图7。

表3 吉安县不同日降雨量临界值与发生地质灾害关系统计结果(2008—2016年)

图5 吉安县地质灾害与降雨量关系统计方式1(2008—2016年)Fig. 5 Stastistical chart 1 of relationship between geological disasters and rainfall in Ji’an County (2008—2016)

图6 吉安县地质灾害与降雨量关系统计方式2(2008—2016年)Fig. 6 Stastistical chart 2 of relationship between geological disasters and rainfall in Ji’an County (2008—2016)

图7 吉安县地质灾害与降雨量关系统计方式3(2008—2016年)Fig. 7 Stastistical chart 3 of relationship between geological disasters and rainfall in Ji’an County (2008—2016)

三种不同统计方法拐点都明显(图5、图6、图7),但从发展趋势来看,方式1(图5)与方式2(图6)中后段线斜率较前段线段斜率小,反映当降雨量增大时,地质灾害发生的速率反而比前期的小。方式3(图7)按日降雨量<70mm、70～110 mm、≥110 mm统计时,线段斜率随着降雨量增大而斜率增大,反映当降雨量增大时,发生地质灾害现象明显增多,这更符合实际情况,本文采用方式3划分方案。

3.3 地质灾害气象预警等级

参考中国地质环境监测院《地质灾害气象预报预警技术要求》[24]和《江西省突发地质灾害应急预案》(2013年)[25],预警等级由弱到强分为4个等级:四级(蓝色)预警、三级(黄色)预警、二级(橙色)预警、一级(红色)预警。

4 研究成果与认识

4.1 地质灾害气象预警降雨量临界值确定

上文分析表明,当日降雨量达到50 mm、70 mm、110 mm时,发生地质灾害的数量有一个明显的变化,可作为一个降雨量临界值,区内>150 mm或160 mm日降雨量时发生的地质灾害数量资料不多,参考江西省地质灾害气象预警降雨量临界值方案,确定吉安县地质灾害气象预警降雨量临界值方案(表4)。

表4 江西省与吉安县地质灾害气象预警降雨量临界值对比

4.2 降雨量临界值合理性分析

对比分析吉安县地质灾害气象预警降雨量临界值与江西省地质灾害气象预警降雨量临界值,每个降雨强度等级条件下,吉安县降雨量临界值均比江西省的高10 mm。与全省比,吉安县境内基本缺失岩浆岩。资料表明,江西省境内岩浆岩类分布面积广[7],岩浆岩强风化,产物呈砂土状,结构松散,其凝聚力小,遇雨极易失稳,岩浆岩地质灾害发育强度大于其他岩性[23,26]。吉安县缺失岩浆岩,其降雨量临界值比江西省的气象预警降雨临界值高,符合实际情况。另外,从表3及吉安县地质灾害与降雨量关系图看,日降雨量50 mm、70 mm、110 mm临界阈值明显。因此,吉安县地质灾害气象预警降雨量临界值是可信的。

4.3 地质灾害气象预警区划方案与预警响应

综合吉安县崩塌、滑坡、泥石流地质灾害易发程度分区和降雨量临界阈值,将吉安县地质灾害气象预警划分四级(表5):一级预警(红色)、二级预警(橙色)、三级预警(黄色)、四级预警(蓝色)。

表5 吉安县地质灾害气象风险预警等级划分标准

(1)崩塌、滑坡、泥石流地质灾害高易发区内,当气象预报日降雨量达到4级降雨强度阈值(50 mm)时,进入地质灾害气象四级(蓝色)预警,当日降雨量达到3级降雨强度阈值(70 mm)时,进入地质灾害气象三级(黄色)预警,当日降雨量达到2级降雨强度阈值(110 mm)时,进入地质灾害气象二级(橙色)预警,当日降雨量达到1级降雨强度阈值(160 mm)时,进入地质灾害气象一级(红色)预警。

(2)崩塌、滑坡、泥石流地质灾害中易发区内,当气象预报日降雨量达到3级降雨强度阈值(70 mm)时,进入地质灾害气象四级(蓝色)预警,当日降雨量达到2级降雨强度阈值(110 mm)时,进入地质灾害气象三级(黄色)预警,当日降雨量达到1级降雨强度阈值(160 mm)时,进入地质灾害气象二级(橙色)预警。

(3)崩塌、滑坡、泥石流地质灾害低易发区内,当气象预报日降雨量达到2级降雨强度阈值(110 mm)时,进入地质灾害气象四级(蓝色)预警,当日降雨量达到1级降雨强度阈值(160 mm)时,进入地质灾害气象三级(黄色)预警。

(4)崩塌、滑坡、泥石流地质灾害不易发区内不需进行地质灾害预警。

不同崩塌、滑坡、泥石流地质灾害易发区,根据本区未来24 h降雨量预报值按表5进入相应地质灾害气象预警状态。当发布地质灾害气象风险预警时,红色、橙色、黄色、蓝色不同预警级别区按照《江西省突发地质灾害应急预案》(2020年)[27]规定分级启动相应的应急响应行动。在预警区内要加强对地质灾害隐患的监测、排查、宣传等工作,及时采取避让等措施,划定危险区,设立警示标志,组织受威胁群众疏散,保护人民生命财产安全,最大限度减少地质灾害造成的损失。

5 意义

吉安县低山、丘陵区占全县总面积的82.6%,崩塌、滑坡、泥石流地质灾害占全县地质灾害的77.3%,是地质灾害防治的主要对象。随着社会经济的发展,人类工程活动如建房切坡、交通切坡等现象增加,地质灾害隐患数量呈上升态势。

目前,吉安县还没有针对本县自己的一套地质灾害气象预警方案,仅是按照江西省自然资源厅与江西省气象台联合发布的地质灾害气象预警进行地质灾害预报。本文地质灾害气象预警方案的提出,填补了吉安县县域地质灾害气象预警的空白,同时为吉安县地质灾害防治提供了具有针对性和可操作性的方案。

6 问题与建议

1∶5万地质灾害调查表明,多数地质灾害发生时间难以查明,灾害发生时间大多数仅凭群众回忆,准确性低,有确切发生时间的地质灾害数量太少。崩塌、滑坡、泥石流地质灾害的发生除受日降雨量直接影响以外,也与多天过程降雨量密切相关,本文由于灾害与多天过程降雨量的精确时空资料不足,难以分析其相关性,没有研究地质灾害多天过程降雨量临界值。上述情况影响了地质灾害气象预警的精度和准确度。

今后应及时开展地质灾害应急调查,系统收集灾害发生时间、发生前多天过程降雨量和当日24 h降雨量等数据信息;另一方面,收集地质灾害气象预警发布后地质灾害与降雨量资料,做好预警校验,加强地质灾害与降雨强度的关联性分析。为以后地质灾害气象预警降雨量临界值修正积累准确资料,进一步提高地质灾害气象预警准确率。

在进行地质灾害气象预警时同时应注意前期降雨量资料的分析,前期降雨量对地质灾害的发生是有较大影响的,相同的预报降雨量,有没有前期降雨,发生地质灾害的几率是不同的。对于前期降雨量的有效性,单九生等[22]提出了日综合雨量概念,前期每天降雨量删减系数为0.75,前2天是前1天的0.75,前3天是前2天的0.75,时间越长,删减系数越大,在前10天时,删减系数为0.056。将前期降雨量乘以删减系数之和求得前期有效降雨量。