突发地质灾害中应急数值模拟技术应用浅析*

李俊峰　张小趁　刘红岩③　倪天翔　褚宏亮　陈红旗

(①中国地质大学(北京)工程技术学院　北京　100083)

(②华北科技学院　北京　101601)

(③西藏大学工学院　拉萨　850000)

(④中国地质环境监测院(国土资源部地质灾害应急技术指导中心)　北京　100081)



突发地质灾害中应急数值模拟技术应用浅析*

李俊峰①张小趁②刘红岩①③倪天翔④褚宏亮④陈红旗④

(①中国地质大学(北京)工程技术学院北京100083)

(②华北科技学院北京101601)

(③西藏大学工学院拉萨850000)

(④中国地质环境监测院(国土资源部地质灾害应急技术指导中心)北京100081)

以滑坡为例，通过综合分析大量前人研究成果，阐述数值模拟在常规状态地质灾害稳定性分析、滑坡机制分析和滑动过程推演中的应用状况，总结数值模拟技术应用于地质灾害应急的基础，比较常规状态和应急情境下数值模拟技术条件的区别，明确应急数值模拟过程特点，提出应急数值模拟技术关键问题。最后，以峨眉山九里镇王山—抓口寺滑坡灾情为例，结合非连续变形分析方法(DDA)大变形和真实时间模拟的优势，进一步研究了数值模拟在灾情应急过程中启动原因和滑动过程反演的实际应用。结果表明，滑体滑动过程中最大水平位移约87m，峰值速度7m·s-1，滑距约38m，前缘堆积厚度53m，滑坡滑动过程持续时间70s左右，与应急调查结果相符。模拟结果有助于理解滑坡发生机制，有助于确立滑坡应急防治重点，辅助应急决策效果明显。

地质灾害数值模拟应急响应滑坡DDA

0　引　言

突发地质灾害应急是一种非常规防灾减灾行动，具有跨阶段、高要求、大集成、快反应和求实效的特征(刘传正，2010)。突发地质灾害事件应急响应具有社会聚焦性和时间紧迫性，处置决策离不开科学高效的技术支撑，数值模拟技术具有低成本、高效率、且多工况模拟等优势，在公共安全、气象、水资源与环境等应急事件领域得到广泛的应用，其决策辅助功能突出(李俊峰等，2016)。然而，数值模拟技术在地质灾害应急领域现场应用却并不多见，尚未形成一种范式的程序化方法，对灾害应急支撑效果十分有限，其主要原因有6点：①在应急情况下，数值模拟的应用受到技术条件限制；②在数值模拟理论研究与应用实践中，长期以来对时间效率考虑较少，导致常规模拟耗时较长，这与抢险救灾时效性强的特点有差距；③缺乏专门应急情景下的数值模拟技术适应性研究；④突发地质灾害应急工作本身具有技术集成性、方法协调性和追求实效性的特点；⑤灾害体自身的特异性决定了似乎很难找到普遍适用的数值方法和模型；⑥数值分析人员往往在软件模拟中带有较强主观性，降低了结果的可信度。但是，数值模拟技术在常规地质灾害研究领域却相当广泛，结合应急实践经验来看，两者在功能和思路基本相似，大体上包括稳定性分析、变形机制分析和还原或推演灾害过程等3个方面，区别在于常规数值模拟缺乏应急情境适应性，例如，建模、选参及模拟结果表达等方面的特殊性。

刘传正等(2010)在重大地质灾害应急响应技术支撑体系中提出数值模拟技术的需求，肯定其辅助应急决策的作用。张小趁等(2015a)以时间为主线提出突发地质灾害应急过程模式，明确了应急情境下灾害体稳定性判断、情景演变及后果分析、处置方案优化、划定灾害区域和应急措施有效性评定等应急重难点问题。张小趁等(2015b)提出的“情景-应对”模式是符合现场应急实际的，结合数值模拟技术自身优势，满足现场情景多变的现实要求，其在解决上述应急重难点问题上辅助决策功能突出。例如，灾情应急侧重成因判断和减灾措施，可借助数值模拟反演，理解灾害形成原因、危害机制和过程；险情则强调动态监测预警和灾害体应急处置措施的选取，借助数值模拟技术可研判危险性、预测危害范围，并为应急治理现场概念设计提供依据。应急避险是防灾减灾成效和降低损失的关键一环，而实际成功避险案例不多，少数成功案例也具有偶然性和随机性，且多源于群测群防等传统技术手段，设想数值模拟能够实现快速的情景推演和危害范围圈定，将具有直观的实效性，符合科学防灾减灾的需求。

