不同气氛下富油煤受热裂隙演化及热解动力学参数变化

王双明 ，孙 强 ，胡 鑫 ，葛振龙 ，耿济世 ，薛圣泽 ，师庆民

（1.陕西省煤炭绿色开发地质保障重点实验室, 陕西 西安 710054；2.西安科技大学 煤炭绿色开采地质研究院, 陕西 西安 710054；3.西安科技大学地质与环境学院, 陕西 西安 710054；4.山西大同大学 建筑与测绘工程学院, 山西 大同 037003）

0 引 言

受“缺油、少气、相对富煤”的能源资源禀赋特征制约[1]，我国油气长期保持着较高的对外依存。为弥补油气能源不足的现状，煤制油有望成为实现油气补给的一种有效途径。由于我国富油煤资源储量丰富[2]，焦油产率高（Tar,d>7%），这一特殊煤炭资源为煤制油制气提供了先天优势[3]。富油煤是集煤、油、气属性于一体的煤基油气资源，但直接作为燃料的传统方式未能使其高效利用，造成了能源浪费[4]。为响应煤炭低碳、清洁、高效开发转型的国家要求，探索富油煤新的开发利用方式具有现实意义。

热解可以通过升高温度和压力的方式将富油煤转化为更有价值的产品，如煤焦油、煤气和焦炭，从而实现资源的多元化与最大化利用[5]。其次，热解可以有效地分解并去除富油煤中的有害物质，降低对环境的污染，因此热解成为了适合富油煤的高效利用方法。富油煤主要赋存于中低阶煤类中，具有水分高、挥发分高、氧含量高的特点。但相较于传统低阶煤,富油煤具有较大的孔比表面积与孔体积[6]。其含有大量的富氢结构，热解后可生成较多的CH4、H2、CO 和烃类气体与煤焦油[2]。而热解的温度阶段与气体产率等特征与煤化程度、热解温度、煤粒径等因素相关。且在热解高温影响下富油煤的结构、成分、性质等也会发生变化，影响热解的进程和产物分布[7-8]。为推进富油煤热解制油制气的开展，解决热解过程中富油煤的高温响应成为了重要的科学问题[9]。

富油煤原生裂隙较少，作为煤热解油气的运移通道[10]，裂隙越丰富，富油煤提油效率越高[11-12]。通过改造富油煤裂隙，能够提高富油煤的渗透性[13-14]。而热解的高温可以促进富油煤裂隙的扩展与发育[15]，李利峰[16]、李海琪[17]等对褐煤、长焰煤等热解条件下裂隙的演化作了探索，但缺少不同气氛条件下的对比研究。富油煤地下原位热解过程中存在低氧与无氧高温区域，不同的气氛条件会对煤体裂隙的发育产生影响。另一方面，热解过程中富油煤的动力学特征也是一个亟待探索的关键科学问题。毛崎森等[18]、董光顺等[19]、陈美静等[20]利用数值模拟技术对富油煤原位热解过程中温度场分布、传热效率、余热利用效率等进行了研究，但关于热解过程中热导率、热失重、活化能等参数的变化规律有待深入探究。

笔者从探索高温作用下富油煤裂隙演化与热解动力学参数响应角度入手，选取柠条塔矿富油煤，利用体视显微镜和Matlab 研究富油煤表面裂隙变化规律，并利用激光闪射法监测富油煤加热过程中的实时热导率变化，分析富油煤热解失重（TG）、微商热重分析（DTG）、差示扫描量热（DSC）曲线以及转化率α随温度的变化特征，基于热解动力学分析富油煤热解特征与规律。研究结果将为富油煤热解提供重要的理论和实践指导意义。

1 试验样品与方法

1.1 试验样品

采样点位于陕西省神木县柠条塔煤矿，取2-2 煤层进行研究，对煤样进行工业分析和元素分析等测试，结果见表1。

表1 富油煤分析测试结果Table 1 Characterization of tar-rich coal sample

1.2 试验方法

