硫含量对煤自燃特性的影响实验与模拟研究

刘军见 ，唐一博 ，刘洪刚

（1.山西地宝煤炭综合检测中心有限公司，山西 太原 030045；2.太原理工大学 安全与应急管理工程学院，山西 太原 030024）

煤炭是我国重要的能源支柱[1]，在能源结构中占主体地位[2-3]。煤炭的开采通常伴随着安全事故的发生，包括瓦斯爆炸、透水、塌陷、顶板、火灾等[4]，其中煤矿火灾是常见事故之一，多发于采空区。2001—2022 年的数据表明[5-9]，煤矿火灾事故总量虽然有着减少的趋势，但煤矿火灾伤亡数和事故数的比例在不断上升。煤矿火灾发生的主要原因是采空区遗煤自燃。煤自燃是煤炭自身产热大于煤炭向外界散发的热量从而导致热量积聚、温度上升，进而自发燃烧的现象。

煤自燃有着众多影响因素，如孔隙、结构、含水量、元素成分等。硫元素的含量是煤炭自燃的影响因素之一。煤中硫的主要存在形式分为无机硫和有机硫。学者们对于无机硫的研究主要围绕煤中黄铁矿的成分开展。DENG 等[10]考察了不同黄铁矿含量的煤样对标志气体、温度和放热等自燃参数的影响，发现黄铁矿可以非线性地加速自燃过程，且黄铁矿含量为5%的煤样CO 释放率和氧吸附量最大，黄铁矿含量为7%的煤样具有最大的热流率，黄铁矿含量为5%～7%对自燃的影响最为显著；YANG 等[11]通过热分析表征了不同黄铁矿含量的煤在4 种升温速率下的反应过程，计算了不同黄铁矿含量煤的活化能，发现煤样的表观活化能随着黄铁矿含量的增加而降低，即煤中黄铁矿可以促进煤的自燃；DING 等[12]选取褐煤、烟煤、无烟煤3 种不同变质程度的低硫煤煤样来研究黄铁矿对煤自燃的影响，结果发现Fe3+、Fe2+和H+对煤的低温氧化都有催化作用，其中Fe3+的催化能力最强，黄铁矿在干湿条件下对煤自燃都有催化作用，且潮湿环境下催化作用更显著；常绪华等[13]通过对10 余种阻化剂的研究和筛选，得出水玻璃，Ca(OH)2、Na2CO3、NaHCO3、CaCl2、MgCl2等阻化剂对新峪矿高硫煤的氧化有明显的抑制作用。上述研究表明煤中无机硫对煤炭自燃有着促进作用。然而除了无机硫，煤中有机硫也是影响煤炭自燃特性的重要因素。煤中有机硫结构主要有硫醇、硫醚、噻吩等。国内外学者围绕有机硫对煤自燃的影响进行了深入研究。胡军等[14]研究了26 个地区煤样中的有机硫，发现煤中有机硫含量多在1%范围内，高硫煤中主要以无机硫为主，低硫煤中主要以有机硫为主；王雪峰等[15]围绕金属离子对煤中含硫活性基团的影响进行研究，发现Ca2+能与煤中含硫结构形成配合物，形成配合物后分子的稳定性增强，化学活性减弱，起到抑制煤自燃的作用；叶正亮等[16]结合宏观与微观分析发现过硫酸钠可以作为一种防止高硫煤自燃的化学阻化剂；张兰君[17]、ZJANG 等[18]通过程序升温、气相色谱、TG/DSC 对有机硫模型化合物的自燃特性进行了实验分析，发现有机硫官能团中氧化活性最高的3 个官能团为硫醇、硫醚、噻吩；ZHENG 等[19]利用模型化合物，选择含有不同有机硫官能团的代表性有机硫模型化合物进行低温氧化实验和热分析，揭示了它们的氧化和热力学特性，获得了低温氧化过程中有机硫与煤中其他活性基团相互作用的机理；GAO 等[20]通过低温氧化实验,研究了变质程度相近、有机硫含量不同的煤的自燃倾向，发现煤中有机硫的形态主要有硫醇类、硫醚类、噻吩类、砜类，经低温氧化后，煤中硫醇、硫醚、噻吩、甲基、亚甲基、吡啶含量下降，而砜、亚砜、硫酸根、羧基、氮氧化物等含氧基团含量增加，模型化合物低温氧化反应后，部分有机硫以芳香族硫或硫氧化物的形式存在，小部分硫以SO2和H2S 气体的形式在固体氧化产物中逸出。上述研究表明，煤中有机硫的氧化会促进煤自燃的过程，通过抑制有机硫活性可以对煤自燃起到抑制作用。

