ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

基于LBM-DEM细观数值模拟的水力诱导覆盖型岩溶地面塌陷发育过程分析

陶小虎, 叶明, 龚建师, 王赫生, 胡晓雨

陶小虎,叶明,龚建师,等. 基于LBM-DEM细观数值模拟的水力诱导覆盖型岩溶地面塌陷发育过程分析[J]. 中国地质灾害与防治学报,2024,35(1): 124-131. DOI: 10.16031/j.cnki.issn.1003-8035.202207027
引用本文: 陶小虎,叶明,龚建师,等. 基于LBM-DEM细观数值模拟的水力诱导覆盖型岩溶地面塌陷发育过程分析[J]. 中国地质灾害与防治学报,2024,35(1): 124-131. DOI: 10.16031/j.cnki.issn.1003-8035.202207027
TAO Xiaohu,YE Ming,GONG Jianshi,et al. Analysis of the formation process of the covered karst ground collapse induced by groundwater changes based on the coupled LBM-DEM numerical simulation at micro scale[J]. The Chinese Journal of Geological Hazard and Control,2024,35(1): 124-131. DOI: 10.16031/j.cnki.issn.1003-8035.202207027
Citation: TAO Xiaohu,YE Ming,GONG Jianshi,et al. Analysis of the formation process of the covered karst ground collapse induced by groundwater changes based on the coupled LBM-DEM numerical simulation at micro scale[J]. The Chinese Journal of Geological Hazard and Control,2024,35(1): 124-131. DOI: 10.16031/j.cnki.issn.1003-8035.202207027

基于LBM-DEM细观数值模拟的水力诱导覆盖型岩溶地面塌陷发育过程分析

基金项目: 中国地质调查局地质调查项目(DD20221756;DD20230079)
详细信息
    作者简介:

    陶小虎(1987—),男,水力学及河流动力学专业,博士,工程师,主要从事地下水环境研究。E-mail:taoxiaohu_cgs@126.com

    通讯作者:

    王赫生(1984—),男,资源与环境专业,博士,高级工程师,主要从事水工环地质调查研究。E-mail:411099641@qq.com

  • 中图分类号: P642.26

