泊江海子矿围岩蠕变特性解析与本构关系构建

泊江海子矿围岩蠕变特性解析与本构关系构建一、绪论1.1研究背景与目的在煤炭资源开发中，煤矿巷道的稳定性至关重要，而围岩的力学特性是影响巷道稳定性的关键因素之一。泊江海子矿作为重要的煤炭生产基地，其开采活动面临着复杂的地质条件。该矿井田位于鄂尔多斯高原之西北部，区域性地表分水岭“东胜梁”的北侧，属高原半沙漠地貌。井田构造形态总体为一向北西倾斜的单斜构造，倾向N30°-40°W，地层倾角1-3°，区内虽无岩浆岩侵入，也未发现大的断裂和褶皱构造，但围岩主要由砂岩、砂质泥岩、泥岩和煤层组成。其中，砂质泥岩、泥岩质地较软，强度低，遇水软化、强度急剧下降，甚至砂化、泥化，成为松散体；砾岩及各粒级的砂岩均为孔隙式砂泥质、泥质胶结，胶结较疏松，失水后易干裂破碎。随着开采深度和强度的不断增加，泊江海子矿围岩所承受的地应力不断增大，围岩的变形和破坏问题日益突出。在煤矿开采过程中，围岩会受到开挖扰动、地应力、地下水等多种因素的作用，其力学行为表现出明显的时间相关性，即流变特性。蠕变作为岩石流变的一种重要表现形式，是指岩石在恒定应力作用下，随时间发生的不可逆变形。蠕变过程通常包括初始蠕变、稳态蠕变和加速蠕变三个阶段。在初始蠕变阶段，岩石的应变随时间快速增加；进入稳态蠕变阶段，应变增加速度逐渐稳定；而在加速蠕变阶段，应变速度急剧增大，最终导致岩石的破坏。蠕变现象的存在对岩石工程的长期稳定性具有重要影响，可能导致巷道变形、支护结构失效，甚至引发安全事故。岩石流变力学的基本原理为岩石工程的长期稳定性分析提供了理论基础。在实际工程中，综合考虑岩石的蠕变、松弛和长期强度特性，是确保工程安全性和稳定性的关键。近年来，随着科学技术的不断进步，岩石流变力学在工程中的应用日益广泛，涵盖了从地下资源开发到大型基础设施建设的多个领域。在矿山工程中，通过对岩石流变特性的研究，可以预测矿山的变形和破坏过程，为矿山的开采提供安全保障。因此，开展泊江海子矿围岩蠕变试验及本构关系研究具有重要的现实意义。本研究旨在通过对泊江海子矿围岩进行蠕变试验，深入了解其蠕变特性和变形规律，建立准确的本构关系模型，为该矿的巷道支护设计、开采方案优化以及安全生产提供科学依据。具体来说，通过试验获取围岩在不同应力条件下的蠕变数据，分析其蠕变过程中的应变随时间的变化规律、不同阶段的特征以及破坏形式；基于试验结果，建立能够准确描述围岩蠕变行为的本构关系模型，并对模型参数进行辨识和优化，使其能够更好地反映围岩的实际力学行为，从而为工程实践提供可靠的理论支持。1.2国内外研究现状1.2.1岩石流变理论研究岩石流变理论的研究历史悠久，早在19世纪，国外学者就开始关注岩石的流变现象。随着时间的推移，岩石流变理论不断发展，先后出现了弹性理论、塑性理论和粘弹性理论等。这些理论从不同角度对岩石的力学行为进行了描述，为岩石流变学的发展奠定了基础。在弹性理论中，胡克定律的提出为描述岩石的弹性行为提供了基本依据；塑性理论则通过引入屈服准则，对岩石的塑性变形进行了研究；粘弹性理论将岩石视为弹性元件和粘性元件的组合，考虑了时间因素对岩石变形的影响，更能反映岩石在实际工程中的力学行为。近年来，随着分数阶微积分、非局部弹性理论等新的数学模型和计算方法的引入，岩石流变理论得到了进一步的发展。分数阶微积分能够更准确地描述岩石流变过程中的记忆特性和遗传效应，使模型能够更全面地反映岩石在长时间尺度下的变形规律。非局部弹性理论考虑了岩石内部微观结构的相互作用，引入了空间相关性，从而能够更准确地描述岩石在复杂应力状态下的力学行为。这些新的理论和方法的应用，使得岩石流变理论能够更精确地描述岩石的流变行为，为岩石工程的设计和分析提供了更有力的理论支持。在国内，岩石流变理论的研究也取得了显著进展。众多学者结合我国丰富的工程实践，对岩石流变特性进行了深入研究，提出了许多具有创新性的理论和方法。例如，一些学者通过对岩石微观结构的研究，揭示了岩石流变的微观机制，为建立更准确的岩石流变本构模型提供了理论基础。同时，国内学者还在岩石流变理论的工程应用方面进行了大量研究，将岩石流变理论与实际工程相结合，解决了许多工程中的关键问题，推动了岩石流变理论在工程中的广泛应用。1.2.2岩石流变数值研究随着计算机技术的飞速发展，数值模拟方法在岩石流变研究中得到了广泛应用。有限元法、有限差分法、离散元法等数值方法成为研究岩石流变行为的重要工具。有限元法通过将连续的岩石介质离散为有限个单元，将复杂的岩石力学问题转化为线性代数方程组进行求解，能够准确地模拟岩石在复杂应力条件下的变形和应力分布。有限差分法是一种基于差分原理的数值方法，它将岩石的连续介质模型离散为网格，通过对网格节点上的物理量进行差分计算，求解岩石的力学响应，在处理一些具有明显边界条件的岩石流变问题时具有独特的优势。离散元法则适用于模拟岩石的非连续变形和破坏过程，它将岩石视为由离散的颗粒组成，通过考虑颗粒之间的相互作用来模拟岩石的力学行为，能够很好地反映岩石在破裂和失稳过程中的力学特性。国外在岩石流变数值研究方面起步较早，取得了一系列重要成果。许多学者利用数值模拟方法对各种岩石工程问题进行了研究，如隧道开挖、边坡稳定性分析等，为工程设计和施工提供了重要的参考依据。例如，通过数值模拟研究隧道开挖过程中围岩的流变特性，分析围岩的变形规律和应力分布，预测隧道的长期稳定性，为隧道支护设计提供科学依据。在边坡稳定性分析中，利用数值模拟方法考虑岩石的流变特性，评估边坡在长期荷载作用下的稳定性，为边坡的加固和治理提供决策支持。国内在岩石流变数值研究方面也取得了长足的进步。众多科研团队和学者针对我国复杂的地质条件和工程需求，开展了深入的研究工作。通过自主研发和改进数值模拟软件，结合实际工程案例，对岩石流变行为进行了精细化模拟，为我国的岩石工程建设提供了有力的技术支持。例如，在一些大型水利工程、矿山开采工程中，利用数值模拟方法考虑岩石的流变特性，优化工程设计方案，确保工程的长期安全稳定运行。同时，国内学者还在数值模拟方法的创新和发展方面做出了积极贡献，提出了一些新的数值算法和模型，提高了数值模拟的精度和效率。1.2.3岩石流变试验研究岩石流变试验是研究岩石流变特性的重要手段，通过试验可以直接获取岩石在不同应力条件下的应变随时间的变化规律，为岩石流变理论和数值模拟研究提供基础数据。岩石流变试验包括单轴蠕变试验、三轴蠕变试验、松弛试验等。单轴蠕变试验是在恒定的轴向应力作用下，测量岩石试样的轴向应变随时间的变化，能够直观地反映岩石在单轴应力状态下的蠕变特性。三轴蠕变试验则是在恒定的围压和轴向应力作用下，研究岩石在三维应力状态下的流变行为，更符合岩石在实际工程中的受力状态。松弛试验是在恒定的应变条件下，测量岩石试样的应力随时间的变化，用于研究岩石的松弛特性。在国外，岩石流变试验研究已经开展了多年，积累了丰富的经验和大量的试验数据。许多先进的试验设备和技术不断涌现，如电液伺服岩石力学试验系统、高精度应变测量装置等，为岩石流变试验提供了可靠的技术支持。同时，国外学者还在试验方法和数据处理方面进行了深入研究，提出了一些新的试验方案和数据处理方法，提高了试验结果的准确性和可靠性。国内在岩石流变试验研究方面也取得了显著成就。近年来，随着我国对岩石工程建设的重视和投入不断增加，岩石流变试验研究得到了快速发展。许多高校和科研机构建立了先进的岩石力学试验平台，具备了开展各种岩石流变试验的能力。通过大量的试验研究，国内学者对我国不同地区、不同类型岩石的流变特性有了更深入的了解，为我国岩石工程的设计和施工提供了重要的试验依据。同时，国内学者还在试验设备的研发和改进方面做出了努力，提高了试验设备的性能和自动化程度，推动了我国岩石流变试验研究的发展。1.3研究内容与方法1.3.1研究内容泊江海子矿围岩基本力学特性研究：对取自泊江海子矿的围岩岩样进行常规单轴压缩试验和常规三轴压缩试验。在单轴压缩试验中，依据相关标准，利用高精度的岩石力学试验设备，严格控制加载速率，记录岩样在加载过程中的应力-应变数据，从而获取泊江海子矿围岩的单轴抗压强度、弹性模量、泊松比等基本力学参数。通过对试验数据的分析，深入研究围岩在单轴应力状态下的变形特征和破坏模式，为后续的研究提供基础数据和理论依据。在三轴压缩试验中，模拟围岩在实际工程中的三维受力状态，施加不同的围压和轴向压力，全面分析围压对围岩力学特性的影响规律，包括强度变化、变形特征以及破坏机制等方面的变化，为深入理解围岩在复杂应力条件下的力学行为提供支持。泊江海子矿围岩蠕变特性试验研究：开展岩石单轴蠕变试验和三轴蠕变试验。在单轴蠕变试验中，选用合适的蠕变试验仪器，精心制备岩样，严格按照试验标准和步骤，对岩样施加不同等级的恒定轴向应力，持续监测并记录岩样的轴向应变随时间的变化情况。通过对试验数据的深入分析，全面研究单轴应力状态下围岩的蠕变特性，包括蠕变曲线的特征、蠕变各阶段的变形规律、蠕变速度随时间的变化以及不同应力水平对蠕变特性的影响等方面。在三轴蠕变试验中，同样模拟实际工程中的三维应力状态，对岩样施加不同的围压和轴向应力组合，详细分析围压和轴向应力对围岩蠕变特性的综合影响，揭示围岩在复杂应力条件下的蠕变变形机制和规律。同时，对蠕变试验后的岩样进行仔细观察和分析，研究其破坏形式和破坏机理，为进一步理解围岩的蠕变破坏过程提供依据。泊江海子矿围岩蠕变本构模型的建立及参数辨识：基于岩石流变理论，结合泊江海子矿围岩的蠕变试验结果，建立能够准确描述其蠕变行为的本构模型。考虑到岩石的弹性、塑性和粘性等特性，以及蠕变过程中的损伤演化，选择合适的流变元件进行组合，构建合理的本构模型结构。通过对试验数据的非线性拟合和参数优化，辨识出本构模型中的各个参数，使模型能够精确地反映围岩在不同应力条件下的蠕变特性。同时，对建立的本构模型进行严格的验证和对比分析，与其他常用的本构模型进行比较，评估模型的准确性和适用性，确保模型能够为工程实际提供可靠的理论支持。基于蠕变本构模型的数值模拟研究：将建立的围岩蠕变本构模型引入到数值模拟软件中，对泊江海子矿巷道开挖过程进行数值模拟分析。在模拟过程中，充分考虑围岩的力学特性、蠕变特性以及支护结构的作用，真实地模拟巷道开挖后的应力分布、变形演化以及支护结构的受力情况。通过对不同支护方案和开采工艺的模拟分析，研究其对巷道稳定性的影响规律，为巷道支护设计和开采方案的优化提供科学依据和技术支持。同时，将数值模拟结果与现场监测数据进行对比验证，进一步评估本构模型和数值模拟方法的可靠性和准确性。1.3.2研究方法室内试验方法：在实验室环境中，利用先进的岩石力学试验设备，如电液伺服岩石力学试验系统、高精度应变测量装置等，对泊江海子矿围岩岩样进行基本力学特性试验和蠕变特性试验。严格按照相关试验标准和规范，精心制备岩样，精确控制试验条件，确保试验数据的准确性和可靠性。通过室内试验，直接获取围岩在不同应力状态下的力学参数和蠕变数据，为后续的理论分析和数值模拟提供坚实的数据基础。理论分析方法：深入研究岩石流变力学的基本理论，结合室内试验结果，对泊江海子矿围岩的力学行为和蠕变特性进行全面的理论分析。运用弹性力学、塑性力学、粘弹性力学等相关理论，建立合理的力学模型，深入分析围岩在不同应力条件下的变形机制和破坏机理。同时，对建立的蠕变本构模型进行严格的理论推导和分析，明确模型的物理意义和适用范围，为模型的参数辨识和工程应用提供理论依据。数值模拟方法：借助专业的数值模拟软件，如ANSYS、FLAC3D等，将室内试验得到的力学参数和建立的蠕变本构模型应用于数值模拟中。通过数值模拟，对泊江海子矿巷道开挖过程进行全面的模拟分析，深入研究围岩的应力分布、变形演化以及支护结构的受力情况。利用数值模拟方法，可以快速、准确地预测不同条件下巷道的稳定性，为巷道支护设计和开采方案的优化提供科学依据和技术支持。同时，通过与现场监测数据的对比验证，不断优化数值模拟模型和参数，提高模拟结果的准确性和可靠性。1.4技术路线本研究采用理论分析、室内试验与数值模拟相结合的技术路线，深入开展泊江海子矿围岩蠕变试验及本构关系研究，具体如下：资料收集与理论研究：全面收集泊江海子矿的地质资料，包括地层分布、岩石类型、地质构造等，以及矿井开采现状，如开采深度、开采范围、巷道布置等信息。同时，深入研究岩石流变力学的基本理论，系统分析国内外相关研究现状，充分了解岩石流变理论、数值研究和试验研究的最新进展，为后续研究奠定坚实的理论基础。室内试验：在实验室环境下，对泊江海子矿围岩岩样开展基本力学特性试验和蠕变特性试验。利用先进的电液伺服岩石力学试验系统进行常规单轴压缩试验和常规三轴压缩试验，严格按照相关标准和规范，精确控制试验条件，获取围岩的单轴抗压强度、弹性模量、泊松比等基本力学参数，以及围压对围岩力学特性的影响规律。采用高精度应变测量装置进行单轴蠕变试验和三轴蠕变试验，对岩样施加不同等级的恒定应力，持续监测并记录岩样的应变随时间的变化情况，深入研究围岩的蠕变特性，包括蠕变曲线特征、各阶段变形规律以及破坏形式和机理。本构模型建立与参数辨识：基于岩石流变理论和室内试验结果，综合考虑岩石的弹性、塑性和粘性等特性，以及蠕变过程中的损伤演化，选用合适的流变元件进行组合，构建能够准确描述泊江海子矿围岩蠕变行为的本构模型。运用非线性拟合和参数优化方法，对试验数据进行分析处理，精确辨识出本构模型中的各个参数，使模型能够精准地反映围岩在不同应力条件下的蠕变特性。通过与其他常用本构模型进行对比分析，严格验证所建模型的准确性和适用性。数值模拟：将建立的围岩蠕变本构模型引入专业的数值模拟软件，如ANSYS、FLAC3D等，对泊江海子矿巷道开挖过程进行全面的数值模拟分析。在模拟过程中，充分考虑围岩的力学特性、蠕变特性以及支护结构的作用，真实地模拟巷道开挖后的应力分布、变形演化以及支护结构的受力情况。通过对不同支护方案和开采工艺的模拟分析，深入研究其对巷道稳定性的影响规律，为巷道支护设计和开采方案的优化提供科学依据和技术支持。将数值模拟结果与现场监测数据进行对比验证，进一步评估本构模型和数值模拟方法的可靠性和准确性，不断优化模型和参数，提高模拟结果的精度。结果分析与应用：对室内试验结果、本构模型计算结果和数值模拟结果进行综合分析，深入探讨泊江海子矿围岩的力学行为和蠕变特性，总结其变形规律和破坏机制。根据分析结果，为泊江海子矿的巷道支护设计提供合理的参数建议，如支护形式、支护强度、支护时间等；为开采方案的优化提供科学依据，如开采顺序、开采速度、开采方法等，以确保矿井的安全生产和长期稳定运行。同时，对研究成果进行总结和归纳，为类似矿山工程的岩石流变研究提供参考和借鉴。二、泊江海子矿工程地质条件分析2.1矿区地质概况泊江海子矿位于鄂尔多斯高原之西北部，区域性地表分水岭“东胜梁”的北侧，行政区划属鄂尔多斯市东胜区管辖。其地理位置优越，交通较为便利，矿井工业广场布置在泊江海子镇北部，距109国道2.5km处，距乌蒙公路30m。该区域属高原半沙漠地貌，地形较为平坦，有利于矿山的建设和开采活动。从地层分布来看，泊江海子矿地层属于华北地层大区晋冀鲁豫地层区鄂尔多斯地层分区的中生代大型内陆坳陷盆地沉积，出露地层主要有中生界三叠系上统、侏罗系下统和中统、白垩系下统，以及新生界新近系上新统和第四系。地层产状平缓，倾向西偏南，倾角在0-3°之间。在这些地层中，侏罗系下统延安组岩性为蓝灰、浅灰色粉砂岩、粉砂质泥岩夹淡黄色块状中粒砂岩、细砂岩，中、下部夹有煤线和多层可采煤层，底部有一层2-15m的灰白色含砾中粒石英砂岩，是与下伏地层的分层标志，区域出露厚度168m，与下伏三叠系延长组呈假整合或微角度不整合接触。侏罗系中统直罗组岩性为一套灰白、灰黄、灰绿、紫红色中细粒砂岩、粉砂岩夹蓝灰色泥岩和煤线，其底部常含砾石、铁质结核及龟背石，龟背石直径达0.5-0.8m，出露厚度96.31m，与下伏延安组为整合接触。白垩系下统是工作区分布最广的地层，由红色大型斜层理发育的砂岩与杂色中小型斜层理的砂岩、细砂岩及泥岩组成。矿区构造形态总体为一向北西倾斜的单斜构造，倾向N30-40°W，地层倾角1-3°，地层产状沿走向及倾向均有一定变化，但变化不大。虽无岩浆岩侵入，也未发现大的断裂和褶皱构造，但小的构造变动仍对岩体的完整性和力学性质产生一定影响。在局部区域，由于小型断层或节理的存在，岩体的连续性被破坏，导致其强度降低，在矿山开采过程中容易引发围岩的变形和破坏。这些小构造的存在使得地应力分布更加复杂，对巷道的稳定性产生不利影响。泊江海子矿可采煤层有3-1、3-1下、3-2、4-1、4-2、5-1计6层。其中，3-1煤层厚度0.21-7.26m，平均3.87m，为薄-厚煤层，可采性指数Km为0.91，煤厚变异系数γ为50%，面积可采率为95%，结构简单，是主采煤层，连续性好，评为稳定煤层；4-2煤层厚度0-3.06m，平均1.27m，为薄-中厚煤层，Km为0.55，γ为70%，面积可采率为61%，结构简单，连续性较好，评为较稳定煤层；5-1煤层厚度0-4.25m，平均2.80m，为薄-厚煤层，Km为0.95，γ为33%，面积可采率为96%，结构简单，连续性好，评为稳定煤层；3-1下煤层厚度0-1.30m，平均0.75m，为薄煤层，Km为0.56，γ为37%，面积可采率为26%，结构简单，可采区连续性较差，评为不稳定煤层；3-2煤层厚度0-2.07m，平均0.59m，以薄煤层为主，Km为0.38，γ为74%，面积可采率为36%，结构简单，可采区的连续性较差，评为不稳定煤层；4-1煤层厚度0-2.20m，平均0.60m，以薄煤层为主，Km为0.20，γ为60%，面积可采率为22%，结构简单，可采区连续性较差，评为不稳定煤层。稳定和较稳定煤层的资源储量占总资源储量的89%，稳定煤层的资源储量占总资源储量的77%，综合评定本矿煤层稳定程度为较稳定型。目前，泊江海子矿的开采深度和强度不断增加。随着开采深度的增加，地应力增大，围岩所承受的压力也相应增大，这使得围岩的变形和破坏问题更加突出。同时，开采强度的增加导致开采速度加快，对围岩的扰动更加频繁，进一步加剧了围岩的不稳定。在一些深部开采区域，已经出现了巷道变形、底鼓等问题，严重影响了矿山的安全生产和正常运营。因此，深入研究泊江海子矿的工程地质条件，对于解决矿山开采过程中的围岩稳定性问题具有重要意义。2.2围岩岩性特征泊江海子矿围岩主要由白垩系中砂岩、侏罗系砂质泥岩等构成，不同岩石的矿物成分、结构构造及物理力学性质各异，这些特性对围岩的稳定性及力学行为有着关键影响。白垩系中砂岩主要由石英、长石等矿物组成，石英含量通常在50%-70%之间，长石含量约为20%-30%，此外还含有少量的云母、黏土矿物及其他杂质。其结构多为中粒砂状结构，颗粒大小相对均匀，粒径一般在0.5-2mm之间。颗粒之间以孔隙式胶结为主，胶结物主要为泥质、钙质或硅质。这种胶结方式使得砂岩具有一定的孔隙度，其孔隙度一般在10%-20%之间，从而具有较好的透水性能。从构造上看，白垩系中砂岩多呈层状构造，层理较为清晰，这是在沉积过程中由于不同时期的沉积环境差异所导致的。层理的存在使得砂岩在力学性质上表现出各向异性，平行于层理方向的抗压强度和抗拉强度与垂直于层理方向有所不同。侏罗系砂质泥岩的矿物成分则以黏土矿物为主，含量可达50%-70%，常见的黏土矿物有蒙脱石、伊利石、高岭石等。此外，还含有一定量的石英、长石等碎屑矿物，含量约为30%-50%。其结构为泥质结构，颗粒细小，质地细腻。由于黏土矿物的存在，砂质泥岩具有较强的吸水性和可塑性。在结构构造方面，砂质泥岩多呈薄层状或页片状构造，层理发育，页理面较为光滑。这种构造特点使得砂质泥岩在受力时容易沿层理面发生滑动或分离，从而降低其整体强度。在物理力学性质方面，白垩系中砂岩具有较高的强度和弹性模量。其单轴抗压强度一般在30-80MPa之间，弹性模量约为10-30GPa，泊松比在0.2-0.3之间。较高的强度和弹性模量使得砂岩在承受荷载时变形较小，具有较好的承载能力。然而，由于其孔隙度较大，在长期的地下水作用下，可能会发生溶蚀等现象，导致强度降低。侏罗系砂质泥岩的强度相对较低，单轴抗压强度通常在5-20MPa之间，弹性模量约为1-5GPa，泊松比在0.3-0.4之间。较低的强度使得砂质泥岩在受到外力作用时容易发生变形和破坏。而且，砂质泥岩的吸水性强，遇水后会发生膨胀，强度进一步降低，这对巷道的稳定性极为不利。例如，当巷道穿过砂质泥岩地层时，如果没有采取有效的防水和支护措施，砂质泥岩遇水膨胀可能会导致巷道变形、底鼓等问题。白垩系中砂岩和侏罗系砂质泥岩的矿物成分、结构构造和物理力学性质存在显著差异。这些差异使得它们在矿山开采过程中表现出不同的力学行为，对巷道的稳定性产生不同程度的影响。因此，在进行巷道支护设计和开采方案制定时，必须充分考虑围岩的岩性特征，采取针对性的措施，以确保矿山的安全生产和巷道的长期稳定。2.3地应力分布特征地应力是存在于地层中的未受工程扰动的天然应力，也称岩体初始应力、绝对应力或原岩应力，其大小和方向对泊江海子矿围岩稳定性有着至关重要的影响。为了准确掌握矿区的地应力分布情况，采用水压致裂法和套钻孔应力解除法进行地应力测量。水压致裂法通过在钻孔中封隔一小段钻孔，然后向封隔段注入高压流体，从而确定原位地应力，该方法测试周期短，不需要岩石力学参数参与计算，适合工程初勘阶段，且可进行大深度测量。套钻孔应力解除法则是根据解除方式和传感器的安装部位，分为探孔应力解除法、孔底应变解除法和孔壁切割解除法，其中孔壁应变法基于岩石各向同性、均质、连续、线弹性的假设，通过孔壁6个以上不同方向的应变值来计算岩体的三维地应力。在矿区不同位置布置多个测量点，共完成了[X]个钻孔的地应力测量工作。测量结果显示，泊江海子矿地应力大小和方向呈现出一定的规律。最大主应力值在[X1]MPa-[X2]MPa之间，最小主应力值在[X3]MPa-[X4]MPa之间，最大主应力与最小主应力的比值在[X5]-[X6]之间。在方向上，最大主应力方向主要集中在[具体方向区间]，与矿区主要构造方向基本一致。地应力的大小和方向对围岩稳定性有着显著影响。当最大主应力与巷道轴向夹角较小时，巷道周边的应力集中程度相对较低，围岩稳定性较好；而当夹角较大时，巷道周边会出现较大的应力集中，容易导致围岩的破坏和变形。在实际工程中，泊江海子矿部分巷道由于地应力作用，出现了顶板下沉、底板鼓起和两帮收敛等变形现象。在一些地应力较大且方向不利的区域，巷道变形速率明显加快，支护结构承受的压力增大，甚至出现了支护结构损坏的情况。这表明地应力的分布特征是影响泊江海子矿围岩稳定性的重要因素之一，在巷道设计和支护方案制定时，必须充分考虑地应力的大小和方向，以确保巷道的长期稳定和安全生产。三、岩石基本力学特性试验研究3.1试验设备与试样制备本次岩石常规压缩试验采用的是YAW-2000型微机控制电液伺服压力试验机，该设备由主机、液压源、控制器、计算机及打印机等部分组成。主机采用四立柱框架式结构，具有较高的刚度和稳定性，能够保证试验过程中荷载的准确施加。液压源采用高压齿轮泵，可提供稳定的液压动力。控制器采用先进的数字信号处理器（DSP）技术，能够精确控制试验过程中的加载速率和荷载大小。该试验机的最大试验力为2000kN，试验力测量范围为4kN-2000kN，示值相对误差不超过±1%，变形测量范围为0-50mm，分辨率为0.001mm，能够满足岩石常规压缩试验的高精度要求。对于单轴压缩试验，试样取自泊江海子矿不同深度和位置的围岩。为保证试验结果的准确性和代表性，在取样时尽量选取完整、无明显裂隙和风化的岩石。首先，使用金刚石锯将大块岩石切割成尺寸大致为直径50mm、高度100mm的圆柱体粗坯。然后，利用磨石机对粗坯进行精细加工，使试样两端面的平行度偏差不超过0.05mm，直径偏差不超过0.2mm，两端面与轴线的垂直度偏差不超过0.25°。在加工过程中，通过游标卡尺和直角尺对试样的尺寸和垂直度进行实时测量和调整，确保试样符合标准要求。加工完成后的试样表面应光滑、平整，无明显的加工痕迹和损伤。三轴压缩试验的试样制备过程与单轴压缩试验类似，但对试样的密封性和防水性有更高要求。在试样加工完成后，将其套上特制的橡胶膜，橡胶膜应具有良好的弹性和密封性，能够有效防止压力室内的液体渗入试样内部。使用热缩套将橡胶膜紧密固定在试样表面，确保在试验过程中橡胶膜与试样之间不会发生相对滑动。在橡胶膜的两端，使用密封胶进行密封处理，进一步提高试样的防水性能。同时，在试样的顶部和底部放置透水石，以保证在试验过程中试样能够均匀受力，并便于孔隙水的排出。经过严格的制备过程，共获得单轴压缩试验试样[X]组，每组[X]个；三轴压缩试验试样[X]组，每组[X]个。这些试样将用于后续的基本力学特性试验，以获取泊江海子矿围岩的准确力学参数。3.2岩石常规单轴压缩试验3.2.1试验原理与步骤单轴压缩试验是研究岩石力学特性的基本方法之一，其原理基于材料力学中的轴向拉伸与压缩理论。在单轴压缩试验中，将圆柱形岩石试样放置在压力试验机的上下承压板之间，通过试验机以恒定的加载速率对试样施加轴向压力，使试样在无侧向约束的条件下发生轴向压缩变形。根据胡克定律，在弹性阶段，岩石的应力与应变呈线性关系，其比例系数即为弹性模量。当应力达到一定值时，岩石内部开始出现微裂纹，随着应力的继续增加，微裂纹逐渐扩展、贯通，最终导致岩石试样的破坏。在试验开始前，首先对试样进行外观检查，确保试样无明显缺陷。使用游标卡尺在试样的不同部位测量其直径，取平均值作为试样的直径，并测量试样的高度，精确到0.01mm。将测量后的试样放置在压力试验机的下承压板中心位置，调整试验机，使上承压板与试样接触，并确保试样与承压板之间均匀接触，避免偏心受压。试验过程中，按照国家标准《工程岩体试验方法标准》（GB/T50266-2013）的规定，采用位移控制方式，以0.05mm/min-0.1mm/min的加载速率对试样进行加载。在加载过程中，通过试验机自带的数据采集系统实时记录荷载和位移数据，每隔一定时间间隔（如0.1s）采集一次数据，确保数据的连续性和准确性。同时，密切观察试样的变形情况，当试样出现明显的破坏迹象，如出现裂缝、剥落或突然的变形增加时，停止加载，并记录此时的最大荷载。试验结束后，将破坏后的试样从试验机上取出，仔细观察其破坏形态，包括裂缝的方向、数量和分布情况等，并拍照记录。对试验数据进行整理和分析，根据记录的荷载和位移数据，计算试样在不同阶段的应力和应变值，绘制应力-应变曲线。3.2.2试验结果分析通过对单轴压缩试验数据的处理，得到了泊江海子矿围岩的应力-应变曲线，典型的应力-应变曲线如图1所示。从曲线中可以看出，岩石的变形过程可分为以下几个阶段：在OA段，曲线呈近似直线，应力与应变呈线性关系，岩石表现出弹性变形特征，符合胡克定律。此阶段岩石内部的微裂纹尚未开始扩展，变形主要是由于岩石颗粒之间的弹性压缩和晶格的弹性变形引起的。该阶段的斜率即为岩石的弹性模量，通过计算该阶段应力增量与应变增量的比值，可得到弹性模量E。随着应力的增加，进入AB段，曲线开始偏离直线，应力-应变关系呈现非线性，岩石的变形逐渐进入弹塑性阶段。在这个阶段，岩石内部开始出现微裂纹，随着应力的增大，微裂纹不断扩展和连通，导致岩石的变形不再完全可逆，塑性变形逐渐增加。当应力达到B点时，岩石内部的微裂纹大量扩展并相互贯通，形成宏观的破裂面，岩石开始进入破坏阶段。BC段曲线斜率急剧减小，表明岩石的承载能力迅速下降，变形急剧增加，最终导致岩石的完全破坏。B点所对应的应力即为岩石的单轴抗压强度σc，它是衡量岩石强度的重要指标。在岩石达到峰值强度后，虽然其承载能力迅速下降，但仍能承受一定的荷载，曲线进入CD段。此阶段岩石已经破坏，但由于岩石碎块之间的相互摩擦和咬合作用，仍具有一定的残余强度。根据试验数据，计算得到泊江海子矿围岩的弹性模量、泊松比和抗压强度等力学参数，具体结果如表1所示。从表中数据可以看出，不同试样的力学参数存在一定的差异，这主要是由于岩石的非均质性以及试样制备过程中的微小差异导致的。通过对多组试样数据的统计分析，得到了泊江海子矿围岩力学参数的平均值，这些参数将为后续的研究和工程应用提供重要依据。试样编号弹性模量E(GPa)泊松比ν抗压强度σc(MPa)1[X1][X2][X3]2[X4][X5][X6]3[X7][X8][X9]平均值[X10][X11][X12]表1泊江海子矿围岩单轴压缩试验力学参数泊江海子矿围岩在单轴压缩条件下表现出明显的弹塑性变形特征，其弹性模量、泊松比和抗压强度等力学参数对于深入了解围岩的力学行为和工程稳定性分析具有重要意义。3.3岩石常规三轴压缩试验3.3.1试验原理与步骤三轴压缩试验是研究岩石在三维应力状态下力学特性的重要手段，其原理基于摩尔-库仑强度理论。在三轴压缩试验中，将圆柱形岩石试样放置在三轴压力室内，通过压力室中的液体对试样施加恒定的围压σ₃，模拟岩石在地下所受到的侧向压力。然后，通过活塞杆对试样施加轴向压力，使试样在围压和轴压的共同作用下发生变形直至破坏。随着轴向压力的增加，当试样达到极限平衡状态时，其内部某一平面上的剪应力达到岩石的抗剪强度，从而导致岩石发生剪切破坏。试验开始前，首先对试样进行外观检查和尺寸测量，确保试样符合要求。将准备好的试样套上橡胶膜，放入三轴压力室的底座上，安装好压力室的上盖，确保密封良好。连接好围压系统、轴压系统和数据采集系统，检查各系统是否正常工作。试验过程中，先通过围压系统对试样施加预定的围压σ₃，保持围压恒定。然后，按照一定的加载速率，通过轴压系统对试样施加轴向压力σ₁，使试样在三维应力状态下发生变形。加载速率通常根据岩石的类型和试验要求进行选择，一般为0.05mm/min-0.1mm/min。在加载过程中，利用数据采集系统实时记录轴向荷载、轴向位移、围压以及孔隙水压力等数据，每隔一定时间间隔采集一次数据，以获取岩石在加载过程中的应力-应变关系。同时，密切观察试样的变形和破坏情况，当试样出现明显的破坏迹象，如出现裂缝、剥落或突然的变形增加时，停止加载，并记录此时的轴向荷载和轴向位移。试验结束后，关闭围压系统和轴压系统，拆卸压力室，取出破坏后的试样，观察其破坏形态，包括裂缝的方向、数量和分布情况等，并拍照记录。对试验数据进行整理和分析，根据记录的荷载和位移数据，计算试样在不同阶段的应力和应变值，绘制应力-应变曲线。3.3.2试验结果分析通过对三轴压缩试验数据的处理，得到了泊江海子矿围岩在不同围压下的应力-应变曲线，典型的应力-应变曲线如图2所示。从曲线中可以看出，随着围压的增加，岩石的变形特性和强度发生了显著变化。在弹性阶段，应力-应变曲线近似为直线，岩石表现出弹性变形特征。围压的增加使得岩石的弹性阶段延长，弹性模量增大，这表明围压能够增强岩石的抵抗变形能力。在弹塑性阶段，随着应力的增加，岩石内部开始出现微裂纹，塑性变形逐渐增加。围压的增大使得岩石的塑性变形阶段更加明显，岩石的韧性增强，能够承受更大的变形而不发生破坏。当应力达到峰值强度时，岩石内部的微裂纹大量扩展并相互贯通，形成宏观的破裂面，岩石开始进入破坏阶段。围压的增加显著提高了岩石的峰值强度，即岩石的抗压强度随着围压的增大而增大。这是因为围压的存在限制了岩石内部微裂纹的扩展和贯通，使得岩石需要更大的应力才能发生破坏。在岩石达到峰值强度后，进入残余强度阶段。围压的增加使得岩石的残余强度也有所提高，这是因为围压能够抑制岩石碎块之间的相对滑动，使得岩石在破坏后仍能保持一定的承载能力。为了更直观地分析围压对岩石力学性质的影响，对不同围压下的岩石力学参数进行了统计分析，结果如表2所示。从表中数据可以看出，随着围压的增加，岩石的弹性模量、抗压强度和残余强度均呈现增大的趋势，而泊松比则略有减小。这表明围压对岩石的力学性质有着显著的影响，在工程设计和分析中，必须充分考虑围压的作用。围压σ₃(MPa)弹性模量E(GPa)泊松比ν抗压强度σc(MPa)残余强度σr(MPa)5[X1][X2][X3][X4]10[X5][X6][X7][X8]15[X9][X10][X11][X12]表2不同围压下泊江海子矿围岩力学参数泊江海子矿围岩在三轴压缩条件下，围压对其力学性质有着显著的影响。随着围压的增加，岩石的强度和抵抗变形能力增强，韧性提高。这些试验结果为深入了解围岩在复杂应力条件下的力学行为提供了重要依据，对于泊江海子矿的巷道支护设计和开采方案优化具有重要的指导意义。四、岩石蠕变特性试验研究4.1试验方案设计蠕变试验选用TAW-2000型微机控制电液伺服岩石三轴试验机，该设备主要由主机、液压源、控制系统、数据采集系统等部分组成。主机采用高强度框架结构，具备良好的稳定性和刚度，能够确保在试验过程中精确施加荷载并保持稳定。液压源能够提供稳定的液压动力，满足试验所需的加载要求。控制系统采用先进的微处理器技术，可实现对试验过程的精确控制，包括加载速率、加载方式、荷载保持时间等参数的设定。数据采集系统则能够实时采集试验过程中的应力、应变、时间等数据，采集频率可根据试验需求进行调整，确保数据的准确性和完整性。该试验机可实现最大试验力2000kN，轴向位移测量精度为±0.001mm，围压控制精度为±0.1MPa，能够满足岩石单轴和三轴蠕变试验的高精度要求。在单轴蠕变试验中，采用分级加载方式，即先对试样施加一个初始应力水平，待试样的蠕变变形基本稳定后，再施加下一级应力，直至试样破坏。根据前期常规单轴压缩试验得到的岩石单轴抗压强度，确定应力水平分别为单轴抗压强度的30%、40%、50%、60%、70%。这样设置应力水平的目的是全面研究岩石在不同应力条件下的蠕变特性，从较低应力水平下的稳定蠕变阶段，到较高应力水平下逐渐进入加速蠕变阶段，从而更深入地了解岩石蠕变的全过程。试验时间设定为每个应力水平下持续加载24小时，以确保能够充分观测到岩石在该应力水平下的蠕变变形规律。若在24小时内试样出现明显的加速蠕变迹象或破坏，则提前终止试验并记录相关数据。三轴蠕变试验同样采用分级加载方式，先施加预定的围压，保持围压恒定后，再分级施加轴向应力。围压水平根据矿区的实际地应力情况确定，分别选取5MPa、10MPa、15MPa，以模拟不同深度的地应力环境。轴向应力水平则根据常规三轴压缩试验结果，分别设定为不同围压下岩石抗压强度的30%、40%、50%、60%。这种设置方式能够综合考虑围压和轴向应力对岩石蠕变特性的影响，更真实地反映岩石在地下复杂应力状态下的蠕变行为。试验时间在每个应力水平下同样持续加载24小时，密切观察试样在不同应力组合下的蠕变变形随时间的变化情况。如果在加载过程中试样出现异常变形或破坏，及时停止试验并记录相关数据。在试验过程中，使用高精度的位移传感器测量试样的轴向应变和径向应变，位移传感器的精度可达±0.001mm，能够精确捕捉试样在蠕变过程中的微小变形。通过数据采集系统，以每秒1次的频率实时采集并记录应变和时间数据，确保数据的连续性和完整性。同时，密切观察试样的外观变化，如是否出现裂缝、剥落等现象，并及时拍照记录，以便后续对试样的破坏形式和破坏机理进行分析。4.2岩石单轴蠕变试验4.2.1试验过程与数据采集在正式开展单轴蠕变试验前，需对TAW-2000型微机控制电液伺服岩石三轴试验机进行全面检查与调试，确保设备的各项性能指标满足试验要求。仔细检查主机的框架结构是否稳固，液压源的油压是否正常，控制系统的参数设置是否准确，数据采集系统的传感器是否灵敏等。同时，对试验所需的各种辅助设备和工具，如位移传感器、应变片、导线等进行检查和校准，确保其精度和可靠性。将制备好的岩石试样小心放置在试验机的下承压板中心位置，确保试样与承压板紧密接触，且受力均匀。安装好位移传感器，使其测量端与试样的轴向方向垂直，并与试样表面紧密接触，以准确测量试样的轴向变形。连接好数据采集系统的导线，确保信号传输稳定。根据预先设计的试验方案，采用分级加载方式对试样施加轴向应力。首先，通过控制系统设置初始应力水平为单轴抗压强度的30%，然后以缓慢且稳定的速率逐渐增加轴向应力，直至达到预定的初始应力值。在加载过程中，密切关注试验机的运行状态和试样的变形情况，确保加载过程平稳，无异常波动。当应力达到初始应力水平后，立即启动数据采集系统，开始记录试样的轴向应变随时间的变化数据。数据采集系统以每秒1次的频率自动采集并存储应变和时间数据，确保数据的连续性和准确性。在每个应力水平下，持续加载24小时，期间实时监测试样的变形情况。若在24小时内试样出现明显的加速蠕变迹象，如应变随时间急剧增加，或者出现破坏现象，如试样表面出现裂缝、剥落等，则立即停止加载，并记录此时的时间和应变数据。当一个应力水平下的试验结束后，按照同样的方法，逐渐增加轴向应力至下一个预定的应力水平，即单轴抗压强度的40%，重复上述加载、监测和数据采集过程。依次类推，直至完成所有预定应力水平下的试验。在整个试验过程中，除了采集应变和时间数据外，还密切观察试样的外观变化。每隔一定时间间隔，对试样进行拍照记录，以便后续分析试样在不同应力水平下的变形特征和破坏形式。同时，详细记录试验过程中的各种异常情况，如试验机故障、传感器失灵等，以便对试验数据进行合理的分析和处理。4.2.2试验结果分析通过对单轴蠕变试验数据的深入分析，得到了泊江海子矿围岩在不同应力水平下的典型蠕变曲线，如图3所示。从曲线中可以清晰地看出，岩石的蠕变过程可分为三个明显的阶段。在OA段，即初始蠕变阶段，也称为减速蠕变阶段。在这一阶段，岩石的应变随时间迅速增加，但应变增加的速率逐渐减小，曲线斜率逐渐变缓。这是因为在初始加载时，岩石内部存在较多的微裂纹和孔隙，随着应力的作用，这些微裂纹和孔隙逐渐被压实和闭合，岩石的结构逐渐趋于紧密，从而导致应变增加的速率逐渐减小。在这个阶段，岩石的变形主要是由弹性变形和少量的塑性变形组成，弹性变形占主导地位。随着时间的推移，进入AB段，即等速蠕变阶段，也称为稳态蠕变阶段。在该阶段，岩石的应变随时间以较为稳定的速率增加，蠕变曲线近似为一条直线，斜率保持相对稳定。此时，岩石内部的微裂纹和孔隙已经基本压实，岩石的结构相对稳定，变形主要是由于岩石内部的位错运动和晶界滑移等机制引起的，塑性变形逐渐占据主导地位。在等速蠕变阶段，岩石的蠕变速率保持相对恒定，这一阶段的持续时间与岩石的性质、应力水平等因素有关。当应力达到一定程度后，岩石进入BC段，即加速蠕变阶段。在这一阶段，岩石的应变随时间急剧增加，蠕变曲线斜率迅速增大，变形速率显著加快。这是因为随着应力的持续作用和时间的延长，岩石内部的微裂纹不断扩展、贯通，形成宏观的裂缝，岩石的结构逐渐破坏，承载能力急剧下降，最终导致岩石的破坏。在加速蠕变阶段，岩石的变形主要是由塑性变形和断裂变形组成，岩石的破坏过程迅速发展。为了更深入地研究应力水平对蠕变特性的影响，对不同应力水平下的蠕变速率和蠕变极限进行了详细分析。蠕变速率是指单位时间内岩石的应变变化量，它反映了岩石蠕变的快慢程度。通过对试验数据的计算和分析，得到了不同应力水平下的蠕变速率随时间的变化曲线，如图4所示。从图中可以看出，随着应力水平的增加，蠕变速率明显增大。在较低的应力水平下，如单轴抗压强度的30%和40%，蠕变速率相对较小，且在等速蠕变阶段保持相对稳定。然而，当应力水平提高到单轴抗压强度的60%和70%时，蠕变速率显著增大，且在加速蠕变阶段增加的速度更快。这表明应力水平越高，岩石内部的损伤发展越快，蠕变变形越剧烈。蠕变极限是指岩石在一定时间内所能承受的最大应变值，它反映了岩石的长期承载能力。通过对试验数据的分析，得到了不同应力水平下的蠕变极限，如表3所示。从表中数据可以看出，随着应力水平的增加，蠕变极限逐渐减小。在较低的应力水平下，岩石能够承受较大的应变而不发生破坏，表明其长期承载能力较强。然而，当应力水平超过一定限度时，岩石的蠕变极限急剧减小，表明其长期承载能力迅速下降，容易发生破坏。应力水平（%）蠕变极限（με）30[X1]40[X2]50[X3]60[X4]70[X5]表3不同应力水平下的蠕变极限泊江海子矿围岩在单轴蠕变试验中表现出典型的蠕变特性，不同应力水平对蠕变特性有着显著的影响。随着应力水平的增加，蠕变速率增大，蠕变极限减小，岩石的长期稳定性降低。这些试验结果为深入理解泊江海子矿围岩的蠕变行为提供了重要依据，对于该矿的巷道支护设计和开采方案优化具有重要的指导意义。4.3岩石三轴蠕变试验4.3.1试验过程与数据采集在开展三轴蠕变试验之前，对TAW-2000型微机控制电液伺服岩石三轴试验机进行全面检查和调试，确保设备的各项性能指标满足试验要求。仔细检查主机的框架结构是否稳固，液压源的油压是否稳定，控制系统的参数设置是否准确，数据采集系统的传感器是否灵敏等。同时，对试验所需的各种辅助设备和工具，如位移传感器、应变片、导线、橡胶膜、热缩套、透水石等进行检查和校准，确保其精度和可靠性。将制备好的岩石试样套上橡胶膜，使用热缩套紧密固定，确保橡胶膜与试样之间密封良好，防止压力室内的液体渗入试样内部。在试样的顶部和底部放置透水石，以保证在试验过程中试样能够均匀受力，并便于孔隙水的排出。将处理好的试样小心放置在三轴压力室的底座上，安装好压力室的上盖，确保密封良好。连接好围压系统、轴压系统和数据采集系统，检查各系统的连接是否牢固，信号传输是否正常。根据预先设计的试验方案，先通过围压系统对试样施加预定的围压，围压水平分别选取5MPa、10MPa、15MPa，以模拟不同深度的地应力环境。在施加围压时，采用缓慢、稳定的加载方式，加载速率控制在0.1MPa/min左右，以避免因加载过快对试样造成损伤。当围压达到预定值后，保持围压恒定，持续观察围压的稳定性，确保围压在试验过程中波动不超过±0.1MPa。在围压保持稳定后，按照分级加载方式对试样施加轴向应力。轴向应力水平根据常规三轴压缩试验结果，分别设定为不同围压下岩石抗压强度的30%、40%、50%、60%。在施加轴向应力时，同样采用缓慢、稳定的加载方式，加载速率控制在0.05mm/min-0.1mm/min之间。当轴向应力达到预定的第一级应力水平后，立即启动数据采集系统，开始记录试样的轴向应变、径向应变随时间的变化数据。数据采集系统以每秒1次的频率自动采集并存储应变和时间数据，确保数据的连续性和准确性。在每个应力水平下，持续加载24小时，期间实时监测试样的变形情况。若在24小时内试样出现明显的加速蠕变迹象，如应变随时间急剧增加，或者出现破坏现象，如试样表面出现裂缝、剥落等，则立即停止加载，并记录此时的时间和应变数据。当一个应力水平下的试验结束后，按照同样的方法，逐渐增加轴向应力至下一个预定的应力水平，重复上述加载、监测和数据采集过程。依次类推，直至完成所有预定应力水平下的试验。在整个试验过程中，除了采集应变和时间数据外，还密切观察试样的外观变化。每隔一定时间间隔，对试样进行拍照记录，以便后续分析试样在不同应力组合下的变形特征和破坏形式。同时，详细记录试验过程中的各种异常情况，如试验机故障、传感器失灵、压力室泄漏等，以便对试验数据进行合理的分析和处理。4.3.2试验结果分析通过对三轴蠕变试验数据的深入分析，得到了泊江海子矿围岩在不同围压和轴压组合下的典型蠕变曲线，如图5所示。从曲线中可以看出，三轴蠕变曲线同样呈现出明显的阶段性特征。在初始阶段，即OA段，曲线斜率较大，应变随时间迅速增加，但增加速率逐渐减小，这与单轴蠕变试验中的初始蠕变阶段类似。在这个阶段，岩石内部的微裂纹和孔隙在围压和轴压的共同作用下迅速被压实和闭合，导致应变快速增加，但随着压实和闭合过程的进行，应变增加速率逐渐减缓。此阶段的变形主要由弹性变形和少量的塑性变形组成，弹性变形占主导地位。随着时间的推移，进入AB段，曲线斜率趋于稳定，应变随时间以较为稳定的速率增加，进入等速蠕变阶段。在该阶段，岩石内部结构相对稳定，变形主要是由于岩石内部的位错运动和晶界滑移等机制引起的，塑性变形逐渐占据主导地位。与单轴蠕变试验相比，三轴蠕变试验中的等速蠕变阶段持续时间更长，这是因为围压的存在抑制了微裂纹的扩展，使得岩石的结构更加稳定，能够承受更长时间的稳定变形。当应力达到一定程度后，岩石进入BC段，即加速蠕变阶段。在这一阶段，曲线斜率急剧增大，应变随时间急剧增加，变形速率显著加快。这是由于随着应力的持续作用和时间的延长，岩石内部的微裂纹在围压和轴压的共同作用下不断扩展、贯通，形成宏观的裂缝，岩石的结构逐渐破坏，承载能力急剧下降，最终导致岩石的破坏。与单轴蠕变试验相比，三轴蠕变试验中的加速蠕变阶段发展更为迅速，这是因为围压和轴压的共同作用使得岩石内部的损伤发展更快，破坏过程更为剧烈。为了深入研究围压和轴压对蠕变特性的影响，对不同围压和轴压组合下的蠕变速率进行了详细分析。在低围压（如5MPa）下，随着轴压的增加，蠕变速率逐渐增大。当轴压达到一定程度后，蠕变速率急剧增加，表明岩石进入了加速蠕变阶段。在高围压（如15MPa）下，蠕变速率的变化相对较为平缓，轴压的增加对蠕变速率的影响相对较小。这是因为高围压能够抑制微裂纹的扩展，使得岩石的结构更加稳定，从而减缓了蠕变速率的增加。对比不同岩样的蠕变差异，发现即使在相同的围压和轴压条件下，不同岩样的蠕变曲线也存在一定的差异。这主要是由于岩石的非均质性导致的，不同岩样的矿物成分、结构构造以及内部缺陷分布等存在差异，这些差异会影响岩石的蠕变特性。一些岩样中含有较多的黏土矿物，其蠕变变形相对较大，蠕变速率也相对较高；而一些岩样的结构较为致密，矿物成分相对均匀，其蠕变变形相对较小，蠕变速率也相对较低。泊江海子矿围岩在三轴蠕变试验中表现出与单轴蠕变试验相似但又有差异的蠕变特性。围压和轴压对蠕变特性有着显著的影响，不同岩样之间也存在一定的蠕变差异。这些试验结果为深入理解泊江海子矿围岩在复杂应力条件下的蠕变行为提供了重要依据，对于该矿的巷道支护设计和开采方案优化具有重要的指导意义。4.4影响围岩蠕变的因素分析岩性是影响泊江海子矿围岩蠕变特性的重要内在因素之一。不同岩性的围岩，其矿物成分、结构构造以及微观组织结构存在显著差异，这些差异直接决定了围岩的基本物理力学性质，进而对蠕变特性产生重要影响。泊江海子矿围岩主要由白垩系中砂岩和侏罗系砂质泥岩组成。白垩系中砂岩以石英、长石等矿物为主，颗粒间以孔隙式胶结为主，结构相对致密。这种矿物组成和结构使得砂岩具有较高的强度和较低的蠕变变形能力。在蠕变试验中，白垩系中砂岩的蠕变速率相对较低，蠕变变形较小，且在较高应力水平下才会进入加速蠕变阶段。这是因为其矿物颗粒之间的胶结作用较强，能够有效抵抗外力引起的变形，内部微裂纹的产生和扩展相对缓慢。相比之下，侏罗系砂质泥岩以黏土矿物为主，结构较为疏松，颗粒间的连接较弱。黏土矿物具有较大的吸水性和膨胀性，遇水后会发生膨胀，导致岩石结构破坏，强度降低。这种特性使得砂质泥岩的蠕变变形能力较强，蠕变速率较高。在较低的应力水平下，砂质泥岩就可能进入加速蠕变阶段，且蠕变变形随时间的增长较为明显。在实际工程中，当巷道穿过砂质泥岩地层时，由于砂质泥岩的蠕变特性，容易导致巷道变形、底鼓等问题，严重影响巷道的稳定性。应力水平对泊江海子矿围岩蠕变特性有着显著的影响。随着应力水平的增加，围岩的蠕变速率明显增大，蠕变变形也显著增加。在单轴蠕变试验中，当应力水平从单轴抗压强度的30%增加到70%时，蠕变速率逐渐增大，且在较高应力水平下，蠕变曲线更快地进入加速蠕变阶段。这是因为在较高应力水平下，岩石内部的微裂纹更容易产生和扩展，导致岩石的结构损伤加剧，从而使得蠕变变形加速发展。在三轴蠕变试验中，围压和轴压的共同作用也体现了应力水平的影响。随着轴压和围压的增加，围岩的蠕变特性发生明显变化。在低围压下，轴压的增加对蠕变速率的影响更为显著，当轴压达到一定程度时，蠕变速率急剧增加，岩石迅速进入加速蠕变阶段。而在高围压下，由于围压对微裂纹扩展的抑制作用，轴压的增加对蠕变速率的影响相对较小，岩石的蠕变变形相对较为稳定。应力水平的变化还会影响围岩的长期强度和蠕变极限。较高的应力水平会降低围岩的长期强度和蠕变极限，使得围岩在较短的时间内发生破坏。温度是影响岩石蠕变特性的重要外部因素之一，泊江海子矿围岩也不例外。一般来说，随着温度的升高，岩石的蠕变速率增大，蠕变变形增加。这是因为温度升高会使岩石内部的分子热运动加剧，降低岩石内部颗粒之间的相互作用力，从而使得岩石更容易发生变形。在高温环境下，岩石内部的位错运动和晶界滑移等机制更容易发生，导致岩石的蠕变变形加速。虽然泊江海子矿的开采深度相对较浅，地温变化相对较小，但在一些特殊区域，如靠近热源（如深部岩浆活动区域或地下热水区域）或在开采过程中由于机械摩擦等原因产生热量的区域，温度对围岩蠕变特性的影响仍不可忽视。当温度升高时，围岩的蠕变速率可能会显著增加，从而影响巷道的稳定性。在深部开采区域，如果地温较高，围岩的蠕变变形可能会更加明显，对巷道支护结构的压力也会增大，需要采取相应的措施来应对温度对围岩蠕变特性的影响。湿度对泊江海子矿围岩蠕变特性的影响主要体现在对岩石物理性质的改变上。泊江海子矿围岩中的砂质泥岩等岩石具有较强的吸水性，当湿度增加时，岩石中的黏土矿物会吸收水分，发生膨胀，导致岩石结构破坏，强度降低。这种强度降低会使得岩石的蠕变变形能力增强，蠕变速率增大。在实际工程中，地下水的存在会使围岩处于高湿度环境中，从而对围岩的蠕变特性产生显著影响。当巷道周围存在地下水时，地下水会渗入围岩，使围岩的湿度增加，进而加速围岩的蠕变变形。在一些富水区域，巷道的变形速率明显加快，这与湿度对围岩蠕变特性的影响密切相关。湿度的变化还可能导致岩石内部产生干湿循环，进一步加剧岩石的结构破坏和强度降低，从而对围岩的蠕变特性产生更为复杂的影响。五、岩石蠕变本构关系研究5.1岩石蠕变理论基础岩石蠕变是指岩石在恒定应力作用下，其变形随时间不断发展的现象。这一现象在自然界和工程实践中广泛存在，对岩石工程的长期稳定性具有重要影响。岩石蠕变曲线通常呈现出典型的阶段性特征，清晰地反映了岩石在不同阶段的变形行为和内部结构变化。在初始阶段，也称为减速蠕变阶段。当应力施加于岩石时，岩石内部的微结构开始调整。由于岩石内部存在各种缺陷和孔隙，在应力作用下，这些缺陷和孔隙首先被压缩和闭合，导致应变迅速增加。随着时间的推移，微结构的调整逐渐趋于稳定，应变增加的速率逐渐减小，曲线斜率逐渐变缓。在这个阶段，岩石的变形主要是由弹性变形和少量的塑性变形组成，弹性变形占主导地位。随着时间的进一步推移，岩石进入等速蠕变阶段，也称为稳态蠕变阶段。在该阶段，岩石内部的微结构已经基本调整完毕，变形主要是由于岩石内部的位错运动和晶界滑移等机制引起的。这些微观机制的作用使得岩石的应变随时间以较为稳定的速率增加，蠕变曲线近似为一条直线，斜率保持相对稳定。在等速蠕变阶段，岩石的蠕变速率保持相对恒定，这一阶段的持续时间与岩石的性质、应力水平等因素有关。当应力达到一定程度后，岩石进入加速蠕变阶段。在这一阶段，岩石内部的微裂纹不断扩展、贯通，形成宏观的裂缝，岩石的结构逐渐破坏，承载能力急剧下降。随着裂缝的不断发展，岩石的变形速率显著加快，蠕变曲线斜率迅速增大，最终导致岩石的破坏。在加速蠕变阶段，岩石的变形主要是由塑性变形和断裂变形组成，岩石的破坏过程迅速发展。岩石蠕变的机理较为复杂，涉及多个方面的因素。从微观角度来看，岩石是由各种矿物颗粒组成，矿物颗粒之间通过晶界相互连接。在应力作用下，矿物颗粒内部会发生位错运动，位错的滑移和攀移导致矿物颗粒的变形。晶界处的原子排列不规则，能量较高，在应力作用下容易发生晶界扩散和滑动，从而导致岩石的变形。岩石内部还存在各种缺陷，如孔隙、微裂纹等，这些缺陷在应力作用下会不断扩展和连通，进一步加剧岩石的变形。岩石的化学成分和矿物组成对蠕变也有重要影响。不同的矿物具有不同的晶体结构和力学性质，其蠕变特性也各不相同。含有较多黏土矿物的岩石，由于黏土矿物的吸水性和膨胀性，在水分作用下容易发生膨胀和软化，从而导致岩石的蠕变变形增加。而含有较多石英、长石等矿物的岩石，其结构相对稳定，蠕变变形相对较小。岩石的微观结构，如孔隙率、孔隙分布、颗粒大小和形状等，也会影响岩石的蠕变性能。孔隙率较大的岩石，其内部缺陷较多，在应力作用下更容易发生变形，蠕变速率也相对较高。5.2现有蠕变本构模型分析经验流变模型主要基于试验数据，通过数学拟合的方式建立应力、应变与时间之间的关系。这类模型的建立过程相对简单，直接依据试验所获取的蠕变曲线，运用数学函数进行拟合。例如，幂函数模型\varepsilon(t)=\varepsilon_0+At^n，其中\varepsilon(t)为t时刻的应变，\varepsilon_0为瞬时应变，A和n为与岩石性质及应力水平相关的常数。该模型能够较好地描述岩石在初始蠕变阶段和稳态蠕变阶段的部分特性，通过调整参数A和n，可以使模型曲线与试验数据在这两个阶段实现较好的拟合。经验流变模型的优点在于其建立过程简单直接，对试验数据的依赖程度高，能够快速地根据试验结果得到一个初步的模型，适用于对模型精度要求不是特别高，且试验数据较为充足的情况。在一些工程初步设计阶段，需要快速了解岩石的大致蠕变特性，经验流变模型可以提供一个较为便捷的参考。然而，经验流变模型缺乏明确的物理意义，它只是对试验数据的一种数学描述，无法从微观机制上解释岩石的蠕变行为。这使得在实际应用中，对于模型参数的变化难以从物理本质上进行理解和解释，模型的外推能力较差，当应用于与试验条件差异较大的情况时，其预测的准确性难以保证。组合流变模型是将弹性元件（如弹簧，符合胡克定律，代表弹性变形）、粘性元件（如阻尼器，符合牛顿粘性定律，代表粘性变形）和塑性元件（如摩擦片，代表塑性变形）通过串联、并联或串并联等方式组合而成，以模拟岩石的复杂流变行为。马克斯威尔（Maxwell）体由一个弹性元件和一个粘性元件串联而成，其本构方程为\dot{\varepsilon}=\frac{\dot{\sigma}}{E}+\frac{\sigma}{\eta}，其中\dot{\varepsilon}为应变速率，\dot{\sigma}为应力速率，E为弹性模量，\eta为粘性系数。该模型具有瞬时弹性变形和等速蠕变特性，能够描述岩石在恒定应力作用下，先产生瞬时弹性变形，随后变形随时间以恒定速率增加的蠕变行为。开尔文（Kelvin）体则是由一个弹性元件和一个粘性元件并联而成，本构方程为\sigma=E\varepsilon+\eta\dot{\varepsilon}，它具有弹性后效特性，无松弛现象，在恒定应力作用下，应变随时间逐渐增加，卸载后应变逐渐恢复，适用于描述具有衰减蠕变和弹性后效的岩石。组合流变模型的优点是能够通过不同元件的组合，较为直观地反映岩石的弹性、塑性和粘性等多种力学特性，物理意义明确。不同的元件组合可以模拟出不同类型岩石在不同应力条件下的流变行为，具有较强的适应性。然而，组合流变模型也存在一定的局限性。对于复杂的岩石流变行为，需要大量的元件进行组合，这会导致模型结构复杂，参数众多，增加了参数辨识的难度和不确定性。在实际应用中，过多的参数也使得模型的计算量增大，计算效率降低，并且对于一些特殊的岩石流变现象，如岩石的损伤演化等，传统的组合流变模型难以准确描述。改进流变模型是在传统流变模型的基础上，针对其局限性进行改进和优化。一些改进流变模型通过引入新的流变元件或对传统元件进行修正，来提高模型对岩石复杂流变行为的描述能力。在传统的西原模型中，为了更好地描述岩石在加速蠕变阶段的非线性特性，引入了非线性粘壶元件。传统西原模型由虎克体、黏弹性体以及黏塑性体串并联组合而成，难以描述岩石非线性加速蠕变阶段的流变特性。通过在西原模型上串联一个带应变触发的非线性黏壶，提出改进的西原模型，推导岩石在恒应力情况下的三维蠕变本构方程。改进后的模型不仅可以充分反映岩石初期蠕变、稳定蠕变阶段的流变特性，还可以很好地描述岩石加速蠕变阶段的蠕变规律，明显优于传统的西原模型。改进流变模型的优点是能够针对传统模型的不足，通过合理的改进措施，提高模型对岩石复杂流变行为的模拟精度，尤其是对于一些传统模型难以描述的特殊流变现象，如岩石的非线性加速蠕变等，改进模型能够提供更准确的描述。然而，改进流变模型的建立需要对岩石的流变特性有深入的理解和研究，改进措施的合理性需要通过大量的试验和理论分析来验证。一些改进模型可能会因为引入新的参数或元件，导致模型的复杂性增加，参数的物理意义不够明确，在实际应用中需要谨慎对待。引入损伤理论的流变模型将损伤力学与流变学相结合，考虑岩石在蠕变过程中的内部损伤演化对其力学行为的影响。这类模型认为，岩石在蠕变过程中，内部会逐渐产生微裂纹、孔隙等损伤，这些损伤的发展会导致岩石的力学性能劣化，从而影响其蠕变特性。基于损伤理论的流变模型通常通过引入损伤变量来描述岩石的损伤程度，损伤变量可以与岩石的微裂纹密度、孔隙率等微观结构参数相关联。在模型中，损伤变量会影响流变模型中的参数，如弹性模量、粘性系数等，使得这些参数随着损伤的发展而发生变化，从而更准确地描述岩石的蠕变行为。引入损伤理论的流变模型的优点是能够从微观机制上解释岩石蠕变过程中的力学性能劣化现象，考虑了岩石内部结构的变化对蠕变的影响，使模型更加符合岩石的实际力学行为。对于研究岩石在长期荷载作用下的渐进破坏过程，这类模型具有重要的意义。然而，损伤理论本身较为复杂，损伤变量的定义和测量方法尚未统一，不同的损伤理论和损伤变量定义会导致模型的差异较大。在实际应用中，损伤变量的获取往往需要借助微观测试技术，这增加了模型应用的难度和成本。而且，将损伤理论与流变模型相结合时，模型的数学推导和求解过程也会变得更加复杂，对计算能力和理论水平提出了更高的要求。5.3基于试验数据的本构模型建立5.3.1模型选择与建立在众多流变模型中，西原模型由于其能够较好地描述岩石的弹性、塑性和粘性等特性，在岩石蠕变研究中得到了广泛应用。传统西原模型由虎克体、黏弹性体以及黏塑性体串并联组合而成，然而，对于泊江海子矿围岩复杂的蠕变特性，传统西原模型存在一定的局限性，难以准确描述其非线性加速蠕变阶段的流变特性。为了更精确地描述泊江海子矿围岩的蠕变行为，对传统西原模型进行改进。在传统西原模型的基础上，串联一个带应变触发的非线性黏壶，提出改进的西原模型。改进后的模型结构如图6所示，其中，H1为弹性元件（弹簧），代表岩石的瞬时弹性变形，其应力-应变关系满足胡克定律，即\sigma=E_1\varepsilon，E_1为弹性模量；K为开尔文体，由弹性元件H2和粘性元件N1并联组成，代表岩石的黏弹性变形，其本构方程为\sigma=E_2\varepsilon+\eta_1\dot{\varepsilon}，E_2为弹性模量，\eta_1为粘性系数；B为宾汉体，由弹性元件H3、塑性元件Y和粘性元件N2串联组成，代表岩石的黏塑性变形，当应力超过屈服应力\sigma_s时，宾汉体开始产生塑性变形，其本构方程为\sigma=E_3\varepsilon+\sigma_s+\eta_2\dot{\varepsilon}，E_3为弹性模量，\eta_2为粘性系数；新增的非线性黏壶N3用于描述岩石在加速蠕变阶段的非线性特性，其应变速率与应力之间的关系采用非线性函数表示，如\dot{\varepsilon}_3=A\sigma^n，A和n为与岩石性质相关的常数。通过这种改进，新模型不仅能够反映岩石在初始蠕变和稳态蠕变阶段的流变特性，还能准确描述岩石在加速蠕变阶段的非线性变形行为，更符合泊江海子矿围岩的实际蠕变特性。5.3.2模型参数辨识与验证采用非线性拟合方法对改进西原模型的参数进行辨识。将蠕变试验数据作为输入，利用最小二乘法等优化算法，通过不断调整模型中的参数，使模型计算得到的应变值与试验测量的应变值之间的误差最小化。具体过程如下：建立目标函数：定义目标函数F(\theta)，其中\theta=[E_1,E_2,\eta_1,E_3,\sigma_s,\eta_2,A,n]为模型参数向量，目标函数F(\theta)表示模型计算应变\varepsilon_{cal}(\theta)与试验测量应变\varepsilon_{exp}之间的误差平方和，即F(\theta)=\sum_{i=1}^{m}(\varepsilon_{cal}(\theta)_i-\varepsilon_{exp}_i)^2，m为试验数据点的数量。优化算法求解：采用优化算法，如Levenberg-Marquardt算法，对目标函数F(\theta)进行求解。该算法结合了梯度下降法和高斯-牛顿法的优点，能够在保证收敛速度的同时，提高求解的稳定性。在求解过程中，不断调整参数向量\theta的值，使目标函数F(\theta)逐渐减小，直至达到预设的收敛条件。确定参数值：当目标函数F(\theta)收敛到最小值时，此时的参数向量\theta即为辨识得到的模型参数值。通过上述方法，得到了改进西原模型在不同应力水平下的参数值，如表4所示。应力水平（%）E1(MPa)E2(MPa)η1(MPa·s)E3(MPa)σs(MPa)η2(MPa·s)An30[X1][X2][X3][X4][X5][X6][X7][X8]40[X9][X10][X11][X12][X13][X14][X15][X16]50[X17][X18][X19][X20][X21][X22][X23][X24]60[X25][X26][X27][X28][X29][X30][X31][X32]表4不同应力水平下改进西原模型参数值为了验证改进西原模型的准确性和可靠性，将模型计算结果与试验数据进行对比。以应力水平为50%的蠕变试验数据为例，对比结果如图7所示。从图中可以看出，改进西原模型的计算结果与试验数据具有良好的一致性，能够准确地描述泊江海子矿围岩在该应力水平下的蠕变特性。在初始蠕变阶段，模型计算的应变随时间迅速增加，与试验数据相符；在稳态蠕变阶段，模型计算的应变以较为稳定的速率增加，与试验数据的变化趋势一致；在加速蠕变阶段，模型能够准确地捕捉到应变随时间急剧增加的趋势，计算结果与试验数据的偏差较小。通过对不同应力水平下的试验数据进行验证，结果表明改进西原模型能够较好地描述泊江海子矿围岩的蠕变特性，模型参数辨识准确，具有较高的准确性和可靠性，为泊江海子矿巷道稳定性分析和支护设计提供了可靠的理论依据。六、数值模拟与工程应用6.1基于本构模型的数值模拟利用有限元软件ANSYS对泊江海子矿巷道开挖过程进行数值模拟，将前文建立的改进西原蠕变本构模型应用于模拟分析中。该软件具有强大的前处理、求解和后处理功能，能够方便地建立复杂的几何模型、划分网格、施加边界条件和荷载，并对模拟结果进行直观的可视化展示。首先，根据泊江海子矿的实际地质条件和巷道设计参数，建立三维数值模型。模型尺寸根据巷道的影响范围确定，一般取巷道周边一定距离，以保证模型