现阶段高硫煤阻化剂的选取依据主要考虑阻化剂与煤中无机硫(FeS2)的相互作用，而阻化剂对高硫煤中有机硫的抑制作用却鲜有研究。为此，使用程序升温实验系统、综合热分析仪对典型高硫煤和典型低硫煤的氧化过程进行实验分析；使用量子化学模拟对高硫煤常用阻化剂中Na+、Mg2+、Ca2+和2 种煤中典型有机硫结构(苯硫醇、二苯硫醚)的络合作用进行计算分析。

1 研究方法

1）煤样选取。从若干煤样中筛选出硫元素占比差异较大且其余元素占比相近的2 种煤样，分别代表典型低硫煤和典型高硫煤，煤样的元素分析见表1。

表1 典型低硫煤的元素参数Table 1 Elemental parameters of typical low-sulfur coal %

2）程序升温实验系统。程序升温实验系统如图1。测试温度为20～200 ℃，升温速率为1 ℃/min，干空气通入速率为10 L/min，利用气相色谱仪检测各温度下CO 与CO2体积分数。

图1 程序升温实验系统Fig.1 Experimental system of temperature programming

3）热重/差示扫描量热测试(TG/DSC)。使用耐驰STA449F5 同步热分析仪，获得样品的热重(TG)和差示扫描量热(DSC)曲线。升温区间为20～800 ℃，升温速率为10 K/min，保护气(N2/O2)和吹扫气(N2/O2)通入速率分别为30 mL/min 和20 mL/min。

4）分子动力学模拟。使用Gaussian09 软件，使用密度泛函理论(DFT)在B3LYP/6-31G(d)水平上对煤中含硫活性基团(苯硫醇、二苯硫醚)和Na+、Mg2+、Ca2+形成的络合物结构进行计算。对络合物的静电势图和HOMO/LUMO(最高占据轨道/最低非占据轨道)进行计算分析。

2 实验结果

2.1 程序升温实验结果

通过2 种煤样的程序升温实验系统测试结果可知，随着温度不断升高，陆续检测出CO2、CO、CH4等气体，其中CO2和CO 在初期就可测得，二者的生成量可以用来进行煤自燃预测。煤样CO2体积分数变化曲线如图2，煤样CO 体积分数变化曲线如图3。

图2 煤样CO2 体积分数变化曲线Fig.2 CO2 concentration variation curves of coal samples

图3 煤样CO 体积分数变化曲线Fig.3 CO concentration variation curves of coal samples

从图2 可以看出：高硫煤和低硫煤的CO2体积分数曲线趋势基本一致。高硫煤在20～50 ℃时，CO2释放量随温度升高的变化较小，煤与氧气反应较为轻微；50 ℃后CO2体积分数随着温度升高显著增大；当煤样温度达到140 ℃后，CO2体积分数达到爆发式增加，总体符合指数关系。低硫煤相比高硫煤，煤样的特征温度向后推迟，160 ℃后，CO2体积分数才开始爆发式增加。

