目前，绿色矿山建设受到高度重视，且随着矿山开采逐渐转向深部，充填采矿法得到了广泛应用（曹国华等，2024；张长青等，2024）。在充填采矿过程中，充入采场的充填体与围岩组合在一起，形成充填体—围岩复合结构。充填体与围岩共同承压、相互作用和协调变形（Zhao et al，2022；Long et al，2024），导致单一的充填体或岩体力学特征与充—岩组合体的力学特征存在显著差异（Wang et al，2023）。探讨充—岩组合结构的力学与破坏特征及其影响因素对研究采空区的安全和稳定性具有重要意义。
针对充—岩组合体的力学特征，研究人员从充—岩界面特征（卢宏建等，2022a；徐文彬等，2023）、充—岩组合结构（冯国瑞等，2024；Yu et al，2024）和充填体配比（杨孝亮等，2022；张雯等，2024）等方面，探讨了其对充—岩组合体力学与破坏特征的影响规律，丰富了充—岩组合体相互作用的力学理论。应力加载形式也是影响充—岩组合体力学特征的因素之一。李胜等（2022）在单轴压缩条件下，得出组合体的初始损伤、延性破坏特征、破坏后的残余强度与充填体厚度呈正相关关系。龙大愚等（2024）在疲劳荷载条件下，得出灰砂比下降导致组合体体积膨胀增加，割线模量先增加后减小，提出了累积损伤演化模型。为进一步揭示充—岩组合体的损伤破坏规律，数值仿真技术被用于充—岩组合体破坏机理的研究。汪杰等（2020）基于PFC颗粒流软件，得到微裂纹先从充填体内部萌生，之后扩展到围岩中，导致组合体发生整体性破坏。卢宏建等（2022b）利用RFPA3D数值模拟揭示了组合体的宏观剪切—张拉混合裂纹破坏类型。
然而，充填体力学参数变化会导致充填体与岩石力学参数比值的改变，表明当前人们对充—岩组合体力学特征的产生及其影响机理的认识尚不充分，进而说明当前对充—岩组合体破坏过程中裂纹萌生与扩展规律的研究也存在不足。为此，本文通过控制充填体的料浆浓度来改变充填体与岩石的弹性模量、峰值应变和峰值强度的比值，制备了不同充—岩参数比值的试样，开展单轴压缩力学试验，并采用声发射和数字图像（DIC）技术对单轴压缩过程中组合体的破坏过程进行监测。分析弹性模量、峰值应变和峰值强度比值的变化对充—岩组合体相应力学特征变化规律的影响，探究不同比值下充—岩组合体的裂纹扩展规律，进一步完善充—岩组合体的相互作用机理，为矿山充填配比优化提供参考，从而保障地下采场的稳定。
1 组合体单轴压缩试验方案
1.1 试验材料与试样制备
充—岩组合试样制备：先将岩石试样放入内径为50 mm×50 mm×100 mm的长方体模具右侧，再将混合均匀的充填料浆（灰砂比为1∶4，料浆质量浓度为70%、75%和80%）倒入模具中。终凝24 h脱模后立即将组合试样放入PCD-90P型标准恒温恒湿养护箱，在温度为（20±5）℃和相对湿度为（95±5）%的条件下养护28 d。养护期满后将充—岩组合试样取出，得到3组料浆质量浓度分别为70%、75%和80%的充—岩组合试样，每个料浆浓度制备3个平行试样，共9个组合试样进行单轴压缩试验。充—岩组合试样制备过程如图3所示。
1.2 试验设备及方法
试验采用QKX-ZSZ-4000刚性伺服加载系统进行单轴压缩试验，并采用Micro-Ⅱ多通道声发射95探测系统（美国物理声学公司）和高速摄像机对组合试样破坏过程进行监测。试验前先将组合试样进行抛光处理（Yang et al，2022），然后在组合试样的充填体和结合面喷涂白漆。等自然风干后，将碳粉喷洒在组合试样结合面得到均匀的散斑（Lin et al，2021）。最后在充填体左侧与岩石右侧上下各放置2个（共4个）涂抹完耦合剂的声发射探头，并用胶皮进行固定（Chen et al，2023）。试验过程中，选择力控制加载方式，以0.5 kN/s的加载速率向试样施加力，直至应力下降至试样峰值强度的80%，停止试验。试验设备如图4所示。
2 试验结果与分析
2.1 组合体力学特征分析
为研究充—岩组合体中充填体与岩石的力学参数比值对组合体相应力学特征的影响规律，选取了3个具有代表性的力学参数比值作为分析指标，分别是充填体与岩石弹性模量的比值（ ），充填体与岩石峰值应变的比值（ ），充填体与岩石峰值强度的比值（ ）。通过控制充填体料浆质量浓度的大小来改变这些力学参数比值，深入探究力学参数比值变化对组合体相应力学特征的影响规律。
图5 不同试样单轴压缩应力—应变曲线Fig. 5 Uniaxial compression stress-strain curves of different samples |
2.2 弹性模量特征分析
图6揭示了充填体与岩石弹性模量比值 对组合体弹性模量的影响规律。随着弹性模量比 的增大，充—岩组合体的弹性模量呈先减小后增大的趋势，组合体的弹性模量在 为0.072时最高，在 为0.097时最低。组合体的弹性模量受到充填体弹性模量的限制，远低于岩石的弹性模量，接近充填体的弹性模量。当弹性模量比为0.072时，组合体弹性模量大于充填体的弹性模量，当弹性模量比为0.097和0.102时，组合体的弹性模量小于充填体的弹性模量。
2.3 变形特征分析
2.4 强度特征分析
图8揭示了充填体与岩石峰值强度比值 对组合体峰值强度的影响规律。以DIC全场应变演化为基准，观察到组合体中界面作为薄弱处，其先破坏导致组合体的整体峰值强度降低，使得组合体的峰值强度介于充填体与岩石的峰值强度之间，受到充填体峰值强度的限制。与岩石峰值强度相比，更接近于充填体的峰值强度。随着峰值强度比值 的增大，组合体的峰值强度呈先减小后增大的变化规律。当 为0.124时，组合体的峰值强度最高；当 为0.177时，组合体的峰值强度最低。当峰值强度比 增大至0.177时， 的增大促使界面应力集中现象的发生，加剧岩石破坏。岩石发生破坏后，强度较低的充填体在较小荷载增量下便发生破坏，最终引发组合体的整体性破坏，导致组合体整体峰值强度下降。当峰值强度比 继续增大至0.243时，由于充填体强度的增大，可以承担更多的荷载，虽然充填体先达到破坏强度发生破坏，但充填体发生破坏后，岩石强度较高，能够承担较大荷载后再发生破坏，因此组合体峰值强度有了小范围的上升。
综上所述，弹性模量比 为0.072~0.102，峰值应变比 为1.295~1.747，峰值强度比 为0.124~0.243范围内，因此选取 为0.072~0.097， 为1.295~1.387， 为0.124~0.177时，通过充填体与岩石的相互作用，可以较好地平衡变形协调与荷载分配，表现出较好的承载效果。
3 组合体破坏特征分析
3.1 组合体声发射特征分析
为了研究模量比 、应变比 和强度比 对组合体损伤演化规律的影响，采用声发射技术对单轴压缩过程中组合体的充填体侧和岩石侧进行监测，研究组合体内部裂纹萌生与扩展。模量比 、应变比 和强度比 均随料浆浓度的增大而增大，因此下述采用强度比 对组合体的损伤过程进行分析。图9所示为不同强度比 组合体的充填体侧与岩石侧破坏过程中声发射振铃计数、累计振铃计数和应力随时间的变化关系。
由图9可以看出，声发射振铃计数变化与应力的变化规律一致。当应力突降时，声发射振铃计数激增。当 为0.124时，充填体侧与岩石侧的声发射振铃计数聚集在28 s时应力突降处，以及峰值前80%至峰后阶段。只是声发射信号数值上存在差异，当应力达到峰值时，岩石侧的声发射信号数值更高，表明岩石受损更加严重。当 为0.177时，应力存在双峰现象。充填体侧的声发射振铃计数在应力达到第一次峰值时小幅上升，最大值出现在第二次峰值附近，表明充填体破坏主要在第二次峰值时。在应力达到第一次峰值时，岩石就发生破坏，存在较高的声发射振铃计数，随着荷载的继续增加，岩石破坏加剧，声发射振铃计数持续保持着活跃状态。当 为0.243时，随着 的增加，充填体侧与岩石侧的声发射信号明显增多，声发射振铃计数变化与 为0.177时相似，聚集在峰后阶段。当 为0.177~0.243时，充填体侧与岩石侧的声发射振铃计数在应力达到第二个峰值时仍保持显著活跃状态，存在声发射振铃计数激增现象，表明此时组合体存在二次破坏。
组合体的声发射累计振铃计数随着 的增加呈先降低后增加的趋势。当 为0.177时，声发射累计振铃计数最低，表明充填开采过程中，充填体与岩石进行耦合形成组合体时，由于组合体之间的相互作用关系，需要从充—岩组合体角度考虑充填体的最佳配比问题。
3.2 组合体裂纹扩展规律
图10 为0.124的组合体水平方向应变云图Fig. 10 Horizontal strain cloud maps of the composite with of 0.124 |
图11 为0.177的组合体水平方向应变云图Fig. 11 Horizontal strain cloud maps of the composite with of 0.177 |
如图10所示，当 为0.124时，孔隙压密阶段（Ⅰ），组合体中原有孔隙和微裂纹被压密，组合体整体受力均匀，应变随机分布（a点）。当应力增加到b点发生应力降现象，从应变云图中得出该现象源于界面处裂纹1的萌生。弹塑性阶段（Ⅱ），应力增加至c点，裂纹1沿界面持续扩展形成贯通性破坏面。应力达到峰值时（d点），岩石达到破坏强度，开始发生破坏，裂纹2自岩石顶部萌生并向界面处扩展。峰后阶段（Ⅲ），裂纹1和裂纹2继续发育扩展，使得岩石部分发生宏观破坏，应力下降。而充填体部分无明显表观破坏。
如图11所示，当 为0.177时，组合体在孔隙压密阶段（Ⅰ）呈现出与 为0.177时相同的裂纹扩展规律，a点应变随机分布，b点发生应力降，组合体界面出现裂纹1。在弹塑性阶段（Ⅱ），应力再次出现应力降（c点），应变云图显示是裂纹1扩展贯通界面所致。随着应力增大至第一个峰值（d点），岩石部分出现从界面萌生并向下扩展的裂纹2。在一次破坏阶段（Ⅲ），岩石部分裂纹2继续发育扩展。充填体部分出现由裂纹1发育引起的裂纹3。随着应力达到第二个峰值时（f点），裂纹2和裂纹3继续发育。在二次破坏阶段（Ⅳ），裂纹1、裂纹2和裂纹3持续扩展，导致组合体发生整体性破坏。
如图12所示，当 为0.243时，孔隙压密阶段（Ⅰ）应变随机分布（a点）。弹塑性阶段（Ⅱ）发生应力降，界面处出现裂纹1萌生，从而发生破坏（b点）。应力达到第一个峰值时（c点），裂纹1扩展贯通界面。在一次破坏阶段（Ⅲ），充填体部分发生破坏，出现由裂纹1扩展引起的裂纹2和裂纹3萌生（d点）。与应力达到第二个峰值时（e点），岩石部分开始发生破坏，出现了裂纹4。在二次破坏阶段（Ⅳ），裂纹1、裂纹2、裂纹3和裂纹4持续发育扩展，导致组合体发生整体性破坏，界面、充填体和岩石均发生破坏。
综上所述，充—岩组合体中充填体与岩石之间胶结效果较差，界面是最薄弱的地方，导致组合体受到荷载作用时，界面先于岩石和充填体发生破坏，后续岩石与充填体破坏的先后顺序受到模量比、应变比和强度比的复杂作用影响。组合体受到荷载作用时，荷载会按充填体和岩石各自的刚度和截面积进行分配。峰值应变强，能够通过发生变形延缓破坏的发生。当强度比 为0.124时，充填体峰值强度远低于岩石峰值强度，充填体与岩石模量比 为0.072，充填体弹性模量低，承担较小的载荷，而岩石弹性模量高，承担较大的荷载。同时，界面的破坏导致界面出现应力集中，加速岩石破坏。此时充填体与岩石应变比 为1.295，充填体峰值应变较大，可通过塑性变形耗散能量，延缓破坏，即使承载能力下降，仍保持一定的完整性，因此充填体表面无明显损伤。当强度比 为0.177时，充填体与岩石模量比 为0.097，岩石承担的荷载仍超过其破坏强度，且组合体峰值强度的增大促使界面处的应力集中发生，导致岩石先发生破坏。随后充填体单独承担荷载，由于充填体强度较低，不足以完全承载，使得组合体最终发生破坏。当 为0.243， 为0.102时，充填体的弹性模量再次增大，充填体自身强度低，因承担的荷载增多率先达到破坏强度，发生破坏。充填体破坏后，岩石单独承压，岩石的峰值强度大，荷载继续增加，达到破坏强度才开始破坏。
4 结论
（1）组合体的弹性模量、峰值应变和峰值强度均受到充填体的限制，与岩石相应参数相比，组合体参数更接近充填体相应的参数。随着充填体与岩石弹性模量比值 和峰值强度比值 的增大，充—岩组合体的弹性模量和峰值强度呈先减小后增大的变化趋势。随着充填体与岩石峰值应变比值 的增大，组合体的峰值应变呈单调递增的趋势。
（2）声发射振铃计数变化与应力的变化规律一致。当应力突降时，声发射振铃计数激增。随着充填体与岩石弹性模量比值 、峰值应变比值 和峰值强度比值 的增大，声发射累计振铃计数呈先减少后增加的变化趋势。
（3）组合体发生破坏时均是界面先发生贯穿破坏，但岩石与充填体破坏的先后顺序受到模量比、应变比和强度比的复杂作用的影响。当强度比 为0.124时，组合体破坏顺序是界面和岩石。当强度比 为0.177时，组合体破坏顺序是界面、岩石和充填体。当强度比 为0.243时，组合体破坏顺序是界面、充填体和岩石。
（4）在实际采矿过程中，充填体存在最佳配比问题。从力学特性、声发射与DIC的角度对组合体进行分析，模量比 为0.072~0.102，应变比 为1.295~1.747，强度比 在0.124~0.243范围内，选取 为0.072~0.097， 为1.295~1.387， 为0.124~0.177时，通过充填体与岩石的相互作用，可以较好地平衡变形协调与荷载分配，表现出较好的承载效果。
甘公网安备 62010202000672号
