自重作用下覆盖层的抗剪强度与岩溶塌陷高度的数值分析研究.doc

广 东 地 质 GUANGDONG GEOLOGY 覆盖层的抗剪强度与土洞塌陷高度关系的数值分析研究* 郭改梅　陶 波　邱向荣 （广东省地质建设工程集团公司，广州　510080） 摘　要　通过三维快速拉格朗日有限差分法分析了阳春市大河水库移民河西安置区典型土洞塑性区的分布情况。数值分析结果表明，研究区典型土洞覆盖层在自重作用下所引起的塑性区已发展至地面，因此存在较大塌陷的可能。以此为基础，通过数值分析方法，研究了在自重及均布外载联合作用下土洞塌陷的高度与覆盖层抗剪强度之间的关系。研究结果表明：仅在自重作用下，粘聚力与塌陷高度呈指数函数关系，而随内摩擦角的增加，塌陷高度随之呈线性下降；在自重与均布外载联合作用下，对于不同的粘聚力及内摩擦角，塌陷高度与均布外载呈线性关系，而且线性方程由粘聚力及内摩擦角共同决定。 关键词　自重　均布外载　抗剪强度　土洞　塌陷高度　数值分析 覆盖型岩溶指可溶岩石大部分被第四系沉积物覆盖，覆盖率一般在70%以上，仅局部零星见有基岩出露，覆盖层厚度一般小于30 m。覆盖层岩土体是岩溶塌陷的物质基础，其物理力学性质（抗剪强度、孔隙率、含水量、自重及渗透系数等）、结构组成及厚度等均对土洞塌陷的形成产生重要影响。本文采用三维快速拉格朗日有限差分法（FLAC3D）对广东省阳春市大河水库移民春城河西安置区的某一典型岩溶地面塌陷进行了数值分析研究，结果表明在自重作用下，覆盖层所具有的抗塌力（主要由抗剪强度提供）不足以克服致塌力（主要由自重作用产生），因而导致了岩溶地面塌陷；以此为基础，全面研究了仅在自重作用情况及与外载联合作用情况下，抗剪强度与塌陷高度之间的关系。 1 研究区土洞发育概况 研究区位于阳春复向斜岩溶盆地内，地貌以河谷冲积平原为主，地势平坦，地面高程一般为13.6～15.5 m，微向南东倾斜，坡度1～2°。表层被第四系覆盖，厚度5.00～50.00 m；下伏石炭纪黄龙组灰岩，节理裂隙发育，部分被方解石脉充填。在节理裂隙的交叉或密集带，岩溶较发育。覆盖层中土洞发育。据勘查报告，研究区共揭示土洞119个，土洞大小不一，埋深各不相同，且各土洞的覆盖层均较复杂。经综合分析表明，该移民安置区典型土洞高度约为4 m，埋深约13 m①。 本文2008年1月收到，3月改回。 *广东省地质勘查局地质科学研究基金资助项目(NO.2006211)。 ①广东省地质勘查局.广东省阳春市大河水库移民春城河西安置区地质灾害勘查报告,2005 2　数值分析模型与岩土参数 2.1　数值分析模型 本次研究采用美国ITASCA公司开发研制的三维快速拉格朗日有限差分法（FLAC3D）的数值分析软件。该模型（图1）有896个单元，1 850个节点。在模拟土洞所引发的岩溶塌陷过程中，只需将以坐标原点为圆心，以2 m为半径的圆所在范围内的单元置为null即可。圆形土洞顶部距地面13 m。由于这是一个平面应变对称问题，所以只取了一个半圆为研究对象。边界条件如下：在x=0 m及x＝15 m的面上约束x向变形；在y=0 m及y＝0.5 m的面上约束y方向的变形；在z=－15 m的面上约束z向变形，z=15 m的面为自由面。 图1　数值分析模型 2.2 岩土弹塑性物理力学参数 覆盖层主要为粉砂及粉质粘土。据室内试验，粉质粘土、粉砂及其灰岩的弹塑性物理力学参数见表1。 表1 岩土体物理力学参数 岩 性 容重（kN/m3） 弹性模量（MPa） 泊松比 粘聚力（MPa） 内摩擦角（°） 抗拉强度（MPa） 粉质粘土 18.5 9 0.35 0.018 12 0 粉砂 19.4 3.32 0.28 0.016 25 0 灰岩 27.0 20 000 0.23 30 35 1 3 抗剪强度与塌陷高度的关系 3.1 自重作用下抗剪强度与塌陷高度的关系 3.1.1 粘聚力与塌陷高度 由土洞所引发的塌陷与覆盖层的各弹塑性物理力学参数有关，如弹性模量、泊松比、抗剪强度等。在数值分析研究过程中，仅研究塌陷高度与覆盖层抗剪强度（粘聚力与内摩擦角）之间的关系。在覆盖层仅受自重作用的前提下，为研究塌陷高度与粘聚力的关系，只需将内摩擦角固定不变，分别将不同的粘聚力代入数值分析模型中进行计算，然后建立粘聚力与塌陷高度之间的关系曲线，即可分析塌陷高度随粘聚力变化的规律。为了更好地揭示抗剪强度指标与塌陷高度之间的变化规律，覆盖层仅取粉质粘土一层。在研究过程中，将内摩擦角固定为25°，粘聚力的取值分别为0.016 、0.018 、0.020 和0.022 MPa，直至0.084 MPa。 由图2可知，随粘聚力的增加，塌陷高度总体上呈指数函数下降的趋势。总体上可以将曲线分为3段：当粘聚力为0.016～0.022 MPa时，塌陷高度保持为14.3 m不变，这是由于土洞所