汪志林，叶海旺，李子旋，雷 涛

（1.矿物资源加工与环境湖北省重点实验室，湖北 武汉430070；2.武汉理工大学 资源与环境工程学院，湖北 武汉430070）

近年来，受国家基建项目对水泥和砂石骨料需求剧增的影响，我国大规模露天灰岩矿山的数量也越来越多。然而，作为一种典型的沉积岩，在后期地质运动作用下，灰岩中往往分布着两组以上的裂隙，其稳定性更易受到裂隙组分布形式的影响，失稳机制也更加复杂。因此，如何预估含有两组及多组裂隙岩体边坡的破坏过程，揭示其破坏模式，采取有效的工程措施，是近年来相关学者关注的重要课题之一［1－4］。

针对裂隙的不同分布形式对岩质边坡破坏模式的影响成果丰硕［5－14］。但对双组裂隙边坡尚未开展过成体系的研究。究其原因是在传统数值模拟中，依赖于连续介质理论的强度折减难以重现边坡动态失稳过程。

针对上述问题，本文以山西省某灰岩矿山为背景，引入可更有效呈现边坡失稳过程的离散元法，基于MatDEM离散元软件的二次开发功能，建立了含双组裂隙边坡的数值计算模型，研究了双裂隙组灰岩边坡的渐进失稳过程。

1 工程概况

某石灰石矿山位于山西省朔州市怀仁县何家堡乡悟道村南西1 km至楼子口村北西1.2 km处，地理坐标东经112°57′26″～112°59′03″，北纬39°49′08″～39°50′35″，矿区北东-南西长3.23 km，北西-南东宽0.45 km，面积1.454 km2。矿区位于大同盆地西侧洪涛山山脉中部，区内山势陡峭，西高东低，地形高差较大，最低海拔1 175 m，最高海拔1 575 m，相对高差约400 m，沟谷呈近东西向分布，切割深度较大。……