本文通过阐述数值模拟技术在滑坡常规应用研究的基础之上，比较分析常规和应急两种模拟工况下的技术条件异同点，提出应急情景下滑坡数值模拟技术应用的关键技术条件和问题。最后，以2015年6月14日峨眉山九里镇王山—抓口寺滑坡为例，根据现场应急需求，利用DDA能较符合实际模拟大变形和真实时间的优势，探讨了数值模拟技术在应急响应行动的应用实践和效果。结果表明，数值模拟技术对应急决策辅助作用明显。鉴于该项研究尚粗浅，仅供讨论商榷。

1　数值模拟技术在滑坡应急应用中的基础

在地质灾害理论研究和防治工程实践过程中，已经积累了大量数值模拟技术研究成果。尽管尚缺乏应急情境下数值模拟技术应用的专门研究，但是通过常规数值模拟的应用，尤其是单一问题的实践，可以探求数值模拟技术应急条件下应用的技术基础，相关文献分析结果及其对应急数值模拟技术应用的参考意义(表1)。

表1　数值模拟技术在常规地质灾害中应用对应急的意义

Table 1　Numerical simulation technology applied in conventional geo-hazard of the significance for emergency

常规应用文献研究内容应急适宜性临滑研判张雪东等(2005)FLAC3D呷爬库岸滑坡天然和蓄水工况分析数值分析的多工况分析便宜性蔡跃等(2008)UDEC模拟边坡反倾层弯曲倾倒破坏型式,分析边坡稳定性影响因素参数选取和取值可参考成永刚等(2008)Geoslope模拟全断面无防护条件下坡体稳定性,强调关键点的监测充分结合应急监测信息反馈和数值模拟内力分析优势,为应急治理的设计提供依据毕小勇等(2015)Midas-GTS分析边坡稳定性,云图反映滑动面位置,动画反映滑动面形成、贯通和发展过程满足应急处置现场提出的对滑体稳定及其发展过程直观化和可视化的要求机制分析殷跃平等(2008)FLAC3D分析多应急治理工况下滑体变形破坏特征,确立应急堆载和变形控制后施用预应力锚索的应急治理措施辅助应急治理,强化措施的有效性和实用性曹琰波等(2011)UDEC分析唐家山滑坡变形累计和破坏过程,失稳前塑性本构模型和M-C屈服准则,失稳后滑体设为刚性模型岩土体本构模型的选择直接关系到应急数值模拟的效率曹琰波等(2011)UDEC分析唐家山滑坡在自重和地震条件下,滑坡启动破坏发生机制与诱发因素的对应关系数值模拟能动态反映诱发因素对灾害体变形的影响,确定影响灾害体变形失稳的关键诱发因素葛云峰等(2014)3DEC分析鸡尾山变形机制,基于正交试验原理,利用参数敏感性手段分析滑动面节理参数对变形特征的影响的显著性常规情境下参数的敏感性分析有助于应急条件下迅速掌握灾害体关键参数和取值吴韩等(2015)区别以往统一强度参数赋值,通过编制程序实现模型的差异化赋值,结果表明精细化赋值能更好反映研究对象的力学响应在条件时间允许的情况下,可适当考虑精细化赋值过程推演Chenetal.(2003)LFEM方法建立准三维模型反演香港大屿山滑坡滑动过程,结合GIS确立灾害区域前景GIS与数值模拟更有效的结合,如确定灾害区域、堆积体体积等邬爱清等(2006)DDA反演千江坪滑坡启动条件和滑动过程,分析表明启动原因为降雨及库水位上升,得到滑坡最大位移速度、最大滑动错距和持续时间。说明应急数值模拟能辅助成因分析和成灾过程反演或预测Pirullietal.(2007)RASH3D反演方量基本一直而滑动路径区域地貌差别很大的两个岩质滑坡,指出理想模拟结果要求堆积区域和堆积体深度方向土石分布与实际一致堆积体区范围和堆积体深度方向岩土体分布直接关系到应急搜救和应急治理有效性张龙等(2012)PFC3D模拟鸡尾山滑坡滑动全过程,指出滑坡堆积区分布与滑面摩擦因素和滑体强度的影响,滑距和区域地形的关系滑体滑动过程及堆积方面,国内外学者进行大量参数敏感性分析研究。有助于应急条件下对滑坡滑动过程反演和预测的参数选择和取值;为应急数值模拟提供参数类型参考;有利于确定研究重点,也存在参数过多和复杂的不利因素,具体需根据现场状况综合确定影响因素Cascinietal.(2014)准3D耦合SPH模型反演意大利南部降雨诱发滑坡,指出滑床夹带率和侵蚀区范围影响滑坡流通路径、滑坡速度和堆积物高度Katzetal.(2015)二维离散元RICEBAL研究单体滑坡尺寸和堆积体形状控制因素,指出滑坡尺寸受材料峰值强度和斜坡初始角控制,堆积体形状受材料残余摩擦角和斜坡初始角控制Zhangetal.(2014)PFEM模拟中国南部一滑坡,结果显示滑动面的几何形态对于滑体的运动有重要影响,滑面与滑体的摩擦角对堆积体形状有重要影响,而滑体材料密度不论是对于滑动过程还是堆积体轮廓均没有影响杜娟等(2015)利用有限体积法,考虑滑体下表面侵蚀作用和摩擦阻力的变化,建立滑坡-碎屑流运动过程三维计算模型,对比未考虑滑体下表面侵蚀作用的计算结果,发现侵蚀作用对于滑体运动物质总体积和滑距有重要影响。段钊等(2014)2D-Block模拟人工切坡公路滑坡,确定滑坡破坏模式和阶段,预测滑坡失稳威胁距离和堆积体体积灾害区域的划定直接关系应急避险方案的选择

