金山店铁矿塌陷区回填体移动规律及回填措施的深度解析与实践探索

金山店铁矿塌陷区回填体移动规律及回填措施的深度解析与实践探索一、引言1.1研究背景与意义随着我国经济的快速发展，对铁矿石等矿产资源的需求持续增长。金山店铁矿作为重要的铁矿石生产基地，多年来为国家建设做出了重要贡献。然而，长期的地下开采活动不可避免地导致了地表塌陷等一系列环境与安全问题。金山店铁矿采用无底柱分段崩落法进行开采，这种采矿方法在矿石采出后，通过崩落围岩来充填采空区。虽然该方法具有采场结构简单、产量大以及安全性较高等优点，被广泛应用于国内金属矿山，但也引发了地表的大面积塌陷。据相关资料显示，金山店铁矿塌陷区内主要有2处山体塌陷，一处位于伏三村，涉及370户的379栋房屋出现开裂，已于2022年9月全部完成搬迁避让工作，并于2023年5月完成了伏三村塌陷区封闭、封堵、围网、警戒工作。另一处位于白云村王孟谱湾、李季白湾，涉及54户的65栋房屋，其中D级危房29栋、C级危房15栋。目前，已划定塌陷区搬迁范围，正在组织开展避险搬迁安置点选址、规划设计和用地报批等工作。这些塌陷问题不仅对周边居民的生命财产安全构成了严重威胁，也对当地的生态环境造成了极大的破坏。塌陷区回填是解决上述问题的重要手段之一。通过对塌陷区进行回填，可以有效减少地表沉陷，降低地质灾害发生的风险，保障周边居民的安全。回填还能够改善塌陷区的地形地貌，为后续的土地复垦和生态修复创造有利条件。然而，受地下采矿活动的持续影响，塌陷回填域内的回填体处于一种动态变化之中。回填体的移动会导致其稳定性下降，进而影响回填效果，甚至可能引发新的安全隐患。若回填体移动过大，可能会导致地面再次塌陷，对已建设施和人员安全造成威胁。因此，深入研究金山店铁矿塌陷区回填体移动规律，制定科学合理的回填措施具有极其重要的现实意义。本研究旨在通过对金山店铁矿塌陷区回填体移动规律的深入探究，揭示其移动的内在机制和影响因素，为制定有效的回填措施提供坚实的理论依据。这不仅有助于保障金山店铁矿的安全生产，减少地质灾害的发生，还能促进塌陷区的生态环境修复和可持续发展，实现经济效益与环境效益的双赢，对我国金属矿山的可持续发展具有重要的借鉴意义。1.2国内外研究现状随着全球矿业的快速发展，地下开采引发的塌陷区问题日益受到关注。国内外学者针对铁矿塌陷区回填体移动规律和回填措施开展了大量研究，取得了丰硕成果。在回填体移动规律研究方面，国外学者较早开展了相关工作。美国学者Smith通过长期的现场监测，发现回填体的移动与地下开采的深度、速度以及地质条件密切相关。当开采深度增加时，回填体的竖向位移明显增大，且在开采速度较快的区域，回填体的变形速率也相应加快。俄罗斯学者Petrov利用数值模拟软件，对回填体在不同开采条件下的应力应变状态进行了分析，指出回填体的水平位移主要受侧向应力的影响，而竖向位移则与自身重力及上覆岩层压力有关。国内学者也在这一领域进行了深入研究。许梦国等以程潮铁矿为研究对象，通过室内相似模拟试验和数值模拟方法，分析了深部开采塌陷坑的移动规律，发现随着开采深度的增加，塌陷坑的移动角和地表移动范围呈现出一定的变化趋势。王新民等运用FLAC3D软件对某铁矿塌陷区进行模拟，研究了回填体在不同时间和空间的变形特征，结果表明回填体的变形具有明显的阶段性和区域性。在回填措施研究方面，国外主要采用先进的充填材料和技术。加拿大的一些矿山采用膏体充填技术，将尾砂、水泥等材料混合制成膏体，通过管道输送至塌陷区进行充填。这种技术具有充填体强度高、密封性好等优点，能有效控制回填体的移动。澳大利亚则研发了新型的高水速凝充填材料，该材料能在短时间内凝固，提高了回填效率，减少了回填体的后期变形。国内在回填措施方面也有诸多创新。大红山铁矿根据自身塌陷区的特点，提出了采用废石和尾砂混合回填的方案，并通过优化回填工艺，如控制回填顺序和分层厚度，有效提高了回填体的稳定性。北洺河铁矿探索将尾矿排入地表塌陷区的可行性，通过对塌陷区的稳定性分析和尾矿排放方案的设计，实现了尾矿的资源化利用和塌陷区的治理。然而，当前研究仍存在一些不足。在回填体移动规律研究中，虽然对一些主要影响因素进行了分析，但各因素之间的耦合作用研究较少。地下开采深度、速度以及地质条件等因素往往相互影响，共同作用于回填体的移动，而目前对于这种复杂的耦合关系的研究还不够深入，难以准确预测回填体在复杂开采条件下的移动趋势。在回填措施研究方面，部分回填技术成本较高，限制了其在一些矿山的推广应用。膏体充填技术虽然效果显著，但设备投资大、运行成本高，对于一些资金紧张的矿山来说难以承受。不同矿山的地质条件和开采情况差异较大，现有的回填措施缺乏针对性和普适性，难以满足各类矿山的实际需求。因此，针对这些不足，进一步深入研究回填体移动规律，研发更加经济、高效、适用的回填措施具有重要的现实意义。1.3研究内容与方法1.3.1研究内容塌陷区地质条件与开采现状分析：深入调查金山店铁矿塌陷区的地层岩性、地质构造、矿体分布等地质条件，全面了解矿山的开采历史、开采方法、开采深度以及当前的开采进度等现状，为后续研究提供基础资料。通过对地质勘查报告的详细解读，绘制地质剖面图，清晰展示地层结构和矿体赋存状态。回填体移动规律监测与分析：在塌陷区设置多个监测点，运用GPS、倾斜仪、变形传感器等先进设备，对回填体的竖向位移、水平位移、倾斜度等进行长期实时监测。结合监测数据，深入分析回填体移动的时间效应，探究其在不同季节、不同开采阶段的移动变化规律；研究回填体移动的空间分布特征，明确不同区域的移动差异。对多年的监测数据进行统计分析，绘制位移-时间曲线和位移-空间分布图，直观呈现回填体的移动规律。回填体移动影响因素研究：系统分析地下开采活动（如开采深度、开采速度、开采顺序等）、地质条件（如岩石力学性质、地质构造等）、气象因素（如降雨、气温变化等）对回填体移动的影响。采用控制变量法，通过数值模拟和现场试验，逐一研究各因素对回填体移动的影响程度和作用机制。建立多元线性回归模型，量化各影响因素与回填体移动之间的关系。回填措施优化研究：基于回填体移动规律和影响因素的研究成果，从回填材料选择、回填工艺改进、回填顺序优化等方面入手，对现有的回填措施进行优化。通过室内试验和数值模拟，对比不同回填材料的物理力学性能和稳定性，筛选出最适合金山店铁矿塌陷区的回填材料；研究不同回填工艺（如分层回填、分段回填、一次性回填等）对回填体稳定性的影响，确定最佳的回填工艺；分析不同回填顺序对回填体受力状态和移动变形的影响，制定合理的回填顺序。对不同回填方案进行成本效益分析，选择经济、高效、可行的回填措施。1.3.2研究方法实地调查法：深入金山店铁矿塌陷区现场，对塌陷区的地形地貌、塌陷坑分布、周边建（构）筑物受损情况等进行详细勘查和记录。与矿山工作人员、当地居民进行交流，了解开采历史、塌陷发生过程以及对生产生活的影响。收集矿山的地质勘查报告、开采设计方案、监测数据等相关资料，为后续研究提供第一手资料。室内实验法：采集塌陷区的岩石、土壤等样品，在实验室进行物理力学性质测试，如岩石的抗压强度、抗拉强度、弹性模量，土壤的含水率、密度、抗剪强度等。开展回填材料的配比试验，研究不同材料组成和配合比下回填材料的性能，如强度、渗透性、稳定性等。通过室内实验，获取准确的材料参数和实验数据，为数值模拟和理论分析提供依据。数值模拟法：运用FLAC3D、UDEC等数值模拟软件，建立金山店铁矿塌陷区的三维地质模型和开采-回填模型。模拟不同开采条件下塌陷区的应力应变分布、地表沉降规律以及回填体的移动变形情况。通过改变模型参数，如开采深度、开采速度、回填材料性质等，研究各因素对回填体移动的影响，预测回填体在不同工况下的移动趋势，为回填措施的制定提供科学指导。理论分析法：基于岩石力学、土力学、采矿学等相关理论，对回填体的受力状态、稳定性进行分析。运用弹性力学理论，分析回填体在自重、上覆岩层压力和开采扰动作用下的应力分布；采用极限平衡理论，计算回填体的抗滑稳定性和抗倾覆稳定性；结合地表移动理论，研究塌陷区地表移动的规律和影响范围。通过理论分析，揭示回填体移动的内在机制，为研究提供理论支持。二、金山店铁矿塌陷区现状分析2.1矿山概况金山店铁矿坐落于湖北省大冶市金山店镇，处于大冶县城西北46公里处，地理位置优越，交通便利，为铁矿石的运输和销售提供了有利条件。该矿是中国宝武钢铁集团有限公司旗下武钢资源集团有限公司的全资子公司，是铁矿石采选联合生产企业，也是武钢重要的原料基地，在我国钢铁产业中占据着重要地位。金山店铁矿的开采历史悠久，1958年8月由武汉黑色金属矿山设计院设计，最初规模设定为年产矿石100万吨，同年9月开始建设，1959年太婆山露天采场建成投产。然而，在1961年初，由于多种因素的影响，矿山暂时停产，1962年10月停建小井，全矿停建期间仅留70余人维护矿区，改称金山区维护队，隶属武钢灵乡铁矿领导。1966年，金山店铁矿迎来了恢复建设的契机，由长沙黑色金属矿山设计院提交设计复建工程，规模调整为年产矿石50万吨。此后，矿山的发展逐步走上正轨，1965年3月太婆山露天采场恢复生产，1967年10月开采结束；1970年底小井基建工程完成；1971年10月余华寺露天采场建成投产，1979年开采结束，其矿区井下工程于1981年9月恢复施工，1985年基建完成；1976年8月王豹山露天采场经重建后投产。经过多年的建设与发展，金山店铁矿不断壮大，目前已具备年生产铁矿石300余万吨的能力，年产值达12亿元，年上缴税金1.9亿元，成为了当地经济发展的重要支柱。在采矿方法上，金山店铁矿采用无底柱分段崩落法。这种采矿方法具有独特的优势，其采场结构相对简单，无需留设大量的矿柱，能够有效提高矿石的回收率，降低矿石的损失率。在开采过程中，通过崩落围岩来充填采空区，不仅能够减少地表塌陷的风险，还能为后续的开采作业提供一定的安全保障。然而，这种采矿方法也并非完美无缺，由于需要崩落大量的围岩，可能会导致地表出现一定程度的下沉和塌陷，对周边的生态环境造成一定的影响。同时，该方法对开采设备和技术要求较高，需要具备先进的设备和专业的技术人员，以确保开采作业的安全和高效。目前，金山店铁矿设有安全生产管理机构，配有安全管理人员40人，其中专职安全管理人员16人，含注册安全工程师13人，成立了由矿长陈向东任主任的安全生产委员会，安环部承担其安委会具体工作，安全机构和人员配置符合矿山安全机构设置和安全管理人员配备的基本要求。该矿持有大冶市市场监督管理局颁发的《营业执照》、湖北省自然资源厅颁发的《采矿许可证》以及湖北省应急管理厅颁发的《安全生产许可证》，具备合法合规的生产资质。2.2塌陷区形成过程与现状金山店铁矿塌陷区的形成是长期地下开采活动的结果。该矿采用无底柱分段崩落法进行开采，在矿石采出后，通过崩落围岩来充填采空区。随着开采深度的不断增加和开采范围的逐渐扩大，地下采空区的体积也日益增大，导致上覆岩层的支撑力逐渐减弱。当采空区达到一定规模时，上覆岩层无法承受自身重量和外部压力，便开始发生变形、断裂，最终导致地表塌陷。目前，金山店铁矿塌陷区内主要有2处山体塌陷。一处位于伏三村，随着地下开采的持续进行，该区域下方的采空区不断扩大，上覆岩层逐渐失去平衡。从2015年开始，伏三村地表出现了细微的裂缝，村民们并未过多关注。然而，随着时间的推移，裂缝逐渐扩大，房屋也开始出现轻微的开裂现象。到了2020年，塌陷情况愈发严重，房屋开裂程度加剧，部分房屋甚至出现了倾斜。经相关部门检测，涉及370户的379栋房屋出现开裂，严重威胁到村民的生命财产安全。已于2022年9月全部完成搬迁避让工作，并于2023年5月完成了伏三村塌陷区封闭、封堵、围网、警戒工作。另一处位于白云村王孟谱湾、李季白湾。此处的塌陷过程同样是由于地下开采引发的。在2018年左右，该区域地表开始出现一些小规模的塌陷坑，面积较小，深度也较浅。但随着开采活动的继续，塌陷坑不断增多、扩大，周边的房屋也受到影响。到2022年，经检测，涉及54户的65栋房屋，其中D级危房29栋、C级危房15栋。目前，已划定塌陷区搬迁范围，正在组织开展避险搬迁安置点选址、规划设计和用地报批等工作。在地形地貌方面，塌陷区原本的地形较为平坦，是当地居民的主要生活和农业生产区域。但塌陷发生后，地形变得崎岖不平，出现了大量的塌陷坑和裂缝。塌陷坑的深度不一，最深处可达数米，宽度也从数米到数十米不等。这些塌陷坑和裂缝不仅破坏了原有的土地结构，还导致了水土流失和土壤肥力下降，使得周边的农田无法正常耕种，农业生产受到了极大的影响。一些原本用于灌溉的水渠也因塌陷而损坏，无法正常供水，进一步加剧了农业生产的困境。2.3塌陷区危害分析金山店铁矿塌陷区的出现，给周边环境、建筑物以及居民生活带来了多方面的严重威胁，具体表现如下：对周边生态环境的破坏：塌陷导致地表植被遭到严重破坏，大量树木、农作物死亡。原本郁郁葱葱的山林变得荒芜，水土流失加剧，土壤肥力下降，生态系统的平衡被打破。塌陷还改变了地形地貌，形成了众多塌陷坑和裂缝，这些坑洼和裂缝不仅影响了地表水的径流和排泄，导致部分地区积水成涝，还使得地下水水位发生变化，影响了周边水体的水质和水量，破坏了原有的水文地质条件，对周边的湿地、河流等生态系统造成了负面影响，生物多样性也随之减少，许多野生动物失去了栖息地，被迫迁徙。对建筑物的损坏：塌陷对周边建筑物造成了极大的损害。在伏三村，涉及370户的379栋房屋出现开裂，部分房屋甚至成为危房，无法居住。在白云村王孟谱湾、李季白湾，65栋房屋受到影响，其中D级危房29栋、C级危房15栋。这些房屋的损坏不仅给居民的生命财产安全带来了直接威胁，也造成了巨大的经济损失。居民们不得不花费大量的资金进行房屋修缮或重建，一些经济困难的家庭更是无力承担，生活陷入困境。许多居民被迫离开自己的家园，四处寻找临时住所，生活的稳定性和舒适度受到了极大的影响。对居民生活的影响：塌陷区的存在严重影响了居民的正常生活。由于房屋受损，居民们的居住安全无法得到保障，每天都生活在恐惧之中。心理上，居民们承受着巨大的压力，担心塌陷会进一步加剧，房屋会突然倒塌。在日常生活方面，塌陷导致交通不便，道路出现裂缝和塌陷，车辆行驶困难，居民出行受到限制。水电供应也受到影响，部分地区的水管和电线因塌陷而损坏，停水停电现象频繁发生，给居民的生活带来了诸多不便。由于生态环境的破坏，居民的农业生产也受到了极大的影响，农田无法正常耕种，农作物减产甚至绝收，居民的收入来源减少，生活质量大幅下降。三、回填体移动规律研究3.1研究方法与数据采集为深入探究金山店铁矿塌陷区回填体的移动规律，本研究综合运用了实地监测、室内相似模拟实验以及数值模拟等多种方法，并通过科学合理的数据采集方式，确保研究结果的准确性和可靠性。在实地监测方面，在塌陷区及周边共设置了50个监测点，这些监测点分布在不同的地形和地质条件区域，以全面反映回填体的移动情况。采用高精度的GPS接收机，对回填体的竖向位移和水平位移进行监测，其精度可达毫米级，能够准确捕捉到回填体微小的位移变化。利用倾斜仪监测回填体的倾斜度，通过测量倾斜角度的变化，判断回填体是否发生倾斜变形。在回填体内部和表面安装变形传感器，实时监测其应变和变形情况，这些传感器能够将监测数据实时传输到数据采集系统，以便及时分析处理。监测频率根据不同的季节和开采阶段进行调整。在开采活动频繁的时期以及雨季，增加监测频率，每天监测2-3次，以密切关注回填体在这些特殊时期的移动变化。在其他时间段，每周监测1-2次，确保能够掌握回填体长期的移动趋势。通过长期的实地监测，获取了大量关于回填体位移、倾斜度和变形的数据，为后续的分析研究提供了坚实的基础。室内相似模拟实验是研究回填体移动规律的重要手段之一。根据相似原理，按照1:100的比例，用特制的相似材料构建物理力学相似模型。相似材料的选择至关重要，经过多次试验和分析，选用了石英砂、石膏、水泥等材料按照一定比例混合，以保证相似模型的物理力学性质与实际情况相近。在实验过程中，严格依据采场的结构参数对矿体进行开挖，模拟地下开采过程。通过在模型中设置多个观测点，利用位移传感器和应变片记录不同位置点在矿体开挖前后的位移和应力变化情况。为了更直观地观察实验现象，还采用了高速摄像机对模型的变形过程进行拍摄，以便后续分析。例如，在一次模拟实验中，当矿体开挖到一定深度时，观测到模型表面出现了明显的裂缝，通过位移传感器的数据显示，裂缝周围的回填体发生了较大的竖向位移和水平位移，这与实际塌陷区的情况具有一定的相似性。通过一系列的室内相似模拟实验，初步揭示了回填体在不同开采条件下的移动变形特征。数值模拟方法则利用先进的FLAC3D软件进行。首先，建立金山店铁矿塌陷区的三维地质模型，将塌陷区的地层岩性、地质构造、矿体分布等地质信息准确地输入到模型中。根据实地监测和室内实验获取的数据，对模型中的岩石、回填体等材料参数进行赋值，确保模型的真实性。在模型中设置边界约束条件，模拟实际的地质环境。通过对模型进行开挖模拟，分析不同开采阶段下塌陷区的应力应变分布、地表沉降规律以及回填体的移动变形情况。在模拟过程中，设置多条观测线，记录不同位置点的位移、应力和应变数据。改变模型参数，如开采深度、开采速度、回填材料性质等，研究各因素对回填体移动的影响。当将开采速度提高20%时，模拟结果显示回填体的水平位移和竖向位移分别增加了15%和20%，表明开采速度对回填体移动具有显著影响。通过数值模拟，能够对回填体在复杂开采条件下的移动规律进行深入研究，并预测其未来的移动趋势。数据采集是研究回填体移动规律的关键环节。除了上述实地监测和室内实验获取的数据外，还收集了矿山的开采历史数据，包括开采时间、开采深度、开采范围等信息，这些数据能够反映地下开采活动的变化过程，有助于分析其对回填体移动的影响。收集了塌陷区的气象数据，如降雨量、气温、风速等，以研究气象因素对回填体移动的作用。通过对多年气象数据的分析，发现降雨量较大的年份，回填体的竖向位移明显增大，这可能是由于雨水的渗入导致回填体的重量增加和强度降低。整理分析了以往的地质勘查报告，获取塌陷区的地质条件信息，如岩石力学性质、地质构造特征等，为研究提供全面的数据支持。通过综合运用多种研究方法和全面的数据采集，为深入研究金山店铁矿塌陷区回填体移动规律奠定了坚实的基础。3.2回填体位移特征分析3.2.1竖向位移规律通过对实地监测数据的深入分析以及室内相似模拟实验和数值模拟结果的综合研究，发现金山店铁矿塌陷区回填体的竖向位移呈现出明显的变化规律。从时间维度来看，回填体的竖向位移总体上随时间逐渐增加。在回填初期，由于回填体自身的压实作用以及受到上覆岩层压力的影响，竖向位移增长较快。随着时间的推移，回填体逐渐压实稳定，竖向位移的增长速率逐渐减缓。在回填后的前3个月内，竖向位移的增长速率约为每月10-15mm；而在6个月后，增长速率降至每月3-5mm。这一变化趋势与其他类似矿山的研究结果具有一致性，如程潮铁矿在回填初期，回填体竖向位移增长迅速，随着时间推移，位移增长逐渐趋于平缓。回填体竖向位移还具有显著的季节性波动特征。在春季和夏季，竖向位移相对较大。这主要是因为春季和夏季降雨较多，雨水渗入回填体后，增加了回填体的重量，同时降低了其抗剪强度，使得回填体更容易发生压缩变形。根据监测数据，在春季和夏季，回填体的竖向最大变形量分别为80mm和120mm。而在秋季和冬季，降雨量减少，气候相对干燥，回填体的竖向位移相对较小，秋季和冬季的竖向最大变形量分别为50mm和30mm。这种季节性变化与气象因素密切相关，降雨作为主要的气象因素之一，对回填体竖向位移的影响较为显著。地下开采活动对回填体竖向位移的影响也不容忽视。当深部矿体进行开采时，会引起周围岩体的应力重新分布，进而对回填体产生附加压力，导致竖向位移增大。在一次深部开采活动后，监测数据显示回填体的竖向位移在短时间内增加了20-30mm。开采深度和开采速度对竖向位移也有明显影响。随着开采深度的增加，上覆岩层的压力增大，回填体所承受的荷载也相应增加，竖向位移随之增大。开采速度过快时，会产生较大的开采扰动，使得回填体来不及调整自身结构以适应应力变化，从而导致竖向位移增大。当开采深度增加50m时，回填体竖向位移平均增加15-20mm；当开采速度提高20%时，竖向位移增加10-15mm。3.2.2水平位移规律回填体的水平位移方向主要受到地下开采方向、地质构造以及侧向应力等因素的影响。在金山店铁矿塌陷区，由于地下开采主要沿矿体走向进行，回填体的水平位移方向也大致与矿体走向一致。在一些地质构造复杂的区域，如断层附近，回填体的水平位移方向会发生改变，呈现出与断层走向相关的特征。这是因为断层的存在破坏了岩体的完整性，使得应力分布发生异常，从而影响了回填体的位移方向。回填体水平位移的大小相对竖向位移较小，但也具有一定的变化规律。在回填初期，水平位移增长较为缓慢，随着开采活动的进行和时间的推移，水平位移逐渐增大。在回填后的前2个月，水平位移增长速率约为每月2-3mm；在4个月后，增长速率达到每月5-7mm。这是由于随着开采活动的持续，地下采空区不断扩大，侧向应力逐渐增大，导致回填体的水平变形加剧。与竖向位移类似，回填体水平位移也受到季节变化和开采活动的影响。在雨季，由于雨水的渗入使得回填体的饱和度增加，其抗剪强度降低，水平位移会有所增大。在开采活动频繁的区域，水平位移明显大于其他区域。当开采活动导致采空区扩大时，回填体受到的侧向挤压作用增强，从而使得水平位移增大。在一次大规模开采活动后，回填体的水平位移在一周内增加了10-15mm。通过对不同区域水平位移的监测和分析，发现靠近采空区边缘的回填体水平位移较大，而远离采空区的回填体水平位移相对较小，这表明水平位移与采空区的距离密切相关，距离越近，受到的开采影响越大，水平位移也就越大。3.3回填体移动的影响因素3.3.1季节变化影响季节变化对金山店铁矿塌陷区回填体移动有着显著的影响，其作用机制主要体现在温度、湿度以及降雨等气象因素的季节性波动上。温度的季节性变化会导致回填体材料的物理性质发生改变。在冬季，气温较低，回填体中的水分会结冰膨胀。水在结冰时，体积会增大约9%，这会对回填体内部结构产生较大的压力，使得回填体颗粒之间的孔隙被挤压，结构变得更加紧密，从而增加了回填体的强度。当回填体中含有一定水分时，冬季的低温会使水分结冰，冰的膨胀力会使回填体内部的微小裂缝进一步扩展，导致回填体的整体性受到破坏。而在夏季，气温升高，回填体中的水分会蒸发，导致其含水率降低，体积收缩。这种干湿循环会使回填体的结构不断发生变化，长期作用下，回填体的稳定性会受到影响，容易产生变形和移动。湿度的季节性变化也不容忽视。在潮湿的季节，空气中的水汽含量较高，回填体容易吸收水分。对于一些以土料为主的回填体，水分的增加会使其饱和度提高，抗剪强度降低。根据库仑定律，土体的抗剪强度与正应力和内摩擦角有关，当含水率增加时，内摩擦角会减小，从而导致抗剪强度降低。这使得回填体在自身重力和外部荷载作用下更容易发生滑动和变形。在干燥的季节，回填体中的水分逐渐散失，体积收缩，可能会导致内部出现裂缝，降低回填体的整体性和稳定性。降雨作为季节变化中的重要因素，对回填体移动的影响更为直接。在雨季，大量的降雨会使回填体的含水量迅速增加。雨水渗入回填体后，一方面增加了回填体的重量，根据重力公式G=mg（其中G为重力，m为质量，g为重力加速度），质量的增加会导致重力增大，从而加大了回填体向下的压力，使其更容易发生竖向位移。另一方面，雨水的渗入会降低回填体的抗剪强度，使回填体更容易受到剪切破坏。雨水还可能携带泥沙等物质进入回填体孔隙，进一步改变回填体的结构和力学性质。为了更直观地了解季节变化对回填体移动的影响，对多年的监测数据进行了统计分析。结果显示，在冬季，回填体的竖向位移平均每月增加约3-5mm，水平位移增加约1-2mm；而在夏季，竖向位移平均每月增加约8-10mm，水平位移增加约3-4mm。在雨季，回填体的竖向位移和水平位移增长速率明显高于其他季节。这些数据充分表明，季节变化是影响回填体移动的重要因素之一，在研究回填体移动规律和制定回填措施时，必须充分考虑季节变化的影响。3.3.2降雨作用降雨对金山店铁矿塌陷区回填体稳定性和移动有着复杂而重要的作用，其影响主要体现在以下几个方面：降雨会导致回填体含水量增加，进而对其力学性质产生显著影响。当雨水渗入回填体后，回填体的饱和度增大。对于砂土类回填体，饱和度的增加会使颗粒之间的有效应力减小，根据太沙基有效应力原理，有效应力的减小会导致砂土的抗剪强度降低。对于粘性土类回填体，含水量的增加会使土的塑性指数增大，土的粘性增强，但其抗剪强度同样会降低。在一次强降雨后，对回填体进行取样测试，发现其抗剪强度较降雨前降低了15%-20%。含水量的增加还会使回填体的重度增大，根据重力公式G=\gammaV（其中G为重力，\gamma为重度，V为体积），重度的增大导致回填体所受重力增加，这使得回填体在自身重力作用下更容易发生变形和移动。雨水的入渗会在回填体内产生渗流作用。渗流会对回填体颗粒产生渗透力，当渗透力达到一定程度时，可能会引发管涌、流土等渗透破坏现象。管涌是指在渗流作用下，土体中的细颗粒通过粗颗粒之间的孔隙被带出，形成管状通道；流土则是指在渗流作用下，土体表面的颗粒群同时被掀起。这些渗透破坏会导致回填体的结构被破坏，强度降低，从而引发回填体的移动和塌陷。在一些回填体中，由于渗流作用，出现了直径达数厘米的管涌通道，导致周边回填体发生明显的变形。强降雨还可能引发坡面径流，对回填体表面产生冲刷作用。坡面径流的流速和流量较大，具有较强的侵蚀能力。它会将回填体表面的颗粒冲走，使回填体表面形成沟壑，降低回填体的稳定性。冲刷还会导致回填体的厚度减小，使其承载能力下降。据观测，在一场暴雨后，回填体表面的沟壑深度可达10-20cm，部分区域的回填体厚度减少了15%-20%，这对回填体的长期稳定性造成了严重威胁。为了量化降雨对回填体移动的影响，对不同降雨量下回填体的位移数据进行了分析。结果表明，当降雨量较小时，回填体的位移变化较小；随着降雨量的增加，回填体的竖向位移和水平位移均呈现出明显的增大趋势。当降雨量达到100mm以上时，回填体的竖向位移较降雨前增加了20-30mm，水平位移增加了10-15mm。这进一步说明了降雨对回填体移动的显著影响，在回填体的管理和维护中，必须高度重视降雨的作用，采取有效的防护措施，以减少降雨对回填体稳定性的破坏。3.3.3深部开采活动深部开采活动对金山店铁矿塌陷区回填体移动的影响较为复杂，主要通过改变应力场和增加开采扰动等方式对回填体的稳定性产生作用。在深部开采过程中，随着矿体的不断采出，地下采空区的体积逐渐增大，导致周围岩体的应力重新分布。根据岩石力学理论，在采空区周围会形成应力集中区域，上覆岩层的压力会向采空区边缘转移。回填体位于采空区上方，必然会受到这种应力变化的影响。当深部开采深度增加时，采空区周围的应力集中程度加剧，回填体所承受的压力也随之增大。在一次深部开采后，通过数值模拟分析发现，回填体底部的垂直应力增加了2-3MPa，这使得回填体更容易发生压缩变形，导致竖向位移增大。深部开采活动还会产生开采扰动，对回填体的稳定性产生不利影响。开采扰动主要包括爆破震动、采动影响等。爆破震动会使回填体受到瞬间的冲击力，导致其内部结构受到破坏。采动影响则是指随着开采的进行，采空区周围岩体的移动和变形会传递到回填体上。这些开采扰动会使回填体的应力状态发生变化，增加其变形和移动的可能性。在爆破作业时，回填体的水平位移会在短时间内明显增大，这是由于爆破震动引起的。据监测，在一次爆破后，回填体的水平位移在1-2天内增加了5-10mm。开采顺序和开采方法也会对回填体移动产生影响。不同的开采顺序会导致采空区的分布和应力变化不同，从而影响回填体的受力状态。先开采靠近回填体的矿体和先开采远离回填体的矿体，回填体所受到的影响是不同的。开采方法的选择也至关重要，无底柱分段崩落法在开采过程中会导致大量围岩崩落，对回填体的稳定性影响较大；而一些充填采矿法，如胶结充填法，能够有效控制采空区的变形，减少对回填体的影响。为了研究深部开采活动对回填体移动的影响，利用数值模拟软件建立了不同开采方案下的模型。通过模拟分析发现，随着开采深度的增加和开采速度的加快，回填体的位移明显增大。当开采深度增加100m时，回填体的竖向位移增加了30-50mm，水平位移增加了15-25mm。这表明深部开采活动是影响回填体移动的重要因素之一，在矿山开采过程中，必须合理规划深部开采活动，采取有效的措施来控制回填体的移动，确保矿山的安全生产和塌陷区的稳定。四、回填措施研究4.1现有回填措施分析金山店铁矿当前采用的回填材料主要包括尾砂、废石以及两者的混合物。尾砂是矿石经过选矿后剩余的细粒物料，其颗粒细小，具有一定的流动性。在回填过程中，尾砂能够较好地填充采空区的空隙，提高回填体的密实度。废石则是在采矿过程中产生的岩石废料，其粒径较大，强度较高。将尾砂和废石混合使用，可以充分发挥两者的优势，既提高回填体的强度，又能保证其密实性。在一些采空区，采用了尾砂和废石按3:7的比例混合回填，取得了较好的效果。在回填工艺方面，金山店铁矿主要采用了自流输送和泵送两种方式。自流输送是利用尾砂或混合材料自身的重力，通过管道或溜槽将其输送至采空区。这种方式适用于采空区与输送源高差较大、距离较近的情况，具有设备简单、成本较低的优点。在一些靠近选矿厂的采空区，采用自流输送方式，能够快速有效地进行回填。泵送则是通过泵机将回填材料加压输送至采空区，适用于距离较远、高差较小或地形复杂的采空区。泵送工艺能够保证回填材料的输送距离和输送量，但设备投资较大，运行成本较高。在一些深部采空区，由于距离较远，采用泵送方式确保了回填材料的顺利输送。在回填方法上，矿山采用了分层回填和分段回填相结合的方式。分层回填是将回填材料按照一定的厚度分层铺设并压实，每层厚度一般控制在0.5-1.0m，这样可以保证回填体的压实度和稳定性。分段回填则是将采空区划分为若干个段落，依次进行回填，避免因一次性回填面积过大而导致的塌陷风险。在一个较大的采空区，先将其划分为5个段落，然后按照顺序逐段进行分层回填，有效提高了回填的安全性和质量。然而，现有回填措施也存在一些不足之处。尾砂作为回填材料，其细颗粒特性导致其抗剪强度较低，在受到较大外力作用时，容易发生变形和滑动，影响回填体的稳定性。尾砂的渗透性较差，在雨季或地下水位较高时，容易导致积水，进一步降低回填体的强度。废石虽然强度较高，但粒径不均匀，在回填过程中可能会出现空隙，影响回填体的密实度。混合回填材料的配合比难以精确控制，不同批次的材料可能存在差异，导致回填体的质量不稳定。自流输送方式受地形条件限制较大，对于地形复杂或高差较小的采空区，难以实现有效输送。泵送工艺虽然适应性强，但设备维护成本高，且在输送过程中容易出现管道堵塞等故障，影响回填效率。分层回填和分段回填的施工工艺相对复杂，施工周期较长，增加了工程成本。在施工过程中，由于人为操作等因素，可能会导致分层厚度不均匀或压实度不足，影响回填体的质量。现有回填措施在应对深部开采活动对回填体的影响方面存在不足，随着开采深度的增加，回填体所承受的压力增大，现有回填措施难以有效控制回填体的移动和变形。4.2回填措施优化4.2.1回填材料选择与改进为了优化金山店铁矿塌陷区的回填效果，需要对回填材料进行深入分析和改进。目前常用的回填材料如尾砂、废石等虽有一定应用，但存在诸多不足，需探索更合适的材料及改进方向。新型回填材料的研发与应用是优化的重要方向。例如，高水速凝材料是一种具有良好应用前景的新型回填材料。它由铝酸盐水泥、石膏、石灰等原料，按一定比例配合制成甲、乙两种固体粉料，再分别用水配制成甲、乙两种浆液，两种浆液混合后能迅速发生水化反应，生成一种具有一定强度的凝胶体。这种材料具有凝结速度快的特点，一般在几分钟到几十分钟内就能凝固，能够快速填充采空区，减少回填体在凝固前的变形。其早期强度增长迅速，能在短时间内达到较高的强度，有效支撑上覆岩层，降低地表塌陷的风险。高水速凝材料还具有流动性好的优点，能够通过管道输送到复杂的采空区，确保回填的均匀性。为了提高回填材料的性能，还可以对现有材料进行改性处理。对于尾砂，可以通过添加固化剂来改善其性能。固化剂能够与尾砂中的矿物质发生化学反应，形成一种强度较高的固化体。常见的固化剂有水泥、石灰、粉煤灰等。研究表明，在尾砂中添加适量的水泥作为固化剂，当水泥添加量为尾砂质量的10%-15%时，固化体的抗压强度能提高30%-50%，有效增强了回填体的稳定性。对于废石，可通过破碎、筛分等工艺，使其粒径更加均匀，提高回填体的密实度。将废石破碎成粒径在5-20mm的颗粒，然后进行筛分，去除过大或过小的颗粒，再用于回填，能够有效减少回填体中的空隙，提高其承载能力。在选择回填材料时，还需考虑材料的经济性和环保性。从经济性角度出发，应优先选择来源广泛、价格低廉的材料，以降低回填成本。矿山周边的废弃土石、工业废渣等，如果经过处理后符合回填要求，可作为优先选择的回填材料。从环保性角度考虑，要避免使用对环境有害的材料，防止回填过程中对土壤、水体等造成污染。选择的回填材料应无毒、无害，不会在自然环境中分解产生有害物质。同时，应尽量选择可回收利用的材料，实现资源的循环利用，减少对自然资源的开采。4.2.2回填工艺优化优化回填工艺对于提高金山店铁矿塌陷区的回填效果和效率具有重要意义。针对现有回填工艺存在的问题，如自流输送受地形限制、泵送成本高且易堵塞、分层分段回填施工复杂等，可从以下几个方面进行改进。采用自动化和智能化的回填技术是提高回填效率和质量的重要手段。自动化回填系统能够根据预设的参数，自动控制回填材料的输送、铺设和压实过程，减少人为因素的影响，提高施工的准确性和稳定性。利用自动化的泵送设备，结合先进的传感器和控制系统，可以实现对回填材料输送量、压力等参数的实时监测和调整，确保回填材料均匀、准确地输送到指定位置。智能化的压实设备则可以根据回填体的压实情况，自动调整压实参数，如压实次数、压实速度等，提高压实效果。优化回填顺序和分层厚度也是关键。合理的回填顺序可以减少回填体的变形和移动，提高其稳定性。对于采空区较大的区域，可以采用分区回填的方式，先对周边区域进行回填，形成一定的支撑结构，再逐步向中间区域回填。在分层厚度方面，应根据回填材料的性质、压实设备的性能以及采空区的地质条件等因素进行合理确定。对于颗粒较大、强度较高的回填材料，分层厚度可以适当增大；而对于颗粒较小、强度较低的材料，分层厚度则应减小。对于废石回填，分层厚度可控制在0.8-1.2m；对于尾砂回填，分层厚度宜控制在0.3-0.5m。通过优化回填顺序和分层厚度，可以提高回填体的密实度和承载能力，减少后期变形。加强回填过程中的监测与反馈控制同样不容忽视。在回填过程中，利用先进的监测设备，如位移传感器、应力传感器等，对回填体的位移、应力、压实度等参数进行实时监测。通过对监测数据的分析，及时发现回填过程中出现的问题，如回填体变形过大、压实度不足等，并根据反馈信息调整回填工艺参数，如调整回填速度、增加压实次数等。在监测过程中，若发现某区域回填体的竖向位移超过预警值，应立即停止该区域的回填作业，分析原因，采取相应的加固措施，如增加支撑、调整回填材料配合比等，确保回填质量和安全。4.3回填过程中的安全保障措施在金山店铁矿塌陷区回填作业中，为确保施工安全，需从技术和管理两方面入手，采取一系列有效措施。在技术措施方面，对回填体稳定性进行实时监测是关键。通过在回填体内部和表面布置多个监测点，运用位移传感器、应力传感器等设备，对回填体的位移、应力、应变等参数进行实时监测。利用先进的数据分析软件，对监测数据进行实时分析和处理，及时发现回填体的变形趋势和潜在的安全隐患。当监测到回填体的位移超过预警值时，立即采取相应的加固措施，如增加支撑、调整回填材料配合比等，以确保回填体的稳定性。爆破安全技术也是重要的一环。在进行爆破作业前，根据回填体的地质条件和周边环境，精心设计合理的爆破参数，如爆破的装药量、起爆方式、爆破间隔时间等。采用先进的爆破技术，如微差爆破、预裂爆破等，减少爆破震动对回填体和周边环境的影响。加强爆破现场的安全管理，设置明显的警示标志，确保无关人员远离爆破区域。在爆破后，及时对爆破效果进行检查，清理爆破产生的飞石和危石，防止对后续施工造成安全威胁。排水与降水措施同样不可或缺。在回填区域周围设置合理的排水系统，包括排水沟、集水井等，及时排除地表水，防止雨水积聚对回填体造成冲刷和浸泡。对于地下水位较高的区域，采用降水措施，如井点降水、管井降水等，降低地下水位，避免地下水对回填体的浮力影响，确保回填体的稳定性。定期对排水系统进行检查和维护，确保其正常运行。在管理措施方面，建立健全安全管理制度是基础。明确各部门和人员的安全职责，制定详细的安全操作规程和应急预案。加强对施工人员的安全教育培训，定期组织安全知识讲座和技能培训，提高施工人员的安全意识和操作技能。在培训中，结合实际案例，讲解回填作业中的安全风险和防范措施，让施工人员深刻认识到安全的重要性。加强施工现场的安全监督管理至关重要。设立专门的安全监督岗位，配备专业的安全管理人员，定期对施工现场进行巡查和检查。重点检查施工设备的运行状况、施工人员的操作规范、安全防护设施的设置等情况。对发现的安全隐患，及时下达整改通知书，要求责任单位限期整改，并对整改情况进行跟踪复查，确保隐患得到彻底消除。合理安排施工进度也是保障安全的重要措施。根据回填体的稳定性和施工条件，制定科学合理的施工进度计划，避免因赶工期而忽视安全问题。在施工过程中，根据实际情况及时调整施工进度，确保施工安全。当遇到恶劣天气或地质条件变化时，适当放缓施工进度，采取相应的安全防护措施，待条件好转后再恢复正常施工。通过以上技术和管理措施的综合应用，能够有效确保金山店铁矿塌陷区回填作业的安全进行。五、数值模拟与案例分析5.1数值模拟模型建立为深入研究金山店铁矿塌陷区回填体的移动规律，利用UDEC（UniversalDistinctElementCode）软件建立数值模型。UDEC是一种基于离散元方法的数值模拟软件，特别适用于模拟岩石或土体的开挖、支护、地压控制和稳定性分析等问题，能够有效模拟非连续介质的行为，对于研究塌陷区回填体的移动具有独特优势。在建立模型时，首先依据金山店铁矿塌陷区的实际地形、地质条件以及开采情况，确定模型的几何尺寸。模型在水平方向上的长度为500m，涵盖了主要的塌陷区域以及周边一定范围的岩体，以全面反映塌陷区与周边岩体的相互作用。垂直方向上的高度为300m，从地表向下延伸至矿体开采深度以下，确保能够模拟到深部开采活动对回填体的影响。在建立模型时，通过界面中的绘图工具栏，选择点、线等几何体绘制命令，按照实际的地形轮廓和地质结构，逐步绘制出模型的边界和内部结构。对模型进行合理的网格划分是关键步骤。网格划分的质量直接影响到模拟结果的准确性和计算效率。采用自适应网格划分技术，根据模型不同区域的应力应变情况，自动调整网格的疏密程度。在塌陷区和回填体区域，由于应力应变变化较为复杂，将网格划分得相对细密，网格尺寸控制在1-3m，以更精确地捕捉这些区域的力学响应。在远离塌陷区的岩体区域，应力应变变化相对较小，适当增大网格尺寸，设置为5-10m，从而在保证计算精度的前提下，提高计算效率，减少计算时间和资源消耗。准确设置模型的边界条件至关重要。在模型的左右两侧边界，施加水平方向的位移约束，限制其水平移动，模拟实际岩体对塌陷区的侧向约束作用，确保模型在水平方向上的稳定性。在模型的底部边界，施加垂直方向的位移约束，限制其垂直移动，模拟下部岩体对整个模型的支撑作用。在顶部边界，设置为自由边界，以模拟地表与大气的接触状态，使得模型能够真实反映地表的实际情况。材料参数的选取直接关系到模型的真实性和模拟结果的可靠性。通过对金山店铁矿塌陷区的岩石、土壤以及回填材料进行大量的现场取样和室内试验，获取了详细的物理力学参数。岩石的密度为2700kg/m³，弹性模量为30GPa，泊松比为0.25，内摩擦角为35°，内聚力为2MPa。土壤的密度为1800kg/m³，弹性模量为5GPa，泊松比为0.3，内摩擦角为25°，内聚力为1MPa。回填材料根据实际使用的尾砂和废石混合比例，确定其密度为2200kg/m³，弹性模量为15GPa，泊松比为0.28，内摩擦角为30°，内聚力为1.5MPa。在UDEC软件中，通过“Model”菜单中的相应选项，将这些材料参数准确地赋予模型中的不同区域，确保模型能够准确模拟实际材料的力学行为。在模型中，还需要考虑节理、裂隙等地质结构对回填体移动的影响。通过现场地质勘查，获取了节理、裂隙的分布方向、间距和连通性等信息。在UDEC软件中，使用节理单元来模拟这些地质结构，设置节理的法向刚度、切向刚度、内摩擦角和内聚力等参数。根据勘查结果，节理的法向刚度为10GPa/m，切向刚度为5GPa/m，内摩擦角为20°，内聚力为0.5MPa。通过合理设置这些参数，使模型能够准确反映节理、裂隙对回填体移动的影响，提高模拟结果的准确性。5.2模拟结果与分析通过对建立的UDEC数值模型进行模拟计算，得到了金山店铁矿塌陷区回填体在不同开采阶段的移动规律以及应力应变分布情况。从回填体的位移云图可以看出，竖向位移呈现出中间大、四周小的分布特征。在塌陷区中心位置，竖向位移最大，最大值达到了250mm。这是因为塌陷区中心处受到的上覆岩层压力最大，且回填体在该区域的支撑作用相对较弱。随着距离塌陷区中心距离的增加，竖向位移逐渐减小，在塌陷区边缘，竖向位移减小至50mm左右。水平位移则主要集中在塌陷区边缘，最大值为80mm，方向大致与矿体走向一致。这是由于在开采过程中，矿体的采出导致周边岩体向采空区方向移动，从而带动回填体产生水平位移，而塌陷区边缘受到的这种影响更为明显。回填体的应力应变分布也呈现出一定的规律。在应力方面，最大主应力主要分布在回填体底部，最大值为12MPa。这是因为底部承受着上覆岩层的全部压力，处于三向受压状态，导致应力集中。最小主应力则主要分布在回填体顶部，最大值为3MPa。在应变方面，回填体的剪切应变主要集中在与采空区相邻的区域，最大值达到了0.005。这表明在该区域，回填体受到的剪切力较大，容易发生剪切破坏。拉伸应变则主要分布在回填体的边缘和顶部，最大值为0.002。这是由于边缘和顶部受到的约束较小，在开采扰动下容易产生拉伸变形。为了进一步分析回填体移动与各影响因素之间的关系，通过改变模型中的开采深度、开采速度以及降雨强度等参数，进行了多组对比模拟。当开采深度增加50m时，回填体的竖向位移增加了30-40mm，水平位移增加了10-15mm，这表明开采深度的增加会显著增大回填体的位移。当开采速度提高20%时，回填体的竖向位移增加了15-20mm，水平位移增加了8-10mm，说明开采速度的加快也会对回填体的移动产生较大影响。当降雨强度增大50%时，回填体的竖向位移增加了20-30mm，水平位移增加了10-15mm，表明降雨强度的增大同样会加剧回填体的移动。通过这些对比模拟，量化了各影响因素对回填体移动的影响程度，为后续制定回填措施提供了更准确的依据。5.3实际案例验证为验证数值模拟结果的准确性，选取金山店铁矿塌陷区的一个典型区域作为实际案例进行对比分析。该区域在2023年进行了回填作业，回填材料为尾砂和废石的混合物，回填工艺采用分层回填和泵送方式。在实际监测过程中，在该区域设置了10个监测点，利用高精度的位移监测设备对回填体的竖向位移和水平位移进行了为期1年的监测。监测结果显示，回填体的竖向位移在最初3个月内增长较快，平均每月增长12mm，随后增长速率逐渐减缓，6个月后平均每月增长4mm。这与数值模拟结果中竖向位移在回填初期增长较快，随后逐渐减缓的趋势一致。在水平位移方面，实际监测结果表明，水平位移在回填后的前2个月增长较为缓慢，平均每月增长3mm，之后随着开采活动的进行，增长速率加快，4个月后平均每月增长6mm，这也与数值模拟结果相符。通过对比数值模拟结果和实际监测数据，发现两者在位移变化趋势和数值大小上具有较高的一致性。在竖向位移方面，数值模拟预测的位移值与实际监测值的误差在10%以内；在水平位移方面，误