稻秆的热解特性及热解模型的建立

1、Zhaochuang JiangenchenChenliangguangXiongleiming (Engineering college in south China Agricural University ,guangzhou, China，510642)Email：ecjiang，zhao020057 , Abstract: By using STA449C Jupiter®thermogravimetric apparatus, effect of the heating rateand ratio of adding catalyst for rice straw wer

2、e studied, which showed that as the heating rates increased,weight loss was little,weight loss rate changed heavily,;after adding catalyst, retention rate increased( yield of carbon rised), activation energy of pyrolysis lowered,as ratio of adding catalyst increased, activation energy of pyrolysis r

3、ised, activation energy lowered,then rised .By Malek method , it could be determined that pyrolysis process of rice straw accorded with the random nucleation and later growth. By estimating trial data,pyrolysis reaction series of rice straw was 0.4, pyrolysis reaction series of rice straw plus 10% c

4、atalyst was 0.46, pyrolysis reaction series of rice straw plus 15% catalyst was 0.42.By estimating pyrolytic reaction of rice straw, which were 184.07 kJ/mol, 158.70 kJ/mol and 172.60 kJ/mol, andin accord with Coats and Redfern method . Correctness and applicability of the mechanism model were furth

5、er confirmed by pyrolysis of rice straw segment model calculation, reliability of the mechanism model was vetified by Kinetics software.Keywords: pyrolysis,catalyst, mechanism,kinetic,pyrolysis parameters,稻秆的热解特性及热解模型建立赵创蒋恩臣 陈亮广熊磊明华南农业大学工程学院，广州 ，中国，510642Email：ecjiang，zhao020057 , 摘要：该文采用耐驰STA 449C

6、Jupiter®同步热重分析仪，研究了升温速率和添加催化剂比例对稻秆的的热解过程的影响。结果表明，随着升温速率的增大，失重量变化不明显，但失重速率变化显著；添加催化剂后，热解的残留率增大，即焦炭产率上升，随着催化剂添加量的增加，活化能先降低后增加。用Malek法确定其热解机理模型为随机成核和随后成长，稻秆热解的反应级数n为0.4，加入10%催化剂的稻秆的反应级数n为0.46，加入剂15%催化剂的稻秆的反应级数n为0.42。通过稻秆热解反应公式拟合的方法得到的主热解阶段活化能与Coats 和Redfern 法的计算结果相符，其活化能依次为184.07 kJ/mol、158.70 kJ/

7、mol和172.60 kJ/mol。通过稻秆热解的分段模型计算验证其热解机理的正确性及适用性，用Kinetics软件验证其模型的可靠性。关键词：热解，催化剂，反应机理，动力学，1引言对化石燃料急剧的需求和环境问题的担忧，促使人们对可再生能源之一的生物质能广泛关注，其在整个能源结构中占有重要地位。可以将生物质通过不同的方式和机理转化成生物质燃料，例如，H2，CH4，乙醇或者气化燃气等1-4。燃气可以用作动力装置燃料，输出电力或动力，国外在生物质发电方面已经取得一些成就，国内发展较晚，但是发展较快5。热解是生物质热化学转化技术之一，近些年受到了人们的广泛关注6-12。热解动力学的研究不仅有助于了解

8、生物基金项目：广东省科技厅国际合作项目（2009B050600008）质热解过程，而且也将为生物质热转化技术所涉及的优化产物组成、反应条件及反应器设计选型提供重要的理论依据13。当前，国内外对于生物质催化热解反应动力学开展了大量的研究14-19，而基于热重分析对稻秆展开的关于催化热解的研究比较少。因此本文以稻秆为原料，在热重热差分析仪器上进行试验，分析添加催化剂前后稻秆的热解特性和动力学，通过稻秆热解的分段模型计算验证其热解机理的正确性及适用性，用Kinetics软件验证热解模型的可靠性。2试验原料和方法2.1 试验原料试验所用稻秆来自华南农业大学的试验田，经研碎为粉末状，在烘箱干燥40h，试

9、验所用秸秆均过80目标准筛。本文借鉴煤的工业分析标准进行测量，其中水分采用标准GB/T211-1996中的空气干燥法测量，灰分采GB/T212-2001中的快速灰分法测量，挥发分采用GB/T212-2001中的方法测量；采用美国EAI公司EA-2400型快速元素分析仪分析秸秆的主要元素。结果列于表1。Table.1 Proximate analysis ,elemental analysis and quantity of rice straw表1稻秆的工业分析和元素分析/%水分/%灰分/%挥发份/%固定炭/%C/%H/%N/%S/%O/%5.3515.1567.4517.440.327.64

