Experimental research on low and strength characteristics of coal slime based cemented filling material
-
摘要:
为研究煤泥基充填材料的流动与强度特性,运用正交试验设计了以煤泥、高水、水泥、粉煤灰为组合的胶结充填材料,采用极差、方差分析及多元回归分析等方法,研究了煤泥水百分浓度(A)、高水添加比例(B)、水泥添加比例(C)、粉煤灰添加比例(D)等4个关键因素对充填材料的初凝时间、扩散度及抗压强度、抗剪强度的影响机制。试验结果表明:① 煤泥基充填材料的初凝时间为52~67 min,扩散度为42.0~76.5 cm,抗压强度为0.076~0.247 MPa,抗剪强度为0.033~0.139 MPa。对于初凝时间,影响显著性有B>C>A>D;对于扩散度,有A>D>B>C;对于抗压强度,有D>C>A>B;对于抗剪强度,有C>D>A>B。② 拟合出了A、B、C、D与流动性指标及强度指标的数学关系式,初凝时间与A、B成二次函数关系、与C、D均呈指数关系;扩散度与A、B、C均呈指数关系,与D呈三角函数关系;抗压强度与A、D均呈指数关系,与B、C呈二次多项式关系;抗剪强度与A、B呈二次多项式关系、与C、D均呈指数关系;在A、B、C、D等4个因素与4个参数独立拟合方程的基础上,建立4因素共同影响下煤泥基充填材料流动性及强度特性参数预计模型,并进行了验证。
Abstract:To study the flow and strength characteristics of coal slime based backfill materials, orthogonal tests were used to design one cemented backfill material combining coal slime, high-water-content materials, cement, and fly ash. By using range, analysis of variance, and multi-variate regression analysis, influences of four key factors on the initial setting time, bleeding ratio, diffusivity, compressive strength, and shear strength of the backfill material were studied. These four factors include coal slurry water mass per cent concentration (A), high water material mass / slurry water mass (B), cement mass / slurry water mass (C) and fly ash mass / slurry water mass (D). Results show that initial setting time, diffusivity, compressive strength, and shear strength of the backfill material are 52−67 min, 42.0−76.5 cm, 0.076−0.247 MPa, and 0.033−0.139 MPa, respectively, all distributed in wide ranges. For initial setting time, the four factors are listed in a descending order asB>C>A>Daccording to their influences; for the compressive strength, the four factors are ranked asD>C>A>B; for the shear strength, the four factors are listed in an order asC>D>A>B. Mathematical equations betweenA,B,C,Dand initial setting time, diffusivity, compressive strength, and shear strength were gotten, initial setting time andA,Bpresents quadratic function relation, andC,Dpresent exponential relation; diffusivity andA,B,Cpresents quadratic function relation, andDpresents trigonometric relation; compressive strength andA,Dpresents exponential relation, andB,Cpresents quadratic function relation; shear strength andA,Bpresents quadratic function relation, andC,Dpresents exponential relation.One prediction model for characteristic parameters, namely, fluidity and strength, of the coal slime based backfill material under action of various factors was established and calibrated.
-
Keywords:
- coal slime /
- backfill material /
- fluidity /
- compressive strength /
- shear strength
-
0. 引 言
我国坚硬顶板赋存煤层占全国赋存的30%左右 [1-2],受开采深度与围岩作用共同影响,易造成“悬而不跨”等现象。坚硬顶板一旦大面积突然破断,极易引发冲击地压、飓风、煤与瓦斯突出等煤岩动力灾害[3-6],严重制约工作面的正常回采,造成巨大的经济损失和人员伤亡。
爆破技术、水力压裂技术、静态膨胀致裂技术是矿井处理坚硬顶板的主要手段[7-9]。近年来,针对矿山坚硬顶板弱化问题,国内外学者取得了丰富的研究成果。于斌等[10]开展了坚硬顶板复合爆破定向造缝地面试验,揭示了相邻孔间裂缝扩展延伸特征。黄志增等[11]分析了深孔爆破对岩梁应力分布的影响作用,揭示了爆破切顶前后岩梁内部应力分布规律。HE等[12]基于复变函数法,建立了水力压裂破裂压力判据。康红普等[13]基于水力压裂工业性试验,获得了顶板水力压裂裂缝空间展布特征,建立了煤矿井下定向钻孔水力压裂成套技术。ZHAO等[14]提出了一种密集线性多孔控制无声开裂拆除剂定向压裂方法,分析了裂缝起裂扩展过程及定向压裂机理。谢雄刚等[15-16]经煤层静态膨胀致裂现场试验研究指出,选用合理的布孔参数可实现突出煤层的开裂。但在实际应用中,爆破技术无法保证围岩质量,致裂过程不够绿色安全,且易受爆炸物管理政策约束[17-18]。水力压裂技术受施工工艺限制,在卸压范围以及连续性方面存在不足[19]。静态破碎技术使用过程需要严格控制配比,且水化反应过程不受人为干预,易出现喷口等问题使其应用受限[20-21]。
套筒压裂技术通过体积膨胀产生的膨胀力来压裂岩石,能够在围岩环境较差的条件下进行施工[22],因其操作简单、可控性更强、无需封孔操作等优点而备受学者们的青睐。郑志涛等[23]基于PFC3D与室内试验研究了不同影响因素对套筒压裂裂纹扩展方向的影响,获得了定向致裂最优法。李丁[24]采用控制变量法开展了不同条件的套筒致裂试验,获得了围压、孔间距、布孔角度等参量对裂缝起裂扩展的影响规律。胡善超等[25]基于弹性力学、线弹性断裂力学推导了多孔套筒定向压裂机制。LI等[26]对比了套筒压裂与水力压裂的致裂效果,认为套筒压裂相较水力压裂所产生的裂缝更直。
上述学者针对套筒压裂取得了一系列重要的研究成果,为进一步研究套筒压裂机理提供了理论与试验基础。但以往的研究主要集中在岩石宏观裂缝扩展规律等方面,压裂过程中细观裂隙及能量演化规律的研究还不多见。为此,开展不同围压作用下套筒压裂室内试验以及PFC3D数值模拟试验,基于弹性力学推导套筒压裂力学机制,结合声发射特征信息探究围压对套筒压裂过程中岩石细观裂隙及能量演化过程的影响。
1. 试验与数值模拟方案
1.1 试件制作
天然岩体中富含节理与裂隙,较难用于制备大量尺寸与力学性质相近的试样。根据试验设备条件与本试验的研究目的,采用水泥砂浆作为类岩石材料开展室内试验研究。如图1所示,试样尺寸为100 mm×100 mm×100 mm,浇筑时在模具中心预埋实心圆柱以形成岩样孔眼,孔眼直径28 mm,孔深100 mm。套筒压裂试验与静态膨胀致裂试验均通过体积膨胀开裂岩石,二者压裂原理相同,两个试验本身具有一定的相似性,参考文献[27],确定水泥、砂子、石膏、水质量配比为2∶1∶2∶2,共进行5组常规单轴压缩试验,确定水泥砂浆标准试样力学性质如下:单轴抗压强度11.64 MPa,弹性模量10.22 GPa。试件自然养护48 h后脱模,在恒温恒湿条件下((20±1) ℃,湿度≥95%)养护30 d达到设计强度。试验前对试件进行密度测试,挑选表面完整、密度接近的试样作为研究对象。
1.2 试验方案
1.2.1 膨胀力校核试验
套筒壁采用柔性材料,在压裂过程中,注液系统将压裂液输送至密封套筒内部,引起套筒壁径向膨胀。该过程存在压力损耗,因此注液系统所提供的注液压力并不等于套筒所产生的膨胀力。
膨胀力是描述套筒压裂性能的重要指标,参考TANG等[28]所使用的静态膨胀剂膨胀力测定方法,采用无缝钢管以及静态应变电阻仪对套筒压裂过程中的膨胀力进行测定与校核。试验装置如图2所示,采用Q235型号钢管作为套筒外接管,钢管管径稍大于套筒,壁厚2 mm,弹性模量206 GPa,泊松比为0.3。采用DH3816静态应变电阻仪用于测量外接管轴向、径向应变。在套筒安全压力范围内($ \leqslant $35 MPa),设定注液速率为40 mL/min,试验过程钢管应变信号可被实时采集。套筒在钢管内所产生的径向膨胀力,可由式(1)进行计算[21]。
$$ P = \frac{{E({K^2} - 1)}}{{2 - \mu }}{\varepsilon _\theta } $$ (1) 式中:P为钢管所受套筒产生的径向膨胀压力,MPa;E为钢管的弹性模量,GPa;K为钢管内外半径比;εө为钢管径向应变;μ为钢管泊松比。
1.2.2 套筒压裂试验
套筒压裂试验是在不同应力条件下以恒定流量向套筒内注入压裂液,注液期间压裂液不会外流,通过套筒膨胀来致裂岩石的过程。本试验采用山东科技大学RLJW-2000型微机控制岩石伺服试验机,如图3所示,该试验机由轴向加载系统、剪切系统、孔隙水压系统构成。其中,轴向加载、剪切系统均采用液压伺服控制,用以双向围压加载。孔隙水压系统通过油水转换装置为套筒提供压裂液,同时对试验过程中的水压进行实时监测与记录。
为保证加载过程中的稳定性,减少加载过程对岩石的影响,轴向应力σ1、侧向应力σ2均以0.01 MPa/s加载至目标值,待围压稳定后再开启孔隙水压系统。1号试件出现明显裂隙时随即停止注液,2、3、4号试件注液压力达到设定值后再停止注液,并卸载水压与围压。试验选用清水做为压裂液,所施加的注液速率、外部应力条件见表1。
表 1 外部应力条件及注液速率Table 1. External stress conditions and liquid injection rate试件编号 σ1/
MPaσ2/
MPaσ1-σ2/
MPa注液速率/
(mL·min−1)终止应力/
MPa1 0 0 0 40 5.3 2 2 1.0 1.0 40 35 3 3 1.5 1.5 40 35 4 4 2.0 2.0 40 35 1.3 数值模拟方案
颗粒流离散元法(PFC)能够有效处理非连续介质问题,逼近真实材料的力学响应特征,被广泛应用在岩土工程等领域[29]。为揭示套筒压裂作用下岩石力学特性及演化规律,选用PFC3D开展数值模拟研究。
1.3.1 膨胀力校核数值模拟
基于膨胀力测定试验,建立如图4所示的数值模型。模型高100 mm,外径28 mm,内径16 mm。由于钢管形变量较小,可选用刚性墙体近似模拟,套筒则是选用具有相同大小的一组颗粒进行模拟,颗粒间采用接触黏结模型。为保证加载过程中颗粒间接触不会破坏,参考以往柔性加载的数值案例,将颗粒的法向、切向黏结强度设置为一个较大值[30-31],其余参数依据膨胀力校核试验所获取的膨胀力曲线进行标定。套筒膨胀过程压裂液不会外流,套筒注液压力可通过对套筒内边界施加等效集中力进行模拟。套筒膨胀力通过钢管的累计法向接触力表征,等效集中力依据实际注液压力进行拟合。
1.3.2 套筒压裂数值模拟
采用平行黏结模型建立与岩石试样具有相同尺寸的数值模型,通过单轴压缩与巴西劈裂室内试验标定岩石模型强度参数。图5a为套筒压裂模型的局部透视图,显示了岩样孔眼与套筒位置关系。为提高数值模拟精度,套筒颗粒相较于岩石颗粒更小,膨胀过程接触更加均匀。套筒边界处所施加的等效集中力依据不同应力条件下所获取的注液压力进行拟合。图5b为模型的俯视图,模型通过伺服机制沿x方向与y方向设置不同应力σx与σy,z方向设为自由边界。
2. 试验结果
2.1 模型标定结果
2.1.1 膨胀力校核结果
注液压力修正结果如图6a所示,可以看出,注液压力与套筒膨胀力随时间变化趋势一致,套筒膨胀力大致呈现指数式分布。受管壁材料性质的影响,膨胀力较注液压力存在压力损失,且随注液压力的增大而增大,最大损失量为9.79 MPa,约占最大总注液压力的28%。依据注液压力损失率可将套筒膨胀过程划分为3个阶段,即充填阶段(阶段Ⅰ)、增压阶段(阶段Ⅱ)、快速增压阶段(阶段Ⅲ)。在充填阶段,压裂液充填套筒内腔,注液压力损失率较高,膨胀力增速极为缓慢;在增压阶段,套筒内腔趋于饱和,膨胀力开始上涨,注液压力损失率开始快速下跌;在快速增压阶段,注液压力以及套筒膨胀力近乎线性增加,注液压力损失率趋于稳定。
膨胀力与注液压力关系如图6b所示,两者呈线性关系,其拟合曲线为
$$ y = 0.678x - 0.33 $$ (2) 式中,y为套筒膨胀力,MPa;x为注液压力,MPa。
PFC中所需的模拟时间与实际试验所需的时间并不一致,因此需要确定合理的循环数目以保证整个模拟过程的稳定。依据套筒注液压力曲线,拟合5组不同循环数目(1×104、2×104、3×104、4×104、5×104)的等效集中力曲线,所测得的膨胀力变化曲线如图7a所示。可以看出,不同循环数目下膨胀力曲线变化趋势一致,所测得的最大膨胀力大致相等,最大差异值为0.1 MPa。为定量描述膨胀力的离散情况,定义相邻循环间隔的膨胀力差为应力波动,不同循环数目下的应力波动如图7b所示。可以看出,随着循环数目的增加应力波动逐渐减小。当循环数为1×104时波动最为明显,最大波动为0.19 MPa;当循环数目为5×104时,最大应力波动为0.03 MPa,相较于前者下降约84.2%。增加循环数目会显著延长计算时间,循环数为5×104时基本符合数值模拟要求。
套筒细观参数与实际宏观参数存在一定的对应关系,基于膨胀力校核曲线采用试错法逐步调整,所得细观参数见表2。数值模拟膨胀力校核曲线如图8a所示,均可根据注液压力损失率大致划分为3个阶段。注液压力最大损失量9.69 MPa,误差仅为1%。数值模拟与室内试验的膨胀力与注液压力关系如图8b所示,两者均呈线性关系,且吻合度较高,因此该参数能够反映套筒膨胀过程。
表 2 校正后的套筒模型参数Table 2. Corrected sleeve model parameters接触类型 模型参量 参数 线性接触黏结模型 颗粒半径下限/ mm 0.6 颗粒粒径比 1.0 密度/(kg·m−3) 2000 有效模量/GPa 0.95 法向黏结强度/N 10300 切向黏结强度/N 10300 刚度比 1.1 摩擦因数 0.58 2.1.2 岩石微观参数标定结果
颗粒强度以及接触参数直接影响模型宏观破坏形态,为确定岩石细观参数,采用平行黏结模型建立尺寸为ø50 mm×100 mm的圆柱体试件。综合考虑颗粒总数以及计算机计算效率,确定颗粒半径。用于参数标定以及套筒压裂试验的数值模型均符合模型最小尺度颗粒数标准[32-33]。基于单轴压缩试验,采用试错法对岩石模型细观参数进行标定。对于平行黏结模型,黏结抗拉强度、黏聚力强度以及黏结弹性模量、黏结刚度比分别决定了材料的强度性质与变形性质,在试错过程中通常表现的更为敏感,应作为主要调试对象[34]。经多次校核,所得最优细观参数见表3。
表 3 校正后的岩石模型参数Table 3. The corrected rock model parameters接触类型 模型参量 参数 线性模型 颗粒半径下限/mm 1.5 颗粒粒径比 1.66 密度/(kg·m−3) 2500 弹性模量/GPa 1.1 刚度比 1.76 摩擦因数 0.577 平行黏结模型 弹性模量/GPa 1.95 刚度比 1.63 抗拉强度/MPa 2.2 黏聚力/MPa 3.8 摩擦角/(°) 35 室内试验与PFC模拟结果对比如图9与表4所示,可以看出:①岩石应力-应变曲线基本吻合,由于PFC建模初期消除了颗粒间的不平衡力,因此无法模拟岩石的压密阶段;②宏观破坏形态基本一致,均沿轴向产生贯通岩石表面的宏观裂缝,呈现脆性破坏特征;③峰值强度、弹性模量高度一致,误差仅为1%与2%。由此可知,该参数能够捕获岩石整体力学特性,后续模拟研究将基于此参数进行。
表 4 试验结果对比Table 4. comparison of experimental results参量 室内试验 数值模拟 单轴压缩强度/MPa 11.64 11.76 弹性模量/GPa 10.22 10.43 不同应力条件下的套筒注液压力曲线如图10所示,从图中可以看出,各试验组注液压力随注液时间变化规律相似。随着注液压力的增加,岩石表面产生宏观裂缝,钻孔体积出现扩增,因此试验组2、3、4号注液压力出现明显应力波动。而在试验组1号中,当岩石出现宏观裂缝时随即停止注液,因此注液压力并未表现出波动性。
根据注液压力曲线变化特征可将套筒压裂试验分为2个阶段。其中,$a \to b$为注液阶段,$b \to c$为卸压阶段。岩石的破坏主要集中在ab阶段,因此在数值模拟研究中,仅对ab段注液压力曲线进行拟合,bc段不做研究。
2.2 岩石破坏形态
套筒压裂过程伴随着剧烈声响,岩石破坏表现出脆性破坏特征,试验结束后岩石被压裂为多个部分。不同围压条件下岩石物理试验与数值模拟破坏形态对比如图11所示,可以看出,物理试验与数值模拟所获取的岩石裂缝形态基本一致。无围压条件下,1号试件仅受套筒膨胀力作用,形成了3条互成一定夹角的主裂缝,岩石试件被拉裂为3个部分。对于不同应力差值的2、3、4号试件,均形成了贯穿岩石表面的裂缝网络。其中,2号岩石试件裂缝数目最多,裂纹分布最复杂。随着围压的增大,同一注液压力作用下裂缝数目逐渐减少,裂缝网络趋于简单。受外部应力及岩石材料性质影响初生裂缝表现出一定的偏转,但其扩展路径最终趋向于大主应力方向。随着套筒膨胀力的增加,初生裂缝开度逐渐增大,裂缝面岩块呈现碎块状剥落,同时岩石在小主应力方向也形成了一定长度的次生裂缝。
颗粒间通过平行黏结接触互相关联,接触的力学行为是离散元计算方法的核心,当平行黏结接触在外力作用下发生破坏时即认为产生了微裂纹,可依据接触失效类型将裂纹划分为张拉型裂纹与剪切型裂纹。图12以两向应力差1.5 MPa时为例,分别截取等效集中力加载至0.37Pmax、0.48Pmax、0.81Pmax、Pmax时微观裂纹的分布状况。其中,红色圆盘代表张拉型裂纹,绿色圆盘代表剪切型裂纹。
在套筒压裂初期,岩石内部无微裂缝产生,这是由于等效集中力被套筒壁与岩石所消耗,整个模型处于弹性阶段;当套筒膨胀力达到0.37Pmax时,岩石发生起裂,裂纹主要集中在套筒周围,且基本沿大主应力方向产生;当套筒膨胀力达到0.48Pmax时,模型沿大主应力方向形成了贯通表面的宏观裂缝,断裂面以张拉破坏为主,仅存在少量剪切型裂纹;当套筒膨胀力达到0.81Pmax时,模型在边界处形成剪切裂纹,断裂面处裂纹进一步发育,孔眼位置裂纹进一步萌生;当套筒膨胀力达到Pmax时,等效集中力达到预设值,模型沿小主应力方向形成贯通表面的宏观裂缝,裂缝网络整体呈现“十字形”分布。
图13为不同围压条件下PFC数值模型套筒压裂过程中内部裂纹数量及不同种类裂纹所占比例统计变化。可以看出,无围压条件下,1号试件仅受套筒膨胀力作用,裂隙类型以张拉型裂纹为主,剪切裂纹极少,仅占裂纹总数的3.38%。有围压条件下,2、3、4号试件内部张拉型裂纹占比分别为93.05%、90.76%、88.49%,可以看出,套筒压裂过程中试件内部接触以张拉破坏为主,岩石内部张拉作用更为明显。随着两向应力差的逐渐增大,岩石内部裂纹总数及张拉型裂纹占比均逐渐减小,4号试件相较2号试件裂纹总数下降了23.76%,这是由于外部应力限制了套筒的膨胀,岩石破坏更加困难。套筒压裂过程中,岩石剪切裂纹占比随两向应力差的增大而增大,最大提升4.56%,这表明同一膨胀力作用下,两向应力差越大,岩石内部剪切破坏越明显。
2.3 声发射模拟结果
声发射法因其可监测岩石材料内部裂缝的动态演化特征,揭示不同变形阶段的岩石破损机制为而被广泛应用。在数值模拟过程中,一旦接触力超过其相应的黏结强度,接触便会断裂,即对应一次声发射事件。根据此原理,可采用FISH语言将模拟过程中每一时刻(时步)的声发射数目实时记录,并在黏结键破坏的地方生成一个裂隙标记,来实现套筒压裂过程中声信号空间分布特征的监测[35]。值得注意的是,该标记不参与力学计算,并随后续计算不断更新其相应的位置。
无围压条件下累计裂缝数目与声发射关系如图14a所示,可以看出,在致裂初期几乎没有声发射信号,整个岩石模型处于弹性阶段。随着套筒膨胀力的逐渐增大,岩石开始损伤,但AE振铃强度较弱且相对分散。当AE振铃次数达到峰值时,裂纹累计数目曲线出现瞬时激增,该过程所产生的裂纹数目约占总裂纹数目的72.19%,内部裂纹发生大规模扩展并形成宏观裂隙。整体来看,无围压条件下声发射特征可分为平静期与活跃期,且呈现单峰值分布,内部裂纹仅出现一次大规模扩展。在致裂末期,累计裂纹数目继续增加,但增速趋于平缓,AE振铃次数显著降低,该过程试件断裂为多个部分,并进行着摩擦运动。
有围压条件下累计裂缝数目与声发射关系如图14b、14c、14d所示,可以看出,其声发射特征与无围压条件相似,均存在着平静期与活跃期。但不同点在于,有围压条件下的声发射特征存在多个峰值,岩石内部裂隙发育以及沟通交替存在。结合图12裂隙演化过程可知,试件内部的裂纹贯通并不是一次性完成的,而是随着套筒膨胀力的增加表现出一定的周期性。压裂过程中,累计裂纹曲线出现多次瞬时激增,近似呈现阶梯状上升趋势,其中2号试件声发射活动表现最为剧烈,共出现3次瞬时激增。
结合套筒压裂过程声发射特征可知,套筒压裂技术能够有效弱化岩体,压裂过程存在多个破坏阶段,由于压裂过程压裂液不会外流,其可控性更强,压裂过程可对膨胀力实现人为干预。
3. 讨 论
3.1 岩石破坏力学机制
为揭示双向围压条件下无限大弹性体中套筒膨胀力学机制,建模时做以下假设[36]:①围岩岩体为各向同性,线弹性,无蠕变或粘性行为;②钻孔周围围岩应力状态为静水压力;③钻孔沿深度方向围岩性质相同,将计算简化为平面应变问题;④套筒所产生的膨胀压均布于孔壁。
建立如图15a所示的套筒压裂二维力学模型,模型分沿x方向与y方向施加有远场应力σH与σV,钻孔半径为a,钻孔内壁施加均匀的膨胀力P0,钻孔周围一点n距孔心距离为ρ,与远场应力σH方向夹角记为φ。
根据叠加原理,图15a的应力分布可转为图15b与图15c两种受力状态叠加,即:
$$ \left\{ \begin{gathered} \sigma _{\text{φ} } ^{\mathrm{a}} = \sigma _{\text{φ} } ^{\mathrm{b}} + \sigma _{\text{φ} } ^{\mathrm{c}} \\ \tau _{\text{ρφ} }^{\mathrm{a}} = \tau _{{\text{ρφ} } }^{\mathrm{b}} + \tau _{{\text{ρφ} } }^{\mathrm{c }} \\ \end{gathered} \right. $$ (3) 式中,$ \sigma _{\text{φ} } ^{\mathrm{a}} $、$ \sigma _{\text{φ} } ^{\mathrm{b}} $、$ \sigma _{\text{φ} } ^{\mathrm{c}} $分别为图15a、15b、15c中n点处的环向应力,$ \tau _{{\text{ρφ} } }^{\mathrm{a}} $、$ \tau _{{\text{ρφ} }}^{\mathrm{b }}$、$ \tau _{{\text{ρφ} } }^{\mathrm{c }}$分别为图15a、15b、15c中n点处的切向应力。
对于图15b,双向围压条件下单孔无限大板钻孔周围环向应力、切向应力在极坐标下的基尔斯解[37]为
$$\begin{split}& \sigma _{\text{φ} } ^{\mathrm{b}} = - \frac{1}{2}({\sigma _{\mathrm{H}}} + {\sigma _{\mathrm{V}}})\left[ {1 + {{(\frac{a}{\rho })}^2}} \right] + \\& \frac{1}{2}({\sigma _{\mathrm{H}}} - {\sigma _{\mathrm{V}}})\left[ {1 + 3{{\left(\frac{a}{\rho }\right)}^4}} \right]\cos \;2\varphi \end{split} $$ (4) $$ \tau _{{\text{ρφ} } }^{\mathrm{b}} = \frac{1}{2}({\sigma _{\mathrm{H}}} - {\sigma _{\mathrm{V}}})\left[ {1 - {{\left(\frac{a}{\rho }\right)}^2}} \right]\left[ {1 + 3{{\left(\frac{a}{\rho }\right)}^2}} \right]\sin\; 2\varphi $$ (5) 对于图15c,根据无限大弹性体中圆筒受力模型,圆筒内均匀受压,外半径趋于无穷大,对于这一问题的拉梅解[37]为:
$$ \left\{ \begin{gathered} \sigma _{\text{φ} }^{\mathrm{c}} = {\left(\frac{a}{\rho }\right)^2}{P_0} \\ \tau _{{\text{ρφ} } }^{\mathrm{c}} = 0 \\ \end{gathered} \right. $$ (6) 为简化分析,引入两向应力差Δσ,其计算方法为Δσ = σH−σV。引入特征参量λ作为距离ρ与钻孔半径a的比值,计算方法为λ = ρ/a,其中$\lambda \geqslant 1$。综合式(3)—式(6),可得图15a的应力分量表达式为:
$$ \left\{\begin{split} & \sigma _{\text{φ} } ^{\mathrm{a}} = - \frac{1}{2}(\Delta \sigma + 2{\sigma _V})\left[ {1 + \frac{1}{{{\lambda ^2}}}} \right] +\\ & \frac{1}{2}\Delta \sigma \left[ {1 + 3\frac{1}{{{\lambda ^4}}}} \right]\cos \;2\varphi + \frac{1}{{{\lambda ^2}}}{P_0} \\ \; & \tau _{{\text{ρφ} }}^{\mathrm{a}} = \frac{1}{2}\Delta \sigma \left[ {1 - \frac{1}{{{\lambda ^2}}}} \right]\left[ {1 + \frac{3}{{{\lambda ^2}}}} \right]\sin\; 2\varphi \\ \end{split} \right. $$ (7) 取σV = 1 MPa,Δσ = 2 MPa,φ = 0°作为基本参数,综合式(6)易知该工况条件下n点出环向应力解析式为:
$$ \sigma _{\text{φ} }^ * = - 2\left[ {1 + \frac{1}{{{\lambda ^2}}}} \right] + \left[ {1 + 3\frac{1}{{{\lambda ^4}}}} \right] + \frac{{{P_0}}}{{{\lambda ^2}}} $$ (8) 式中,$ \sigma _{\text{φ} }^ * $为所设工况n点处环向应力。
以P0 = 10 MPa为例,远场应力以及套筒膨胀力所产生的环向应力随参量λ的变化规律如图16所示。从图16可以看出,环向应力呈现对称式分布,靠近钻孔的位置张拉作用更加明显,随着λ的增加环向应力迅速衰减并趋于稳定。受远场应力作用,套筒膨胀力作用范围有限,环向应力存在受拉区间(−λ0 , λ0)以及受压区间(−∞ , −λ0)∪(λ0 , ∞)。其中参量λ0为受拉区间临界值,与套筒膨胀力P0存在一定协调关系,即套筒膨胀力P0越大,环向应力受拉区间范围越大。
当远场应力与套筒膨胀力共同作用下位置点处环向应力大于0时,该点受力状态为受拉,即受拉区间临界值应满足:
$$ \sigma _{\text{φ} } ^ * \geqslant 0 $$ (9) 综合式(8)与式(9)可知受拉区间临界值λ0表达式为
$$ {\lambda _0} = \frac{{\sqrt {2{P_0} - 3 + \sqrt {4P_0^2 - 12{P_0} + 33} } }}{2} $$ (10) 受拉区间临界值λ0随套筒膨胀力P0变化规律如图17所示。从图17可以看出,曲线增速随着套筒膨胀力的增加逐渐变缓。经拟合,参量λ0与套筒膨胀力大致呈现对数分布。当套筒膨胀力从10 MPa增加至40 MPa时,其有效致裂范围约为6倍孔径,相较于前者受拉区间范围提升了79.56%。
由式(7)可知,环向应力与应力差Δσ、参量λ、夹角φ、膨胀力P0密切相关, 切向应力仅与应力差Δσ、参量λ、夹角φ密切相关,而与套筒膨胀力P0无关。不同应力差作用下环向应力随参量λ变化规律如图18a所示。从图18a可以看出,同一膨胀力作用下,随着应力差的逐渐增大,不同位置点处的环向应力出现不同程度的变化。钻孔周围环向应力改变幅度相较其他位置更为明显,即随着应力差的逐渐增大,钻孔周围受围压限制程度更大,裂隙减少现象将更为突出,这与图11的数值模拟结果相一致。不同应力差作用下切向应力随夹角变化规律如图18b所示。从图18b中可以看出,无围压条件下,切向应力恒定为0,岩石内部不受剪切作用。有围压条件下,切向应力随着夹角的变化呈现一定的周期性,当夹角取值为45°或135°时,切向应力取得最大值。同一夹角下,切向应力随着应力差的增大而增大。这一理论结果可对图13模拟结果进行解释,应力差的增加使得岩石内部承受了更大的切向应力,因此岩石内部剪切型裂纹占比将有所上升。无围压条件下,岩石内部仅存在张拉作用,少量的剪切型裂纹可能来自裂缝面间的摩擦。
3.2 套筒压裂过程中岩石能量演化规律
假设岩石套筒压裂过程与外界没有热交换产生,基于能量守恒定律,岩石的破坏行为是输入能、弹性能、耗散能累计与转化的结果[38]。如图19所示,在PFC能量计算中,应变能主要包括平行黏结应变能与线性应变能,耗散能包括滑动摩擦能、阻尼能、缓冲器消耗能、动能等[39]。在图5所示的数值模型中,等效集中力对岩石所做的功可视为系统的总输入能,输入能又可转化为套筒弹性能以及岩石弹性能。随着等效集中力的加载,岩石内部产生微裂隙,进而形成耗散能。由于等效集中力加载过程符合准静态过程,因此耗散能主要包括摩擦能与阻尼能。
不同围压条件下套筒压裂能量演化规律如图20所示,从图20可以看出,不同岩石的能量演化规律相似。在压裂初期,岩石内部各类型能量几乎为零,该过程输入能增加极为缓慢。随着输入能的不断增加,套筒与钻孔紧密接触,岩石进入弹性阶段。该过程中,输入能一方面被套筒自身膨胀所消耗,另一方面通过套筒壁间接转移至岩石内部。当能量积聚至极限承载值时,岩石进入破坏阶段,耗散能近乎线性增加,宏观裂隙面迅速形成。
在套筒压裂末期,岩石内部发生宏观破坏,钻孔体积出现扩容,套筒将进一步膨胀。由于PFC中能量的计算遵循力−位移定律,因此输入能、弹性能均出现不同程度激增。受围压限制,不同应力条件下岩石内部的输入能大小并不相等,同一膨胀力作用下,围压越大岩石的输入能越小,4号岩石相较2号岩石输入能约下降47.92%。对比不同应力条件下岩石能量占比可以发现,套筒压裂过程中耗散能占比均大于弹性能,且不同应力条件下能量平均增长速率各不相同。能量平均增长速率计算方法为某一阶段内所产生的某一能量与该阶段所用循环数之比,定义耗散能产生时所对应的循环数为初始循环数,不难看出,套筒压裂末期耗散能增量相较弹性能更多,且曲线平均增长速率相较弹性能更快。对于1、2、3、4岩石,耗散能平均增长速率约为弹性能的1.48、1.27、1.33、1.55倍,能量释放更迅速。
4. 结 论
1)套筒注液过程存在压力损失,依据注液压力损失率可将套筒注液过程大致划分为3个阶段。在快速增压阶段(阶段Ⅲ)中,压力损失率趋于稳定,注液压力与膨胀力变化趋势基本一致,两者呈线性关系。
2)套筒压裂过程岩石以张拉破坏为主,随着应力差的增大,岩石剪切破坏越明显。围压对裂缝网络的形成有着抑制作用,并主导着裂缝的扩展路径。有围压条件下,初生裂缝总是沿着大主应力方向扩展,次生裂缝沿着小主应力方向产生,裂缝网络大致呈现“十字型”分布。
3)声发射信号特征可分为平静期与活跃期,内部裂隙发育以及沟通交替存在。无围压条件下,声发射特征呈单峰值分布;有围压条件下,声发射特征呈现多峰值分布,压裂过程存在多个破坏阶段。
4)套筒膨胀力作用范围有限,靠近钻孔的位置张拉作用更明显。环向应力存在受拉区间,其作用范围与套筒膨胀力呈现对数关系;切向应力随夹角φ的变化呈现一定的周期性,同一夹角下,切向应力随应力差的增大而增大。
5)同一套筒膨胀力作用下,围压越大岩石的输入能越小,岩石开裂越困难;宏观裂缝的形成引发了钻孔扩容现象,岩石的输入能、弹性能、耗散能均出现上升,但耗散能释放更为迅速,其能量占比以及平均增长速率均高于弹性能。
-
表 1 因素水平
Table 1 Factor levels
水平组数 A/% B C D 1 40.00 0.030 0 0.040 0 0.040 0 2 33.75 0.037 5 0.080 0 0.080 0 3 28.75 0.045 0 0.120 0 0.120 0 表 2 正交试验结果
Table 2 Summary of orthogonal test results
试验组数 影响因素 试验结果 A/% B C D 初凝时间/min 扩散度/cm 抗压强度/MPa 抗剪强度/MPa 1 40.00 0.030 0 0.040 0 0.040 0 57.000 42.000 0.107 0.056 2 40.00 0.037 5 0.080 0 0.080 0 56.000 50.800 0.161 0.095 3 40.00 0.045 0 0.120 0 0.120 0 52.000 40.000 0.234 0.369 4 33.75 0.030 0 0.080 0 0.120 0 58.000 57.800 0.247 0.109 5 33.75 0.037 5 0.120 0 0.040 0 67.000 61.300 0.086 0.043 6 33.75 0.045 0 0.040 0 0.080 0 53.000 62.200 0.076 0.052 7 28.75 0.030 0 0.120 0 0.080 0 62.000 76.500 0.115 0.139 8 28.75 0.037 5 0.040 0 0.120 0 55.000 67.200 0.102 0.036 9 28.75 0.045 0 0.080 0 0.040 0 52.000 68.300 0.105 0.033 表 3 试验数据极差分析
Table 3 Ranged analysis of test data
因素 初凝时间/min 扩散度/cm A B C D A B C D k1 55.000 59.000 55.000 58.667 44.267 58.767 57.133 57.200 k2 59.333 59.333 55.333 57.000 60.433 59.767 58.967 63.167 k3 56.333 52.333 60.333 55.000 70.667 56.833 59.267 55.000 R 4.333 7.000 5.333 3.667 26.400 2.933 2.133 8.167 因素 抗压强度/MPa 抗剪强度/MPa A B C D A B C D k1 0.168 0.157 0.095 0.099 0.173 0.101 0.048 0.044 k2 0.136 0.116 0.171 0.117 0.068 0.058 0.079 0.095 k3 0.107 0.138 0.145 0.194 0.069 0.151 0.184 0.171 R 0.061 0.040 0.076 0.095 0.105 0.093 0.136 0.127 注:k为因素各个水平下指标总和的平均数;R为极差,表示因子对结果的影响幅度。 表 4 试验数据方差分析
Table 4 Analysis of test data ANOVA
变异来源 偏差平方和 自由度 方差 F值 Fa 显著水平 A 29.5319 2 14.766 1.187 F0.01(2,4)=18 — B 93.5319 2 46.7659 3.759 F0.05(2,4)=6.944 o C 53.5319 2 26.766 2.151 F0.1(2,4)=4.325 o 误差e1 20.2324 2 10.1162 0.813 F0.25(2,4)=2 修正误差e2 49.7643 4 12.4411 总和 196.828 注:Fa为F分布临界值;***表示影响极为显著;**表示影响较为显著;*表示影响一般显著;o表示有显著影响;—表示无显著影响。 表 5 试验数据方差分析
Table 5 Analysis of test data ANOVA
变异来源 偏差平方和 自由度 方差 F值 Fa 显著水平 A 1063.0560 2 531.528 99.466 F0.01(2,4)=18 *** B 13.3581 2 6.679 1.250 F0.05(2,4)=6.944 — D 107.1497 2 53.5749 10.026 F0.1(2,4)=4.325 ** 误差e1 8.0169 2 4.0085 0.750 F0.25(2,4)=2 修正误差e2 21.3750 4 5.3438 总和 1191.5810 表 6 试验数据方差分析
Table 6 Analysis of test data ANOVA
变异来源 偏差平方和 自由度 方差 F值 Fa 显著水平 A 0.0055 2 0.0027 2.250 F0.01(2,2)=99 — C 0.0090 2 0.0045 3.750 F0.05(2,2)=19 o D 0.0154 2 0.0077 6.417 F0.1(2,2)=9 o 误差e1 0.0025 2 0.0012 1.000 F0.25(2,2)=3 总和 0.0320 表 7 试验数据方差分析
Table 7 Analysis of test data ANOVA
变异来源 偏差平方和 自由度 方差 F值 Fa 显著水平 A 0.0218 2 0.010 9 1.1010 F0.01(2,6)=10.925 — C 0.0303 2 0.015 2 1.5354 F0.05(2,6)=5.143 — D 0.0245 2 0.012 3 1.2424 F0.1(2,6)=3.463 — 误差e1 0.0132 2 0.006 6 0.6667 F0.25(2,6)=1.762 修正误差e2 0.0594 6 0.009 9 总和 0.09 表 8 拟合关系
Table 8 Fitting relationship
因素 指标 拟合关系式 拟合度R2 A 初凝时间 $ {t_1} = - 1\;150{A^2} + 778.5A - 72.47 $ 0.99 扩散度 $ {S_1} = 114.866 - 13.355{{\text{e}}^{A/0.24}} $ 0.99 抗压强度 $ {\sigma _{Y1}} = 0.33 - 0.497{{\text{e}}^{ - A/0.36}} $ 0.99 抗剪强度 $ {\sigma _{J1}} = {\text{15}}{\text{.22}}{A^2} - 9.537A + 1.553 $ 0.99 B 初凝时间 $ {t_2} = - 65\;190{B^2} + 4\;444B - 15.67 $ 0.99 扩散度 $ {S_2} = 59.89{{\text{e}}^{ - {{{\text{[(}}B - 0.035\;63{\text{)/}}0.040\;91{\text{]}}}^2}}} $ 0.99 抗压强度 $ {\sigma _{Y2}} = 551.1{B^2} - 42.53B + 0.936\;3 $ 0.99 抗剪强度 $ {\sigma _{J2}} = 1\;215{B^2} - 87.77B + 1.641 $ 0.99 C 初凝时间 $ {t_3} = {\text{54}}{\text{.976\;7}} + 0.023\;3{{\text{e}}^{(C - 0.04)/0.014\;7}} $ 0.99 扩散度 $ {S_3} = 59.328\;44 - 13.483\;77{{\text{e}}^{ - C/0.022\;05}} $ 0.99 抗压强度 $ {\sigma _{Y3}} = - 31.875{C^2} + 5.725C - 0.083 $ 0.99 抗剪强度 $ {\sigma _{J3}} = {\text{0}}{\text{.034\;95}} + 0.013\;1{{\text{e}}^{(C - 0.04)/0.032\;87}} $ 0.99 D 初凝时间 $ {t_4} = 67 - 6.944\;4{{\text{e}}^{D/0.219\;4}} $ 0.99 扩散度 $ {S_4} = 63.21\sin \left( {11.94D + 0.653\;6} \right) $ 0.99 抗压强度 $ {\sigma _{Y4}} = {\text{0}}{\text{.093\;8}} + 0.005\;5{{\text{e}}^{(D - 0.04)/0.027\;5}} $ 0.99 抗剪强度 $ {\sigma _{J4}} = - {\text{0}}{\text{.062\;8}} + 0.106\;8{{\text{e}}^{(D - 0.04)/0.101\;94}} $ 0.99 表 9 校验方案设计
Table 9 Fitting relationship
测试件编号 .煤泥水百分
浓度A/%高水添加
比例B水泥添加
比例C粉煤灰添加
比例DTEST-1 35.00 0.035 0 0.100 0 0.100 0 TEST-2 30.00 0.028 0 0.050 0 0.050 0 TEST-3 39.00 0.025 0 0.050 0 0.100 0 表 10 校验结果分析
Table 10 Fitting relationship
试件编号 材料 初凝时
间/min扩散
度/cm抗压强
度/MPa抗剪强
度/MPaTEST-1 预测值 57.290 60.200 0.173 0.056 试验值 60.010 63.980 0.159 0.062 误差/% 4.53 5.91 8.74 9.68 TEST-2 预测值 55.900 69.900 0.105 0.095 试验值 52.060 63.900 0.096 0.086 误差/% 7.38 9.39 9.69 10.47 TEST-3 预测值 50.910 51.330 0.241 0.370 试验值 48.330 53.710 0.263 0.411 误差/% 5.34 4.43 8.52 9.88 -
[1] 郑丽辉,周跃进,张雷鸣,等. 矸石聚合物充填材料承载特性试验研究[J]. 煤炭科学技术,2022,50(9):222−229. ZHENG Lihui,ZHOU Yuejin,ZHANG Leiming,et al. Experimental research on load-bearing characteristics of gangue polymer backfilling material[J]. Coal Science and Technology,2022,50(9):222−229.
[2] 林海飞,翟雨龙,李树刚,等. 基于正交设计的“固-气”耦合相似材料力学与渗透特性试验研究[J]. 煤炭学报,2016,41(3):672−679. LIN Haifei,ZHAI Yulong,LI Shugang,et al. Research on mechanics and permeability characteristics of similar material of solid-gas coupling based on orthogonal design[J]. Journal of China Coal Society,2016,41(3):672−679.
[3] 常庆粮,袁崇亮,王永忠,等. 膏体充填综采台阶煤壁稳定性半凸力学分析[J]. 中国矿业大学学报,2022,51(1):46−55. doi: 10.3969/j.issn.1000-1964.2022.1.zgkydxxb202201005 CHANG Qingliang,YUAN Chongliang,WANG Yongzhong,et al. Semi-convex mechnical analysis on stability of step coal wall in fully mechanized mining with paste filling[J]. Journal of China University of Mining and Technology,2022,51(1):46−55. doi: 10.3969/j.issn.1000-1964.2022.1.zgkydxxb202201005
[4] 黄 鹏,张吉雄,郭宇鸣,等. 深部矸石充填体黏弹性效应及顶板时效变形特征[J]. 中国矿业大学学报,2021,50(3):489−497. HUANG Peng,ZHANG Jixiong,GUO Yuming,et al. Viscoelastic effect of deep gangue backfill body and time-dependent deformation characteristics of roof in deep mining[J]. Journal of China University of Mining and Technology,2021,50(3):489−497.
[5] 秦海忠,李 佳. 氟石膏基新型膏体充填材料配比试验研究[J]. 有色金属工程,2022,12(9):144−153. doi: 10.3969/j.issn.2095-1744.2022.09.020 QIN Haizhong,LI Jia. Study on the proportion of a new fluorgypsum-based paste filling material[J]. Nonferrous Metals Engineering,2022,12(9):144−153. doi: 10.3969/j.issn.2095-1744.2022.09.020
[6] 刘铁军,万擎宇,许高锋. 某铅锌矿全尾膏体充填工艺的选择及分析[J]. 有色冶金设计与研究,2022,43(4):1−4. doi: 10.3969/j.issn.1004-4345.2022.04.001 LIU Tiejun,WAN Qingyu,XU Gaofeng. Selection and analysis of unclassified tailings paste filling process for a lead-zinc mine[J]. Nonferrous Metals Engineering and Research,2022,43(4):1−4. doi: 10.3969/j.issn.1004-4345.2022.04.001
[7] 轩召军,常庆粮. 膏体充填开采底板破坏机理及分区治理技术[J]. 煤矿安全,2022,53(8):210−217. XUAN Zhaojun,CHANG Qingliang. Study on floor failure mechanism of paste backfill mining and division area control technology[J]. Safety in Coal Mines,2022,53(8):210−217.
[8] 余伟健,万 幸,刘芳芳,等. 红土膏体充填材料及其物理特性试验研究[J]. 煤炭科学技术,2021,49(2):61−68. YU Weijian,WAN Xing,LIU Fangfang,et al. Experimental study on red clay paste backfilling material and its physical characteristics[J]. Coal Science and Technology,2021,49(2):61−68.
[9] 孙希奎. “三下”采煤膏体充填开采技术研究[J]. 煤炭科学技术,2021,49(2):61−68. SUN Xikui. Research on paste backfilling mining technology of coal mining under buildings,water bodies and railways[J]. Coal Science and Technology,2021,49(2):61−68.
[10] 王永岩,于卓群,崔立桩. 不同含水率膏体充填材料的单轴压缩试验研究[J]. 煤炭科学技术,2022,50(6):219−224. WANG Yongyan,YU Zhuoqun,CUI Lizhuang. Experimental stsudy on compressive behavior of cemented paste backfilling material with different water contents[J]. Coal Science and Technology,2022,50(6):219−224.
[11] 郝 喆. 小窑采空区地质灾害综合评价研究[J]. 地质与资源,2019,28(4):394−400. doi: 10.3969/j.issn.1671-1947.2019.04.014 HAO Zhe. Comprehensive assessment on the geological disasters in small pit goaf[J]. Geology and Resources,2019,28(4):394−400. doi: 10.3969/j.issn.1671-1947.2019.04.014
[12] 冯国瑞,侯水云,梁春豪,等. 复杂条件下遗煤开采岩层控制理论与关键技术研究[J]. 煤炭科学技术,2020,48(1):144−149. FENG Guorui,HOU Shuiyun,LIANG Chunhao,et al. Basic theories and key technologies of rock strata control for residual coal resources mining under complex conditions[J]. Coal Science and Technology,2020,48(1):144−149.
[13] 冯国瑞,李庆东,戚庭野,等. 超声波法评价矸石膏体的静态抗离析性能[J]. 采矿与安全工程学报,2020,37(6):1231−1237. FENG Guorui,LI Qingdong,QI Tingye,et al. Assessing static stability of cemented coal gangue backfill with ultrasonic pulse velocity method[J]. Journal of Mining and Safety Engineering,2020,37(6):1231−1237.
[14] 毋林林,康天合,尹 博,等. 粉煤灰膏体充填材料水化放热特性的微量热测试与分析[J]. 煤炭学报,2015,40(12):2801−2806. WU Linlin,KANG Tianhe,YIN Bo,et al. Microcalorimetric test and analysis of hydration heat of fly ash paste-filling material[J]. Journal of China Coal Society,2015,40(12):2801−2806
[15] 刘鹏亮,张华兴,崔 锋,等. 风积砂似膏体机械化充填保水采煤技术与实践[J]. 煤炭学报,2017,42(1):118−126. LIU Pengliang,ZHANG Huaxing,CUI Feng,et al. Technology and practice of mechanized backfill mining for water protection with aeolian sand paste-like[J]. Journal of China Coal Society,2017,42(1):118−126.
[16] 邵小平,陶叶青,刘二帅,等. 陕北浅埋煤层似膏体充填条带开采参数研究及应用[J]. 煤炭科学技术,2021,49(7):63−70. SHAO Xiaoping,TAO Yeqing,LIU Ershuai,et al. Study and application of paste-like filling mining parameters of shallow buried coal seam in Northern Shaanxi[J]. Coal Science and Technology,2021,49(7):63−70.
[17] 朱建明,马中文,徐金海,等. 小煤窑采空区充填材料及复采技术研究[J]. 金属矿山,2012(3):10−14. doi: 10.3969/j.issn.1001-1250.2012.03.003 ZHU Jianming,MA Zhongwen,XU Jinhai,et al. Study on goaf filling material of small coal mine and repeated mining[J]. Metal Mine,2012(3):10−14. doi: 10.3969/j.issn.1001-1250.2012.03.003
[18] 周保精,徐金海,吴 锐,等. 特厚煤层小窑采空区充填复采技术研究与应用[J]. 采矿与安全工程学报,2012,29(3):317−321. ZHOU Baojing,XU Jinhai,WU Rui,et al. Research on filling repeated mining technology in small pit goaf in extremely thick coal seam and its application[J]. Journal of Mining and Safety Engineering,2012,29(3):317−321.
[19] 李永刚. 综放工作面过小窑采空区注浆充填技术应用[J]. 河北煤炭,2009,(5):3−4,17. LI Yonggang. Application of the grouting filling technology in goaf on fully mechanized top coal cavingmining face through the large coalpit[J]. Hebei Coal,2009,(5):3−4,17.
[20] 白锦文,崔博强,戚庭野,等. 关键柱柱旁充填岩层控制基础理论[J]. 煤炭学报,2021,46(2):424−438. BAI Jinwen,CUI Boqiang,QI Tingye,et al. Fundamental theory for rock strata control of key pillar-side backfilling[J]. Journal of China Coal Society,2021,46(2):424−438.
[21] 彭广虎,祁传西,赵孔友. 基于CFB发电技术的煤泥及矸石利用途径探讨[J]. 煤炭工程,2017,49(S1):56−60,66. PENG Guanghu,QI Chuanxi,ZHAO Kongyou. Discussion on utilization approach of coal slime and gangue based on CFB power technology[J]. Coal Engineering,2017,49(S1):56−60,66.
[22] 任瑞晨,胡秀明,李彩霞,等. 动力煤煤泥处理及利用现状[J]. 贵州大学学报(自然科学版),2021,38(5):1−6. REN Ruichen,HU Xiuming,LI Caixia,et al. Treatment and utilization of steam coal slime[J]. Journal of Guizhou University(Natural Sciences),2021,38(5):1−6.
[23] 王瑞斌. 煤泥燃烧资源化利用现状及展望[J]. 洁净煤技术,2021,27(S2):23−26. WANG Ruibin. Present situation and prospect of coal slurry combustion as a resource[J]. Clean Coal Technology,2021,27(S2):23−26.
[24] 王栋达. 煤泥基预充填胶结材料力学特性研究及应用[D]. 湘潭:湖南科技大学,2022. WANG Dongda. Mechanical properties of coal slime based pre filling cementitious material and its application[D]. Xiangtan:Hunan University of Mining and Technology,2022.
[25] 孙恒虎,刘文永. 高水固结充填采矿[M]. 北京:机械工业出版社,1998. SUN Henghu,LIU Wenyong. High-water cemented material backfilling mining[M]. Beijing:China Machine Press,1998.
-
期刊类型引用(2)
1. 毕银丽,武超,彭苏萍,田乐煊,张延旭. 西部煤矿区微生物修复促进植物水分高效利用策略. 煤炭学报. 2024(02): 1003-1010 . 百度学术
2. 毕银丽,柯增鸣,高学江. 露天煤矿区不同生态措施对紫穗槐水分利用特征的影响. 水土保持学报. 2024(02): 190-196 . 百度学术
其他类型引用(1)