将采回的原煤切割制成约20 mm×20 mm×20 mm 的块状煤样并分为2 组，一组样品置于欠氧气氛的BLMT-1400 ℃高温节能箱式炉（图1a），另一组样品置于N2气氛的OTF-1200X 石英管式马弗炉中（图1b），均以5 ℃/min 的升温速率进行加热，达到预设温度后炉内恒温2 h，恒温完成后炉内降温至室温。焦油蒸汽与脂肪烃在炉内冷凝（图1c），煤样经过不同温度热处理后产生不同程度破裂（图1d）。

图1 试验仪器与过程Fig.1 Experimental instrument and process

煤样冷却至室温后，利用体视显微镜（图1e）对比了不同温度处理后的表面微观图，并利用数字图像处理（Matlab）技术对显微图像进行分析。数字图像处理通过计算机对图像进行去除噪声、增强、复原、分割、提取特征等方法和手段，将图像信号转换成数字信号并对其进行处理，提取图像中所包含的特征与信息，将图片信息可视化，并量化这些特征与信息[21-23]。

将煤样磨成直径小于0.5 mm 的粉末（图1f），取用10 mg 的粉末煤样装入热重仪的小坩埚，放入TGA/DSC1/1 600 型热重及同步热分析仪中（图1g），设定加温速率为10 ℃/min，最高温度为1 200 ℃。测试环境为氩气气氛，流量60 mL/min，进行热稳定性测试。单盘配有2 对Pt-Pt/Rh 热电偶用于实时测试温度。可以输出TG，DSC 曲线以及DTG 曲线，可同时测量热流以及重量变化。

取部分原煤和热处理后的煤样进行粉碎、研磨和筛分，并将得到的不同条件下的煤粉（d<0.075 mm）压 制 成Ø12.7 mm×1 mm 的 圆 薄 片 样 品（图1h，图1i）。利用LFA457 激光闪射仪（图1j）进行热导率测试，试验过程中以100 mL/min 速率的氮气吹拂样品。

2 试验结果与分析

2.1 煤样裂隙特征

使用Matlab 数字图像技术，将体视显微镜获取的煤样表面显微图处理为灰度等级只有0 和1 的二值矩阵图像，黑色裂隙部分用0 来表示，其余部分用1 表示[24]。因此，可以定量分析煤岩表面裂隙随温度的发育过程与程度，见表2 与表3。

表2 欠氧气氛热处理后煤样表面裂隙发育Table 2 Development of surface cracks in coal samples after heat treatment under air atmosphere

表3 N2 气氛热处理后煤样表面裂隙发育Table 3 Development of surface cracks in coal samples after heat treatment under N2 atmosphere

在室温至150 ℃之间，脱吸附水与吸附气阶段煤体裂隙变化不明显，而在150～325 ℃，随着温度上升，CO、CO2、CH4等挥发分热解气体快速逸出[25]。在欠氧气氛中，煤样表面与氧气发生氧化反应，产生压力梯度，突破孔隙，致使裂隙数目明显增多，并逐渐贯通[26]。逸出的挥发分在煤块附近燃烧，进一步提高了煤表面的温度。300 ℃后煤样表面受高温影响产生琉璃色斑，至325 ℃时部分煤体棱角或表面粗糙处发生燃烧（图2）。

图2 高温处理后煤样表面现象Fig.2 Surface morphology of coal samples after high-temperature treatment

在N2气氛中，煤样在200～350 ℃已完成脱水，并开始析出少量热解气，但热解程度较低，煤样表面仅产生少量细小裂纹。从400 ℃开始热解气逸散加剧，煤样表面裂隙数量和裂隙宽度明显增加，煤体掉渣。随着加热温度的进一步提高，煤体裂隙继续增加并相互贯通，部分裂隙加深至内部，煤块破碎分裂严重，且镜煤的裂隙发育更强烈（图3）。

图3 N2 气氛500 ℃下暗煤与镜煤高温处理后裂隙发育状况Fig.3 Development of cracks in dark coal and specular coal after heat treatment at 500 ℃ under N2 atmosphere