10、31.1620.48649.64试验所用的催化剂是450煅烧的膨润土，热稳定性比较好。其成分以蒙脱石为主，MgO、CaO 和Fe2O3等一些成分在一些研究过程中有催化裂解的作用。其成分分析列于表2。Table 2 Main composition of swell soil表2 膨润土的主要成分成分Al2O3FeOSiO2Fe2O3MgOCaOK2O含量16.54%0.26%50.95%1.36%4.65%2.26%0.47%2.2 试验方法热重试验在德国耐驰公司STA 449C Jupiter®同步热分析仪上完成。以10、20和30·min-1从室温升温到800, 加热炉

11、通入氮气保护气，流速为20 ml·min-1，样品量约为5mg，实验参比是Pt，记录热解曲线，分析它的热解特性。3 试验结果分析3.1稻秆热解曲线的分析稻秆在不同升温速率下的失重(TG)和失重速率(DTG)曲线如图1和图2所示。从图1和图2可以看出稻秆的热解可以分为三个阶段：失水预热阶段、主热解阶段和炭化阶段。生物质的化学组成主要有中性洗涤剂溶解物、半纤维素、纤维素和木质素组成，一般农业秸秆类物质主要成分是后三种物质，升温速率对富含纤维素的物质有很大的影响。随着升温速率增大，热解起始温度高，主热解最大峰向高温区转移，且残留率越大；随着升温速率减小，热解起始温度低，物料吸热较慢，主热解

12、最大峰向低温转移，残留率较低，热解较充分9、20-23。分析认为这主要是由于升温速率增加，导致颗粒内部的热解气体未能及时扩散，影响了内部热解的进行13。表3列出了稻秆的热解特性。Fig1 TG curve of rice straw pyrolysis Fig2 DTG curve of rice straw pyrolysis图1 稻秆热解的热重曲线图2 稻秆热解的失重速率曲线Table.3 pyrolysis speciality and characteristic parameters of main pyrolysis rice straw表3稻秆的热解特性和主热解阶段的特征参数热解

13、步骤升温速率/minTe/Tf/Tp/（da/dt）max /%.min-1失重/%失水预热阶段103018047-0.557.762030186681.607.353030203682.027.57主热解阶段101803883206.5552.972018639533213.2053.613020340833718.5053.68炭化阶段10388800-18.6220395800-14.9230408800-10.44注：Te主热解起始温度，Tf主热解结束温度，Tp主热解最大失重速率对应的温度以图1和图2中升温速率10/min曲线为例进行分析，在180以前为失水预热阶段，样品失去自由水，一

14、些不稳定成分开始发生分解反应，成分变化不大5，但内部发生重组，例如结聚和玻璃化转变9、25。随着升温速率的增大，对失重没有大的影响，平均失重7.75%，失重速率变化明显，出现失水滞后现象。180388为稻秆的主热解阶段，此阶段样品发生复杂的化学变化，失重主要是半纤维素、纤维素和木质素等物质的分解5、8。随着升温速率的增大，失重量变化不明显，样品平均失重为53.42%，但失重速率变化显著。388以后为炭化阶段，样品失重速率随升温速率的增大而升高。随着温度升高，木炭中的挥发物被排出，提高固定炭含量5，炭化阶段是热解反应中最耗能的阶段。3.2稻秆催化热解曲线的分析添加催化剂（膨润土）后，在20/mi

15、n的升温速率下，稻秆热解的TG线和DTG线如图3和图4所示，其热解曲线的趋势和未添加催化剂一样，分为三个阶段，添加的催化剂比例分别为10%和15%。 从图3可以看出，添加10%催化剂的残留率为35.72%，添加15%催化剂的残留率为43.75%，其明显比未添加催化剂的残留率增大，即焦炭产量增加，说明添加催化剂有助于焦炭产量的增加。当添加比例越大时，催化作用明显，但热传质阻力增大，影响热解挥发物的溢出。从图4可以看出，随着添加催化剂的比例增加，样品在主热解阶段的失重速率逐渐减小，添加10%催化剂的最大失重速率为10.95%/min，20%的最大失重速率为10.01%/min，最大热解峰变逐渐变窄

16、，与热传质阻力有关；热解时间缩短，这认为是与添加催化剂相关，催化剂降低了热解反应的活化能，稻秆甚至会直接热解。从图4可以看出，120前，添加催化剂量对热解过程的影响不明显，120332之间，催化剂量对热解过程影响较明显，332后，催化剂量对热解过程影响不显著；未添加催化剂的热解过程中，在300左右有一个肩峰，在添加10%催化剂的情况下肩峰不明显，但随着添加催化剂比例增加，其肩峰又变明显。由此推论，添加催化剂比例越大，催化作用是先强后弱。Fig 3 TG curve of rice straw pyrolysis after adding catalyst图 3 添加催化剂后稻杆热解的热重曲线F

