冲击倾向煤样在不同应变率载荷作用下破裂的声发射及碎片分形特征

吕鹏飞，路康斌，曹树斌，丛日盛

（1.内蒙古科技大学 矿业与煤炭学院，内蒙古 包头 014010；
2.国家能源集团 神东煤炭经销中心，陕西 神木 719315；
3.辽宁省化工地质勘查院有限责任公司，辽宁 锦州 121007）

声发射是指煤、岩石类材料体在外载荷作用下产生微裂纹而导致材料体以弹性波形式释放应变能的现象[1-2]。声发射通过对岩石破坏过程的动态连续监测，可细观分析受载岩石内部裂纹生成、发展、动态演化及损伤累积过程，揭示地下工程中的煤岩体破裂特征，对煤矿冲击地压、巷道底鼓、冒顶、金属矿岩爆等灾害的研究具有重要的理论指导意义。

许多学者在岩石类材料受载破裂声发射特征方面进行了研究。曹安业等[3]对砂岩、泥岩在不同加载速率下的声发射特征进行了研究，认为随着加载应变率增大，声发射能量释放形式由孤震型向震群型转变。张黎明等[4]基于单轴压缩试验分析了大理岩受载破坏过程的声发射参量变化特征，发现声发射频率随加载速率增减呈明显的波动，大理岩峰值强度前频率突增，峰值强度后频率突降。张国凯等[5-7]利用声波、声发射一体化装置研究了花岗岩在单轴加载下的应力波传播、裂纹扩展与声发射演化之间的关系，认为初始加载阶段波速与声发射幅值增长趋势一致，但增速较慢，声发射能量在花岗岩宏观破裂前发生突发性增长，并将此作为煤岩体灾害性破坏预判的前兆信息。陈宇龙等[8-9]基于力学实验获取了岩石应力-应变曲线4 个阶段的声发射特征，认为在初始压密和弹性阶段，声发射撞击数少、能量低、振幅小，较少破裂事件产生；
在应变硬化阶段，撞击数骤增、能量高、振幅大，大量破裂事件产生；
在应变软化阶段撞击数骤减、能量低、振幅小，少量破裂事件产生。姚强岭等[10-11]研究了含水砂岩受载的声发射特征，在干燥状态下，含水砂岩声发射计数峰值与应力峰值位置基本一致；
随着含水率增加，声发射计数峰值滞后于应力峰值，且含水率越高，滞后现象越明显。

岩石类材料的声发射特征可清晰地反馈岩石破裂内部细观结构演化及损伤情况，而岩石类材料受载发生宏观破裂程度的评价及相关性分析亦是地下工程研究中的重点问题。分形理论[12]作为描述无规则事件的前沿理论，多次被应用于定量描述煤岩破坏程度，分析煤岩破坏与其他参数的关系。刘少虹等[13-15]采用一维过分离式霍普金森压杆（Split Hopkinson Pressure Bar，SHPB）动静组合实验分析了试件破坏碎片分维与应力波关系，认为组合煤岩试件碎片分维随应力波能量的增大而增大，随静态载荷的增大呈先增大后减小的规律。郑克洪[16]研究了煤矸石破碎粒度的分形表达式，并与威布尔分布结果进行了对比，认为采用分形分布表示煤矸石破碎粒度能更好地表达煤矸石破碎粒度的分布规律。纪杰杰等[17]探究了冲击荷载作用下岩石破碎分形特征，认为岩石动态抗压强度随应变率增大而增大，二者满足乘幂函数关系。

上述研究对岩石受载破坏过程的声发射特征和试件碎片分形特征进行了相关分析，并取得了一定成果，但针对冲击倾向煤样在不同应变率单向受载条件下破坏程度的定量描述及与加载应变率的定量关系的研究较少。本文通过冲击倾向煤样在单轴静态载荷基础上施加不同应变率动态载荷加载实验，并采用分形理论分析试件碎片与载荷应变率的定量关系，以及载荷应变率与试件破坏模式的对应关系等，来研究冲击倾向煤样在不同应变率载荷作用下破裂的声发射及碎片分形特征。

不同岩石工程中对载荷状态的界定不同，应变率是界定载荷状态的唯一指标[3]。在煤矿动静态载荷应变率的界定中，当应变率大于10-3s-1时属于动态载荷[18-19]，大于10-5s-1小于10-3s-1时属于准动态载荷，小于10-5s-1时属于静态载荷，见表1。

表1 煤矿载荷的应变率界定标准Table 1 Strain rate definition criterion of loadings in coal mine

2.1 试验系统

试验系统如图1 所示。选用MTS-C64.106 电液伺服系统作为压力加载系统，额定承载能力为1 000 kN，动作器位移速度为0.5～90 mm/min，位移分辨率为0.2 μm。声发射系统选用集成PCI-2 卡型声发射监测设备，可同时实现特征参数提取和波形处理，频率范围为103～3×107HZ，采样率为4×107Hz。试验中，声发射门槛值设定为40 dB，受载试件采用6 个传感器控制，距离试件顶底端10 mm 分别布置3 个传感器且呈180°布控。

