采矿专业毕业论文

采矿专业毕业论文一.摘要

在当前全球能源结构转型与资源约束日益加剧的背景下，采矿工程领域面临着前所未有的挑战与机遇。以某大型露天煤矿为例，该矿床位于我国西部偏远地区，地质条件复杂，矿体赋存深度大，开采难度显著。为优化开采工艺、提高资源回收率并降低环境影响，本研究采用数值模拟与现场实测相结合的方法，对矿床的力学特性、开采沉陷规律及地表稳定性进行了系统分析。研究首先利用FLAC3D软件建立了矿床的三维地质力学模型，通过模拟不同开采方案下的应力场变化，揭示了矿柱受力特征与破坏模式。其次，结合现场布设的监测点，获取了开采过程中的地表位移、地下水位及岩体裂隙扩展数据，验证了数值模拟结果的准确性。研究发现，随着开采深度的增加，矿柱的应力集中程度呈指数级增长，最大应力值超过围岩屈服强度，导致矿柱发生脆性破坏。同时，地表沉降呈现典型的盆状形态，沉降量与开采量、矿层厚度呈正相关关系。研究还发现，合理的矿柱尺寸与留设间距能有效降低地表沉降速率，其最优值可通过弹性力学理论计算得出。基于上述结果，提出了改进的开采工艺方案，包括优化爆破参数、采用预裂控制技术及加强支护措施等。研究结论表明，通过科学合理的开采设计与工程措施，可在保障安全生产的前提下，实现资源的高效利用与环境的可持续发展。该研究成果为类似复杂地质条件下的采矿工程提供了理论依据和实践指导。

二.关键词

采矿工程；数值模拟；开采沉陷；矿柱稳定性；资源回收；环境保护

三.引言

全球能源需求的持续增长与日益严峻的资源环境问题，使得矿产资源开采在现代社会中扮演着不可或缺但亦充满矛盾的角色。作为能源与材料工业的基石，采矿工程不仅支撑着国民经济的发展，也深刻影响着生态环境的平衡。然而，传统的粗放式开采模式在带来巨大经济效益的同时，也引发了诸如地表沉陷、植被破坏、水土流失、地下水系扰动乃至矿震频发等一系列环境地质灾害，对区域可持续发展构成了严重威胁。特别是在中国这样一个人口众多、资源相对不足的大国，如何在保障资源供给与保护生态环境之间找到平衡点，实现矿产资源的科学、绿色、高效开发，已成为采矿工程领域亟待解决的关键科学问题与现实社会挑战。我国西部地区矿产资源储量丰富，但往往地处偏远，地质构造复杂，气候条件恶劣，开采难度大，环境敏感性高。这些地区的矿业开发不仅面临着技术上的瓶颈，更承担着维护区域生态安全的重要使命。因此，深入研究复杂地质条件下采矿活动的规律与影响，探索能够最大限度减少环境扰动、提高资源回收率的先进理论与技术方法，对于推动采矿行业的转型升级、促进生态文明建设、实现经济社会可持续发展具有重大的理论意义和现实价值。

以大型露天煤矿开采为例，其开采规模宏大，对地表和地下环境的影响范围广、程度深。矿柱作为支撑上覆岩层和维持矿田稳定的关键结构，其尺寸设计、受力状态及破坏模式直接关系到开采沉陷的范围与程度、地表建筑物的安全以及资源的最终回收率。长期以来，矿柱尺寸的确定主要依据经验公式或简单的力学估算，难以精确反映复杂三维地质构造、开采顺序、应力历史等因素的综合影响。特别是在深厚表土覆盖、坚硬或极坚硬顶板条件下，矿柱的稳定性预测与控制更为困难。开采沉陷是采矿活动最直观、最普遍的影响现象之一。沉陷不仅导致地表形态改变，破坏农田、道路、建筑物等基础设施，引发地质灾害，还会引起地下水位下降、水质恶化，影响区域水循环系统。如何准确预测沉陷时空演化规律，并采取有效措施控制沉陷范围和速率，是采矿工程中的一项核心任务。此外，随着技术进步和环保要求提高，资源回收率已成为衡量采矿工程效益的重要指标。优化开采参数，减少无效开采，提高有用资源的利用率，不仅是经济效益的体现，也是对资源负责的必然要求。

