关于突发地质灾害中应急数值模拟技术应用浅析论文

Word第第页关于突发地质灾害中应急数值模拟技术应用浅析论文突发地质灾难应急是一种特别规防灾减灾行动，具有跨阶段、高要求、大集成、快反应和求实效的特征。突发地质灾难大事应急响应具有社会聚焦性和时间紧迫性，处置决策离不开科学高效的技术支撑，数值模拟技术具有低本钱、高效率、且多工况模拟等优势，在公共平安、气象、水资源与环境等应急大事领域得到广泛的应用，其决策帮助功能突出。

然而，数值模拟技术在地质灾难应急领域现场应用却并不多见，尚未形成一种范式的程序化方法，对灾难应急支撑效果非常有限，其主要缘由有6点:

①在应急状况下，数值模拟的应用受到技术条件限制；

②在数值模拟理论讨论与应用实践中，长期以来对时间效率考虑较少，导致常规模拟耗时较长，这与抢险救灾时效性强的特点有差距；

③缺乏特地应急情景下的数值模拟技术适应性讨论；

④突发地质灾难应急工作本身具有技术集成性、方法协调性和追求实效性的特点；

⑤灾难体自身的特异性确定了好像很难找到普遍适用的数值方法和模型；

⑥数值分析人员往往在软件模拟中带有较强主观性，降低了结果的可信度。

但是，数值模拟技术在常规地质灾难讨论领域却相当广泛，结合应急实践阅历来看，两者在功能和思路基本相像，大体上包括稳定性分析、变形机制分析和还原或推演灾难过程等3个方面，区分在于常规数值模拟缺乏应急情境适应性，例如，建模、选参及模拟结果表达等方面的特别性。

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

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

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

2.1技术条件对比

文献分析说明数值模拟在常规和应急状态讨论内容具有全都性，即灾难体稳定性分析、变形机理分析和破坏过程推演等3方面。现场应急处置过程中，应急调查快速查明灾难体周边工程地质条件，是现场应急模拟几何模型和物理力学参数等基本信息的主要数据来源，是建立初步数值模型的基础。此外，现场动态应急监测实时把握灾难进展趋势，为验证数值模型的正确性和修正参数取值供应拟合数据，以便灾难体关键部位的快速确立。应急治理阶段，数值模拟便利考虑支护结构与灾难体的耦合作用，动态反映处置措施的有效性，满意灾难应急求实效的客观要求。然而常规数值模拟应用非直面应急情景，必定与应急数值模拟技术条件存在差异。

2.2关键技术问题

应急决策需要快速研判滑坡稳定性、推演可能的滑动过程，不要求模拟结果非常精确，但必需是牢靠的，这是平安法则所确定的。分析突发地质灾难应急响应技术流程，说明数值模拟技术在实际应急处置过程可以发挥帮助应急决策作用，应急响应启动之后便可介入。表3从应急建模、数值方法选择和参数选取及结果表达等3个方面，提出数值模拟技术在应急过程中应用的关键技术方法。

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

3典型滑坡应急案例实证讨论

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

3.1应急建模

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

3.2结果解读与实证分析

3.2.1滑坡启动缘由

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

3.2.2滑动过程推演

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

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

4结语

突发地质灾难应急现场对数值模拟技术提出了对时间、模拟结果和牢靠度的更高要求。实际灾难应急是多学科、多领域、多部门和多层次的综合性防灾减灾活动。灾难应急现场对数值模拟技术的应用提出了实际需求，虽然常规数值模拟在灾难领域的讨论存在非基于时间和应急模式模拟