总体上看，滑坡险情灾害体稳定性分析和灾情二次滑动可能性判定是其应急响应启动后面临的首要问题。目前主要通过两种途径进行灾害体稳定性分析：极限平衡法和数值极限分析法。目前，极限平衡法在工程设计和可行性分析应用较广，也是规范推荐的一种方法；数值极限分析法自提出起已获得广泛认可和应用，具有一定可靠性，避免了传统方法需假定滑动面的不足，便于内力分析和考虑结构共同作用对稳定性的影响(郑颖人，2012)。其次，实际应急防治过程中，要求遵循先稳后治的原则，如何选择关键时机、关键部位(张小趁等，2015c)或工程措施，数值模拟能直观反映灾害过程，方便获取各阶段内力变化图，可辅助理解灾害机理和确立防治重点(李俊峰等，2016)。最后，现场应急经验告诉我们，应急最为棘手的不是稳定性验算，而是如何反演或预测滑距、滑速和影响范围(杨海清等，2015)。这方面，应用较多、效果较好的有流固耦合和离散元两类方法。

2　滑坡应急数值模拟关键技术

2.1技术条件对比

文献分析表明数值模拟在常规和应急状态研究内容具有一致性，即灾害体稳定性分析、变形机理分析和破坏过程推演等3方面。现场应急处置过程中，应急调查快速查明灾害体周边工程地质条件，是现场应急模拟几何模型和物理力学参数等基本信息的主要数据来源，

表2　数值模拟应急与常规技术条件比较

Table 2　Comparison of emergent and conventional condition in numerical simulation

对比内容应急数值模拟常规数值模拟研究目的突发地质灾害应急处置救援,防灾减灾科学研究研究对象突发地质灾害,灾情和险情两种情况常规地质灾害工况分析具体可能的几种固定工况通常会考虑多种可能工况,对比分析监测数据需现场进行应急监测,由于突发地质灾害的随机性和不确定性,基本没有既有监测数据研究对象一般已进行或可进行较精确的监测,可以根据需要安排监测内容物理力学参数获取通常无试验数据,只能通过现场快速量测、已有附近工程资料、经验值和专家意见确定参数取值可通过室内外试验,如钻探、平硐、现场原位测试等手段获取较精确取值几何模型构建现场即时获取,如照片、3S技术(GPS、RS、GIS)、目测画示意图等手段,且以2D为主可以通过后期加工获取较为准确的几何模型,且3D模拟更能反映实际边界范围圈定尚未有相关方面研究,然其范围圈定直接影响几何模型大小及工作量目前见于文献的主要有两种方式,其一,张倬元等(2009)对应力影响域的描述确定边界范围;其二,郑颖人等(2006)规定有限元模拟中坡脚边界为坡高的1.5倍,坡顶边界为坡高的2.5倍,且上下边界总高不低于2倍坡高,计算精度较好本构模型必须选用已有常规和较为成熟的本构模型,在尽可能准确的描述岩土体性质的基础上,要求参数少且易获取无特殊要求,根据分析目的,选用逼近真实状况的本构模型,一般不考虑参数获取难易程度,若没有合适的本构模型,则推导和修正现有模型,但通常具有特异性时间要求在保证一定的可靠度的情况下,越快越好没有时间要求,以研究目的为导向