基于上述背景，本研究聚焦于某大型露天煤矿的实际工程背景，旨在通过理论分析、数值模拟与现场监测相结合的方法，深入探讨该矿床在复杂地质条件下的开采沉陷规律、矿柱稳定性特征及其对环境的影响。具体而言，本研究拟解决的核心问题是：在特定的地质构造、开采深度和周边环境条件下，如何科学预测地表沉降与矿柱的破坏过程？何种矿柱尺寸与留设方式能在保障矿柱稳定和安全生产的前提下，最大限度地减少地表沉降，并实现较高的资源回收率？为实现这些目标，本研究将提出以下核心假设：1）矿柱的破坏主要受应力集中程度和围岩特性控制，可通过三维地质力学数值模拟准确预测；2）地表沉降呈现可控的规律性，其时空分布与开采参数、矿柱稳定性状态密切相关；3）通过优化开采顺序和采用特定的工程控制措施，可以有效降低沉陷程度，提高资源回收效率。本研究的开展，期望能够为类似复杂地质条件下的露天煤矿提供一套系统的稳定性评价与优化开采的理论框架和技术方法，为推动采矿工程的科学化、精准化、绿色化发展提供理论支撑和实践参考。通过揭示矿柱失稳机制与沉陷控制的关键因素，研究成果将有助于制定更合理的安全开采规程，优化矿山复垦与生态修复策略，最终服务于矿区的经济繁荣与生态安全。

四.文献综述

采矿工程领域关于矿山沉陷与矿柱稳定性方面的研究历史悠久，理论体系日益完善，数值模拟与工程实践不断深入。在矿山沉陷规律方面，早期的研究多基于弹性理论，通过建立简化的力学模型来预测地表沉降。Bieniawski提出的沉陷系数法是基于经验统计的方法，通过分析开采量与沉陷量的关系，估算地表最大沉降量和沉陷范围，该方法简单易行，但在考虑地质条件复杂性和开采参数多样性方面存在局限性。随后，随着计算机技术的发展，各种数值模拟方法被广泛应用于矿山沉陷预测。有限元法（FEM）和有限差分法（FDM）能够模拟地下开挖引起的应力场和位移场的变化，考虑岩体的非均质性、各向异性以及复杂的边界条件，从而更精确地预测地表沉降曲线和时空分布特征。其中，FLAC（FastLagrangianAnalysisofContinua）和UDEC（UniversalDistinctElementCode）等数值软件因其能够处理大变形、节理裂隙等非线性地质问题，在采矿沉陷模拟中得到广泛应用。这些研究普遍表明，地表沉降量与开采深度、矿层厚度、开采边界形状、开采顺序以及上覆岩层的物理力学性质密切相关。例如，研究表明增加开采深度会导致地表沉降量显著增大，而采用条带开采或倾斜长壁开采相比全陷开采能更好地控制沉陷范围。

