代碧波 孙丽军 赵兴东

（1.东北大学采矿地压与控制研究中心，辽宁沈阳，110004；2.金属矿山安全与健康国家重点实验室，安徽马鞍山，243000；3.中钢集团马鞍山矿山研究总院股份有限公司，安徽马鞍山，243000；4.武汉科技大学资源与环境工程学院，湖北武汉，430081）

缓倾斜薄矿体的非充填法回采会形成大量的采空区，这些空区群的稳定性直接关系到井下生产的安全［1-2］。故而对空区群实际三维形态的测定和稳定性评价至关重要［3］。

采空区的实际三维形态和准确的空间位置是影响空区稳定性数值模拟结果准确性的关键因素。过去，在对空区稳定性进行数值模拟时，往往将采空区的形态视为空区形成之前（即开采前）进行采场设计时确定的规则形态，部分则是依据传统的测量方法测定的简化空区形态进行模拟计算［3-5］。然而，采空区的实际形态往往在实际开采过程中，因测量、爆破等因素的影响而呈不规则状，加之采空区形成后受爆破震动的影响，以及空区围岩的冒落、片帮，致使空区的实际形态和开采设计中确定的相差甚远、这就势必严重影响空区稳定性数值模拟结果的精度和可靠性［6-8］。因此，如何精确地确定空区的实际形态和空间位置，并以此为基础建立空区群稳定性数值模拟模型进行数值模拟计算分析，是空区群稳定性数值模拟研究中应加以解决的一个重要问题［4］。

针对上述问题，本项目开展以精密探测获得的空区实际三维形态和空间位置……