Experimental study on confined compression deformation and breakage characteristics for different types of coal particles
在煤矿采掘活动中，煤体在不同的应力环境下产生破裂破碎和重新压实，而构造煤本身就是在地质条件下形成的煤颗粒的集合体，强度低、黏结性弱［1-4］，因此不同类型煤颗粒体在加载条件下的破碎特性、粒径迁移以及孔隙率演化等对研究煤与瓦斯突出［5-7］和瓦斯抽采［8-9］等具有重要的研究意义。
ZHONG 等［10］通过单个煤颗粒的单轴压缩试验，指出粒径与煤颗粒的强度、煤颗粒的杨氏模量均呈负相关。 DONG 等［11］对近似球形的单个煤颗粒进行大量的单轴压缩试验，提出煤颗粒在单轴加载过程依次表现出弹性与弹塑性，煤颗粒的抗压强度、当量杨氏模量都与粒径成负幂指数关系，且构造煤与原生煤在强度上存在较大差异。 WANG 等［12］也对大量的单个煤颗粒进行单轴压缩试验，利用试验数据分析对比各传统的强度理论，指出其中一种较少被关注的理论较为适用于煤颗粒。 由于单个煤颗粒无法反映出煤体的裂隙，其并不能完全表征煤的力学性质，而构造煤多为弱黏结甚至无黏结，因此可以用煤颗粒群进行替代，进一步反映构造煤的性质。GUO 等［13］将构造煤视为无黏结颗粒，系统地研究了煤颗粒双轴剪切特性，结果表明构造煤颗粒只能存储少量的弹性势能，大量外界做功被耗散，通过研究煤颗粒群揭示了构造煤部分特性。 马占国等［14］对饱和煤岩颗粒进行了侧限压缩试验，其指出饱和煤岩颗粒压缩过程的应力与应变成指数关系，且粒径对应变速率具有较大的影响，以及含水率对应变影响较大，但是没有深入研究煤颗粒的侧限加载过程。
通过对构造煤和原生煤的煤颗粒进行侧限压缩试验，通过研究加载过程的应力应变曲线和切线模量，分析煤颗粒侧限压缩过程中的变形特性。 采用Hardin 的相对破碎率模型［15］和染色标定法［16］，分析煤颗粒群破碎特征并探讨构造煤与原生煤的煤颗粒之间的差异。
试验所采用的煤样分别来自平煤十三矿、任楼矿、寺河矿和红沙梁矿，其中平煤十三矿和任楼矿煤样为构造煤样，寺河矿和红沙梁矿均为原生结构煤煤样。 将4 个煤样破碎并筛分出大颗粒（粒径为3～4 mm）、中颗粒（粒径为2 ～3 mm）和小颗粒（粒径为1～2 mm），按照表1 对4 种煤样都配置6 种不同级配的煤颗粒群试样，6 种级配试样的初始级配曲线如图1 所示。
表1 试样组分
Table 1 Sample components
级配编号各粒径煤颗粒质量／g 1～2 mm2～3 mm3～4 mm试样总质量／g S1480048 S2048048 S3004848 M116161648 M220161248 M32416848
图1 初始粒径级配
Fig.1 Initial particle size grading
试验采用图2所示的YAW4206 微机控制电液伺服压力试验机进行试样加载并记录其应力应变曲线。 试验时将煤颗粒放入图3 中内径50 mm、高100 mm的缸筒里，并将压头置于煤颗粒之上。 为避免加载过程中细小煤颗粒漏出所造成的误差，在缸筒底部进行密封处理。 试验时控制下压板以0.1 mm／s 的速率向上移动，记录应力应变曲线，各试样所加载的最大应力见表2。
图2 YAW4206 试验机
Fig.2 YAW4206 Experimental machine
图3 侧限压缩装置
Fig.3 Confined compression device
表2 试验方案
Table 2 Test scheme
试验压力／MPa级配任楼矿平煤十三矿红沙梁矿寺河矿S15、10、15、205、10、15、201010 S210101010 S35、10、15、205、10、15、201010 M11010（染色）10（染色）10（染色）M21010（染色）10（染色）10（染色）M31010（染色）10（染色）10（染色）
试样加载完后对其进行筛分，记录各粒径区间上煤颗粒的质量，并绘制其级配曲线。 对染色试样组而言，还需要对筛分后的煤颗粒进行颜色识别。
试验中的侧限压缩装置为刚度远大于煤颗粒的钢制品，可近似认为加载过程中该装置无形变，只认为试样发生轴向的形变。 从图4 所示的煤颗粒加载应力应变曲线可以看出，初始时刻曲线存在一定幅度且较为频繁的波动，随后曲线较为平缓斜率小，试验后期曲线较陡斜率较高，因而可将试样加载过程分为3 个阶段：滑移阶段、破碎阶段、压固阶段。
图4 煤颗粒加载应力应变曲线
Fig.4 Stress-strain curve of coal particles under loading
煤颗粒加载初期为滑移阶段，曲线呈现明显的波动性，此阶段松散的煤颗粒在轴向应力作用下发生相对滑移填充孔隙，伴随少量煤颗粒破碎，整体呈现抗压能力较弱、易变形特性，该过程中外界输入能量主要用于颗粒摩擦耗散；随后进入破碎阶段，该阶段煤颗粒更为密集，难以发生相对滑移，颗粒间接触力随加载增强，破碎能阈值较低的颗粒首先破碎，破碎后的小煤颗粒再次滑移填充孔隙，颗粒间接触力重新分配，使得试样更加密实，因而宏观表现为试样变形量较大，该过程中外界输入能量主要被颗粒破碎和摩擦消耗［17］；压固阶段出现在加载后期，煤颗粒经过破碎滑移后结构更为稳定，宏观上表现为应变硬化，此阶段外界能量较多地存储在煤颗粒中［18］。
对比构造煤与原生煤试样的应力应变曲线可以发现，破碎阶段构造煤曲线的斜率明显低于原生煤，表明该阶段构造煤抗变形能力更弱。 在图5 所示的煤颗粒加载后的状况中，图5a 和图5b 等构造煤颗粒在应力10 MPa 下已经出现了明显黏结现象，而图5c 和图5d 等原生煤颗粒十分松散且没有黏结成块现象，表明构造煤颗粒在应力作用后被快速压实，呈现塑性变形，具有流变特性，存在结构咬合［19］甚至黏聚力［20］，且并不能存储较多的能量。
图5 加载后的煤颗粒
Fig.5 Loaded coal particles
试验结束后卸压瞬间，原煤试样出现明显回弹，但构造煤并没有观测到相关现象，表明在压固阶段，构造煤试样外界输入能量主要被塑性变形和破碎消耗，而原煤试样中能量存储为弹性势能在卸荷瞬间释放。
切线模量是应力相对于应变变化率，能反映物体瞬时的变形能力［21］，其计算公式为式（1）：
式中：Et 为切线模量；σ 为轴向应力；ε 为试样轴向应变；F 为轴向载荷；s 为压头轴向位移；l 为试样初始高度。
由于加载过程中煤颗粒滑移与破碎现象以及试验机非完全连续加载均会导致试样加载曲线出现波动，并对计算切线模量产生一定影响，因此对数据进行等间隔取样处理，图6 为处理过后S1 级配试样的切线模量随应力的变化情况。 在滑移阶段和破碎阶段的前期（应力低于2 MPa），原生煤试样的切向模量均高于构造煤试样，进一步印证了构造煤颗粒本身的抵抗变形能力较弱。 破碎阶段末端构造煤与原生煤试样的切线模量的差距逐渐减少，进入压实阶段后等于甚至超过原生煤样，说明构造煤试样经过破碎阶段后更早进入应变硬化阶段，硬化特性更明显。
图6 切线模量-应力曲线
Fig.6 Tangential modulus-stress curve
图7 为2 个不同矿的构造煤的S1、S3 级配试样经不同应力加载后的颗粒级配曲线。 显示了应力的增加，试样的级配曲线不断向左上角移动，趋近其最终级配曲线［22］。 试样级配在应力5 MPa 作用下发生了明显变化，主要对应颗粒滑移阶段和破碎阶段，但增加应力后，级配曲线变化减缓，该过程对应固化阶段，颗粒破碎量相对较少，因此级配变化幅度降低。
图7 试验前后级配曲线
Fig.7 Gradation curves before and after experiment
平煤十三矿S1、S3 和任楼S3 三个试样在应力为15 MPa 和20 MPa 时，级配曲线接近重合，表明3种试样在15 MPa 时已形成稳定级配和力学结构且具有很强的抗变形能力，因此只有极少颗粒出现破碎。
此外，试验发现除S1 级配其余级配试样的级配曲线在加载后均出现了“阶梯”形状，如图7 中S3试样加载后级配曲线所示。 显示了破碎后粒径2.5～3.0 mm 的煤颗粒极少，出现“粒径丢失”的情况，表明侧限加载过程中粒径2.5 ～3.0 mm 的煤颗粒不稳定容易破碎。
颗粒群破碎与初始空隙率、颗粒形貌、粒径级配及颗粒本身力学特性密切关联，是一个复杂的力学过程。 颗粒群的级配曲线变化能较好地表征其破碎情况，许多学者通过级配曲线积分量化得到破碎程度。 采用了Hardin 的相对破碎率模型［15］进行试样破碎的量化分析。 图8 为相对破碎率计算示意，称破碎量与残余破碎势之和为整体破碎势Bp ，则关于相对破碎率计算公式如式（2）—式（4）：
式中：Br 为相对破碎率，无量纲；Bt 为破碎量，为图8 中破碎前后粒径级配曲线和横轴所围成的面积，无量纲；Bp 为整体破碎势，为图8 中初始级配曲线与竖轴所围成的面积，无量纲；d 为粒径，Hardin 认为粒径小于0.074 mm 的颗粒非常稳定，难以发生破碎［15］；bp 为对粒径进行无量纲化处理后的值；p 为小于某粒径的质量占比。
图8 相对破碎率计算示意
Fig.8 Diagram of relative breakage rate calculation
通过计算试样相对破碎率可以得到图9 所示的各因素与相对破碎率关系，图9a 显示了应力对相对破碎率有较大影响，但是当应力较高时，影响程度较弱，这也进一步印证了煤颗粒在压固阶段发生的破碎较少。
从图9b 中可以看出，单一级配系列（S1、S2、S3）的试样相对破碎率近似随着粒径的增加而增加。 这是因为，粒径较大的颗粒自然堆积后容易产生较大的孔隙，且难以通过颗粒滑移填充，故易产生较强的接触力，更容易达到颗粒破碎阈值导致颗粒破碎，破碎后的颗粒填充孔隙后形成抗荷载能力更强的结构；反之对于粒径较小的颗粒，自然堆积后产生孔隙较小，产生少量破碎后即可填充孔隙形成稳定抗压结构，因此，单一级配系列的试样相对破碎率近似随粒径增大而增大。
图9 各因素对试样相对破碎率的影响
Fig.9 Influence of various factors on relative breakage rate of samples
由于混合级配系列（M1、M2、M3）的试样中煤颗粒粒径分布相对更均匀，初始孔隙较低，故混合级配系列试样相对破碎率总体较单一级配系列更低。煤颗粒的相对破碎率受到粒径级配和煤种共同影响。 混合级配系列中红沙梁矿试样的相对破碎率受级配影响极小；任楼矿和寺河矿煤颗粒破碎率随大粒径颗粒减少而增大；平煤十三矿的单一级配和混合级配系列试样均未表现出明显规律，这与构造煤的塑性、弱黏结和力学性质不均匀有关。
粒径3～4 mm 为大粒径；2～3 mm 为中粒径；1～2 mm 为小粒径。 试验后染色试样组的煤颗粒在一定粒径区间内尺寸差距较小，可认为均匀混合的各染色煤颗粒质量百分比与其在一个面上的面积百分比相等［23］。 因此可通过图10 所示的颜色识别步骤可以得到中粒径上绿色与红色煤颗粒的质量。
图10 染色煤颗粒颜色识别过程
Fig.10 Color recognition process of dyed coal particles
式中：pb-m 为大粒径煤颗粒破碎至中粒径的占比；mb-m 为大粒径煤颗粒破碎至中粒径的质量；ms，b 为大粒径煤颗粒存活的质量；m0，b 为其初始质量。
对颜色识别后的数据采用式（5）和式（6）处理可得到图10 所示的大、中粒径煤颗粒破碎特征。 中粒径煤颗粒的绝对破碎率均在50%以上且均高于大粒径破碎率，与前文中粒径2.5 ～3.0 mm 颗粒缺失对应，可能是因为煤颗粒结构纹理特性导致了该粒径区间煤颗粒不稳定，极易破碎，因而2～3 mm 煤颗粒破碎率更高。
图11a 中相同条件下构造煤大、中粒径绝对破碎率均明显高于原生煤，也进一步印证了图9b 中相同条件下构造煤相对破碎率均高于原生煤。 此外，图11b 中原生煤试样大粒径煤颗约50%以上破碎至中粒径范围，而构造煤试样的大粒径煤颗粒发生破碎后在中粒径范围内比例更小。 这是因为破碎为中粒径颗粒强度不足以支撑其载荷，因此更多大粒径构造煤破碎为小粒径甚至更小的颗粒。 综合以上结果，表明试验条件相同，构造煤更容易发生破碎且破碎程度更高。
图11 大、中粒径煤颗粒的破碎特征
Fig.11 Breakage characteristics of large and medium size coal particles
在应力加载下，煤颗粒逐渐压实并产生破碎，尤其是构造煤颗粒，本身在经历了地质时期应力破坏破碎后，呈现塑性特征。 而在此次试验中，再次加载破碎，使得构造煤颗粒更容易破碎至更小粒径，甚至呈现流变特征，使得加载后煤体孔隙率更低，如图12 所示。 但外界加载做功主要消耗在煤颗粒摩擦、破碎流变等阶段，构造煤颗粒本身存储的弹性能很低。 在现场条件下，受采掘影响，构造煤也存在二次破碎重新压实作用，改造煤颗粒孔隙结构，导致构造煤区域渗透率更低，具备了保存吸附瓦斯条件［24］。在卸压边界条件下，复杂构造应力叠加采动应力，使得构造煤体产生局部大变形，为转化富集游离瓦斯能提供条件，同时也为破碎抛出弱黏结的构造煤提供能量和物质基础［25-26］。 而原煤颗粒表现出较难破碎，且具有较高的弹性能存储特性。 进一步对不同类型煤颗粒的细观力学、颗粒破碎动力学的深入研究为煤岩动力灾害预测预报提供基础。
图12 煤颗粒加载前后变形破碎结构示意
Fig.12 Schematic diagram of coal particles
1）煤颗粒侧限加载过程可按照应力应变曲线特性分为滑移阶段、破碎阶段和压固阶段。 滑移阶段原生煤和构造煤应力应变曲线特征相同，但构造煤抗变形能力弱，破碎阶段原生煤曲线斜率高于构造煤，构造煤抗变形能力仍弱于原生煤，相比原生煤变形更大，但压固阶段其硬化趋势明显强于原煤，抗变形能力显著增强。
2）煤颗粒粒径级配曲线大幅变化主要集中在破碎阶段，后续继续增大压力后粒径级配曲线变化不再显著，破碎率变化也较小。 煤颗粒破碎后出现中部分粒径区间颗粒极不稳定，易破碎，出现粒径缺失现象。 单一粒径级配试样大颗粒含量高对应破碎率更高，均匀粒径级配较单一级配整体破碎率更低。
3）在相同级配、载荷条件下，构造煤大粒径、中粒径破碎率均明显高于原生煤样，整体破碎率也高于原生煤，破碎过程中，大粒径原生煤更多破碎为中粒径，而构造煤大粒径则更多破碎为小粒径及以下的煤颗粒。
［1］ 王恩营，刘明举，魏建平.构造煤成因-结构-构造分类新方案［J］.煤炭学报，2009，34（5）：656-660.WANG Enying，LIU Mingju，WEI Jianping.New genetic texturestructure classification system of tectonic coal［J］. Journal of China Coal Society，2009，34（5）：656-660.
［2］ 李 明.构造煤结构演化及成因机制［D］.徐州：中国矿业大学，2013.LI Ming.Structure evolution and deformation mechanism of tectoni⁃cally deformed coal［D］.Xuzhou： China University of Mining and Technology，2013.
［3］ 琚宜文， 姜 波， 侯泉林，等.构造煤结构-成因新分类及其地质意义［J］.煤炭学报， 2004，29（5）：513-517.JU Yiwen， JIANG Bo， HOU Quanlin，et al.The new structure-ge⁃netic classification system in tectonically deformed coals and its ge⁃ological significance［J］.Coal Science and Technology， 2004，29（5）：513-517.
［4］ 张玉贵，张子敏，曹运兴. 构造煤结构与瓦斯突出［J］.煤炭学报，2007，32（3）：281-284.ZHANG Yugui， ZHANG Zimin， CAO Yunxing. Deformed coal structure and control to coal-gas outburst［J］.Journal of China Coal Society，2007，32（3）：281-284.
［5］ 张友谊，黄战峰，余贵军，等.地质条件影响下的煤与瓦斯突出特征及防治对策［J］.煤炭科学技术，2013，41（12）： 39-42，45.ZHANG Youyi，HUANG Zhanfeng，YU Guijun，et al.Coal and gas outbursts characteristics and countermeasures under the influence of geological conditions［J］.Coal Science and Technology， 2013，41（12）：39-42，45.
［6］ GRAY I.Outburst risk determination and associated fac-tors［M］.Brisbane： Australian Coal Research Ltd，2015.
［7］ 程远平，雷 杨.构造煤和煤与瓦斯突出关系的研究［J］.煤炭学报，2021，46（1）：180-198.CHENG Yuanping，LEI Yang.Causality between tectonic coal and coal and gas outbursts［J］.Journal of China Coal Society，2021，46（1） ： 180-198.
［8］ LIU Wei，WU Deyao，XU Hao，et al.A permeability evolution model of coal particle from the perspective of adsorption deforma⁃tion［J］. Energy Science ＆ Engineering，2021，9 （4）：577-787.
［9］ 董 骏. 基于等效物理结构的煤体瓦斯扩散特性及应用［D］.徐州：中国矿业大学，2018.DONG Jun. Gas diffusion properties of coal mass based on equivalent physical structure［D］.Xuzhou： China University of Mining and Technology，2018.
［10］ SHAN Zhong，FELIX BAITALOW，PETR NIKRITYUK，et al.The effect of particle size on the strength parameters of German brown coal and its chars［J］.Fuel，2014，125：200-205.
［11］ DONG Jun，CHENG Yuanping，HU Biao，et al.Experimen tal study of the mechanical properties of intact and tectonic coal via compression of a single particle［J］.Powder Technology，2018，325：412-419.
［12］ WANG Chenghao，CHENG Yuanping，HE Xinxin， et al. Size effect on uniaxial compressive strength of single coal particle un⁃der different failure conditions［J］.Powder Technology，2019，345： 169-181.
［13］ GUO Jingjie，LI Wei，JIAO Yang，et al.A discrete numerical simu⁃lation of micromechanics for dense coal granular systems： Impli⁃cation for coal and gas outbursts［J］.Powder Technology，2021，392：448-458.
［14］ 马占国，郭广礼，陈荣华，等.饱和破碎岩石压实变形特性的试验研究［J］.岩石力学与工程学报，2005，24（7）： 1139-1144.MA Zhanguo， GUO Guangli， CHEN Ronghua， et al. An experimental study on the compaction of water-saturated overbroken rock［J］.Chinese Journal of Rock Mechanics and Engi⁃neering，2005，24（7）：1139-1144.
［15］ HARDIN B O. Crushing of soil particles［J］.Journal of Geotech⁃nical Engineering，1985，111（10）：1177-1192.
［16］ 于际都，沈超敏，刘斯宏.染色石膏颗粒一维压缩破碎与粒径迁移［J］ .岩石力学与工程学报，2020，39（5）：1071-1079.YU Jidu，SHEN Chaomin，LIU Sihong. Breakage and migration of dyed gypsum particles under one-dimensional compression［J］.Chinese Journal of Rock Mechanics and Engineering， 2020，39（5）：1071-1079.
［17］ 李 振. 老空区破碎煤岩体变形与渗流特性研究及在煤层气抽采中的应用［D］.太原：太原理工大学，2018.LI Zhen. Research on the deformation and seepage charac -teristics of crushed coal and rock and its application in abandoned coal mine gob methane extraction ［ D］. Taiyuan： Taiyuan University of Technology，2018.
［18］ LI Meng，ZHANG Jixiong，ZHOU Nan，et al. Effect of par-ticle size on the energy evolution of crushed waste rock in coal mines［J］. Rock Mechanics and Rock Engineering，2017，50（5）：1347-1354.
［19］ 李 震，王宏强，高雨航，等.沙柳生物质颗粒致密成型特性的离散元仿真［J］.锻压技术，2020，45（3）：152-158.LI Zhen， WANG Hongqiang， GAO Yuhang， et al. Discrete element simulation on dense forming characteristics for Salix bio⁃mass particles ［J］. Forging ＆ Stamping Technology， 2020，45（3）：152-158.
［20］ 张 帅. 构造煤型煤物理力学性质影响因素分析及制作方法试验研究［D］.徐州：中国矿业大学，2020.ZHANG Shuai.Analysis of factors influencing physical and me⁃chanical properties of tectonically deformed coal｀s moulded coals and experimental research on manufacturing methods ［ D］.Xuzhou： China University of Mining and Technology，2020.
［21］ 冯梅梅，吴疆宇，陈占清，等.连续级配饱和破碎岩石压实特性试验研究［J］.煤炭学报，2016，41（9）：2195-2202.FENG Meimei，WU Jiangyu，CHEN Zhanqing，et al. Experimental study on the compaction of saturated broken rock of con⁃tinuous gradation ［J］. Coal Science and Technology， 2016，41（9）：2195-2202.
［22］ EINAV，I.Breakage mechanics：part Ⅰ： Theory［J］. Journal of the Mechanics and Physics of Solids，2006，55（6）：1274-1297.
［23］ 彭 宇，丁选明，肖 杨，等.基于染色标定与图像颗粒分割的钙质砂颗粒破碎特性研究［J］.岩土力学，2019，40（7）：2663-2672.PENG Yu， DING Xuanming， XIAO Yang，et al.Study of particle breakage behaviour of calcareous sand by dyeing tracking and particle image segmentation method ［ J ］. Rock and Soil Mechanics，2019，40（7）：2663-2672.
［24］ 胡 彪，程远平，王 亮.原生结构煤与构造煤孔隙结构与瓦斯扩散特性研究［J］.煤炭科学技术，2018，46（3）：103-107，24.HU Biao，CHENG Yuanping，WANG Liang. Study on porous structure and gas diffusion characteristics of primary structure coal and tectonic coal ［J］. Coal Science and Technology， 2018，46（3）：103-107，24.
［25］ 涂庆毅.构造煤表观物理结构及煤与瓦斯突出层裂发展机制研究［D］.徐州： 中国矿业大学，2019.TU Qingyi.Study on apparent physical structure of tectonic coal and spallation development mechanism of coal and gas outburst［D］.Xuzhou： China University of Mining and Technology，2019.
［26］ 刘彦伟，刘明举.粒度对软硬煤粒瓦斯解吸扩散差异性的影响［J］.煤炭学报，2015，40（3）：579-587.LIU Yanwei，LIU Mingju. Effect of particle size on difference of gas desorption and diffusion between soft coal and hard coal［J］.Journal of China Coal Society，2015，40（3）： 579-587.