由图3 可以看出：高硫煤和低硫煤的CO 浓度曲线趋势也基本一致。高硫煤在20～110 ℃时，CO 释放量随温度的升高变化较小，煤与氧气反应较为轻微；110 ℃后CO 体积分数随着温度升高显著增加；当煤样温度达到160 ℃后，CO2体积分数达到爆发式增加，总体符合指数关系。低硫煤相比高硫煤，煤样的特征温度向后推迟，170 ℃后，CO2体积分数才开始爆发式增加。

综上，典型高硫煤在低温氧化过程中各温度阶段的CO2、CO 标志气体释放量大于典型低硫煤，说明高硫煤相比低硫煤具有更高的自燃倾向。

2.2 热分析(TG/DSC)结果

对2 种煤样进行热分析（TG/DSC），得到煤样的TG 和DSC 曲线，根据煤样氧化过程的特点定义5 个特征温度点 （T1～T5）。典型高硫煤和典型低硫煤氧化过程的特征温度点见表2。

表2 典型高硫煤和典型低硫煤氧化过程的特征温度点Table 2 Characteristic temperature points of typical high-sulfur coal and typical low-sulfur coal ℃

1）T1干裂温度。煤样由于脱水导致质量达到极小值时对应的温度，此时DTG 数值上为0。

2）T2失重温度。煤样吸氧增重，质量达到极大值时即将受热分解质量下降时对应的温度，此时DTG 数值上为0。

3）T3着火温度。煤样的着火温度点，过T4做TG 曲线的切线，该切线与过T2的水平线交于一点，该点则为T3着火温度。

4）T4最大失重温度。煤样燃烧最剧烈时对应的温度，对应于DTG 曲线的峰值。

5）T5燃尽温度。煤样停止燃烧且质量保持不变时的温度。

典型高硫煤的TG 曲线图如图4，典型低硫煤的TG 曲线图如图5。

图4 典型高硫煤的TG 曲线图Fig.4 TG curves of typical high-sulfur coal

图5 典型低硫煤的TG 曲线图Fig.5 TG curves of typical low-sulfur coal

从图4 可以看出，煤样从室温到800 ℃的氧化过程被划分成脱水失重阶段、吸氧增重阶段、热解失重阶段、燃烧阶段、燃尽阶段5 个阶段；随着温度上升，煤样中的水分蒸发，质量下降，当温度达到189 ℃时，煤样温度达到极小值；随着温度继续上升，煤样吸收空气中氧气进入吸氧增重阶段，质量略微上升；温度达到365 ℃时，煤样进入热解失重阶段，煤中的芳香环等有机结构受热分解，释放出CO2、CO 等气体并放出大量热量；温度达到516 ℃时，煤样达到着火点，标志着煤样进入到氧化最剧烈的燃烧阶段；温度达到593 ℃时，煤样的燃烧最为剧烈；温度达到751 ℃时，煤样质量不再随温度的升高而变化，煤样进入燃尽阶段。

由图5 可以看出：煤样从室温到800 ℃的氧化过程被划分成脱水失重阶段、吸氧增重阶段、热解失重阶段、燃烧阶段、燃尽阶段5 个阶段；随着温度上升，煤样中的水分蒸发，质量下降，当温度达到213 ℃时，煤样温度达到极小值；随着温度继续上升，煤样吸收空气中氧气，进入吸氧增重阶段，质量略微上升；温度达到372 ℃时，煤样进入热解失重阶段，煤中的芳香环等有机结构受热分解，释放出CO2、CO 等气体并放出大量热量；温度达到523 ℃时，煤样达到着火点，标志着煤样进入到氧化最剧烈的燃烧阶段；温度达到604 ℃时，煤样的燃烧最剧烈；当温度达到766 ℃时，煤样质量不再随温度的升高而变化，煤样进入燃尽阶段。