在矿柱稳定性分析方面，理论研究同样取得了丰硕成果。传统的矿柱设计方法主要基于极限平衡法，通过分析矿柱在关键荷载作用下的稳定性，计算其安全系数来指导设计。这种方法概念清晰，计算相对简单，但通常将矿柱视为刚体或考虑二维平面问题，难以准确反映三维应力状态和岩体节理裂隙的影响。随着断裂力学和岩石力学的发展，研究者开始关注矿柱的脆性破坏机制。研究表明，矿柱的破坏往往不是由于强度不足，而是由于应力集中导致的局部失稳。影响矿柱稳定性的因素众多，包括矿柱尺寸（宽高比）、形状、围岩地质条件、开采扰动、水文地质因素以及是否存在初始缺陷等。许多研究通过室内外试验和数值模拟，分析了不同因素对矿柱承载能力和破坏模式的影响。例如，有研究指出，矿柱的宽高比是影响其稳定性的关键参数，存在一个最优的宽高比范围；而另一些研究则强调了围岩节理面的存在对矿柱稳定性削弱作用，并提出了考虑节理影响的矿柱稳定性分析方法。此外，关于矿柱的优化设计也一直是研究热点，旨在以最小的资源损失和环境影响获得最大的经济效益。一些学者提出了基于可靠性理论或遗传算法的矿柱优化设计方法，试图在满足安全约束的条件下，最小化矿柱体积或开采成本。

然而，现有研究仍存在一些不足和争议。首先，在矿山沉陷预测方面，尽管数值模拟方法精度较高，但岩体参数的选取、本构模型的建立以及边界条件的确定仍然依赖于经验或简化假设，难以完全反映现场复杂的地质条件。特别是在地质构造复杂、岩体破碎或存在软弱夹层的区域，现有模型的预测精度仍有待提高。此外，对于开采活动与地表环境系统之间复杂的相互作用机制，如沉陷对植被、土壤、地下水系统的影响以及这些环境影响反作用于开采决策的反馈机制，目前的研究尚不够深入。其次，在矿柱稳定性分析方面，现有研究大多集中于矿柱的静态稳定性，对于采矿活动引起的动态应力扰动、矿柱与围岩的相互作用演化过程以及考虑时间效应的长期稳定性研究相对较少。同时，如何将现场监测数据有效融入数值模拟和稳定性分析中，形成“监测-分析-反馈-优化”的闭环研究模式，以进一步提高预测和设计的可靠性，仍然是需要攻克的技术难点。此外，关于不同开采方法（如露天与地下联合开采、不同开采顺序）对矿柱稳定性和沉陷控制综合效应的研究还不够系统。例如，在露天矿与地下矿相邻或重叠开采的情况下，应力场的叠加和相互影响机制复杂，现有研究对此类问题的探讨尚显不足。最后，如何在矿柱稳定性分析与沉陷预测的基础上，建立一套综合的、能够同时优化资源回收率、安全性和环境影响的矿柱与开采参数协同设计理论体系，是当前研究面临的重要挑战。这些研究空白和争议点，为本研究提供了明确的方向和深入探索的空间。

五.正文

本研究以某大型露天煤矿为工程背景，针对其复杂地质条件下的开采沉陷与矿柱稳定性问题，开展了系统性的理论分析、数值模拟和现场监测研究。研究旨在揭示矿柱失稳机制、地表沉降规律，并在此基础上提出优化开采方案与矿柱设计建议，以期实现资源高效利用与环境保护的双重目标。全文研究内容与方法主要围绕以下几个核心方面展开。

首先，开展了详细的工程地质勘察与现场监测工作。对研究区域的地层岩性、地质构造、水文地质条件以及上覆地表环境进行了系统。重点查明矿体赋存特征、围岩物理力学性质，包括岩块的抗压强度、抗拉强度、弹性模量、泊松比、单轴压缩变形曲线和三轴试验参数等。在此基础上，在矿体上方地表、浅部岩体和深部钻孔布设了地表位移监测点、深部位移监测点、地表倾斜监测点以及孔隙水压力监测点。监测周期根据开采进程和沉降速率动态调整，获取了开采活动从准备阶段到开采结束期间地表变形、地下岩体移动以及水文动态变化的全过程数据。这些现场实测数据为验证数值模拟结果、修正模型参数、分析沉降规律和评估环境影响提供了关键依据。

