论文部分内容阅读
区域滑坡的预报技术作为一项非工程措施,是滑坡减灾工作的重要措施之一。预报单元的选取是区域滑坡预报的基础,主要有五种:栅格单元(Grid cell)、地域单元(Terrain unit)、特定条件单元(Unique-condition unit)、斜坡单元(Slope unit)、地貌单元(Topological unit)。其中,栅格单元与斜坡单元应用最为广泛。斜坡单元比栅格单元更具物理意义,且预报结果的指导性更为明确。然而,目前斜坡单元的提取以地表水文过程分析方法（以下称传统方法）为主,提取的斜坡单元内部存在坡度突变的问题,无法满足滑坡稳定分析要求的均一性假定。并且,基于力学概念的滑坡预警报模型主要依靠安全系数的大小给出滑坡发生与否的确定性预报结果,但区域范围岩土体力学参数的不确定性却导致目前力学预报模型不具有发布滑坡确定性预报的能力。为解决斜坡单元内部坡度突变问题,本文基于滑坡稳定分析的均一性假定,利用地形学与形态影像学,建立了新的斜坡单元提取方法MIA-HSU(Morphological Image Analysis-Homog eneous Slope Unit)。为解决物理模型不具有发布确定性预报结果能力的问题,利用MIA-HSU方法提取的斜坡单元建立了基于物理框架的区域滑坡概率分析模型。基于建立的区域滑坡概率分析模型进行了三峡库区奉节县不同前期降水等级下的区域滑坡评估与预报。通过C＃和Fortran混合编程实现了奉节县区域滑坡预报系统的构建和系统应用。通过以上研究,取得了以下研究成果: (1)基于滑坡分析的均一性基本假定,提出一种新的斜坡单元提取方法MIA-HSU。MIA-HSU方法将斜坡单元定义为三维空间中连续、均质、闭合的小区域。每个小区域具有均一的坡度和坡向,小区域的边界即为地形特征出现微小起伏变化之处。采用形态影像学技术从原始DEM中提取出小区域边界,采用主元分析法提取每个小区域的拟合平面,根据向量相似度准则合并小区域以形成新的斜坡单元。 (2)云南蒋家沟流域提取结果表明,MIA-HSU方法提取的斜坡单元解决了传统方法提取的斜坡单元存在的坡度突变的问题,所提取的斜坡单元具有更为均一的坡度和坡向。通过地面雷达扫描和野外滑坡台地测量等手段,分别比较传统方法和MIA-HSU方法提取的斜坡单元对浅层与深层滑坡地貌特征的反映能力。蒋家沟流域提取结果表明,传统方法提取的斜坡单元不能反映浅层滑坡的地貌形态。MIA-HSU方法提取的的斜坡单元不仅能够满足坡度均一性需求,而且能够准确反映浅层滑坡的地貌形态。奉节县新铺滑坡与寂静滑坡提取结果表明,传统方法提取的斜坡单元无法反映深层滑坡的地貌特征,MIA-HSU方法提取的斜坡单元能够反映深层滑坡中的滑坡台地和坡度较为均一的区域,地貌学意义更加清晰。 (3)针对MIA-HSU方法提取的斜坡单元,研究并建立了基于蒙特卡洛方法的滑坡稳定性概率分析模型。利用计算机图形学和地形学建立斜坡单元计算剖面提取方法。然后,利用蒙特卡洛思想建立斜坡单元临界滑面搜索方法。采用概率密度函数反映输入力学参数的不确定性,假定土体在一定的边界范围内服从均匀分布。以塑限和液限下的力学参数作为上下边界,采用蒙特卡洛法在边界范围内随机取值n次,计算n次取值中安全系数小于1所占的比重P,利用比重P来建立滑坡发生概率与输入参数不确定性之间的关联,最终建立起基于物理机制的滑坡概率分析模型。 (4)进行了奉节县不同前期降水等级下的区域滑坡评估工作。首先,对研究区进行高密度的野外取样和室内实验以提高下垫面力学参数取值精度。然后将前期降水分为7个等级:0mm、30mm、60mm、90mm、120mm、150mm、170mm。每个降水等级分别进行重现期5年、10年、20年、50年、100年的区域滑坡评估。评估结果表明,当前期降水小于120mm时,滑坡的发生明显受到降雨频率的控制,随着重现期的增加,潜在不稳定坡体沿着境内长江、朱衣河、梅溪河两侧逐渐向境内纵深持续扩展。当前期降水大于120mm时,随着重现期的增加,不稳定斜坡单元数目增加趋势趋于平缓,降雨频率对滑坡的激发作用不再明显。 (5)进行了奉节县区域滑坡预报工作。根据重庆市气象局提供的2017年3月29-4月7日逐日的雷达估测降水数据和4月8日雷达预报降水数据,对2017年4月8日降水诱发的区域滑坡事件进行预报。预报结果表明,在野外考察查明的45个滑坡点中,预报成功的滑坡点为41个,漏报率8.8％。采用C＃和Fortran混合编程,结合多线程并行编程思想构建了奉节县区域滑坡预报系统。利用构建的区域滑坡预报系统对2017年10月3日江北部分地区发生的区域滑坡灾害进行预报,在奉节县国土局提供的21个滑坡灾害点中,共成功预报滑坡点19个,漏报率9.5％。 本文的创新性主要体现在: (1)利用形态学等技术手段建立了一种新的斜坡单元提取方法MIA-HSU。该方法提取的斜坡单元具有均一的坡度和坡向,解决了传统方法提取的斜坡单元因内部存在坡度突变而难以满足滑坡稳定分析均一性基本假定的问题。 (2)建立了斜坡单元计算剖面提取方法和斜坡单元临界滑面搜索方法,采用概率密度函数反应输入力学参数的不确定性,建立了基于斜坡单元的区域滑坡稳定性概率分析模型,避免了使用不确定的岩土力学参数输入而获得确定性预报结果的问题。