相比高硫煤，低硫煤的特征温度点均有不同程度提高，煤样的干裂温度升高了24 ℃、煤样的失重温度升高了7 ℃、煤样的着火温度升高了7℃、煤样的最大失重温度升高了11 ℃、煤样的燃尽温度升高了15 ℃。说明高硫煤相比低硫煤有更高的煤自燃倾向。

典型高硫煤和典型低硫煤的DSC 曲线图如图6。

图6 高硫煤和低硫煤的DSC 曲线图Fig.6 DSC curves of high-sulfur coal and low-sulfur coal

由图6 可以看出：随着温度升高，2 种煤样的放热量逐渐增大；温度在500 ℃附近，煤样达到着火点，放热量增幅较大；温度达到600 ℃附近，煤样放热量最大，分别为6.914、5.992 mW/mg。

通过以上分析可以得出：高硫煤在脱水失重、吸氧增重、热解失重和燃烧4 个阶段的放热量相比低硫煤更高，发生煤自燃的风险更高。

煤炭氧化反应的活化能是影响煤炭火灾发生和发展的关键因素之一，活化能是指在化学反应过程中必须克服的能量障碍，也是能使反应发生的最小能量。活化能越小，代表氧化过程越容易发生。煤样热解失重和燃烧阶段的活化能可以通过Coats-Redfern 积分(式（1）)求得。由于煤炭氧化反应属于初级化学反应，故下式中g(a)取-ln(1-a)。

式中：g(a)为煤氧反应的模型函数；a为煤样的转化率，%；T为热力学温度，K；A为指前因子，min-1；R为通用气体常数，7.314 J/(mol·K)；β为升温速率，K/min；E为活化能，kJ/mol。

使用ln[-ln(1-a)/T2]对1/T作图并进行曲线拟合，根据拟合直线的斜率为-E/R即可求出活化能。根据上述方法，得到的煤样热解失重和燃烧阶段的活化能大小见表3。

表3 高硫煤和低硫煤在热解阶段和燃烧阶段的活化能Table 3 Activation energy of high-sulfur coal and low-sulfur coal in pyrolysis stage and combustion stage

通过表3 可以得出：热解阶段中高硫煤的活化能大小为75.112 kJ/mol，低硫煤的活化能大小为82.689 kJ/mol；燃烧阶段中高硫煤的活化能大小为32.767 kJ/mol，低硫煤的活化能大小为38.582 kJ/mol；在热解阶段和燃烧阶段，高硫煤的活化能相比低硫煤分别小9.16%和15.07%。说明高硫煤的氧化过程更容易发生，需要更加注重高硫煤自燃灾害的防治。

2.3 阻化剂对煤中有机硫结构的影响

利用量子化学模拟软件对Ca2+、Na+、Mg2+与二苯硫醚和苯硫醇的络合作用进行计算。结果表明， Ca2+、Na+、Mg2+可以与二苯硫醚和苯硫醇形成络合物，计算得到络合物的分子结构如图7。络合物的键长、键角见表4、表5。

图7 有机硫结构与络合物的几何构型Fig.7 Geometric configurations of organic sulfur structures and complexes

表4 络合物键长表Table 4 Bond length of complexes

表5 络合物键角表Table 5 Bond angles of complexes

从表4 可以看出：3 种离子与2 种有机硫结构形成的配位键键长中，Na+与二苯硫醚和苯硫醇形成的配位键键长最小，分别为0.327 09、0.325 52 nm。说明Na+与2 种有机硫结构形成的络合物结构稳定性更强。

有机硫结构与络合物的静电势投影图如图8，图8 中：红色代表负电性，亲电相关；蓝色代表正电性，亲核相关。

图8 有机硫结构与络合物的静电势图Fig.8 Electrostatic potential diagrams of organic sulfur structures and complexes