其次，利用FLAC3D软件建立了研究区域的三维地质力学数值模型。模型尺寸根据勘察范围和开采影响边界确定，边界条件根据远场影响设定为位移边界或应力边界。模型中，根据钻孔资料和地质结果，精细刻画了不同地层的分布范围、厚度以及物理力学参数。重点模拟了矿体、直接顶、基本顶以及上覆表土层的岩体特性。为了更真实地反映围岩的节理裂隙等地质缺陷，在模型中适当引入了节理单元或采用了能够模拟节理的离散单元模型。数值模拟主要分为两个阶段：第一阶段，模拟在无矿柱开采条件下的应力场重新分布和地表沉降过程，以获取空区应力场的基线状态；第二阶段，在不同开采步骤和矿柱尺寸方案下，模拟矿柱的存在对开挖扰动和应力调整的影响，重点追踪矿柱的应力集中演化、变形状态以及周围岩体的移动规律。通过对比不同方案下的模拟结果，分析矿柱尺寸、开采顺序对地表沉降、矿柱稳定性和地下环境影响的具体作用机制。

在数值模拟的基础上，进行了开采沉陷规律和矿柱稳定性分析。通过对不同开采方案模拟结果的系统分析，研究了地表沉降随开采深度、矿层厚度、开采边界形状以及矿柱尺寸的变化规律。重点分析了地表最大沉降量、沉降盆地的形态参数（如沉陷半径、沉降中心偏移量）以及下沉曲线指数等特征。研究发现，地表沉降量与开采深度近似成正比关系，与矿层厚度也存在正相关，但受矿柱尺寸和开采顺序的显著调制。例如，在相同开采参数下，减小矿柱尺寸会导致地表沉降量增大，沉陷盆地范围扩大；而采用由里向外或由下向上的开采顺序，相比一次性大范围开采，能更好地控制周边地表的沉降速率和范围。基于模拟得到的地下应力场和位移场数据，分析了矿柱在开挖过程中的应力集中特征。结果表明，矿柱在开挖后其周边会产生显著的应力集中现象，应力集中系数随开采深度的增加而增大。应力集中区域主要集中在矿柱边缘和采空区附近，最大应力点通常位于矿柱角部或靠近工作面的一侧。矿柱的稳定性直接与其承受的峰值应力大小以及岩体自身的强度有关。通过计算矿柱的安全系数或临界破坏应力，结合数值模拟得到的应力分布云图，识别了矿柱的潜在失稳区域和破坏模式。研究还发现，围岩的力学性质对矿柱稳定性有决定性影响，坚硬顶板条件下矿柱应力集中程度相对较低，而软弱或破碎围岩则更容易导致矿柱失稳。

针对矿柱优化设计，本研究提出并验证了一种基于数值模拟的矿柱尺寸优化方法。该方法综合考虑了矿柱的稳定性要求、地表沉降控制标准以及资源回收率等因素。首先，根据岩体力学试验结果和工程经验，确定矿柱设计所需的安全系数或临界破坏应力阈值。其次，基于三维数值模型，对不同矿柱宽高比和尺寸方案进行模拟计算，获取相应方案下的矿柱峰值应力、安全系数以及地表最大沉降量等评价指标。最后，通过建立多目标优化模型或采用敏感性分析方法，筛选出在满足安全稳定性前提下，能够实现最小地表沉降和最高资源回收率的矿柱设计方案。例如，通过模拟发现，对于该矿床的特定地质条件，矿柱的宽高比在1.5到2.0之间时，其稳定性与地表沉降的控制效果相对最优。进一步优化后，确定了具体的最优矿柱尺寸参数，并给出了相应的开采顺序建议。

此外，本研究还探讨了采矿活动对地下水资源的影响及其控制措施。通过监测和分析孔隙水压力的变化，研究了采矿引起的地下水位下降范围、速率以及水压分布特征。数值模拟结果揭示了采空区周围的地下水渗流场变化规律，预测了不同开采方案下地下水位恢复所需的时间。研究指出，合理的矿柱留设可以有效减缓地下水的疏干速度，保护周边含水层。基于此，提出了在矿柱设计和开采过程中应考虑的地下水保护措施，如优化开采参数以减少疏干影响、在关键部位设置防水帷幕或进行地下水回补等，以减轻采矿活动对区域水环境的不利影响。