是建立初步数值模型的基础。此外，现场动态应急监测实时掌握灾害发展趋势，为验证数值模型的正确性和修正参数取值提供拟合数据，以便灾害体关键部位的迅速确立。应急治理阶段，数值模拟方便考虑支护结构与灾害体的耦合作用，动态反映处置措施的有效性，满足灾害应急求实效的客观要求。然而常规数值模拟应用非直面应急情景，必然与应急数值模拟技术条件存在差别，表2(李俊峰等，2016)从研究目的、对象、工况等9方面分析比较了两者技术条件的差异。

2.2关键技术问题

应急决策需要快速研判滑坡稳定性、推演可能的滑动过程，不要求模拟结果十分精确，但必须是可靠的，这是安全法则所决定的(李俊峰等，2016)。分析突发地质灾害应急响应技术流程(刘传正等，2010；张小趁等，2015a)，表明数值模拟技术在实际应急处置过程可以发挥辅助应急决策作用，应急响应启动之后便可介入。表3从应急建模、数值方法选择和参数选取及结果表达等3个方面，提出数值模拟技术在应急过程中应用的关键技术方法。

针对灾害应急情景模式，使数值模拟技术应用更加高效，有必要分类建立灾害类型数据库，库内包含常见典型应急案例数值模型建立、方法选择和主要参数类型及取值等基本类，以便应急处置过程中直接调用，在此基础上作出特异性修改；同时为规范化建模，有必要探索半自动和自动化应急建模技术，最终实现参数输入、模型输出的现场应急快速建模，也有利于数值模拟技术更为广泛的应用，当然也要充分降低数值模拟主观性因素的影响，增加模拟的准确性。

表3　数值模拟技术应急应用关键技术方法

Table 3　Key technical method in numerical simulation technology of emergency

技术方法内容内容具体描述备 注应急建模地质模型以现场应急调查为基础,现场(或远程会商)应急专家意见为参考,是数值模型建立的基础和修正的根据专家经验与现场实际相结合几何模型其概化基于现场或历史地勘资料的二维主剖面,边界区域圈定(表2)边界条件以位移边界条件为主,适当考虑构造应力边界,特殊工况如地震则需设置底部黏滞边界数值模型确立合理边界条件,根据现场经验、历史资料和参数反演综合确立物理力学参数,本构模型要求简单而较为准确保证可靠性前提,尽可能简化模型表达数值方法选择稳定性分析及变形机制分析连续介质数值方法和非连续介质方法均可,连续方法为佳,如FLAC3D、Geo-studio和ANSYS等有限元适合处理非线性问题不适用模拟非连续大变形问题,内力分析优势明显;FLAC方法在处理非线性大变形问题具有优势,且在多工况分析和支护型式效果分析应用广泛滑动过程推演非连续介质数值方法为主,如DDA、UDEC、PFC等离散元法适合处理非均质、非连续和大变形问题,滑动过程推演模拟优势明显,重视局部失稳导致计算过程失败和耗时问题参数选取与结果表达参数类型岩土体参数与数值方法和本构模型直接相关,但存在基本物理参数,应急现场经验取值和工程类比是参数主要来源结构面包括节理面、岩层面和滑动面结构面参数结果表达 据分析目的不同,内力机制分析包括应力云图、位移云图和塑性区云图等,运动过程包括位移时间关系图、速度时间关系图和堆积范围图等地质灾害应急处置是多学科、多领域、多部门和多层次的综合性防灾减灾活动,要重视将模拟结果以图片和动画的形式输出,使非专业人士也能及时了解和掌握险情或灾情动态

3　典型滑坡应急案例实证研究

