麦秆热解机理研究

第42卷第4期vo.42,No.42012年4月JOURNAL OF UNIVERSITY OF SCIENCE AND TECHNOLOGY OF CHINAApr.2012文章编号:02532778(2012)04031807麦秆热解机理研究李小民,李永平,邓权威,林其钊中国科学技术大学热科学与能派工程系安合屁230026)摘要:采用热重法对不同加热速率、氮气气氛条件下的麦秆失重行为进行了研究研究发现,麦秆失重主要由水分析出、有机物热解和无机物分解等3个阶段组成,其中在无机物分解阶段还伴随有轻微的炭化;此外,随着加热速率增加,有杋物分解的初始温度、終了溫度、峰值失重速率对应溫度等相应增加,利用双外推法推断了麦秆中有杌物热解的最概然机理函数,发现麦秆热解可以用二级反应机理进行模化,同时,确定了麦秆热解的活化能、指前因子等反应动力学参数.对利用双外推法得到的模型计算结果与实验结果进行了对比,结果表明,两者之间符合较妤,但仍存在一定偏差.对引起偏差的原因进行了分析关键词:生物质;热解;热重分析;双外推法中图分类号:TK6文献标识码:Adoi:10.3969/.isn.0253-2778.2012.04.009引用格式:李小民李永平邓权威,等,麦秆热解机理研究[门].中国科学技术大学学报,2012,42(4):318324Li Xiaomin, Li Yongping, Deng Quanwei, et al Research on mechanisms of pyrolysis of wheat straw[]. Journal of University of Science and Technology of China, 2012, 42(4): 318-324.Research on mechanisms of pyrolysis of wheat strawLI Xiaomin, LI Yongping, DENG Quanwei, LIN QizhaoCUhniuersity of Sciene and Technology of China. De partment of Thermal Science and Energy Engineering, Hefei 230026, China)Abstract: a thermogravimetric experiment was performed to investigate the mass loss of wheat straw ininertial ambience of nitrogen. The results show that, in the process of mass loss, there exist three stages:preheat drying, pyrolysis and thermal decomposition of inorganic components, accompanying mildcarbonization in the third stage. At the same time, the characteristic parameters, such as onsettemperature, final temperature and peak temperature, increase along with the increase of heating rate.Furthermore, the dual extrapolation method was utilized to obtain appropriate kinetic mechanisms tomodel the pyrolysis of wheat straw. It was concluded that the mechanisms of pyrolysis for wheat strawfollow second-order reaction. In addition, the parameters of chemical kinetics, such as activation energyand frequency factor, were obtained as well. The preliminary comparisons between the model and thexperimental data validate the reliability of the kinetic mechanisms as a whole. Finally, the causes of thediscrepancy were analyzed.Key words: biomass; pyrolysis; thermo-gravimetric analyse; dual extrapolation method收稿日期:201107-18修回日期:2011-10-14基金项目:国家重点基础研究发展(973)计划(2010CB22730)资助YHS中国煤化工作者简介:李小民,男,1971年生博士生,研究方向:秸秆生物质利用CNMHG通讯作者林其钊博上/教授,Emai:qlin@ use edu en第4期麦秆热解机理研究0引言1实验方法热解是生物质热化学转换的初始反应,对热1.1实验原料化学转换装置的效率、运行水平等具有重要影本文利用热重法对麦秆热解进行了研究.为了响口.热重法是生物质热解研究的常用方法14.消除颗粒粒径对反应的影响,采用的试样粒度为刘乃安等2利用热重法对林木的热解进行了研100~120pm麦秆试验样品的主要性质参数如表1究,认为可以使用二级反应动力学模型对其进行所示模化. Nassar等{3利用热重法研究了稻壳热解,使1.2实验仪器与实验条件用的是准一级反应模型.张晓东等利用热重法使用日本岛津的热分析装置DTG60H得到了研究了玉米秸秆的热解,采用的反应机理为一级麦秆样品的热重分析(TGA)曲线和微分热重分析反应,孔晓英等利用热重法对麦秆和蔗渣的热(DTG)曲线;采用DSC60得到了差热分析(DTA)解进行了比较研究,分析了加热速率对麦秆和蔗线,参比物为aAl2O為.实验时加热炉由室温加热渣热解特性的影响. Vuthaluru利用热重法对麦到1273K,升温速率分别为5,10,20,40K/min,加秆/煤混合物的热解特性进行了研究,分析了不同热炉通入的氮气流速为50mL/min,每次实验使用比例混合物的热解特性众所周知,在热解动力学的样品约为10mg为了消除系统误差,每个实验结研究中,由于反应活化能和指前因子的补偿效应,束后作一个相同条件下的空白实验使很多机理函数都能有一个良好的线性结果81这使研究得到的同一物质的动力学参数出入很2结果与讨论大. Vuthaluru!°在麦秆热解研究中,当反应级数取2.1麦秆热解失重行为0.8时,计算得到的活化能为80.2kJ/mol,指前因实验得到的4种加热速率条件下的热重分析曲子为1.6×10min-;而当反应级数取0.52时,计线如图1所示,微分热重曲线如图2所示,差热分析算得到的活化能为1148kJ/mol,指前因子为曲线如图3所示4.7×10°min-.这是由于选择的机理函数与实际从图1、图2和图3可以看出,麦秆失重主要经发生的动力学过程之间存在差异造成的,故选择最历3个阶段首先是水分析出阶段,主要是麦秆中自概然机理函数非常重要.最概然机理函数的推断方法很多,如双外推法、 Satava法、 bagchi法、 Malek法等.王明峰等利用 Malek法推断了玉米秸秆解的机理函数,认为玉米秸秆热解服从随机成核K/随后生长机理任庚坡等0利用双外推法对元宝山褐煤的热解机理进行了推断,认为其热解服从 Anti-Jander的三维扩散机理.本文首先利用热重法,对麦秆热解进行了实验研究然后利用双外推法,确定了麦秆中有机物热解的最概然机理函数最后,将模型计算结果与实验结果进行了对比,验证了模型结温度/K果的可靠性,并对造成二者之间存在的偏差的原因图1热重(TG曲线进行了分析Fig 1 Weight loss curve from thermogravimetry analysi表1麦秆元素分析与工业分析Tab. 1 Ultimate analysis and proximate analysis of wheat straw质量分数/%低热值/(kJ·kg1)EalSaf FCar VarTHE中国煤化工CNMHG中国科学技术大学学报第42卷应增加.这表明提高加热速率也会导致挥发份析出滞后,原因也是材料内部传热需要一定时间.此外,0.05如图2所示,当加热速率为10K/min和20K/minv-0.15时,DTG曲线上出现了台阶,这是由于麦秆中纤维-0.20素的热解初始温度较高引起的麦秆是由半纤维素、纤维素和木质素组成的混合物.Yang等1)曾对各组分的独立热解过程进行了研究结果表明,木质素0.450.50和半纤维素的热解过程为放热过程,其中前者的热055400600解初始温度约为160℃,后者约为220℃;纤维素的温度/K热解为吸热过程,且其热解初始温度最高,约为图2微分热重DG)曲线315℃.根据图3,当加热速率为5K/min和10K/Fig 2 Derivative thermogravimetry curve of wheat strawmin时,热解过程总体表现为吸热过程,这表明半纤维素和木质素热解放出的热量不足以维持纤维素热解所需热量,而须补充外加热量以维持热解的进行第三阶段是无机物的分解,并伴随有轻微的炭化过程根据廖翠平等的研究,麦秆中的无机元素主要是K,CaMg,Al,由这些元素的组成的无机物分解为吸热反应由图3的DTA曲线可以看出,当加热速率为20和40K/min时,放热分别在707K和746K达到峰值,表明此时木质素热解已趋结束,而无机物分解开始占据优势.图2的DTG曲线中的4000001200小失重峰就对应着某种无机物的分解,如CaCO3在温度/K950K左右的分解.此外,该阶段还存在轻微的炭化图3差热分析(DTA)曲线过程.由图3的DTA曲线可以看出,当加热速率分Fig 3 DtA curve of wheat straw at different heating rates别为20和40K/min时,仍存在放热峰.考虑到麦由水、结晶水和吸附水的析出.水分析出阶段的主要秆中的无机物分解是吸热反应,故放热峰的存在特征参数如表2所示,随着加热速率的增加,水分析定是由某种其他物质发生放热反应造成的.根据出的终了温度T1、峰值失重速率对应温度Tm均相Yang等的研究结果,木质素热解为放热反应,热应增加这表明提高加热速率将引起水分析出滞后,解温度范围较宽(约为160~900℃),因此该阶段这是由于材料内部传热需要一定时间造成的.其次的放热是由木质素热解造成的进一步,结合图2的是麦秆中的有机物热解阶段,其失重曲线的主要特DTG曲线,该阶段的失重速率很小,在DTG曲线上征参数如表2所示,随着加热速率的增加,挥发份析近似为一条直线,表明在无机物分解阶段尽管存在出的终了温度T、峰值失重速率对应温度Tmx均相炭化,但已经不明显表2麦秆热解失重曲线的特征参数Tab. 2 The characteristic parameters of the TG-DTG curves of pyrolysis of wheat straw热解步骤B(K·min-1初始温度终了温度峰值失重速率对应温度峰值失重速率失重率Ti/K T/KTmax/K/(s-1)/%0.0149水分析出200.02654.800.04404.27中国煤化工259有机物热解CNMHG388第4期麦秆热解杌理研究2.2热解反应机理研究最概然机理函数推断有双外推法、 Satava法、2.2.1双外推法原理Bgch法Maek法等,其中,双外推法是指将加多相反应的动力学方程式为热速率和转化率a双双外推为零以得到样品在热da/dt= k(T)f(a).(1)平衡状态下的E*及原始状态下的E~然后根据式中,a为转化率其定义为a=[m-m(t)]/(m-两者之间的差异大小而确定最概然机理函数其基m),mo为试样初始时刻的质量,m(t)为试样在t本思想是认为固体样品在一定加热速率的温度场中时刻时的质量,m为试样热解完成后的最终质量;t的受热过程是非定温过程样品本身的热传导造成为反应时间;k(T)为反应速率常数;f(a)为反应机了样品本身以及样品与温度场之间处于非热平衡状理函数确定速率常数k(T),需要选择合适的f(a)态,故在此基础上得到的反应机理及动力学参数与形式.多相反应动力学模型的主要研究方法如表3真实情况有偏离:加热速率越大,偏离就越大显然,所示(13如果将加热速率外推到零,就可以获得样品处于热表3非均相反应动力学模型研究方法平衡状态时的动力学参数,它将反映过程的真实状Tab. 3 The model methods of kinetics况同时,当样品在不同转化率时,其表观活化能也on heterogeneous reaction将呈现规律性的变化.因此,如果获得转化率为零时序号其本思想开究方法的动力学参数,就可认为是体系处于原始状态时的I基于 Arrhenius定律的表达式微分法参数0.双外推法的步骤是:积分法首先,由 Cotas-Redfern积分式:使用其他形式的k(T表达式HE公式:k(T)=BTBH公式:Ce"rn[G(a)/T]=ln(AR/配)一E/(RT).(3)直接使用式(1)NPK方法固定加热速率由ln[G(a)/T]与1T之间的直线关系,求出反应的表观活化能和指前因子.在某固表3中方法Ⅱ和方法Ⅲ的优点是其拟合公式与定加热速率下,选择多个线性关系极佳的G(a),由数据符合较好,但诸如B、m、C、D等参数却缺乏明式(3)确定动力学参数进而由下述两式外推而分别确的物理意义而方法I中,反应活化能E和指前得到E和Ap:因子A等动力学参数的物理意义却分别有碰撞理E=a+b1B+c1P+d1P,E-0=a1,(4论和活化络合物理论的诠释,故在此使用方法InA=a2+b2+c2p+d?,(lnA)p-=a2进行研究从数学处理方法上来看,方法I可分为积(5)分法和微分法较之微分法受限于实验数据噪卢的然后,根据 Ozawa积分式:影响积分法更具优势1.由式(1)积分得到lgB=lgAE/R·G(a)]-2.315-0.4567E/(Rr),G(a)=T(2)当转化率a固定,则G(a)一定这样由g8与1/T式中,G(a)为机理函数;B=dT/dt为加热速率;A为的直线关系,求出一定a时的表观活化能然后按照指前因子;E为活化能;R为通用气体常数从式(2)下述方程将a外推到零而得到无任何副反应干扰、可以看出,其右端的积分部分没有解析表达式,因此体系处于原始状态下的E:积分法的难点在于其温度积分的求解.温度积分式E=a3+b2a+c3a2+d3a3,E0=a3.(7)的求解主要有解析法和数值方法,其中解析法的使将选定的几个G(a)式的E与E-相比较,其中用较多,如 Coats- Redfern, Flynn-Wall-Ozawa等".相同或相近者即为该反应的最概然机理函数众所周知,由于反应活化能和指前因子的补偿效应,2.2.2推断过程与结果使很多机理函数都可能呈现良好的线性结果双外推法的第一次外推是根据式(3)进行的要这使得在相同实验条件下,不同研究者得到的同一求所选择的机理函数G(a)使ln[G(a)/7]~1/T物质的反应动力学参数出入很大显然,这是由于选之间的关系曲线具有极佳的直线性目前的异相反择的机理两数形式与实际发生的动力学过程之间的应机中国煤化工相界面模型扩散差异造成的,故选择最概然机理函数对增强模型机理的通用性显得尤其重要模型xCNMHG于试算时后3种反应级数模型,如中国科学技术大学学报第42卷表4所示能E比较发现麦秆热解的机理函数为二级反衰4采用的机理函数应,相应的动力学参数为E=1.17×10°J/mol和Tah. 4 The several mechanisms adopted in this paperA=1.27×103s反应级数nf(a)2.3结果可靠性验证与误差分析利用双外推法,得到麦秆中有机物热解为二级n=2,3,4(1-(1-a)-)/(1-n)(1-c)n反应机理活化能E为1.17×105J/mol指前因子根据实验得到的加热速率分别为5,10,20,A为1.27×108s-1.为了验证该模型的可靠性,将40K/min的热重曲线,由式(3)得到的动力学参数,其与4种加热速率下的实验结果进行了比较.二级反应的机理函数为即活化能E和指前因子A如表5所示同时表5也G(a)=(1-a)1-1给出了由公式(4)和(5)外推得到的E和A-,该这样由式(3)得到的麦秆热解理论计算式为次外推在表5中由“外推1”标记另一个外推是根据 Ozawa公式和式(7)进行的首先根据一定的转化率a得到不同加热速率下e邓p(-的体系温度,进而根据式(6)得到该转化率时的活化式中,E=117×10°J/mol;A=1.27×108s';R=能计算中所选定的转化率,一定转化率、不同加热8.314J/(molK);B/(K·s-1)为加热速率.图4速率时的体系温度,以及根据 Ozawa公式计算得到给出了按照式(9)得到的理论曲线与实验数据的比的活化能如表6所示.然后将得到的活化能根据式较,其中加热速率B分别为5,10,20和40K/(7)外推得到E.最终得到E0=111236.12从图4可以看出,理论曲线与实验数据符合较J/moL将E与表4给出的一次外推得到的活化好.但二者之间仍有一定偏差并且不同加热速率下表5 otas-Redfen公式得到的对应于不同机理函数的动力学参数Tab. 5 The parameters of kinetics from Cotas Redfern formula corresponding to several mechanisms反应级数np(K·min-1)E/(kJ·m-1)A/(s-1)E/(kJ·ml-1)A/(s1)E/(kJ·mol-1)A/(s1)E/(kJ·mol-1)A/(s-1)8.86E+041.66E+051.11E+055,62E+071,40E+057.42E+101.71E+051.93E+8.04E+042.24E+041.1E+054.49E+071.51E+056.05E+111.93E+051.71E+168.09E+043.05E+041.14E+058.70E+071.57E+051.82E+122.03E+058.20E+168.68E+041.40E+051.19E+052.67E+081.60E+053.36E+122.05E+058,79E+外推11.05E+071.14E+071.17E+051.27E+081.21E+051.95E+091.32E+054.07E+10表60wa公式得到的不同转化率时的活化能Tab, 6 The parameter of activation energy from Ozawa formula at different conversion ratio加热速率/(K·min-1)E/(kJ·mol-1)KT528.73555.61571.161.17E+050.20576.26592.051.21E+05578.43593.43608.811.23E+05576.2422.910.50601.26617.68634.07124E+050.6059277124E+05600.18618.06635.411.22E+05609.72中国煤化工汁5648.29666.HaCNMHG第4期麦秆热解机理研究B: 5 K/minB: 10 K/'min0.88哥测量测量值理论曲线理论曲线00450500550600650700750450500550600650700750温度/K温度/Kp: 20 K/minB:40▲测最值测量值理论曲线理论曲线50500550600650700750550600650700750温度/K温度/K图4理论曲线与实验数据的对照Fig 4 The preliminary comparisons between the model and the experimental data产生偏差的大小与位置也不相同.现对此进行简单分析出、有机物热解和无机物分解等阶段组成,其中分析:①热重实验中影响麦秆热解的因素既有动力在无机物分解阶段还伴随有轻微的炭化过程学方面的作用,也有传热、传动量、传质等方面的作(Ⅱ)根据得到的理论计算式与实验结果的对用,单纯动力学模型研究尤法对其作出准确描述结比情况,表明采用二级反应机理模化麦秆中的有机果造成计算值与测量值之间的偏差.②麦秆是由半物热解具有一定的可靠性纤维素、纤维素以及木质素组成的混合物.首先,实(Ⅲ)麦秆中有机物热解的活化能为1.17×105际热解过程中随着加热速率发生变化,由于材料内J/mol,指前因子为1.27×108s部传热需要一定时间,导致各组分的热解初始温度、终了温度、峰值温度等均会发生变化;其次,各组分参考文献( References)的析出产物也不同如半纤维素热解主要析出CO、[1] Vuthaluru H B. Thermal behavior of coal/biomass纤维素主要析出CO、木质素热解主要析出CH4和blends during co-pyrolysis [J ]. Fuel ProcessingH2u1],同时,挥发份中各组分的扩散系数也不相Technology,2003,85:141-155.同;第三热解过程中半纤维素等3种组分之间还会[2] Liu Naian,, Wang haihui, Xia dunhuang,ctal存在相互影响结果使测量值与计算值之间出现的Modelling the thermal decomposition of wood and偏差大小与位置也不相同.leaves [J]. Journal of University of Science andT3结论中国煤化工解动力学模型研究CNMHG81):4048(I)麦秆在惰性气氛下的失重过程主要由水[3] nassatmn. Thermal analysis of kinetics of324中国科学技术大学学报第42卷bagasse and rice straws]. Energy Sources, 1999, 21: [10] Ren Gengpo, Zhang Chaoqun, Wei Lihong, et al.131-137Pyrolysis characteristics of super-fine pulverized lignite[4] Zhang Xiaodong, Xu Min, Sun rongfeng, et al. Studyand its pyrolysis mechanism []. Journal of Engineeringon the kinetics of corn stalk pyrolysis by TG-DTGfor Thermal Energy and Power, 2005, 20(4):analysis [J]. Jourmal of Fuel Chemistry and407-410Technology,2006,31(1):123-125任庚坡张超群魏砾宏,等超细褐煤粉的热解特性张晓东,许敏,孙荣峰等.玉米秸热解动力学研究[及其热解机理[J].热能动力工程,2005,20(4):407燃料化学学报,2006,31(1):123-125.[5] Kong Xiaoying, Ma Longlong, Wu Chuangzhi, et aL. [11] Yang Haiping, Yan Rong, Chen Hanping, et al.Comparsion of thermogravimetry characteristicsCharacteristics of hemicellulose, cellulose and ligninwheat and basage [J]. Renewable Energypyrolysis[J.Fuel,2007,86;1781-1788112(6):12-14.[12] Liao Cuiping, Wu Chuangzhi, Yanyongjie, et al,孔晓英,马隆龙,吴创之,等.麦草与蔗渣基本热解特性Chemical elemental characteristics of biomass fuels in的比较[],可再生能源,2003,112(6):12-14China [J]. Biomass and Bioenergy, 2004, 27: 119-[6] Vuthaluru H B. Thermal behaviour of coal/biomass130blends during co-pyrolysis [J]. Fuel Process[13] Chen Haixiang, Liu Naian. Approximation expressionsTechnology,2003,85:141-155for the temperature integral]. Progress in Chemist[7] Wang Mingfeng, Jiang Enchen, Zhou Ling. Kinetic2008,20(7/8):1015-1020analysis of corn stalk pyrolysis [J]. Transactions of the陈海翔,刘乃安温度积分近似式研究[J化学进展,CSAE,2009,25(2):204207.2008,20(7/8):1015-1020王明峰蒋恩臣,周岭.玉米秸杆热解动力学分析[].[14] Ran jingyu, Zeng Yan, Zhang li, et al Pyrolysis and农业工程学报,2009,25(2):204-207kinetics characteristics of typical crop biomass [J].[8] Liu Naian, Wang Binghong, Fan Weicheng. KineticJournal of Chongqing University, 2009, 32(1):76-compensation effect in biomass thermal decompositionJ]. Fire Safety Science, 2002, 11(2): 63-69.冉景煜,曾艳张力,等.几种典型农作物生物质的热解[9]胡荣祖,史启祯热分析动力学LM北京:科学出版及动力学特性[]重庆大学学报,2009,32(1):7681社,2001:113-115中国煤化工CNMHG