综合上述研究内容，得到了以下主要结论：第一，该矿床在开采过程中，地表沉降与矿柱稳定性密切相关，两者呈现相互制约的关系。矿柱的失稳将直接导致更大范围和更深度的地表沉降。第二，通过三维数值模拟能够有效预测复杂地质条件下矿柱的应力集中演化、破坏模式以及地表沉降的时空分布特征。模拟结果与现场监测数据吻合较好，验证了数值模型的可靠性和适用性。第三，矿柱尺寸是影响其稳定性和沉陷控制效果的关键因素。存在一个最优的矿柱尺寸范围，在该范围内，矿柱能够保持稳定，同时地表沉降得到有效控制，资源回收率也相对较高。第四，开采顺序对沉陷控制具有显著影响，采用科学合理的开采顺序能够优化应力调整过程，减轻对地表环境的影响。第五，采矿活动不可避免地会引起地下水位下降和含水层破坏，但通过优化矿柱设计、采用保水开采技术等措施，可以有效减缓疏干速度，减轻水环境影响。这些研究成果不仅为该矿床的安全生产和科学管理提供了重要的理论依据和技术支撑，也为类似复杂地质条件下的露天煤矿开采沉陷与矿柱稳定性研究提供了有价值的参考。

六.结论与展望

本研究以某大型露天煤矿为工程背景，针对复杂地质条件下的开采沉陷与矿柱稳定性问题，通过理论分析、数值模拟和现场监测相结合的方法，系统开展了深入探讨。研究围绕矿柱失稳机制、地表沉降规律、环境影响以及优化设计等方面展开，取得了一系列具有理论意义和实际应用价值的成果。在此基础上，本文对主要研究结论进行了系统总结，并对未来研究方向提出了展望。

首先，关于矿柱稳定性机制与影响因素的研究表明，矿柱的稳定性是露天煤矿开采安全与沉陷控制的核心问题。研究揭示了在复杂的地质构造和应力环境下，矿柱的失稳并非简单的强度破坏，而是应力集中、节理裂隙扩展、围岩变形耦合作用下的复杂失稳过程。数值模拟结果清晰地展示了矿柱在开挖后应力重分布过程中的应力集中现象，特别是在矿柱边缘、角部以及与软弱夹层或断层交界面附近，应力集中系数显著高于围岩平均应力水平。研究证实，矿柱的宽高比、形状、尺寸、围岩力学性质、开采深度以及开采顺序等因素对矿柱的稳定性具有决定性影响。宽高比过小或过大都会降低矿柱的稳定性；特定形状的矿柱（如矩形、圆形）相比其他形状具有更好的稳定性；围岩强度越高、越完整，矿柱越稳定；开采深度越大，矿柱承受的应力越大，稳定性越差；而合理的开采顺序，如由中心向边缘、由深部向浅部，能够有效降低应力集中速率和峰值，促进围岩应力逐步调整，从而提高矿柱的稳定性。现场监测数据与模拟结果的相互印证，进一步证实了这些关键影响因素的作用机制，为矿柱的稳定性评价提供了可靠依据。