2015年6月14日，位于峨眉山市九里镇兴阳村九沙河右岸发生大型顺层岩质滑坡，滑坡体总体积约500×104m3，27户共52人安全受到威胁，滑体前缘冲进九沙河，形成厚度达50m左右的堰塞坝体，汛期对下游存在重大安全隐患(图1)。应急现场迫切需要对滑坡发生机理和滑动过程进行分析，以便迅速确定滑坡成因和应急治理重点。在应急调查观测基础上，经综合考虑，结合非连续大变形分析(石根华，1997)真实时间模拟和大变形问题处理的优势，对滑坡滑动成因和滑动过程进行推演，辅助应急决策。

图1　滑坡地质剖面图Fig.1　Geological profile of landslide

3.1应急建模

峨眉山九里滑坡灾情发生后，迫切需要对灾害发生过程进行反演，确立防治重点。通过滑坡历史资料收集和现场应急调查，确定潜在滑动面和基岩出露情况，滑体被两组节理切割破碎，凝灰岩滑带泥化现象严重，充分利用三维激光扫描技术非接触式测量的技术优势(褚宏亮等，2015)，获得滑坡表面云图，结合历史卫星遥感图，推测原始地面(图1)，将滑坡划分为滑体和滑床两个部分，滑床为一个整体，初步建立滑坡二维几何模型(图2)。滑坡物理力学参数取值综合考虑经验取值、反演分析和既有资料。本工程以历史资料为主，辅以专家意见，利用数值分析参数反演的优势综合确定符合现场实际的物理力学参数。根据计算模型简化结果，将滑坡参数分为岩体和节理参数两类(表4)。模型的左右边界和底部边界设为固定边界，上部自由边界以模拟滑体在自重作用下的位移。根据现场应急调查，滑坡出露基岩位置为滑动面，故将基岩以下划分为一个块体；单元尺寸趋于均一，以便充分反映降雨弱化滑体物理力学参数对滑体稳定状况的影响，滑体部分块体系统包含块体单元306个，布置前中后3个监测点(图3)。

图2　滑坡几何模型示意图Fig.2　Geometry model of landslide

图3　滑坡DDA数值计算模型Fig.3　DDA numerical model of landslide

表4　滑坡材料力学参数取值

Table 4　Mechnical parameters of landslide

材料参数密度/kg·m-3弹性模量/GPa泊松比黏聚力/kPa内摩擦角/(°)抗拉强度/MPa岩体250010.25———节理———0～300～300

图4　天然状态滑坡关键块体稳定系数和滑距与时间关系曲线　Fig.4　Natural state of landslide stability coefficient of a key block and slip curve relationship with timea.天然状态关键块体稳定性系数-时间曲线；b.天然状态监测块体滑距-时间曲线

图5　滑坡滑动全过程位移、速度与时间关系曲线Fig.5　Relation curve of displacement，velocity with time during whole process of landslide slidinga.滑坡过程水平位移-时间曲线；b.滑坡过程滑距-时间曲线；c.滑坡过程速度-时间曲线

3.2结果解读与实证分析

3.2.1滑坡启动原因

通过现场应急调查、参数反演和历史数据收集，滑体天然状态下黏聚力约30kPa，内摩擦角约30°。天然状态模拟结果显示，由于滑体前缘存在髙10余米的临空面，滑体前缘部分块体发生滑动，11#监测块体显示滑距最大为0.65m，滑体缓慢蠕动后，关键块体稳定系数保持在1.8左右，说明滑坡整体稳定(图4)。应急调查表明，现场水文地质条件简单，主要为第四系松散地层孔隙潜水和基岩裂隙水，滑坡处于地下水位以上，且灾情发生前该区域多日连续降雨。综合模拟结果及实际情况初步判定造成灾情的主要诱因是降雨弱化节理岩体物理力学参数，临空面为滑坡发生提供有利条件。

3.2.2滑动过程推演

根据前述对滑坡发生诱因的初步判断，拟合现场滑体变形和关键标志物的运动情况，通过数值参数反演确定实际滑体物理力学参数估计值。根据DDA模拟结果，滑坡滑动分为启动、快速滑动和堆积3个阶段，启动阶段滑体前缘块体首先滑出，滑体后缘出现拉裂面；快速滑动阶段，滑体表层松动，部分块体飞起，滑体滑动过程中最大水平位移约87m，速度峰值最高达7m·s-1；堆积停止阶段，随着滑坡势能的摩擦消耗，以及滑体前缘出现反压，滑体中下部趋于稳定，表层碎石土局部发生缓慢滑动，滑距约为38m，前缘堆积体厚度53m，持续时间70s左右(图5，图6)。