图1 试验系统Fig.1 Experimental system

2.2 试件加工

煤样采集于兖州煤业股份有限公司东滩煤矿63上05 工作面。选取具有弱冲击危险煤层中无明显裂隙发育的大块原煤，钻取直径为50 mm 的圆柱体，截锯成高度为100 mm 的试件，将两端打磨平整，确保两端不平整度范围为±0.5 mm，断面与轴线不垂直度范围为±0.15°。共加工24 个试件。试件参数见表2，试件成品如图2 所示。

图2 试件成品Fig.2 Test samples

表2 试件基本参数Table 2 Basic parameters of samples

2.3 试验方案

为研究冲击倾向煤样在不同应变率载荷作用下的破坏程度，先采用静态载荷加载试件，当加载到5 kN 后保持20 s，再采用不同应变率载荷控制试件，直到试件破坏。载荷应变率分别为8.33×10-6，4.17×10-5，8.34×10-5，1.67×10-4，3.34×10-4，1.00×10-3，2.50×10-3，5.00×10-3s-1，8 次为一组试验，共进行3 组试验，试验结果取平均值。

2.4 试验结果

不同应变率载荷加载作用下的应力-应变曲线如图3 所示。可看出在静态载荷基础上，试件破坏的峰值强度随着应变率载荷的升高而逐渐上升，说明高应变率载荷触发试件破坏的峰值强度更大。在冲击地压煤矿生产中，高应变率触发的冲击地压强度更大。

图3 不同应变率载荷加载作用下的应力-应变曲线Fig.3 Stress-strain curves under loading modes with different strain rates

3.1 声发射振铃计数和能量特征分析

选取应变率为8.33×10-6，1.67×10-4，2.50×10-3s-13 种条件下监测到的声发射振铃计数、能量幅值做图，如图4 所示。可看出声发射振铃计数和能量幅值都随时间增加而增大，且二者的变化趋势一致；
当应变率增大，加载时间减少，振铃计数减少，但能量幅值升高，说明应变率增大有利于高能震动产生；
随着应变率增大，振铃计数与能量幅值的增长模式由缓慢增长逐渐发展为急速增长，最终发展为突然增长，说明煤矿中的高应变率载荷更易触发高能动力事件。

图4 不同应变率载荷加载下声发射振铃计数、能量特征Fig.4 Ringing count and energy characteristics of AE under different strain rate loading

3.2 能量输入与声发射参数关系分析

煤体声发射现象及内部损伤受外界能量影响，煤体在加载之后产生裂纹，同时释放能量，并以应力波的形式传递，在传递到试件表面时可被声发射探头接收，每个探头可监测能量输入速率、振铃计数、撞击事件数[13-14]。因此通过上述3 种监测值对比分析能量输入与试件内部变化的关系，如图5 所示。一维情况下，能量输入为

式中：Q为总能量输入，J；
V为体积，cm3；
σ为应力，Pa；
ε为应变。

由图5 可看出：加载应变率较低时能量输入较均衡，应变率较高时能量输入不稳定，突增突降现象明显；
随着应变率增大，能量输入经历缓增-急增-突增的转变，这与声发射监测到的能量特征一致。能量输入速率与振铃计数、撞击事件数增长趋势基本相同，尤其表现在能量输入速率达到最大时，在微小范围内振铃计数与撞击事件数也达到最大，此时试件内部裂纹扩展和损伤程度最大，因此认为能量输入速率与试件内部损伤呈正相关；
试件达到峰值载荷时，裂纹扩展主要表现为裂隙的贯通和宏观破坏，内部损伤急速增加；
在能量输入速率峰值之后的短时间内，振铃计数与撞击事件数迅速降低，说明试件正在完全失去承载能力。对比不同应变率加载的振铃计数与撞击事件数发现，应变率越大，单位时间内振铃计数与撞击事件数越大，说明应变率增大促进了试件内部裂隙的快速扩展与撞击事件的形成。

图5 能量输入与声发射参数关系Fig.5 Relationship between AE parameter and energy input

3.3 声发射相关常数分析

为深入揭示试件受载内部损伤及破裂程度，采用古登堡公式处理声发射数据。

式中：N为震动频次；
a为与震动强弱正相关的常数；
b为与高低能震动数之比负相关的常数；
E为震级。

选取试验中应变率为8.33×10-6，1.67×10-4，2.50×10-3s-13 种条件下接收的20 组lgE与lgN的声发射数据进行处理，求解a，b，结果见表3。可看出a随应变率增大而减小，b随应变率增大而增大。这是因为应变率增大，加载时间缩短，监测到的声发射频次减少，试件裂纹扩展活动性减弱，但高能声发射事件增加，高低能声发射比例增大，所以b增大，试件内部损伤加剧，煤矿生产中更易触发高能动力事件。