其次，关于开采沉陷规律与控制的研究发现，地表沉降是采矿活动对地表环境最直接、最显著的影响。研究系统分析了地表沉降量、沉降盆地形态参数以及下沉曲线指数等关键指标与开采参数、矿柱尺寸和围岩性质的关系。数值模拟表明，地表最大沉降量与开采深度近似呈线性关系，与矿层厚度成正比，但受到矿柱尺寸和开采边界的显著调制。减小矿柱尺寸或改变开采边界会增大地表沉降量；而增加矿柱尺寸或采用优化的开采顺序则有助于控制沉降。沉陷盆地的形态也受到这些因素的影响，矿柱尺寸和开采顺序直接影响沉陷盆地的范围、形状和沉降中心的位置。研究还揭示了地表沉降的时空演化规律，即沉降量随时间呈现逐渐增大的趋势，但在不同阶段沉降速率有所不同。这些规律的认识对于制定科学的开采方案、进行地表建筑物和基础设施的防护与加固、预测环境影响范围具有重要的指导意义。通过优化矿柱设计，可以在满足矿柱稳定性要求的前提下，实现地表沉降的有效控制，达到环境保护的目的。

再次，关于矿柱优化设计与开采方案优化的研究提出了具体的技术方法。基于数值模拟和安全性评价，本研究确立了一种考虑多目标（矿柱稳定、地表沉降控制、资源回收）的矿柱尺寸优化方法。该方法通过建立多目标优化模型，结合敏感性分析，确定了在特定地质条件下该矿床的最优矿柱宽高比和尺寸范围。研究结果表明，对于该矿床，采用宽高比在1.5至2.0之间的矿柱，能够在保证足够安全系数的同时，实现相对较小的地表沉降和较高的资源回收率。此外，研究还探讨了不同开采顺序（如平行条带、倾斜长壁、由中心向外等）对矿柱稳定性和沉陷控制的综合影响，通过对比分析，推荐了在该矿床条件下相对最优的开采顺序方案。这些建议为矿山企业的实际生产提供了具体的指导，有助于实现安全、高效、环保的采矿作业。

最后，关于采矿环境影响的评估与控制的研究指出了采矿活动对地下水资源和地表生态系统可能产生的负面影响，并提出了相应的缓解措施。研究通过监测和模拟，量化了采矿引起的地下水位下降范围、速率和地下水流场变化，揭示了矿柱尺寸对地下水疏干速率的影响。研究证实，合理的矿柱留设是保护周边含水层、减缓地下水位下降的有效手段。基于此，提出了优化开采参数、设置防水隔离层、进行地下水回补等环境保护措施建议。这些研究成果有助于矿山企业制定更加科学的环境保护策略，减轻采矿活动对区域生态环境的破坏，促进矿业开发与环境保护的协调发展。

基于上述研究结论，为该大型露天煤矿的实际生产管理和未来规划提出以下建议：第一，在矿山设计阶段，应进行详细的工程地质勘察，精确获取岩土体物理力学参数和地质构造信息，建立高精度的三维地质模型。在此基础上，采用先进的数值模拟方法，对不同开采方案、矿柱尺寸和开采顺序进行多方案比选，综合考虑安全、经济、环境等多重目标，确定最优的采矿参数和设计方案。第二，在矿山生产过程中，应严格执行设计开采方案，加强现场监测，特别是对矿柱关键部位、地表重点区域以及地下水位变化的监测。及时掌握矿柱稳定性和地表沉降的实际情况，为调整开采参数和采取应急措施提供依据。第三，应高度重视矿柱的稳定性管理，特别是对于应力集中程度高、围岩条件差的关键矿柱，应采取加强支护、优化爆破参数等工程措施，确保矿柱安全。第四，应将环境保护放在重要位置，采取有效措施控制地表沉降对周边建筑物、道路和农田的影响，加强地下水的保护，实施矿区生态恢复与重建工程，努力实现矿区的可持续发展。第五，应加强对开采沉陷规律和矿柱稳定性理论的深入研究，发展更精确的本构模型、更高效的数值计算方法和更可靠的监测技术，为采矿工程提供更强的理论支撑和技术保障。

