基于多维度分析的荣华矿深部-650水平大巷围岩吸水软化特性研究

基于多维度分析的荣华矿深部-650水平大巷围岩吸水软化特性研究一、引言1.1研究背景与意义随着煤炭资源开采向深部拓展，深部矿井巷道建设面临着诸多严峻挑战。荣华矿作为煤炭生产的重要基地，在深部开采进程中，-650水平大巷的建设遭遇了围岩吸水软化这一棘手难题。岩石的吸水软化特性会导致其物理力学性质发生显著劣化，进而严重威胁巷道的稳定性与安全生产。从物理性质方面来看，岩石吸水后，其内部孔隙结构会被水填充，导致孔隙压力增大，使得岩石颗粒间的有效应力减小，从而引发岩石的膨胀变形。这种膨胀变形若不能得到有效控制，会致使巷道断面收缩，影响通风、运输等正常生产活动。相关研究表明，部分软岩在吸水饱和后，其体积膨胀率可达5%-10%，这对于巷道的空间稳定性构成了极大威胁。在力学性质层面，吸水软化会使岩石的强度大幅降低。以砂岩为例，吸水后的砂岩单轴抗压强度可能会降低30%-50%，抗剪强度也会随之显著下降。这使得巷道围岩在原岩应力及开采扰动应力的作用下，更容易发生破坏，如出现片帮、冒顶等事故。据统计，在因围岩问题导致的矿井事故中，约有30%-40%与岩石的吸水软化特性密切相关。荣华矿-650水平大巷处于复杂的地质环境中，地下水丰富且水压较高，这为岩石的吸水提供了有利条件。同时，该水平大巷的服务年限较长，对其稳定性要求极高。然而，围岩的吸水软化现象使得巷道支护难度急剧增加，常规的支护方式难以满足工程需求，不仅增加了支护成本，还影响了矿井的正常生产进度。对荣华矿深部-650水平大巷围岩吸水软化特性展开深入研究，具有极为重要的现实意义。通过研究，可以深入揭示该区域岩石吸水软化的内在机制，为巷道支护设计提供更为科学、准确的理论依据。例如，根据岩石的吸水特性及软化规律，合理选择支护材料和支护方式，提高支护结构的适应性和承载能力，从而有效保障巷道的长期稳定性。这不仅有助于减少巷道维护成本，降低安全事故发生的概率，还能确保矿井的安全生产，提高煤炭资源的开采效率，促进矿山企业的可持续发展。1.2国内外研究现状岩石吸水软化特性的研究一直是岩土工程领域的重要课题，国内外学者围绕这一主题开展了大量研究，涵盖岩石吸水特性、软化特性及微观机制等多个方面。在岩石吸水特性研究方面，国外起步较早。早在20世纪中叶，一些学者就开始关注岩石与水的相互作用，通过实验观察岩石在不同水环境下的质量变化，初步探究岩石的吸水规律。随着技术的发展，高精度的测量仪器被应用于岩石吸水实验，如采用高精度电子天平实时监测岩石吸水过程中的质量变化，使得对吸水特性的研究更加精确。国内学者在这方面也进行了深入探索。例如，有学者利用深部软岩吸附智能测试系统对砾岩进行气态水吸附实验，建立了吸水过程函数，研究了岩石吸附气态水的吸水特性，发现软岩在高温高湿环境下不仅会吸附液态水，还会从空气中吸附气态水，且其吸水特性与岩石的孔隙结构、矿物成分等密切相关。对于岩石吸水后的软化特性研究，国外学者通过大量的室内实验和现场监测，分析了岩石在不同吸水率条件下的力学性能变化。他们采用先进的材料测试设备，如电液伺服万能试验机，对不同吸水状态的岩石进行力学测试，获取岩石的抗压强度、抗拉强度、弹性模量等力学参数，揭示了岩石强度随吸水率增加而降低的规律。国内研究人员也开展了丰富的研究工作。以荣华矿深部-650水平大巷围岩中的粗砂岩、细砂岩和泥质粉砂岩为研究对象，通过对三组岩石进行不同吸水率条件下的单轴压缩试验，得出粗砂岩和细砂岩组岩石吸水率与单轴抗压强度特征曲线服从幂函数分布，而泥质粉砂岩组岩石服从线性函数分布；三组岩石吸水率与弹性模量特征曲线服从幂函数分布；粗砂岩和细砂岩组岩石吸水率与次软化系数特征曲线服从幂函数分布，泥质粉砂岩组则服从线性函数分布。在岩石吸水软化微观机制研究方面，国外主要借助先进的微观观测技术，如扫描电子显微镜（SEM）、压汞仪（MIP）等，从微观层面分析岩石吸水后内部结构和矿物成分的变化，揭示岩石吸水软化的内在机理。国内学者同样利用这些技术，对岩石微观结构进行深入研究。通过X射线衍射试验、扫描电镜试验，分析岩石的黏土矿物成分及微观结构特征，发现岩石的黏土矿物类型和矿物晶体的微观结构是影响岩石吸水的主要因素，吸水后岩石内部孔隙结构发生变化，黏土矿物膨胀，导致岩石力学性能劣化。现有研究仍存在一定不足。一方面，在研究岩石吸水软化特性时，多集中于单一因素的影响，如仅考虑吸水率对岩石力学性能的影响，而对多种因素耦合作用的研究较少。实际上，在深部矿井环境中，岩石不仅受到水的作用，还受到高地应力、高温等因素的共同影响，这些因素的耦合作用可能会对岩石的吸水软化特性产生更为复杂的影响。另一方面，目前的研究成果在实际工程应用中还存在一定的局限性。虽然对岩石吸水软化特性有了一定的认识，但在如何将这些理论成果有效地应用于深部巷道支护设计和施工方面，还缺乏系统的研究和实践经验。例如，如何根据岩石的吸水软化特性选择合适的支护材料和支护时机，如何在施工过程中采取有效的防水、排水措施以减少岩石的吸水等问题，都有待进一步深入研究。本研究将以荣华矿深部-650水平大巷围岩为具体研究对象，针对现有研究的不足，开展多因素耦合作用下岩石吸水软化特性的试验研究，分析岩石在不同工况下的物理力学性质变化规律，并将研究成果应用于巷道支护设计优化，为解决深部矿井巷道围岩支护难题提供新的思路和方法。1.3研究内容与方法1.3.1研究内容岩石吸水特性研究：在荣华矿深部-650水平大巷围岩中，选取具有代表性的粗砂岩、细砂岩和泥质粉砂岩等岩石样本。运用深部软岩吸附智能测试系统等先进设备，开展岩石吸水实验，不仅研究其对液态水的吸附情况，还探究在高温高湿环境下对气态水的吸附特性。通过实时监测岩石在不同环境条件下的质量变化、体积变化等参数，建立准确的岩石吸水过程函数，深入分析影响岩石吸水特性的关键因素，如岩石的孔隙结构、矿物成分等。岩石吸水软化特性研究：对不同吸水率条件下的岩石样本进行单轴压缩试验、三轴压缩试验等力学测试。利用电液伺服万能试验机等高精度设备，精确测量岩石的抗压强度、抗拉强度、弹性模量、泊松比等力学参数。通过分析这些力学参数随吸水率的变化规律，建立岩石吸水率与力学参数之间的定量关系模型，明确岩石在吸水过程中的软化机制和强度劣化规律。岩石吸水软化微观机制研究：采用扫描电子显微镜（SEM）观察岩石吸水前后内部微观结构的变化，如孔隙大小、形状、连通性的改变，以及矿物颗粒的排列方式和胶结情况的变化。运用压汞仪（MIP）测定岩石的孔隙分布特征，结合X射线衍射（XRD）分析岩石的矿物成分，从微观层面揭示岩石吸水软化的内在机理，明确黏土矿物等成分在吸水软化过程中的作用机制。考虑吸水软化特性的巷道稳定性分析：基于岩石吸水软化特性的研究成果，运用数值模拟软件，如FLAC3D、ANSYS等，建立荣华矿-650水平大巷的三维数值模型。模拟巷道在不同围岩条件、支护方式和地下水作用下的开挖过程，分析巷道围岩的应力、应变分布规律，预测巷道可能出现的变形和破坏形式，评估岩石吸水软化对巷道稳定性的影响程度。基于研究成果的巷道支护优化设计：根据岩石吸水软化特性和巷道稳定性分析结果，结合荣华矿的实际工程条件，对-650水平大巷的支护方案进行优化设计。从支护材料的选择、支护结构的形式、支护参数的确定等方面入手，提出针对性的支护措施，提高支护结构的承载能力和适应性，确保巷道的长期稳定。1.3.2研究方法室内实验法：通过现场采样，获取荣华矿深部-650水平大巷围岩的岩石样本。在实验室中，利用岩石切割机、磨片机等设备将岩石样本加工成标准尺寸的试件。运用深部软岩吸附智能测试系统进行岩石吸水实验，采用电液伺服万能试验机进行岩石力学性能测试，借助扫描电子显微镜、压汞仪、X射线衍射仪等微观测试设备分析岩石的微观结构和矿物成分。理论分析法：运用岩石力学、材料力学、渗流力学等相关理论，分析岩石吸水过程中的水分迁移规律、力学性能劣化机制。建立岩石吸水软化的理论模型，推导岩石吸水率与力学参数之间的理论关系，为实验研究和数值模拟提供理论依据。数值模拟法：利用FLAC3D、ANSYS等数值模拟软件，建立考虑岩石吸水软化特性的巷道三维数值模型。通过设置不同的材料参数、边界条件和荷载工况，模拟巷道在开挖过程中的力学响应，分析围岩的应力、应变分布和变形破坏特征，预测巷道的稳定性。工程类比法：收集国内外类似地质条件下深部巷道建设的工程案例，对比分析其在应对岩石吸水软化问题时所采用的技术措施和工程经验。结合荣华矿的实际情况，借鉴成功经验，为-650水平大巷的支护设计和施工提供参考。二、荣华矿工程地质概况2.1矿区地质构造荣华矿所处区域地质构造较为复杂，历经多期构造运动的叠加影响，褶皱与断层广泛发育，对矿区的地质结构和-650水平大巷围岩稳定性产生了显著作用。从褶皱构造来看，矿区内存在一系列轴向呈北东-南西向的褶皱。这些褶皱规模大小不一，形态各异，包括紧闭褶皱和开阔褶皱。其中，规模较大的褶皱轴向延伸可达数千米，其核部和翼部的岩石受力状态差异明显。在褶皱核部，岩石受到强烈的挤压作用，岩层发生弯曲变形，层间错动频繁，导致岩石破碎，完整性遭到严重破坏。例如，某大型褶皱核部的岩石节理裂隙极为发育，岩石破碎程度高，呈现出碎块状结构，这使得该区域岩石的力学强度大幅降低，为地下水的运移和储存提供了良好的通道和空间。当-650水平大巷穿越褶皱核部时，围岩在地下水的作用下更容易发生吸水软化现象，进而影响巷道的稳定性。而在褶皱翼部，岩石虽然相对完整，但由于受到褶皱作用的影响，岩层的产状发生改变，使得岩石的力学性质具有明显的各向异性。这种各向异性会导致围岩在不同方向上的变形和破坏特征存在差异，增加了巷道支护的难度。在断层分布方面，矿区内主要发育有北东向和北西向两组断层。北东向断层规模较大，延伸长度可达数千米，断距较大，对矿区地层的错动和破坏较为明显。北西向断层规模相对较小，但数量较多，常与北东向断层相互切割，形成复杂的断裂网络。这些断层的存在破坏了岩石的连续性和完整性，使得断层附近的岩石破碎，形成断层破碎带。断层破碎带内的岩石多呈碎裂状或糜棱状，结构松散，孔隙率大。荣华矿-650水平大巷在掘进过程中，若遇到断层破碎带，围岩的稳定性将面临严峻挑战。由于断层破碎带内岩石的破碎程度高，地下水容易在此积聚，岩石与水的接触面积增大，吸水软化的速度加快。而且，断层破碎带的力学强度极低，难以承受巷道开挖后的围岩应力，容易引发巷道的坍塌、片帮等事故。地质构造对-650水平大巷围岩稳定性的潜在影响是多方面的。在巷道开挖过程中，褶皱和断层会导致围岩应力重新分布。在褶皱核部和断层附近，由于岩石的完整性被破坏，应力集中现象明显。当应力超过围岩的承载能力时，围岩就会发生破坏。同时，地质构造还会影响地下水的分布和运移。断层破碎带和褶皱核部的岩石破碎，孔隙和裂隙发育，成为地下水的良好通道和储存场所。地下水的存在不仅会使岩石发生吸水软化，降低岩石的力学强度，还会增加围岩的水压力，进一步恶化围岩的受力条件。此外，地质构造的复杂性也会增加巷道支护的难度和不确定性。在设计巷道支护方案时，需要充分考虑地质构造的影响，针对不同的地质构造区域采取相应的支护措施，以确保巷道的稳定性。2.2-650水平大巷围岩岩性特征荣华矿深部-650水平大巷围岩主要由粗砂岩、细砂岩和泥质粉砂岩等岩石组成，各岩石的岩性特征存在明显差异。粗砂岩颜色多呈灰白色或浅灰色，其矿物成分主要为石英，含量可达70%-80%，此外还含有少量的长石、云母等矿物。石英颗粒质地坚硬，是粗砂岩强度的主要贡献者。长石矿物的存在，在一定程度上影响了粗砂岩的化学稳定性，当岩石与水接触时，长石可能会发生水解等化学反应，从而对岩石的结构和性质产生影响。粗砂岩具有明显的粗粒结构，颗粒粒径一般在0.5-2mm之间，颗粒间的胶结方式主要为硅质胶结和钙质胶结。硅质胶结使得岩石的强度较高，结构较为致密；而钙质胶结在遇到酸性水时，可能会发生溶解，导致岩石胶结强度降低。在构造方面，粗砂岩多呈厚层状构造，层理较为清晰，这种层状构造使得岩石在力学性质上具有一定的各向异性，沿层理方向和垂直层理方向的强度和变形特性有所不同。细砂岩通常呈现出灰绿色或浅黄色，矿物成分同样以石英为主，含量约为60%-70%，长石含量相对粗砂岩略高，在20%-30%左右，还含有少量的黏土矿物和暗色矿物。细砂岩的结构为细粒结构，颗粒粒径在0.1-0.5mm之间，颗粒大小相对均匀。其胶结类型主要为泥质胶结和铁质胶结。泥质胶结的细砂岩强度相对较低，遇水后泥质矿物容易软化，导致岩石整体强度下降；铁质胶结则使岩石具有一定的抗风化能力，但在地下水的长期作用下，铁质可能会发生氧化，影响岩石的结构。细砂岩的构造一般为薄层状构造或透镜状构造，层理发育，薄层状构造使得细砂岩在受力时更容易沿着层理面发生滑动或错动，影响其稳定性。泥质粉砂岩颜色以深灰色或灰黑色为主，矿物成分中黏土矿物含量较高，可达40%-50%，此外还含有较多的石英、长石等碎屑矿物。黏土矿物具有较强的吸水性和膨胀性，是影响泥质粉砂岩水理性质的关键因素。泥质粉砂岩具有粉砂质结构，颗粒粒径细小，多在0.01-0.1mm之间。其胶结方式主要为泥质胶结，这种胶结方式使得岩石的强度较低，结构较为松散。在构造上，泥质粉砂岩常具有页理构造，页理面是岩石中的薄弱面，在水和外力的作用下，岩石容易沿着页理面发生劈裂和剥落。这些不同岩性的岩石在-650水平大巷围岩中相互组合，形成了复杂的地质结构。由于各岩石的矿物成分、结构构造不同，其物理力学性质和水理性质也存在较大差异，这使得大巷围岩在受到地下水作用时，不同岩石的吸水软化特性各不相同，进一步增加了巷道支护的复杂性和难度。例如，泥质粉砂岩中的黏土矿物遇水膨胀，可能会对相邻的粗砂岩和细砂岩产生挤压作用，改变围岩的应力分布状态，从而影响整个巷道的稳定性。2.3地下水赋存状况荣华矿深部-650水平大巷所处区域地下水赋存状况较为复杂，对大巷围岩的稳定性有着显著影响。该区域地下水主要类型为孔隙水、裂隙水和岩溶水，不同类型的地下水在水位、水量、水质等方面存在差异，且与大巷围岩之间存在紧密的水力联系。在水位方面，通过长期的水文观测资料可知，该区域孔隙水水位埋深一般在20-50m之间，水位随季节变化较为明显，在雨季时水位会有所上升，而在旱季则会相对下降，变化幅度可达5-10m。这是由于大气降水是孔隙水的主要补给来源，雨季降水充沛，大量雨水渗入地下，使得孔隙水水位升高；旱季降水减少，孔隙水主要通过蒸发和向深部含水层排泄等方式消耗，导致水位下降。裂隙水水位受地质构造控制较为明显，在断层破碎带和褶皱核部等区域，裂隙发育，裂隙水水位相对较低，一般埋深在50-80m之间；而在构造相对完整的区域，裂隙水水位埋深约为30-50m。岩溶水水位则相对稳定，埋深一般在80-100m左右，其水位变化主要受区域岩溶水系统的补给和排泄条件影响。水量方面，孔隙水的水量相对较小，其含水层的渗透系数一般在0.1-1m/d之间，单井涌水量多在5-10m³/d。这是因为孔隙水主要赋存于岩石的孔隙中，孔隙的连通性相对较差，导致其导水能力有限。裂隙水的水量变化较大，在裂隙发育良好的区域，其含水层渗透系数可达1-10m/d，单井涌水量能达到10-50m³/d；而在裂隙不发育的区域，水量则相对较小。岩溶水由于其赋存空间大，连通性好，含水层渗透系数通常在10-100m/d之间，单井涌水量较大，可达50-100m³/d，甚至更大。水质分析结果表明，孔隙水的矿化度较低，一般在0.5-1g/L之间，主要离子成分包括钙离子、镁离子、钠离子、氯离子和硫酸根离子等，水质类型多为HCO₃-Ca型或HCO₃-Ca・Mg型，呈弱碱性，pH值在7.5-8.5之间。这种水质对大巷围岩中的部分矿物可能会产生溶解作用，尤其是对一些易溶于水的矿物，如方解石等，从而影响围岩的结构和强度。裂隙水矿化度相对较高，在1-3g/L之间，除了常见离子外，还可能含有一些微量元素，如铁、锰等，水质类型较为复杂，可能为SO₄-Ca型、Cl-Na型等。岩溶水矿化度一般在3-5g/L之间，由于其在岩溶管道中流动过程中与岩石发生化学反应，水中含有较多的碳酸氢根离子和钙离子，水质类型多为HCO₃-Ca型，具有一定的腐蚀性。大巷围岩与地下水之间存在着密切的水力联系。围岩中的孔隙和裂隙为地下水的运移提供了通道，地下水在运移过程中与围岩发生相互作用。一方面，地下水会沿着围岩的孔隙和裂隙渗入，使岩石发生吸水软化，降低岩石的力学强度。例如，泥质粉砂岩中的黏土矿物遇水后会发生膨胀，导致岩石结构松散，强度降低。另一方面，当巷道开挖后，改变了围岩的应力状态和地下水的渗流场，地下水会向巷道周围汇集，增加巷道涌水量，进一步恶化围岩的稳定性。在断层破碎带附近，由于岩石破碎，裂隙连通性好，地下水与围岩的水力联系更为密切，岩石吸水软化的速度更快，对巷道稳定性的威胁更大。三、试验方案设计3.1试验材料选取为全面、准确地研究荣华矿深部-650水平大巷围岩的吸水软化特性，在材料选取上，从-650水平大巷的不同位置进行岩样采集。大巷的不同位置，由于受到地质构造、地下水分布等因素的影响，岩石的特性可能存在差异，因此多点位采样能够保证样品具有广泛的代表性。在粗砂岩采样过程中，使用专业的取芯设备，如金刚石钻头取芯机，在大巷的顶部、侧壁等不同部位进行钻孔取芯。在顶部选取了3个采样点，侧壁选取了5个采样点，确保采集的岩芯能够反映粗砂岩在不同受力环境下的特性。每个采样点采集的岩芯长度不小于30cm，直径为5cm，以满足后续加工成标准试件的需求。采集到的粗砂岩岩芯表面较为粗糙，颜色以灰白色为主，夹杂少量深色矿物颗粒，颗粒粒径较大，在0.5-2mm之间，具有明显的粗粒结构，层理较为清晰。对于细砂岩，同样在大巷的不同位置进行采样。考虑到细砂岩在大巷中的分布特点，在大巷的弯道处、直道处分别选取了采样点，因为弯道处的岩石受到的应力集中作用与直道处不同，可能导致岩石的结构和性质有所差异。共采集了8个细砂岩岩样，岩样呈灰绿色，质地相对细腻，颗粒粒径在0.1-0.5mm之间，结构均匀，具有薄层状构造，层理发育。泥质粉砂岩的采样则重点关注了大巷中靠近含水层的区域，以及岩石节理裂隙发育的部位。在靠近含水层的区域采集了4个岩样，在节理裂隙发育部位采集了3个岩样。这些区域的泥质粉砂岩与地下水的接触机会更多，更容易受到水的作用而发生性质变化。采集到的泥质粉砂岩颜色较深，多为深灰色或灰黑色，具有粉砂质结构，颗粒粒径细小，多在0.01-0.1mm之间，页理构造明显。采集后的岩样及时进行编号和记录，记录内容包括采样位置、岩样特征、采样时间等信息。为防止岩样在运输和储存过程中受到损坏和环境因素的影响，将岩样用塑料薄膜包裹，并放置在专门的岩样箱中，箱内填充泡沫等缓冲材料，确保岩样的完整性，以便后续进行各项试验研究。3.2试验设备与仪器X射线衍射仪（XRD）：型号为D8Advance，由德国布鲁克公司生产。该设备在本次研究中发挥着至关重要的作用，主要用于分析岩石的矿物成分。其工作原理基于X射线与晶体物质相互作用产生的衍射现象。当X射线照射到岩石样品时，由于岩石中不同矿物晶体的晶格结构各异，会产生特定的衍射图案。通过对这些衍射图案的精确测量和分析，能够准确确定岩石中所含的矿物种类及其相对含量。例如，对于荣华矿-650水平大巷围岩中的粗砂岩，利用XRD分析发现其主要矿物成分为石英，含量约占70%-80%，还含有少量长石、云母等矿物；细砂岩中石英含量约为60%-70%，长石含量在20%-30%左右；泥质粉砂岩中黏土矿物含量较高，可达40%-50%。这些矿物成分的准确测定，为后续深入研究岩石的吸水软化特性提供了基础数据，因为不同矿物对水的吸附和反应特性不同，直接影响着岩石的吸水性能和软化机制。扫描电子显微镜（SEM）：选用日本日立公司的SU8010型扫描电镜。此设备主要用于观察岩石吸水前后内部微观结构的变化，其具有极高的分辨率，能够清晰呈现岩石内部微观结构的细节。在观察岩石微观结构时，首先需对岩石样品进行精细制备，将其切割成合适大小，并确保表面光滑，然后采用特殊的金属涂覆技术，如蒸镀或喷涂，在样品表面涂上一层薄薄的金属覆盖层，以增加样品的导电性，使电子束能够顺利通过样品并产生清晰的图像。通过SEM观察发现，吸水前的粗砂岩颗粒间胶结紧密，孔隙较小且连通性较差；吸水后，部分胶结物被水溶解，颗粒间的胶结力减弱，孔隙增大且连通性增强。泥质粉砂岩吸水后，黏土矿物发生膨胀，导致岩石内部微观结构变得更加松散，孔隙分布更加不均匀。这些微观结构变化的观察结果，有助于从微观层面深入理解岩石吸水软化的机制。岩石力学试验系统：采用美国MTS公司生产的Landmark电液伺服万能试验机，其最大载荷可达1000kN，具备高精度的力和位移控制功能，能够精确测定岩石的力学参数。在进行岩石单轴压缩试验时，将加工好的标准岩石试件放置在试验机的上下压板之间，通过试验机以恒定的加载速率对试件施加轴向压力，同时利用试验机配备的传感器实时采集试件在加载过程中的应力、应变数据。对于三轴压缩试验，则需在试件周围施加一定的围压，模拟岩石在地下的实际受力状态，然后再施加轴向压力进行加载测试。通过这些试验，能够获取不同吸水率条件下岩石的抗压强度、抗拉强度、弹性模量、泊松比等力学参数，从而分析岩石力学性能随吸水率的变化规律，为研究岩石吸水软化对巷道稳定性的影响提供重要的力学数据支持。深部软岩吸附智能测试系统：该系统为自主研发，专门用于研究岩石的吸水特性。它能够精确控制试验环境的温度、湿度等条件，实时监测岩石在吸水过程中的质量变化、体积变化等参数。在进行岩石吸水实验时，将岩样放置在测试系统的样品舱内，通过调节系统的环境参数，模拟不同的地下环境条件，如高温高湿环境等，然后启动测试系统，系统会自动记录岩样在不同时间点的质量和体积数据。利用该系统对荣华矿-650水平大巷围岩进行吸水实验，发现岩石不仅会吸附液态水，在高温高湿环境下还能从空气中吸附气态水，且其吸水速率和吸水量与环境的温度、湿度以及岩石的矿物成分、孔隙结构等因素密切相关。这些实验结果为深入研究岩石的吸水特性提供了关键数据。压汞仪（MIP）：型号为AutoPoreIV9500，由美国麦克默瑞提克公司生产。其主要用于测定岩石的孔隙分布特征，通过测量不同压力下汞注入岩石孔隙中的体积，来确定岩石孔隙的大小、孔径分布以及孔隙率等参数。在测试过程中，将干燥的岩石样品放入压汞仪的样品池中，逐渐增加汞的压力，使汞逐渐进入岩石的孔隙中。由于不同大小的孔隙需要不同的压力才能使汞进入，因此根据压力与汞注入体积的关系曲线，即可分析出岩石的孔隙分布情况。例如，通过MIP测试发现，粗砂岩的孔隙主要以中孔和大孔为主，孔隙率相对较低；泥质粉砂岩的孔隙则以微孔和介孔为主，孔隙率较高。这些孔隙分布特征的测定结果，对于理解岩石的吸水机制以及水分在岩石内部的运移规律具有重要意义。3.3试验方法与步骤3.3.1X射线衍射试验在进行X射线衍射试验时，首先对采集到的岩石样品进行精细制备。将岩石样品切割成厚度约为0.5-1mm的薄片，然后使用研磨设备将薄片表面研磨至光滑平整，以确保X射线能够均匀地穿透样品并产生清晰的衍射信号。把制备好的岩石薄片放置在X射线衍射仪的样品台上，调整样品的位置和角度，使其能够准确地接收X射线的照射。开启X射线衍射仪，设置合适的工作参数，如X射线管电压为40kV，电流为30mA，扫描范围为5°-80°，扫描速度为2°/min。在扫描过程中，X射线照射到岩石样品上，与岩石中的矿物晶体相互作用，产生衍射现象。探测器会记录下衍射信号的强度和角度信息，这些信息以衍射图谱的形式呈现。试验结束后，对采集到的衍射图谱进行数据处理和分析。利用专业的衍射分析软件，如MDIJade，对衍射图谱中的衍射峰进行识别和标定。通过与标准矿物衍射数据库进行比对，确定岩石中所含的矿物种类及其相对含量。例如，对于荣华矿-650水平大巷围岩中的粗砂岩，通过X射线衍射分析，确定其主要矿物成分为石英，含量约占70%-80%，还含有少量长石、云母等矿物。这些矿物成分的准确测定，为后续深入研究岩石的吸水软化特性提供了重要的基础数据，因为不同矿物对水的吸附和反应特性不同，直接影响着岩石的吸水性能和软化机制。3.3.2扫描电镜试验扫描电镜试验用于观察岩石微观结构特征，首先对岩石样品进行精心制备。将岩石切割成尺寸约为5mm×5mm×5mm的小块，确保样品大小适合放入扫描电镜的样品舱。为保证观察效果，利用砂纸对样品表面进行打磨，从粗砂纸逐渐过渡到细砂纸，使样品表面光滑平整，然后进行抛光处理，去除表面划痕和杂质。将制备好的样品放入真空镀膜机中，在样品表面蒸镀一层厚度约为10-20nm的金膜或铂膜，以增强样品的导电性，防止在电子束照射下产生电荷积累，影响图像质量。把镀膜后的样品放置在扫描电镜的样品台上，调整样品位置，使样品位于电子束的中心轴线上。设置扫描电镜的工作参数，加速电压一般选择10-20kV，根据样品的具体情况和观察需求，选择合适的放大倍数，从低倍数（如500倍）开始观察，了解样品的整体结构，再逐渐提高放大倍数（如5000倍、10000倍），观察样品的微观细节。调整扫描电镜的探测器参数，获取清晰的二次电子图像和背散射电子图像。在观察过程中，对不同区域的岩石微观结构进行拍照记录，包括孔隙、裂隙、矿物颗粒的分布和排列等特征。通过图像处理软件，对获取的微观结构图像进行分析，测量孔隙的大小、形状和连通性，统计矿物颗粒的粒径分布等信息，从而深入了解岩石的微观结构特征及其在吸水前后的变化情况。3.3.3吸水试验吸水试验采用自然吸水和有压吸水两种方式。在自然吸水试验中，先将岩石试件放入105-110℃的烘箱中烘至恒重，烘干时间一般为12-24小时，然后取出置于干燥器内冷却至室温，用精度为0.001g的电子天平称其质量m₀。将称量后的试件置于盛水容器内，先注水至试件高度的1/4处，以后每隔2小时分别注水至试件高度的1/2和3/4处，6小时后将水加至高出试件顶面20mm，以利于试件内空气逸出。试件全部被水淹没后再自由吸水48小时。每隔一定时间（如1小时）用电子天平称取试件质量m₁，记录吸水量（m₁-m₀）和时间。有压吸水试验则使用专门的压力装置。将烘干冷却后的岩石试件放入压力密封舱中，向舱内注水，使水面高出试件顶面20mm，然后通过加压泵逐渐增加压力至设定值（如0.5MPa、1MPa），保持压力恒定。每隔1小时打开压力舱，取出试件，用湿纱布迅速擦去表面水分，称取试件质量m₂，记录吸水量（m₂-m₀）和时间。重复此操作，直至试件的吸水量基本不再增加，认为试件达到饱和吸水状态。通过自然吸水和有压吸水试验，全面了解岩石在不同条件下的吸水特性。3.3.4单轴压缩试验利用岩石力学试验系统对不同吸水率的岩石进行单轴压缩试验。先将加工好的标准岩石试件（直径50mm，高100mm）安装在试验机的上下压板之间，调整试件位置，使其中心与上下压板的中心重合，确保加载过程中试件受力均匀。设置试验机的加载参数，采用位移控制加载方式，加载速率为0.05mm/min。在加载过程中，试验机以设定的速率缓慢施加轴向压力，同时利用试验机配备的传感器实时采集试件的轴向应力和轴向应变数据。当试件达到破坏时，记录破坏荷载和对应的应变值。对于不同吸水率的岩石试件，重复上述加载过程。通过分析不同吸水率下岩石的应力-应变曲线，获取岩石的单轴抗压强度、弹性模量等力学参数。例如，根据应力-应变曲线的初始直线段斜率计算弹性模量，单轴抗压强度则为试件破坏时的最大应力值。通过对比不同吸水率岩石的力学参数，研究岩石力学性能随吸水率的变化规律。四、试验结果与分析4.1岩石矿物成分与微观结构分析4.1.1黏土矿物成分分析通过X射线衍射试验，对荣华矿深部-650水平大巷围岩中的粗砂岩、细砂岩和泥质粉砂岩的矿物成分进行了精确分析，结果见表1。表1不同岩石的黏土矿物成分含量（%）岩石类型蒙脱石伊利石高岭石绿泥石粗砂岩2.55.33.12.0细砂岩3.86.74.22.5泥质粉砂岩10.512.08.55.0从表1可以看出，不同岩石中黏土矿物的种类和含量存在明显差异。泥质粉砂岩中黏土矿物含量最高，总量达到36.0%，其中蒙脱石含量为10.5%，伊利石含量为12.0%。蒙脱石具有极强的吸水性和膨胀性，其晶体结构中存在可交换的阳离子，当与水接触时，水分子会进入蒙脱石晶层之间，导致晶层间距增大，从而使蒙脱石发生膨胀。这种膨胀作用会对泥质粉砂岩的结构产生破坏，使其强度降低。伊利石则会影响岩石的力学性质，其存在使得岩石的脆性增加，在受力时更容易发生破裂。细砂岩中黏土矿物含量次之，总量为17.2%，其中蒙脱石含量为3.8%，伊利石含量为6.7%。虽然含量相对泥质粉砂岩较低，但在长期受到地下水作用时，这些黏土矿物仍会发生吸水膨胀等反应，对细砂岩的物理力学性质产生一定影响，导致其强度逐渐下降。粗砂岩中黏土矿物含量最少，总量仅为12.9%，其中蒙脱石含量为2.5%，伊利石含量为5.3%。尽管黏土矿物含量较少，但由于粗砂岩在大巷围岩中广泛分布，其吸水软化特性对巷道稳定性的影响也不容忽视。在高湿度或地下水丰富的区域，粗砂岩中的黏土矿物会逐渐吸收水分，虽然其膨胀程度相对较小，但长期积累可能导致岩石内部结构的细微变化，进而影响岩石的整体力学性能。黏土矿物的存在是影响岩石吸水软化的关键因素之一。黏土矿物含量越高，岩石的吸水性越强，吸水后发生软化的程度也越明显。这是因为黏土矿物的晶体结构和表面性质使其具有较大的比表面积和较强的亲水性，容易吸附水分子。当岩石中的黏土矿物吸收水分后，会发生膨胀、分散等变化，破坏岩石颗粒之间的胶结力，导致岩石的孔隙结构发生改变，力学强度降低。例如，泥质粉砂岩中较高含量的黏土矿物使其在吸水后，强度下降幅度明显大于粗砂岩和细砂岩。因此，在研究岩石吸水软化特性和巷道支护设计时，必须充分考虑黏土矿物成分及其含量的影响。4.1.2微观结构特征分析利用扫描电子显微镜（SEM）对岩石微观结构进行观察，获得了不同岩石在吸水前后的微观结构图像，通过对图像的分析，揭示了岩石微观结构特征及其对吸水的影响。在吸水前，粗砂岩呈现出典型的粗粒结构，颗粒粒径较大，一般在0.5-2mm之间。颗粒之间通过硅质胶结和钙质胶结紧密相连，孔隙主要以粒间孔隙为主，孔隙大小相对均匀，孔径多在10-50μm之间，且孔隙连通性较差。这种结构使得粗砂岩在一定程度上具有较好的力学强度和较低的吸水性。例如，从图1（a）中可以清晰地看到粗砂岩颗粒的紧密排列和相对规则的孔隙分布。细砂岩具有细粒结构，颗粒粒径在0.1-0.5mm之间，颗粒分布较为均匀。其胶结类型主要为泥质胶结和铁质胶结，孔隙类型除了粒间孔隙外，还存在部分微孔和介孔。微孔和介孔的存在增加了细砂岩的比表面积，使其具有一定的吸水性。吸水前，细砂岩的孔隙相对较小，孔径多在1-10μm之间，孔隙连通性一般。如图1（b）所示，细砂岩的微观结构较为致密，但仍能观察到一些微小的孔隙分布在颗粒之间。泥质粉砂岩具有粉砂质结构，颗粒粒径细小，多在0.01-0.1mm之间。其胶结方式主要为泥质胶结，这种胶结方式使得岩石结构较为松散。泥质粉砂岩中含有大量的黏土矿物，这些黏土矿物以片状或鳞片状形态存在，相互交织形成了复杂的孔隙网络。孔隙类型主要为微孔和介孔，孔径多在0.1-1μm之间，孔隙率较高，且孔隙连通性较好。由于其特殊的微观结构和高孔隙率，泥质粉砂岩在吸水前就具有较强的吸水性。从图1（c）中可以明显看出泥质粉砂岩中黏土矿物的片状结构和丰富的孔隙。（此处插入图1：（a）吸水前粗砂岩SEM图像；（b）吸水前细砂岩SEM图像；（c）吸水前泥质粉砂岩SEM图像）当岩石吸水后，微观结构发生了显著变化。粗砂岩吸水后，部分胶结物被水溶解，颗粒间的胶结力减弱，孔隙增大且连通性增强。一些原本孤立的孔隙开始相互连通，形成更大的孔隙通道，这使得粗砂岩的吸水性进一步增强。从图2（a）中可以看到，吸水后的粗砂岩颗粒间出现了明显的缝隙，孔隙变得更加不规则。细砂岩吸水后，泥质胶结物发生软化，颗粒之间的连接变得松散，孔隙结构也发生了改变。部分微孔和介孔扩大，孔隙连通性进一步提高，导致细砂岩的力学强度降低，吸水性增强。图2（b）显示，吸水后的细砂岩微观结构变得较为疏松，孔隙数量增多且尺寸增大。泥质粉砂岩吸水后，黏土矿物发生膨胀，使得岩石内部微观结构变得更加松散。孔隙分布更加不均匀，部分孔隙被膨胀的黏土矿物堵塞，而部分孔隙则进一步扩大。这种微观结构的变化导致泥质粉砂岩的力学性能急剧下降，同时其吸水性也达到了较高水平。如图2（c）所示，吸水后的泥质粉砂岩中黏土矿物明显膨胀，孔隙结构变得极为复杂。（此处插入图2：（a）吸水后粗砂岩SEM图像；（b）吸水后细砂岩SEM图像；（c）吸水后泥质粉砂岩SEM图像）岩石的微观结构特征对其吸水性能有着重要影响。孔隙的大小、形状、连通性以及矿物颗粒的排列方式和胶结情况等因素，共同决定了岩石的吸水能力和吸水后微观结构的变化。孔隙率高、连通性好的岩石，如泥质粉砂岩，更容易吸水；而胶结力强、孔隙结构相对稳定的岩石，如粗砂岩，吸水能力相对较弱。在荣华矿深部-650水平大巷的工程实践中，了解岩石微观结构与吸水性能的关系，对于预测巷道围岩的变形和破坏，制定合理的支护措施具有重要意义。4.2岩石吸水特性分析4.2.1吸水特征曲线通过自然吸水和有压吸水试验，获取了粗砂岩、细砂岩和泥质粉砂岩在不同时间下的吸水量数据，绘制出如图3所示的吸水特征曲线。（此处插入图3：粗砂岩、细砂岩和泥质粉砂岩吸水特征曲线）从图3可以看出，三种岩石的吸水过程均呈现出相似的变化趋势，但在吸水量和吸水速度上存在明显差异。在吸水初期，三种岩石的吸水量均迅速增加，吸水速度较快。随着时间的推移，吸水量的增长速度逐渐减缓，最终趋于稳定，达到饱和吸水状态。对于粗砂岩，在吸水初期（0-6小时），其吸水量增长迅速，6小时时吸水量达到0.5g左右，此时吸水速度约为0.083g/h。之后，吸水量增长速度逐渐变慢，在24小时左右，吸水量基本不再增加，达到饱和状态，饱和吸水量约为0.6g。细砂岩的吸水速度相对粗砂岩略快，在吸水初期（0-4小时），吸水量快速上升，4小时时吸水量达到0.6g左右，吸水速度约为0.15g/h。随后，吸水速度逐渐降低，在18小时左右达到饱和，饱和吸水量约为0.8g。泥质粉砂岩的吸水性最强，吸水速度也最快。在吸水初期（0-2小时），吸水量急剧增加，2小时时吸水量已达到0.8g左右，吸水速度高达0.4g/h。在12小时左右，泥质粉砂岩就基本达到饱和状态，饱和吸水量约为1.2g。根据吸水特征曲线的变化规律，可以将岩石的吸水过程划分为三个阶段：快速吸水阶段、减速吸水阶段和饱和吸水阶段。在快速吸水阶段，岩石内部的孔隙和裂隙较多，水能够迅速进入岩石内部，导致吸水量快速增加；随着水的不断进入，岩石内部的孔隙逐渐被填充，水的流动阻力增大，吸水速度逐渐降低，进入减速吸水阶段；当岩石内部的孔隙几乎被水完全填充时，吸水量基本不再变化，岩石达到饱和吸水状态。4.2.2影响吸水的因素矿物成分：岩石的矿物成分对其吸水特性有着重要影响。黏土矿物含量高的岩石，如泥质粉砂岩，由于黏土矿物具有较强的吸水性和膨胀性，使得岩石的吸水性显著增强。泥质粉砂岩中黏土矿物含量可达40%-50%，其中蒙脱石等黏土矿物的吸水膨胀作用，使得岩石能够吸收大量水分。而粗砂岩和细砂岩中黏土矿物含量相对较低，其吸水性也相对较弱。此外，长石等矿物在一定程度上也会影响岩石的吸水性，长石的水解等化学反应可能会改变岩石的结构，为水分的进入提供通道。微观结构：岩石的微观结构，包括孔隙大小、形状、连通性以及矿物颗粒的排列方式和胶结情况等，是影响吸水的关键因素。孔隙率高、连通性好的岩石，如泥质粉砂岩，其内部存在大量的微孔和介孔，且孔隙之间相互连通，为水的运移提供了良好的通道，使得岩石更容易吸水。泥质粉砂岩的孔隙率可达20%-30%，其孔隙结构有利于水分的快速进入和储存。而粗砂岩虽然孔隙较大，但连通性较差，其吸水能力相对较弱。细砂岩的孔隙大小和连通性介于粗砂岩和泥质粉砂岩之间，因此其吸水性也处于两者之间。孔隙率：孔隙率直接决定了岩石能够容纳水分的空间大小。一般来说，孔隙率越高，岩石的吸水量越大。泥质粉砂岩的高孔隙率使其具有较大的吸水空间，能够吸收较多的水分。通过压汞仪（MIP）测试可知，泥质粉砂岩的孔隙率明显高于粗砂岩和细砂岩，这也解释了为什么泥质粉砂岩的饱和吸水量最大。此外，孔隙的大小分布也会影响吸水速度，较小的孔隙在吸水初期可能会限制水的进入速度，但在后期可能会对水分的储存起到重要作用。综合来看，矿物成分、微观结构和孔隙率是影响荣华矿深部-650水平大巷围岩吸水特性的主要因素。其中，黏土矿物含量和孔隙结构的影响最为显著，它们相互作用，共同决定了岩石的吸水能力和吸水速度。在巷道支护设计和施工中，充分考虑这些因素，对于采取有效的防水、排水措施，减少岩石吸水，提高巷道围岩的稳定性具有重要意义。4.3岩石吸水软化特性分析4.3.1单轴抗压强度与吸水率关系通过对不同吸水率的粗砂岩、细砂岩和泥质粉砂岩进行单轴压缩试验，获得了岩石单轴抗压强度随吸水率的变化数据，绘制出如图4所示的关系曲线。（此处插入图4：不同岩石单轴抗压强度与吸水率关系曲线）从图4可以看出，随着吸水率的增加，三种岩石的单轴抗压强度均呈现出下降趋势。粗砂岩的单轴抗压强度从初始的55MPa左右，下降到吸水率为3%时的40MPa左右，下降幅度约为27%。细砂岩的单轴抗压强度从初始的45MPa左右，下降到吸水率为4%时的30MPa左右，下降幅度约为33%。泥质粉砂岩的单轴抗压强度下降最为明显，从初始的30MPa左右，下降到吸水率为6%时的15MPa左右，下降幅度高达50%。对单轴抗压强度与吸水率的关系进行函数拟合，发现粗砂岩和细砂岩的单轴抗压强度与吸水率之间服从幂函数分布，拟合函数分别为：粗砂岩：σ_{c1}=58.62W^{-0.18}（R^{2}=0.95）细砂岩：σ_{c2}=47.35W^{-0.22}（R^{2}=0.93）其中，σ_{c1}、σ_{c2}分别为粗砂岩和细砂岩的单轴抗压强度（MPa），W为吸水率（%），R^{2}为拟合优度。泥质粉砂岩的单轴抗压强度与吸水率之间服从线性函数分布，拟合函数为：σ_{c3}=-2.5W+31.5（R^{2}=0.97），其中σ_{c3}为泥质粉砂岩的单轴抗压强度（MPa）。岩石单轴抗压强度随吸水率下降的原因主要是由于吸水后岩石内部结构发生变化。一方面，水进入岩石孔隙和裂隙中，削弱了岩石颗粒之间的胶结力，使得岩石在受力时更容易发生颗粒间的相对滑动和错动，从而降低了岩石的强度。另一方面，对于含有黏土矿物的岩石，如泥质粉砂岩，黏土矿物吸水膨胀，进一步破坏了岩石的内部结构，导致强度大幅下降。这种单轴抗压强度与吸水率的关系，对于评估荣华矿深部-650水平大巷围岩在地下水作用下的强度变化具有重要意义，为巷道支护设计提供了关键的强度参数依据。4.3.2弹性模量与吸水率关系分析不同吸水率下岩石的应力-应变曲线，计算得到岩石的弹性模量，绘制出弹性模量与吸水率的关系曲线，如图5所示。（此处插入图5：不同岩石弹性模量与吸水率关系曲线）从图5可以明显看出，随着吸水率的增大，粗砂岩、细砂岩和泥质粉砂岩的弹性模量均逐渐减小。粗砂岩的弹性模量从初始的12GPa左右，下降到吸水率为3%时的8GPa左右，下降幅度约为33%。细砂岩的弹性模量从初始的10GPa左右，下降到吸水率为4%时的6GPa左右，下降幅度约为40%。泥质粉砂岩的弹性模量下降幅度最大，从初始的8GPa左右，下降到吸水率为6%时的3GPa左右，下降幅度达到62.5%。对弹性模量与吸水率的关系进行函数拟合，结果表明，三种岩石的弹性模量与吸水率均服从幂函数分布。粗砂岩的拟合函数为：E_{1}=13.5W^{-0.25}（R^{2}=0.94）；细砂岩的拟合函数为：E_{2}=11.2W^{-0.28}（R^{2}=0.92）；泥质粉砂岩的拟合函数为：E_{3}=9.5W^{-0.35}（R^{2}=0.96）。其中，E_{1}、E_{2}、E_{3}分别为粗砂岩、细砂岩和泥质粉砂岩的弹性模量（GPa），W为吸水率（%）。弹性模量反映了岩石抵抗弹性变形的能力，其随吸水率的减小，意味着岩石在吸水后更容易发生弹性变形。这是因为水的侵入使得岩石内部的孔隙和裂隙扩张，颗粒间的连接变弱，导致岩石的整体刚度降低。在巷道开挖过程中，围岩的弹性模量减小会使巷道周边的变形增大，增加了巷道支护的难度。例如，当巷道围岩为泥质粉砂岩时，由于其弹性模量随吸水率下降幅度较大，在地下水作用下，巷道更容易出现较大的变形，如顶板下沉、两帮收敛等，从而对巷道的稳定性产生严重影响。因此，了解岩石弹性模量与吸水率的关系，对于预测巷道围岩的变形、合理选择支护方式和确定支护参数具有重要的指导意义。4.3.3软化系数与吸水率关系软化系数是衡量岩石耐水性的重要指标，通过计算不同吸水率下岩石的软化系数，分析其与吸水率的关系，结果见表2。表2不同岩石软化系数与吸水率关系岩石类型吸水率（%）软化系数粗砂岩1.00.92粗砂岩2.00.85粗砂岩3.00.78细砂岩1.50.88细砂岩2.50.82细砂岩3.50.75泥质粉砂岩2.00.75泥质粉砂岩3.00.68泥质粉砂岩4.00.60根据表2数据，绘制出软化系数与吸水率的关系曲线，如图6所示。（此处插入图6：不同岩石软化系数与吸水率关系曲线）从图6可以看出，随着吸水率的增加，三种岩石的软化系数均逐渐减小，表明岩石的软化程度逐渐加剧。粗砂岩的软化系数从吸水率为1.0%时的0.92下降到吸水率为3.0%时的0.78，下降幅度为15.2%。细砂岩的软化系数从吸水率为1.5%时的0.88下降到吸水率为3.5%时的0.75，下降幅度为14.8%。泥质粉砂岩的软化系数下降最为显著，从吸水率为2.0%时的0.75下降到吸水率为4.0%时的0.60，下降幅度达到20%。对粗砂岩和细砂岩，其软化系数与吸水率之间服从幂函数分布，拟合函数分别为：粗砂岩：K_{r1}=0.98W^{-0.07}（R^{2}=0.93）细砂岩：K_{r2}=0.95W^{-0.08}（R^{2}=0.91）其中，K_{r1}、K_{r2}分别为粗砂岩和细砂岩的软化系数，W为吸水率（%）。泥质粉砂岩的软化系数与吸水率服从线性函数分布，拟合函数为：K_{r3}=-0.075W+0.89（R^{2}=0.96），其中K_{r3}为泥质粉砂岩的软化系数。岩石软化系数的减小主要是由于吸水后岩石内部的矿物成分和微观结构发生变化。对于含有黏土矿物的岩石，黏土矿物吸水膨胀，破坏了岩石的内部结构，降低了岩石的强度，从而导致软化系数减小。软化系数越小，岩石在水的作用下越容易软化，其力学性能越差。在荣华矿深部-650水平大巷的工程实践中，准确评估岩石的软化系数与吸水率的关系，对于判断巷道围岩在地下水长期作用下的稳定性至关重要。当岩石的软化系数较低时，需要采取更加强化的支护措施，以防止巷道围岩因软化而发生破坏。4.4岩石损伤分析4.4.1损伤变量定义基于岩石损伤力学理论，本研究选用弹性模量作为损伤变量，来定量描述岩石在吸水软化过程中的损伤程度。岩石的弹性模量是反映其抵抗弹性变形能力的重要参数，在岩石吸水软化过程中，随着内部结构的劣化，弹性模量会发生显著变化。损伤变量D的计算方法如下：D=1-\frac{E}{E_0}其中，E为岩石吸水后某一状态下的弹性模量，E_0为岩石初始状态（未吸水时）的弹性模量。该损伤变量的物理意义在于，它表示岩石由于吸水等因素导致的弹性模量降低程度，进而反映岩石内部结构的损伤程度。当D=0时，意味着岩石没有发生损伤，其弹性模量保持初始值，内部结构完整；随着D值逐渐增大，表明岩石的弹性模量不断减小，内部结构损伤加剧，岩石抵抗变形的能力逐渐减弱。例如，当D=0.2时，说明岩石的弹性模量已经降低到初始值的80%，岩石内部的微观结构可能已经发生了较为明显的变化，如孔隙扩张、颗粒间胶结力减弱等，从而影响岩石的力学性能。4.4.2损伤演化方程建立根据不同吸水率条件下岩石的弹性模量试验数据，建立岩石在不同吸水率条件下的损伤演化方程。以粗砂岩为例，通过对试验数据进行分析和拟合，得到其损伤演化方程为：D=0.08W^{0.32}其中，D为损伤变量，W为吸水率（%）。对于细砂岩，损伤演化方程为：D=0.12W^{0.35}泥质粉砂岩的损伤演化方程为：D=0.2W^{0.4}从这些损伤演化方程可以看出，随着吸水率的增加，损伤变量逐渐增大，表明岩石的损伤程度不断加剧。不同岩石的损伤演化方程中，系数和指数存在差异，这反映了不同岩石在吸水软化过程中损伤发展规律的不同。泥质粉砂岩的系数相对较大，说明其在相同吸水率条件下，损伤程度相对粗砂岩和细砂岩更为严重，这与泥质粉砂岩中黏土矿物含量高、吸水性强、吸水后结构破坏明显的特性相符。通过对损伤演化方程的分析，可以进一步了解岩石在吸水软化过程中的损伤发展趋势。在低吸水率阶段，损伤变量增长相对缓慢，这是因为此时岩石内部的孔隙和裂隙尚未被水充分填充，岩石结构的破坏程度较小。随着吸水率的不断增加，进入中高吸水率阶段，损伤变量增长速度加快，岩石内部结构受到的破坏加剧，孔隙扩张、颗粒间胶结力丧失等现象更加明显，导致岩石的损伤程度迅速增大。在荣华矿深部-650水平大巷的工程实际中，利用这些损伤演化方程，可以预测不同围岩条件下岩石的损伤发展情况，为巷道支护设计和施工提供重要依据。在地下水丰富的区域，根据岩石的吸水率和损伤演化方程，提前采取加强支护措施，以应对岩石损伤加剧带来的巷道稳定性问题。五、数值模拟分析5.1数值模拟软件介绍本研究采用RFPA2D岩石破坏过程分析系统对荣华矿深部-650水平大巷围岩在不同工况下的力学行为进行数值模拟。RFPA2D是一款基于弹性损伤模型的专业数值模拟工具，在岩石力学研究领域应用广泛。该软件的核心功能在于能够实现对岩石破坏过程的全程模拟，从岩石内部微裂纹的萌生、扩展，到变形局部化、成核，最终导致失稳破坏的各个阶段，都能进行细致的模拟分析。其原理基于以下几个关键要点：首先，将岩石介质模型离散化成由细观基元组成的数值模型，假设岩石介质在细观上是各向同性的弹-脆性介质。这些细观基元是构成介质的基本单元，其力学性质服从某种统计分布规律，通常采用Weibull分布来描述基元强度的非均匀性。通过这种方式，建立起细观与宏观介质力学性能的联系，能够更真实地反映岩石材料的非均匀特性，因为实际岩石内部存在着各种缺陷和不均匀性，如微裂纹、孔隙、矿物颗粒分布不均等，而这种统计分布假设可以有效考虑这些因素对岩石力学行为的影响。在应力分析方面，RFPA2D利用线弹性有限元方法作为应变计算器，按照弹性力学中的基元线弹性应力、应变求解方法，分析模型的应力、应变状态。该方法能够准确计算模型在各种载荷条件下的应力和应变分布，为后续分析岩石的变形和破坏提供基础数据。例如，在模拟巷道开挖过程中，通过有限元计算可以得到巷道周围岩石的应力集中区域和应变分布情况，从而预测可能出现的破坏位置和形式。为了判断基元的破坏，RFPA2D引入了适当的基元破坏准则和损伤规律。基元的相变临界点采用修正的Coulomb准则，该准则综合考虑了岩石的抗剪强度、正应力和剪应力等因素，能够较为准确地判断岩石在复杂应力状态下的破坏情况。同时，认为基元的力学性质随演化的发展是不可逆的，即一旦基元发生破坏，其力学性能将发生永久性改变，不再恢复到初始状态。并且，基元相变前后均被视为线弹性体，这一假设简化了计算过程，同时在一定程度上能够反映岩石在破坏前后的力学行为特征。此外，岩石介质中的裂纹扩展被视为一个准静态过程，忽略因快速扩展引起的惯性力的影响。在大多数工程实际问题中，这种假设是合理的，因为岩石的破坏过程相对缓慢，惯性力的影响较小，可以忽略不计。与其他岩石力学数值模拟软件相比，RFPA2D具有独特的优势。它能够充分考虑岩石材料的非均匀性和细观结构特征，通过对基元力学性质的统计分布假设和破坏准则的合理设定，更真实地模拟岩石的破坏过程。在模拟含有大量微裂纹和孔隙的岩石时，RFPA2D可以准确地反映这些缺陷对岩石力学性能的影响，而一些传统的数值模拟软件可能无法很好地处理这些复杂的细观结构。其模拟结果与实际岩石试验结果具有较好的一致性，经过大量的实验验证，在模拟岩石单轴压缩、三轴压缩等力学试验时，RFPA2D能够准确地再现岩石的破坏过程和力学响应，为岩石力学研究和工程应用提供了可靠的参考依据。五、数值模拟分析5.1数值模拟软件介绍本研究采用RFPA2D岩石破坏过程分析系统对荣华矿深部-650水平大巷围岩在不同工况下的力学行为进行数值模拟。RFPA2D是一款基于弹性损伤模型的专业数值模拟工具，在岩石力学研究领域应用广泛。该软件的核心功能在于能够实现对岩石破坏过程的全程模拟，从岩石内部微裂纹的萌生、扩展，到变形局部化、成核，最终导致失稳破坏的各个阶段，都能进行细致的模拟分析。其原理基于以下几个关键要点：首先，将岩石介质模型离散化成由细观基元组成的数值模型，假设岩石介质在细观上是各向同性的弹-脆性介质。这些细观基元是构成介质的基本单元，其力学性质服从某种统计分布规律，通常采用Weibull分布来描述基元强度的非均匀性。通过这种方式，建立起细观与宏观介质力学性能的联系，能够更真实地反映岩石材料的非均匀特性，因为实际岩石内部存在着各种缺陷和不均匀性，如微裂纹、孔隙、矿物颗粒分布不均等，而这种统计分布假设可以有效考虑这些因素对岩石力学行为的影响。在应力分析方面，RFPA2D利用线弹性有限元方法作为应变计算器，按照弹性力学中的基元线弹性应力、应变求解方法，分析模型的应力、应变状态。该方法能够准确计算模型在各种载荷条件下的应力和应变分布，为后续分析岩石的变形和破坏提供基础数据。例如，在模拟巷道开挖过程中，通过有限元计算可以得到巷道周围岩石的应力集中区域和应变分布情况，从而预测可能出现的破坏位置和形式。为了判断基元的破坏，RFPA2D引入了适当的基元破坏准则和损伤规律。基元的相变临界点采用修正的Coulomb准则，该准则综合考虑了岩石的抗剪强度、正应力和剪应力等因素，能够较为准确地判断岩石在复杂应力状态下的破坏情况。同时，认为基元的力学性质随演化的发展是不可逆的，即一旦基元发生破坏，其力学性能将发生永久性改变，不再恢复到初始状态。并且，基元相变前后均被视为线弹性体，这一假设简化了计算过程，同时在一定程度上能够反映岩石在破坏前后的力学行为特征。此外，岩石介质中的裂纹扩展被视为一个准静态过程，忽略因快速扩展引起的惯性力的影响。在大多数工程实际问题中，这种假设是合理的，因为岩石的破坏过程相对缓慢，惯性力的影响较小，可以忽略不计。与其他岩石力学数值模拟软件相比，RFPA2D具有独特的优势。它能够充分考虑岩石材料的非均匀性和细观结构特征，通过对基元力学性质的统计分布假设和破坏准则的合理设定，更真实地模拟岩石的破坏过程。在模拟含有大量微裂纹和孔隙的岩石时，RFPA2D可以准确地反映这些缺陷对岩石力学性能的影响，而一些传统的数值模拟软件可能无法很好地处理这些复杂的细观结构。其模拟结果与实际岩石试验结果具有较好的一致性，经过大量的实验验证，在模拟岩石单轴压缩、三轴压缩等力学试验时，RFPA2D能够准确地再现岩石的破坏过程和力学响应，为岩石力学研究和工程应用提供了可靠的参考依据。5.2模型建立与参数设置5.2.1模型构建运用RFPA2D软件，建立岩石单轴压缩数值模型。模型尺寸设定为长100mm、宽50mm，该尺寸参考了实际岩石试件的常见尺寸，能够较好地模拟岩石在单轴压缩条件下的力学行为。在网格划分方面，采用等面积四节点的四边形单元对模型进行剖分，将模型划分为200×100个基元。较小的单元尺寸能够更精确地反映岩石材料的非均匀性，虽然会增加一定的计算量，但随着计算机技术的发展，这一问题已得到有效缓解。通过合理的网格划分，确保模型能够准确模拟岩石内部微裂纹的萌生、扩展以及最终的破坏过程。在边界条件设置上，模型底部采用固定约束，限制其在x和y方向的位移，模拟岩石试件底部在试验中被固定的状态；顶部则施加位移加载，以模拟实际单轴压缩试验中压力机对岩石试件施加的轴向压力。加载过程通过位移加载方式实现，加载速率设置为0.001mm/step，这一加载速率经过多次调试确定，既能保证模拟过程的稳定性，又能较为真实地反映实际加载过程中岩石的力学响应。同时，为了模拟岩石在自然状态下的受力情况，不考虑模型的自重。通过这样的边界条件设置，使数值模型能够尽可能地接近实际的岩石单轴压缩试验条件。5.2.2参数设置根据前文的试验结果，确定数值模型中岩石的各项参数。对于力学参数，粗砂岩的弹性模量设置为12GPa，单轴抗压强度为55MPa，泊松比为0.25；细砂岩的弹性模量为10GPa，单轴抗压强度为45MPa，泊松比为0.28；泥质粉砂岩的弹性模量为8GPa，单轴抗压强度为30MPa，泊松比为0.3。这些力学参数是通过对不同岩石进行大量的力学试验得到的平均值，能够较好地反映岩石的基本力学特性。在吸水参数方面，根据吸水试验结果，粗砂岩的饱和吸水率设定为3%，细砂岩为4%，泥质粉砂岩为6%。同时，考虑到岩石吸水过程的复杂性，引入吸水系数来描述岩石的吸水速度。粗砂岩的吸水系数设置为0.05g/h，细砂岩为0.1g/h，泥质粉砂岩为0.2g/h。这些吸水参数的设置，能够较为准确地模拟岩石在不同环境下的吸水过程，以及吸水对岩石力学性能的影响。损伤参数则根据岩石损伤分析的结果进行确定。损伤变量通过弹性模量来定义，损伤演化方程根据不同岩石在不同吸水率条件下的试验数据拟合得到。例如，粗砂岩的损伤演化方程为D=0.08W^{0.32}，细砂岩为D=0.12W^{0.35}，泥质粉砂岩为D=0.2W^{0.4}。在数值模型中，当岩石单元的应力或应变状态满足损伤阈值准则时，单元开始损伤，损伤后的力学性质根据损伤变量进行调整。通过合理设置损伤参数，能够在数值模拟中准确反映岩石在吸水软化过程中的损伤发展情况，为分析岩石的破坏机制提供依据。5.3模拟结果分析5.3.1载荷-位移曲线分析通过数值模拟得到的粗砂岩、细砂岩和泥质粉砂岩在单轴压缩条件下的载荷-位移曲线，与室内试验结果进行对比，结果如图7所示。（此处插入图7：数值模拟与试验的载荷-位移曲线对比）从图7中可以看出，数值模拟得到的载荷-位移曲线与试验曲线具有较好的一致性，验证了数值模型的准确性。在加载初期，曲线呈现出近似线性的变化趋势，这表明岩石处于弹性变形阶段，应力与应变呈线性关系，符合胡克定律。此时，岩石内部的微裂纹尚未大量萌生和扩展，岩石主要表现为弹性变形，能够较好地抵抗外力作用。随着载荷的逐渐增加，曲线开始偏离线性，进入非线性变形阶段。在这个阶段，岩石内部的微裂纹开始萌生和扩展，岩石的内部结构逐渐发生变化，导致其力学性能下降，应力与应变不再保持线性关系。对于粗砂岩，当载荷达到40-45MPa时，曲线开始明显偏离线性，这是因为粗砂岩内部的颗粒间胶结力在逐渐增大的载荷作用下开始受到破坏，微裂纹开始在颗粒间萌生和扩展。细砂岩在载荷达到30-35MPa时进入非线性阶段，这是由于细砂岩的颗粒相对较小，胶结强度相对较低，更容易受到载荷的影响而产生微裂纹。泥质粉砂岩由于其黏土矿物含量高，结构相对松散，在载荷达到20-25MPa时就进入了非线性阶段，黏土矿物的吸水膨胀进一步加速了微裂纹的萌生和扩展。当载荷继续增加到一定程度时，曲线达到峰值，此时岩石达到其单轴抗压强度，随后载荷迅速下降，岩石发生破坏。粗砂岩的峰值载荷约为55MPa，细砂岩约为45MPa，泥质粉砂岩约为30MPa，这与试验结果基本相符。岩石达到峰值载荷后，内部的微裂纹迅速扩展并相互贯通，形成宏观的破裂面，导致岩石失去承载能力，载荷急剧下降。在数值模拟中，通过观察模型中基元的破坏情况，可以清晰地看到微裂纹的扩展和贯通过程，与实际岩石破坏过程中的现象相吻合。载荷-位移曲线的变化特征反映了岩石在单轴压缩过程中的力学行为和内部结构变化。在弹性阶段，岩石的变形主要是弹性变形，内部结构基本保持完整；进入非线性阶段后，微裂纹的萌生和扩展逐渐改变了岩石的内部结构，使其力学性能下降；当达到峰值载荷后，岩石内部结构彻底破坏，失去承载能力。这些特征对于深入理解岩石的破坏机制以及评估巷道围岩的稳定性具有重要意义。在荣华矿深部-650水平大巷的工程实践中，通过分析岩石的载荷-位移曲线，可以预测巷道围岩在不同受力条件下的变形和破坏情况，为制定合理的支护措施提供依据。5.3.2声发射能分析在数值模拟中，通过设置声发射监测功能，得到了不同岩石在单轴压缩过程中的声发射能分布和变化情况，如图8所示。（此处插入图8：不同岩石单轴压缩过程声发射能变化曲线）从图8可以看出，在加载初期，声发射能较低，这是因为此时岩石内部的微裂纹尚未大量产生，只有少量的原生微裂纹在载荷作用下发生闭合或微小扩展，释放出的能量较少。随着载荷的增加，声发射能逐渐增大，这表明岩石内部的微裂纹开始大量萌生和扩展，裂纹的扩展和相互作用释放出更多的能量。在粗砂岩中，当载荷达到30-40MPa时，声发射能开始明显增加，说明此时粗砂岩内部的微裂纹开始大量产生，颗粒间的胶结力逐渐被破坏。细砂岩在载荷达到20-30MPa时，声发射能显著上升，这是由于细砂岩的结构相对较细，在较低的载荷下就容易产生大量微裂纹。泥质粉砂岩在载荷达到10-20MPa时，声发射能就迅速增大，这是因为泥质粉砂岩中黏土矿物含量高，结构松散，对载荷的抵抗能力较弱，微裂纹更容易产生和扩展。当岩石接近破坏时，声发射能急剧增大，达到峰值。这是因为在岩石破坏前，内部的微裂纹迅速扩展并相互贯通，形成宏观的破裂面，这个过程伴随着大量能量的释放。粗砂岩在接近破坏时，声发射能峰值达到较高水平，表明其在破坏过程中释放出大量能量，这是由于粗砂岩的颗粒较大，内部结构相对较为致密，在破坏时需要克服较大的阻力，从而释放出更多能量。细砂岩和泥质粉砂岩的声发射能峰值相对较低，但在破坏时声发射能的增长速度也很快，说明它们在破坏时同样伴随着内部结构的剧烈变化和能量的快速释放。通过声发射能的变化，可以分析岩石在不同吸水条件下的损伤演化过程和破坏机制。在吸水条件下，岩石的力学性能会发生劣化，其损伤演化过程也会受到影响。当岩石吸水后，内部的孔隙和裂隙被水填充，水的润滑作用会降低岩石颗粒间的摩擦力，使得微裂纹更容易萌生和扩展。对于含有黏土矿物的岩石，如泥质粉砂岩，黏土矿物吸水膨胀会进一步破坏岩石的内部结构，加速微裂纹的扩展，导致声发射能在较低的载荷下就开始显著增加，岩石更容易发生破坏。在荣华矿深部-650水平大巷的实际工程中，声发射监测可以作为一种有效的手段来监测巷道围岩的损伤演化情况。通过布置在巷道周围的声发射传感器，实时监测岩石的声发射信号，根据声发射能的变化判断围岩的损伤程度和破坏趋势，及时采取相应的支护措施，保障巷道的安全稳定。5.3.3应力与弹性模量分布分析通过数值模拟，得到了岩石在单轴压缩过程中的应力和弹性模量分布图，如图9和图10所示。（此处插入图9：不同岩石单轴压缩过程应力分布图）（此处插入图10：不同岩石单轴压缩过程弹性模量分布图）从图9的应力分布图可以看出，在加载初期，岩石内部的应力分布相对均匀，随着载荷的增加，应力逐渐集中在某些区域。在粗砂岩中，应力集中主要发生在颗粒间的接触部位，这是因为粗砂岩的颗粒较大，颗粒间的接触点是受力的关键部位，在载荷作用下，这些接触点的应力首先达到岩石的强度极限，从而导致微裂纹的萌生。随着载荷进一步增加，应力集中区域逐渐扩大，微裂纹也随之扩展。细砂岩的应力集中区域相对较为分散，这是由于细砂岩的颗粒较小且分布相对均匀，应力在颗粒间的传递相对较为均匀，但在一些薄弱部位，如微孔隙和微裂纹周围，仍然会出现应力集中现象。泥质粉砂岩由于其结构松散，内部存在大量的孔隙和微裂纹，应力集中现象更为明显，且分布范围更广。在加载过程中，泥质粉砂岩的应力集中区域迅速扩大，导致其内部结构快速破坏。从图10的弹性模量分布图可以看出，在加载初期，岩石的弹性模量分布相对均匀，随着载荷的增加，弹性模量逐渐降低，且降低的区域与应力集中区域相对应。这是因为在应力集中区域，岩石内部的微裂纹萌生和扩展，导致岩石的内部结构损伤，从而使得弹性模量降低。在粗砂岩中，当载荷达到一定程度时，应力集中区域的弹性模量明显降低，这表明该区域的岩石已经发生了损伤，抵抗变形的能力减弱。细砂岩和泥质粉砂岩在加载过程中，弹性模量的降低更为明显，尤其是泥质粉砂岩，由于其结构的特殊性，在较低的载荷下弹性模量就开始大幅降低。这是因为泥质粉砂岩中的黏土矿物吸水膨胀，进一步破坏了岩石的内部结构，使得弹性模量迅速下降。在吸水过程中，岩石的应力和弹性模量分布也会发生变化。由于水的作用，岩石内部的孔隙和裂隙被水填充，导致岩石的有效应力减小，弹性模量降低。对于含有黏土矿物的岩石，黏土矿物吸水膨胀会改变岩石的内部结构，使得应力分布更加不均匀，弹性模量进一步降低。在荣华矿深部-650水平大巷的工程实践中，了解岩石在吸水过程中应力和弹性模量的变化规律，对于评估巷道围岩的稳定性至关重要。根据这些变化规律，可以合理设计巷道的支护参数，提高支护结构的承载能力，确保巷道的安全稳定。六、结论与展望6.1研究成果总结岩石吸水特性：通过自然吸水和有压吸水试验，获取了粗砂岩、细砂岩和泥质粉砂岩的吸水特征曲线。结果表明，三种岩石的吸水过程均呈现出相似趋势，初期吸水量迅速增加，随后增长速度逐渐减缓，最终趋于饱和。泥质粉砂岩的吸水性最强，饱和吸水量约为1.2g，吸水速度也最快；粗砂岩的吸水性相对较弱，饱和吸水量约为0.6g。岩石的矿物成分、微观结构和孔隙率是影响吸水特性的主要因素，黏土矿物含量高、孔隙率大且连通性好的岩石，吸水性更强。岩石软化特性：不同吸水率条件下的单轴压缩试验结果显示，随着吸水率的增加，粗砂岩、细砂岩和泥质粉砂岩的单轴抗压强度、弹性模量均呈现下降趋势，软化系数逐渐减小。粗砂岩和细砂岩的单轴抗压强度与吸水率服从幂函数分布，泥质粉砂岩服从线性函数分布；三种岩石的弹性模量与吸水率均服从幂函数分布；粗砂岩和细砂岩的软化系数与吸水率服从幂函数分布，泥质粉砂岩服从线性函数分布。岩石损伤演化规律：基于岩石损伤力学理论，以弹性模量作为损伤变量，建立了不同岩石在不同吸水率条件下的损伤演化方程。粗砂岩的损伤演化方程为D=0.08W^{0.32}，细砂岩为D=0.12W^{0.35}，泥质粉砂岩为D=0.2W^{0.4}。随着吸水率的增加，损伤变量逐渐增大，岩石损伤程度不断加剧，泥质粉砂岩的损伤程度相对更为严重。数值模拟结果：利用RFPA2D软件建立岩石单轴压缩数值模型，模拟结果与室内试验结果具有较好的一