图6　滑坡滑动过程变形图Fig.6　Deformation pattern of landslide sliding processa.0s滑体变形情况；b.3.5s滑体变形情况；c.21s滑体变形情况；d.70s滑体变形情况

综上所述，王山—抓口寺滑坡破坏模式为以后部推移为主、前部牵引为辅的“强降雨-入渗-岩体软化-推拉”破坏模式。模拟结果很好的为确立应急治理重点和理解滑坡发生机制提供了参考，辅助应急效果明显。

4　结　语

突发地质灾害应急现场对数值模拟技术提出了对时间、模拟结果和可靠度的更高要求。实际灾害应急是多学科、多领域、多部门和多层次的综合性防灾减灾活动。灾害应急现场对数值模拟技术的应用提出了实际需求，虽然常规数值模拟在灾害领域的研究存在非基于时间和应急模式模拟等诸多应急工况适应性不足问题，但是也为数值模拟技术应急应用提供了大量素材。应急过程中，应急调查为数值模拟提供了基础数据，应急监测作为验证建模准确性和模拟结果的合理性的重要手段。通过峨眉山九里镇王山—抓口寺滑坡灾情的数值模拟技术应用分析，直观的反映了数值模拟在灾害应急的极大优势，进一步说明了数值模拟技术对应急决策辅助效果明显。

Bi X Y，Yan T J，Lu J.2015.Application of Midas-GTS(SRM)to 2D stability analysis of slope[J].Journal of Natural Disasters，24(1)：170～176．

Cai Y，Mitani Y，Esaki T.2008.Numerical analysis of stability for an anti-dip stratified rock slope[J].Chinese Journal of Rock Mechanics and Engineering，27(12)：2517～2522．

Cao Y B，Dai F C，Xu C，et al.2011.Discrete element simulation of deformation and movement mechanism for Tangjiashan landslide[J].Chinese Journal of Rock Mechanics and Engineering，30(S1)：2878～2887．

Cascini L，Cuomo S，Pastor M，et al.2014.SPH run-out modelling of channelised landslides of the flow type[J].Geomorphology，214(2)：502～513．

Chen H，Lee C F.2003.A dynamic model for rainfall-induced landslides on natural slopes[J].Geomorphology，51(4)：269～288．

Cheng Y G.2008.Numerical simulation and monitoring analysis of bedding landslide[J].Chinese Journal of Rock Mechanics and Engineering，27(S2)：3746～3752.

Chu H L，Yin Y P，Cao F，et al.2015.Research on deformation monitoring of large collapses and landslides based on 3D laser scanning technology[J].Hydrogeology & Engineering Geology，42(3)：128～134．

Du J，Yin K L，Wang J J.2015.Simulation of three-dimensional movement of landslide-debris flow based on finite volume method[J].Chinese Journal of Rock Mechanics and Engineering，34(3)：480～488．

Duan Z，Tang H，Dang Q，et al.2014.Discrete lelement simulation of landslide induced by slope excavation[J].Journal of Chang’an University(Natural Science Edition)，34(5)：49～55．

Ge Y F，Tang H M，Xiong C R，et al.2014.Effect of sliding plane mechanical parameters on landslide stability-a case study of Jiweishan rockslide in Wulong，Chongqing[J].Chinese Journal of Rock Mechanics and Engineering，33(S2)：3873～3884．

Katz O，Morgan J K，Aharonov E，et al.2015.Controls on the size and geometry of landslides：Insights from discrete element numerical simulations[J].Geomorphology，220(3)：104～113．

Li J F，Zhang X C，Liu H Y，et al.2016.Numerical simulation technology analysis of geological disaster emergency[J].The Chinese Journal of Geological Hazard and Control，27(2)：154～161．

Liu C Z，Chen H Q，Han B，et al.2010.Technical support system of emergency response for serious geo-hazards[J].Geological Bulletin of China，29(1)：147～156．

Pirulli M，Mangeney A.2007.Results of Back-Analysis of the Propagation of Rock Avalanches as a Function of the Assumed Rheology[J].Rock Mechanics & Rock Engineering，41(1)：59～84．

Shi G H.1997.Numerical Manifold Method and Discontinuous Deformation Analysis[M].Translated by Pei J M.Beijing：Tsinghua University Press．

Wen J J，Wang L.2015.Numerical study on failure mechanism of a loess landslide[J].Subgrade Engineering，179(2)：184～186．