Matlab 获取的灰度值可以量化煤样表面裂隙率（图4），在欠氧气氛下150 ℃以下煤裂隙率增加缓慢，而在N2气氛下，400 ℃以下煤裂隙率增加幅度均较低，且同一温度下的裂隙率明显低于欠氧气氛，因此N2环境提高了裂隙率的拐点温度，降低了裂隙的增加程度。但随着加热温度超过拐点，2 种气氛下煤的裂隙率均进入快速增长阶段。

图4 不同温度处理后煤样表面裂隙率与质量损失率Fig.4 Surface crack rate and quality loss rate of coal samples treated at different temperatures

质量损失率规律与裂隙率相似（图4），同一温度下欠氧气氛的质量损失率明显高于N2气氛，拐点温度为200 ℃，而N2气氛拐点温度为300 ℃。超过拐点温度后欠氧气氛的质量损失率的增长速率明显高于N2气氛。

2.2 煤样热导率

激光闪射法在以往的研究中被广泛用于不同煤样实时热导率的测试[27]。图5 展示了笔者及其他研究中不同类型煤利用激光闪射法得到的实时热导率λ。由于种类和开采矿区不同，煤样在持续升温过程中的实时热导率也存在显著区别，如图5a 所示。对各煤样实时热导率根据最低温度值进行归一化分析（图5b），结果表明，在室温至350 ℃范围内的持续升温过程中，热导率与温度成正相关。

图5 不同煤样的实时热导率曲线[26-30]Fig.5 Real time thermal conductivity curves of different coal samples[26-30]

将经过不同温度热解后的煤样热导率和原煤在实时热导率测试中的数据进行对比（图6），发现热解后煤样与实时测试数据随温度的变化规律相反。这是由于在实时热导率测试的持续升温过程中，煤样体积不断膨胀，导致热导率不断升高[28]。而对于热解后的煤样，煤样膨胀造成的热导率升高现象在热解后的冷却阶段逐步消失，热解温度的升高造成了声子振动幅度增大，声子平均自由程减小，导致煤样热导率降低[29]。

图6 不同温度热解后煤样热导率与原煤实时热导率变化Fig.6 Variation in thermal conductivity of coal samples and real-time thermal conductivity of pyrolized coal at different temperatures

图7 展示了经过不同温度氧化和热解处理后煤样热导率的变化。随温度的升高，氧化煤样的热导率在100 ℃和275 ℃时分别达到最大值0.108 5 W/(m·K)和最小值0.083 W/(m·K)，热解煤样的热导率始终随温度升高而减小，从常温下的0.097 W/(m·K)降至350 ℃下的0.084 W/(m·K)。

图7 氧化和热解后煤样的热导率Fig.7 Thermal conductivity of coal samples after oxidation and pyrolysis

利用激光闪射法测定热导率的过程中，煤样经过研磨并重新压制成样，煤样孔隙率远小于原煤，其热导率数据（图7）可近似视为固相煤热导率λS。煤高温处理中（后）存在大量孔隙和裂隙[30-31]。应用式（1）所示的多孔介质热导率模型更为符合实际情况[32]：

式中：λP为煤的有效热导率；φC为沿孔隙流管截面的孔隙体积比；λa为孔隙流管部分有效热导率。

考虑煤样孔隙中空气对热导率的影响，利用固-流二相模型计算了煤中各相物质在串联分布下的最小有效热导率λmin和并联分布下的最大有效热导率λmax2 种极限分布模式下的热导率公式[33]（如式（2）和式（3）所示）。

式中，λf为空气的热导率。而煤样在实际中的固-流相分布模式总是介于2 种极限模式之间，利用式（4）计算得到氧化和热解处理后煤样的有效热导率，并利用式（5）计算出热导率变化率τ，结果如图8 所示。煤样的有效热导率小于激光闪射法测得的值，热导率变化率在氧化条件下随氧化温度的升高显著增大，在热解条件下变化较小。

