基于差示扫描量热技术的生物质热解两步连续反应模型研究
- 格式:pdf
- 大小:236.50 KB
- 文档页数:5
物理化学学报(WuliHuaxueXuebao)ActaPhys.鄄Chim.Sin.,2006,22(7):786耀790
基于差示扫描量热技术的生物质热解两步连续反应模型研究陈海翔刘乃安*范维澄(中国科学技术大学火灾科学国家重点实验室,合肥230026)
摘要采用差示扫描量热分析仪对我国的一种生物质试样在空气气氛中进行了实验,发现试样从常温到923K高温的低速升温过程中,经历了两步明显的放热过程.对放热机理的分析表明,第一步主要是由半纤维素和纤维素的有氧热解过程控制,第二步放热过程则受木质素热解和炭的氧化反应的共同作用.采用等转化率方法和优化计算方法,对热解过程的动力学模型进行了研究,结果表明,两步连续反应机理可用于描述生物质在空气气氛中热解的放热动力学.
关键词:生物质,热解,动力学模型,差示扫描量热中图分类号:O643
Two鄄stepConsecutiveReactionModelofBiomassThermalDecompositionbyDSC
CHEN,Hai鄄XiangLIU,Nai鄄An*FAN,Wei鄄Cheng(StateKeyLaboratoryofFireScience,UniversityofScienceandTechnologyofChina,Hefei230026,P.R.China)
AbstractThethermaldecompositionofonekindofbiomassinairhasbeeninvestigatedbydifferentialscanningcalorimetry(DSC)analyzer.Theresultsindicatethattheheatingprocessofsamplesfromambienttemperatureto923Katlowheatingratesshowtwoobviousexothermicpeaks.Accordingtothedecompositionmechanism,thefirstexothermicstepisattributedtooxidativedegradationofhemicelluloseandcellulose,andthesecondexothermicstepisattributedtolignindegradationandcharoxidation.Thereactionmodelofthestudiedbiomassthermaldecompositionhasbeenstudiedbyiso鄄conversionalmethodsandoptimizationcomputation.Theresultssuggestthatthetwo鄄stepconsecutivereactionmodelissuitabletodescribetheexothermofbiomassthermaldecompositioninair.
Keywords:Biomass,Thermaldecomposition,Kineticmodel,DSC
生物质的热解是火灾安全科学领域的重要研究课题.据统计,我国的森林资源人均少于1.1伊103m2,然而1950到1990年间森林火灾年均毁掉8.2伊109m2[1].为了有效减少森林火灾造成的损失,我国发展了林火蔓延模型和防火隔离带等技术.这些技术需要生物质热解机理研究的支撑.深入研究生物质热解的动力学规律,是理解生物质着火燃烧过程的基础.基于生物质热解动力学的研究,可以进一步完善林火蔓延的物理模型,发展筛选防火树种的新技术.同时,生物质的热解也是化学工程科学领域的重要研究课题.用生物质通过热解制取生物燃油前景广阔,是解决我国石油资源不足的有效途径.近几十年来,人们广泛研究了生物质热解过程,特别是通过热重分析(TG)建立了各种生物质热解失重过程的动力学模型.对文献的总结可以看出,生物质表观热解失重动力学模型大致可以分为两类:一类模型将生物质看成单一组分,使用单一组分在全局温度区间的固鄄气反应动力学来描述其热解失重
[Article]www.whxb.pku.edu.cnReceived:November21,2005;Revised:February21,2006.*Correspondent,E鄄mail:liunai@ustc.edu.cn;Tel:+86551鄄3601668.国家自然科学基金(50323005,50576090),973项目(2001CB409600),教育部新世纪优秀人才支持计划,安徽省优秀青年科技基金(2004鄄2005)和中国科学技术大学青年基金资助鬁EditorialofficeofActaPhysico鄄ChimicaSinica
July786No.7刘乃安等:基于差示扫描量热技术的生物质热解两步连续反应模型研究
过程[2鄄3];另一类模型将生物质失重过程看成是它的各种赝组分在全局温度区间上分别热解过程的加权叠加[4鄄6].Liu等[7鄄8]分析了它们的缺陷,针对生物质DTG曲线的双峰特点,发展了一种“双组分分阶段反应模型”,对两阶段分别进行建模.基于TG技术的生物质热解机理的研究仍在蓬勃发展,而差示扫描量热法(differentialscanningcalorimetry,简称DSC)作为另一类热分析技术,在生物质热解过程的研究中却往往只是作为TG分析的一种补充,DSC图谱也只作定性的说明.本研究将运用DSC技术,对我国的一种经济树种进行热解机理研究,探讨空气气氛中生物质热解的基本过程,建立基于DSC图谱的生物质热解动力学模型.1实验实验采用德国耐驰公司的STA409C热分析仪,可以同时记录生物质热解过程中的TG和DSC图谱.实验试样是安徽祁门山区的油茶树.该树种在安徽舒城等地广泛种植,有很高的经济价值.研究它的热解特性,可以对其防火性能进行一些评估.实验时,将干燥后的试样用粉碎机研磨后,用筛子选取粒径小于300滋m的颗粒.选用10mg左右的样品用于热分析实验.加热炉由室温加热到923K,高于这个温度样品的质量变化不再明显.升温速率控制为5、10、15、25K·min-1,通入的空气流速控制为42mL·min-1.热分析仪可以记录升温过程中试样的TG和DSC图谱.每种条件下的实验进行两次,以保证可重现性.2结果与分析2.1DSC曲线的定性分析油茶树试样在热分析仪中热解过程的DSC曲线见图1.很明显,图中DSC曲线呈现了两个放热峰.对应于5、10、15、25K·min-1的升温速率,第一个放热峰温度分别为590.6、600.3、615.6、618.2K;第二个放热峰温度分别为693.3、700.7、714.5、732.6K.这两个放热峰可以与TG实验分析对比来加以解释.在TG实验中,生物质在563~613K和713~773K温度区间会呈现较大的质量损失过程,在DTG曲线上表现为两个峰值[7鄄8].第一步失重过程主要是半纤维素和纤维素的有氧热解过程控制,第二步失重过程则受木质素热解和炭的氧化反应的共同作用.这两个失重过程会伴随产生热流变化.DSC分析仪正是记录了此热流信号,在DSC曲线上就表现为两个放热峰.图1还表明随着升温速率的提高,伴随热解过程的热量流率也增大,而且第二峰逐渐比第一峰尖锐.这说明升温速率的增加有利于炭的氧化反应的发生.升温速率增大,生物质在低温段停留的时间短,这不利于挥发气体的产生而有利于炭的生成,从而使高温区炭的氧化反应比较剧烈,释放更多热量.2.2热解动力学分析方法
通常单步固体热降解过程可以表述为S(solid)寅R(solid)+V(gas)
有下面的反应动力学方程:d琢f(琢)=A茁exp(-Ea/RT)dT(1)
其中琢是固体热解转化率,茁是升温速率(K·min-1),Ea是活化能(kJ·mol-1),A是指前因子(s-1),R是气体
常数(8.314J·mol-1·K-1),T是绝对温度,f(琢)是反应
的机理函数.基于上式发展了一系列分析方法,其中等转化率方法可以在不知道机理函数的情形下计算活化能,被广泛使用.这类方法的代表是(I)Friedman法[9]:
ln(茁d琢dT)=ln[Af(琢)]-EaRT(2)
对不同升温速率的曲线,在等转化率处,由ln(茁d琢/dT)对1/T作图,用最小二乘法拟合数据,由斜率可
以求得Ea
.
(II)OFW法[10鄄11]:
ln茁艿-5.333-ln[RAEa琢0乙1f(琢)d琢]-1.052EaRT(3)
对不同升温速率的曲线,在等转化率处,由ln茁对1/T
图1油茶树试样在空气中不同升温速率下热解过程的DSC曲线
Fig.1DSCcurvesofthedecompositionprocessofoil鄄teawoodsamplesinairatdifferentheatingrates
787ActaPhys.鄄Chim.Sin.(WuliHuaxueXuebao),2006Vol.22作图,用最小二乘法拟合数据,由斜率可以求得Ea.本研究使用这两种方法对全局温度区间的DSC曲线进行分析,计算活化能.为此,先从原始数据中截取453~853K之间,每隔0.5K的数据点.起始温度选取453K是考虑在此温度下样品的水分可以完全蒸发,而热解反应几乎没有进行.终止温度选取853K是因为从图1中可以发现此时的DSC曲线已经开始呈水平状,可以认为反应已经结束.热解的转化率由下式计算:琢=tti乙[DSC(t)-Baseline(t)]dttfti乙[DSC(t)-Baseline(t)]dt(4)其中ti是反应起始的时刻,tf是反应结束的时刻.由于DSC信号受仪器基点漂移的影响比较严重,因此减去直线型的基线Baseline进行校准[12].样品热解过程的转化率曲线及其变化率曲线见图2.图3是等转化率方法计算得到的活化能对转化率的曲线图.Friedman法计算的活化能在琢<0.2时几乎保持在120kJ·mol-1,然后减小到94.5kJ·mol-1(琢抑0.35),再很快增加.而OFW法计算的活化能先减小到约110kJ·mol-1(琢抑0.41),然后增加.两种计算结果有一定的差别,但是均显示样品的热解反应比较复杂,存在多步的反应过程.而且,都可以在琢抑0.4附近把活化能曲线分为前后两部分,这是
DSC曲线存在两个峰的必然结果,0.4是第一个DSC峰下面积占总面积的比值的近似值.在下面的
热解模型研究中,将据此采取两步反应的热解模式.2.3热解动力学模型研究
引言中提到研究者建立了各种生物质热解失重过程的动力学模型,而几乎没有见到对热解过程中热流变化进行建模的研究.本研究提出两步连续反应的机理模式用于DSC模拟.单组分和多组分全局区间的热解失重动力学模型没有考虑热解过程是连续的,存在重要的中间产物炭.由于DSC两峰交叠比较明显,分阶段模型存在不易确定分界点的问题.两步连续反应模型虽然存在动力学优化计算的难题,但揭示了热解反应过程,在聚合物热解中得到了很多运用[13鄄14].