从图8（a）可以看出：二苯硫醚中H 原子区域为蓝色，亲核相关；苯环上的C 原子和S 原子区域为红色，亲电相关。从图8（b）～图8（d）中可以看出，二苯硫醚与3 种离子形成络合物后，红色区域与蓝色区域的颜色变淡，说明络合物的形成减弱了二苯硫醚的亲电反应能力和亲核反应能力。从图8（e）可以看出：苯硫醇中H 原子区域为蓝色，亲核相关；苯环上的C 原子和S 原子区域为红色，亲电相关。从图8（f）～图8（h）中可以看出：形成络合物后，S 原子区域的红色明显变淡；H 原子区域的蓝色也一定程度变淡；说明络合物的形成降低了苯硫醇的亲电反应能力，减弱了苯硫醇的亲核反应能力。综上，Na+、Mg2+、Ca2+与2 种有机硫形成络合物后降低了有机硫的化学活性。

在1 个分子中，HOMO/LUMO 轨道决定了分子的主要性质。有机硫结构与络合物的HOMO 图如图9，有机硫结构与络合物的LUMO 图如图10。

图9 有机硫结构与络合物的HOMO 图Fig.9 HOMO of organic sulfur structures and complexes

图10 有机硫结构与络合物的LUMO 图Fig.10 LUMO of organic sulfur structures and complexes

从图9（a）可以看出：二苯硫醚的HOMO 轨道主要由苯环上的C 原子和S 原子提供；从图9（b）～图9（d）可以看出：形成络合物后，络合物的HOMO 轨道由金属原子提供，苯环上的C 原子和S 原子对HOMO 轨道不做贡献；从图9（e）可以看出：苯硫醇的HOMO 轨道主要由苯环上的C 原子和S 原子提供；从图9（f）～图9（h）可以看出：形成络合物后，络合物的HOMO 轨道主要由金属原子提供，小部分由S 原子提供。

从图10（a） 可以看出：二苯硫醚的LUMO 轨道主要由2 个苯环上的C 原子和S 原子提供；从图10（b）可以看出：二苯硫醚与Na+形成络合物后，提供LUMO 轨道的原子相同，但红色和绿色部分体积减小；从图10（c）可以看出：二苯硫醚与Mg2+形成络合物后，2 个苯环上的C 原子和S 原子对LUMO 轨道的贡献减少，LUMO 轨道主要由Mg 原子提供；从图10（d）可以看出：二苯硫醚与Ca2+形成络合物后，LUMO 轨道主要由Ca2+提供；从图10（e）可以看出：苯硫醇的LUMO 轨道主要由苯环上的C 原子提供；从图10（f）可以看出：苯硫醇与Na+形成络合物后，C 原子的贡献降低，Ca 原子的贡献较多；从图10（g）～图10（h）可以看出：苯硫醇与Mg2+、Ca2+形成络合物后，LUMO轨道由Mg 原子和Ca 原子提供，苯环无贡献。

综上所述，二苯硫醚和苯硫醇的轨道主要由S 原子和苯环上的C 原子提供，具有较高的化学活性，容易被氧气氧化；与金属离子形成络合物后，S 原子和苯环上的C 原子对分子轨道的贡献下降，降低了2 种有机硫结构的化学活性。

3 结 语

1)低硫煤相比高硫煤，CO2、CO 煤自燃标志气体释放量较少，煤样氧化过程中特征温度点向高温方向推移，说明高硫煤相比低硫煤发生煤自燃的倾向更高。

2)热解阶段和燃烧阶段，高硫煤的活化能相比低硫煤分别小9.16%和15.07%。说明高硫煤的氧化过程更容易发生，需要更加重视高硫煤自燃灾害的防治。

3)Na+、Mg2+、Ca2+与2 种有机硫形成络合物后降低了有机硫结构的化学活性，对煤自燃起到抑制作用。

4)Na+与二苯硫醚和苯硫醇形成的配位键键长较短，形成的络合物结构稳定性更强。从络合作用的角度， Na2CO3相比Ca(OH)2和MgCl2对预防高有机硫煤自燃的效果更好。