图8 考虑孔隙中气氛影响的煤样有效热导率Fig.8 Effective thermal conductivity of coal samples considering influence of atmosphere in pores

2.3 氩气环境中煤样热分析特征

图9 为煤样的TG、DTG、DSC 曲线以及转化率α的变化。在25～600 ℃的热解温度区间，可分为25～300 ℃、300～500 ℃和500～600 ℃ 3 个阶段。在阶段Ⅰ，主要是吸附水析出，升温速率对TG曲线无明显影响，发生吸热反应，DSC 曲线呈下降趋势。阶段Ⅱ中，有机质开始分解，同时部分矿物析出结构水，发生脱羟基的吸热反应[34-36]；DTG 曲线表明煤样在热解过程最大失重速率位于435 ℃左右；随着升温速率的增大，DTG 曲线波谷逐渐向高温区移动。阶段Ⅲ，煤样中有机物已基本分解完毕，无机矿物开始分解并放热[37-38]，DSC 曲线上升。图9d 表明：转化率随温度的升高而增大：阶段Ⅰ（α<0.2），随温度的升高，转化率呈缓慢上升趋势；阶段Ⅱ（0.2<α<0.8）转化率迅速上升；阶段Ⅲ（0.8<α<1.0）转化率呈现缓慢上升。

图9 升温速率对TG、DTG 和DSC 和转化率的影响Fig.9 Effect of heating rate on TG, DTG, DSC and conversion rate

2.4 煤样热解动力学特征

在热反应过程中，t时刻的质量转化率α为

式中：m0为样品的初始质量；mt为t时刻样品的质量；m∞为样品的最终质量。

采用Coats-Redfern 积分法[39]，进行热解动力学分析[40]，可得式（7）：

式中：A为频率因子，min-1；E为活化能，kJ/mol；R为摩尔气体常数，8.314 J/(mol·K)，T为绝对温度，K；β为升温速率，K/min。

热反应机理函数G(α)可取-ln(1-α)[41]，由于式（7）可化简为式（8），各温度段ln[G(α)/T2]与1/T关系如图10 所示。

图10 各温度段ln[G(α)/T2]与1/T 的线性拟合结果汇总Fig.10 Linear fitting of ln [G(α)/ T2] with 1/T in different temperature ranges

定义参数b为式（9），根据式（10）和式（11）计算活化能E和频率因子A[42]。

不同加热速率下E和A随温度的变化如图11所示，在阶段Ⅰ，随着温度的升高，E和A逐渐增大；阶段Ⅱ，随着有机质开始热解，E和A继续增加，在450 ℃左右达到峰值；阶段Ⅲ，有机质基本被分解，E和A降低。

图11 活化能和频率因子随温度的变化Fig.11 Variation of activation energy and frequency factor with temperature

3 结 论

1）随着加热温度的升高，富油煤样表面裂隙数目明显增多，并逐渐贯通。相比N2气氛，同温度下欠氧环境中煤样裂隙率更高。高温下热解水与热解气体逐步脱离煤体逸散，使富油煤样质量变化率持续增加，欠氧环境中增幅更显著。

2）室温至350 ℃范围内的持续升温过程中，热导率与温度成正相关，热解后煤样与实时测试数据随温度的变化规律相反。富油煤样氧化后的热导率在100 ℃和275 ℃时分别达到最大值和最小值。

3）根据富油煤样的TG、DTG、DSC 曲线以及α可将热解温度分为区间25～300 ℃、300～500 ℃和500～600 ℃ 3 个阶段。热解动力学表明热解阶段Ⅰ和阶段Ⅱ中，E和A随温度的升高逐渐增大，在450 ℃左右达到峰值，但在阶段Ⅲ中呈降低趋势。

致谢 感谢西安科技大学侯恩科教授和王生全教授对本文给予的大力支持和帮助，感谢博士研究生袁士豪和张鹤对本文图表的编辑工作。