Analysis of the formation process of the covered karst ground collapse induced by groundwater changes based on the coupled LBM-DEM numerical simulation at micro scale

  • 摘要:

    文章以水力驱动的覆盖型岩溶地面塌陷为背景,基于离散元方法和格子Boltzmann方法,采用Bouzidi插值反弹边界格式和动量交换法,建立一种可以从细观角度模拟覆盖型岩溶塌陷的二维格子Boltzmann方法—离散元方法流固耦合模型。在此基础上对承压水下降引起的覆盖型岩溶塌陷数值模拟开展了探索性研究。模拟结果表明:承压水位下降工况中地下水主要对隔水层岩溶开口处的颗粒产生影响,对土洞周围土体产生向下的作用力;土体颗粒的剥落容易造成土颗粒原位置和上方位置处水压的陡降,从而造成较强的水力坡降,使得地下水对内部颗粒作用力陡增,容易引起上方颗粒在地下水作用力和重力作用下失稳,导致从土体颗粒失稳至土层塌陷逐渐加速。研究成果对进一步从细观尺度进行水力驱动的覆盖型岩溶地面塌陷的发育过程与特征研究具有理论及实际意义。

    Abstract:

    To investigate the hydraulic characteristics and soil behaviors during the formation of covered karst ground collapse induced by the groundwater changes, a 2D fluid-solid coupling model was developed based on discrete element method and lattice Boltzmann method. This model utilizes the linearly interpolated bounce-back scheme of Bouzidi and the momentum exchange method, allowing for the simulation of the formation of covered karst ground collapse from a microscopic perspective. Using the fluid-solid coupling model, an exploratory study was conducted to simulate the formation of covered karst ground collapse induced by a decrease in the hydraulic head of confined aquifers. Simulation results indicate that when the water level of a confined aquifer declines, the groundwater flow mainly affects the particles located above a cave opening and produces a downward force on the surrounding soil. When soil particles spall, the hydraulic heads at the positions of the spalled soil particles drop sharply. This results in a significant increase in the hydraulic gradient, causing the groundwater force on internal particles to sharply increase as well. As a result, the upper particles lose stability due to the combined force of groundwater dragging and gravitational force, which can lead to a gradual acceleration process of collapse. The research results provide valuable insights into the understanding of covered karst ground collapse formation induced by the groundwater changes.

  • 覆盖型岩溶地面塌陷的形成是自然环境变迁的一个组成部分,而人类活动对岩溶发育地区自然平衡状态的干扰增多,大大加快了这一地质灾害发生的进程,由此引发的地面塌陷问题日趋频繁、后果严重,已成为岩溶地区可持续发展的一大障碍[12]。在众多覆盖型岩溶塌陷形成过程中,由地下水位变化引起的岩溶塌陷居多,地下水流提供了最主要的动力。地下水过量开采、水库蓄水等引起的地下水位变化为主要诱发因素[36]

    长期以来,国内外学者对地下水位变化引起的覆盖型岩溶地面塌陷问题一直十分关注并开展了卓有成效的研究。在数学解析方面,基于渗透力学、极限平衡理论、牛顿摩擦定律等方法建立简单力学模型,获得宏观规律,确定岩溶塌陷的机制及其临界条件的研究[79]。在物理实验方面,覆盖型岩溶塌陷模型试验研究从单一致塌因素研究到多因素组合、从简单的二维定性研究到三维定量分析以及从真实土到透明土可视化分析都体现了试验研究的巨大进步[1012]。在数值模拟方面,基于连续介质理论,在岩溶塌陷临界条件、主控因素、塌陷模式等内容的研究上取得了较好的成果[1314]。然而由于受到当前科学技术的局限、土洞发育以及塌陷过程的非连续、非均匀、流固耦合等特性,现有的物理实验和数值模拟方法很难对水力诱导覆盖型岩溶塌陷进行系统的岩溶塌陷发育和塌陷过程的细观研究,对各影响因子之间的定量关系研究也受到了限制[15]

    离散元方法(discrete element method, DEM)是根据离散物质本身的离散特性建立数值模型,与连续介质理论方法相比,DEM可以从细观角度完成大变形分析,更好地描述结构破坏过程[1617]。格子Boltzmann方法(lattice Boltzmann method,LBM)是一种新兴的计算流体学方法,能够从介观尺度上描述流体的运动和模拟多相多尺度问题。与传统流体力学计算方法相比,LBM在程序实现上更为简单,能够实现与大量移动颗粒的完全耦合,很好地模拟多相多尺度耦合问题,尤其是在模拟复杂边界条件的多孔介质水流中,更能体现出其相比于传统流体力学计算方法的优势[18]。因此,为了能够从细观上对覆盖性岩溶塌陷模拟研究,本文将采用DEM模拟土体颗粒运动,选择库仑定律来判别土体颗粒之间的黏性接触是否受到破坏,采用LBM中D2Q9模型[19]模拟土体孔隙间水流运动,应用Bouzidi等[20]反弹边界格式和动量交换法实现LBM与DEM双向耦合,建立覆盖型岩溶塌陷二维流固耦合细观数值模型,并对承压水下降工况下的覆盖性岩溶塌陷过程进行探索研究。

    本节主要介绍LBM-DEM双向流固耦合的方法,主要包括固体对流体的作用、流体对固体的作用、耦合模型时间步长同步。

    为了计算固体颗粒与流体之间的相互作用的准确性,特别是流体对固体颗粒的作用力的计算,本文采用具有二阶精度且更为简单的线性插值的Bouzidi插值反弹边界格式(图1),假设每个时步末节点A$ - \alpha $方向的分布函数是由虚拟流体节点D处的粒子向右迁移,经过点C反弹运动至节点A所得,其中一个时间步长所迁移长度为一个空间步长,即$ \left| {DC} \right| + \left| {AC} \right| = \left| {AB} \right| $。取$ q{\text{ = }}\left| {AC} \right|/\left| {AB} \right| $${f_\alpha }(r,t)$为在$ r $位置离散速度在$ {\alpha } $方向上粒子密度分布函数。当$ q < \dfrac{1}{2} $时,虚拟节点D处于流体节点EA之间,可由流体节点EA处的分布函数插值得到虚拟节点D的分布函数;当$ q \geqslant \dfrac{1}{2} $时,虚拟节点D处于流体节点A与固体节点B之间,可以通过节点A处的碰撞后的分布函数$ f_\alpha ^ + (r,t) $$ f_{ - \alpha }^ + (r,t) $进行线性插值得到点A处下一时步的$ - \alpha $方向的分布函数。

    图  1  Bouzidi插值反弹边界示意图
    Figure  1.  Schematic view of Bouzidi interpolation bounce-back diagram

    动量交换法能够精确有效且便捷的处理流体流动对固体的力学作用的计算[21]。因此本文将采用动量交换法评估二维流体对固体的力和力矩。在评估流体对固体颗粒的作用力和力矩时,固体颗粒中的分布函数不参与计算。

    $$ \left\{ \begin{gathered} {F_{\sigma x}}({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t){\text{ = [}}f_\alpha ^ + ({{\boldsymbol{r}}_{\rm{f}}},t){e_{\alpha x}}{-}{f_{ - \alpha }}({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t){e_{{-}\alpha x}}]{\kern 1pt} \frac{{{\delta}{x^2}}}{{{\delta}t}} \\ {F_{\sigma y}}({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t){\text{ = [}}f_\alpha ^ + ({{\boldsymbol{r}}_{\rm{f}}},t){e_{\alpha y}}{-}{f_{ - \alpha }}({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t){e_{{-}\alpha y}}]{\kern 1pt} \frac{{{\delta}{x^2}}}{{{\delta}t}} \\ \end{gathered} \right. $$ (1)
    $$ {{\boldsymbol{F}}_\sigma }({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t){\text{ = }}{F_{\sigma x}}({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t){\text{ + }}{F_{\sigma y}}({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t) $$ (2)
    $$ {{\boldsymbol{T}}_\sigma }({{\boldsymbol{r}}_{\text{f}}},t + {\delta}t) = {{\boldsymbol{F}}_\sigma }({{\boldsymbol{r}}_{\rm{f}}},t + {\delta}t) {{\boldsymbol{r}}_\sigma } $$ (3)

    式中:$ {{{r}}_{\rm{f}} } $−颗粒边界流体节点位置;

    $ {{{F}}_\sigma } $−流体节点处流体对固体颗粒的作用力;

    $ {{{T}}_\sigma } $−流体节点处流体对固体颗粒的作用的力矩;

    $ {\delta}x $−LBM中格子步长;

    $ {{{e}}_\alpha } $−速度方向;

    $ {\delta}t $−时间步长;

    $ {f_\alpha }({{r}},t) $−在$ {{r}} $位置离散速度在$ {\alpha } $方向上粒子密度 的分布函数;

    $ {{{r}}_\sigma } $−作用力点到颗粒圆心的矢量。

    在耦合模型中,涉及到两套方法的时间步长,作如下考虑:

    当LBM时间步长$ {\delta}t $小于或等于DEM预估时间步长$ \delta t $时,DEM时间步长取LBM时间步长,$\delta t{{ = \delta }}t$

    当LBM时间步长$ {\delta}t $大于或等于DEM预估时间步长$ \Delta t $时,令${{n}} = {\text{ceil}}\left(\dfrac{{\delta t}}{{\delta t}}\right)$,将DEM时间步长更新为$\Delta t{\text{ = }} \dfrac{{\delta t}}{{{n}}}$

    实际覆盖层土体中,流体能够穿过由多个颗粒形成的多孔介质,即使相互紧密接触的颗粒之间也存在渗透路径。然而在二维LBM-DEM耦合模型中,三维颗粒会被简化成二维圆形,相互紧密接触的圆形之间不存在任何通道,阻止了流体的流动,这与实际情况相悖。此外,三维颗粒与二维圆形的受力也不相同。因此在二维LBM-DEM耦合模型中需要对颗粒进行处理。

    Cui[22]从流体等效拖曳力出发,将三维球体颗粒等效为圆柱体,从而得到二维圆形颗粒的水力半径$ {{{r}}_{\rm{h}}} \approx 0.8R $,其中$ R $为颗粒的实际半径,该式适用于雷诺数较大($ {Re} > 1\;000 $)的情况。本文将推导较小雷诺数的情况。

    当雷诺数较小时,由斯托克斯阻力公式[23]可知三维中流体作用于圆球的力为:

    $$ F = 6\text{π} \mu UR $$ (4)

    同样,二维中,流体作用于圆盘的力为:

    $$ {F_{2{\text{D}}}} = 3\text{π} \mu {U_{2{\text{D}}}} $$ (5)

    式中:$ U $−实体颗粒相对于流体的速度;

    $ F $−三维实体颗粒收到流体作用力;

    $ {U_{2{\text{D}}}} $−二维圆形颗粒相对于流体的速度;

    $ {F_{2{\text{D}}}} $−二维圆形颗粒收到流体的作用力。

    因此当球体颗粒在静水中自由沉降时,阻力应与速度成正比增大,逐渐达到匀速运动,作用于球体的重力、浮力以及阻力保持平衡,将二维和三维计算中的平衡速度进行等效,则有:

    $$ {R_{\rm{h}}}{\text{ = }}0.816R $$ (6)

    本文建议LBM-DEM耦合模型模拟多孔介质时水力半径因子取:

    $$ {R_{\rm{h}}}{{ = 0.8\sim}}0.816R $$ (7)

    由于流体是在孔隙介质中流动,并对细小的孔隙介质产生力的作用,为保证流体运动和流体作用力达到精度要求,需要选择细小颗粒影响下的格子网络最优步长[24]。定义颗粒直径与格子步长比$ \;\beta = D/{\delta}x $,其中D为颗粒直径。本文通过平板间流体对圆形阻碍物作用力的收敛性分析,以选择计算合理的$ \beta $。数值模拟取平板长度为0.15 m,宽度为0.01 m,流体为20 °C下的水,其密度为1×103 kg/m3,动力黏滞系数μ为1.010×10−6 kPa·s。左端流入压强3.4 Pa,右端流出压强为3.3 Pa。颗粒直径为1 mm,位于x=0.070 m、y=0.005 m处。对3种不同直径D(1.0 mm、1.2 mm、1.5 mm)的颗粒,计算分析不同β下流体对颗粒的作用力,计算结果如图2所示。

    图  2  颗粒所受到水流的水平作用力随β的变化
    Figure  2.  The variation of horizontal force on solid particle introduced by water flow with different β

    颗粒受力计算值随着β的增大而逐渐收敛。当β≥8时,流体对颗粒作用力与参考值(最小步长0.05 mm工况)比较,相对误差可控制在0.5%以内。因此,对覆盖层多孔介质中进行流固耦合计算时,为确保计算的准确性和收敛性,在进行流体计算时,模型中经过水力半径因子缩放后的最小颗粒直径与格子步长比需满足β≥8。

    图3所示,覆盖型岩溶塌陷的流固耦合仿真计算共包括4个部分:DEM初始设置、LBM初始设置、塌陷模型初始设置、塌陷模型计算。

    图  3  覆盖型岩溶塌陷的建模流程
    Figure  3.  Modeling flowchart of covered karst ground collapse formation

    (1)DEM初始设置:锁定LBM程序;确定模型的范围和边界,形成模拟空间;模型的左右边界、下边界为Wall边界,模型内颗粒不能穿越此类边界。岩溶开口在前期设置中假设为Wall边界,在后期修改为开放边界。上边界为开放边界,即颗粒可以自由出入该边界。在空间内按照四边形规则排布或者随机生成颗粒,并赋予颗粒力学参数(杨氏模量和内摩擦角,为了颗粒间紧密,此时不设黏性键);引入局部阻尼[2526],在重力作用下将颗粒压实,经过足够长时间的循环,消除土体颗粒内部的不平衡力,让其达到相对稳定状态。

    (2)LBM初始设置:在DEM初始设置完成后,激活LBM计算程序,锁定DEM程序(固体颗粒位置固定);按照模拟范围生成离散网格;施加水力边界条件,计算固体颗粒间的水的流动情况,直至达到稳定;打开DEM程序,修改DEM程序中的时间步长,保证DEM和LBM之间能够实时的同步交换,在水力和重力共同作用下,运行足够步数循环,使得土体达到相对稳定状态。

    (3)塌陷模型初始设置:对已经生成的模型中的颗粒之间存在的接触按照参数黏性力学参数配置黏性键,经过足够长时间的循环,消除土体颗粒内部的不平衡力,完成颗粒流模型的黏性设置(黏性键只配置一次,后续计算中新生成的接触或者颗粒黏性键破坏后的接触不再拥有黏性属性);去除模型中岩溶开口处的墙体,改为开放边界。开口处的颗粒由于失去墙体的顶托力在重力作用和应力释放影响下产生细微变形,少部分颗粒向下脱落,并形成平衡拱,运行足够步数,使得颗粒达到相对稳定状态。至此,完成对覆盖型岩溶塌陷模型的初始状态设置。一般情况下,拆除岩溶开口处的墙体(Wall)后的计算属于动力计算,离散元系统中的阻尼应设置为0,但由于岩溶开口的形成一般是非常缓慢的扩展过程,而模型中为了方便,将岩溶开口处的墙体直接拆除,容易造成模型系统内颗粒的流失,与实际状态不符。因此,在此步骤下可以仍然考虑继续采用阻尼来避免这一问题。

    (4)塌陷模型计算:模型完成初始状态设置后,确保去掉阻尼(设为0);根据计算需求修改流体边界条件;进一步进行耦合计算至稳定或者塌陷。(图4

    图  4  覆盖型岩溶塌陷LBM-DEM模型
    Figure  4.  Schematic view of LBM-DEM model of covered karst ground collapse

    在实际中黏土颗粒非常小,模拟数量巨大,普通计算机能力有限而不能有效模拟。颗粒放大法是目前较为可行的处理方式,已广泛应用于工程研究中,该方法将颗粒放大,降低模型中的颗粒数目,使原物理模型问题能在合理有效的时间内解决[2730]。因此本文假设模型中的颗粒为多个黏性土颗粒结合的土体,以增加颗粒的半径。如图5所示,模型空间范围取15.9 cm × 8.0 cm,岩溶开口大小为0.6 cm,颗粒直径取0.15 cm,对应LBM中颗粒缩放后的直径为0.12 cm(水力半径因子取0.8),LBM空间步长取$ {\delta}x $=0.0001 m,即β为12,生成127.2万个网格,1 274 391个格点。固体颗粒与流体具体参数主要参照文献[17]、[22]、 [25](表1)。在LBM中水头的计算可以用等效水压代替,令隔水层上方水压为Pt,下方岩溶开口处水压力为Pb,其余边界为水量平衡边界或不透水边界。

    表  1  计算模型参数
    Table  1.  Summary of simulation model parameters
    参数名称
    固体(DEM)密度/(kg·m−32700
    杨氏模量/Pa100×106
    泊松比0.3
    内摩擦角/(°)20
    法向黏聚力/(N·m−1150
    切向黏聚力/(N·m−1150
    时间步长/s2.5×10−5
    重力加速度/(m·s−29.8
    流体(LBM)密度/(kg·m−31000
    运动黏滞系数/(m2·s−11.01×10−6
    时间步长/s0.0001
    空间步长/m0.0001
    下载: 导出CSV 
    | 显示表格
    图  5  承压水位下降工况覆盖型岩溶塌陷过程
    Figure  5.  The process of covered karst ground collapse formation with declining confined water level

    为了方便模拟,本次模拟工况为初始时刻模型中各处水头相等,等效成水压为933.33 Pa,Pb以降速50 Pa/s下降500 Pa的情况,模拟结果如图5所示。

    (1)覆盖型岩溶塌陷形成过程

    在初始时刻,模型中各处的水头相等,水流处于相对静止状态,对土体仅有向上的浮力作用。图5(a)为模型完成初始状态设置后的颗粒分布情况。

    随着承压水的水位下降,岩溶开口处的水头低于土层中的水头,土层中孔隙水向下流动,对土体颗粒造成向下的力。由于水位下降初期水流对岩溶开口上方土颗粒的作用力较小,对颗粒未造成较大影响,随着承压水位的继续下降,孔隙水对土颗粒的作用力逐渐增加,岩溶开口上方土体颗粒逐渐向洞口处移动,并脱离内侧相邻土体颗粒,见图5(b)。如图5(c)所示,随着水力的作用,岩溶上方失稳土体颗粒增多,土洞不断扩展,但是此时黏性土层表面并不会出现明显的下沉。土洞不断扩大,在自身重力以及水力作用下,表层土体受到破坏并突然塌陷进下方的洞穴,见图5(d),而形成覆盖型岩溶塌陷,此后坠入土洞的土体在水力作用下逐渐通过岩溶开口流失,见图5(e),并最终形成落水洞。

    比较图5(b)(c)(d)(e),土洞扩展由初期土体颗粒运动呈现单颗粒向下剥落,逐渐演变成块状失稳剥落,块体内的土体颗粒之间仅发生极小的相对位移甚至不发生相对运动,发生失稳破坏的块状构造不断变大,土洞扩展加速,最终导致土洞上方土层失稳。隔水层中的土洞从开始扩展到土洞顶部达到极限厚度是一个逐渐加速破坏的过程,顶板失稳塌陷是土洞扩展的特例。

    (2)塌陷过程中的水力特性

    在承压水位下降的影响下,岩溶开口处的颗粒优先发生失稳,发生较大运动,从而影响附近的水压分布。如图6所示为塌陷过程中岩溶开口附近高程水压随时间变化曲线,图6y表示为塌陷过程中岩溶开口中心线不同高度(底部y=0 mm)。发生颗粒剥落时,颗粒向下运移,原位置和上方附近水压有陡降,原因来自两方面:一方面是颗粒向下运动,原位置需要得到周围地下水补偿,从而引起周围水压释放;另一方面,颗粒失稳后会在重力作用下加速下降,这个速度一般会比水流速度大得多,会对颗粒附近的局部范围的水流造成扰动,加速水流,从而导致颗粒上方附近的水压降低。这使得上方颗粒附近形成较大的水压差,从而加大上方颗粒塌陷的可能。

    图  6  塌陷过程中岩溶开口中心线不同高程处水压随时间变化曲线
    Figure  6.  The variation of water pressure at different elevations of karst opening centerline during covered karst ground collapse formation

    (3)地下水对颗粒作用力分析

    对岩溶开口中心线上7个颗粒(0377、0588、0799、1854、3331、4808、6285,位置从下至上排序)受到地下水垂向作用力随时间变化作出曲线,如图7所示,其中0377、1854、3331、4808、6285五个颗粒在垂向上等距,0377、0588、0799为垂向上相邻的颗粒。

    图  7  塌陷过程中岩溶开口中心线颗粒受地下水垂向作用力随时间变化曲线
    Figure  7.  The variation of vertical force on the particles at karst opening centerline during covered karst ground collapse formation

    图7可以看出,岩溶开口附近的土体颗粒受力失稳剥落后,导致内部颗粒受到地下水的作用力明显增大。例如,颗粒0377在地下水垂向作用力为−0.073 2 N(负值表示方向向下)时失稳剥落,导致颗粒0588受到的地下水垂向作用力陡增至−0.144 N从而失稳,进一步导致颗粒0799上地下水垂向作用力突变为−0.129 2 N。颗粒一旦开始下降后,受到向下的水流作用力逐步减小。模型表层的颗粒因水流实现贯通,造成所受到的水流作用力紊动。

    由此不难推断出,当岩溶开口上方土洞周围存在强度较弱、甚至对土洞内部产生黏性拉力土体时,该部分土体可以在较小水力作用下失稳,但其失稳后会造成较强的瞬间水流作用力,容易使得土洞边沿土体失稳。土洞发展过程中,对土体颗粒的作用力会越来越大,从而导致剥落力呈增大趋势。

    (4)位移变化分析

    图8所示为塌陷过程中岩溶开口中心线不同高程处颗粒的垂向位移随时间变化曲线。从图8中可以看出,岩溶开口中心线上土体颗粒从下向上逐级剥落,高程越高的土体颗粒位移变化较为一致,这说明岩溶开口上方土层从下向上,以颗粒剥落逐步演变为块状剥落,直至盖层整体失稳,整个过程是加速过程。

    图  8  塌陷过程中岩溶开口中心线颗粒的垂向位移随时间变化曲线
    Figure  8.  The variation of vertical displacement of the particles at karst opening centerline during covered karst ground collapse formation

    (1)承压水位下降情况下地下水主要对岩溶开口处土层的颗粒产生向下的作用力。在克服土颗粒之间的黏结作用后,岩溶开口处土洞附近的颗粒在重力和水力作用下剥落。土体颗粒的剥落容易造成土颗粒原位置和上方位置处水压的陡降,从而造成较强的水力坡降,使得地下水对内部颗粒作用力陡增,容易进一步引起上方颗粒在地下水作用力和重力作用下失稳,从而使得承压水位下降引起的塌陷从土体颗粒失稳至土层塌陷过程加速。

    (2)本文所建立的流固耦合模型对进一步从细观上研究水力驱动的覆盖型岩溶地面塌陷的发育特征以及发育的临界条件、定量分析各影响因子之间的关系具有理论及实际意义。

  • 图  1   Bouzidi插值反弹边界示意图

    Figure  1.   Schematic view of Bouzidi interpolation bounce-back diagram

    图  2   颗粒所受到水流的水平作用力随β的变化

    Figure  2.   The variation of horizontal force on solid particle introduced by water flow with different β

    图  3   覆盖型岩溶塌陷的建模流程

    Figure  3.   Modeling flowchart of covered karst ground collapse formation

    图  4   覆盖型岩溶塌陷LBM-DEM模型

    Figure  4.   Schematic view of LBM-DEM model of covered karst ground collapse

    图  5   承压水位下降工况覆盖型岩溶塌陷过程

    Figure  5.   The process of covered karst ground collapse formation with declining confined water level

    图  6   塌陷过程中岩溶开口中心线不同高程处水压随时间变化曲线

    Figure  6.   The variation of water pressure at different elevations of karst opening centerline during covered karst ground collapse formation

    图  7   塌陷过程中岩溶开口中心线颗粒受地下水垂向作用力随时间变化曲线

    Figure  7.   The variation of vertical force on the particles at karst opening centerline during covered karst ground collapse formation

    图  8   塌陷过程中岩溶开口中心线颗粒的垂向位移随时间变化曲线

    Figure  8.   The variation of vertical displacement of the particles at karst opening centerline during covered karst ground collapse formation

    表  1   计算模型参数

    Table  1   Summary of simulation model parameters

    参数名称
    固体(DEM)密度/(kg·m−32700
    杨氏模量/Pa100×106
    泊松比0.3
    内摩擦角/(°)20
    法向黏聚力/(N·m−1150
    切向黏聚力/(N·m−1150
    时间步长/s2.5×10−5
    重力加速度/(m·s−29.8
    流体(LBM)密度/(kg·m−31000
    运动黏滞系数/(m2·s−11.01×10−6
    时间步长/s0.0001
    空间步长/m0.0001
    下载: 导出CSV
  • [1]

    ZALASIEWICZ J,WILLIAMS M,STEFFEN W,et al. The new world of the anthropocene[J]. Environmental Science & Technology,2010,44(7):2228 − 2231.

    [2]

    PRICE S J,FORD J R,COOPER A H,et al. Humans as major geological and geomorphological agents in the Anthropocene:The significance of artificial ground in Great Britain[J]. Philosophical Transactions Series A,Mathematical,Physical,and Engineering Sciences,2011,369(1938):1056 − 1084.

    [3] 李前银. 再论岩溶塌陷的形成机制[J]. 中国地质灾害与防治学报,2009,20(3):52 − 55. [LI Qianyin. Further study on formation mechanism of karst collaps[J]. The Chinese Journal of Geological Hazard and Control,2009,20(3):52 − 55. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2009.03.018

    LI Qianyin. Further study on formation mechanism of karst collaps[J]. The Chinese Journal of Geological Hazard and Control, 2009, 20(3): 52-55. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2009.03.018

    [4] 余政兴,金福喜,段选亮. 河床透-阻型岩溶塌陷形成机理[J]. 中国地质灾害与防治学报,2020,31(2):57 − 66. [YU Zhengxing,JIN Fuxi,DUAN Xuanliang. Formation mechanism of karst collapse with unconfined aquifer-aquitaed system in riverbed[J]. The Chinese Journal of Geological Hazard and Control,2020,31(2):57 − 66. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2020.02.08

    YU Zhengxing, JIN Fuxi, DUAN Xuanliang. Formation mechanism of karst collapse with unconfined aquifer-aquitaed system in riverbed[J]. The Chinese Journal of Geological Hazard and Control, 2020, 31(2)57-66(in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2020.02.08

    [5] 车增光,刘洪,喻永祥. 苏州市金庭镇蒋东岩溶塌陷地质条件及形成机理研究[J]. 华东地质,2021,42(1):85 − 92. [CHE Zengguang,LIU Hong,YU Yongxiang. Study on geological conditions and formation mechanism of karst collapse in Jiangdong Village,Jinting Town,Suzhou City[J]. East China Geology,2021,42(1):85 − 92. (in Chinese with English abstract)

    Che Zengguang, LIU Hong, YU Yongxiang. Study on geological conditions and formation mechanism of karst collapse in Jiangdong Village, Jinting Town, Suzhou City[J]. East China Geology, 2021, 42(1): 85-92. (in Chinese with English abstract)

    [6] 缪世贤,黄敬军,武鑫,等. 徐州岩溶地质调查及其发育特征分析[J]. 水文地质工程地质,2017,44(2):172 − 177. [MIAO Shixian,HUANG Jingjun,WU Xin,et al. Karst geological survey and analysis of its development characteristics in Xuzhou[J]. Hydrogeology and Engineering Geology,2017,44(2):172 − 177. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.2017.02.26

    MIAO Shixian, HUANG Jingjun, WU Xin, et al. Karst geological survey and analysis of its development characteristics in Xuzhou[J]. Hydrogeology and Engineering Geology, 2017, 44(2)172-177(in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.2017.02.26

    [7] 万志清,秦四清,李志刚,等. 土洞形成的机理及起始条件[J]. 岩石力学与工程学报,2003,22(8):1377 − 1382. [WAN Zhiqing,QIN Siqing,LI Zhigang,et al. Formation mechanism and initial condition of soil cavity[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22(8):1377 − 1382. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-6915.2003.08.028

    WAN Zhiqing, QIN Siqing, LI Zhigang, et al. Formation mechanism and initial condition of soil cavity[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(8): 1377-1382. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-6915.2003.08.028

    [8] 王滨,贺可强. 岩溶塌陷临界土洞的极限平衡高度公式[J]. 岩土力学,2006,27(3):458 − 462. [WANG Bin,HE Keqiang. Study on limit equilibrium height expression of critical soil cave of karst collapse[J]. Rock and Soil Mechanics,2006,27(3):458 − 462. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2006.03.024

    WANG Bin, HE Keqiang. Study on limit equilibrium height expression of critical soil cave of karst collapse[J]. Rock and Soil Mechanics, 2006, 27(3): 458-462. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2006.03.024

    [9] 陶小虎,赵坚,WANG Xiaoming,等. 地下水位变化对透-阻型岩溶塌陷影响的分析[J]. 中国岩溶,2017,36(6):777 − 785. [TAO Xiaohu,ZHAO Jian,WANG Xiaoming,et al. Analysis of seepage effect on the formation of sinkhole in unconfined aquifer-aquitard system caused by groundwater changes[J]. Carsologica Sinica,2017,36(6):777 − 785. (in Chinese with English abstract) DOI: 10.11932/karst2017y50

    TAO Xiaohu, ZHAO Jian, WANG Xiaoming, et al. Analysis of seepage effect on the formation of sinkhole in unconfined aquifer-aquitard system caused by groundwater changes[J]. Carsologica Sinica, 2017, 36(6): 777-785. (in Chinese with English abstract) DOI: 10.11932/karst2017y50

    [10] 孙金辉. 覆盖型岩溶塌陷临界参数模型试验与数值模拟研究[D]. 成都: 西南交通大学, 2011: 7 − 35

    SUN Jinhui. Study on critical parameterin cover karst collapseby model experimentand numerical simulation[D]. Chengdu: Southwest Jiaotong University, 2011: 7 − 35. (in Chinese with English abstract)

    [11]

    LEI Mingtang, GAO Yongli, JIANG Xiaozhen, et al. Experimental study of physical models for sinkhole collapses in Wuhan, China[C]//Sinkholes and the Engineering and Environmental Impacts of Karst. San Antonio, Texas, USA. Reston, VA: American Society of Civil Engineers, 2005: 91-102.

    [12]

    MAHMOUD A. Application of digital image cross correlation to study sinkhole collapse[J]. ISRN Soil Science,2013:1 − 6.

    [13]

    Y F,ZHOU. The mechanism of soil failures along cracks subjected to water infiltration[J]. Computers and Geotechnics,2014,55:330 − 341. DOI: 10.1016/j.compgeo.2013.09.009

    [14]

    ALI S M. Modelling the effect of void migration underneath landfill liner system [D]. University of Nottingham, 2003.

    [15] 金晓文,陈植华,曾斌,等. 岩溶塌陷机理定量研究的初步思考[J]. 中国岩溶,2013,32(4):437 − 446. [JIN Xiaowen,CHEN Zhihua,ZENG Bin,et al. Preliminary thinking of quantitative research on the mechanism of karst collapse[J]. Carsologica Sinica,2013,32(4):437 − 446. (in Chinese with English abstract)

    JIN Xiaowen, CHEN Zhihua, ZENG Bin, et al. Preliminary thinking of quantitative research on the mechanism of karst collapse[J]. Carsologica Sinica, 2013, 32(4): 437-446. (in Chinese with English abstract)

    [16]

    CUNDALL P A. A computer model for simulating progressive large scale movements in blocky rock systems[M]. In Proc Int Symp Rock Fracture. Nancy; ISRM. 1971: 2 − 8.

    [17] 周健, 贾敏才. 土工细观模型试验与数值模拟[M]. 北京: 科学出版社, 2008

    ZHOU Jian, JIA Mincai. Meso-model test and numerical simulation of geotechnical engineering[M]. Beijing: Science Press, 2008. (in Chinese)

    [18] 陈松贵. 宾汉姆流体的LBM-DEM方法及自密实混凝土复杂流动研究[D]. 北京: 清华大学, 2014: 44 − 50

    CHEN Songgui. Development of LBM-DEM for Bingham suspensions with application to self-compacting concrete[D]. Beijing: Tsinghua University, 2014: 44 − 50. (in Chinese with English abstract)

    [19]

    QIAN Y H,D'HUMIÈRES D,LALLEMAND P. Lattice BGK models for navier-stokes equation[J]. Europhysics Letters (EPL),1992,17(6):479 − 484. DOI: 10.1209/0295-5075/17/6/001

    [20]

    BOUZIDI M,FIRDAOUSS M,LALLEMAND P. Momentum transfer of a Boltzmann-lattice fluid with boundaries[J]. Physics of Fluids,2001,13(11):3452 − 3459. DOI: 10.1063/1.1399290

    [21]

    MEI Renwei, YU Dazhi, SHYY W, et al. Force evaluation in the lattice Boltzmann method involving curved geometry[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 2002, 65(4 Pt 1): 041203.

    [22]

    CUI Xilin. Numerical simulation of internal fluidisation and cavity evolution due to a leaking pipe using the coupled Dem-lbm technique[D]. Birmingham, West Midlands, UK: University of Birmingham, 2013.

    [23] 陈玉璞, 王惠民. 流体动力学[M]. 2版. 北京: 清华大学出版社, 2013.

    CHEN Yupu, WANG Huimin. Fluid dynamics[M]. 2nd ed. Beijing: Tsinghua University Press, 2013. (in Chinese)

    [24]

    LOMINÉ F,SCHOLTÈS L,SIBILLE L,et al. Modeling of fluid-solid interaction in granular media with coupled lattice Boltzmann/discrete element methods:Application to piping erosion[J]. International Journal for Numerical and Analytical Methods in Geomechanics,2013,37(6):577 − 596. DOI: 10.1002/nag.1109

    [25]

    LHL A,GN B,ZB A,et al. Hydro-mechanical modeling of sinkhole occurrence processes in covered karst terrains during a flood[J]. Engineering Geology,2019,260:105249. DOI: 10.1016/j.enggeo.2019.105249

    [26]

    HENTZ S. Modélisation d’une structure en béton armé soumise à un choc par la méthode des eléments discrets [D]. Grenoble: Université Grenoble 1 - Joseph Fourier, 2003.

    [27]

    MIKIO,SAKAI. Study on a large-scale discrete element model for fine particles in a fluidized bed[J]. Advanced Powder Technology,2012,23(5):673 − 681. DOI: 10.1016/j.apt.2011.08.006

    [28] 汤志刚,蔡承刚,王艳红,等. 基于光纤传感的石膏矿地面塌陷监测预警系统[J]. 中国地质灾害与防治学报,2022,33(5):93 − 101. [TANG Zhigang,CAI Chenggang,WANG Yanhong,et al. Monitoring and warning system for ground subsidence of gypsum mine based on fiber sensing[J]. The Chinese Journal of Geological Hazard and Control,2022,33(5):93 − 101. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.202108034

    TANG Zhigang, CAI Chenggang, WANG Yanhong, et al. Monitoring and warning system for ground subsidence of gypsum mine based on fiber sensing[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(5)93-101(in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.202108034

    [29] 何军,刘磊,黎清华,等. 隐伏岩溶区地下空间探测技术方法研究—以武汉市为例[J]. 水文地质工程地质,2020,47(6):47 − 56. [HE Jun,LIU Lei,LI Qinghua,et al. Techniques for detecting underground space in hidden karst region:Taking Wuhan as an example[J]. Hydrogeology & Engineering Geology,2020,47(6):47 − 56. (in Chinese with English abstract)

    HE Jun, LIU Lei, LI Qinghua, et al. Techniques for detecting underground space in hidden karst region: taking Wuhan as an example[J]. Hydrogeology & Engineering Geology, 2020, 47(6): 47-56. (in Chinese with English abstract)

    [30] 查甫生,刘从民,苏晶文,等. 铜陵市朝山地区岩溶塌陷形成条件与地面稳定性评价分析[J]. 地质论评,2020,66(1):246 − 254. [ZHA Fusheng,LIU Congmin,SU Jingwen,et al. Formation conditions of karst collapse and evaluation of ground stability in Chaoshan area of Tongling City[J]. Geological Review,2020,66(1):246 − 254. (in Chinese with English abstract) DOI: 10.16509/j.georeview.2020.01.018

    ZHA Fusheng, LIU Congmin, SU Jingwen, et al. Formation conditions of karst collapse and evaluation of ground stability in Chaoshan area of Tongling City[J]. Geological Review, 2020, 66(1): 246-254. (in Chinese with English abstract) DOI: 10.16509/j.georeview.2020.01.018

  • 期刊类型引用(13)

    1. 徐文正,卢书强,林振,周王敏. 联合InSAR与神经网络的范家坪滑坡形变监测及预测研究. 水文地质工程地质. 2025(02): 150-163 . 百度学术
    2. 刘鹏. 新一轮找矿突破战略行动背景下地质调查发展路径. 世界有色金属. 2024(04): 93-95 . 百度学术
    3. 刘伯元,王婷婷,陈慧玲. 刍议数字时代地勘行业参与生态文明建设路径. 世界有色金属. 2024(04): 182-184 . 百度学术
    4. 岳磊,刘昌义,丛晓明,唐彬元,付江涛,邢光延,雷浩川,赵吉美,吕伟涛,胡夏嵩. 基于InSAR技术的夏藏滩滑坡区地表变形监测与分析. 水文地质工程地质. 2024(03): 158-170 . 百度学术
    5. 王东仓,何景媛. 创新金融模式下生态保护修复路径探析. 世界有色金属. 2024(06): 147-149 . 百度学术
    6. 黄宝华,周利霞,孔祥侨. 基于PS-InSAR的建筑及道路动态沉降安全监测. 山东交通学院学报. 2024(02): 53-59 . 百度学术
    7. 刘鹏. 创新视角下的生态保护修复探讨. 世界有色金属. 2024(08): 174-176 . 百度学术
    8. 刘栋永. 地勘行业战略转型背景下科学事业单位项目成本核算探究. 世界有色金属. 2024(08): 181-183 . 百度学术
    9. 陈一超,赵文博. 基于地质调查转型发展背景下创新技术应用探讨. 世界有色金属. 2024(08): 188-191 . 百度学术
    10. 亓国涛,王荣,高晓琦. 地勘单位参与生态修复行业高质量发展的路径探析. 世界有色金属. 2024(10): 178-181 . 百度学术
    11. 王婷婷,刘伯元,孔萌. 浅析矿山生态修复治理存在问题及治理方向. 世界有色金属. 2024(14): 102-104 . 百度学术
    12. 张婉娜,高乃昌. 数据可视化背景下生态修复规划实现路径. 世界有色金属. 2024(14): 190-192 . 百度学术
    13. 马梓程,陈思,熊忠招,谢菲,孙天成,沈鹏,张思琦,田晓月,刘飞. 基于机器学习协同长时序主被动遥感数据的地质灾害易发性评价. 资源环境与工程. 2023(05): 587-598 . 百度学术

    其他类型引用(4)

图(8)  /  表(1)
计量
  • 文章访问数:  1283
  • HTML全文浏览量:  1051
  • PDF下载量:  213
  • 被引次数: 17
出版历程
  • 收稿日期:  2022-07-20
  • 修回日期:  2023-03-05
  • 录用日期:  2023-04-16
  • 网络出版日期:  2023-04-24
  • 刊出日期:  2024-01-31

目录

/

返回文章
返回