基于FLAC3D的西藏吉隆镇热玛村后山崩塌稳定性分析

多吉次旦，边巴次仁，普布扎西，格桑多吉，米玛扎西

（西藏自治区地质矿产勘查开发局第六地质大队， 中国拉萨 851400）

1 危岩基本特征

热玛村后山崩塌位于日喀则市吉隆县吉隆镇热玛村，坡脚处为热玛村聚居区。崩塌后缘高程约3365m，前缘以崩塌可能影响范围为界，前缘高程3115m，相对高差约250m。勘查区域斜坡宽度约350m，纵向长约160m，面积约5.6×104m2，斜坡坡向165°～175°，斜坡坡度65°～75°。根据现场勘查危岩位置分布、变形破坏特征及节理裂隙分布特征，将勘查区危岩分为11个危岩体［1］（见表1）。

表1 吉隆镇热玛村后山崩塌危岩特征

2 模型构建

本文计算模型的构建采用实地勘测、无人机测绘、三维激光扫描、崩塌调查表等一手资料确定了边界、岩性分带、产状等要素收集后，使用Rhinoceros griddle三维地质建模软件建立数值计算模型并完成单元体划分，克服了专业岩土分析软件FLAC3D在前处理复杂网格建立中的不足后，通过Rhinoceros-FLAC3D接口程序，将Rhinoceros建立的计算模型数据导入FLAC3D计算，充分发挥FLAC3D在岩土数值计算方面的优势。最终计算模型单元体采用四面体单元（见图1、图2），共有1378815个单元，257001个节点。模型x轴方向长263m，y方向宽160m，z轴方向高334m，底部标高为3119，顶部最高处标高为3453m。

图1 斜坡模型总视图

图2 斜坡模型正视图

根据计算结果及8个典型剖面中的分布情况，分析位移分布情况，判断模型中崩塌面的稳定性问题、不稳定情况以及可能发生崩塌破坏的区域。为考虑能全面反映模型内部应力应变分布情况和模型不同区域重要性，选取2处典型剖面，即剖面1和3（见图3）。……