Wu A Q，Ding X L，Li H Z，et al.2006.Numerical simulation of startup and whole failure process of Qianjiangping landslide using discontinuous deformation analysis method[J].Chinese Journal of Rock Mechanics and Engineering，25(7)：1297～1303．

Wu H, Jian W X, Xiao D, et al.2015.Numerical modeling method for slope stability under variable distribution of strength parameters [J].Journal of Engineering Geology, 23(1)：135～139.

Xie J J，Li Y Q.2011.Rapid evaluation of dam-break hazard based on numerical simulation of flood-routing[J].Journal of Catastrophology，26(2)：31～34．

Yang H Q，Lan Y F，Zeng Y Y，et al.2015.A three-dimensional model with spring-deformable-blocks for post-failure motion of landslides[J].Chinese Journal of Rock Mechanics and Engineering，34(3)：528～536．

Yin Y P，Li T Q，Tang J.2008.Landslide reactivation and emergency stabilization on Danba county town in Sichuan province[J].Chinese Journal of Rock Mechanics and Engineering，27(5)：971～978．

Zhang L，Tang H M，Xiong C R，et al.2012.Movement process simulation of high-speed long-distance Jiweishan landslide with PFC3D[J].Chinese Journal of Rock Mechanics and Engineering，31(S1)：2601～2611．

Zhang X C，Chen H Q.2015a.Geological disaster emergency technology：process model[J].Journal of Catastrophology，30(4)：149～155．

Zhang X C，Chen H Q.2015b.Study on emergent prevention of sudden geological disaster based on scenario-response model[J].Yangtze River，46(23)：29～33．

Zhang X C，Chen H Q.2015c.Geological disaster emergency technical：emergency monitoring[J].Journal of Engineering Geology，23(S1)：355～361．

Zhang X D，Chen J P，Huang R Q，et al.2005.A study of deformation features of gapa landslide using FLAC3D[J].Rock and Soil Mechanics，26(1)：131～134．

Zhang X，Krabbenhoft K，Sheng D，et al.2014.Numerical simulation of a flow-like landslide using the particle finite element method[J].Computational Mechanics，55(1)：167～177．

Zheng Y R.2012.Development and application of numerical limit analysis for geotechnical materials[J].Chinese Journal of Rock Mechanics and Engineering，31(7)：1279～1316．

Zheng Y R，Zhao S Y，Deng C J，et al.2006.Development of finite element limit analysis method and its applications in geotechnical engineering[J].Engineering Sciences，8(12)：39～61．

毕小勇，闫天俊，鲁杰.2015.Midas-GTS(SRM)在边坡二维稳定性分析中的运用[J].自然灾害学报，24(1)：170～176．

蔡跃，三谷泰浩，江琦哲郎.2008.反倾层状岩体边坡稳定性的数值分析[J].岩石力学与工程学报，27(12)：2517～2522．

曹琰波，戴福初，许冲，等.2011.唐家山滑坡变形运动机制的离散元模拟[J].岩石力学与工程学报，30(增1)：2878～2887．

成永刚.2008.顺层滑坡数值模拟与监测分析[J].岩石力学与工程学报，27(增2)：3746～3752．

褚宏亮，殷跃平，曹峰，等.2015.大型崩滑灾害变形三维激光扫描监测技术研究[J].水文地质工程地质，42(3)：128～134．

杜娟，殷坤龙，王佳佳.2015.基于有限体积法的滑坡-碎屑流三维运动过程模拟分析[J].岩石力学与工程学报，34(3)：480～488．

段钊，唐皓，党琪，等.2014.边坡开挖诱发滑坡离散元模拟[J].长安大学学报(自然科学版)，34(5)：49～55．

葛云峰，唐辉明，熊承仁，等.2014.滑动面力学参数对滑坡稳定性影响研究——以重庆武隆鸡尾山滑坡为例[J].岩石力学与工程学报，33(增2)：3873～3884．

李俊峰，张小趁，刘红岩，等.2016.突发地质灾害应急数值模拟浅析[J].中国地质灾害与防治学报，27(2)：154～161．

刘传正，陈红旗，韩冰，等.2010.重大地质灾害应急响应技术支撑体系研究[J].地质通报，29(1)：147～156．

石根华.数值流形方法与非连续变形分析[M].裴觉民，译.北京：清华大学出版社．

文建军，王雷.2015.某黄土滑坡破坏机理数值研究[J].路基工程，(2)：184～186．

