ISSN 1003-8035 CN 11-2852/P

    降雨和地下水位下降作用下岩溶塌陷的渐进破坏机制与数值模拟研究

    Progressive Failure Mechanism and Numerical Simulation of Covered Karst Collapse under Rainfall Infiltration and Groundwater Drawdown

    • 摘要: 覆盖型岩溶塌陷的数值模拟普遍存在岩体强度参数固定不变、溶洞形态效应缺乏量化等不足,导致渐进破坏过程刻画失真。针对上述问题,本文旨在揭示降雨入渗及地下水位下降作用下岩溶塌陷的渐进破坏机制,并量化溶洞形态对塌陷敏感性的控制效应。以深圳坪山咸水湖可溶岩分布区为研究对象,基于FLAC数值模拟平台,采用内置FISH语言实现岩体强度参数的动态折减,模拟溶洞围岩拉张破坏的发生与发展过程;对比分析狭长型(2m×8m)与椭圆型(6m×8m)溶洞在真空吸蚀作用下的塌陷响应。研究成果如下:(1)实现了对溶洞围岩拉张破坏单元的动态追踪,有效解决了传统固定参数模拟难以刻画渐进破坏的难题;(2)揭示了降雨入渗条件下岩溶塌陷的五阶段演化规律,通过0.5h、0.75h、1h三个典型时刻的位移场演变分析,发现溶洞结构在入渗后期(1h阶段)呈现突变性失稳特征,最终塌陷发生于后期阶段的渗透作用;(3)量化了溶洞形态效应,狭长型溶洞形成顶板塑性贯通区的临界吸力(80kPa)仅为椭圆型溶洞(160kPa)的50%,表明其对真空吸蚀作用更为敏感;(4)系统阐明了岩溶塌陷的力学演化链:顶板初始塑性破坏→拉张裂隙扩展→破坏区垂向延伸→土体塑性化发展→塑性区全面贯通→竖向位移突变致塌,建立了基于顶板位移速率的塌陷预警判据。本研究成果可为岩溶区工程地质勘察、灾害监测预警系统优化及防治工程设计提供理论支撑与技术参考。

       

      Abstract: Numerical simulations of covered karst collapse commonly suffer from deficiencies such as fixed rock mass strength parameters and insufficient quantification of cave morphology effects, leading to a distorted representation of the progressive failure process. To address these issues, this study aims to reveal the progressive failure mechanism of karst collapse under rainfall infiltration and groundwater drawdown, and to quantify the controlling effect of cave morphology on collapse susceptibility. Using the soluble-rock area of Xianshui Lake in Pingshan, Shenzhen, as a case study, the built-in FISH scripting language in the FLAC numerical simulation platform was employed to dynamically degrade rock mass strength parameters, simulating the progressive failure process of tensile failure in the surrounding rock mass of karst caves. A comparative analysis was conducted on the collapse responses of a narrow elongated cave(2 m × 8 m)and an elliptical cave(6 m × 8 m) under vacuum suction erosion. The results are as follows: (1) Dynamic tracking of tensile failure units in the surrounding rock mass of karst caves was achieved, effectively resolving the challenge of depicting progressive failure in conventional fixed-parameter simulations;(2) The five-stage evolution mechanism of karst collapse under rainfall infiltration was revealed. Displacement field evolution analysis at 0.5 h, 0.75 h, and 1 h demonstrated abrupt instability characteristics during the late infiltration phase (1 h stage), with final collapse occurring through late-stage seepage; (3) The morphological effects of karst caves were quantified. The critical suction force (80 kPa) required to form a through-going plastic zone in the roof of elongated caves was only 50% of that for elliptical caves (160 kPa), indicating greater sensitivity to vacuum suction erosion; (4) The mechanical evolution chain was clarified: initial roof plastic failure → tensile crack propagation → vertical failure zone extension → soil plasticization development → full-section plastic zone penetration → collapse induced by sudden vertical displacement. An early-warning criterion based on roof displacement rate was established. This study provides theoretical support and technical reference for engineering geological surveys, disaster monitoring system optimization, and prevention engineering design in karst areas.

       

    /

    返回文章
    返回