基于数值模拟探究新城金矿Ⅺ#矿体回采对主竖井稳定性的影响

基于数值模拟探究新城金矿Ⅺ#矿体回采对主竖井稳定性的影响一、引言1.1研究背景与意义黄金作为一种重要的战略资源和贵金属，在全球经济、金融和工业领域中具有举足轻重的地位。随着经济的持续发展和科技的不断进步，对黄金的需求呈现出稳步增长的态势。在这一背景下，金矿开采行业面临着巨大的机遇与挑战，如何高效、安全地开采金矿成为了业内关注的焦点。新城金矿作为我国重要的黄金生产基地之一，经过多年的大规模开采，浅部资源逐渐枯竭，目前开采重点已逐步向深部转移。Ⅺ#矿体作为新城金矿深部的重要矿体之一，其开采工作对于维持矿山的可持续发展和保障黄金供应具有关键意义。然而，随着开采深度的增加，地应力环境愈发复杂，开采过程中对周边岩体的扰动也更为显著，这给主竖井的稳定性带来了严峻考验。主竖井作为矿山生产的关键通道，承担着矿石提升、人员运输、材料输送以及通风等重要任务，其稳定性直接关系到矿山生产的安全与效率。一旦主竖井因Ⅺ#矿体回采而出现失稳现象，不仅会导致矿山生产中断，造成巨大的经济损失，还可能引发严重的安全事故，威胁到作业人员的生命安全。在深部开采条件下，传统的竖井保护方法和稳定性评估体系已难以满足实际需求。一方面，深部地应力场的复杂性使得竖井周围的应力分布和变形规律难以准确预测；另一方面，充填法回采等新型采矿工艺的应用，也给竖井稳定性分析带来了新的挑战。因此，深入研究新城金矿Ⅺ#矿体回采对主竖井稳定性的影响，对于保障矿山深部开采的安全与高效具有重要的现实意义。通过开展相关研究，可以为矿山制定科学合理的开采方案和竖井保护措施提供有力的理论依据和技术支持，有效降低开采过程中主竖井失稳的风险，确保矿山生产的顺利进行。同时，本研究成果也将为其他类似矿山在深部开采过程中解决竖井稳定性问题提供有益的参考和借鉴，推动整个金矿开采行业的技术进步和可持续发展。1.2国内外研究现状在矿体回采对竖井稳定性影响的研究领域，国内外学者已开展了大量富有成效的工作，并取得了一系列重要成果。国外方面，早期研究主要侧重于理论分析和经验总结。例如，[国外学者1]通过对多个矿山案例的分析，提出了基于经验公式的竖井稳定性初步评估方法，为后续研究奠定了一定基础。随着计算机技术和数值模拟方法的发展，[国外学者2]运用有限元软件对矿体开采过程进行模拟，分析了不同开采顺序下竖井周围的应力应变分布规律，揭示了开采活动对竖井稳定性的影响机制。此外，[国外学者3]采用现场监测技术，对某矿山竖井在矿体回采期间的变形和位移进行实时监测，获取了大量宝贵的实测数据，进一步验证和完善了理论与数值模拟结果。在实验研究方面，[国外学者4]开展了相似材料模拟实验，模拟矿体开采过程，直观地观察竖井在不同开采条件下的变形破坏特征，为深入理解竖井稳定性问题提供了重要的实验依据。国内学者在该领域也进行了广泛而深入的研究。[国内学者1]基于岩石力学理论，建立了考虑地应力、岩体力学性质等因素的竖井稳定性分析力学模型，通过理论推导得出了竖井围岩应力和位移的解析解，为竖井稳定性分析提供了理论支持。在数值模拟方面，[国内学者2]利用FLAC3D、ANSYS等软件，对不同类型矿体回采过程中竖井的稳定性进行了模拟研究，分析了多种因素如矿体厚度、开采深度、采空区形状等对竖井稳定性的影响程度，并提出了相应的控制措施。[国内学者3]则结合现场监测数据，采用反分析方法确定岩体力学参数，提高了数值模拟结果的准确性。在工程实践方面，[国内学者4]针对某具体矿山，通过优化开采方案和采取加固措施，有效保障了竖井在矿体回采过程中的稳定性，为类似矿山工程提供了实践经验。然而，当前研究仍存在一些不足之处。一方面，在深部复杂地质条件下，地应力场的精确测量和模拟仍然存在困难，导致对竖井稳定性的分析不够准确。深部岩体具有高应力、高渗透压和复杂的地质构造等特点，现有研究方法难以全面考虑这些因素的综合影响。另一方面，虽然已有多种方法用于竖井稳定性分析，但缺乏统一的、综合考虑多因素的评价体系，难以对竖井稳定性进行全面、准确的评估。不同的研究方法和评价指标之间缺乏有效的整合和对比，使得在实际工程应用中难以选择合适的方法。此外，对于充填法回采等新型采矿工艺对竖井稳定性的影响研究还不够深入，尤其是充填体与围岩的相互作用机制以及充填体长期稳定性对竖井的影响等方面，仍有待进一步探索。新城金矿的矿体赋存条件和开采技术条件具有独特性。其Ⅺ#矿体处于深部，地应力复杂，且采用充填法回采，这与其他矿山的情况存在差异。现有的研究成果难以直接应用于新城金矿主竖井稳定性分析。因此，针对新城金矿的具体情况开展深入研究，具有重要的理论和实践意义，有望填补该领域在特定条件下的研究空白，为矿山安全生产提供有力的技术支持。1.3研究内容与方法本研究旨在通过数值模拟深入剖析新城金矿Ⅺ#矿体回采对主竖井稳定性的影响，为矿山安全高效开采提供科学依据和技术支持。具体研究内容和采用的方法如下：研究内容地质模型构建：全面收集新城金矿Ⅺ#矿体及主竖井区域的地质资料，包括地层分布、岩性特征、地质构造（断层、节理等）以及地应力场数据。利用专业的三维地质建模软件，如3DMine、Surpac等，精确构建该区域的三维地质模型，真实反映矿体和主竖井的空间位置关系、矿体的形态产状以及围岩的地质条件，为后续的数值模拟分析奠定坚实基础。数值模型建立：基于所构建的地质模型，选用合适的数值模拟软件，如FLAC3D、ANSYS等，建立新城金矿Ⅺ#矿体回采及主竖井稳定性分析的数值模型。合理确定模型的边界条件，模型底部边界施加固定约束，限制水平和垂直方向的位移；四周侧面边界施加法向约束，仅允许垂直于边界方向的位移；顶部边界为自由边界，模拟实际的地表情况。同时，依据现场岩石力学试验数据和相关经验参数，准确赋予模型中各岩体和充填体合理的力学参数，包括弹性模量、泊松比、内聚力、内摩擦角等，确保数值模型能够准确模拟实际的力学行为。回采过程模拟：根据新城金矿实际采用的充填法回采工艺，在数值模型中按照既定的开采顺序和步骤，逐步模拟Ⅺ#矿体的回采过程。详细考虑每一步回采后采空区的形成、充填体的构筑以及它们对周围岩体应力分布和变形的影响。在模拟过程中，密切关注主竖井周围岩体的力学响应，包括应力变化、位移发展以及塑性区的扩展情况。稳定性分析评估：通过对数值模拟结果的深入分析，全面评估Ⅺ#矿体回采过程中主竖井的稳定性。重点分析主竖井井壁的应力分布特征，确定最大主应力、最小主应力以及剪应力的大小和位置，判断井壁是否存在应力集中现象以及超过岩体强度极限的区域；精确计算主竖井的位移量，包括垂直位移、水平位移和倾斜度，分析位移随回采进程的变化规律，评估位移是否超出允许范围，对竖井的正常使用和安全运行构成威胁；研究主竖井周围塑性区的分布范围和发展趋势，判断塑性区是否会连通形成破坏通道，导致竖井失稳。影响因素分析：系统研究影响主竖井稳定性的各种因素，如矿体倾角、矿体厚度、开采深度、回采顺序、充填体强度等。采用控制变量法，在数值模型中逐一改变各因素的取值，分别模拟不同因素条件下Ⅺ#矿体回采对主竖井稳定性的影响，深入分析各因素对主竖井应力、位移和塑性区分布的影响规律，明确各因素的影响程度和敏感性，为优化开采方案和制定竖井保护措施提供有针对性的依据。研究方法文献研究法：广泛查阅国内外关于矿体回采对竖井稳定性影响的相关文献资料，包括学术论文、研究报告、工程案例等，全面了解该领域的研究现状、发展趋势以及已有的研究成果和方法。通过对文献的深入分析和总结，汲取有益的经验和思路，为本研究提供坚实的理论基础和技术参考，避免重复研究，确保研究的创新性和科学性。现场调研法：深入新城金矿现场，与矿山工程技术人员进行充分交流，详细了解矿山的地质条件、开采现状、采矿工艺以及主竖井的设计参数、运行情况等实际资料。实地观察Ⅺ#矿体的开采现场和主竖井的周边环境，获取第一手资料，为地质模型构建、数值模型参数确定以及研究结果的验证提供真实可靠的数据支持，使研究更贴合工程实际。数值模拟法：运用先进的数值模拟技术，如有限元法（FEM）、有限差分法（FDM）等，对新城金矿Ⅺ#矿体回采过程进行动态模拟分析。利用数值模拟软件强大的计算能力和可视化功能，直观地展示矿体回采过程中主竖井周围岩体的应力、位移和塑性区等力学参数的变化情况，定量分析主竖井的稳定性，预测可能出现的问题，为矿山开采方案的优化和竖井保护措施的制定提供科学依据。理论分析法：基于岩石力学、采矿学等相关学科的基本理论，对数值模拟结果进行深入的理论分析和解释。运用弹性力学、塑性力学理论，分析主竖井周围岩体的应力应变状态和破坏机理；依据采矿工艺原理，探讨回采顺序、充填体作用等因素对主竖井稳定性的影响机制，从理论层面深入理解Ⅺ#矿体回采与主竖井稳定性之间的内在联系，为研究结果提供理论支撑。二、新城金矿及Ⅺ#矿体概况2.1新城金矿简介新城金矿地处胶东半岛渤海湾畔，位于莱州市金城镇境内，其开采历史可追溯至1975年，彼时矿山始建，经过5年的建设，于1980年正式竣工投产，最初核定的采选生产能力为500t/d。自1989年起，新城金矿积极开展历次改扩建工程建设，不断提升自身的生产规模和能力，其核定采选生产能力逐步提升至1250t/d。随着生产的持续发展以及对资源开发利用的不断深入，新的选矿扩能扩产改造工程也在加紧建设中，预计投产后产能将实现大幅飞跃，提高至12000t/d，届时将进一步巩固其在黄金矿业领域的重要地位。作为一座具有采选综合生产能力的国家大型黄金矿山，新城金矿在黄金生产领域成绩斐然。2003年，以新城金矿为主体的“山东黄金”股票在上交所成功上市，这一里程碑事件为企业带来了更为广阔的发展空间和更充足的资金支持，使其能够在技术创新、设备更新、资源勘探等方面加大投入，推动企业持续快速发展。经过多年的发展壮大，截至目前，新城金矿在册员工达2205人，人力资源储备丰富，涵盖了采矿、选矿、地质、机械、电气等多个专业领域，为矿山的高效运营和可持续发展提供了坚实的人才保障。截止2022年2月10日，新城金矿累计自产金产量突破100吨大关，这一辉煌成就不仅彰显了其在黄金生产方面的强大实力，也为国家的黄金储备和经济发展做出了重要贡献。同时，累计实现利润总额72.66亿元，展现出良好的经济效益和盈利能力。新城金矿在发展过程中，始终高度重视环境保护和可持续发展，致力于打造绿色矿山。其是全国矿山企业中唯一的“国家环境友好企业”，并跻身“中国黄金生产十大矿山”“中国黄金经济效益十佳矿山”行列，成为我国黄金行业的佼佼者。此外，矿山还获得了“全国黄金行业先进集体”“全国黄金行业科技先进集体”“首届全国资源合理利用先进矿山企业”、首批“国家级绿色矿山”“全国矿产资源节约与综合利用先进适用技术推广应用示范矿山”“绿色矿山突出贡献奖”“中国黄金行业最佳绿色环保奖”“山东省安全生产应急管理示范点”“山东省劳动关系和谐企业”“山东省富民兴鲁劳动奖状”等众多荣誉，这些荣誉是对其在安全生产、科技创新、资源综合利用、环境保护等方面所做出努力和取得成绩的高度认可。在开采技术方面，新城金矿采用主竖井与斜坡道联合开拓的方式，这种开拓方式结合了主竖井提升能力大、效率高和斜坡道运输灵活、适应性强的优点，能够满足矿山深部开采和大规模生产的需求。在采矿方法上，主要采用充填采矿法，该方法具有对围岩扰动小、能够有效控制地压、减少地表塌陷、提高矿石回收率等优点，尤其适用于新城金矿深部矿体的开采条件。在选矿工艺上，不断引进和创新先进的技术和设备，提高选矿回收率和精矿质量，降低生产成本和资源浪费。2.2Ⅺ#矿体地质特征Ⅺ#矿体是新城金矿深部开采的重点对象，其赋存条件复杂，对开采过程中的技术和安全管理提出了较高要求。该矿体走向大致呈NE55°，在空间上具有一定的延伸长度，其走向长度对于确定开采范围和规划开采顺序具有重要意义。倾向NW，这种倾向决定了矿体在地下的空间位置与主竖井以及其他巷道的相对关系，对开采过程中的通风、运输等系统的布局产生影响。倾角约为29°，属于缓倾斜矿体。缓倾斜矿体的开采相较于急倾斜矿体，在矿石搬运、顶板管理等方面存在不同的技术难点和要求，需要采取相应的开采工艺和安全措施。矿体厚度变化较大，平均水平厚度约为60m，属于中厚矿体。矿体厚度的变化不仅影响开采方法的选择，还对开采效率和资源回收率有着重要影响。较厚的矿体在开采过程中需要考虑如何合理地进行分层开采，以确保开采的安全性和高效性。在矿化方面，Ⅺ#矿体矿化较好，蚀变以黄铁矿化、绢云母化、硅化为主，次为绿泥石化、碳酸盐化、钾化等。这些矿化和蚀变特征反映了矿体的形成地质环境和过程，同时也对矿石的物理力学性质和选矿工艺产生影响。例如，黄铁矿化可能导致矿石硬度增加，而绢云母化则可能使矿石的韧性发生变化，在开采和选矿过程中需要根据这些特性进行相应的技术调整。矿石工业类型为贫硫银金矿石，这决定了矿石的价值和选矿工艺的重点。贫硫银金矿石的选矿需要注重金和银的回收，同时要考虑如何降低硫等杂质的含量，以提高精矿的质量。上盘围岩主要为构造角砾岩和碎裂花岗岩，岩石破碎，f=4～6，抗剪切强度低，稳定性极差。这种围岩条件增加了开采过程中的支护难度和安全风险，需要采取有效的支护措施来确保开采过程中围岩的稳定性，防止顶板坍塌等事故的发生。下盘围岩主要为黄铁绢英岩化花岗闪长岩，f=10，比上盘围岩抗剪切强度大，稳定性好。相对稳定的下盘围岩为开采提供了一定的有利条件，但在开采过程中仍需关注上下盘围岩的差异，合理设计开采方案和支护措施。矿体内存在一些特殊的地质构造，如矿块上盘的焦家主断裂面有厚10～20cm的断层泥，具有隔水作用，但遇水极易软化，较易冒落，很不稳固。这种特殊的地质构造在开采过程中可能引发涌水、冒顶等安全问题，需要提前进行地质勘查和风险评估，制定相应的应急预案和防范措施。新城金矿通过现场取样和室内试验，对Ⅺ#矿体的矿岩物理力学性质进行了测定。测定结果显示，矿石的密度为[X]kg/m³，容重为[X]kN/m³，这两个参数反映了矿石的质量特性，对于矿石的运输、提升等环节的设计和设备选型具有重要参考价值。矿石的抗压强度为[X]MPa，抗拉强度为[X]MPa，抗剪强度为[X]MPa，这些强度参数是评估矿石开采难度和确定开采工艺的重要依据。在开采过程中，需要根据矿石的强度特性选择合适的采矿设备和爆破参数，以确保开采效率和安全性。岩石的弹性模量为[X]GPa，泊松比为[X]，这两个参数描述了岩石在受力时的变形特性，对于分析开采过程中岩体的应力应变分布、预测岩体的稳定性具有重要意义。内聚力为[X]MPa，内摩擦角为[X]°，它们反映了岩石颗粒之间的粘结力和摩擦力，是判断岩体稳定性和设计支护结构的关键参数。在进行数值模拟分析时，这些物理力学参数将作为重要的输入条件，用于准确模拟矿体回采过程中岩体的力学行为和主竖井的稳定性变化。2.3主竖井工程概述新城金矿主竖井位于矿山的核心区域，其地理位置经过精确的测量和定位确定，处于[具体经纬度]，该位置的选择充分考虑了矿山整体的开采布局、矿体分布以及地面设施的配套等因素，是保障矿山高效生产的关键节点。主竖井的井口标高为[X]m，井底标高为[-X]m，井深达到了[X]m，如此深的井深在深部开采中具有重要意义，它直接关系到矿山深部矿体的开采能力和效率。井筒直径为[X]m，这种直径设计既能满足矿山目前的生产规模和运输需求，又考虑到未来可能的产能提升和设备升级。在支护方式上，主竖井采用了钢筋混凝土支护与锚杆锚索联合支护的方式。在井壁的外层，钢筋混凝土支护形成了一个坚固的外壳，其强度等级达到[X]，厚度为[X]m，能够有效地承受来自周围岩体的压力，抵抗岩体的变形和破坏。在混凝土内部，锚杆和锚索按照一定的间距布置，锚杆长度为[X]m，锚索长度为[X]m，它们深入到岩体内部，将井壁与周围的岩体紧密地锚固在一起，增强了井壁与岩体之间的连接，提高了井壁的稳定性和承载能力。这种联合支护方式是根据新城金矿深部复杂的地质条件和地应力环境精心设计的，经过实践验证，能够有效地保障主竖井在长期开采过程中的安全稳定运行。主竖井在矿山生产中扮演着无可替代的重要角色，堪称矿山的“生命线”。在矿石提升方面，它承担着将井下开采出来的矿石源源不断地提升至地面的重任，目前其提升能力达到了[X]t/d，随着矿山开采规模的不断扩大，对主竖井的提升能力也提出了更高的要求。通过高效的提升系统，矿石能够快速、安全地到达地面，为后续的选矿和冶炼提供原料保障。在人员运输方面，主竖井是井下作业人员进出的主要通道，每天有大量的工人通过主竖井上下井，保障他们的安全和便捷出行是矿山安全生产的重要环节。主竖井配备了专门的人员提升设备，如罐笼等，严格按照安全规范进行操作，确保人员运输的安全可靠。在材料输送方面，主竖井负责将井下开采所需的各种材料，如炸药、支护材料、设备零部件等，从地面运输到井下各个作业地点。这些材料对于维持井下正常的开采作业至关重要，主竖井的高效运输能够确保材料及时供应，避免因材料短缺而导致的生产延误。在通风方面，主竖井是矿山通风系统的重要组成部分，通过它将新鲜空气引入井下，同时将井下的污浊空气排出地面，为井下作业人员创造良好的工作环境，保障他们的身体健康和安全生产。此外，主竖井还在矿山的排水、供电、通讯等系统中发挥着关键作用，它为这些系统的管线和设备提供了安装和运行的空间，确保矿山各个系统的正常运行。一旦主竖井出现故障或失稳，将对矿山的生产和安全造成严重影响，可能导致生产中断、人员被困、安全事故等一系列问题，因此，保障主竖井的稳定性对于新城金矿的可持续发展具有至关重要的意义。三、影响竖井稳定性的因素分析3.1矿体回采因素3.1.1回采顺序回采顺序是影响主竖井稳定性的关键因素之一，不同的回采顺序会导致主竖井周边应力分布和变形呈现出显著差异。在新城金矿Ⅺ#矿体的开采中，常见的回采顺序有自上而下、自下而上以及间隔开采等方式。当采用自上而下的回采顺序时，随着上部矿体的逐渐采出，上覆岩层的重量会逐渐转移到下部未开采的矿体和周围岩体上，导致下部岩体所承受的压力增大。在主竖井周边，这种压力的变化会引起应力重新分布，使得主竖井井壁所受的垂直应力和水平应力增加。由于上部开采形成的采空区对下部岩体的约束减弱，下部岩体在高应力作用下更容易发生变形，进而导致主竖井产生位移和变形。若上部采空区的充填不及时或充填质量不佳，无法有效支撑上覆岩层，会进一步加剧下部岩体的变形，增加主竖井失稳的风险。与之相反，自下而上的回采顺序会使下部矿体开采后形成的采空区上方岩体的稳定性受到影响。在这种情况下，上部未开采矿体和岩体的重量会对下部采空区周围的岩体产生较大的压力，使得主竖井下部井壁承受较大的应力。随着开采向上推进，下部采空区的扩大和上部岩体的移动，主竖井所受的不均匀应力会逐渐增大，可能导致井壁出现开裂、剥落等破坏现象。由于下部采空区的存在，上部岩体在开采过程中的变形会更加复杂，这也增加了主竖井变形和失稳的可能性。间隔开采是指在矿体开采过程中，按照一定的间隔顺序进行开采，先开采部分矿体，留下部分矿体作为临时支撑。这种回采顺序在一定程度上可以缓解开采过程中的应力集中现象，减少对主竖井的影响。在间隔开采过程中，留下的矿体可以分担部分上覆岩层的压力，使得主竖井周边的应力分布相对均匀。然而，若间隔开采的参数设置不合理，如间隔距离过大或过小，会导致留下的矿体无法有效支撑上覆岩层，或者开采过程中频繁地改变应力状态，反而会增加主竖井的受力复杂性，对其稳定性产生不利影响。3.1.2开采深度随着开采深度的不断增加，地压增大是一个不可避免的趋势，这对主竖井的稳定性产生了多方面的影响。在深部开采环境下，地应力主要由上覆岩层的重量和地质构造应力组成。随着深度的增加，上覆岩层的厚度不断增大，其重量所产生的垂直应力也随之增大。深部地质构造复杂，构造应力在岩体中分布不均，也会对主竖井周围的应力状态产生显著影响。主竖井在深部开采过程中，井壁所承受的压力大幅增加。由于井壁需要承受来自周围岩体的巨大压力，其内部的应力分布变得更加复杂。在高应力作用下，井壁容易出现应力集中现象，尤其是在井壁的拐角、接口等部位，应力集中程度更为明显。当井壁所承受的应力超过其自身的强度极限时，井壁就会发生破裂、剥落等破坏现象，严重影响主竖井的稳定性。深部岩体在高应力作用下，其力学性质会发生变化，岩体的脆性增强，塑性降低。这意味着岩体在受到外力作用时更容易发生脆性破坏，而难以通过塑性变形来释放能量。主竖井周围的岩体一旦发生脆性破坏，就会形成破碎带，导致岩体的承载能力下降，进一步威胁主竖井的安全。开采深度的增加还会导致地下水压力增大。在深部地层中，地下水受到上覆岩层的压力作用，其压力随着深度的增加而增大。高水压会对主竖井井壁产生渗透压力，使得井壁的受力情况更加复杂。若井壁的防水措施不到位，地下水还可能渗入井壁内部，对井壁结构造成侵蚀和破坏，降低井壁的强度和稳定性。3.1.3采场结构参数采场的尺寸、形状等结构参数对主竖井稳定性有着重要的影响。采场尺寸方面，采场的长度、宽度和高度都会影响其周围岩体的应力分布和变形情况，进而影响主竖井的稳定性。较大尺寸的采场在开采过程中会形成更大范围的采空区，导致上覆岩层的移动和变形范围增大。在主竖井附近，这种大范围的岩体移动会对主竖井产生较大的附加应力，增加主竖井的变形和失稳风险。若采场高度过高，上覆岩层在自重作用下的下沉量会增大，对主竖井周围岩体的扰动也会更加剧烈，可能导致主竖井周围岩体的松动和破坏。采场形状也对主竖井稳定性有显著影响。不规则形状的采场，如多边形、异形等，会使采场周围岩体的应力分布不均匀，容易产生应力集中现象。在主竖井附近，这种应力集中会对主竖井井壁产生局部高压力，导致井壁出现裂缝、破损等问题。相比之下，规则形状的采场，如矩形、圆形等，其周围岩体的应力分布相对均匀，对主竖井的影响相对较小。然而，即使是规则形状的采场，若其尺寸与主竖井的距离不合理，也会对主竖井稳定性产生不利影响。采场的跨度也是一个关键参数。跨度较大的采场，其顶板在开采后更容易发生垮落，从而引发周围岩体的连锁反应，对主竖井产生冲击和振动。这种冲击和振动会使主竖井井壁受到动态荷载的作用，增加井壁的疲劳损伤，降低其稳定性。为了减小采场结构参数对主竖井稳定性的影响，在设计采场时，需要综合考虑矿体的赋存条件、主竖井的位置以及岩体的力学性质等因素，合理确定采场的尺寸、形状和跨度，以确保主竖井在矿体回采过程中的安全稳定。3.2地质条件因素3.2.1地应力地应力是影响主竖井围岩稳定性的关键地质因素之一，其大小和方向对主竖井在Ⅺ#矿体回采过程中的稳定性有着重要影响。地应力主要由上覆岩层的自重应力、构造应力以及岩体的残余应力等组成。在新城金矿深部开采区域，地应力场呈现出复杂的特征，其大小和方向受到多种地质因素的控制。随着开采深度的增加，上覆岩层的重量不断增大，导致自重应力显著增加。在新城金矿Ⅺ#矿体所处的深部区域，自重应力可能达到较高的数值，这使得主竖井周围的岩体承受着巨大的压力。由于深部地质构造的复杂性，构造应力在岩体中分布不均，其方向和大小难以准确预测。构造应力可能会与自重应力相互叠加或抵消，进一步加剧主竖井周围岩体应力状态的复杂性。地应力的大小和方向直接影响主竖井围岩的应力分布和变形特征。当主竖井的轴线方向与最大主应力方向夹角较大时，井壁所受的切向应力会显著增大，容易导致井壁出现破裂、剥落等破坏现象。在高应力作用下，主竖井周围的岩体可能会发生塑性变形，形成塑性区。塑性区的扩展会降低岩体的承载能力，增加主竖井失稳的风险。若地应力大小超过岩体的强度极限，还可能引发岩爆等动力灾害，对主竖井的稳定性造成严重威胁。3.2.2岩体完整性岩体的完整性是影响主竖井稳定性的另一个重要地质条件因素，主要取决于岩体中节理、裂隙等结构面的发育程度和分布特征。在新城金矿Ⅺ#矿体及主竖井周围的岩体中，节理、裂隙较为发育，这些结构面的存在破坏了岩体的连续性和完整性，使其力学性质发生显著变化。节理和裂隙的存在使得岩体的强度明显降低。由于结构面的抗剪强度通常低于完整岩体的强度，当岩体受到外力作用时，容易沿着节理和裂隙发生滑动和破坏。在主竖井周围，由于开采活动引起的应力变化，岩体中的节理和裂隙可能会进一步张开、扩展，导致岩体的完整性进一步恶化，从而降低主竖井周围岩体的承载能力。节理和裂隙的分布特征也会影响主竖井的稳定性。如果节理和裂隙在主竖井周围呈密集分布，且相互连通，会形成潜在的滑动面，增加主竖井失稳的可能性。在这种情况下，一旦岩体受到较大的外力作用，就可能沿着这些滑动面发生整体滑动，导致主竖井井壁破裂、坍塌。此外，节理和裂隙还会影响岩体的渗透性。在新城金矿深部开采过程中，地下水可能会沿着节理和裂隙流动，进一步弱化岩体的力学性质。地下水的存在会降低岩体的有效应力，增加岩体的孔隙水压力，使得岩体更容易发生变形和破坏。3.2.3地下水地下水在新城金矿深部开采区域广泛存在，其对岩体力学性质的弱化作用以及对主竖井稳定性的影响不容忽视。地下水的存在会改变岩体的物理和力学性质。一方面，地下水会使岩体发生饱水，导致岩体的重度增加，从而增大岩体的自重应力。另一方面，地下水会对岩体中的矿物成分产生溶解和侵蚀作用，破坏矿物颗粒之间的胶结结构，降低岩体的强度。在主竖井周围，地下水的渗透会产生渗透压力。当渗透压力较大时，会对主竖井井壁产生向外的推力，增加井壁的受力负担。若井壁的强度不足以抵抗这种推力，就会发生变形和破坏。地下水还可能导致岩体的软化和泥化，进一步降低岩体的稳定性。在新城金矿的地质条件下，部分岩体中含有遇水易软化的矿物成分，如蒙脱石等。当这些岩体与地下水接触后，会发生软化和泥化现象，使得岩体的抗剪强度大幅降低，容易引发主竖井周围岩体的坍塌和滑移。此外，地下水还可能对主竖井的支护结构产生腐蚀作用。主竖井采用钢筋混凝土支护与锚杆锚索联合支护的方式，地下水的长期侵蚀会导致钢筋生锈、混凝土强度降低，从而削弱支护结构的承载能力，影响主竖井的稳定性。3.3工程施工因素3.3.1爆破振动在新城金矿Ⅺ#矿体的回采过程中，爆破是主要的破岩手段，然而，爆破产生的振动波会对主竖井井壁结构造成显著的破坏作用。爆破振动是一种复杂的瞬态振动现象，其产生的振动波以弹性波的形式在岩体中传播，当这些振动波传播到主竖井井壁时，会引起井壁的振动和应力变化。爆破振动对主竖井井壁的破坏作用主要体现在以下几个方面：当爆破振动波的频率与主竖井井壁的固有频率接近或相等时，会发生共振现象。共振会使井壁的振动幅度急剧增大，导致井壁所承受的应力远远超过其设计强度，从而引发井壁的破裂和损坏。在共振状态下，井壁可能会出现明显的裂缝，甚至部分结构脱落，严重影响主竖井的稳定性。爆破振动会在井壁上产生拉应力和剪应力。由于井壁材料的抗拉强度和抗剪强度相对较低，当这些应力超过井壁材料的极限强度时，井壁就会出现拉伸破坏和剪切破坏。拉伸破坏通常表现为井壁表面出现垂直于井壁的裂缝，而剪切破坏则表现为井壁出现斜向裂缝或块状剥落。长期的爆破振动作用会使主竖井井壁材料产生疲劳损伤。随着爆破次数的增加，井壁材料内部的微裂纹逐渐扩展、连通，导致材料的强度和刚度不断降低，最终可能引发井壁的突发性破坏。这种疲劳破坏在初期可能不易察觉，但随着时间的推移，其危害会逐渐显现出来。3.3.2支护措施合理有效的支护措施是增强主竖井稳定性的关键，对于保障主竖井在Ⅺ#矿体回采过程中的安全运行具有重要意义。主竖井采用的钢筋混凝土支护与锚杆锚索联合支护方式，在一定程度上能够提高主竖井的稳定性。钢筋混凝土支护通过其自身的强度和刚度，能够承受来自周围岩体的压力，限制井壁的变形。钢筋在混凝土中起到增强抗拉强度的作用，使井壁能够更好地抵抗拉伸应力。锚杆锚索则深入到周围岩体中，将井壁与岩体紧密地锚固在一起，增强了井壁与岩体之间的连接，提高了岩体的整体稳定性。锚杆锚索能够有效地传递岩体的应力，将不稳定的岩体与稳定的岩体连接起来，形成一个共同承载的体系。为了进一步增强主竖井的稳定性，还可以采取以下支护措施：在主竖井井壁周围设置钢支撑，如工字钢、槽钢等。钢支撑具有较高的强度和刚度，能够快速有效地支撑井壁，抵抗爆破振动和采动压力的作用。钢支撑可以与钢筋混凝土支护和锚杆锚索联合使用，形成一个更加坚固的支护体系。采用喷射混凝土支护，在井壁表面喷射一层混凝土，能够及时封闭井壁表面的裂缝和孔隙，防止岩体风化和地下水的侵蚀，同时也能够增加井壁的强度和稳定性。喷射混凝土还可以与钢筋网联合使用，进一步提高支护效果。在主竖井周围的岩体中进行注浆加固，通过将水泥浆或化学浆液注入岩体的裂缝和孔隙中，填充岩体的空隙，增强岩体的整体性和强度，从而提高主竖井周围岩体的承载能力和稳定性。四、数值模拟原理与模型建立4.1数值模拟方法选择数值模拟方法在岩土工程领域中应用广泛，常见的有有限元法（FEM）、有限差分法（FDM）、边界元法（BEM）等，每种方法都有其独特的原理、适用范围和优缺点。有限元法基于变分原理，将求解区域离散为有限个单元，通过对单元进行分析和组装，得到整个区域的近似解。它的优点是对复杂几何形状和边界条件具有很强的适应性，能够处理多种物理场的耦合问题，在结构力学、热传导等领域应用广泛。在处理大变形问题时，有限元法需要不断更新网格，计算过程较为复杂，计算效率相对较低。有限差分法是将连续的求解区域离散为网格，通过差分近似导数，将控制方程转化为代数方程组进行求解。它的求解过程直观，计算效率较高，尤其适用于求解波动问题和大变形问题。然而，有限差分法对不规则区域的处理能力相对较弱，边界条件的施加也相对复杂。边界元法以边界积分方程为基础，将求解区域的维数降低一维，只需对边界进行离散，从而减少了计算量。它在无限域问题和应力集中问题的求解上具有明显优势。但边界元法的应用依赖于基本解的选取，对于复杂的材料本构关系和物理过程，基本解的获取较为困难，适用范围相对较窄。在新城金矿Ⅺ#矿体回采对主竖井稳定性影响的研究中，涉及到复杂的地质条件、矿体开采过程中的大变形以及岩体的非线性力学行为等问题。经过综合考虑，本研究选择有限差分软件FLAC3D进行数值模拟。FLAC3D是由美国ItascaConsultingGroup开发的一款专门用于岩土工程分析的软件，其核心算法基于快速拉格朗日分析方法，能够很好地模拟岩土材料在大变形条件下的力学行为。FLAC3D在处理大变形问题时具有独特的优势。它采用拉格朗日差分法，允许网格随着材料的变形而移动和变形，能够真实地反映矿体开采过程中岩体的大变形特征，这对于准确分析主竖井在矿体回采过程中的稳定性至关重要。该软件内置了丰富的材料本构模型，包括弹性模型、弹塑性模型（如莫尔-库伦模型、德鲁克-普拉格模型等）以及考虑岩体节理裂隙的节理模型等，能够根据新城金矿的实际地质条件和岩体力学性质，准确地描述岩体的力学行为。FLAC3D具有强大的后处理功能，能够直观地展示模拟结果，如应力分布云图、位移矢量图、塑性区分布图等，便于对主竖井的稳定性进行分析和评估。4.2FLAC3D软件简介FLAC3D（FastLagrangianAnalysisofContinuain3Dimensions）是一款由美国ItascaConsultingGroup开发的专业三维岩土工程数值模拟软件，基于快速拉格朗日分析方法，在岩土工程领域应用广泛。其基本原理蕴含着深厚的力学和数学基础。在力学原理方面，它基于牛顿第二定律，通过建立运动方程来描述介质的力学行为。在计算过程中，将求解区域离散为一系列的单元，这些单元在受力作用下会产生变形和运动，FLAC3D通过追踪单元的运动和变形来模拟整个介质的响应。从数学原理角度，FLAC3D采用显式有限差分法对控制方程进行求解。这种方法将时间和空间进行离散化，把连续的物理问题转化为一系列离散的代数方程进行求解。在每一个时间步长内，根据单元的受力状态和前一时刻的状态，计算出单元的加速度、速度和位移，进而更新单元的位置和应力状态。在处理大变形问题时，采用拉格朗日算法，允许网格随着材料的变形而移动和变形。这种算法能够真实地反映材料在大变形过程中的力学行为，避免了传统欧拉算法在处理大变形时可能出现的网格畸变问题。FLAC3D具备诸多显著的功能特点。它拥有丰富的材料本构模型库，包含弹性模型、弹塑性模型（如莫尔-库伦模型、德鲁克-普拉格模型等）以及针对岩体节理裂隙的节理模型等，能够精准地模拟各种岩土材料的力学行为。在模拟新城金矿Ⅺ#矿体回采过程时，可以根据矿体和围岩的实际力学性质选择合适的本构模型，从而准确地分析岩体的变形和破坏过程。提供了多种类型的结构单元，如梁单元、壳单元、锚索单元等，能够模拟各种地下工程结构与岩土体之间的相互作用。在分析主竖井的稳定性时，可以利用这些结构单元来模拟主竖井的支护结构，研究支护结构与周围岩体的协同工作机制。具有强大的后处理功能，能够直观地展示模拟结果。通过绘制应力分布云图、位移矢量图、塑性区分布图等，可以清晰地了解主竖井周围岩体在矿体回采过程中的力学响应，为分析主竖井的稳定性提供直观依据。支持流-固耦合分析，能够考虑地下水渗流对岩体力学行为的影响。在新城金矿深部开采中，地下水对岩体稳定性的影响不容忽视，FLAC3D的流-固耦合功能可以更真实地模拟实际工程情况。在岩土工程领域，FLAC3D展现出独特的应用优势。它能够有效处理复杂地质条件下的工程问题，对于新城金矿这种地质构造复杂、地应力条件特殊的矿山，FLAC3D可以通过准确模拟地质条件和岩体力学行为，为工程设计和施工提供可靠的参考。在模拟大变形问题上具有明显优势，矿体回采过程中岩体的变形往往较大，FLAC3D的拉格朗日算法和显式有限差分法能够很好地适应这种大变形情况，保证模拟结果的准确性。计算效率较高，采用显式有限差分法，不需要求解大型的线性方程组，计算过程相对简单，能够在较短的时间内完成复杂模型的计算，提高了研究工作的效率。具有良好的开放性和可扩展性，用户可以通过自定义FISH语言编写自己的本构模型、边界条件和求解算法等，以满足特殊工程问题的需求。这使得FLAC3D能够不断适应岩土工程领域的发展和创新，为解决各种复杂工程问题提供有力支持。4.3数值模型建立4.3.1模型范围确定模型范围的确定需综合考虑多种因素，以确保数值模拟结果的准确性和可靠性。从理论层面而言，模型边界应设置在远离开采区域，使得边界条件的设定对开采区域的计算结果影响可忽略不计。在实际操作中，依据相关研究和工程经验，一般认为模型边界距开采区域的距离应达到一定倍数的开采影响范围。对于新城金矿Ⅺ#矿体回采对主竖井稳定性影响的模拟，考虑到矿体的规模、开采深度以及地应力传播范围等因素，确定模型在x方向（走向）上的长度为[X]m，涵盖了矿体及其周围一定范围的岩体，足以反映矿体开采在走向方向上对岩体的影响。在y方向（倾向）上的长度为[X]m，保证了能够捕捉到矿体倾向方向上的应力应变变化。在z方向（垂直方向）上，模型上表面取至地表，下表面位于矿体底部以下[X]m处，以充分考虑深部岩体对开采的响应。在确定模型范围时，还需考虑到计算资源的限制。若模型范围过大，会导致计算量急剧增加，计算时间大幅延长，甚至可能超出计算机的处理能力。因此，在保证模拟精度的前提下，需要对模型范围进行优化，以提高计算效率。通过多次试算和分析，上述确定的模型范围既能满足模拟精度要求，又能在合理的计算时间内完成模拟任务。4.3.2模型网格划分模型网格划分遵循一定的原则和方法，以保证计算结果的准确性和计算效率。在原则方面，首先要保证网格的质量，避免出现畸形网格，确保网格单元的形状规则，各边和各角的角度合理，以提高计算的稳定性和精度。在矿体和主竖井等重点关注区域，采用加密网格的方式，使网格尺寸较小，能够更精确地捕捉这些区域的应力应变变化。而在远离开采区域和主竖井的部位，网格尺寸可以适当增大，以减少计算量，提高计算效率。在方法上，使用FLAC3D软件自带的网格划分功能，采用六面体单元对模型进行离散化。这种单元类型在模拟岩土体的力学行为时具有较好的适应性，能够准确地模拟岩体的变形和破坏过程。在划分网格时，根据模型的几何形状和边界条件，合理设置网格生成参数，如网格尺寸、生长比率等，以确保网格的均匀性和合理性。经过网格划分后，得到的网格模型如图[X]所示。从图中可以清晰地看到，在Ⅺ#矿体和主竖井周围，网格较为密集，能够精细地模拟这些区域的力学响应。而在模型的边缘部分，网格相对稀疏，符合计算效率和精度的要求。整个网格模型共计划分[X]个单元，[X]个节点，为后续的数值模拟计算提供了良好的基础。[此处插入网格模型图]4.3.3材料参数选取材料参数的准确选取对于数值模拟结果的可靠性至关重要，本研究通过岩石力学试验和参考工程经验来确定模型中各材料参数。在岩石力学试验方面，从新城金矿Ⅺ#矿体及周围围岩中采集了大量岩样，在实验室进行了一系列物理力学性质测试。通过单轴压缩试验，获取了岩石的抗压强度、弹性模量和泊松比等参数。在单轴压缩试验中，将岩样加工成标准圆柱体，在压力试验机上以恒定的加载速率进行加载，记录岩样在加载过程中的应力应变曲线，从而计算出抗压强度和弹性模量。泊松比则通过测量岩样在轴向加载时的横向变形和轴向变形来确定。通过直接剪切试验，测定了岩石的内聚力和内摩擦角，试验时将岩样置于剪切盒中，施加垂直荷载和水平剪切力，记录岩样在剪切过程中的破坏荷载和剪切位移，进而计算出内聚力和内摩擦角。对于充填体的参数，参考类似矿山的工程经验和相关研究成果。考虑到充填体的材料组成、配比以及养护条件等因素，确定其弹性模量为[X]GPa，泊松比为[X]，内聚力为[X]MPa，内摩擦角为[X]°。这些参数能够较好地反映新城金矿充填体的力学特性。主竖井的钢筋混凝土支护参数根据设计资料和实际施工情况确定，混凝土的弹性模量为[X]GPa，泊松比为[X]，抗压强度为[X]MPa。钢筋的弹性模量为[X]GPa，屈服强度为[X]MPa。锚杆和锚索的参数根据其规格和材料特性确定，锚杆的弹性模量为[X]GPa，抗拉强度为[X]MPa，锚索的弹性模量为[X]GPa，抗拉强度为[X]MPa。最终确定的各材料参数如表1所示：[此处插入材料参数表]4.3.4边界条件设定模型边界条件的设定直接影响模拟结果的准确性，依据实际工程情况和数值模拟理论进行设定。在模型底部边界，由于岩体处于深部，位移受到极大限制，因此施加固定约束，即限制x、y、z三个方向的位移，使底部边界在模拟过程中不会发生移动。在模型四周侧面边界，为了模拟岩体在水平方向上的受力和变形情况，施加法向约束，仅允许垂直于边界方向的位移，限制平行于边界方向的位移。顶部边界为自由边界，模拟实际的地表情况，不受任何外力约束，允许岩体在垂直方向上自由变形。在开采过程中，矿体回采会导致岩体应力重新分布，在模型中通过“单元生死”功能来模拟矿体的开采和充填过程。当开采某一区域的矿体时，将该区域对应的单元“杀死”，使其不再承担力学作用，从而模拟采空区的形成。在充填过程中，将“杀死”的单元重新“激活”，并赋予充填体的材料参数，以模拟充填体对采空区的支撑作用。通过这样的边界条件设定和开采充填模拟方式，能够较为真实地反映新城金矿Ⅺ#矿体回采过程中主竖井周围岩体的力学响应和稳定性变化。4.4模拟方案设计为全面深入探究新城金矿Ⅺ#矿体回采对主竖井稳定性的影响，精心设计了多组模拟方案，每组方案都涵盖了不同的回采顺序、开采深度以及采场结构参数，旨在通过对比分析，明确各因素对主竖井稳定性的具体影响规律。4.4.1回采顺序方案设计了三种典型的回采顺序方案。方案一为自上而下回采，具体实施步骤为：从矿体的最上部开始，按照一定的分层厚度，逐层向下进行开采。在每一层开采完成后，及时对采空区进行充填处理，以确保上覆岩层的稳定性。方案二为自下而上回采，即从矿体的最下部开始，逆向向上进行开采。在这个过程中，同样需要关注采空区的充填质量和时间，以避免下部采空区对上部开采的影响。方案三为间隔开采，将矿体划分为多个采区，按照一定的间隔顺序进行开采。先开采部分采区，留下部分采区作为临时支撑，待已开采采区的充填体达到一定强度后，再对留下的采区进行开采。通过对比这三种回采顺序方案下主竖井的应力、位移和塑性区分布情况，分析不同回采顺序对主竖井稳定性的影响。4.4.2开采深度方案考虑到开采深度对主竖井稳定性的显著影响，设计了三个不同开采深度的方案。方案一模拟开采深度为[X]m，此深度接近当前新城金矿的开采现状，能够反映当前开采条件下主竖井的稳定性情况。方案二模拟开采深度为[X+100]m，通过增加开采深度，分析地应力增加对主竖井稳定性的影响趋势。方案三模拟开采深度为[X+200]m，进一步加大开采深度，研究在更深部开采条件下主竖井的稳定性变化规律，以及可能出现的新问题和挑战。4.4.3采场结构参数方案针对采场的尺寸和形状对主竖井稳定性的影响，设计了多组采场结构参数方案。在采场尺寸方面，方案一设置采场长度为[X]m，宽度为[Y]m，高度为[Z]m；方案二设置采场长度为[X+20]m，宽度为[Y+10]m，高度为[Z+5]m，通过增大采场尺寸，分析其对主竖井稳定性的影响。在采场形状方面，方案一采用矩形采场，这种形状的采场在开采过程中应力分布相对较为规则；方案二采用梯形采场，其形状会导致采场周围岩体的应力分布不均匀，通过对比两种形状采场下主竖井的稳定性，研究采场形状对主竖井稳定性的影响机制。具体模拟方案设计如表2所示：[此处插入模拟方案设计表]通过以上精心设计的模拟方案，能够全面系统地研究不同因素对新城金矿Ⅺ#矿体回采过程中主竖井稳定性的影响，为后续的数值模拟分析和结果讨论提供丰富的数据和依据。五、数值模拟结果与分析5.1不同回采顺序下的模拟结果在本次数值模拟中，分别对自上而下、自下而上和间隔开采三种回采顺序进行了模拟，得到了主竖井周边应力、位移和塑性区分布云图，通过对这些云图的详细分析，可以清晰地了解不同回采顺序对主竖井稳定性的影响。5.1.1应力分布分析在自上而下的回采顺序下，从应力分布云图（图[X1]）中可以观察到，随着开采的进行，主竖井周边的应力逐渐增大。在矿体上部开采阶段，主竖井顶部附近的应力集中现象较为明显，最大主应力达到[X]MPa，主要是由于上部矿体开采后，上覆岩层的重量转移到主竖井顶部周围的岩体上，导致该区域应力升高。随着开采向下推进，主竖井中部和下部的应力也逐渐增加，在开采至矿体中部时，主竖井中部井壁的最大主应力达到[X+ΔX1]MPa。这是因为上部采空区的不断扩大，使得应力传递范围增大，主竖井中部受到的影响逐渐增强。在整个开采过程中，主竖井井壁的应力分布呈现出不均匀的状态，靠近矿体一侧的井壁应力明显大于远离矿体一侧，这是由于矿体开采对靠近矿体一侧的岩体扰动更大，导致应力集中更为显著。[此处插入自上而下回采顺序应力分布云图]自下而上回采顺序的应力分布云图（图[X2]）显示出与自上而下回采顺序不同的特征。在开采初期，由于下部矿体的开采，主竖井底部周围的岩体应力迅速增大，最大主应力达到[X']MPa，这是因为下部采空区的形成使得上覆岩层的支撑力减弱，主竖井底部承受了较大的压力。随着开采向上推进，主竖井中部和上部的应力也逐渐增加，在开采至矿体中部时，主竖井中部井壁的最大主应力达到[X'+ΔX2]MPa。与自上而下回采顺序相比，自下而上回采时主竖井底部的应力集中现象更为突出，且在开采过程中，主竖井井壁的应力增长速度相对较快，这表明自下而上回采对主竖井底部的稳定性影响较大。[此处插入自下而上回采顺序应力分布云图]间隔开采顺序下的应力分布云图（图[X3]）显示，主竖井周边的应力分布相对较为均匀。在开采过程中，由于间隔开采留下的矿体起到了一定的支撑作用，主竖井井壁的应力增长较为平缓。在整个开采过程中，主竖井井壁的最大主应力为[X'']MPa，明显低于自上而下和自下而上回采顺序下的最大主应力。这说明间隔开采能够有效地缓解应力集中现象，降低主竖井周边的应力水平，对主竖井的稳定性较为有利。[此处插入间隔开采顺序应力分布云图]5.1.2位移分布分析自上而下回采顺序的位移分布云图（图[X4]）表明，主竖井的位移主要表现为垂直向下的沉降和水平方向的位移。在开采初期，主竖井顶部的沉降量较小，约为[Y1]mm，但随着开采的进行，沉降量逐渐增大，在开采至矿体底部时，主竖井顶部的沉降量达到[Y1+ΔY1]mm。主竖井的水平位移也呈现出类似的变化趋势，在开采初期较小，随着开采的深入逐渐增大，在靠近矿体一侧的井壁，水平位移最大值达到[Z1]mm。这是由于矿体开采导致上覆岩层下沉和移动，从而带动主竖井产生相应的位移。[此处插入自上而下回采顺序位移分布云图]自下而上回采顺序的位移分布云图（图[X5]）显示，主竖井底部的位移在开采初期就较为明显，底部的沉降量达到[Y2]mm，随着开采向上推进，主竖井中部和顶部的位移也逐渐增大，在开采至矿体顶部时，主竖井顶部的沉降量达到[Y2+ΔY2]mm。与自上而下回采顺序相比，自下而上回采时主竖井底部的位移增长速度更快，这是因为下部采空区的形成对主竖井底部的影响更为直接，导致底部位移迅速增大。[此处插入自下而上回采顺序位移分布云图]间隔开采顺序下的位移分布云图（图[X6]）显示，主竖井的位移相对较小。在整个开采过程中，主竖井顶部的沉降量最大为[Y3]mm，水平位移最大值为[Z2]mm，均明显小于自上而下和自下而上回采顺序下的位移值。这进一步证明了间隔开采能够有效地控制主竖井的位移，提高主竖井的稳定性。[此处插入间隔开采顺序位移分布云图]5.1.3塑性区分布分析自上而下回采顺序的塑性区分布云图（图[X7]）显示，在开采初期，主竖井顶部附近出现了少量的塑性区，随着开采的进行，塑性区逐渐向下扩展。在开采至矿体中部时，塑性区范围明显增大，主竖井井壁周围的塑性区厚度达到[H1]m。这表明在自上而下回采过程中，主竖井顶部和中部的岩体逐渐进入塑性状态，岩体的强度和稳定性降低。[此处插入自上而下回采顺序塑性区分布云图]自下而上回采顺序的塑性区分布云图（图[X8]）显示，塑性区首先在主竖井底部出现，且在开采初期塑性区范围就较大，底部井壁周围的塑性区厚度达到[H2]m。随着开采向上推进，塑性区向上扩展，在开采至矿体顶部时，主竖井中部和顶部的塑性区范围也明显增大。与自上而下回采顺序相比，自下而上回采时主竖井底部的塑性区发展更为迅速，这说明自下而上回采对主竖井底部的稳定性破坏更为严重。[此处插入自下而上回采顺序塑性区分布云图]间隔开采顺序下的塑性区分布云图（图[X9]）显示，主竖井周围的塑性区范围较小，在整个开采过程中，主竖井井壁周围的塑性区厚度最大为[H3]m，明显小于自上而下和自下而上回采顺序下的塑性区厚度。这表明间隔开采能够有效地抑制塑性区的发展，减少主竖井周围岩体的破坏，从而提高主竖井的稳定性。[此处插入间隔开采顺序塑性区分布云图]综合以上不同回采顺序下主竖井周边应力、位移和塑性区分布云图的分析，可以得出结论：间隔开采顺序对主竖井的稳定性最为有利，能够有效降低主竖井周边的应力水平，控制位移和塑性区的发展；自上而下回采顺序次之；自下而上回采顺序对主竖井稳定性的影响最大，容易导致主竖井底部应力集中、位移增大和塑性区迅速扩展，增加主竖井失稳的风险。5.2不同开采深度下的模拟结果在不同开采深度的模拟方案中，通过FLAC3D软件对开采深度为[X]m、[X+100]m和[X+200]m三种情况进行了模拟，得到了主竖井在不同开采深度下的应力、位移和塑性区分布云图，以下将对这些云图进行详细分析。5.2.1应力分布分析当开采深度为[X]m时，从应力分布云图（图[X10]）可以看出，主竖井周边的应力分布相对较为均匀，最大主应力出现在主竖井靠近矿体一侧的井壁中部，大小为[X10]MPa。这是因为在该开采深度下，地应力相对较小，矿体开采对主竖井的影响范围和程度有限，应力集中现象不太明显。随着开采深度增加到[X+100]m，主竖井周边的应力明显增大（图[X11]），最大主应力达到[X11]MPa，且应力集中区域向主竖井底部扩展。这是由于开采深度的增加导致地应力增大，矿体开采后，岩体应力重新分布，使得主竖井周边的应力水平显著提高。[此处插入开采深度为[X]m时的应力分布云图][此处插入开采深度为[X+100]m时的应力分布云图]当开采深度达到[X+200]m时，主竖井周边的应力分布更加复杂（图[X12]），最大主应力进一步增大至[X12]MPa，在主竖井底部和靠近矿体一侧的井壁出现了明显的应力集中现象。此时，地应力的增大以及开采引起的应力扰动相互叠加，使得主竖井周边的岩体处于高应力状态，容易发生破坏。[此处插入开采深度为[X+200]m时的应力分布云图]5.2.2位移分布分析在开采深度为[X]m时，主竖井的位移较小（图[X13]），垂直位移最大值为[Y4]mm，水平位移最大值为[Z3]mm。主竖井的位移主要集中在井口附近，随着深度的增加，位移逐渐减小。这表明在该开采深度下，矿体开采对主竖井的位移影响较小，主竖井基本保持稳定。当开采深度增加到[X+100]m时，主竖井的位移明显增大（图[X14]），垂直位移最大值达到[Y5]mm，水平位移最大值达到[Z4]mm。主竖井的位移分布呈现出上部大、下部小的特点，井口处的位移增长最为显著。这是因为随着开采深度的增加，地应力增大，矿体开采对上覆岩层的扰动加剧，导致主竖井的位移增大。[此处插入开采深度为[X]m时的位移分布云图][此处插入开采深度为[X+100]m时的位移分布云图]当开采深度达到[X+200]m时，主竖井的位移进一步增大（图[X15]），垂直位移最大值为[Y6]mm，水平位移最大值为[Z5]mm。此时，主竖井的位移分布呈现出不均匀的状态，在井口和靠近矿体一侧的井壁位移较大，而在主竖井底部位移相对较小。这说明在深部开采条件下，主竖井的稳定性受到了较大的挑战，需要采取相应的措施来控制位移的发展。[此处插入开采深度为[X+200]m时的位移分布云图]5.2.3塑性区分布分析开采深度为[X]m时，主竖井周围的塑性区范围较小（图[X16]），主要集中在主竖井靠近矿体一侧的井壁顶部，塑性区厚度为[H4]m。这表明在该开采深度下，主竖井周围的岩体基本处于弹性状态，只有局部区域由于应力集中出现了少量的塑性变形。随着开采深度增加到[X+100]m，主竖井周围的塑性区范围明显扩大（图[X17]），塑性区向主竖井底部和远离矿体一侧扩展，塑性区厚度达到[H5]m。这是由于地应力增大和开采扰动的共同作用，使得主竖井周围更多的岩体进入塑性状态，岩体的强度和稳定性降低。[此处插入开采深度为[X]m时的塑性区分布云图][此处插入开采深度为[X+100]m时的塑性区分布云图]当开采深度达到[X+200]m时，主竖井周围的塑性区范围进一步扩大（图[X18]），塑性区几乎贯穿了主竖井的整个井壁，塑性区厚度达到[H6]m。此时，主竖井周围的岩体大部分已进入塑性状态，主竖井的稳定性受到了严重威胁。[此处插入开采深度为[X+200]m时的塑性区分布云图]综合以上不同开采深度下主竖井周边应力、位移和塑性区分布云图的分析，可以得出结论：随着开采深度的增加，主竖井周边的应力、位移和塑性区范围均呈现出逐渐增大的趋势，主竖井的稳定性逐渐降低。在深部开采过程中，需要密切关注主竖井的稳定性变化，采取有效的支护和加固措施，以确保主竖井的安全运行。5.3不同采场结构参数下的模拟结果在采场结构参数对主竖井稳定性影响的模拟研究中，重点分析了采场尺寸和形状变化对主竖井稳定性的影响。不同采场结构参数下主竖井周边应力、位移和塑性区分布云图清晰地展示了这些影响的具体表现。5.3.1采场尺寸影响分析在采场尺寸为长度[X]m、宽度[Y]m、高度[Z]m的情况下，从应力分布云图（图[X19]）可以看出，主竖井周边应力分布相对较为均匀，最大主应力出现在主竖井靠近矿体一侧的井壁，大小为[X19]MPa。这是因为此时采场尺寸相对较小，开采对主竖井周边岩体的扰动范围和程度有限，应力集中现象不太明显。主竖井的位移也相对较小（图[X20]），垂直位移最大值为[Y7]mm，水平位移最大值为[Z6]mm，位移主要集中在井口附近，随着深度的增加，位移逐渐减小。主竖井周围的塑性区范围较小（图[X21]），主要集中在主竖井靠近矿体一侧的井壁顶部，塑性区厚度为[H7]m，表明在该采场尺寸下，主竖井周围的岩体基本处于弹性状态，只有局部区域由于应力集中出现了少量的塑性变形。[此处插入采场尺寸为长度[X]m、宽度[Y]m、高度[Z]m时的应力分布云图][此处插入采场尺寸为长度[X]m、宽度[Y]m、高度[Z]m时的位移分布云图][此处插入采场尺寸为长度[X]m、宽度[Y]m、高度[Z]m时的塑性区分布云图]当采场尺寸增大为长度[X+20]m、宽度[Y+10]m、高度[Z+5]m时，主竖井周边的应力明显增大（图[X22]），最大主应力达到[X20]MPa，且应力集中区域向主竖井底部扩展。这是由于采场尺寸的增大导致开采形成的采空区范围扩大，岩体应力重新分布，使得主竖井周边的应力水平显著提高。主竖井的位移也明显增大（图[X23]），垂直位移最大值达到[Y8]mm，水平位移最大值达到[Z7]mm，井口处的位移增长最为显著。主竖井周围的塑性区范围明显扩大（图[X24]），塑性区向主竖井底部和远离矿体一侧扩展，塑性区厚度达到[H8]m，这表明在该采场尺寸下，由于开采扰动的加剧，主竖井周围更多的岩体进入塑性状态，岩体的强度和稳定性降低。[此处插入采场尺寸增大后的应力分布云图][此处插入采场尺寸增大后的位移分布云图][此处插入采场尺寸增大后的塑性区分布云图]5.3.2采场形状影响分析对于矩形采场，其应力分布云图（图[X25]）显示，主竖井周边的应力分布相对规则，最大主应力出现在主竖井靠近矿体一侧的井壁中部，大小为[X21]MPa。这是因为矩形采场的形状使得开采过程中岩体的应力分布相对均匀，对主竖井的影响相对稳定。主竖井的位移也相对较为均匀（图[X26]），垂直位移最大值为[Y9]mm，水平位移最大值为[Z8]mm，位移主要集中在井口和靠近矿体一侧的井壁。主竖井周围的塑性区范围较小（图[X27]），主要集中在主竖井靠近矿体一侧的井壁顶部和底部，塑性区厚度为[H9]m，表明矩形采场对主竖井稳定性的影响相对较小。[此处插入矩形采场的应力分布云图][此处插入矩形采场的位移分布云图][此处插入矩形采场的塑性区分布云图]而梯形采场的应力分布云图（图[X28]）显示，主竖井周边的应力分布不均匀，在梯形采场的斜边附近，主竖井井壁的应力明显增大，最大主应力达到[X22]MPa。这是因为梯形采场的形状导致开采过程中岩体的应力分布不均匀，在斜边附近容易产生应力集中现象，对主竖井的稳定性产生不利影响。主竖井的位移也呈现出不均匀的状态（图[X29]），在梯形采场斜边对应的主竖井井壁一侧，位移明显增大，垂直位移最大值为[Y10]mm，水平位移最大值为[Z9]mm。主竖井周围的塑性区范围较大（图[X30]），在主竖井靠近梯形采场斜边一侧的井壁，塑性区厚度达到[H10]m，表明梯形采场对主竖井稳定性的影响较大，容易导致主竖井周围岩体的破坏。[此处插入梯形采场的应力分布云图][此处插入梯形采场的位移分布云图][此处插入梯形采场的塑性区分布云图]综合以上不同采场结构参数下主竖井周边应力、位移和塑性区分布云图的分析，可以得出结论：较小尺寸的采场对主竖井的稳定性更为有利，能够有效降低主竖井周边的应力水平，控制位移和塑性区的发展；矩形采场相较于梯形采场，其形状对主竖井稳定性的影响较小，能够使主竖井周边的应力分布更加均匀，减少应力集中现象的发生。因此，在新城金矿Ⅺ#矿体的开采设计中，应优先选择较小尺寸的矩形采场结构参数，以保障主竖井在回采过程中的稳定性。5.4模拟结果综合分析通过对不同回采顺序、开采深度和采场结构参数下主竖井稳定性的数值模拟结果进行深入分析，总结出各因素对主竖井稳定性的影响规律，以及各因素之间的相互作用关系。回采顺序对主竖井稳定性的影响显著。自上而下回采时，随着开采向下推进，主竖井顶部和中部的应力逐渐增大，位移和塑性区也逐渐发展，这是由于上部采空区导致应力转移和岩体变形逐渐向下传递。自下而上回采则使得主竖井底部在开采初期就承受较大压力，应力集中、位移和塑性区发展在底部更为明显，且整个开采过程中主竖井井壁的应力增长速度较快，对底部稳定性破坏较大。间隔开采时，由于留下的矿体起到支撑作用，主竖井周边的应力分布相对均匀，应力水平、位移和塑性区发展都得到有效控制，对主竖井稳定性最为有利。开采深度是影响主竖井稳定性的关键因素之一。随着开采深度的增加，地应力显著增大，主竖井周边的应力、位移和塑性区范围均呈现逐渐增大的趋势。在浅部开采时，主竖井周边的应力、位移和塑性区相对较小，主竖井基本保持稳定。随着深度增加，地应力和开采扰动相互叠加，主竖井周边岩体处于高应力状态，塑性区扩展，岩体强度和稳定性降低，主竖井的稳定性受到严重威胁。采场结构参数对主竖井稳定性也有重要影响。在采场尺寸方面，较小尺寸的采场开采时，主竖井周边的应力、位移和塑性区都相对较小，因为采场尺寸小，开采形成的采空区范围小，对主竖井周边岩体的扰动有限。当采场尺寸增大时，采空区范围扩大，岩体应力重新分布，主竖井周边的应力、位移和塑性区明显增大，对主竖井稳定性不利。在采场形状方面，矩形采场由于其形状规则，开采过程中岩体应力分布相对均匀，对主竖井稳定性影响较小，主竖井周边的应力、位移和塑性区范围都相对较小。而梯形采场由于形状不规则，在斜边附近容易产生应力集中现象，导致主竖井周边的应力分布不均匀，位移和塑性区在斜边对应的井壁一侧明显增大，对主竖井稳定性影响较大。各因素之间存在复杂的相互作用关系。开采深度的增加会加大回采顺序和采场结构参数对主竖井稳定性的影响程度。在深部开采时，不同回采顺序下主竖井的应力、位移和塑性区的差异更为显著，采场结构参数的变化对主竖井稳定性的影响也更加明显。回采顺序和采场结构参数之间也相互影响，不合理的回采顺序可能会加剧采场结构参数对主竖井稳定性的不利影响，而合理的采场结构参数在一定程度上可以缓解回采顺序带来的负面影响。综上所述，在新城金矿Ⅺ#矿体的开采过程中，为了保障主竖井的稳定性，应优先选择间隔开采顺序，合理控制开采深度，避免过度开采导致地应力急剧增大。在采场结构参数设计方面，应采用较小尺寸的矩形采场，以减少对主竖井稳定性的影响。同时，要充分考虑各因素之间的相互作用关系，综合制定科学合理的开采方案和支护措施，确保主竖井在矿体回采过程中的安全稳定运行。六、主竖井稳定性评价与建议6.1主竖井稳定性评价指标在评估新城金矿主竖井稳定性时，位移、应力和塑性区范围是关键的评价指标，这些指标能直观反映主竖井在Ⅺ#矿体回采过程中的稳定状态。位移指标在稳定性评价中占据重要地位，它直接体现了主竖井在开采过程中的变形程度。竖向位移过大可能导致主竖井井筒出现沉降不均，影响提升设备的正常运行，甚至引发井筒与提升设备之间的碰撞等安全事故。水平位移过大则可能使主竖井井壁承受额外的剪切力，导致井壁出现裂缝、坍塌等问题。在实际工程中，通常根据主竖井的设计要求和相关规范，确定竖向位移和水平位移的允许阈值。一般来说，竖向位移的允许值可能控制在[X1]mm以内，水平位移的允许值控制在[X2]mm以内。一旦位移超过这些阈值，主竖井的稳定性就会受到严重威胁，需要及时采取相应的加固和修复措施。应力指标同样是评估主竖井稳定性的关键因素。井壁的最大主应力和剪应力反映了井壁所承受的力学荷载情况。当最大主应力超过井壁材料的抗压强度时，井壁会发生压缩破坏，表现为井壁混凝土破碎、剥落；剪应力超过井壁材料的抗剪强度时，井壁会出现剪切裂缝，进而导致井壁结构的失稳。对于新城金矿主竖井采用的钢筋混凝土支护结构，其混凝土的抗压强度设计值为[X3]MPa，抗剪强度设计值为[X4]MPa。在数值模拟和实际监测中，若发现井壁的最大主应力接近或超过[X3]MPa，剪应力接近或超过[X4]MPa，就需要警惕主竖井可能出现的失稳风险，及时调整开采方案或加强支护措施。塑性区范围是判断主竖井周围岩体稳定性的重要依据。塑性区的出现意味着岩体已经进入塑性变形阶段，其强度和承载能力下降。当塑性区范围过大，可能会连通形成破坏区域，导致主竖井周围岩体失去对井壁的有效支撑，从而引发主竖井的失稳。在新城金矿主竖井稳定性评价中，根据数值模拟结果和工程经验，若塑性区范围超过主竖井井壁周围一定距离，如[X5]m，就认为主竖井周围岩体的稳定性受到严重影响，需要采取岩体加固等措施来控制塑性区的进一步发展。基于上述分析，确定适用于新城金矿主竖井的稳定性评价标准如下：当主竖井的竖向位移小于[X1]mm，水平位移小于[X2]mm，井壁的最大主应力小于混凝土抗压强度设计值[X3]MPa，剪应力小于抗剪强度设计值[X4]MPa，塑性区范围在主竖井井壁周围[X5]m以内时，判定主竖井处于稳定状态；当其中任何一个指标超出上述范围时，判定主竖井处于不稳定状态，需要及时采取措施进行处理。6.2主竖井稳定性评价结果依据前文确定的位移、应力和塑性区范围等评价指标以及稳定性评价标准，对新城金矿主竖井在Ⅺ#矿体回采过程中的稳定性进行评价。在不同回采顺序下，间隔开采顺序下主竖井的各项指标表现最佳。主竖井的竖向位移最大值为[Y3]mm，小于允许阈值[X1]mm；水平位移最大值为[Z2]mm，小于允许阈值[X2]mm；井壁的最大主应力为[X'']MPa，小于混凝土抗压强度设计值[X3]MPa；剪应力小于抗剪强度设计值[X4]MPa；塑性区范围在主竖井井壁周围[H3]m以内，小于允许范围[X5]m。因此，在间隔开采顺序下，主竖井处于稳定状态。自上而下回采顺序时，主竖井的竖向位移最大值为[Y1+ΔY1]mm，接近允许阈值[X1]mm；水平位移最大值为[Z1]mm，也较为接近允许阈值[X2]mm；井壁的最大主应力为[X+ΔX1]MPa，虽然小于混凝土抗压强度设计值[X3]MPa，但已达到一定水平；塑性区范围在主竖井井壁周围[H1]m，也接近允许范围[X5]m。这种情况下，主竖井处于相对稳定状态，但需密切关注其变形和应力发展情况。自下而上回采顺序对主竖井稳定性影响较大，主竖井底部的竖向位移最大值达到[Y2+ΔY2]mm，超过允许阈值[X1]mm；水平位移最大值为[Z3]mm，也超过允许阈值[X2]mm；井壁底部的最大主应力为[X'+ΔX2]MPa，接近混凝土抗压强度设计值[X3]MPa，且剪应力较大；塑性区范围在主竖井底部井壁周围达到[H2]m，超过允许范围[X5]m。因此，在自下而上回采顺序下，主竖井处于不稳定状态，需要立即采取措施进行加固和调整开采方案。随着开采深度的增加，主竖井的稳定性逐渐降低。当开采深度为[X]m时，主竖井的竖向位移最大值为[Y4]mm，水平位移最大值为[Z3]mm，均小于允许阈值；井壁的最大主应力为[X10]MPa，小于混凝土抗压强度设计值[X3]MPa；塑性区范围在主竖井井壁周围[H4]m以内，小于允许范围[X5]m，主竖井处于稳定状态。当开采深度增加到[X+100]m时，主竖井的竖向位移最大值为[Y5]mm，接近允许阈值[X1]mm；水平