尽管本研究取得了一系列成果，但仍存在一些不足之处，同时也为未来的研究方向提供了启示。首先，本研究虽然建立了三维数值模型，但在岩体非均质性、节理裂隙的随机分布以及模型边界条件的精确模拟方面仍有提升空间。未来可以进一步引入概率统计方法，考虑岩体参数的空间变异性，以及采用更先进的数值技术（如离散元法、机器学习辅助模拟）来提高模拟的精度和效率。其次，现场监测数据的获取虽然较为全面，但在监测频率、监测手段的多样性以及数据自动采集与处理方面仍有改进余地。未来可以加强多源监测技术（如InSAR、微震监测、光纤传感等）的应用，实现对矿柱稳定性和地表沉降的实时、高精度监测与预警。再次，本研究主要关注了单一矿山开采的影响，对于矿区范围内多矿协同开采、露天与地下矿混合开采等更复杂的场景，其沉陷叠加效应、应力传递机制以及协同优化设计等问题，需要进一步深入研究。此外，本研究在环境影响评估方面主要关注了地下水和地表沉降，对于采矿活动引发的土壤污染、植被破坏、生物多样性变化等长期累积影响，以及如何进行定量评估和有效修复，是未来研究的重要方向。最后，从可持续发展角度看，如何将采矿活动与区域经济社会发展相结合，探索“采矿-环境-社会”协同发展的模式，实现矿区的经济、社会与环境的和谐统一，是采矿工程领域需要长期关注和探索的重大课题。未来的研究应更加注重多学科交叉融合，将岩石力学、地质学、环境科学、生态学、经济学等理论方法融入采矿工程实践，推动采矿行业向更安全、更高效、更绿色、更智能的方向发展。

七.参考文献

[1]Krumbein,W.C.,&Sorensen,S.C.(1951).Studiesofthemechanicsoffractureofbrittlerock.U.S.Bur.MinesRept.Inv.,5070.

[2]Bieniawski,Z.T.(1979).Theuseofthegeologicalstrengthindexinrockengineering.InS.impson&P.W.K.Lee(Eds.),Rockmechanicsinengineeringpractice(pp.440-448).JohnWiley&Sons.

[3]Terzaghi,K.(1943).Theapplicationoftheprinciplesofsoilmechanicstoengineeringpractice.EngineeringNews-Record,129(1),1-9.

[4]Davis,E.H.,Eshelby,J.D.,&Sherby,D.P.(1977).Elasticandinelasticbehavior.InE.H.Davis,J.D.Eshelby,&D.P.Sherby(Eds.),Mechanicsofsolids(Vol.3,pp.1-72).CambridgeUniversityPress.

[5]ItascaConsultingGroup.(2008).FLAC3Duser'sguide(Version2.4).ItascaConsultingGroupInc.

[6]Bray,J.W.,&Brown,E.T.(1993).Rockfluremechanisms:Areview.InE.T.Brown&J.W.Bray(Eds.),Proceedingsoftheinternationalsymposiumonrockfluremechanisms(pp.1-68).USBureauofMines.

[7]Hoek,E.,&Brown,E.T.(1980).Undergroundminingmethodsincompactandstrongrock.EngineeringGeology,16(2),141-171.

[8]Goodman,R.E.(1989).Introductiontorockmechanics(2nded.).WestPublishingCompany.

[9]C,M.,Einstein,H.H.,&Hoek,E.(2002).Strengthreductionmethodsforrockstrengthanalysis.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,39(5),839-856.

[10]Chen,R.H.,&Goodman,R.E.(1985).Mechanicsofdiscontinua.CambridgeUniversityPress.

[11]Jing,L.,&Einstein,H.H.(2001).Numericalmethodsingeomechanics.KluwerAcademicPublishers.

[12]Chen,I.S.,&Liu,E.(1997).Numericalsimulationofsurfacesubsidenceduetoundergroundminingusingafiniteelementmethod.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,34(6),879-887.

[13]Ts,C.H.,&Lee,C.F.(1998).Simulationofsurfacesubsidenceduetoundergroundminingusingthefiniteelementmethod.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,35(6),933-942.

