油页岩主要由藻类低等浮游生物经腐化作用和煤化作用而生成,是一种固体可燃矿产[1-2].由于石油资源日益紧缺,油页岩的开发利用已愈来愈受到重视[3-4].在世界范围内,油页岩利用途径主要包括提炼页岩油、直接燃烧和发电等,利用过程中均会导致粉尘产生.同煤尘一样,也存在火灾爆炸危险.
油页岩开发利用在我国起步较晚,有关其着火、爆炸研究的文献较少.在着火理论方面,王国金等[5]采用化学反应动力学及传热学理论建立了油页岩颗粒在流化床中燃烧的非稳态渐进扩散模型,并开发了数值求解方法.姜秀民等[6-7]利用热重分析测试仪、热显微镜和热天平等分析得到了氧浓度、粒径以及升温速率对油页岩着火的影响,依据传热学理论建立了油页岩颗粒在流化床中低温段均相燃烧以及高温段多相燃烧数学模型.
据统计,德国近1/3的粉尘爆炸事故是由于设备热表面温度过高而导致堆积的粉尘层着火[8].本文首先建立粉尘层着火的理论模型,并利用专用测试装置,测试研究油页岩粉尘层的着火规律,然后利用测试结果计算得到粉尘层燃烧动力学参数值,进而得到相应厚度粉尘层的临界着火温度,验证模型对粉尘层着火问题的适用性,以便于生产过程中预测和评估堆积油页岩粉尘的着火危险性,为安全开发利用这种资源奠定基础.
1 粉尘层着火的理论模型生产过程中粉尘层着火一般由于设备热表面温度过高而导致,且粉尘径向面积较大,而厚度较薄,因此建模过程中将粉尘层视为无限大平板,忽略粉尘层径向的传热,而仅考虑粉尘沿厚度方向的能量交换.粉尘层底部受恒温热源加热,热量以热传导的方式垂直向上传递至粉尘层,粉尘层顶部只考虑自然对流换热.
假设粉尘层着火前反应物质量无损失,反应为遵从Arrhenius定律的一级反应.由于自燃过程比较缓慢,可以假设温度在粉尘层内部的分布是稳态的,依Fourier导热定律[9]粉尘在恒温热板上的一维稳态能量守恒方程为
(1) |
式中:λ为导热系数,依文献[10],4种油页岩样品统一取导热系数为0.2 W/(m·K);A为指前因子,s-1;ΔHR为反应热,J/kg;ρ为粉尘堆积密度,kg/m3;E为反应活化能,J/mol;R为普适气体常数,8.314 J/(mol·K);T为粉尘温度,K.
设TA为环境温度,则粉尘温度可表示为T=TA+ΔT,其中ΔT为粉尘层与环境间的温差.引入无量纲温度变量θ,设
(2) |
其中,δ为无量纲加热速率:
(3) |
由Thomas边界条件[11],粉尘层内部温度先增加后降低,在粉尘层底部至半高处的某点处达到最大值,公式表示为毕奥数:Bi=hr/λ,其中h为对流换热系数,
设TS为粉尘层与环境接触表面的温度,则由Thomas边界条件:
(4a) |
(4b) |
(4c) |
对于边界条件(4c),粉尘层顶部温度的变化由热板向上导热引起,即
(5) |
将以上各方程联立可得无量纲加热速率:
(6) |
由式(2)知,δ最大时的温度为粉尘层理论临界着火温度.因此,将以上各式联立计算,对于非线性方程,利用C语言编程,采用Newton迭代法[13]计算,多次迭代至|xk+1-xk| < ε(本文取ε=10-5),可以求得不同r下的各参数值.
2 油页岩粉尘层最低着火温度实验研究 2.1 实验样品实验样品分别取自辽宁抚顺、吉林桦甸、山东龙口及广东茂名四个油页岩主要产区现场,分别记为FS, HD, LK, MM.实验时先将块状油页岩原矿石破碎,再利用球磨机粉碎研磨,并在105 ℃的烘干箱中烘干1 h,除去样品中的水分,然后利用200目筛子进行筛分,取筛下物作为实验样品.实验样品工业分析、元素分析以及粒径分析结果如表 1所示.
根据国家标准GB/T16430-1996和ASTM2021-09,测试装置如图 1所示.
盛粉环置于热板上,用于盛装待测样品,分别采用高度为2,5,12.5和20 mm,直径100 mm环体进行测试.温控模块确保热板在测试过程中保持恒定的预设温度,测温热电偶以及温度采集模块用于测试并记录粉尘层在受热过程中温度的变化.
2.3 实验结果及分析4种油页岩粉尘的测试结果如图 2所示,随粉尘层厚度的增大,最低着火温度逐渐降低;相同厚度不同产地油页岩粉尘层的最低着火温度不同,由高到低依次为FS>HD>LK>MM.依牛顿冷却定律,粉尘层顶部对流散热系数与粉尘层厚度有关,粉尘层厚度越小,对流换热系数越大,从而散热率越大,粉尘层着火所需时间和温度也就越大.导致各产区油页岩不同着火温度的原因很复杂,但油页岩成分应该是主要因素,对比四者的挥发分含量(表 1)高低依次为LK>MM>HD>FS,与粉尘层最低着火温度的变化规律基本相反,即挥发分含量越高的油页岩粉尘层的最低着火温度越低.这一规律对于MM和LK略有不同,MM油页岩含硫量大于LK,而物质中硫含量的增加会降低着火温度.
由式(6)可知,
将以上计算结果代入式(6),利用C语言计算可以得到理论临界着火温度TP数据.油页岩粉尘着火温度理论分析结果与实验值对比见表 3.
由表 3可见,实验值与理论值随粉尘层厚度的变化规律是一致的,由于理论分析忽略了粉尘层径向的传热,导致预测结果略低于实测数据,理论值与实验值误差在10%以内,该理论模型可用于分析工业生产过程中不同厚度粉尘层的着火问题.
4 结论本文利用化学反应动力学、传热学以及Thomas不对称热自燃理论建立了油页岩粉尘层着火的理论模型,并利用专用装置测试了2,5,12.5和20 mm厚度条件下四种油页岩粉尘样品的最低着火温度和变化规律,主要结论如下:
1)不同地区油页岩最低着火温度随挥发分以及硫分的升高而降低;同一地区油页岩粉尘层最低着火温度随厚度的增加而降低.
2)利用建立的油页岩粉尘层临界着火温度理论模型,求解得到了油页岩粉尘层燃烧表观活化能.
3)理论着火温度与实验测试结果相比误差在10%以内,为预测生产过程粉尘层着火危险性提供了新的方法.
[1] | Akash B A, Jaber J O. Characterization of shale oil as compared to crude oil and some refined petroleum products[J]. Energy Sources , 2003, 25 : 1171–1182. DOI:10.1080/00908310390233612 |
[2] | Burnham A K. Oil evolution from a self-purging reactor:kinetics and composition at 2 ℃/min and 2 ℃/h[J]. Energy Fuels , 1991, 5 : 205–214. DOI:10.1021/ef00025a034 |
[3] | Lin L X, Zhang C, Li H J, et al. Pyrolysis in indirectly heated fixed bed with internals:the first application to oil shale[J]. Fuel Processing Technology , 2015, 138 : 147–155. DOI:10.1016/j.fuproc.2015.05.023 |
[4] | Taciuk W. Does oil shale have a significant future?[J]. Oil Shale , 2013, 30 (1) : 1–5. DOI:10.3176/oil.2013.1.01 |
[5] |
王国金, 王剑秋, 李术元, 等.
油页岩颗粒燃烧释放二氧化硫的非催化气-固反应通用模型[J]. 燃料化学学报 , 1996, 24 (1) : 66–73.
( Wang Guo-jin, Wang Jian-qiu, Li Shu-yuan, et al. General noncatalyzed gas-soild reaction model for sulphur dioxide release during oil shale particle combustion[J]. Journal of Fuel Chemistry and Technology , 1996, 24 (1) : 66–73. ) |
[6] |
姜秀民, 韩向新, 刘德昌, 等.
油页岩着火机理的研究[J]. 发电设备 , 2002 (5) : 1–4.
( Jiang Xiu-min, Han Xiang-xin, Liu De-chang, et al. Research on the ignition mechanism of oil shale[J]. Power Equipment , 2002 (5) : 1–4. ) |
[7] |
韩向新, 姜秀民, 崔志刚.
油页岩颗粒燃烧的高温段多相燃烧模型[J]. 哈尔滨工业大学学报 , 2007, 39 (9) : 1470–1472.
( Han Xiang-xin, Jiang Xiu-min, Cui Zhi-gang. Heterogeneous combustion model for the high-temperature stage of oil shale particle combustion[J]. Journal of Harbin Institute of Technology , 2007, 39 (9) : 1470–1472. ) |
[8] | Ulrich K, Willi H.Hazards arising from electrical devices surrounded by deposits of flammable dusts[C]//Proceedings of 1994 International Colloquium on Safety Science and Technology.Shenyang:Northeastern University Press, 1994. |
[9] | Serth R W, Lestina T G. In process heat transfer[M]. Boston: Academic Press, 2014 . |
[10] |
王秋雯.吉林桦甸油页岩热物理和电物理性质实验研究[D].长春:吉林大学, 2011.
( Wang Qiu-wen.Experiment on thermal and electrical physical properties of oil shale in Jinlin Huadian area [D].Changchun:Jilin University, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10183-1011100930.htm ) |
[11] | Yuan C M, Huang D Z, Li C, et al. Ignition behavior of magnesium powder layers on a plate heated at constant temperature[J]. Journal of Hazardous Materials , 2013, 246 : 283–290. |
[12] | Bowes P C. Self-heating evaluation and controlling the hazards[M]. New York: Elsevier, 1984 . |
[13] | Zienkiewicz O C, Taylor R L, Fox D. Finite element method for solid & structural mechanics[M]. Oxford: Butterworth-Heinemann, 2014 . |