ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD核心库期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

鄂西山区崩塌落石运动特征及危险性分析以远安县瓦坡崩塌区为例

黄维, 艾东, 胡胜华, 许汇源, 寇磊, 周宝, 袁晶晶, 夏冬生

黄维,艾东,胡胜华,等. 鄂西山区崩塌落石运动特征及危险性分析−以远安县瓦坡崩塌区为例[J]. 中国地质灾害与防治学报,2022,33(6): 37-43. DOI: 10.16031/j.cnki.issn.1003-8035.202107023
引用本文: 黄维,艾东,胡胜华,等. 鄂西山区崩塌落石运动特征及危险性分析−以远安县瓦坡崩塌区为例[J]. 中国地质灾害与防治学报,2022,33(6): 37-43. DOI: 10.16031/j.cnki.issn.1003-8035.202107023
HUANG Wei, AI Dong, HU Shenghua, et al. Characteristics of rockfall trajectory and hazard assessment in western Hubei Province: A case study of the Wapo collapse area in Yuan’an County[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(6): 37-43. DOI: 10.16031/j.cnki.issn.1003-8035.202107023
Citation: HUANG Wei, AI Dong, HU Shenghua, et al. Characteristics of rockfall trajectory and hazard assessment in western Hubei Province: A case study of the Wapo collapse area in Yuan’an County[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(6): 37-43. DOI: 10.16031/j.cnki.issn.1003-8035.202107023

鄂西山区崩塌落石运动特征及危险性分析——以远安县瓦坡崩塌区为例

基金项目: 湖北省地质局科研项目(KJ2022-24);湖北省自然资源厅科研计划项目(ZRZY2022KJ01)
详细信息
    作者简介:

    黄 维(1992- ),男,湖北荆门人,博士,工程师,主要从事地质灾害防治工作。E-mail:huangwei@cug.edu.cn

    通讯作者:

    艾 东(1986- ),男,湖北松滋人,硕士,高级工程师,主要从事地质灾害防治工作。E-mail:aidong216@163.com

  • 中图分类号: P642.21

Characteristics of rockfall trajectory and hazard assessment in western Hubei Province: A case study of the Wapo collapse area in Yuan’an County

  • 摘要: 以远安县瓦坡崩塌区为例,通过地质调查、野外测绘、无人机航拍,建立了瓦坡崩塌区三维模型,基于Rockfall Analyst(RA)分析软件,实现了瓦坡崩塌区大量崩塌落石三维空间下运动路径、高度、能量等要素模拟,探究了崩塌区落石的三维运动特征,开展了崩塌区落石危险性评估,为崩塌落石的防治提供科学参考。研究结果表明:崩塌区危岩主要破坏形式为倾倒式,目前处于欠稳定状态;模拟落石三维运动轨迹与已有落石点基本重合,说明此次模拟结果与实际情况较为吻合;落石运动过程中以碰撞弹跳、自由飞落为主,落石主要集中在崩塌区下部冲沟及公路内,部分落石达到居民区,在崩塌落石区下部公路、冲沟及崩塌区右侧危险性较高;崩塌落石防治工程建议采用危石孤石清除+被动防护网,在公路内侧、斜坡下方分别设置5 m高和3 m高抗冲击力2 000 kJ的被动防护网,可有效拦截落石。
    Abstract: Taking the Wapo collapse area in Yuan’an County as an example, a three-dimensional model was established by drone aerial photography, geological survey and field mapping. Based on Rockfall Analyst (RA) analysis software, a large number of three-dimensional rockfalls in the Wapo collapse area were studied. The motion path, height, energy and other factors in the space were simulated to explore the three-dimensional motion characteristics of rockfalls. The risk assessment of rockfalls in the collapsed area was carried out to guide the prevention and control of the rockfalls. The research results showed that the main failure form of the dangerous rock was the toppling, and it was in an under-stable state at present; the three-dimensional motion trajectory of the simulated rockfall basically coincided with the existing rockfall point which indicating the simulation results were in good agreement with the actual situation; The mode of motion of the rocks were collision, bounce and free fall. Rockfalls were mainly concentrated in the gullies and roads but a few rocks reached residential areas. Therefore, the roads, gullies and the right side of the collapse area were in high risk area; Dangerous rock mass and boulders removal and passive protective nets were adopt as the prevention engineering,two passive protection nets which were 5 m high and 3 m high respectively were set on the inner side of the road and below the slope, and a net with an impact resistance of 2 000 kJ could effectively intercept falling rocks.
  • 崩塌、滑坡地质灾害是一个发生频率高、分布范围广、破坏形式多样、灾害损失巨大的灾害类型。早在20世纪40年代,国外已开始了针对滑坡预警模型的研究,比如1985年美国地质调查局(USGS)和美国气象服务中心(NWS)联合在旧金山海湾地区正式建立了泥石流预警系统,并于1986年2月12—21日在旧金山海湾地区的一次特大暴雨灾害中用于滑坡预报并得到检验[1]。我国香港是世界上最早研究降雨与滑坡关系的地区,于1984年建立了滑坡灾害预警系统[2],国内对地质灾害预警的研究则始于20世纪90年代末,以中国地质环境监测院为代表,在消化吸收国内外学者关于降雨引发区域地质灾害机理等方面的研究成果,并结合中国内地的降雨引发地质灾害分析,以及具体预警业务研究与实践,将区域地质灾害气象预警方法划分为3大类,即隐式统计预报(临界降雨判据法)、显式统计预报和动力预报模型方法[3]。隐式统计预警方法是把除降雨量以外的其他地质环境因素隐含在降雨参数中,仅用单一临界降雨量预警指标来建立某一个地区的预警判据模型来进行预警,对于地质环境条件比较单一的一些小区域是比较适用的,且只涉及一个或一类参数(临界降雨量),易于推广[23],自自然资源部(原国土资源部)、中国气象局于2003年起联合开展地质灾害气象风险预警工作以来,在国内得到广泛应用,2003—2009年,河北省唐山市、福建省泉州市安溪县、山西省太原市、福建省厦门市、河南省郑州市、浙江省杭州市、广东省梅州市等陆续启动市县级地质灾害气象风险预警工作[4]

    临界雨量阈值模型的研究主要包括两种思路,一是基于地质环境条件及气象因素分析,确定临界雨量阈值;二是基于历史地质灾害信息和降雨数据统计,得到经验性的临界雨量阈值。基于条件分析的建模过程复杂、参数难以确定、计算量大,而基于统计的建模过程获取研究数据相对简单,预警精度也较高,从而得到了广泛的应用[45]。基于统计的临界降雨阈值研究又分为降雨强度临界值和降雨过程临界值两个研究方向,降雨强度临界值研究仅根据降雨强度(包括小时降雨强度mm/h和日降雨强度mm/d)进行预警,如中国香港地区;降雨过程临界值的研究则以整个降雨过程中的积累雨量为预警阈值,如重庆市、四川雅安、浙江宁波、山东青岛、云南龙陵县、甘肃省白龙江等地区[310],通过统计分析地质灾害与前期降雨量关系建立降雨阈值模型。

    珠海市地形地质条件较为复杂,局部构造较发育,雨水充沛和降雨时空分布不均匀,近年来极端降雨频发,加上人类工程活动强烈,地质灾害每年均有发生[1112]。2015年建立基于临界降雨阈值模型的地质灾害气象风险预警信息系统(以下简称“预警系统”)并投入使用,对降雨可能引发地质灾害危险区域进行预警,为保障人民群众生命和财产安全发挥了重要支撑作用。新形势对地质灾害防治工作提出了更高要求,加强临界雨量阈值模型研究,提高地质灾害预警精度显得十分迫切和重要。预警系统基于最新的地质灾害调查成果进行了地质灾害预警区划,将引发地质灾害的基础因子(如地形地貌、地层岩性、断裂构造、人类活动等)综合纳入地质灾害预警区划工作中,基于易发性评价结果分级细化预警单元。本文以2013—2021年珠海市发生的342起突发性地质灾害调查资料和对应降雨数据为依据,分析研究突发性地质灾害与降雨的关系,并对降雨诱发地质灾害的特征雨量进行了综合分析,通过拟合曲线,建立基于有效降雨量阈值的E-D预警模型,按细化的预警单元逐一确定预警雨量阈值。经预警效果检验,预警准确率可达到91.5%,单次预警命中率可达94.7%,可为预警系统的阈值模型改进提供技术依据,也为类似花岗岩地区崩塌、滑坡地质灾害预警技术研究提供借鉴方法。

    珠海地貌形态呈哑铃型,地势北高南低,被北东、北西向断裂切割成断块式隆升与沉降的地貌单元,形成了断块隆升山地与沉降平原,山丘台地主要出露燕山期花岗岩,坡度较和缓;平原海岸堆积作用强烈,发育有广阔的冲积海积平原,珠江口外岛屿众多。出露地层较简单,除广泛发育第四系外,在东北部和中西部零星出露有古生代的寒武系、泥盆系和中生代的侏罗系。区内岩浆活动强烈,侵入岩极为发育,大小岩体约80个,出露总面积550.78 km2,占陆域面积的32.4%,其风化岩土体构成山体边坡的主要物质组成。

    2013—2021年,珠海市共发生342起崩塌、滑坡地质灾害,其中滑坡69个,崩塌273个,规模均为小型,绝大多数分布在平原和丘陵交界处(316个,占比92.4%),少数分布在丘陵区域,并与人类活动密切相关,主要发生在公路沿线两侧山体或削坡建房后山陡坡地段。崩塌、滑坡地质灾害分布和主要特征详见图12

    图  1  珠海市崩塌、滑坡地质灾害分布图(2013—2021年)
    Figure  1.  Distribution map of landslides and collapse geological hazards in Zhuhai City (2013—2021)
    图  2  珠海市崩塌、滑坡地质灾害特征图
    Figure  2.  Characteristic map of landslides and collapse geological hazards in Zhuhai City

    珠海市突发性地质灾害主要特点:以崩塌、滑坡为主,其发生与降雨密切相关,具有“小雨小滑、大雨大滑、无雨不滑”的特点,属典型的降雨型崩塌和滑坡,主要发育于斜坡面第四系残坡积层的浅表层小型崩滑破坏,规模虽然不大,但多发生在人员密集区,威胁人民群众生命和财产安全。刘艳辉等[10]通过研究地质灾害与降雨雨型的关系,将诱发地质灾害的降雨雨型分为3种类型:台风降雨、持续强降雨和局地暴雨。珠海市降雨诱发地质灾害属于典型的台风降雨型,地质灾害基本随着降雨过程的发生而瞬时发生。

    珠海市雨量充沛,降雨量在时间上主要分布于汛期(每年4—10月),占年降雨量的80%~90%;年平均暴雨10~11次,均集中在汛期,伴随台风影响雨强大,一次降雨过程可达200~300 mm,易引发崩塌、滑坡地质灾害。

    (1)崩塌、滑坡与月均降雨量的关系

    地质灾害数量与月平均降雨量的关系详见表1

    表  1  崩塌、滑坡灾害数量与月均降雨量关系表
    Table  1.  Relationship between the number of collapses and landslide hazards and average monthly rainfall
    月份 1 2 3 4 5 6 7 8 9 10 11 12
    平均降雨量/mm 26.9 57.9 84.7 199.6 298.1 390.0 319.3 336.7 220.8 71.7 44.3 30.6
    地质灾害数量/个 1 0 5 15 162 56 14 61 24 3 1 0
    下载: 导出CSV 
    | 显示表格

    整体来看,342起崩塌、滑坡地质灾害集中发生在5—9月份,占当年地质灾害数量的92.7%,正是每年的汛期,降雨集中,降雨量和灾害数量的相关度较高,且地质灾害的发生多与强降雨同步,滞后时间很短,这亦是广东地区突发性地质灾害主要特点[13];尤于5月份为最,占比达47.4%,与每年前汛期“龙舟水”雨势强、历时长的降雨密切相关,其中2013年、2014年和2022年最为典型,当年5月的降雨量明显偏多,当月发生的地质灾害数量占比分别高达89.2%、90.9%和89.3%;6、7月份出现降雨量大但灾害数量相对较少的情况,可能系5月份“龙舟水”期间的持续暴雨已诱发了较多的灾害,因此后续存在一个灾害间歇期。

    (2)崩塌、滑坡与当日降雨量的关系

    342个崩塌、滑坡与当日降雨量对应关系详见表2

    表  2  地质灾害数量与当日降雨量关系表
    Table  2.  Relationship between the number of geological hazards and daily rainfall
    当日降雨量/mm00.1~9.910.0~24.925.0~49.950.0~99.9100.0~249.9≥250
    降雨量等级无雨小雨中雨大雨暴雨大暴雨特大暴雨
    地质灾害数量/个1213646810688
    灾害数量占比/%0.30.63.818.719.931.025.7
    下载: 导出CSV 
    | 显示表格

    (3)崩塌、滑坡与前期累计降雨量的关系

    对于当日降雨量小于50 mm的80个崩塌、滑坡灾害点,与前4天(包含当天)累计降雨量对应关系详见表3

    表  3  崩塌、滑坡灾害数量与前4天累计降雨量关系
    Table  3.  Relationship between the number of collapse and landslide hazards and cumulative rainfall over the previous four days
    前4天累计降雨量/mm 0~50 50~100 100~150 150~200 200~250 250~300 300~350 ≥350
    地质灾害数量/个 8 27 18 3 2 4 16 2
    灾害数量占比/% 10.0 33.7 22.5 3.8 2.5 5.0 20.0 2.5
    下载: 导出CSV 
    | 显示表格

    统计结果表明:

    ①珠海市崩塌、滑坡的发生与当日降雨量呈明显正相关关系,342起崩塌、滑坡灾害均由降雨诱发,其中262个灾点发生的当日降雨量超过50 mm,占比76.6%,194个当日降雨量超过100 mm,占比56.7%,当日降雨量越大,发生灾害的数量越多。

    ②崩塌、滑坡的发生还与前4天累计降雨量相关,当日降雨量小于50 mm发生的80个灾害点中,前4天累计降雨量大于50 mm的有72个,占比90.0%。说明崩塌、滑坡的发生与前4天累计降雨量相关度较高,当前4天累计降雨量达到50 mm以上时,发生崩塌、滑坡灾害的概率大。

    ③上述统计结果与珠海市山体覆盖土层普遍较薄,崩塌、滑坡多发生于浅层,基本与强降雨同步,与雨水入渗土体达到饱和状态的时间较短的规律是相吻合的。

    地质灾害预警主要基于地质条件要素与气象因素的耦合来实现,是对不同孕灾条件的预警单元,分析研判在不同降雨条件下发生地质灾害的可能性。因此,区域滑坡降雨预警阈值模型的构建,应当考虑地质环境的差异,区分不同孕灾条件下的降雨阈值[8]

    林春泽等[14]将降雨诱发滑坡临界值模型归结为日降雨量(或降雨强度)和前期降雨量模型,其基本方法是采用统计技术对历史滑坡和降雨资料进行分析,取其统计意义上的临界点作为降雨诱发滑坡的临界值,适用于地质环境模式比较单一的区域预警。

    根据珠海市“台风降雨型”地质灾害特点,当日降雨量是诱发崩塌、滑坡灾害的主要因素,其次为前4天累计降雨量。因此,降雨预警阈值模型的构建应基于降雨强度为主要预警判据,降雨强度可以是某一给定时期内的累计降雨量。根据《地质灾害气象风险预警规范(DZ/T 0449—2023)》,采用累计雨量-持时临界降水阈值模型累计雨量-持时临界降水阈值模型(E-D)适合珠海实际。

    以地质灾害详细调查和评价、监测数据以及成灾规律研究成果等资料为基础,基于致灾地质环境条件差异和地质灾害发育现状结合镇区界线,把珠海市划分为124个预警单元,按照地质灾害易发程度不同分为高易发、中易发、低易发和不易发四个等级(图3)。342个地灾点中,发生在高易发区的有233个,占比68.1%,发生在中易发区的有106个,占比30.1%,发生在低易发区的有3个,占比0.8%。

    图  3  珠海市地质灾害预警区划图
    Figure  3.  Geological hazard early warning zoning map of Zhuhai City

    从降雨历时和降雨强度两个方面研究地质灾害发生时的降雨特征,需采用雨场分割的方法将其分割成不同时段。根据珠海市降雨诱发地质灾害特点,本文以日降雨量为0或者连续2天降雨量小于2 mm作为雨场分割点,分割点之后的降雨时间和降雨量不计入有效降雨历时和有效降雨量。据此对342个地灾点的降雨数据进行雨场分割,计入分割点前的有效降雨历时为1~10天不等,大部分为前4天,占比72.2%,与前述灾害发生和前4天累计降雨量相关度较高基本吻合。

    由于排泄和蒸发作用,前期降雨的影响会随时间减小。前期降雨在扣除地表径流、蒸发等损失后,对崩塌、滑坡发生具有影响的这部分降雨称为前期有效降雨。通过分析前期有效降雨量建立降雨型滑坡模型,最先由Crozier等[15]引入降雨引发的滑坡预报研究,其以滑坡发生当日降雨之前10天的有效降雨量与日降雨量来确定引发滑坡的降雨阈值,浙江省每日预报和短期预报考虑前7日累计有效雨量[16]

    有效降雨量计算公式如下:

    Rc=R0+kR1+k2R2++kiRi++knRn (1)

    式中:Rc——有效降雨量/mm;

    R0——当日降雨量/mm,

    Ri——灾害发生前i天降雨量/mm;

    k——降雨衰减系数(根据相关研究,珠海市降雨 入渗后的表层土含水率衰减系数为0.77~ 0.87之间,平均值0.82,因此,珠海市降雨衰 减系数可取0.82);

    n——有效降雨天数,根据上述雨场分割确定的各地灾点的有效降雨历时为1~10d不等。

    根据各地灾点的有效降雨历时,按式(1),可以计算出342个地灾点的有效降雨量。将342个地灾点样本在直角坐标系上绘制散点图,可以得出地质灾害事件与有效降雨历时和有效降雨量的关系(图4)。

    图  4  地质灾害事件与有效降雨量、降雨历时散点图
    Figure  4.  Scatter diagram of geological hazard events, effective rainfall, and rainfall duration

    根据有效降雨量计算公式,当日降雨量在有效降雨量中得到了充分反映。累计雨量-持时临界降水阈值模型的基本公式如下:

    E=c+αDβ (2)

    式中:E——累计降雨量/(mm·d−1);

    D——降雨持续时间/d;

    α、β——统计得出的参数;

    c——长历时降雨阈值的渐进值,c≥0,根据珠海实际,地质灾害基本随着降雨过程的发生而瞬时发生,主要与降雨量相关,降雨时间的增长对雨量阈值影响甚小,故取c = 0。

    利用散点图,以342个地灾点为基数,采用OLS法的线性模型分别求解10%、15%、30%、55%、75%和90%概率下诱发地质灾害的E-D阈值曲线(图5)。

    图  5  E-D阈值曲线拟合图
    Figure  5.  E-D threshold curve fitting diagram

    E-D阈值曲线中,发生地质灾害的概率为10%的阈值曲线为E =59.37D−0.185,概率为15%的阈值曲线为E =81.51D−0.185,概率为30%的阈值曲线为E= 147.95D−0.185,概率为55%的阈值曲线为E=258.67D−0.185,概率为75%的阈值曲线为E= 347.25D−0.185,概率为90%的阈值曲线为E=413.65D−0.185

    E-D阈值曲线中,根据不同概率的阈值曲线在有效降雨量坐标轴上的截距,可求取不同概率的有效降雨量阈值(表4)。

    表  4  不同概率的有效降雨量阈值
    Table  4.  Effective rainfall thresholds for different probabilities
    发生概率/% 10 15 30 55 75 90
    阈值模型 E =59.37D−0.185 E =81.51D−0.185 E= 147.95D−0.185 E=258.67D−0.185 E= 347.25D−0.185 E=413.65D−0.185
    有效降雨量阈值/mm 59.37 81.51 147.95 258.67 347.25 413.65
    下载: 导出CSV 
    | 显示表格

    根据珠海市降雨诱发崩塌、滑坡灾害的特点和预警经验,按照发生地质灾害概率为10%~15%、30%、55%和75%~90%分别确定四级、三级、二级和一级预警等级的有效降雨量阈值符合珠海实际。

    E-D阈值曲线的拟合是基于已发生的地质灾害数据,据此得出的不同预警等级有效降雨量阈值是针对孕灾条件好的预警区域,即对应于地质灾害高易发区。而地质灾害中易发区和低易发区,诱发灾害的有效降雨量要比高易发区大,从发生在高、中、低易发区的灾点占比可以反映。根据运行多年的珠海市预警系统的降雨量阈值设置和预警经验,以高易发区各预警等级的降雨量阈值为基准,中易发区、低易发区依次降低一个等级,如高易发区一级预警等级的降雨量阈值对应为中易发区的二级预警等级降雨量阈值和低易发区的三级预警等级降雨量阈值,依次类推;四级预警属于关注级别,高、中易发区的阈值按发生地质灾害的概率为10%取值,低易发区的阈值按发生地质灾害的概率为15%取值。综合上述原则并结合珠海实际,确定不同易发程度分区各预警等级的有效降雨量阈值见表5

    表  5  不同易发程度分区各预警等级的有效降雨量阈值
    Table  5.  Effective rainfall thresholds for different early warning levels in various susceptibility zones /mm
    易发程度 预警级别
    四级预警 三级预警 二级预警 一级预警
    高易发区 59.37 147.95 258.67 347.25
    中易发区 59.37 258.67 347.25 413.65
    低易发区 81.51 347.25 413.65
    下载: 导出CSV 
    | 显示表格

    目前地质灾害气象风险预警主流检验方法主要分为两类,第一类为单指标法,以准确率等单个指标综合评估预警效果,第二类为多指标法,评价指标包括命中率、空报率和漏报率[1718]。本文根据历史发生的地质灾害点是否落入预警区范围内从而计算预警准确率以及单次地质灾害预警的命中率、漏报率和空报率对预警雨量阈值模型进行预警效果评价[19]。各评价指标计算公式如下:

    预警准确率:

    p=mn×100% (3)

    式中:p——预警准确率;

    m——落入预警区的地质灾害点数;

    n——总的地质灾害点数。

    命中率:

    P=NANA+NB×100% (4)

    式中:P命中——命中率;

    NA——预警区内灾害点数;

    NA+NB——研究区范围内灾害点总数。

    漏报率:

    P=NBNA+NB×100% (5)

    式中:P漏报——漏报率;

    NB——预警区外灾害点数;

    NA+NB——同上。

    空报率:

    P=SSAS×100% (6)

    式中:P空报——空报率;

    SSA——某级别预警区内没有灾害发生的预警 单元面积;

    S——预警区总面积。

    根据上述确定的预警雨量阈值,对124个预警单元的临界雨量阈值进行更新,以2021—2023年珠海市发生的153个崩塌、滑坡地质灾害事件所对应的实际降雨数据为依据,通过对每个预警单元进行叠加计算并比较分析,判定发生地质灾害的可能性,从而验证预警准确率,计算结果详见表6

    表  6  地质灾害预警准确率计算表
    Table  6.  Accuracy calculation table for geological hazard early warnings
    预警等级无预警落入四级预警区落入三级预警区落入二级预警区落入一级预警区合计
    高易发区地灾点数量/个81325121674
    占比/%10.817.633.816.221.6100
    中易发区地灾点数量/个5228211066
    占比/%7.633.312.131.815.2100
    低易发区地灾点数量/个01021013
    占比/%076.915.47.70100
    合计地灾点数量/个1345353426153
    占比/%8.529.422.922.217.0100
    下载: 导出CSV 
    | 显示表格

    根据表6结果,总的预警准确率为91.5%,其中落入三、二、一级预警区占比为62.1%;高易发区的预警准确率为89.2%,落入三、二、一级预警区占比为71.6%。

    以2022年5月12日为例进行单次地质灾害预警效果检验。当日珠海市普降大暴雨,局部特大暴雨,当日平均降雨量192.9 mm,叠加5月11日暴雨—大暴雨影响,全市范围当天共发生38起崩塌、滑坡灾害。预警结果显示,5月12日当次地质灾害预警的命中率为94.7%,且全部落入三、二、一级预警区,漏报率为5.3%;空报率为68.6%,空报情况主要出现在四级预警区,且预警区面积也最大,如仅以三、二、一级预警区来计算,空报率为16.2%。地质灾害点分布和预警信息发布图详见图6,预警准确率验证数据详见表7

    图  6  2022年5月12日预警信息发布和地质灾害分布图
    Figure  6.  Distribution map of geological hazards and early warning information issed on May 12,2022
    表  7  2022年5月12日暴雨引发地质灾害预警准确率
    Table  7.  Accuracy of geological hazard Early warning triggered by rainstorm on May 12,2022
    预警等级 预警区外 落入四级预警区 落入三级预警区 落入二级预警区 落入一级预警区 合计
    地灾点数量/个 2.0 0 2.00 22.00 12.00 38.00
    占比/% 5.3 0 5.30 57.90 31.50 100.00
    发布预警区面积/km2 418.86 73.94 112.97 64.70 670.47
    有灾点的预警区面积/km2 0 37.20 108.92 64.70 210.82
    下载: 导出CSV 
    | 显示表格

    按照将预警区划分为网格单元来评价命中率和空报率更为客观。网格单元按1km×1km进行划分,有地灾点落入的网格单元数量与网格单元总数的比值即为命中率,无地灾点落入的网格单元数量与网格单元总数的比值即为空报率。根据统计结果,总的命中率为5.4%,空报率为94.6%,详见表8。空报率高的主要原因系四级预警网格单元数量大,但发生地质灾害概率较低。若以124个预警单元为基数进行计算,则总的命中率为29.0%,空报率为71.0%。

    表  8  2022年5月12日暴雨引发地质灾害命中率
    Table  8.  Hit rate of geological hazards triggered by rainstorm on May 12, 2022
    预警等级 四级预警区 三级预警区 二级预警区 一级预警区 合计
    网格单元数量/个 418 74 113 65 670
    有地灾点的网格单元数量/个 0 2 22 12 36
    占比(命中率)/% 0 2.7 19.5 18.5 5.4
    无地灾点的网格单元数量/个 418 72 91 53 634
    占比(空报率)/% 0 97.3 80.5 81.5 94.6
    下载: 导出CSV 
    | 显示表格

    (1)珠海市突发性地质灾害以崩塌、滑坡为主,其发生与当日降雨量和前4天累计降雨量密切相关,基本与强降雨同步,属典型的降雨型浅表层崩塌、滑坡。

    (2)根据降雨型浅表层崩塌滑坡的特点,采用累计雨量-持时临界降水阈值模型确定预警雨量阈值符合珠海实际。当样本数量足够时,可按不同易发程度分区分别建立E-D阈值曲线,确定不同等级的预警雨量阈值,将进一步提高预警精度。

    (3)命中率和空报率是预警效果评价的两个不同方向指标,要提高预警命中率,须尽量增加预警区域,但同时也提高了空报率。需进一步做好预警区划和预警雨量阈值的分析研究工作,努力提高预警精准度,避免过度预警同时减少空报漏报。

  • 图  1   研究区位置图

    Figure  1.   The geographical location of the study area

    图  2   1号危岩体基本情况

    Figure  2.   Basic information of No.1 dangerous rock mass

    图  3   飞行阶段示意图

    Figure  3.   The schematic diagram of the rockfall fly

    图  4   碰撞阶段示意图

    Figure  4.   The schematic diagram of the rockfall collision

    图  5   滚动过程示意图

    Figure  5.   The schematic diagram of rockfall roll

    图  6   研究区坡面分类

    Figure  6.   Classification of slope surface in the study area

    图  7   崩塌落石三维运动轨迹

    Figure  7.   Three-dimensional trajectory of rockfall

    图  8   落石点2运动轨迹

    Figure  8.   Trajectory of rockfall 2

    图  9   危险性分区图

    Figure  9.   Hazard zoning map of the study area

    图  10   设置防护网后落石运动轨迹

    Figure  10.   Rockfall movement trajectory after setting the protective net

    表  1   三维运动轨迹模拟基本参数

    Table  1   Basic parameters of three-dimensional motion trajectory simulation

    坡面类型法向恢复
    系数
    切向恢复
    系数
    动摩擦角
    /(°)
    最大落石
    质量/kg
    初速度
    /(m·s−1
    基岩(植被
    发育区)
    0.370.83301000 1
    居民生活区0.30.8225
    碎石区0.250.8020
    下载: 导出CSV
  • [1] 黄俊光,张帅. 基于动力有限元法的落石运动轨迹研究[J]. 长江科学院院报,2021,38(11):73 − 79. [HUANG Junguang,ZHANG Shuai. Research on trajectory of rockfall based on dynamic finite element method[J]. Journal of Yangtze River Scientific Research Institute,2021,38(11):73 − 79. (in Chinese with English abstract) DOI: 10.11988/ckyyb.20200720
    [2] 叶四桥,陈洪凯,唐红梅. 危岩落石防治技术体系及其特点[J]. 公路,2010,55(7):80 − 85. [YE Siqiao,CHEN Hongkai,TANG Hongmei. Rockfall mitigation techniques and its characteristics[J]. Highway,2010,55(7):80 − 85. (in Chinese with English abstract)
    [3] 王军义,梁风,彭雄武,等. 基于GIS技术的单体崩塌危险范围评价方法研究[J]. 工程地质学报,2021:1 − 12. [WANG Junyi,LIANG Feng,PENG Xiongwu,et al. Study on the assessment method of single collapse risk range based on GIS technology[J]. Journal of Engineering Geology,2021:1 − 12. (in Chinese with English abstract)
    [4] 何宇航,裴向军,梁靖,等. 基于Rockfall的危岩体危险范围预测及风险评价:以九寨沟景区悬沟危岩体为例[J]. 中国地质灾害与防治学报,2020,31(4):24 − 33. [HE Yuhang,PEI Xiangjun,LIANG Jing,et al. Risk assessment and range prediction of dangerous rockmass based on rockfall:A case study of the Xuangou collapse[J]. The Chinese Journal of Geological Hazard and Control,2020,31(4):24 − 33. (in Chinese with English abstract)
    [5] 孙敬辉,石豫川. 重庆甑子岩崩塌落石动力学特征及危险性分区[J]. 中国地质灾害与防治学报,2019,30(3):6 − 11. [SUN Jinghui,SHI Yuchuan. Dynamics and hazard zoning of collapse and rockfall in Zengziyan,Chongqing[J]. The Chinese Journal of Geological Hazard and Control,2019,30(3):6 − 11. (in Chinese with English abstract)
    [6] 胡厚田. 崩塌与落石[M]. 北京: 中国铁道出版社, 1989

    HU Houtian. Collapse and falling rocks [M]. Beijing: China Railway Publishing House, 1989. (in Chinese)

    [7] 唐红梅,易朋莹. 危岩落石运动路径研究[J]. 重庆建筑大学学报,2003,25(1):17 − 23. [TANG Hongmei,YI Pengying. Research on dangerous rock movement route[J]. Journal of Chongqing Architecture University,2003,25(1):17 − 23. (in Chinese with English abstract)
    [8]

    PALMA B,PARISE M,REICHENBACH P,et al. Rockfall hazard assessment along a road in the Sorrento Peninsula,Campania,southern Italy[J]. Natural Hazards,2012,61(1):187 − 201. DOI: 10.1007/s11069-011-9899-0

    [9]

    LAN H X,MARTIN C D,ZHOU C H,et al. Rockfall hazard analysis using LiDAR and spatial modeling[J]. Geomorphology,2010,118(1/2):213 − 223.

    [10] 王学良,张路青,张中俭,等. 石经山藏经洞坡体滚石灾害危险性分析[J]. 岩土力学,2012,33(1):191 − 196. [WANG Xueliang,ZHANG Luqing,ZHANG Zhongjian,et al. Rockfall hazard analysis of slope at sutra caves of Shijing Mountain[J]. Rock and Soil Mechanics,2012,33(1):191 − 196. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2012.01.030
    [11]

    CHEN G Q,ZHENG L,ZHANG Y B,et al. Numerical simulation in rockfall analysis:A close comparison of 2-D and 3-D DDA[J]. Rock Mechanics and Rock Engineering,2013,46(3):527 − 541. DOI: 10.1007/s00603-012-0360-9

    [12]

    LAN H X,MARTIN C D,LIM C H. RockFall analyst:A GIS extension for three-dimensional and spatially distributed rockfall hazard modeling[J]. Computers & Geosciences,2007,33(2):262 − 279.

    [13] 仉义星,兰恒星,李郎平,等. 综合统计模型和物理模型的地质灾害精细评估:以福建省龙山社区为例[J]. 工程地质学报,2019,27(3):608 − 622. [ZHANG Yixing,LAN Hengxing,LI Langping,et al. Combining statistical model and physical model for refined assessment of geological disaster:A case study of Longshan community in Fujian Province[J]. Journal of Engineering Geology,2019,27(3):608 − 622. (in Chinese with English abstract)
    [14] 王栋,王剑锋,李天斌,等. 西南山区某铁路隧道口高位落石三维运动特征分析[J]. 地质力学学报,2021,27(1):96 − 104. [WANG Dong,WANG Jianfeng,LI Tianbin,et al. Analysis of three-dimensional movement characteristics of rockfall:A case study at a railway tunnel entrance in the southwestern mountainous area,China[J]. Journal of Geomechanics,2021,27(1):96 − 104. (in Chinese with English abstract)
    [15] 刘洪江,兰恒星. “5·12”震后都江堰—汶川公路崩塌灾害模拟及危险性评价[J]. 资源科学,2012,34(2):345 − 352. [LIU Hongjiang,LAN Hengxing. Rockfall disaster simulation and risk assessment on the Dujiangyan-Wenchuan highway after “5·12” earthquake[J]. Resources Science,2012,34(2):345 − 352. (in Chinese with English abstract)
  • 期刊类型引用(7)

    1. 蔡小东. 暴雨条件下赣南山区高速公路高边坡稳定性研究. 价值工程. 2025(04): 139-141 . 百度学术
    2. 赖波,江金进,江山,江宁,赵风顺. 珠海市南水镇金龙边坡变形破坏特征及稳定性评价. 城市地质. 2024(01): 20-28 . 百度学术
    3. 王玉刚,袁伟. 基于遗传算法的地下暗挖工程边坡稳定性分析. 大众标准化. 2023(15): 103-105 . 百度学术
    4. 林逸晖,李广涛,杨天雨,乔登攀,王俊,张希,赵怀军. 基于PCA-RF的边坡稳定性预测. 化工矿物与加工. 2023(12): 59-65 . 百度学术
    5. 王猛,何德伟,贾志宏,胡至华. 基于多源遥感数据的高位滑坡特征分析——以广元市利州区荣山镇岩窝村滑坡为例. 中国地质灾害与防治学报. 2023(06): 57-68 . 本站查看
    6. 刘顺青,蔡宇宸,程涛,周萍,王旭畅. 降雨入渗条件下下蜀土边坡稳定性分析. 地质灾害与环境保护. 2022(02): 32-37 . 百度学术
    7. 郭世兴. 关于崩滑地质灾害抢险治理工程的斜坡稳定性分析. 世界有色金属. 2021(21): 203-204 . 百度学术

    其他类型引用(2)

图(10)  /  表(1)
计量
  • 文章访问数:  5254
  • HTML全文浏览量:  2374
  • PDF下载量:  486
  • 被引次数: 9
出版历程
  • 收稿日期:  2021-07-28
  • 修回日期:  2021-12-26
  • 录用日期:  2021-12-26
  • 网络出版日期:  2022-11-10
  • 刊出日期:  2022-12-21

目录

/

返回文章
返回