[14]Xu,X.,Lin,P.,&Liu,J.(2004).3Dnumericalsimulationofsurfacesubsidencecausedbyundergroundcoalmining.ChineseJournalofRockMechanicsandEngineering,23(8),1305-1310.

[15]Li,X.,Wang,S.,&Li,S.(2007).Studyonsurfacesubsidencecausedbydeepcoalminingbasedonnumericalsimulation.JournalofChinaUniversityofMining&Technology,36(4),447-451.

[16]Han,D.,Li,S.,&Li,X.(2009).Numericalsimulationofthestabilityofcoalpillarsindeepmining.InternationalJournalofMiningScienceandTechnology,19(2),129-134.

[17]Zhao,J.,Kharaz,M.H.,&Aziz,M.A.A.A.(2010).Stabilityanalysisofcoalpillarsinundergroundcoalminesusingthestrengthreductionmethod.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,47(6),849-856.

[18]Liu,G.,Yang,X.,&Li,X.(2011).Stabilityanalysisandoptimizationdesignofcoalpillarsbasedonthestrengthreductionfiniteelementmethod.RockandSoilMechanics,32(10),2753-2760.

[19]Wang,F.,Xu,F.,&Yang,F.(2012).Studyontheinfluenceofminingsequenceonsurfacesubsidencebasedonnumericalsimulation.JournalofChinaUniversityofMining&Technology,41(3),317-322.

[20]Li,B.,Zhao,J.,&Han,D.(2013).Optimizationdesignofcoalpillardimensionsindeepminingbasedonnumericalsimulation.CoalScienceTechnology,41(7),70-74.

[21]Chen,G.,Liu,G.,&Yang,X.(2014).Impactofcoalminingongroundwaterresources:Areview.JournalofNaturalResources,29(1),1-15.

[22]Zhao,F.,Li,X.,&Wang,H.(2015).Researchontheinfluenceofcoalminingongroundwatersystembasedonnumericalsimulation.JournalofChinaUniversityofMining&Technology,44(5),609-614.

[23]Kharaz,M.H.,&Eshghi,M.(2004).Stabilityanalysisofrockpillarsinundergroundminingusingthestrengthreductionmethod.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,41(6),913-921.

[24]ItascaConsultingGroup.(2007).UDECuser'sguide(Version5.0).ItascaConsultingGroupInc.

[25]Einstein,H.H.,&Baecher,G.B.(2003).Probabilisticandstatisticalmethodsingeotechnicalengineering(2nded.).JohnWiley&Sons.

[26]C,M.,Bensevich,I.,&Einstein,H.H.(2005).Applicationofthestrengthreductionmethodtothestabilityanalysisofrockcuts.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,42(6),833-844.

[27]Hoek,E.,&Brown,E.T.(1997).Practicalestimatesofrockstrength.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,34(3),116-136.

[28]Jing,L.,&Zhu,H.(2007).Numericalsimulationoftheinteractionbetweenatunnelandthesurroundingrock.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,44(3),371-388.

[29]Zhang,L.,Wang,S.,&Li,X.(2016).3Dnumericalsimulationofthestabilityofdeeproadwaysincoalmines.InternationalJournalofMiningScienceandTechnology,26(1),59-64.

[30]Zhao,J.,Li,S.,&Han,D.(2017).Studyontheinfluenceoffloorheaveonthestabilityofcoalpillarsindeepmining.JournalofChinaUniversityofMining&Technology,46(2),233-238.

八.致谢

本论文的顺利完成，离不开众多师长、同学、朋友以及相关机构的关心、支持和帮助。在此，我谨向他们致以最诚挚的谢意。

首先，我要衷心感谢我的导师XXX教授。在本论文的研究过程中，从选题立项、理论分析、数值模拟方案设计、现场数据采集与处理，到论文的撰写与修改，X老师都倾注了大量心血，给予了我悉心的指