1. 引言
近些年来,随着人口的不断增加和工业化的发展,人们在享受福利的同时,资源的储量也在不断地消耗。据统计,2022年原煤产量达到45.6亿吨,2023年第一季度原煤的产量持续增加,同比增长达到5.5% [1] 。然而有80%的煤炭都被直接燃烧,这导致煤炭资源的严重浪费,并且也会带来一定的环境问题。因此如何解决煤炭的燃烧问题,实现煤炭的高效清洁利用是重中之重 [2] [3] 。煤炭的干馏是当前煤炭加工的最重要的工艺之一,表现为在隔绝氧气或者在惰性气体的气氛下加热到较高的温度,在热解的过程中发生一系列的变化。例如:煤炭的焦化、气化和液化。因此,对煤炭热解特性的研究对于实现煤炭的高效清洁利用具有一定的指导意义 [4] [5] 。
煤炭在热解的过程中可分为挥发性物质的析出以及半焦状残余物的形成,其中挥发性物质包括:水分、气体以及焦油的产生 [6] [7] [8] 。许多学者通过对升温过程中煤的失重来分析来研究煤炭的热解过程,将煤炭的整个热解过程分为几个阶段:煤的脱水、煤炭的解聚和分解以及热缩聚反应 [9] [10] 。其中影响煤热解的因素有,热解终止温度、热解升温速率、样品粒径等。其中热解升温速率对于其热解的影响较为显著 [11] [12] [13] 。通过采用非等温热重法研究不同升温速率的影响下,热解终止温度和热解速率对于热解产物的影响,以及活化能和频率因子在热解过程中的变化规律 [14] 。即样品的升温速率越快,煤中挥发分物质的析出越快,活化能在最大热解失重温度时表现为最大。常娜等 [15] 根据低变质神府煤在不同热解升温速率下的热解失重率进行分析发现:煤热解失重率会随着热解升温速率的提升而升高,但当升温速率达到20℃/min后,热解失重率会趋于平稳,整体失重变化较小。在不同热解升温速率下,李改改等 [16] 采用3种热解动力学模型对不同热解升温速率下的富油煤进行分析,发现活化能在热解升温速率的提升下其平均活化能也会随之增加,并不止单单局限于最大热解失重温度,随着热解升温速率的提升,其热解温度逐渐向高温度区域移动。
基于以上研究现状,采用Coats-Redfern方法,对不同地区、不同煤层、不同煤种的3种煤炭(邯郸峰峰矿区2#肥煤、唐山开滦矿区9#气煤、唐山开滦矿区12#气煤)开展不同升温速率下的热解失重实验,分析其热解过程,探究热解过程中升温速率的改变对热解失重率的影响,进一步划分不同中等变质煤在不同升温速率下的热解失重率关系及最大失重率与升温速率的变化规律。通过热解失重率及失重速率曲线建立热解动力学模型,分析活化能及频率因子在不同升温速率下的变化规律。基于热解过程变化规律及动力学模型参数归纳中等变质煤热解规律特性,为各地煤炭原位注热改性开采提供参考。
2. 样品与实验
2.1. 煤样分析
实验所选用的煤炭为邯郸峰峰矿区2#肥煤、唐山开滦矿区9#气煤、唐山开滦矿区12#气煤。首先基于GB/T212-2008、GB/T214-2007、GB/T476-2008、GB/T19227-2008标准,将煤样使用研钵进行研磨,使煤样粒径 < 0.2 mm,对于研磨后煤样进行工业分析以及元素分析,为后续不同升温速率下热重实验提供基础物质数据。工业分析与元素分析结果如表1所示,其煤质特点表现为低水分,挥发性高。

Table 1. Results of coal industry analysis and elemental analysis
表1. 煤工业分析及元素分析结果
2.2. 不同升温速率下热重实验
不同升温速率下的热重实验使用德国耐驰公司STA 449-F3同步热分析仪。其中煤样所选粒径均 < 0.2 mm。每组实验所选用样品数量为10 ± 1 mg。在实验开始前,首先将煤样置于真空干燥箱中,在110℃真空环境中干燥12小时。干燥后煤样在氮气流量为50 ml/min的气氛下,分别以10℃、20℃、30℃/min的升温速率升高至终止温度。通过以上条件对于三种煤样进行检测分析。
3. 实验结果
3.1. 不同煤热解特性
图1为不同煤在10℃/min升温速率下重量随温度变化曲线(TG)以及其重量变化速率随温度的变化曲线(DTG)。通过TG和DTG曲线可知,不同煤样的热解过程主要可以分为以下3个阶段:当热解温度低于300℃时,煤的失重主要由于煤间结合水以及孔隙中吸附的甲烷,二氧化碳等少量气体。在此过程中峰峰矿区2#肥煤失重率为7.3%。随着热解温度的升高,煤的热解失重速率逐渐提升,当温度到达482℃时,此时峰峰矿区2#肥煤失重速率达到峰值,在这个阶段峰峰矿区2#肥煤失重率为14.23%。在300℃~600℃区间,煤热解失重速率降低,峰峰矿区2#肥煤失重率为26.05%。而后伴随温度的升高,煤热解失重曲线逐渐趋于平稳。

Table 2. Weight loss of coal in different temperature ranges (%) at heating rate of 10℃/min
表2. 10˚C /min升温速率下不同温度区间煤热解失重率(%)
比较3种煤在10℃/min升温速率下,在25℃~600℃温度区间内变化规律,发现开滦矿区9#气煤在235℃~320℃温度区间内,煤热解失重速率曲线(DTG)有一个小型峰值产生,此过程可能由于煤中包括较多的水以及吸附气体。在这个阶段热解失重率为2.6%,根据热解失重曲线可知在324℃~552℃区间内煤热解速率变化最为剧烈,热解速率最大值在458℃,而在整个温度区间内煤的总体失重率为33.03%,与峰峰矿区2#肥煤热解失重曲线较为相似。通过对比开滦矿区9#气煤与开滦矿区12#气煤热解失重曲线可知:在25℃~600℃区间内开滦矿区12#气煤热解失重率高于开滦矿区9#气煤,据(表1)煤的工业分析与元素分析可知,此过程的变化是由于开滦矿区12#气煤中硫元素含量高于9层煤,在热解区间内挥发性物质受热分解所导致。再根据(表2)所示,不同煤样在10℃/min的升温速率下,在到达600℃时,比较峰峰矿区2#肥煤、开滦矿区9#气煤与开滦矿区12#气煤热解失重曲线与失重速率曲线可知:其整体失重率较为相似。通过(图1)可知3种煤最大热解反应速率区间均在330℃~560℃。
(a) 10℃/min峰峰矿区2#肥煤TG-DTG热解失重曲线 (b) 10℃/min开滦矿区9#气煤TG-DTG热解失重曲线
(c) 10℃/min开滦矿区12#气煤TG-DTG热解失重曲线
Figure 1. TG-DTG weight loss curves of different coal samples at 10˚C/min
图1. 10℃/min不同煤样TG-DTG失重曲线
3.2. 不同升温速率对于煤热解特性影响
在煤的热解过程中,升温速率是最重要的因素之一。在煤的整个热解的过程中,可以将其分为3个阶段,如(图2)为3种煤在10℃/min、20℃/min、30℃/min热解速率下的TG曲线图,从图中可以看出,煤热解的失重会随着热解速率的提升而降低,在30℃/min升温速率下的失重率最低,在10℃/min升温速率下的失重率最高。在不同升温速率的影响下,TG曲线形状基本一致,但随着升温速率的提升,其整体失重率也随之下降。也就是说在保证失重率相同时,所对应的热解温度也会提升,升温速率的提升不利用煤的充分热解。
煤热解特征参数的求取可以反应热解过程中不同阶段的性质,通过比较不同升温速率下中等变质煤的起始失重温度,可以直观地反映原煤在热解过程中分解特性,通过中等变质煤在反应过程中的最大失重温度来判断其大分子结构之间的稳定性。通过比较相同煤样的不同升温速率下最大失重温度来判断煤的热解反应特性,最大失重温度越高,即更加利于煤的热解,反之抑制。
根据热重实验得出热解的TG与DTG曲线,从而分析热解过程中煤的特征参数,中等变质煤的特征参数定义如下:Ts为热解初始温度,Tmax为热解过程中最大失重温度,Te为剧烈热解温度,其值为Te = 2Tmax − Ts,(dw/dt)max为最大热解失重速率,此点为TG曲线上反应最为剧烈的点,也就是热解反应速率曲线上的峰值。其具体数值如下表3所示,可以发现,随着升温速率的改变,其初始析出温度始终维持在一个稳定的范围内,受到升温速率的影响较小。而最大失重温度与剧烈热解温度均随着升温速率的提升,向后偏移,表明在热解过程中升温速率的提升不利于煤的充分热解。

(a) 峰峰矿区2#肥煤不同升温速率下TG-DTG曲线
(b) 开滦矿区9#气煤不同升温速率下TG-DTG曲线
(c) 开滦矿区12#气煤不同升温速率下TG-DTG曲线
Figure 2. TG-DTG curves of different coal samples with different heating rates
图2. 不同煤样不同升温速率的TG-DTG曲线

Table 3. Pyrolysis characteristics of raw coal
表3. 原煤热解特征参数
3.3. 热解动力学参数及分析
研究煤的热解动力学对于了解热解机理,实现资源的清洁高效利用,具有深刻的意义。通过对于不同煤样在不同升温速率下的热重实验数据进行动力学分析,通过对比活化能、频率因子等动力学参数,来探究升温速率对于煤热解过程中的影响。求解热解动力学参数的方法主要有两种,分别为微分法和积分法,积分法更加侧重于整体数据的分析,不单单局限于微分法的局部数据分析,且积分法处理数据较为简单,对原始数据造成的误差较小,基于以上两点,本文采用Coats-Redfern算法来进行热解动力学参数的计算。
总反应速率方程 [6] :
(1)
(2)
其中:a为煤在热解过程中的转化率,%;t为煤热解反应过程时间,min;A为频率因子,min−1,n为热解反应级数;E为活化能,KJ/mol;R为气体摩尔常数,数值为8.314 J (mol∙K);w0为样品初始质量,g;w为样品在热解过程中任意时刻的质量,g;wf为热解结束后样品的剩余质量,g。
在非等温热重试验中,升温速率为固定数值的情况下,热解温度T与时间t成线性关系,可以表示为:T = T0 + βt,β为升温速率,K/min。
若将T = T0 + βt带入方程(1),可得:
(3)
公式(3)中f(a)为机理函数,假设热解反应过程为n级反应,则f(a) = (1 − α)n,对公式(3)进行取对数变化得:
(4)
(5)
其中,公式(4)、(5)为Coats-Redfern方法的通用表达式,通常可以公式(4)右边第一项lnAR/βE看作常数,因此
可以对1/T做出一条直线,其斜率为−E/R,截距为ln(AR/βE),从而可以计算出原煤热解过程中的活化能E和指前因子A的大小。所选温度区间皆为最大失重温度周围,可以发现随着升温速率的提升,煤样的活化能都变现为上升的趋势,根据表4、表5可知,峰峰矿区2#肥煤平均活化能为45.75 KJ/mol,开滦矿区9#气煤的平均活化能为47.06 KJ/mol,开滦矿区12#气煤的平均活化能为57.94 KJ/mol。

Table 4. 3 Dynamic parameters of different coal heating rates
表4. 3种煤不同升温速率下动力学参数

Table 5. 3 Dynamic parameters of different coal heating rates (continued)
表5. 3种煤不同升温速率下动力学参数(续)
基于3种中等变质煤不同升温速率下的动力学参数,可知不同煤样不同升温速率下的热解动力学方程方程如下表6所示:

Table 6. Kinetic equation of pyrolysis reaction
表6. 热解反应动力学方程
4. 结论
1) 从3种中等变质煤在不同升温速率下的TG、DTG曲线可以得知,不同升温速率下的中等变质煤在热解过程中所表现的热解特性不同,大致可分为3个阶段:脱气脱水阶段;热解活跃阶段;二次脱气阶段。且从DTG曲线可以得知,峰峰矿区2#肥煤最大热解速率集中在480℃~490℃之间,开滦矿区9#气煤、开滦矿区12#气煤最大热解速率均集中于450℃~475℃之间。
2) 根据中等变质煤不同升温速率的TG、DTG对比可知,随着热解升温速率的提升,热解失重量以及热解失重速率峰值均向后移动,即在热解失重量相同的情况下,热解升温速率越高,其所需的温度越高。
3) 热解动力学参数的计算采用Coats-Redfern算法来对中等变质煤热解最大速率阶段进行动力学分析,得到不同煤样在不同热解升温速率下的最大热解速率阶段的动力学参数。由动力学参数可知:在固定煤样的情况下,随着热解升温速率的情况下,中等变质煤的活化能和频率因子会随着热解升温速率的增大而升高。
基金项目
河北省重点研发计划项目(21374104D);河北省自然科学基金(E2019402361、E2020402075);河北省教育厅科学研究项目(ZD2021309)。
NOTES
*通讯作者。