茂名油页岩的热解特性

东北电力大学学报第26卷第2期Journal Of Northeast Dianli UniversityVol, 26, No.2006年4月Natural Science editionApr.2006文章编号:1005-299X2006)2-0073-06茂名油页岩的热解特性柏静儒王擎胡爱娟孙佰仲李少华东北电力大学能源与机械工程学院,吉林吉林132012)摘利用热重分析仪在非等温条件下对茂名两个矿区的四种油页岩进行了热解试验研究,试验中以高纯氮气作为载气采用不同的升温速率(10、20、40、50、100℃/min)加热到900℃。粒径范围00.2mm。由失重曲线可知油页岩的热解主要集中在200~600℃的温度区间。最后根据试验数据建立了热解动力学模型利用阿累尼乌斯直接作图法求得表观活化能E和频率因子(A)等动力学参数。关键词油页岩热重分析动力学中图分类号:TQ533.6文献标识码:A油页岩是矿物机体中含有固体可燃有机质的细粒沉积岩,属于腐泥类煤由水生物等浮游植物、低等浮游动物经腐败作用和煤化作用而生成1中国是世界上油页岩资源丰富的国家之一约为2万亿吨已探明的油页岩储量为315.67亿吨仅次于美国、俄罗斯、巴西、爱沙尼亚、扎伊尔列世界第六位且主要集中于茂名、桦甸和抚顺等地的55个矿区[2。本文在非等温条件下对茂名油页岩的热解特性进行了试验研究为油页岩的干馏和直接燃烧提供了理论依据。在世界能源需求不断增长的今天油页岩作为一种替代的能源资源寻求其有效开发与经济利用的途径对于缓解能源供需矛盾推动社会的发展具有重大的现实意义。1试验部分1.1试验样品样品采用广东茂名金塘矿区一层、二层、三层、旧营矿区五层的油页岩分别记以S1、S2、S3、S4试样经曆煤机研磨后再经玛瑙硏钵研曆而成粒径范围为0~0.2mm。期间对油页岩未进行仼何筛分处理,以保证油页岩的原始特性。其工业分析及元素分析的结果如表1所示。油页岩具有含碳量低、灰份高、发热值低的特性属高灰份低热值燃料疽就地就近综合利用表1油页岩的工业分析和元素分析工业分析元素分析编Q。m/kJ/kgCm/%H/%N/%On/%S。/%1.5426.9770.590.9015,68.5710.220.63S22.1321.1073.193.5838800.7422.1668.0316.231.2试验设备和试验条件试验采用美国 Perkin elmer公司生产的 Pyris i tga热重分析仪以99.99%的高纯氮气作为载收稿日期2006-01-0874东北电力大学学报第26卷气流量80mL/min。氮气主要用来维持炉內的惰性气氛同时及时将油页岩热解产生的挥发性产物带离样品从而尽可能降低二次反应的影响样品重量控制在10mg左右。升温速率(β)分别采用10℃/min20℃/minA0℃/min50℃/min,l00℃/min热解终温为900℃。2结果和分析2.1热解特性图1~4分别给出了试样S、S2、S3、S4在不同升温速率下的热解失重(TG服微商失重(DTG油线。油页岩的热解经历了以下几个不同的阶段从室温至200℃左右为低温区的失重主要是水分的析出其中包括内部水、粘土矿物的层间水3通常在该段的质量损失被用来测量样品的含水率。这阶段的曲线表现为TG曲线下滑,DTG曲线出现个峰。在200~600℃的失重区间这一阶段的热解是油页岩热解的主要部分主要是碳氢化合物的分解导致TG曲线大的滑坡和DTG曲线上的尖峰。一般认为油母岩质裂解为油、气和半焦产物的过程是分两步进行的即油母岩质先裂解为焦油焦油再裂解获得最终的产物。但由图1~4上看油页岩两步反应不明显,其与巴基斯坦4和约旦5油页岩的热解有相似的结果。其他的学者对土耳其61和摩洛哥7油页岩进行试验得到两步反应结果。也许是与油页岩的种类有关。不管曲线的形状表现的是步反应还是两步反应油页岩热解的真实机理是相当复杂的它包含了一系列的串行和并行反应600℃以上油页岩S1、S2失重约3-5%而油页岩S3S4只发生了轻微的失重。一般认为此温度阶段是由于诸如方解石、白云石和富铁白云石等碳酸盐的裂解。同时由于CO3(碳酸盐热解产物)的存在在高温下其与剩余碳发生反应如汇+CO2→C0)产生一氧化碳也使油页岩产生一定的失重油页岩S3S在这一阶段的热解量不明显,可能是由于碳酸盐矿物质含量低的缘故。1--10C/min3-40℃/min54-50℃/min5-100℃/min1--10C/min2-20℃/min3-40 C/min4-50℃/min5-100℃/min温度/t温度/℃(a)TG曲线(b)DTG曲线图1油页岩S1在不同升温速率下的TG和DTG曲线1-10C/min2--20C/min54-50℃/min5-100℃/min8x水涩2--20C/min40℃/min4-50℃/min-1801002003004050060070080090010020030040050060070080090温度/℃温度/t第2期柏静儒等茂名油页岩的热解特性08058新兴10℃/min2-20℃/min-123-40℃min4--50C/min70-10020345060700809020温度!t温度/t(a)TG曲线b)DTG曲线图3油页岩S在不同升温速率下的TG和DTG曲线1-10C/min2-20℃/min3-40℃/in4-50℃/min24681-10 C/min2-20℃/in3-40℃/min68070010020030040050060070080090温度/℃温度/℃(a)TG曲线(b)DTG曲线图4油页岩S2在不同升温速率下的TG和DTG曲线2.2升温速率对热解的影响由DTG曲线可以看岀随升温速率提高,出现最大失重速度的温度升高,表现为DTG峰向右移动这是由于升温速率的变化带来的传热速率的改变。升温速率高使样品在特定的温度下暴露的时间短另外化学动力学控制因素也在热分解反应过程中起作用。由DTG曲线可知,升温速率高时脱一定挥发份的温度范围也较宽。油页岩热解特性参数包括热解产物初析温度(t)最大热解速度(2(表示热解发生的剧烈程度)和所对应的失重峰值温度(t2、t4)油母岩质的热解终止温廚(t),以及d)所对应的温度区间·!称为半峰宽(表示油页岩热解产物释放的集中程度)具体数值见表2图2~5给出了各参数的命名法则。升温速率对油页岩热解初始温度(t1)失重峰值温度(t2和t4)及油母岩质的热解终止温度(t3)都有影响。这是因1为升温速率不同热量至外向内传递的速度就不同升温速率直2接影响锅壁与试样、外层试样与內部试样间的传热和温度梯度,-4升温速率慢,试样有充分的时间接收热量使起始分解温度和终5止温度均降低。失重峰值温度(t2)代表了整个油页岩大分子结构7的平均稳定程度。温度越高表明体系结合越紧密在热解过程中-8R不易破坏整个网络结构反之,说明油页岩的活性较高,每一部分010020030040050060070080090076东北电力大学学报第26卷温度下升温速率越低热解越充分挥发分析岀越多余重越少。但对总失重量影响不大。却明显缩短了到达指定失重所需要的时间因而提高了失重速度完成整个热解的时间明显缩短了表2油页岩热解特性參数试样升温速率温度/℃失重速率/w%/minT⊥/℃-11.274.1245352645252456018.942.00000004395155718,093动力学分析3.1动力学模型的建立油页岩热重分析的研究主要针对失重最为剧烈的阶段对该阶段的分析计算可构建热解的表观反应动力学模型并求解主要的反应动力学参数。根据质量守恒定律试验升温速率恒定用热重法分析油页岩受热失重过程常采用下述简单的动力学模型来模拟8]k/(a)(1)式中a为温度为时的转化率其值为i=0二0其中momr、m分别为样品的初始质量温度为T时的质量、反应终止时的质量。反应速率常数h可由阿累尼乌斯( Arrhenius)定律:h=x(-B)求得对于均一的反应动力学常有如下方程:a)=(1-ay这里n为反应级数对于给RT定的升温速率P有如下关系沪=濒率因子A和活化能E以及反应级数n需通过热重曲线的分析计算求取假定不同的升温速率下求解过程中上述参数不变。将a、k以及β的值代入(1)式中然后分离变量可得第2期柏静儒等茂名油页岩的热解特性7732动力学参数的计算运用阿累尼乌斯直接作图法求解动力学参数对式(2)进行整理并将(a)的表达式代入,可得dt(1-ayβRT(3)0对(3)式两边取对数可得:叫(12)y2]=m-k(4)-30b对于油页岩的热解反应在一定的升温速率下活化能E和频率1251301351401451015160165因子A是定值而反应级数n需要试验数据来计算确定。对于上述反图6不同n值所作的Y-X曲线应根据试验TG及DTG曲线可得不同温度时对应的a及值并对n取不同的值进行试算直到所取得n值可使(4)式的左端Y=26(d对右端的函数图像接近直线即可此时的n就是24该油页岩的反应级数由直线-斜率可求出活化能由直线截距量2ln2可求出频率因子。以油页岩S在=50时的数据为例图6所6nA=5072+0.136E示对不同n值(n=0.5、1、1.5、2)所作得的Y-X曲线表明当nE(J/mol)1时的线性关系最好因此油页岩的热解可视为一级反应。动力参数计算结果见表3。图7油页岩热解动力学参数间的补偿效应可见升温速率对活化能的影响不明显。而且样品的种类、所取温表3不同升温速率下油页岩的E和A值度范围、所用的数值求解方法试验升温速率B温度范围活化能E频率因子条件等因素都会造成活化能的差试样/ki me异。据文献资料[910]活化能值393~4578.34E+11151.422.48E的变化范围也很大。但是解的活化B=40388-481147.493.18E+11能和频率因子之间存在着补偿效393-4954.07E+11应即活化能大的情况下,频率因子=100也较大这一对参数之间是线性互3.88E+09补的如图7所示对这一线性关系380~493127.99的最简单的解释可由对阿累尼乌斯125.8921.451.59E+10公式中的反应速率常数k的变形得382-462129.64到384E122.93In a= Ink+RT(5)388-4971.75E+10由于补偿效应的存在,不能仅仅用活化能或频率因子的大小来评价油页岩的热解反应性。不同的油388~496页岩有不同的补偿效应意味着各反343-4956.E+07351-4951.E+0878东北电力大学学报第26卷4结论1)茂名油页岩在200~600℃温度范围内的热解反应最激烈在此期间挥发份几乎全部析出。第三阶段是碳酸盐热解的结果茂名油页岩由于含量低的缘故此阶段不甚明显2)随着升温速率的增大热滞后现象的增加以及化学反应动力学因素使油页岩热解的初始温度失重峰值温度及热解终止温度均升高表现为DˆG峰向髙温侧移动,3)根据试验数据建立热解动力学模型得到油页岩的热解反应属于一级反应。利用阿累尼乌斯直接作图法求得表观活化能E)和频率因子(A)等动力学参数。升温速率对活化能影响不大,但活化能与频率因子间存在着明显的补偿效应。参考文献[]刘柏谦油页岩及其流化床燃烧M]长春洁林科学技术出版社,999[2]于海龙姜秀民桦甸油页岩热解特性的研究J]燃料化学学报.2001295)450-453[3] A, K Burnham. B E, Huss. M. F Singleton. Pyrolysis Kinetics for Green River Oil-shale from the Saline Zond J ] Fuel 1983 6x 1): 199[4] Paul T. Williams Nasir Ahmad Invetigation of Oil-shale Pyrolysis Processing Conditions Using Thermogravimeteic Analysis[ J ]. ApplieEnergy.20006(2):13-133[5] JOJaber S D Probert. Non-isothermal Thermogravimetry and Decomposition Kinetics of Two Jordanian Oil Shales under Different Pro-essing Condition J ]. Fuel Processing Technology. 2000 6357-70[610. Murat Do? an , B Zuhtu U ysal Non-isothermal Pyrolysis of Three Turkish Oil Shale[ J ] Fuel. 1996( 12): 1424-1428[7] Deepak S. Thakur, H. Eric Nuttal Kinetics of Pyrolysis Moroccon Oil Shale by Thermogravimetry J ]. Ind Eng, Chem. Res,Is1987,26(7)1351-13568]陈镜泓李传儒热分析及其应M].北京科学出版社,1985[9] Mustafa Versan Kk, M. Reha Pamir, Comparative Pyrolysis and Combustion Kineties of Oil Shales. Journal of Analy tical and Applied Pyrolysis J1.200055:185-194[10]M. C. Torrente, M.A. Galan. Kinetics of the Thermal Decomposition of Oil Shale from Puertollano( Spain I J ].Fuel.200180327-334The pyrolysis Characteristics of Maoming oil ShalesBaI Jing-ru , WANG Qing HU Ai-juan SUN Bai-zhong LI Shao-huaNortheast Dianli University Jilin 132012 chinaAbstract: Four oil shale samples from Maoming deposits in China have been non-isothermally pyrolysed on athermogravimetric analyser( TGA). The experimental studies were performed at different heating rate( 1020 40 50,100.C/min )using nitrogen as the purge gas. A range of particle size of 0-0.2 mm was used toheat to 900 C. The weight loss data were analysed that the main pyrolysis of oil shale took place in one region and the major weight loss occurred in the range of 200-600C. On the basis of experiment data a pyfactor( A)were obtained by the direct Arrhenius plot method arent activation energy(E)and frequencyrolysis kinetic model was proposed. The kinetic parameters ofKey words oil shale ,thermogravimetric analysis kinetics