邬爱清，丁秀丽，李会中，等.2006.非连续变形分析方法模拟千将坪滑坡启动与滑坡全过程[J].岩石力学与工程学报，25(7)：1297～1303．

吴韩，简文星，肖迪，等.2015.岩土体力学参数空间分布下的滑坡稳定性数值计算方法[J].工程地质学报，23(1)：135～139．

杨海清，蓝一凡，曾酉源，等.2015.滑体失稳后运动过程三维弹簧变形块分析模型[J].岩石力学与工程学报，34(3)：528～536．

殷跃平，李廷强，唐军.2008.四川省丹巴县城滑坡失稳及应急加固研究[J].岩石力学与工程学报，27(5)：971～978．

张龙，唐辉明，熊承仁，等.2012.鸡尾山高速远程滑坡运动过程PFC～(3D)模拟[J].岩石力学与工程学报，31(增1)：2601～2611．

张小趁，陈红旗.2015a.突发地质灾害应急技术：过程模式[J].灾害学，30(4)：149～155．

张小趁，陈红旗.2015b.基于情景应对模式的突发地质灾害应急防治探讨[J].人民长江，46(23)：29～33．

张小趁，陈红旗.2015c.突发地质灾害应急技术：应急监测模式[J].工程地质学报，23(S1)：355～361．

张雪东，陈剑平，黄润秋，等.2005.用FLAC-3D分析呷爬滑坡的变形特征[J].岩土力学，26(1)：131～134．

张倬元，王士天，王兰生，等.工程地质分析原理[M].北京：地质出版社.

郑颖人.2012.岩土数值极限分析方法的发展与应用[J].岩石力学与工程学报，31(7)：1297～1316．

郑颖人，赵尚毅，邓楚键，等.2006.有限元极限分析法发展及其在岩土工程中的应用[J].中国工程科学，8(12)：39-61.

APPLICATION OF NUMERICAL SIMULATION TECHNOLOGY TO EMERGENT RESCUING OF GEOLOGICAL DISASTER

LI Junfeng①ZHANG Xiaochen②LIU Hongyan①③NI Tianxiang④CHU Hongliang④CHEN Hongqi④

(①School of Engineering & Technology，China University of Geosciences(Beijing)，Beijing100083)

(②North China Institute of Science and Technology，Beijing101601)

(③School of Engineering，Tibet University，Lhasa850000)

(④China Institute for Geo-Environmental Monitoring(Geological Disaster Emergency Technical Guidance Center of MLR)，Beijing100081)

This paper is based on the analysis of landslide and a large number of predecessor research results.It summarizes the works of conventional numerical simulation in the analysis of stability，landslide mechanism analysis and reproducing sliding process.It presents the basic conditions of numerical simulation applicable to geo-hazard emergent rescuing.At the same time，it puts forward the state differences between conventional and emergency numerical analyses.It clearly figures out the characteristics of emergency numerical simulation, and the key problems about the application of numerical simulation to geo-hazard emergent rescuing.It gives an application example of the Wangshan-Zhuakousi landslide on the Mount Emei，Sichuan province.It uses the numerical method of discontinuous deformation analysis(DDA)that has outstanding advantages in the analysis of large deformation in real time.The numerical simulation in the process and reason of the landslide have been taken into a further study.The results show，the sliding process is successfully reproduced and the maximum horizontal displacement of the sliding body reached by DDA mode is 87meters.The peak speed is 7m·s-1.The results show that sliding distance is about 38meters and the leading edge thickness is 53meters.The whole duration of the landslide obtained by DDA method is about 70 seconds.The numerical results are consistent with the results of emergent investigation.It’s concluded that numerical simulation in geo-hazard emergency is good for the comprehension of the process of landslide，finding the key point and auxiliary emergency decision．

Geo-hazard，Numerical simulation，Emergency response，Landslide，DDA

10.13544/j.cnki.jeg.2016.04.011

2016-02-29;

2016-06-15.

国家自然科学基金(41162009)，国土资源部公益性行业科研专项“突发地质灾害应急响应支撑关键技术研究”(201211055)，中国地质大学(北京)基本科研业务费专项资金(2-9-2014-019)资助.

李俊峰(1991-)，男，硕士生，主要从事岩土工程方面的科研工作.Email:yookeelee@163.com

简介：刘红岩(1975-)，男，教授，主要从事岩土工程方面的教学和科研工作.Email:lhyan1204@126.com

P642.22

A