表3 不同应变率加载下a，bTable 3 The a and b values under different strain rate loading

3.4 试件动态破坏模式与损伤特性

采用高速摄像机监拍不同冲击载荷下试件动态破坏过程，如图6 所示。可看出加载应变率越大，试件破坏越严重；
随着加载应变率增大，试件破坏形式由剪切破坏向劈裂破坏再向爆裂破坏发生转变；
当加载应变率较低时，主要是试件上半部分破坏；
应变率增大后由试件中部逐渐向下半部分延伸破坏；
整个破坏过程几乎不发生明显的黏性变形，而以脆性破坏为主，所以由动态载荷引发的煤体裂纹扩展主要是脆性扩展行为。

图6 试件动态破坏过程Fig.6 Dynamic failure process of samples

4.1 碎片粒径特征

采用声发射可较好地定量分析煤样在静动组合加载作用下的细观破坏特征，采用分形理论来定量描述煤样加载破坏程度。为便于筛分粒径，试验中采用筛分粒径分别为1，2，5，10，20，40 mm 的筛子对冲击碎片进行筛分，并称出对应粒径煤屑的质量，对比分析不同粒径碎片的质量分数，见表4。可看出冲击碎片粒径主要集中在10～40 mm。应变率较低时，碎片粒径集中在10～20 mm 和20～40 mm，且此时一般不产生大粒径（≥40 mm）的碎片。粒径为10～20 mm 且应变率低于1.67×10-4s-1时的碎片随应变率的增大而增多。小粒径碎片增多，则碎片总表面积增大，消耗能量增大，因此这种扰动载荷通常不引起试件很大强度的破坏。当应变率较高时，大粒径碎片随应变率增大先增后减。大粒径碎片增多，则碎片总表面积减少，消耗的能量减少，而应变率升高又会使储存的弹性能增多，试件破坏强度增强。当应变率达到最大时，大粒径碎片减少，因此推测其存在一个产生最大冲击强度的应变率极值，需要进一步探索。

表4 不同加载条件下冲击碎片质量百分比Table 4 Quality percentage of coal fragment under different loadings

4.2 碎片分维与冲击载荷应变率关系

试验中，试件破碎程度不同的现象在工程中可表射为不同程度的冲击事故。分形理论是一种广泛应用于定量描述无规则事件的理论，采用分形理论计算试件碎片分维，定量评价碎片破碎程度具有合理性。

为计算在不同应变率载荷加载条件下冲击碎片分维，可对ln（Mmax/M）与lnU进行拟合，拟合方程为

式中：d为冲击碎片质量分维；
Mmax为最大块度的碎片质量，g；
M为试件冲击碎片质量，g；
U为质量不小于M的碎片数量。

当M分别为1，2，3，5，10，15，20 g 时，计算冲击碎片分维，结果见表5。可看出随着应变率增大，冲击碎片分维先增大后减少，相关性系数越大，表示拟合效果越好。

表5 不同加载条件下冲击碎片分维计算结果Table 5 Fractal dimension results of rock burst fragments under different loadings

依据表5 数据，对碎片分维进行拟合，拟合后的方程为

不同载荷应变率与碎片分维关系如图7 所示。可看出拟合曲线呈二次函数形状，随着应变率的提高，冲击碎片质量分维先增大后减小，这与岩石材料[20]存在一定区别。从冲击动能角度出发，冲击碎片动能与质量和冲击速度呈正比，但质量增大带来的后果是速度减小，因此推测存在一个极大值，能够使冲击动能达到最大，抛物线极大值对应应变率所产生最大冲击动能的冲击载荷应变率，试验结果显示约为2.8×10-3s-1。

图7 不同载荷应变率与碎片分维关系Fig.7 Relationship between fractal dimension of fragments and different loads strain rate

（1）随着加载应变率增大，声发射总体数量减少，高能声发射事件增多，声发射振铃计数和能量幅值经历缓增-急增-突增的一致转变过程；
试件受载的能量输入速率与声发射振铃计数、内部撞击事件数增长趋势基本一致；
声发射a值随加载应变率增大而减小，b值随加载应变率增大而增大。

（2）随着加载应变率增大，原煤试件的破坏模式会发生剪切破坏-劈裂破坏-爆裂破坏的转变。当加载应变率较低时，主要是试件上半部分破坏，应变率增大后由试件中部逐渐向下半部分延伸破坏，原煤试件在动载应变率作用下的破坏过程主要是裂纹的脆性扩展行为。

（3）原煤试件受载碎片分维与加载应变率呈二次函数关系，即存在应变率极值使碎片破坏程度达到最大，试验结果显示该值为2.8×10-3s-1。