17、ig 4 DTG curve of rice straw pyrolysis after adding catalyst图 4 添加催化剂后稻杆热解的失重速率曲线4主热解阶段反应机理模型确定生物质热解机理函数有多种方法，例如，王明峰等8采用Malek法确定玉米秸秆热解机理函数。本文采取Malek法选择主热解段的机理函数。由于此法已是被人熟知的方法，在此不具体列出来。若试验曲线与某标准曲线重叠或者接近，或试验数据点全部落在某一标准曲线上，则判定该标准曲线所对应的和就是最概然动力学机理函数。把稻秆热解试验曲线与41种动力学机理函数的标准曲线进行对比，结果表明，稻秆热解试验曲线与J-M-A方程的标准

18、曲线最接近，如图5所示，因此推断稻秆的热解最概然动力学机理为随机成核和随后成长，其机理函数为：（5.1）（5.2）稻秆的热解反应可表示为：（5.3）将试验数据带入式5.3进行拟合，可得到反应动力学参数。稻秆热解动力学参数为：lgA=25.0561 s-1，E=184.07 kJ/mol，n=0.4； 添加10%的催化剂后的热解动力学参数为： lgA=18.2632s-1，E=158.71 kJ/mol，n=0.46； 添加15%的催化剂后的热解动力学参数为： lgA=25.8765 s-1，E=172.60 kJ/mol，n=0.42。a无催化剂b 催化剂10%c 催化剂 15%Fig5 St

19、andard curves and experiment curves of vs图5 标准曲线和试验曲线5 机理模型的验证及适用性计算结果表明，本文通过稻秆热解的分段模型计算出升温速率为10/min时不同温度对应的样品量，并与试验结果进行对比。计算方法如下：稻秆主热解阶段反应的积分形式为：（5.1）将的Doyle近似式带入式5.1，整理得：（5.2）任意时刻样品质量的百分比为：（5.3）将稻秆主热解阶段的E、A、n分别带入式5.1和5.2，取=10/min，给定T值可计算出，然后根据式5.3计算T对应的。10/min稻秆主热解阶段试验结果与计算结果对比如图6，稻秆的热解试验数据与计算结果吻合

20、，其平均误差在2%以内，说明热重分析的结果能正确反映稻秆的热解反应。Fig. 6 Comparison of counted data with experiment data ofmain pyrolysis stage of rice straw图6稻秆主热解阶段的试验数据与计算数据对比本文了通过Kinetics软件预测稻秆在升温速率为20/min下的热解数据，试验数据和预测数据见图7。结果表明，在6.50min内完成主热解阶段反应，该结果与试验数据相吻合。说明所建的模型的应用范围不仅仅局限在建模所用的升温速率范围内，而且对于其他升温速率仍具有较高的可靠性。a 主热解段预测数据b 主热解段

21、试验数据Fig.7 Comparison of simulated data with experiment data of main pyrolysis stage of rice straw图720/min稻秆主热解阶段的试验数据与预测数据对比6动力学分析常用的确定动力学参数的积分求解方法有Doyle 和Zsako 法、 Coats和Redfern 法、 MacCallum 和Tammer 和Ozawa法等。本文选用Coats 和Redfern 法积分公式计算了稻秆热解的反应活化能和指前因子。（6.1）由于式6.1中近似等于1,则是1/T的线性函数,其斜率为-E /R,而截距中包含指前因子

22、A。由此可从斜率、截距求得E、A,计算结果列于表4。对表4进行分析，加入催化剂前后，用Coats和Redfern法计算的活化能值与拟合求解的结果一致，确定了稻杆热解动力学参数。随着升温速率的增大，活化能呈波动状态，一方面由于秸秆中含有半纤维素、半纤维素和木质素等，热解过程比较复杂，影响活化能，另一方面由于升温速率增大影响热传质，必然影响热解过程。加催化剂并没有改变稻秆的热解机理，而是对热解反应进行了优化，添加适当比例的催化剂，降低了热解反应的活化能，缩短主热解阶段的时间。添加10%的催化剂后，稻秆热解的活化能明显降低，与热解曲线展现的结果一致；随着催化剂比例的增加，活化能增大，是由于催化剂比例

23、增加，导致传质阻力增大，影响热解反应。相关系数r 几乎都在0.991 之间，说明得到的结果是可靠的。Table 4 Constant of rice straw pyrolysis表 4 稻秆热解动力学常数项目升温速率/min催化剂催化剂比例/%热解温度范围/E/kJ/molA/s相关系数r稻秆10180-388184.843.0850E+150.99820186-395179.433.0825E+150.99730203-408187.719.7681E+150.996稻秆+催化剂20膨润土10187-397158.704.0745E+130.99715188-399172.596.3531

