福建省地质灾害气象预警有效降雨模型研究_陈春利.pdf

引用格式：陈春利，方志伟，2023.福建省地质灾害气象预警有效降雨模型研究 J.地质力学学报，29（1）：99110.DOI：10.12090/j.issn.1006-6616.2022090Citation:CHENCL，FANGZW，2023.ResearchonaneffectiverainfallmodelforgeologicaldisasterearlywarninginFujianProvince,ChinaJ.JournalofGeomechanics，29（1）：99110.DOI：10.12090/j.issn.1006-6616.2022090福建省地质灾害气象预警有效降雨模型研究陈春利，方志伟CHENChunli，FANGZhiwei中国地质环境监测院（自然资源部地质灾害技术指导中心），北京100081China Institute of Geo-Environment Monitoring(Consultative Centre of Geo-Hazard Mitigation,MNR),Beijing 100081,ChinaResearch on an effective rainfall model for geological disaster early warning in Fujian Province,ChinaAbstract:Rainfallisoneofthemostimportantexternalfactorsinducinggeologicaldisasters,especiallyprominentinChinasmountainousregions.Manygroup-occurringgeologicaldisastersinducedbyheavyrainfalloccurredinsoutheastChinasmiddleandlowmountainousregions,causingmanycasualtiesandpropertylosses.TakingFujianProvinceasanexample,wecarefullyexaminedtheactualrainfalldataofhistoricalgeologicaldisastercasesandconductedcorrelationstudiesontypicalrainfallprocessesinducinggroup-occurringgeologicaldisastersbasedonquantitativerainfallestimationwithcounty-leveladministrativedistrictsasthestatisticalunit.Wealsocarriedoutapartialcorrelationanalysisbasedonquantitativeindicatorsofthegeologicalenvironmentforverification.TheresultshowsthattheoccurrenceofgeologicaldisastersinFujiancorrelateswithadequaterainfallwithinthreedays,anditisreducedbyareductioncoefficientof0.79daily.Basedonthatfinding,weestablishedaneffectiverainfallmodelforFujianProvinceandverifieditthroughfieldmonitoringandearlywarningapplications.ApplyingthisrainfallmodeltogeologicaldisasterearlywarninginFujiancanreducethewarningarea,lowerthewarninglevel,shortenthewarningduration,andimprovetheaccuracyofgeologicaldisasterearlywarningbymaintainingthehitratio.Theresultsofthestudycanhelptocharacterizetherainfallaccuratelyandcanprovideascientificbasisforthereasonableassessmentofrainfallfactorsinregionalgeologicaldisasterearlywarning.Keywords:group-occurringgeologicaldisaster；southeasternChina；rainfall；effectiverainfallmodel；geologicaldisasterearlywarning摘 要：降雨是诱发地质灾害最主要的外部因素之一，在中国东南沿海中低山丘陵区强降雨诱发的地质灾害集群发生，造成了大量人员伤亡和财产损失。以福建省为例，深入挖掘多年历史地质灾害案例的降雨实况资料，以县级行政区为统计单元，基于定量化降雨估测对诱发群发性地质灾害的典型降雨过程进行相关性研究，并基于地质环境的量化指标进行了偏相关分析验证，建立了福建省有效降雨模型并通过现场监测和预警应用进行了验证。研究表明：福建地区地质灾害的发生与灾害发生 3 日内降雨相关，且以 0.79 为折减系数逐日折减，据此建立了福建地区有效降雨模型；将福建地区有效降雨模型应用于地质灾害气象预警，可以在保证地质灾害命中率的前提下，减少预警区面积、降低预警等级、缩短预警持续时间，提高地质灾害气象预警的精准性。研究结果有助于准确把握降雨特征，可为区域地质灾害气象预警中降雨因子的合理评估提供科学依据。基金项目：国家自然科学基金项目（42207213）；中国地质调查局地质环境监测项目（121210140001500）；中国地质调查局地质调查项目（0001212018CC60020）ThisresearchisfinanciallysupportedbytheNationalNaturalScienceFoundationofChina(GrantNo.42207213),theGeologicalEnvironmentMonitoringProgramofCGS(GrantNo.121210140001500),andtheGeologicalSurveyProgramoftheChinaGeologicalSurvey(GrantNo.0001212018CC60020).第一作者：陈春利（1987），女，博士，高级工程师，主要从事地质灾害成灾机理及预警预报等方面的研究。E-mail：chen_收稿日期：20220715；修回日期：20221114；责任编辑：范二平第29卷第1期地质力学学报Vol.29No.12023年2月JOURNALOFGEOMECHANICSFeb.2023关键词：群发性地质灾害；东南地区；降雨；有效降雨计算模型；地质灾害预警预报中图分类号：P694文献标识码：A文章编号：10066616（2023）01009912DOI：10.12090/j.issn.1006-6616.20220900引言降雨是诱发地质灾害最主要的外部因素之一，强降雨诱发的地质灾害以群发性、突发性为显著特征，严重危害公民生命财产安全（刘传正，2004；周平根，2004）。在充分考虑地形地貌、地层岩性和力学性质等内在因素的条件下，研究地质灾害与降雨事件的关联性，开展区域地质灾害的预警预报，是防范降雨诱发地质灾害的有效途径，并已取得显著成 效（王 雁 林，2005；李 欣，2006；Baum and Godt，2010；魏 平 新 和 李 秀 娟，2015；陈 春 利 等，2016；Piciulloetal.，2018；方然可等，2021）。国内外学者通过地质灾害案例的统计分析总结诱发地质灾害的降雨特征，主要集中在地质灾害与雨型、降雨周期、降雨持时、暴雨持续时间等时空分布规律研究（刘艳辉等，2009,2011），以及降雨诱发地质灾害的阈值研究，包括临界累计雨量阈值(24h、1d、3d、7d、15d 累计雨量等)和临界降雨强度阈值（谢 守 益 等，1995；柳 源，1998；Gladeetal.,2000；Aleotti,2004）。实际上，由于地表径流和蒸腾作用，渗入岩土体的降雨量小于气象记录实况降雨量，因此近年来在降雨诱发滑坡的研究中逐步引入有效降雨量概念(Crozier,1999；李铁锋和丛威青，2006；简文彬等，2013；吴益平等，2014)。在以往的研究中，国内外学者通过折减系数模拟土壤含水量的衰减 过 程，提 出 了 前 期 有 效 降 雨 模 型（CrozierandEyles，1980；Crozier，1986；陈景武，1990；刘传正等，2009），其 中 Crozier 提 出 的 模 型 应 用 最 为 广 泛（CrozierandEyles，1980；Crozier，1986），许多地方的对比检验效果都比较理想（李铁锋和丛威青，2006；吴益平等，2014）。在中国东南沿海和西部山区等地区，岩体风化差异大、坡残积层发育，持续强降雨过程特征明显，降雨诱发的群发性地质灾害群集发生，易造成大量人员伤亡和财产损失，是降雨诱发群发性地质灾害气象预警的重点研究区(韦方强等，2010；刘艳辉和刘丽楠，2016；贺凯等，2016；熊炜等，2018；闫金凯等，2020；李壮等，2020；韩帅等，2022，张浩韦等，2022)。相关学者针对东南沿海地区的研究得到了台风影响第一天地质灾害对当日雨量更敏感，第二天地质灾害表现为受到前期有效雨量和当日雨量共同影响的特征（刘艳辉和刘丽楠，2016），以及福建山区的 13 天降雨是导致发生滑坡的重要降雨区间（伍宇明等，2014）等重要结论，但针对有效降雨模型的研究较少。文章以福建省中低山丘陵区为例，研究地质灾害与有效降雨量的相关性，针对性地研究岩土体入渗降雨量对地质灾害的诱发作用，提出福建省有效降雨模型，为区域地质灾害气象预警中降雨因子的合理评估提供依据。1研究区概况福建省地处华南褶皱系的东部，地势西北高、东南低（图 1），山地丘陵占陆地面积的 90%左右，山坡坡度一般 1030，局部大于 30。地貌类型复杂，按成因可分为断块侵蚀、断块剥蚀、侵蚀溶蚀和断块堆积 4 种基本类型。地层发育较齐全，除志留系、中下泥盆统和下第三系缺失外，自太古界至第四系均有出露。岩性以岩浆岩为主，约占全省陆地面积的 64%；变质岩主要分布在西北地区，面积约占陆地面积 22，主要分布在福建西北地区；沉积岩主要分布于闽西南地区，占全省陆地面积14。区内北北东北东东、北西东西及南北走向断裂十分发育，其中南平宁化（北东）断裂带与N图例高程市级地名水系0200 m200500 m500800 m8002162 m0 25 50100km图1福建省地势图Fig.1TopographicmapofFujian100地质力学学报https:/2023政和大埔（北北东）构造岩浆带将福建省割切为闽西北、闽西南及闽东 3 个地体（韦德光等，1997；图 2）。地表附近多数为高角度断层和节理，花岗岩风化壳广泛分布，坡麓地区残积、残坡积物构成的坡积裙大量发育，为降雨条件下发生浅层滑坡提供了物质基础（何永金，1995；王涛等，2003；李欣，2006）。构造(断裂)带海岸线闽西北地体闽东地体闽西地体南图 例省级行政中心地级市行政中心省级界南平三明龙岩漳州厦门泉州莆田福州宁德0 20 40kmN图2福建省构造分区图（韦德光等，1997）Fig.2TectoniczonesinFujian（Weietal.,1997）气候分区显示福建省地处南亚热带和中亚热带海洋性季风气候区，温暖湿润，雨量丰沛，大多数地区年降雨时间介于 120200 天之间，年暴雨时间介于 3.37.5 天，内陆多于沿海，集中于春、夏两季，春季多梅雨，夏季多台风暴雨，10 月至翌年 2 月为旱 季。据 气 象 观 测 资 料 统 计，全 省 年 均 降 雨 量10362048mm，大于 1800mm 的多雨中心一处在周宁、寿宁、宁德、柘荣、屏南一带，另一处在武夷山、邵武、光泽一带（图 3）。区内水系发育，河流众多，河网密度达 0.1m/km2，水系走向多呈南东东方向，山地性河流居多，河床比降大，水量丰富，流量和水位季节变化显著。地下水赋存类型主要为基岩裂隙水、碳酸盐岩岩溶裂隙水和松散孔隙水。根据福建省地质灾害详细调查成果，全省地质灾害易发区面积达 8.8104km2，占全省陆地面积的75%。由于福建地区濒临西太平洋，降雨量大且相对集中，台风暴雨多，气候、地质环境条件决定了福建省为降雨诱发群发性地质灾害多发区之一(黄俊宝，2013；叶龙珍等，2015；狄靖月等，2019；闫金凯等，2020)。据全国地质灾害通报（中国地质环境监测院，20112020），20112020 年，福建省共发生地质灾害 5343 起，造成 190 人死亡失踪，直接经济损失 5.4 亿元。滑坡、崩塌、泥石流是主要灾害类型，集中发生在 59 月，地质灾害发生数量与月均降雨量耦合性较好（图 4）。台风暴雨、强降雨作用下群发性地质灾害多发（表 1），在福建北部、中西部等地区降雨诱发的地质灾害集群出现（图 5）。0501001502002503000500100015002000250030003500400045001 月3 月5 月7 月9 月11 月降雨量/mm灾害数量/起地质灾害月均降雨图4地质灾害与月平均降雨量关系图（20102019 年）Fig.4Relationship between monthly rainfall and geologicaldisaster(20102019)2福建地区有效降雨模型降雨是引发地质灾害的重要因素，地质灾害的发生不仅与当日雨量密切相关，也存在一定的滞后效应，且随着降雨时间的前移，前期降雨的影响逐2000 mm1800 mm1600 mm1400 mm1200 mm1800 mm2000 mm1600 mm降雨等值线省级界地级市行政中心省级行政中心图 例0 20 40km海岸线南平三明龙岩漳州厦门泉州莆田福州宁德N图3福建省年降雨等值线图Fig.3ContourmapofannualrainfallinFujian第1期陈春利，等：福建省地质灾害气象预警有效降雨模型研究101渐减弱。1980 年，Crozier 在降雨诱发的滑坡预报研究中引入 API 方法，提出用于计算滑坡事件的前期有效降雨模型（公式（1）；Crozier，1986）。许多研究者采用 Crozier 提出的前期有效降雨模型或者对其修改后计算前期有效降雨量（李长江等，2011），但是各自所取的累积降雨影响时段（n）和有效降雨系数（k）差异较大（表 2）。Pa0=ni=1kiRi（1）公式中：Pa0为前期有效降雨量，mm；i 为计算日前期降雨时间，天；Ri为计算日前第 i 日的日降雨量，mm；k 为有效降雨系数；n 为前期有效降雨时间，天。在地质灾害气象风险预警中，以前期有效降雨和当日降雨作为地质灾害诱发因子，建立有效降雨模型是实现预警的基础，而确定研究区累积降雨影响时段（n）和有效降雨系数（k）是准确计算有效降雨的关键，但由于测试条件限制目前多取经验参数。文中在 Crozier 提出的前期有效降雨模型通式的基础上，考虑群发性地质灾害与降雨的耦合作用，通过相关性分析确定研究区累积降雨影响时段（n）和有效降雨系数（k），进而提出有效降雨模型。通过筛选近 10 年福建地区的降雨诱发地质灾害案例，保留有准确发灾时间记录的地质灾害案例460 起，逐日匹配空间分辨率为 5km 的定量化降雨估测数据（QPE），形成 460 个有灾前逐日历史降雨数据的地质灾害样本。以县级行政区为统计单元，当日 8 时至次日 8 时为统计时间段，若该时段内发生地质灾害，则记录为 1 次诱发地质灾害的降雨事件，筛选的 460 起地质灾害统计为 178 个诱发地质灾害的降雨事件，其中发生 3 起以上地质灾害的降雨 事 件 48 起，最 大 规 模 为 单 日 群 发 地 质 灾 害 41起。对持续性降雨诱发群发性地质灾害的过程进行追溯，形成 28 个未引发地质灾害的降雨事件。利用 SPSS 软件分析 206 个降雨事件中地质灾害发生数量与日降雨量的相关性(表 3)，发现福建地区地质灾害发生与前 1 日降雨相关系数最高，其次为当日降雨和前 2 日降雨，前 3 日以及之前的日降雨量与地质灾害数量在 0.05 和 0.0.1 水平上均不满足显著相关。同时考虑内动力条件对地质灾害的控制作用，采用地质灾害潜势度指标（刘艳辉等，2008）量化地质环境条件，并开展偏相关分析（表 4），同样得到地质灾害发生数量与当日降雨、前 1 日降雨和前 2 日降雨相关性较显著，与前 3 日及其以前降雨量相关性不显著的结论。因此，该区域有效降雨计滑坡崩塌泥石流省级界地级市行政中心省级行政中心0 20 40km图 例南平三明龙岩漳州厦门泉州莆田福州宁德N海岸线图5福建省降雨诱发地质灾害分布图（20102019 年）Fig.5Distribution of the rainfall-induced geological disasters inFujian（20102019）表1福建省典型群发性地质灾害发生数量统计表（20102019 年）Table1Group-occurringgeologicaldisastersinFujian（20102019）年份降雨时间降雨情况地质灾害数量/起典型案例20106月13日27日全省大部分县市降雨量超过100mm，59个县（市、区）的733个乡镇降雨量超过250mm，26个县（市、区）的158个乡镇降雨量超过500mm，以建宁樱桃岭的875mm最大1574南平市“614”滑坡灾害（造成24人死亡失踪）20118月29日9月3日台风“南玛都”在福建登陆，福建东部局部地区出现了大暴雨8520128月2日3日台风“苏拉”在福建登陆，北部局部区域和西部大部分区域出现暴雨、大暴雨4120155月18日21日中南部地区出现大范围的大暴雨，局地特大暴雨10720165月7日9日西部的三明、南平两市发生的特大暴雨，多处雨量站点监测实况雨量大于150mm/d，成片区出现大暴雨量级及以上等级的降雨43福建泰宁“58”特大型泥石流灾害（造成36人死亡失踪，直接经济损失600万元20197月6日9日福建大部有暴雨，局地有大暴雨，累计雨量100250mm22102地质力学学报https:/2023算时间为 3 日内，即 n 取 2。在确定有效降雨影响时段为 3 日内的基础上，采用有效降雨模型通式（公式（1）分别计算有效降雨系数在 01 范围内对应的有效降雨量与地质灾害发生数量的相关性，发现相关系数呈先升后降的趋势（图 6），在有效降雨系数 k=0.79 时地质灾害发生数量与有效降雨量相关性最强（表 5）。根据上述分析，确定福建地区有效降雨系数为 0.79，累计降雨影响时段为 3 日内，提出福建地区有效降雨计算模型：P=2i=00.79iRi（2）00.10.20.30.40.50.600.10.20.30.40.50.60.70.80.91.0相关系数有效降雨系数图6有效降雨系数相关系数变化曲线Fig.6Variation curve of effective rainfall coefficient withcorrelationcoefficient表2地质灾害有效降雨量计算模型一览表（CrozierandEyles，1980；李长江等，2011；唐红梅等，2013；狄靖月等，2019；刘艳辉等，2019）Table2Summaryoftheeffectiverainfallcalculationmodelsforgeologicaldisasters(CrozierandEyles,1980;Lietal.,2011;Tangetal.,2013;Dietal.,2019；Liuetal.,2019)时间提出者计算方法衰减系数 计算时间/天参数说明研究区1980CrozierMJEylesRJPa0=10i=1KiRi0.8410Pa0为前期有效降雨量，mm；i为计算时间，d；K为衰减系数；Ri为第i日的日降雨量，mm美国(渥太华)1985濑微克美Pa0=14i=1iRi0.5i/T14iPa0为前期有效降雨量，mm；i为计算时间，d；为衰减系数；Ri为计算日前第i日的日降雨量，mm；T为半衰期日本1991谭万沛Pa0=14i=1iRi0.5i/214中国(四川、云南、西藏4条泥石流沟)1993谢正伦等Pa0=7i=1iRi0.5i/T7中国(台湾)2000中国科学院水利部成都山地灾害与环境研究所Pa0=20i=1KiRi0.820Pa0为前期有效降雨量，mm；i为计算时间，d；K为衰减系数；Ri为第i日的日降雨量，mm中国(云南)2001王礼先等Pa0=20i=1KiRi0.815中国(北京)2004詹钱登等Pa0=7i=1iRi0.87iPa0为前期有效降雨量，mm；i为计算时间，d；为衰减系数；Ri为计算日前第i日的日降雨量，mm中国(台湾)2005ZezerePa0=30i=1iRi0.930葡萄牙2013唐红梅等Pa0=4i=1KiRi0.64Pa0为前期有效降雨量，mm；i为计算时间，d；K为衰减系数；Ri为第i日的日降雨量，mm中国(重庆)2019刘艳辉等Ti=mx=1kxTix0.847Ti为有效雨量，mm；i为预警单元编号；m为有效降雨时间，d；k为衰减系数；Tix为前x天的降雨量，mm；x为预警当日向前计算的时间中国(四川)2019狄靖月等Pa0=10i=0KiRi0.815Pa0为前期有效降雨量，mm；i为计算时间，d；K为衰减系数；Ri为第i日的日降雨量，mm中国(东南地区)表3地质灾害发生数量与日降雨量相关性分析Table3Correlationbetweenthenumberofgeologicaldisastersanddailyrainfall日降雨量/mm当日前1日前2日前3日前4日前5日前6日前7日前8日前9日相关系数0.290*0.463*0.137*0.0660.0460.0570.0490.0510.0580.085注：*在0.01水平（双侧）上显著相关；*在0.05水平（双侧）上显著相关表4地质灾害发生数量与日降雨量偏相关分析Table4Partialcorrelationbetweenthenumberofgeologicaldisastersanddailyrainfall日降雨量/mm当日前1日前2日前3日前4日前5日前6日前7日前8日前9日偏相关系数0.293*0.464*0.137*0.0670.0480.0570.0490.0520.0590.087注：*在0.01水平（双侧）上显著相关；*在0.05水平（双侧）上显著相关第1期陈春利，等：福建省地质灾害气象预警有效降雨模型研究103公式中：P 为有效降雨量，mm；i 为计算日前期降雨时间，天；Ri为计算日前第 i 日的日降雨量，mm。3现场监测验证为验证基于相关性分析建立的有效降雨模型的适用性，通过在福建北部地区开展的降雨入渗现场监测结果辅以验证。监测点位于福建北部宁德市蕉城区石后乡大岭村，在斜坡表层坡积残积土中开 挖 监 测 井 开 展 降 雨 入 渗 现 场 监 测，监 测 期 自2020 年 9 月2021 年 8 月，监测指标包括降雨量、地下水位以及斜坡表层不同深度土体体积含水率（图7）。根据已有研究结果，福建地区降雨诱发地质灾害以浅层土质滑坡为主，通过地下水位监测可开展滑坡预测(何永金,1995；许建聪等,2005;王智磊等,2011；黄俊宝，2013；伍宇明等，2014)，通过斜坡地下水位的上升消落与降雨序列的响应验证有效降雨模型。降雨监测显示，监测期内降雨时间为 197 天，其中出现暴雨 11 次，大雨 8 次，降雨集中出现在 5 月8 月，2021 年 8 月累计降雨量为全年最多，达到了645mm。一次集中强降雨导致短期内地下水位突涨，降雨停止后地下水位急速下降；受斜坡整体水文环境影响，年尺度斜坡地下水位呈现典型的线性下降段（9 月次年 2 月）和上升段（5 月8 月），与旱季雨季交替一致（图 8）。y=0.0115x500.24R2=0.921y=0.0435x+1938.1R2=0.805602040608010012014001234567892020-9-12020-10-42020-11-62020-12-92021-1-112021-2-132021-3-182021-4-202021-5-232021-6-252021-7-292021-8-31降雨量/mm水位距地表深度/m日期图8降雨地下水位关系图Fig.8Diagramoftherainfallandgroundwaterlevel选择与前、后序降雨有间歇期，且引发地下水位显著波动的降雨过程，通过地下水位的上升消散验证有效降雨作用周期。2021 年 5 月 8 日降雨 66mm，导致地下水位上涨 0.93m，该日前 20 日内最大日降雨量 13.3mm，且前 2 日无降雨，当日降雨停止后无降雨直至 3 日后（5 月 12 日）降雨 24mm，但未引起地下水位波动，可排除前序、后续降雨干扰。此时段地下水位呈线性上升，根据趋势模拟无降雨影响下地下水位，降雨后第 3 天实测地下水位与趋势性模拟水位重合，此后走势基本一致（图 9），由此判断单次降雨影响时段为降雨 3 日内，与上文提到的有效降雨影响时段吻合。进 一 步 选 择 降 雨 集 中 的 2021 年 7 月 26 日 8 月 31 日进行验证，由于此前 25 日内无显著降雨（7 月 2 日7 月 26 日）可排除前序降雨影响。自7 月 26 日期逐日分别计算 3 日内累计降雨与有效降雨（公式（2），可见有效降雨和 3 日累计降雨均与地下水位上涨波动趋势一致，但采用公式（2）计算的有效降雨量更能及时反映地下水位消落，如 8 月7 日、8 月 21 日，前一日峰值后地下水位出现消落，012345687深度/m9水位计雨量计数据中心土壤温湿度计监测井图7降雨入渗现场监测示意图Fig.7Schematicdiagramshowingthesitemonitoringofrainfallinfiltration表5地质灾害发生数量与有效降雨量相关性分析Table5Correlationbetweenthenumberofgeologicaldisastersandeffectiverainfall有效降雨系数k=0.77k=0.78k=0.79k=0.80k=0.81k=0.82k=0.83k=0.84相关系数0.53404*0.53408*0.53409*0.53406*0.53401*0.53392*0.53380*0.53365*注：*在0.01水平（双侧）上显著相关104地质力学学报https:/2023计算的有效降雨量也呈现峰值后的大幅下降，而3 日累计降雨量仍保持不变或仅出现小幅下降（图10），可见上文提出的有效降雨模型更能及时反映斜坡稳定性变化趋势。0100200300400500600700800012345677 月 26 日7 月 28 日7 月 30 日8 月 1 日8 月 3 日8 月 5 日8 月 7 日8 月 9 日8 月 11 日8 月 13 日8 月 15 日8 月 17 日8 月 19 日8 月 21 日8 月 23 日8 月 25 日8 月 27 日8 月 29 日8 月 31 日降雨量/mm水位距地表深度/m日期有效降雨累计降雨地下水位图10地下水位降雨量响应关系图Fig.10Responsediagramofthegroundwaterlevelandrainfall4预警应用验证建立有效降雨模型的相关性分析采用福建地区近 10 年（20112019 年）的样本，为进一步验证模型适用性，将该有效降雨模型应用于 2020、2021 年的地质灾害气象预警计算。对强降雨诱发的地质灾害进行筛选，选择 2021 年 6 月 28 日7 月 3 日福建南平地区的降雨致灾过程进行验证分析。此次降雨过程共持续 5 天，累计降雨量超过 300mm，期间在南平市降雨诱发地质灾害 4 起（图 11）。此次降雨过程前 4 日无强降雨，7 月 2 日降雨结束，无前、后序降雨叠加影响，便于开展针对性的分析验证。将研究区划分为精度为 5km 的预警网格单元，采用相匹配的定量化降雨估测数据为输入因子，计算有效降雨，代入地质灾害气象预警模型进行预警反演（公式（3）、（4），并与目前我国大部分地区采用的有效降雨模型（公式（5）；温铭生等，2019）的预警结果进行对比验证。Wj=QjPj（3）Tj=|红色预警Wj Wr橙色预警Wo Wj Wr黄色预警Wy Wj Wo不预警Wj Wy（4）Pj=6i=00.84iRi（5）公式中，j 为网格预警单元序号；Qj为第 j 个网格预警单元的地质环境条件量化指标；Wj为第 j 个网格预警单元的预警指数；Pj为第 j 个网格预警单元的诱发因子指数，文中采用有效降雨模型计算并进行归一化后进行表达，mm；Wr为预警指数的红色预警临界值；Wo为预警指数的橙色预警临界值；Wy为预警指数的黄色预警临界值；Tj为第 j 个网格预警单元的预警等级；i 为计算日前期降雨时间，天；Ri为计算日前第 i 日的日降雨量，mm。预警结果显示，此次降雨过程引发地质灾害4 起（表 6），在相同的预警模型及预警指标的前提下，采用文章提出的有效降雨模型（公式（2），以下简称“新模型”）与采用目前应用较多的有效降雨模型（公式（5），以下简称“原模型”）计算有效降雨均命中了全部地质灾害（图 12，图 13）。但采用新模型较采用原模型的预警时间短、预警等级低、预警区面积更小（表 7）。尤其在持续降雨的中后期，有效降雨采用新模型的预警区面积显著减小，预警区面0204060801001201401601802006 月 22 日6 月 23 日6 月 24 日6 月 25 日6 月 26 日6 月 27 日6 月 28 日6 月 29 日6 月 30 日7 月 1 日7 月 2 日7 月 3 日7 月 4 日7 月 5 日降雨量/mm日期图112021 年 6 月 28 日7 月 3 日福建南平地区降雨分布图Fig.11Rainfall distribution in Nanping,Fujian from June 28 toJuly3,2021y=0.0435x+1938.1R2=0.8056 0204060806.577.588.55 月 8 日5 月 9 日5 月 10 日5 月 11 日5 月 12 日5 月 13 日降雨量/mm水位距地表深度/m日期降雨实测水位趋势性水位图9地下水位消落曲线Fig.9Fluctuationcurveofthegroundwaterlevel第1期陈春利，等：福建省地质灾害气象预警有效降雨模型研究105积平均减小 55%（图 14）。6 月 28 日，采用新模型和原模型评估有效降雨的预警结果均为局部橙色预警、大范围黄色预警，两个模型的预警区面积较接近，实际发生的地质灾害均位于黄色预警区内。7 月 1 日，采用新模型评估有效降雨的预警结果为局部橙色预警、大范围黄色预警，而原模型的预警结果为局部红色预警、大范围橙色和黄色预警，实际发生的地质灾害均位于黄色预警区内，采用新模浦城县武夷山市政和县建瓯市延平区顺昌县邵武市光泽县松溪县建阳区南平市图 例县级行政中心市级行政中心黄色预警区橙色预警区地质灾害点0 10 20kmNb浦城县武夷山市政和县建瓯市延平区顺昌县邵武市光泽县松溪县建阳区南平市图 例县级行政中心市级行政中心黄色预警区橙色预警区地质灾害点a0 10 20kmNa原模型预警结果b新模型预警结果图126 月 28 日预警结果Fig.12WarningresultsofJune28(a)Resultsoftheoriginalmodel;(b)Resultsofthenewmodel浦城县武夷山市政和县建瓯市延平区顺昌县邵武市光泽县松溪县建阳区南平市图 例县级行政中心市级行政中心黄色预警区橙色预警区地质灾害点b0 10 20kmN浦城县武夷山市政和县建瓯市延平区顺昌县邵武市光泽县松溪县建阳区南平市图 例县级行政中心市级行政中心黄色预警区橙色预警区地质灾害点a0 10 20kmNa原模型预警结果；b新模型预警结果图137 月 1 日预警结果Fig.13WarningresultsofJuly1(a)Resultsoftheoriginalmodel;(b)Resultsofthenewmodel表6发灾情况一览表Table6Listofgeologicaldisasters发生日期发生时间发灾地点灾害类型灾害规模/m3造成人员伤亡/人直接经济损失/万元主要诱因2021-6-2813时南平市建阳区书坊乡松坑村南木溪10号泥石流300002特大暴雨2021-6-2812时南平市建阳区建阳区书坊乡饶坝村滑坡8000310特大暴雨2021-7-19时南平市浦城县九牧镇中垄村崩塌5000.2降雨2021-7-110时南平市邵武市水北镇三都村下洪源组滑坡60000.1暴雨106地质力学学报https:/2023型的预警区面积为原模型的 56%。结合地质灾害气象风险预警等级划分标准（中国地质灾害防治工程行业协会，2018），采用新模型的预警结果更合理。同时，采用新模型评估有效降雨在 7 月 3 日至 7 月5 日未发布预警，缩短了持续预警时长，结合实际地质灾害发生情况分析该时间段不预警更合理，避免了预警空报。可见，在保证预警命中率的前提下，采用新模型评估有效降雨进行预警计算显著减少空报面积，提升了预警精准性。5结论（1）福建地区降雨诱发地质灾害受当日降雨和前期降雨诱发作用显著，通过筛选诱发地质灾害的典型降雨事件，通过相关性分析确定了福建地区的有效降雨系数为 0.79，降雨影响期为 3 日内，提出了福建地区有效降雨模型；（2）通过现场降雨入渗监测验证，有效降雨模型可通过降雨量因子模拟降雨条件下斜坡浅表层地下水位的雨后上涨消落的趋势，验证了福建地区有效降雨模型的合理性；（3）将福建地区有效降雨模型应用于地质灾害气象预警，结果显示，采用文中提出的有效降雨模型进行预警分析可以在保证地质灾害命中率的前提下，减少预警区面积、降低预警等级、缩短预警持续时间，提高地质灾害气象预警的精准性。ReferencesALEOTTIP，2004.Awarningsystemforrainfall-inducedshallowfailuresJ.EngineeringGeology，73（3-4）：247-265.BAUMRL,GODTJW，2010.Earlywarningofrainfall-inducedshallowlandslidesanddebrisflowsintheUSAJ.Landslides，7（3）：259-272.CHENCL,HEK,WENMS,etal.，2016.Trendpredictionresearchofgeolo-gical hazard in the Yangtze economic zone based on gray systemtheoryJ.JournalofGeomechanics，22（3）：706-713.（inChinesewithEnglishabstract）CHENJW,1990.Studyondebr