杨 宁 ，尹贤刚 ，林 尧

（1.长沙矿山研究院有限责任公司，湖南 长沙 410012；2.金属矿山安全技术国家重点实验室，湖南 长沙 410012；3.国家金属采矿工程技术研究中心，湖南长沙410012；4.洛阳坤宇矿业有限公司，河南 洛阳 471700）

0 引言

采空区稳定性分析可通过现场观测、模型试验以及数值模拟等手段进行[1-5]。第一种方法即在现场建立观测点，施工时工作人员进行追踪，察看其应力-应变随施工而变化的情况，该法可记录现场最真实的资料，总结变化规律得出的结论可指导后续工作，同时能够作为反分析的基础数据源，因此该法在有些情况下显得尤为重要。第二种方法需靠大量的资源以及复杂的结构条件去支撑其实现与真实情况一致的简易模型，通过对施工步骤进行模拟得知其发生改变的动态过程。该法取得的效果视觉冲击力大，但在试验过程中所受约束因素太多，不易实现。随着科学技术的进步，渐渐被方案三代替。数值模拟法是通过电脑中输入的各项数据模拟实际情况开展的一项试验。适应能力强，可实现多种复杂模型的模拟工作，可随时改动，变通性强，成本也比较低，因此受到广大学者和生产企业的青睐[6-8]。

研究结合某矿地下采空区的实际情况进行三维数值模拟，主要用于对采空区稳定性及其地压活动规律进行模拟研究，从而对该矿在开采扰动下的矿柱、上下盘围岩的应力、位移状态等进行模拟分析。

1 工程概况

该矿矿体形态为层状、似层状与透镜状，倾角为25°～88°，厚度为1～108 m。……