24、E+140.9977结论1）对稻秆进行热解分析，稻秆的热解反应分为三个阶段：失水预热、主热解和炭化三个阶段，随着升温速率的增大，每个阶段的失重量变化不明显，但失重速率变化显著；添加催化剂前后，热解趋势保持一致。2）添加催化剂前后稻秆热解反应的活化能波动，未添加催化剂时，E为184.84187.71 kJ/mol；添加10%的催化剂后，受到催化剂的影响，活化能降到158.70 kJ/mol，再增加催化剂的比例，活化能又增加到172.59 kJ/mol。3）用Malek法确定稻杆的热解机理为J-M-A方程，为随机成核和随后成长。通过稻秆热解的分段模型计算验证其热解机理对建模升温速率以外的其它升温速

25、率具有正确性及适用性。References(参考文献)· 1R. Lemus and R. Lal, “Bioenergy crops and carbon sequestration,” Critical Reviews in Plant Science vol. 24, Feb. 2005, pp. 121, doi: 10.1080/07352680590910393.· 2S.H. Beis, Ö. Onay and Ö.M. Koçkar, “Fixed-bed pyrolysis of safflower seed: influen

26、ce of pyrolysis parameters on product yields and compositions,” Renewable Energy，vol. 26, May 2002, pp. 21-32, doi:10.1016/S0960-1481(01)00109-4.· 3P.T. Williams and S. Besler, “The influence of temperature and heating rate on the slow pyrolysis of biomass,” Renewable Energy， vol. 7, Mar. 1996,

27、 pp. 233-250, doi:10.1016/0960-1481(96)00006-7.· 4P.T. Williams and N. Nugranad, “Comparison of products from the pyrolysis and catalytic pyrolysis of rice husks,” Energy，vol. 25, Jun. 2000，pp. 493-513, doi: 10.1016/S0360-5442(00)00009-8.· 5Yuanzhenhong,Wuchuangzhi,Malonglong.Principlesand

28、 techniques of biomass applicationM.BeiJing:Chemical Industry Press.2005.· 6Xiushuangning,Yiweiming,Libaoming. Research Progress of Pyrolysis Kinetics of BiomassJ.Biomass Chemical Engineering, 2007，41(2)：4753.· 7Zhouling,Zhoufujun,Jiangenchen,Wangmingfeng.Pyrolysis characteristics and dyna

29、mics of different components of cotton stalkJ,Transactions of the Chinese Society of Agricultural Engineering, 2009，25(8)：220225.· 8Wangmingfeng,Jiangenchen,Zhouling.Kineticanalysisofcornstalk pyrolysisJ.Transactions of the Chinese Society of Agricultural Engineering,2009,25(2)：204207.· 9Z

30、haoming,Wuwenquan. Pyrolysis Kinetics of Rice StrawJ.Transactions of the Chinese Society of Agricultural Engineering, 2002,18(1):107- 110.· 10Ayhan D. Partly chemical analysis of liquid fraction of flash pyrolysisproducts from biomass in the presence of sodium carbonateJ. EnergyConversion and M

31、anagement，2002，43 (14)：18011809· 11Qu Wenwen, Xia Hongying,eng Jinhui ,Zhang Libo ,hang Zhengyong ,Yang Kunbin. Pyrolysis characteristics and kinetic analysis of walnut shellJ.Transactions of the Chinese Society of Agricultural Engineering, 2009，25(2)：194198.· 12Haykiri-Acma H, Yaman S, Ku

32、cukbayrak S. Effect of heating rate on thepyrolysis yields of rapeseedJ. Renewable Energy，2006，31 (6)：803810· 13Liuyaoxin,Chenwei,Lhongwei. STUDY ON BEHAVIOR AND KINETICS OF CATALYTIC PYROLYSIS FOR RICE HUSKJ.Thermal Power Generation, 1002-3364(2011)04-0053-04.· 14Jayeeta Chattopadhyay, Ch

33、ulho Kim, Raehyun Kim, et a.lThermogravimetric study on pyrolysis of biomass with Cu/Al2O3catalysts J. Journal of Industrial and EngineeringChemistry, 2009, 15(1): 72-76.· 15Wangshurong,Liaoyanfang,Luozhongyang. STUDY OF POTASSIUM CHLORIDE CATALYSIS MECHANISM ON CELLULOSE PYROLYSIS KINETICSJ. A

34、CTA ENERGIAE SOLARIS SINICA, 2005, 26(4): 452-457.· 16Xiaojun,Shenlaihong,Zhengmin. TG-FTIR study on catalytic qyrolysis of biomassJ.JOURANAL OF FUEL CHEMISTRY AND TECHNOLOGY, 2007, 35(3): 280-284.· 17LIAO Yan-fen,A Xiao-qian,UN Yong-ming. Kinetics of Pyrolysis and Metal-salt Catalized Pyrolysis of WoodJ.Chemistry and Industry of Forest Products, 2008, 28(5): 45-50)· 18Yanghaiping,Chenhanping,Dushenglei. Influence of Alkali