基于煤岩蠕变试验的深井条带煤柱长期稳定性解析与展望

基于煤岩蠕变试验的深井条带煤柱长期稳定性解析与展望一、引言1.1研究背景与意义煤炭在我国能源结构中占据着举足轻重的地位，是我国的主体能源。我国“富煤、贫油、少气”的资源禀赋特点，决定了煤炭在未来相当长的一段时间内仍将是保障国家能源安全和经济社会稳定发展的关键能源。从储量来看，我国已探明的煤炭储量在世界煤炭储量中占比可观，可采量位居世界前列，产量更是长期稳居世界第一。在一次能源消费构成中，煤炭一直占据着较高的比例，尽管近年来随着新能源和可再生能源的快速发展，煤炭消费比重有所下降，但在2023年仍达到55.3%。煤炭不仅广泛应用于电力、钢铁、化工等传统行业，为国民经济的发展提供了不可或缺的动力支持，而且在能源供应体系中起着“压舱石”和“稳定器”的作用，对确保国家能源安全和产业链供应链稳定至关重要。然而，随着浅部煤炭资源的日益枯竭，我国煤炭开采逐渐向深部发展，深井开采已成为必然趋势。深井开采面临着一系列严峻的挑战，如矿压显现加剧、巷道维护困难、煤岩破坏过程强化、冲击地压危险性增加、瓦斯压力增高、煤与瓦斯突出危险严重、深热矿井增加、气候条件恶化以及矿井生产费用升高、经济效益下降等问题。在这些挑战中，深井条带煤柱的长期稳定性问题尤为突出。条带煤柱作为深井开采中控制上覆岩层移动和维护开采空间稳定的重要结构，其长期稳定性直接关系到矿井的安全生产和资源的高效开采。如果条带煤柱在长期的地质作用和开采扰动下发生失稳破坏，可能导致顶板垮落、地表塌陷、瓦斯泄漏等严重的安全事故，不仅会造成巨大的经济损失，还可能危及人员生命安全。此外，条带煤柱的稳定性还会影响煤炭资源的回收率，不合理的煤柱设计或煤柱失稳可能导致大量煤炭资源被永久丢弃，降低资源的利用效率。因此，深入研究深井条带煤柱的长期稳定性具有重要的理论和实际意义。从理论方面来看，煤岩蠕变特性是影响条带煤柱长期稳定性的关键因素之一，但目前关于煤岩蠕变特性的研究仍存在诸多不足，对其内在机制的认识还不够深入。通过开展基于煤岩蠕变试验的深井条带煤柱长期稳定性研究，可以进一步揭示煤岩蠕变的力学机制和本构关系，丰富和完善煤岩力学理论体系。从实际应用角度出发，准确掌握深井条带煤柱的长期稳定性规律，能够为煤矿开采设计提供科学合理的依据，优化条带煤柱的尺寸和布局，提高煤炭资源的回收率，同时有效降低矿井开采过程中的安全风险，保障煤矿的安全生产，促进煤炭行业的可持续发展。1.2国内外研究现状1.2.1煤岩蠕变试验研究现状国内外众多学者对煤岩蠕变特性展开了广泛而深入的研究。在蠕变特性方面，大量试验表明煤岩具有明显的蠕变现象，其蠕变过程通常可分为初始蠕变、稳态蠕变和加速蠕变三个阶段。不同学者通过试验得到的煤岩蠕变特性存在一定差异，这主要是由于煤岩的成分、结构、赋存条件以及试验加载方式、加载应力水平等因素的不同所导致。例如，一些研究发现，随着加载应力的增加，煤岩的蠕变变形速率加快，稳态蠕变阶段缩短，加速蠕变阶段提前到来；而煤岩中的矿物成分、孔隙结构等也会对其蠕变特性产生显著影响，含黏土矿物较多的煤岩往往具有更强的蠕变特性。在试验方法上，常用的有单轴蠕变试验、三轴蠕变试验以及不同环境条件下的蠕变试验。单轴蠕变试验操作相对简单，能够初步获取煤岩在轴向应力作用下的蠕变特性；三轴蠕变试验则更能模拟煤岩在实际地下工程中的受力状态，考虑了围压对煤岩蠕变的影响。许多学者开展了不同围压条件下的煤岩三轴蠕变试验，结果表明围压的增大可以抑制煤岩的蠕变变形，提高煤岩的蠕变稳定性。此外，为了研究地下水、温度等环境因素对煤岩蠕变的影响，学者们还进行了相应的环境条件下的蠕变试验。研究发现，水对煤岩有显著的弱化作用，饱水煤岩的蠕变变形明显大于干燥煤岩，这是因为水的存在会使煤岩内部的化学键发生断裂，导致煤岩的强度降低，从而加速蠕变变形；温度的升高也会使煤岩的蠕变特性发生改变，一般来说，温度升高会使煤岩的蠕变变形增大，蠕变速率加快。在模型建立方面，国内外学者提出了多种煤岩蠕变本构模型。经典的蠕变模型如Maxwell模型、Kelvin模型、Burgers模型等，通过弹簧和黏壶的不同组合来描述煤岩的蠕变特性。Maxwell模型能够较好地描述煤岩的瞬时弹性变形和黏滞流动，但无法反映煤岩的稳态蠕变阶段；Kelvin模型则主要用于描述煤岩的稳态蠕变阶段，不能体现瞬时弹性变形；Burgers模型结合了Maxwell模型和Kelvin模型的优点，能够描述煤岩蠕变的三个阶段，但对于某些复杂的煤岩蠕变特性，其描述能力仍显不足。为了更准确地描述煤岩的蠕变特性，许多学者在经典模型的基础上进行了改进和拓展。例如，通过引入损伤变量来考虑煤岩在蠕变过程中的内部损伤演化，建立了损伤蠕变模型；考虑到煤岩的非线性特性，采用分数阶微积分理论建立了分数阶蠕变模型，这些改进后的模型在一定程度上提高了对煤岩蠕变特性的描述精度。尽管国内外在煤岩蠕变试验研究方面取得了丰硕的成果，但仍存在一些不足之处。一方面，目前的研究多集中在实验室条件下的短时间蠕变试验，对于煤岩在实际工程中长期的蠕变特性研究相对较少，而实际工程中煤岩的蠕变过程可能持续数年甚至数十年，短时间的试验结果难以准确反映煤岩的长期力学行为。另一方面，不同地区、不同赋存条件下的煤岩性质差异较大，现有的蠕变模型和研究成果在通用性和普适性方面还存在一定的局限，难以直接应用于各种复杂的工程实际。此外，对于煤岩蠕变过程中的微观机制研究还不够深入，虽然一些研究借助扫描电镜、压汞仪等微观测试手段对煤岩的微观结构变化进行了观察和分析，但对于微观结构变化与宏观蠕变特性之间的定量关系还缺乏深入的理解。1.2.2深井条带煤柱稳定性研究现状条带煤柱稳定性的影响因素是多方面的，主要包括煤柱的尺寸、强度、上覆岩层的载荷以及开采条件等。煤柱尺寸是影响其稳定性的重要因素之一，煤柱宽度过小可能导致煤柱强度不足而发生失稳破坏，而煤柱宽度过大则会降低煤炭资源的回收率；煤柱强度与煤岩的物理力学性质密切相关，煤岩的抗压强度、抗拉强度、内摩擦角等参数都会影响煤柱的承载能力；上覆岩层的载荷大小和分布形式直接决定了煤柱所承受的压力，开采深度的增加、采空区跨度的增大等都会使上覆岩层载荷增大，从而对煤柱稳定性产生不利影响；开采条件如开采顺序、开采速度等也会对条带煤柱的稳定性产生影响，不合理的开采顺序可能导致煤柱受力不均，加速煤柱的破坏。在分析方法上，常用的有条带煤柱载荷计算方法、煤柱强度计算方法以及稳定性评价方法。条带煤柱载荷计算方法主要有经验公式法、理论分析法和数值模拟法。经验公式法是根据大量的工程实践和现场观测数据总结得出的，具有简单实用的优点，但缺乏普遍的理论依据，适用范围有限；理论分析法基于弹性力学、塑性力学等理论，通过建立力学模型来计算煤柱载荷，如基于薄板理论、压力拱理论等的计算方法，但这些理论方法往往对实际工程进行了一定的简化假设，与实际情况存在一定的偏差；数值模拟法则借助有限元、离散元等数值计算软件，能够较为真实地模拟煤柱的受力和变形过程，考虑多种复杂因素的影响，但数值模拟结果的准确性依赖于模型的合理性和参数的选取。煤柱强度计算方法也有多种，如经验公式法、Hoek-Brown强度准则等，不同的计算方法适用于不同的煤岩条件和工程情况。在稳定性评价方面，常用的指标有煤柱安全系数、屈服接近度等，通过这些指标可以对条带煤柱的稳定性进行量化评价。国内外学者针对条带煤柱稳定性开展了大量的研究工作，并取得了一系列的成果。一些学者通过理论分析和数值模拟，研究了不同煤柱尺寸和开采条件下条带煤柱的稳定性规律，提出了合理的煤柱尺寸设计方法和开采方案；部分学者通过现场监测，对条带煤柱的实际受力和变形情况进行了长期观测，验证了理论分析和数值模拟的结果，并为工程实践提供了宝贵的经验。然而，当前对条带煤柱长期稳定性的研究还存在一定的欠缺。一方面，在以往的研究中，大多没有充分考虑煤岩的蠕变特性对条带煤柱长期稳定性的影响，而实际上煤岩的蠕变会导致煤柱的力学性质随时间发生变化，进而影响煤柱的长期稳定性。另一方面，对于条带煤柱在复杂地质条件和开采环境下的长期稳定性研究还不够深入，如深部高地应力、高水压、强采动等条件下条带煤柱的长期稳定性问题，还需要进一步开展系统的研究。1.3研究内容与方法1.3.1研究内容本文以煤岩蠕变试验为基础，深入研究深井条带煤柱的长期稳定性，主要涵盖以下几个方面的内容：煤岩蠕变特性研究：开展不同条件下的煤岩蠕变试验，包括单轴蠕变试验、三轴蠕变试验以及不同环境因素（如水、温度等）作用下的蠕变试验。通过试验获取煤岩的蠕变曲线，分析煤岩在不同应力水平和环境条件下的蠕变特性，包括蠕变变形规律、蠕变阶段划分以及蠕变参数的变化规律等。运用微观测试手段，如扫描电镜（SEM）、压汞仪（MIP）等，观察煤岩在蠕变过程中的微观结构变化，探讨微观结构变化与宏观蠕变特性之间的内在联系，揭示煤岩蠕变的微观机制。深井条带煤柱受力分析：综合考虑上覆岩层载荷、开采深度、采空区尺寸、煤柱尺寸等因素，运用理论分析方法建立深井条带煤柱的力学模型，计算条带煤柱在不同开采阶段所承受的载荷。分析煤柱在长期蠕变过程中的应力重分布规律，考虑煤岩蠕变导致的煤柱力学性质变化，研究煤柱内部应力随时间的演化过程。深井条带煤柱稳定性评价：基于煤岩蠕变特性和条带煤柱受力分析结果，建立考虑蠕变影响的深井条带煤柱稳定性评价指标和方法。采用数值模拟软件，如FLAC3D、ANSYS等，对深井条带煤柱的长期稳定性进行模拟分析，研究不同因素（如煤柱尺寸、开采条件、煤岩蠕变参数等）对条带煤柱稳定性的影响规律。结合工程实际，对某深井条带煤柱的稳定性进行实例分析，验证所建立的稳定性评价方法的有效性和可靠性。工程应用与建议：根据研究成果，为深井条带煤柱的设计和开采提供工程应用建议。优化条带煤柱的尺寸和布局，提出合理的开采顺序和开采速度，以提高条带煤柱的长期稳定性，保障矿井的安全生产。针对煤岩蠕变特性和条带煤柱稳定性，提出相应的监测方案和控制措施，以便及时发现和处理潜在的安全隐患。1.3.2研究方法本文综合运用实验室试验、理论分析、数值模拟和现场监测等方法开展研究，具体如下：实验室试验：通过在实验室进行煤岩蠕变试验，获取煤岩在不同条件下的蠕变数据。采用高精度的试验设备，如微机控制电液伺服岩石三轴试验机、高温高压蠕变试验装置等，确保试验数据的准确性和可靠性。设计合理的试验方案，控制试验变量，全面研究煤岩的蠕变特性。同时，利用微观测试设备对煤岩的微观结构进行分析，为揭示煤岩蠕变机制提供微观依据。理论分析：运用弹性力学、塑性力学、岩石力学等相关理论，对深井条带煤柱的受力状态进行分析。建立条带煤柱的力学模型，推导煤柱载荷和应力的计算公式。结合煤岩蠕变理论，分析煤柱在蠕变过程中的力学行为变化，为条带煤柱稳定性评价提供理论基础。数值模拟：借助数值模拟软件对深井条带煤柱的开采过程进行模拟。建立三维数值模型，考虑煤岩的非线性力学特性、蠕变特性以及开采过程中的各种复杂因素。通过数值模拟，直观地展示条带煤柱在开采过程中的受力、变形和破坏过程，分析不同因素对条带煤柱稳定性的影响，预测条带煤柱的长期稳定性。现场监测：选择典型的深井条带煤柱开采工程现场，布置监测系统，对条带煤柱的实际受力、变形和位移等参数进行长期监测。通过现场监测数据，验证实验室试验、理论分析和数值模拟的结果，及时发现工程实际中存在的问题，并为研究提供实际工程数据支持。同时，根据现场监测结果，对条带煤柱的稳定性进行实时评估，为矿井安全生产提供决策依据。二、煤岩蠕变试验原理与方法2.1蠕变的基本概念与原理蠕变是指材料在恒定应力作用下，应变随时间逐渐增大的现象。这一现象广泛存在于各种材料中，如金属、塑料、岩土、冰等，在土木工程领域也被称为徐变。对于煤岩而言，蠕变特性对其在地下工程中的力学行为有着至关重要的影响。煤岩的蠕变过程通常可分为三个阶段，每个阶段都具有独特的变形特征和机制。在初始蠕变阶段，也称为过渡蠕变阶段，当煤岩受到恒定应力作用时，应变会迅速增加，但增加的速度逐渐减慢，即蠕变应变率下降。这是因为在加载初期，煤岩内部的结构处于调整适应阶段，内部的微裂隙和孔隙开始逐渐闭合或扩展，导致应变快速增加。随着时间的推移，煤岩内部结构逐渐调整稳定，抵抗变形的能力增强，使得应变增加的速度减缓。在这一阶段，煤岩的变形主要是由弹性变形和部分塑性变形组成，弹性变形部分在卸载后可以恢复，而塑性变形部分则是不可恢复的。随着时间的进一步延续，煤岩进入稳态蠕变阶段，也称为定常蠕变阶段。在这个阶段，应变随时间延续而匀速增加，蠕变应变率几乎保持不变。此时，煤岩内部的变形机制达到了一种动态平衡，内部结构的调整和损伤的发展速率基本稳定。虽然煤岩内部仍然存在着微裂隙的扩展和新裂隙的产生，但同时也存在着裂隙的闭合和愈合，使得总体的变形速率保持相对稳定。稳态蠕变阶段持续的时间相对较长，它是衡量煤岩长期稳定性的重要阶段。当煤岩经历了长时间的稳态蠕变后，若应力仍然保持不变，煤岩将进入加速蠕变阶段。在这一阶段，应变随时间延续而加速增加，蠕变应变率急剧上升，直至煤岩发生破裂。加速蠕变阶段的出现标志着煤岩内部结构的严重破坏。随着时间的推移，煤岩内部的微裂隙不断扩展、贯通，形成宏观的裂隙网络，导致煤岩的承载能力迅速下降，变形急剧增大，最终导致煤岩的失稳破坏。煤岩蠕变的产生机制较为复杂，与煤岩的微观结构密切相关。煤岩是一种由多种矿物成分和有机物质组成的复杂材料，其内部存在着大量的微裂隙、孔隙和弱面。在应力作用下，这些微结构会发生变形和破坏。微裂隙的扩展和贯通是煤岩蠕变的重要机制之一。当煤岩受到应力作用时，微裂隙尖端会产生应力集中，随着时间的推移，应力集中导致微裂隙不断扩展。当微裂隙相互贯通时，会形成更大的裂隙，从而导致煤岩的变形增大。此外，煤岩内部矿物颗粒之间的相对滑动和位错运动也是蠕变产生的原因。在应力作用下，矿物颗粒之间的摩擦力和粘结力不足以抵抗外力，使得颗粒之间发生相对滑动和位错。这种微观的运动在宏观上表现为煤岩的蠕变变形。煤岩中的水分和气体也会对蠕变产生影响。水分的存在会使煤岩发生软化，降低其强度，从而加速蠕变变形；气体的吸附和解吸作用会改变煤岩的内部结构和力学性质，进而影响蠕变特性。2.2煤岩蠕变试验设备与方法2.2.1试验设备在煤岩蠕变试验中，MTS伺服控制岩石三轴试验机是一种常用且重要的设备，它在模拟煤岩复杂受力状态以及精确测量煤岩力学响应方面发挥着关键作用。以某型号的MTS伺服控制岩石三轴试验机为例，其工作原理基于电液伺服控制技术。该试验机通过液压系统提供稳定且可精确控制的压力，并借助伺服阀来实现对压力和位移的高精度闭环控制。在试验过程中，可根据试验需求设定不同的加载模式和加载速率，能够模拟煤岩在地下工程中所承受的各种复杂应力状态，包括单轴、三轴等加载条件。该试验机具备轴压、围压和孔隙水压三套独立的闭环伺服控制功能，这使得它能够全面考虑煤岩在实际工程中的受力情况。在进行三轴蠕变试验时，轴压系统可以精确施加轴向压力，模拟上覆岩层对煤岩的垂直作用力；围压系统则能够提供侧向压力，反映煤岩周围岩体的约束作用；孔隙水压系统可以模拟地下水对煤岩的影响，通过控制孔隙水压力的大小，研究水-岩相互作用对煤岩蠕变特性的影响。其主要技术参数也十分出色，试验框架整体刚度可达11.0×109N/m，这保证了试验机在加载过程中的稳定性和准确性，能够有效减少因设备自身变形而带来的测量误差。最大轴向力可达4600kN，能够满足对大多数煤岩试件进行加载的需求，可以模拟深部矿井中煤岩所承受的高应力环境。最大轴向拉力出力为2300kN，可用于研究煤岩在拉伸状态下的力学性能。最大围压可达140MPa，能够模拟深部复杂的地应力条件，深入研究围压对煤岩蠕变特性的影响。最大孔隙水压力同样为140MPa，可实现对饱水煤岩在不同孔隙水压力下蠕变特性的研究。除了MTS伺服控制岩石三轴试验机外，还有其他一些设备也常用于煤岩蠕变试验。例如，微机控制电液伺服岩石三轴试验机，它同样能够模拟岩石在自地表至地下约1000千米埋深围压、常温至200摄氏度温度范围的岩石力学试验，并可进行单向低周循环及用户自行设置的组合波形程序控制等多种控制方式，能模拟岩石在高温高压环境下应力-温度-地下水耦合作用下的三轴试验。MT-C60-D800岩石力学多场耦合试验系统，能开展排水与不排水条件下的静水压力试验、三轴压缩试验、静水压力流变试验、三轴压缩流变试验、三轴压缩松弛试验、比例加卸载试验、侧向与轴向加卸载循环试验、静水压力和三轴应力状态下的液体渗透试验以及静水压力和三轴应力状态下的流体耦合试验等多种试验，为研究煤岩在复杂多场耦合作用下的蠕变特性提供了有力的手段。2.2.2试验方法煤岩蠕变试验的流程较为严谨，包括试件制备、安装、加载和数据采集等多个关键环节。在试件制备阶段，首先需要从煤矿现场采集具有代表性的煤岩样本。采样时应严格遵循相关标准，确保样本能够真实反映现场煤岩的特性。采集回来的煤岩样本需要加工成标准试件，对于常规的煤岩蠕变试验，标准试件的尺寸通常为直径50±2mm，高度为直径的2.0-2.5倍，即高度在100-125mm之间。试件的加工精度要求较高，试件两端面应平行且光滑，其不平行度不得大于0.05mm，试件上、下端直径偏差不得大于0.3mm，以保证在试验过程中试件受力均匀。试件制备完成后，进行安装操作。将加工好的试件小心地安装在试验设备的三轴压力室中。在安装过程中，要确保试件与加载装置紧密接触，避免出现偏心加载的情况。同时，安装好各种测量传感器，如位移传感器用于测量试件的轴向和径向变形，压力传感器用于测量轴压、围压和孔隙水压力等参数。在安装位移传感器时，应将其准确地固定在试件的特定位置，以确保能够精确测量试件的变形。安装压力传感器时，要保证其与压力通道连接紧密，无泄漏现象，以准确测量各方向的压力。加载环节是试验的核心部分。加载方式通常分为分级加载和恒定加载两种。分级加载是按照一定的应力增量逐步增加荷载，每级荷载保持一定的时间，待应变稳定后再施加下一级荷载。例如，首先施加初始应力为煤岩单轴抗压强度的30%，保持24小时，记录这期间的应变随时间变化数据；然后增加应力至单轴抗压强度的40%，同样保持24小时并记录数据，以此类推。这种加载方式可以更全面地研究煤岩在不同应力水平下的蠕变特性，能够获取煤岩在不同应力阶段的蠕变参数。恒定加载则是一次性施加预定的恒定荷载，然后持续监测应变随时间的变化。在进行恒定加载时，需要根据试验目的和煤岩的特性确定合适的恒定荷载值。在整个试验过程中，数据采集至关重要。需要实时采集和记录多个测量参数，包括轴向位移、径向位移、轴压、围压、孔隙水压力以及时间等。通过这些数据可以绘制出煤岩的蠕变曲线，即应变-时间曲线。利用专业的数据采集系统，可以按照设定的时间间隔自动采集数据，确保数据的准确性和完整性。一般数据采集的时间间隔可以根据试验的具体情况设定为1分钟、5分钟或10分钟等。通过对采集到的数据进行分析，可以深入研究煤岩的蠕变特性，如确定蠕变阶段、计算蠕变参数等。2.3煤岩蠕变试验案例分析2.3.1某煤矿煤岩单轴蠕变试验以某煤矿的煤样为研究对象，开展了单轴蠕变试验。在试验前，首先从该煤矿的特定煤层采集了煤样，并将其加工成直径为50mm，高度为100mm的标准试件。为了保证试验结果的可靠性，共制备了3个试件。试验过程中，采用分级加载的方式。首先对试件施加5MPa的初始应力，保持该应力水平，并每隔10分钟记录一次试件的轴向变形。在初始应力作用下，试件的轴向应变随时间迅速增加，在开始的1-2小时内，应变增加较为明显，随后应变增加的速率逐渐减慢，进入初始蠕变阶段的后期。经过约10小时后，应变基本趋于稳定，此时认为试件在该应力水平下的初始蠕变阶段结束。接着，将应力增加至10MPa，同样保持该应力并持续监测轴向变形。在新的应力水平下，轴向应变又开始迅速上升，重复了类似初始应力作用下的变形过程，但应变增加的幅度和速率均大于初始应力阶段。随着时间的推移，应变增加速率逐渐降低，再次进入稳定状态，完成了该应力水平下的初始蠕变阶段。按照同样的方法，依次将应力增加至15MPa、20MPa，并记录相应的应变-时间数据。从试验结果得到的不同应力水平下的蠕变曲线可以看出，随着应力的增加，煤样的初始蠕变阶段应变增加的幅度增大，稳态蠕变阶段的应变率也明显增大。在5MPa应力下，稳态蠕变阶段的应变率约为0.005/h；而在20MPa应力下，稳态蠕变阶段的应变率增大至0.02/h左右。这表明应力对煤岩的蠕变特性有着显著的影响，高应力会加速煤岩的蠕变变形，缩短稳态蠕变阶段的持续时间，使煤岩更容易进入加速蠕变阶段。2.3.2某煤矿煤岩三轴蠕变试验某煤矿进行了煤岩三轴蠕变试验，旨在研究不同围压和轴压条件下煤岩的蠕变特性。试验选取了该煤矿的典型煤样，并将其加工成标准试件。在试验过程中，设置了不同的围压和轴压组合。首先进行了围压为5MPa，轴压分别为10MPa、15MPa、20MPa的蠕变试验。当围压为5MPa，轴压为10MPa时，在加载初期，试件的轴向应变和径向应变都迅速增加，随着时间的推移，应变增加速率逐渐减小，进入稳态蠕变阶段。在稳态蠕变阶段，轴向应变率相对较低，约为0.003/h，径向应变率也较小。当轴压增加到15MPa时，初始蠕变阶段的应变增加幅度明显增大，稳态蠕变阶段的轴向应变率增大至0.006/h左右，径向应变率也有所增加。轴压进一步增加到20MPa时，初始蠕变阶段的应变急剧增加，稳态蠕变阶段的轴向应变率达到0.01/h以上，且加速蠕变阶段提前到来。对比不同围压条件下的蠕变特性发现，围压对煤岩的蠕变具有显著的抑制作用。当围压从5MPa增加到10MPa，轴压保持15MPa时，在相同的时间内，试件的轴向应变和径向应变都明显减小。在围压为5MPa时，经过24小时的蠕变，轴向应变达到0.3%；而围压增加到10MPa时，相同时间内轴向应变仅为0.2%。这是因为围压的增大限制了煤岩内部微裂隙的扩展，增强了煤岩的抵抗变形能力，从而抑制了蠕变变形。通过该煤矿的煤岩三轴蠕变试验，可以清晰地认识到围压和轴压对煤岩蠕变特性的影响规律，为深井条带煤柱稳定性分析提供了重要的试验依据，在实际工程中，合理地控制围压和轴压条件，对于保障煤柱的长期稳定性具有重要意义。三、深井条带煤柱受力分析3.1深井开采的特点与挑战深井开采是指在深部地层中进行的煤炭开采活动，一般将开采深度超过1000m的矿井视为深井。随着煤炭资源的不断开采，浅部煤炭储量逐渐减少，深井开采已成为煤炭行业发展的必然趋势。然而，深井开采与浅部开采相比，面临着诸多特殊的问题和挑战，这些问题对条带煤柱的稳定性产生了深远的影响。地应力是指地壳岩体在各种地质作用下所产生的内应力。随着开采深度的增加，上覆岩层的重量增大，地应力也随之显著升高。研究表明，地应力与开采深度呈近似线性关系，一般每增加100m深度，地应力增加2.5-3.5MPa。在深部矿井中，地应力往往高达数十MPa甚至更高。高地应力会使煤岩的力学性质发生显著变化，导致煤岩的强度降低、脆性增加。在高地应力作用下，条带煤柱更容易发生破坏，其稳定性受到严重威胁。煤柱在高地应力下可能会出现脆性破裂、片帮等现象，导致煤柱的承载能力下降，进而引发顶板垮落等安全事故。深部岩体由于长期受到高地应力的作用，其内部结构往往较为破碎，存在大量的节理、裂隙等弱面。这些破碎的围岩使得条带煤柱的边界条件变得复杂，煤柱与围岩之间的相互作用更加剧烈。破碎的围岩容易向煤柱施加不均匀的压力，导致煤柱受力不均，从而增加了煤柱失稳的风险。破碎围岩的变形和移动也会对煤柱产生附加的作用力，进一步影响煤柱的稳定性。随着开采深度的增加，地温梯度逐渐增大，一般每增加100m深度，地温升高2-3℃。在深部矿井中，地温可能会超过35℃，甚至更高。高温环境会使煤岩的物理力学性质发生改变，如热膨胀、热软化等。煤岩的热膨胀会导致其内部产生应力集中，热软化则会降低煤岩的强度。这些变化都会对条带煤柱的稳定性产生不利影响。高温还会加速煤岩的蠕变过程，使煤柱的变形随时间不断增加，进一步降低煤柱的长期稳定性。深部矿井中，瓦斯含量和瓦斯压力通常较高。瓦斯在煤岩中以吸附态和游离态存在，当煤岩受到开采扰动时，瓦斯会解吸并释放出来。高瓦斯压力会对条带煤柱产生附加的膨胀力，使煤柱内部的应力状态更加复杂。瓦斯的存在还会降低煤岩的强度，增加煤柱发生破坏的可能性。如果瓦斯突然大量涌出，可能会引发煤与瓦斯突出等严重的安全事故，对条带煤柱的稳定性造成极大的破坏。综上所述，深井开采的地应力高、围岩破碎、温度高、瓦斯压力大等特点，给条带煤柱的稳定性带来了严峻的挑战。在深井条带煤柱的设计和开采过程中，必须充分考虑这些因素的影响，采取有效的措施来提高条带煤柱的稳定性，确保矿井的安全生产。3.2条带煤柱的载荷计算方法条带煤柱的载荷计算是研究其稳定性的基础，准确计算条带煤柱所承受的载荷对于合理设计煤柱尺寸、保障矿井安全具有重要意义。目前，常用的条带煤柱载荷计算方法主要有经验公式法、理论计算法和数值模拟法。经验公式法是根据大量的工程实践和现场观测数据总结得出的，具有简单实用的优点。例如，常用的威尔逊（Wilson）公式，该公式基于大量的现场实测数据，考虑了煤柱宽度、采空区宽度、上覆岩层平均重度等因素对煤柱载荷的影响。其表达式为：P=\frac{\gammaH(a+b)}{a}其中，P为煤柱单位面积上的载荷（MPa）；\gamma为上覆岩层平均重度（kN/m³）；H为开采深度（m）；a为煤柱宽度（m）；b为采空区宽度（m）。再如，我国学者根据国内煤矿的实际情况，提出了一些适合我国国情的经验公式。其中一种考虑了煤柱的强度和稳定性的经验公式为：P=k\gammaH式中，k为载荷系数，其取值与煤柱的宽高比、煤岩性质、开采条件等因素有关，一般通过现场实测或经验取值。经验公式法的优点是计算简单、快捷，在工程实践中应用较为广泛。然而，该方法也存在明显的局限性。由于经验公式是基于特定的工程条件和数据总结得出的，缺乏普遍的理论依据，其适用范围有限。不同地区、不同地质条件下的煤矿，煤岩性质、开采条件等差异较大，直接套用经验公式可能会导致计算结果与实际情况偏差较大。经验公式往往难以考虑到一些复杂因素对煤柱载荷的影响，如采动影响、地质构造等。理论计算法是基于弹性力学、塑性力学等理论，通过建立力学模型来计算条带煤柱载荷的方法。基于薄板理论的计算方法，将上覆岩层视为弹性薄板，煤柱视为薄板的支撑。根据薄板的弯曲理论，推导得出煤柱所承受的载荷计算公式。假设上覆岩层为均布载荷作用下的简支薄板，煤柱宽度为a，采空区宽度为b，薄板厚度为h，弹性模量为E，泊松比为\nu，则煤柱所承受的载荷P可通过以下公式计算：P=\frac{q(b+a)^4}{8a^3h^3}\left(1-\nu^2\right)其中，q为上覆岩层单位面积上的均布载荷（MPa）。基于压力拱理论的计算方法认为，在采空区上方会形成压力拱，条带煤柱承受的载荷主要来自压力拱范围内的上覆岩层重量。通过建立压力拱的力学模型，可计算出煤柱所承受的载荷。设压力拱的跨度为L，高度为h_0，上覆岩层平均重度为\gamma，则煤柱所承受的载荷P可表示为：P=\frac{\gammah_0L}{a}理论计算法具有一定的理论依据，能够考虑到煤柱和上覆岩层的力学性质等因素。但该方法在实际应用中也存在一些问题。理论计算法往往对实际工程进行了一定的简化假设，如假设煤岩为均质、连续、各向同性材料，忽略了地质构造、节理裂隙等因素的影响，这些简化假设与实际情况存在一定的偏差，可能导致计算结果的准确性受到影响。理论计算过程通常较为复杂，需要具备较高的力学知识和数学计算能力，在实际工程应用中存在一定的难度。数值模拟法是借助有限元、离散元等数值计算软件，对条带煤柱的开采过程进行模拟，从而计算出煤柱所承受的载荷。以有限元软件ANSYS为例，首先需要建立条带煤柱和上覆岩层的三维有限元模型，定义煤岩的材料参数、边界条件和载荷条件。在模型中，将煤岩视为连续介质，采用合适的单元类型进行网格划分。然后，通过施加开采过程中的位移边界条件，模拟采空区的形成过程。最后，利用软件的求解器进行计算，得到条带煤柱在不同开采阶段的应力分布和载荷大小。离散元软件UDEC则是将煤岩视为由离散的颗粒组成，通过模拟颗粒之间的相互作用来计算煤柱的力学响应。在UDEC模型中，定义颗粒的力学参数、接触模型和边界条件，模拟开采过程中颗粒的运动和相互作用，从而得到条带煤柱的载荷和变形情况。数值模拟法能够较为真实地模拟条带煤柱的受力和变形过程，考虑多种复杂因素的影响，如煤岩的非线性力学特性、采动影响、地质构造等。通过数值模拟，可以直观地观察到煤柱内部的应力分布和变化规律，为煤柱稳定性分析提供详细的信息。然而，数值模拟结果的准确性依赖于模型的合理性和参数的选取。如果模型建立不合理或参数选取不准确，可能会导致模拟结果与实际情况相差较大。数值模拟计算量较大，需要较高的计算机硬件配置和较长的计算时间。综上所述，经验公式法、理论计算法和数值模拟法各有优缺点。在实际工程应用中，应根据具体情况选择合适的计算方法，或者将多种方法相结合，以提高条带煤柱载荷计算的准确性和可靠性。3.3条带煤柱的应力分布规律3.3.1基于理论分析的应力分布条带煤柱的应力分布规律是研究其稳定性的关键内容，运用弹性力学和塑性力学理论可以对条带煤柱的应力分布进行深入分析。在弹性力学中，通常假设煤岩为连续、均质、各向同性的弹性体。对于受均布载荷作用的条带煤柱，可将其简化为平面应变问题进行研究。根据弹性力学的基本原理，在平面应变条件下，煤柱内的应力分量\sigma_x、\sigma_y和\tau_{xy}满足平衡方程和几何方程。假设煤柱宽度为a，高度为h，上覆岩层施加的均布载荷为q，采用直角坐标系，x轴沿煤柱宽度方向，y轴沿煤柱高度方向。根据圣维南原理，在远离煤柱边界的区域，煤柱内的应力分布可近似看作均匀分布。在弹性阶段，煤柱内的垂直应力\sigma_y和水平应力\sigma_x可通过以下公式计算：\sigma_y=q\sigma_x=\frac{\nu}{1-\nu}q其中，\nu为煤岩的泊松比。然而，实际的煤岩并非完全符合弹性力学的假设，在高应力作用下，煤岩会发生塑性变形。此时，需要运用塑性力学理论来分析煤柱的应力分布。基于Mohr-Coulomb屈服准则，当煤岩中的剪应力达到一定程度时，煤岩会发生屈服破坏。在条带煤柱中，随着载荷的增加，煤柱边缘首先进入塑性状态，形成塑性区。假设煤柱的屈服准则为\tau=c+\sigma\tan\varphi，其中c为煤岩的黏聚力，\varphi为内摩擦角。在塑性区内，应力分布满足塑性力学的相关方程。通过求解塑性力学的边值问题，可以得到塑性区内的应力分布表达式。在塑性区内，垂直应力\sigma_y和水平应力\sigma_x的分布不再均匀，且与煤柱的边界条件和屈服准则密切相关。随着塑性区的不断扩展，煤柱的承载能力逐渐下降，当塑性区贯通整个煤柱时，煤柱将发生失稳破坏。从理论分析可知，条带煤柱的应力分布具有明显的特征。在弹性阶段，煤柱内的应力分布相对均匀，垂直应力主要由上覆岩层载荷决定，水平应力与泊松比相关。随着载荷的增加，煤柱边缘进入塑性状态，塑性区逐渐向煤柱内部扩展，应力分布变得不均匀。煤柱的应力分布还受到煤岩的力学性质（如黏聚力、内摩擦角、弹性模量等）、煤柱尺寸（宽度、高度）以及开采条件（上覆岩层载荷、采空区尺寸等）的影响。不同的煤岩力学性质和煤柱尺寸会导致应力分布的差异，从而影响煤柱的稳定性。3.3.2基于数值模拟的应力分布利用FLAC3D等数值模拟软件建立数值模型，可以更加直观、准确地模拟条带煤柱的应力分布，并对比不同开采条件下的应力变化。以某深井条带煤柱开采工程为例，运用FLAC3D软件进行数值模拟。首先，根据实际工程地质条件，建立三维数值模型。模型的尺寸为：长×宽×高=500m×300m×200m，其中条带煤柱的宽度为30m，采空区宽度为50m，开采深度为1200m。模型中煤岩的力学参数根据现场试验和相关资料确定，采用摩尔-库仑本构模型来描述煤岩的力学行为。在模型中，对模型的边界条件进行设置。模型的底面固定，限制其在x、y、z三个方向的位移；模型的侧面施加水平约束，限制其在x和y方向的位移；模型的顶部施加与开采深度对应的垂直应力，以模拟上覆岩层的载荷。在模拟开采过程时，采用分步开挖的方式，逐步形成采空区，观察条带煤柱在不同开采阶段的应力变化。通过数值模拟得到的条带煤柱应力分布云图可以清晰地看到，在开采初期，条带煤柱的应力分布相对均匀，主要承受上覆岩层的垂直压力。随着采空区的逐渐扩大，煤柱两侧的应力集中现象逐渐明显。在煤柱与采空区的交界处，垂直应力和水平应力都显著增大，形成应力集中区域。这是因为采空区的形成导致上覆岩层的压力重新分布，煤柱承担了更多的载荷。对比不同开采条件下的应力变化发现，煤柱宽度对煤柱应力分布有显著影响。当煤柱宽度增加时，煤柱内部的应力集中程度有所降低，应力分布更加均匀。这是因为较宽的煤柱具有更强的承载能力，能够更好地分散上覆岩层的压力。开采深度的增加也会使煤柱的应力显著增大。随着开采深度的增加，上覆岩层的压力增大，煤柱所承受的载荷也相应增加，导致煤柱内部的应力水平升高。采空区宽度的变化也会对煤柱应力分布产生影响。采空区宽度越大，煤柱两侧的应力集中越明显，煤柱的稳定性越差。通过数值模拟还可以分析不同开采顺序对条带煤柱应力分布的影响。采用不同的开采顺序，如顺序开采、跳采等，模拟结果表明，合理的开采顺序可以有效降低煤柱的应力集中程度，提高煤柱的稳定性。顺序开采时，煤柱依次承受采动影响，应力集中较为明显；而跳采时，煤柱的受力相对均匀，应力集中程度降低。综上所述，通过FLAC3D等数值模拟软件建立的数值模型，能够直观地展示条带煤柱在不同开采条件下的应力分布和变化规律。数值模拟结果与理论分析结果相互印证，为深入研究条带煤柱的受力特性和稳定性提供了有力的支持，在实际工程中，可根据数值模拟结果优化条带煤柱的设计和开采方案，提高煤矿开采的安全性和经济性。3.4条带煤柱受力的影响因素条带煤柱的受力受到多种因素的综合影响，深入分析这些因素对于准确掌握煤柱的力学行为和保障矿井安全生产具有重要意义。煤柱宽度是影响其受力的关键因素之一。当煤柱宽度较小时，煤柱的承载能力有限，难以承受上覆岩层的压力。在这种情况下，煤柱内部的应力集中现象较为明显，容易发生破坏。随着煤柱宽度的增加，煤柱的承载能力增强，能够更好地分散上覆岩层的压力，煤柱内部的应力分布相对均匀，应力集中程度降低。然而，煤柱宽度过大也会导致煤炭资源的浪费，降低煤炭回收率。因此，在实际工程中，需要综合考虑煤柱的稳定性和煤炭资源回收率，合理确定煤柱宽度。煤柱高度的变化也会对其受力产生影响。煤柱高度增加，煤柱自身的重量增大，同时上覆岩层对煤柱的压力作用面积也增大。这会导致煤柱内部的应力水平升高，尤其是垂直应力显著增加。过高的煤柱在受到较大的压力时，更容易发生压缩变形和破坏。在设计条带煤柱时，应根据开采深度、煤岩力学性质等因素，合理控制煤柱高度，以保证煤柱的稳定性。开采深度的增加意味着上覆岩层的重量增大，地应力升高。条带煤柱所承受的载荷随开采深度的增加而显著增大。深部开采时，高地应力会使煤柱的力学性质发生改变，煤柱更容易发生脆性破裂和塑性变形。开采深度的增加还会导致煤柱周围的岩体更加破碎，增加了煤柱的支护难度和受力复杂性。在深井开采中，必须充分考虑开采深度对条带煤柱受力的影响，采取有效的支护和加固措施，以提高煤柱的稳定性。不同的开采方式对条带煤柱的受力有着不同的影响。例如，一次采全高开采方式会使煤柱在短时间内承受较大的采动影响，煤柱的应力变化较为剧烈。而分层开采方式则可以减小每次开采对煤柱的影响，使煤柱的应力变化相对平缓。开采顺序也会影响四、基于煤岩蠕变试验的条带煤柱长期稳定性评价4.1煤岩蠕变特性对条带煤柱稳定性的影响煤岩蠕变特性是影响条带煤柱稳定性的关键因素之一，其对条带煤柱稳定性的影响主要体现在强度降低和变形增大两个方面。在长期的蠕变过程中，煤岩内部的微观结构会发生显著变化。煤岩是一种含有大量微裂隙和孔隙的非均质材料，在应力作用下，微裂隙尖端会产生应力集中，随着时间的推移，应力集中导致微裂隙不断扩展。这些微裂隙的扩展和贯通会逐渐破坏煤岩的内部结构，削弱煤岩颗粒之间的连接力，从而使煤岩的强度不断降低。相关研究表明，经过长时间的蠕变后，煤岩的抗压强度可能会降低20%-50%。煤岩强度的降低对条带煤柱的稳定性产生了严重的负面影响。条带煤柱在开采过程中承担着上覆岩层的载荷，其稳定性依赖于自身的强度。当煤岩强度因蠕变而降低时，煤柱的承载能力相应下降。在相同的上覆岩层载荷作用下，强度降低的煤柱更容易进入塑性状态，导致塑性区范围扩大。一旦煤柱的塑性区贯通，煤柱就会失去承载能力，发生失稳破坏。在深部开采中，由于地应力较高，煤岩蠕变导致的强度降低更为明显，条带煤柱的稳定性问题也更加突出。煤岩的蠕变还会导致其变形不断增大。在蠕变的初始阶段，煤岩的变形速率相对较快，随着时间的推移，变形速率逐渐减缓并进入稳态蠕变阶段。然而，即使在稳态蠕变阶段，煤岩的变形仍在持续增加。如果应力超过一定阈值，煤岩将进入加速蠕变阶段，变形速率急剧增大。煤岩变形的增大对条带煤柱的稳定性也带来了诸多不利影响。一方面，煤柱的变形增大会导致其与周围岩体之间的相互作用发生变化。煤柱的过度变形可能会使其与顶板或底板之间出现间隙，从而削弱了煤柱对顶板的支撑作用，增加了顶板垮落的风险。另一方面，煤柱的变形增大还会导致其内部应力重新分布。随着变形的增加，煤柱内部的应力集中现象会更加严重，进一步加速煤柱的破坏。在一些开采实例中，由于煤岩蠕变导致煤柱变形过大，煤柱出现了明显的鼓帮和底鼓现象，最终导致煤柱失稳，引发了顶板大面积垮落事故。综上所述，煤岩蠕变特性导致的强度降低和变形增大对条带煤柱的稳定性有着至关重要的影响。在深井条带煤柱的设计和开采过程中，必须充分考虑煤岩蠕变特性，采取有效的措施来提高条带煤柱的长期稳定性。4.2条带煤柱长期稳定性评价指标与方法4.2.1评价指标位移：条带煤柱的位移是反映其稳定性的重要指标之一。煤柱的位移包括垂直位移和水平位移。垂直位移主要是由于上覆岩层的压力作用导致煤柱的压缩变形，而水平位移则可能是由于煤柱与周围岩体之间的相互作用以及开采扰动等因素引起的。过大的位移可能表明煤柱已经进入塑性变形阶段，其稳定性受到威胁。通过监测煤柱的位移变化，可以及时发现煤柱的变形趋势，预测煤柱的稳定性。在某煤矿的条带煤柱开采过程中，通过在煤柱表面布置位移传感器，实时监测煤柱的位移。当发现煤柱的垂直位移在短时间内急剧增加时，及时采取了加固措施，避免了煤柱的失稳破坏。应力：应力是衡量条带煤柱受力状态的关键指标。煤柱在开采过程中承受着上覆岩层的载荷以及周围岩体的约束作用，其内部应力分布复杂。通过分析煤柱的应力大小和分布情况，可以判断煤柱是否处于稳定状态。如果煤柱内部的应力超过了其强度极限，煤柱就会发生破坏。常用的应力指标包括垂直应力、水平应力和剪应力等。在数值模拟分析中，通过计算煤柱内部的应力分布云图，可以直观地了解煤柱的应力状态，为煤柱稳定性评价提供依据。应变：应变是描述煤柱变形程度的物理量。煤柱在受力过程中会发生弹性应变和塑性应变。弹性应变在卸载后可以恢复，而塑性应变则是不可恢复的。当煤柱的塑性应变达到一定程度时，煤柱的结构会发生破坏，其稳定性降低。通过测量煤柱的应变，可以了解煤柱的变形特征，判断煤柱的稳定性。在实验室试验中，可以通过应变片等设备测量煤岩试件在不同应力条件下的应变，进而分析煤岩的蠕变特性对煤柱应变的影响。塑性区范围：塑性区是指煤柱中发生塑性变形的区域。随着开采的进行和煤岩蠕变的发展，煤柱内部的应力不断变化，当应力超过煤岩的屈服强度时，煤柱就会进入塑性状态，形成塑性区。塑性区的范围大小直接反映了煤柱的破坏程度。如果塑性区范围不断扩大并贯通整个煤柱，煤柱将失去承载能力，发生失稳破坏。在数值模拟中，可以通过设置屈服准则来判断煤柱中塑性区的范围。通过分析塑性区的范围及其变化趋势，可以评估条带煤柱的稳定性。在某深部矿井的条带煤柱稳定性研究中，利用数值模拟软件分析了不同开采阶段煤柱的塑性区范围。结果表明，随着开采深度的增加和开采时间的延长，煤柱的塑性区范围逐渐扩大，当塑性区范围超过煤柱宽度的一定比例时，煤柱的稳定性急剧下降。承载能力：承载能力是条带煤柱稳定性的核心指标。它反映了煤柱能够承受的最大载荷。煤柱的承载能力与煤岩的强度、煤柱的尺寸和形状以及开采条件等因素密切相关。在煤岩蠕变的影响下，煤柱的承载能力会随着时间的推移而逐渐降低。通过计算煤柱的承载能力，并与实际承受的载荷进行比较，可以判断煤柱是否处于稳定状态。常用的计算煤柱承载能力的方法有经验公式法、理论计算法和数值模拟法等。在实际工程中，应根据具体情况选择合适的方法来评估煤柱的承载能力。4.2.2评价方法极限平衡法：极限平衡法是一种经典的条带煤柱稳定性评价方法，它基于极限平衡理论，通过分析煤柱在极限状态下的受力情况来判断其稳定性。该方法假设煤柱处于极限平衡状态时，其内部的应力分布满足一定的屈服准则，如Mohr-Coulomb屈服准则。根据该准则，可以建立煤柱的极限平衡方程，通过求解方程得到煤柱的极限承载能力或安全系数。以某条带煤柱为例，假设煤柱的宽度为b，高度为h，上覆岩层施加的均布载荷为q，煤岩的黏聚力为c，内摩擦角为\varphi。根据Mohr-Coulomb屈服准则，煤柱的极限平衡方程可以表示为：q=\frac{2c\cot\varphi+\gammah}{1-\sin\varphi}\left(1-\frac{b_0}{b}\right)其中，b_0为煤柱的极限平衡区宽度，\gamma为上覆岩层的重度。通过求解该方程，可以得到煤柱在给定条件下的极限承载能力。如果实际作用在煤柱上的载荷小于极限承载能力，则认为煤柱是稳定的。极限平衡法的优点是计算简单、直观，物理意义明确，在工程实践中应用较为广泛。然而，该方法也存在一定的局限性。它假设煤柱是均质、连续的材料，忽略了煤岩的非均质性和内部结构的复杂性，同时也没有考虑煤岩的蠕变特性和时间因素对煤柱稳定性的影响。因此，极限平衡法适用于对条带煤柱稳定性进行初步评估，对于深部开采等复杂条件下的煤柱稳定性评价，其准确性可能受到一定影响。数值模拟法：数值模拟法是利用计算机技术对条带煤柱的开采过程进行模拟分析，从而评估其稳定性的方法。常用的数值模拟软件有FLAC3D、ANSYS、UDEC等。在数值模拟中，首先需要建立条带煤柱和周围岩体的三维模型，定义煤岩的材料参数、边界条件和开采过程。然后，通过数值计算求解模型的力学响应，得到煤柱的位移、应力、应变、塑性区范围等参数。以FLAC3D软件为例，在建立条带煤柱模型时，将煤岩视为连续介质，采用合适的本构模型来描述其力学行为，如摩尔-库仑本构模型、霍克-布朗本构模型等。通过设置边界条件和加载过程，模拟开采过程中煤柱的受力和变形情况。数值模拟法的优点是能够考虑多种复杂因素对煤柱稳定性的影响，如煤岩的非线性力学特性、蠕变特性、开采顺序、地质构造等。通过数值模拟，可以直观地观察到煤柱在开采过程中的力学响应和破坏过程，为煤柱稳定性评价提供详细的信息。然而，数值模拟结果的准确性依赖于模型的合理性和参数的选取。如果模型建立不合理或参数选取不准确，可能会导致模拟结果与实际情况相差较大。数值模拟计算量较大，需要较高的计算机硬件配置和较长的计算时间。因此，数值模拟法适用于对条带煤柱稳定性进行深入研究和详细分析，在实际工程应用中，需要结合现场实测数据对模拟结果进行验证和修正。可靠性分析法：可靠性分析法是基于概率统计理论，考虑煤柱稳定性影响因素的不确定性，对条带煤柱的稳定性进行评估的方法。煤柱的稳定性受到多种因素的影响，如煤岩的力学参数、开采条件、地质构造等，这些因素往往具有不确定性。可靠性分析法通过建立煤柱稳定性的功能函数，并考虑各影响因素的概率分布，计算煤柱的失效概率或可靠度指标。假设煤柱的稳定性功能函数为Z=R-S，其中R为煤柱的承载能力，S为煤柱所承受的载荷。如果Z>0，则煤柱处于稳定状态；如果Z<0，则煤柱发生失效。通过对R和S的概率分布进行分析，可以计算出煤柱的失效概率P_f和可靠度指标\beta。可靠性分析法的优点是能够充分考虑煤柱稳定性影响因素的不确定性，更加客观地评估煤柱的稳定性。它可以为煤柱的设计和开采提供概率意义上的指导，有助于合理确定煤柱的尺寸和安全系数。然而，可靠性分析法需要大量的现场实测数据和统计分析，以确定各影响因素的概率分布，这在实际工程中往往具有一定的难度。该方法的计算过程相对复杂，需要具备一定的概率统计知识。因此，可靠性分析法适用于对条带煤柱稳定性要求较高、对不确定性因素较为关注的工程情况，在实际应用中，需要结合其他评价方法进行综合分析。综上所述，极限平衡法、数值模拟法和可靠性分析法各有优缺点，在条带煤柱长期稳定性评价中，应根据具体情况选择合适的评价方法，或者将多种方法相结合，以提高评价结果的准确性和可靠性。4.3条带煤柱长期稳定性预测模型4.3.1基于蠕变模型的稳定性预测根据煤岩蠕变试验结果，选择合适的蠕变模型是建立条带煤柱长期稳定性预测模型的关键。在众多的蠕变模型中，Burgers模型是一种较为常用且能够较好描述煤岩蠕变特性的模型。Burgers模型由一个Maxwell模型和一个Kelvin模型串联组成，能够综合反映煤岩的瞬时弹性变形、黏弹性变形和黏塑性变形。Maxwell模型由一个弹簧和一个黏壶串联而成，其本构方程为：\dot{\varepsilon}=\frac{\dot{\sigma}}{E_1}+\frac{\sigma}{\eta_1}其中，\varepsilon为应变，\sigma为应力，E_1为弹性模量，\eta_1为黏性系数，\dot{\varepsilon}和\dot{\sigma}分别为应变率和应力率。该模型主要描述煤岩的瞬时弹性变形和黏滞流动，在加载瞬间，弹簧立即产生弹性变形，随后黏壶以恒定的速率产生黏滞流动。Kelvin模型由一个弹簧和一个黏壶并联而成，其本构方程为：\sigma=E_2\varepsilon+\eta_2\dot{\varepsilon}其中，E_2为弹性模量，\eta_2为黏性系数。该模型主要描述煤岩的黏弹性变形，在加载后，弹簧和黏壶共同作用，应变逐渐增加，当应力保持不变时，应变最终趋于一个稳定值。Burgers模型将Maxwell模型和Kelvin模型串联起来，其本构方程为：\dot{\varepsilon}=\frac{\dot{\sigma}}{E_1}+\frac{\sigma}{\eta_1}+\frac{1}{E_2}\left(\sigma-\eta_2\dot{\varepsilon}\right)通过对Burgers模型的本构方程进行积分求解，可以得到在恒定应力\sigma_0作用下的应变-时间关系：\varepsilon(t)=\frac{\sigma_0}{E_1}+\frac{\sigma_0t}{\eta_1}+\frac{\sigma_0}{E_2}\left(1-e^{-\frac{E_2t}{\eta_2}}\right)其中，t为时间。该式中，\frac{\sigma_0}{E_1}表示瞬时弹性应变，\frac{\sigma_0t}{\eta_1}表示黏滞流动产生的应变，\frac{\sigma_0}{E_2}\left(1-e^{-\frac{E_2t}{\eta_2}}\right)表示黏弹性应变。在建立条带煤柱长期稳定性预测模型时，将Burgers模型与条带煤柱的力学模型相结合。考虑条带煤柱在开采过程中所承受的载荷P，根据煤柱的受力分析，得到煤柱内部的应力分布。假设煤柱内部某点的应力为\sigma，将\sigma代入Burgers模型的应变-时间关系式中，即可得到该点的应变随时间的变化规律。通过对煤柱内各点的应变进行分析，可以预测煤柱的变形和稳定性。例如，对于一个宽度为b，高度为h的条带煤柱，上覆岩层施加的均布载荷为q。根据弹性力学理论，煤柱内某点的垂直应力\sigma_y可以表示为：\sigma_y=\frac{q(b+a)}{b}其中，a为采空区宽度。将\sigma_y代入Burgers模型的应变-时间关系式中，得到该点的垂直应变\varepsilon_y(t)随时间的变化规律。通过对煤柱不同位置处的\varepsilon_y(t)进行分析，可以了解煤柱的变形情况。当煤柱的变形超过一定限度时，认为煤柱处于不稳定状态。4.3.2模型验证与分析利用实际工程数据对基于Burgers模型建立的条带煤柱长期稳定性预测模型进行验证。选取某深井条带煤柱开采工程作为研究对象，该工程的开采深度为1200m，煤柱宽度为30m，采空区宽度为50m。在开采过程中，对条带煤柱的位移、应力等参数进行了长期监测。将实际工程的地质条件和开采参数代入预测模型中，计算得到条带煤柱在不同时间的位移和应力分布。将计算结果与现场监测数据进行对比分析。从位移对比结果来看，在开采初期，预测模型计算得到的煤柱位移与现场监测数据基本吻合，随着开采时间的延长，虽然两者之间存在一定的偏差，但变化趋势一致。在开采1年后，预测模型计算得到的煤柱垂直位移为50mm，现场监测数据为55mm，相对误差约为9.1%。从应力对比结果来看，预测模型计算得到的煤柱内部应力分布与现场实测应力分布也具有较好的一致性。在煤柱与采空区交界处，预测模型计算得到的垂直应力为15MPa，现场实测应力为16MPa，相对误差约为6.25%。通过对预测模型计算结果与现场监测数据的对比分析可知，该预测模型能够较好地反映条带煤柱的长期稳定性变化规律，具有五、工程实例分析5.1工程概况某深井煤矿位于[具体地理位置]，该区域地质构造较为复杂，地层呈现出多褶皱、断裂的特征。矿井开采深度达到1200m，属于典型的深井开采。井田内主要可采煤层为3号煤层，煤层厚度平均为4.5m，倾角在10°-15°之间。煤层顶板主要为砂质泥岩，厚度约为5m，其抗压强度为30-40MPa，抗拉强度为2-3MPa；底板为粉砂岩，厚度约为4m，抗压强度为35-45MPa，抗拉强度为2.5-3.5MPa。该煤矿采用条带开采方案，旨在有效控制上覆岩层的移动和变形，减少对地表建筑物的影响，同时提高煤炭资源的回收率。条带开采的基本参数如下：采出条带宽度设计为30m，煤柱宽度设计为40m。采出条带和煤柱沿煤层走向交替布置。在开采过程中，采用综采工艺，配备先进的采煤设备，如大功率采煤机、刮板输送机和液压支架等。采煤机的割煤速度根据煤层硬度和开采条件进行调整，一般控制在3-5m/min；刮板输送机的运输能力为1000-1200t/h，以满足采煤的出煤需求；液压支架的工作阻力为6000-8000kN，能够有效支撑顶板，保障采煤工作的安全进行。5.2煤岩蠕变试验结果分析为了深入研究该煤矿煤岩的蠕变特性，在实验室进行了系统的煤岩蠕变试验。从井下采集具有代表性的煤样后，将其加工成标准试件，试件尺寸严格按照直径50mm，高度100mm进行制备。试验设备选用高精度的MTS伺服控制岩石三轴试验机，确保试验数据的准确性和可靠性。在单轴蠕变试验中，采用分级加载方式。首先施加5MPa的初始应力，保持该应力水平，每10分钟记录一次试件的轴向变形。随着时间的推移，轴向应变迅速增加，在最初的2小时内，应变增加较为明显，随后应变增加速率逐渐减慢，呈现出典型的初始蠕变阶段特征。当应力增加到10MPa时，重复上述过程，此时轴向应变的增加幅度和速率均大于5MPa应力阶段。通过对不同应力水平下的单轴蠕变试验数据进行分析，得到煤样在单轴状态下的蠕变曲线，发现随着应力的增大，煤样的稳态蠕变阶段应变率显著增大，加速蠕变阶段提前到来。在三轴蠕变试验中，设置了不同的围压和轴压组合。当围压为5MPa，轴压分别为10MPa、15MPa、20MPa时，进行蠕变试验。试验结果表明，在相同围压下，随着轴压的增加，煤样的初始蠕变阶段应变增加幅度增大，稳态蠕变阶段的轴向应变率和径向应变率均增大。当轴压为15MPa，围压从5MPa增加到10MPa时，煤样的轴向应变和径向应变在相同时间内明显减小。这充分说明围压对煤岩的蠕变具有显著的抑制作用，围压的增大能够限制煤岩内部微裂隙的扩展，增强煤岩的抵抗变形能力。通过对试验数据的进一步分析，确定了该煤矿煤岩的蠕变参数。在Burgers模型中，弹性模量E_1为10GPa，黏性系数\eta_1为500GPa・h，弹性模量E_2为5GPa，黏性系数\eta_2为300GPa・h。煤岩的长期强度通过试验确定为25MPa，当应力超过该值时，煤岩将进入加速蠕变阶段，变形急剧增大，直至破坏。5.3条带煤柱长期稳定性评价运用前文所述的条带煤柱长期稳定性评价方法，对该煤矿的条带煤柱长期稳定性进行全面评价。首先，采用极限平衡法进行初步评估。根据煤柱的实际尺寸、煤岩的力学参数以及上覆岩层的载荷等条件，利用Mohr-Coulomb屈服准则建立极限平衡方程。计算得到在当前开采条件下，煤柱的极限承载能力为35MPa，而实际作用在煤柱上的平均应力为20MPa。通过计算安全系数，安全系数为1.75，大于1，初步判断煤柱在当前状态下处于稳定状态。利用数值模拟软件FLAC3D建立详细的三维数值模型。模型中充分考虑煤岩的非线性力学特性、蠕变特性以及开采过程中的各种复杂因素。通过模拟条带煤柱在开采过程中的受力和变形情况，得到煤柱的位移、应力、应变和塑性区范围等参数。模拟结果显示，煤柱的最大垂直位移出现在煤柱顶部，约为50mm；煤柱内部的最大应力为22MPa，位于煤柱与采空区交界处；塑性区主要分布在煤柱的边缘部分，塑性区范围占煤柱宽度的20%。从模拟结果来看，煤柱目前处于相对稳定状态，但需要密切关注塑性区的发展情况，防止塑性区进一步扩大导致煤柱失稳。采用可靠性分析法，考虑煤柱稳定性影响因素的不确定性。通过对煤岩力学参数、开采条件等因素进行统计分析，确定其概率分布。建立煤柱稳定性的功能函数，计算得到煤柱的失效概率为0.05，可靠度指标为2.0。这表明煤柱在当前条件下具有较高的可靠性，但仍存在一定的失效风险。综合以上三种评价方法的结果，虽然目前条带煤柱处于相对稳定状态，但由于煤岩的蠕变特性以及开采过程中的不确定性因素，煤柱存在一定的潜在风险。随着开采时间的延长，煤岩蠕变可能导致煤柱强度降低，塑性区范围扩大，从而影响煤柱的长期稳定性。5.4稳定性控制措施与效果分析为了有效提高条带煤柱的长期稳定性，降低潜在风险，提出了一系列针对性的稳定性控制措施。在优化煤柱尺寸方面，通过数值模拟和理论分析，对煤柱宽度进行了优化调整。将原设计的40m煤柱宽度增加到45m，经过重新计算和模拟，煤柱的承载能力得到提高，内部应力分布更加均匀，塑性区范围减小到煤柱宽度的15%。这表明适当增加煤柱宽度能够有效增强煤柱的稳定性。加强支护措施也是提高煤柱稳定性的重要手段。在煤柱与采空区交界处，采用锚索和锚杆联合支护的方式。锚索长度为8m，间距为1.5m；锚杆长度为2.5m，间距为1.0m。通过支护，煤柱边缘的位移得到有效控制，最大垂直位移减小到30mm，煤柱与采空区交界处的应力集中现象得到缓解，应力降低了20%。这说明加强支护能够提高煤柱的稳定性，减少煤柱的变形和破坏。在开采工艺方面，优化开采顺序，采用跳采的方式。先开采间隔的采出条带，待煤柱稳定后再开采剩余的采出条带。通过跳采，煤柱的受力更加均匀，避免了连续开采导致的应力集中。数值模拟结果显示，采用跳采方式后，煤柱内部的最大应力降低到20MPa，塑性区范围进一步减小。这表明优化开采工艺对提高煤柱稳定性具有积极作用。通过实施以上稳定性控制措施，对条带煤柱的稳定性进行了再次监测和评估。结果表明，煤柱的各项稳定性指标得到明显改善，位移、应力和塑性区范围等均在安全范围内。这充分说明所提出的稳定性控制措施是有效的，能够显著提高条带煤柱的长期稳定性，为矿井的安全生产提供了有力保障。在后续的开采过程中，仍需持续对煤柱的稳定性进行监测，及时调整控制措施，确保矿井的安全高效开采。六、结论与展望6.1研究成果总结本研究通过开展煤岩蠕变试验，深入研究了深井条带煤柱的长期稳定性，取得了以下主要研究成果：煤岩蠕变特性研究：通过不同条件下的煤岩蠕变试验，包括单轴蠕变试验、三轴蠕变试验以及不同环境因素作用下的蠕变试验，系统地分析了煤岩在不同应力水平和环境条件下的蠕变特性。试验结果表明，煤岩的蠕变过程具有明显的阶段性，随着应力的增加，煤岩的蠕变变形速率加快，稳态蠕变阶段缩短，加速蠕变阶段提前到来。围压的增大可以抑制煤岩的蠕变变形，提高煤岩的蠕变稳定性。水、温度等环境因素对煤岩蠕变特性也有显著影响，饱水煤岩的蠕变变形明显大于干燥煤岩，温度升高会使煤岩的蠕变变形增大，蠕变速率加快。运用微观测