生物燃料的开发和高效利用已成为人类社会发展进步的动力引擎，当前研究聚焦于其与不同煤质的混燃技术[1]、有机污泥燃料开发以及金属粉尘的燃烧特性[2-3]，在深化其燃烧动力学模型的基础上，探索开发提高其燃烧效率的设备及工艺。但是对于可燃有机粉尘的研究报道较少，GANI等[4]发现生物质燃料中纤维素含量越高，其热分解速率越高，而木质素含量越高，其热解速率越慢；HAYKIRI等[5]发现瓜子壳、菜籽、松果、棉花壳、橄榄油壳等农副产品的点燃温度均在 423～475K之间。MI等[6]对比研究了不同比例竹子和木材混合粉末的燃烧特性；LÓPEZ等[7]采用热重/差式扫描量热仪（TG/DSC）研究了不同木材的燃烧过程，并将其燃烧过程分为脱挥发分和炭化两个阶段。

目前，我国对跨境人民币结算业务与贸易关系的研究不多。研究“一带一路”倡议下两者是否存在相互促进关系的研究更是少之又少。

作为一种生物燃料，谷壳糠粉具备成本低廉、环保无腐蚀、易点燃、锅炉升温快，且燃烧后产物可作为钾肥改善土壤板结等优势。但在食品、农产品加工生产过程中常会产生大量谷壳糠粉等可燃轻质粉尘，如啤酒生产车间，在受限高温空间内易发生燃烧、爆炸等事故，存在严重的安全隐患。查阅文献发现，鲜有针对谷壳糠粉燃烧特性的研究报道，本文以西安某啤酒生产车间的谷壳糠粉为研究对象，通过粉尘层和粉尘云实验探索粒径对其最低着火温度（MIT）的影响，采用哈特曼管和锥形量热仪（CONE）测试不同粒径谷壳糠粉的爆炸下限、爆炸压力和热释放性能，基于TG/DSC研究其燃烧特性和燃烧动力学，为生产车间的安全运行提供一定的理论和实验基础。

1 实验材料和方法

1.1 实验原料

谷壳糠粉样品源于西安某啤酒企业平筛工艺过程除尘系统内的新制粉尘，主要由麦芽壳粉尘、大米粉尘及少量灰分组成，其工业和元素分析结果如表1所示，可见其具有高挥发分、低热值、低固定碳等特点。

利用振动筛将谷壳糠粉按照标准筛网粒度进行筛分，选取分别通过80目（粒径180～1250µm），100目（粒径 154～180µm），120目（粒径 120～154µm），140目（粒径 109～120µm），160目（粒径109～96µm）和180目（80～96µm）的筛余为研究对象，其粒度分布曲线如图1所示。将原料及筛分后不同粒径粉尘分别进行粉尘层、粉尘云、哈特曼管爆炸试验、热释放速率及TG/DSC分析。

1.2 样品表征

参照GB/T 16430—1996《粉尘层最低着火温度测定方法》，研究层厚分别为5mm和10mm时，不同粒径谷壳糠粉所形成粉尘层的 MIT。参照GB/T16429—1996《粉尘云最低着火温度测定方法》测定不同样品所形成粉尘云的 MIT。参照 GB/T 16425—1996《粉尘云爆炸下限浓度测定方法》，以烘干、筛分后得到的不同粒径谷壳糠粉为研究对象，采用1.2L哈特曼管测试其爆炸下限浓度、爆炸压力上升速率、最大爆炸压力。实验条件：环境温度为29℃、点火延迟时间为 8s、喷粉压力为 0.5～0.6MPa。

 

表1 谷壳糠粉的元素分析及工业分析

  

元素分析/%工业分析（干燥条件）/% 低位发热量/MJ·kg-1 C H O N 固定碳 水分 灰分 挥发分44.825.84 44.85 2.48 13.86 7.86 5.50 73.09 16.19

  

图1 谷壳糠粉样品的粒度分布

参照GB/T 16172—2007《建筑材料热释放速率试验方法》标准，称取18.0g谷壳糠粉采用锥形量热仪（ZY6243-I东莞中诺）测试其燃烧时的释热性能，辐照强度为30kW/m2。

汉代服饰是中国传统服饰长河中一颗璀璨的明珠，汉代服饰种类的多样及其发展变化对中国传统服饰的发展有着重要影响，汉代的开放和经济的发达以及汉代博大精深的文化使得汉代服饰也具有丰富的文化内涵，不仅体现了汉代的时代特色，等级制度的森严而且反映出人们审美意识。它丰富的文化内涵使得它在中国传统服饰历史中具有承上启下的作用，我们应该通过研究汉代服饰了解汉代的文化，从小的缩影中窥探整个中国的传统服饰文化，同时应该保护汉代服饰等一系列的传统文化代表。

采用 Mettler(德国)热重-差式扫描量热分析仪分析谷壳糠粉样品的燃烧动力学，空气氛围，温度范围为50～600℃，升温速率为20℃/min。

2 结果与分析

2.1 最低着火温度

另外，对比同一种样品所形成粉尘层和粉尘云的MIT可见：粉尘云的着火温度远高于粉尘层，可能由于在粉尘云状态时，样品均匀地分散在空气中，空气的热导率较低，导致样品颗粒表面的温度上升较慢。

 

表2 谷壳糠粉尘层的最低着火温度

  

样品 MIT /℃5mm 10mm原料 135 130 80目（180～1250µm） 140 135 100目（154～180µm） 135 130 120目（120～154µm） 130 125 140目（109～120µm） 130 130 160目（109～96µm） 130 130 180目（80～96µm） 130 130

 

表3 谷壳糠粉尘云的最低着火温度

  

样品 MIT/℃ 样品 MIT/℃原料 480 140目(109～120µm) 440 80目（180～1250µm） 490 160目(109～96µm) 430 100目(154～180µm) 470 180目(80～96µm) 430 120目(120～154µm) 450

谷壳糠粉尘层及粉尘云的最低着火温度分别见表2、表3。由表2可见：谷壳糠粉尘层厚度相同时，粉尘层MIT随粒度的减小呈略降趋势，粒径减小，比表面积增大，与空气的接触面积增大，导致粒径较小的粉尘更易燃烧；但当样品粒径减小至140目时，其MIT维持在130℃，表明随着粒径的减小，氧气浓度并非影响其燃烧的唯一因素；对于粒径相同的谷壳糠粉，其MIT随粉尘厚度增加而减小，由于粉尘层厚度增加使其内部蓄热性能增强，易使粉饼中样品较充分活化；同时由于粉尘层从底部加热燃烧，开始阶段产生的CO2等非助燃性气体在上升过程中可阻滞顶层表面O2的扩散[8]，饼层越薄，阻滞效果越强，导致MIT随饼层厚度增加而减小。POLKA等[9]对比研究了层厚分别为12.5mm和5mm可燃粉尘的MIT，发现厚度为12.5mm粉尘层的MIT较5mm厚的降低20～30℃。从表3谷壳糠粉尘云的MIT可见：随着样品粒径的减小，粉尘云MIT也随之减小，但当样品粒径减小至160目时，MIT维持在430℃。

2.2 粉尘爆炸下限

从表4可见：随着粉尘粒径增大，其爆炸下限（LEL）浓度随之增大，但最大爆炸压力P和爆炸压力上升速率均随之降低。由于粉尘粒径减小，比表面积增大，与空气接触得更充分，被点火源引爆的可能性增加，同时由于小粒径粉尘的质量较小，导致在爆炸时压力上升速率和最大爆炸压力P随之增大。其结果与文献报道相似，BIDABADI等[10-11]发现铁粉在爆炸时随着粒径减小，其火焰燃烧增速变大。根据有机物在燃烧时的“热解-挥发”理论[12]，当谷壳糠粉的粒径减小时，加速了半纤维素等热解形成可燃性气体的过程，即在火焰前锋扩散到颗粒表面之前已发生热解，形成可燃性气体，其与火焰前锋更易发生燃烧，且在瞬间可释放大量热量，促成爆炸现象的发生。

 

表4 不同粒径粉尘爆炸特征参数

  

粉尘粒径/µm LEL浓度/g·m-3 P/MPa 爆炸压力上升速率/ MPa·s-1 80～96 50～58.33 0.90 9.46 96～109 58.33～66.67 0.89 9.35 109～120 64～71.48 0.88 8.73 154～120 75～83.33 0.88 7.33 180～154 116.67～125 0.81 7.28 1250～180 141.67～150 0.72 6.05

2.3 热释放性能

从图2可见，谷壳糠粉在燃烧过程中其热释放速率（HRR）曲线上基本呈现两个峰，第一个由样品中半纤维素等易燃热解物快速燃烧而形成，同时样品的炭化会导致释热速率逐步减小，由于样品的厚度较小（约 5mm），纤维素、木质素及剩余未燃样品等在持续受热作用下燃烧并产生释热速率的第二个峰值。随着谷粉粒径减小，其燃烧时所需的点火温度降低，但当粒径减小至180目时，由于其蓄热能力增强和燃烧时炭化灰壳的阻滞作用，导致其点火时间TTI延长，最大释热峰的出现时间滞后，由 169s延滞为 307s，两个峰合并为一个峰。结合不同粒径谷壳糠粉燃烧时的烟气温度（如图3所示），180目样品燃烧时在321s时产生的最高烟气温度为59℃，而80目样品在226s产生的最高烟气温度为57℃，即随着谷粉粒径的减小，其烟气温度峰值也发生右移；同时粒径为180目样品的烟气温度曲线也呈现一个较大的峰，而80目和120目样品的烟气温度均呈现明显的“双峰”特征。

式(2)中(dw/dτ)mean为平均燃烧速度，%/min，由式(3)计算，τ表示时间，min；Ti为着火温度，K；Th为燃尽温度，K。燃烧特性指数SN综合考虑了燃料的着火和燃烧稳定性等因素对燃尽的影响，是反应物质着火和燃尽的综合特性指标，其值越大，燃烧特性越大。

  

图2 不同粒径样品的热释放速率

  

图3 不同粒径样品的烟气温度

从表5可见，随着谷壳糠粉粒径的减小，其总放热量、有焰燃烧时间和失重率逐渐减小，其中180目样品的总放热量为307.89kJ，而80目样品的总放热量为376.62kJ，结合它们的失重率结果，可推断180目样品燃烧不够充分导致其放热量减小；另外，随样品粒径减小其点燃时间延长，180目样品的点火时间为93s，而80目样品仅为69s，但单位面积热释放速率峰值和平均热释放速率基本不变。表明增大样品燃烧过程中的氧气浓度（颗粒粒径减小），其燃烧速率和总释热量并未显著增大，一定程度上揭示其燃烧过程并非由扩散过程控制。

 

表5 不同粒径样品的热释放性能

  

样品 总放热量/kJ点火时间/s有焰燃烧时间/s平均热释放速率/ kJ/m2单位面积热释放速率峰值/kJ/m2失重率80目 376.62 69 384 20.01 60.67 90.40%120目 344.47 83 338 22.94 60.21 89.30%180目 307.89 93 300 21.84 61.76 82.90%

2.4 燃烧特性及动力学

2.4.1 燃烧特性

样品的燃烧特性由以下参数确定[13-14]，即可燃特性指数Cr、燃烧特性指数SN，其计算式见式(1)、式(2)。

 

没心情去参加诗歌朗诵会的皇甫一兰下班后，被余科长硬拉进他的私家车绑架到电视台。两人到了三楼刚出电梯口，余科长的手机“哟咿哟咿”地叫起来了，极难听。皇甫一兰早就抗议他换了怪叫声，余科长说我喜欢这女人的喊叫，怎么着？

 

技术推广力量不足。乡镇基层既懂农机又懂农艺技术的人员少，难以满足大面积推广水稻育插秧技术的需求，限制了机插秧技术的大面积推广应用。

 

根据 Coats-Redfern积分法研究其燃烧动力学[19-21]，依据方程式(4)～(6)，采用反应级数n为1，不同粒径谷壳糠粉在 20℃/min升温速率下的热解动力学参数如表7所示，由于3种样品的着火温度均在240～270℃，其燃烧机理属于静态渗透式扩散燃烧，燃烧过程从着火后开始。在282.4～360.1℃之间，主要发生谷壳糠粉的燃烧反应，且对应拟合曲线的相关系数R2均接近于1，表明采用一级反应拟合较为合理；在该阶段样品的活化能E和指前因子A随着粒径的减小而增大，其中粒径为80～96µm样品的活化能达到 51.15kJ/mol，而 180～1250µm样品的活化能为35.35kJ/mol，表明微细谷壳糠粉发生燃烧反应的能量壁垒更高，反应过程由燃烧反应动力学过程控制。

式中，，对ln[G(α)/T2]-1/T作图为一直线，直线的斜率和截距可求得动力学参数 E和A，分别选取反应级数为1、1.5、2、2.5、3试算，根据拟合后的相关系数确定最佳反应级数n。

式(1)中(dw/dτ)max为燃烧反应的最大速率，常用DTGmax表示，%/min；Ti为着火温度，K，常采用TG-DTG方法确定[15]。

2.4.2 燃烧动力学

根据参考文献[16-17]，结合图4(d)中的DTG曲线，可将谷壳糠粉的燃烧过程分为以下4个阶段：水分的蒸发（50～151℃），半纤维素的热解（151～239℃），纤维素的热解（239～335℃），木质素的热解（335～525℃）。由于谷壳样品中的木质素含量较低，约为15%[18]，从图4(d)可见随着谷壳糠粉粒径的减小，其木质素热解对应的失重峰逐渐消失。纤维素是由D-葡糖基组成的直链高分子，较半纤维素更难于热解，随着样品粒径的减小，其蓄热能力增强，导致半纤维素和木质素的热解峰变小，而对应纤维素的热解峰增大。

式(3)中wl为到达燃尽点时的样品失重率，%，VT为升温速率，/min℃。

门只好敞着。呼伦继续冥思苦想，可是一个小时过去，纸上再也没有落下一个字。云梦和丈母娘把一部弱智电视剧看得兴趣盎然，每逢精彩或者逗乐处，两个人都要交换一番意见或者交换一下笑声。当然那笑经过了压抑，可是经过压抑的笑声更加可怕更加刺耳。就像两只同时下蛋的母鸡被人同时捏住嘴巴，然那两只母鸡还在憋足劲儿叫，于是声音变着调子往外挤，拐着弯儿往外钻，拖着尾音往外蹦，分着叉儿往外蹿，咯咯咯咯答，答答答答咯，咯咯咯咯咯咯咯……可怜的呼伦如坐针毡，苦不堪言。

 

表6 不同粒径谷壳糠粉的燃烧特性参数

  

样品/目 TiphmaxlNr/K T/K T/K DTG/%·min-1 w/% S/10-7 C/10-5 平均燃烧速率/%·min-1 80 515 574 789 9.47 73.46 2.43 3.57 5.36 120 533 581 808 11.22 70.11 2.49 3.95 5.11 180 550 583 811 16.63 73.46 3.82 5.49 5.63

 

  

图4 不同粒径样品的热重及其微分曲线TG/DTG

 

表7 不同粒径谷壳糠粉的热解动力学参数

  

颗粒 温度区间/℃ E/kJ·mol-1 A/min-1 R2 80目 约188.7 1.47 0.003 0.99 188.7～226.9 14.80 0.897 0.95 226.9～282.4 20.04 4.269 0.95 282.4～360.1 35.35 242.845 0.98 360～537.3 3.81 0.064 0.84 120目 约188.7 2.12 0.006 0.99 188.7～226.9 12.12 0.469 0.92 226.9～282.4 18.55 3.219 0.99 282.4～360.1 39.86 803.421 0.99 360～537.3 11.64 0.963 0.99 180目 约188.7 0.95 0.002 0.99 188.7～226.9 6.31 0.049 0.87 226.9～282.4 15.28 1.015 0.95 282.4～360.1 51.15 9554.203 0.99 360.1～537.4 3.41 0.061 0.89

式中，α为样品燃烧时的质量转化率，n为反应级数，A为指前因子，即频率因子，min-1；E为活化能，kJ/mol；R为通用气体常数，kJ/(mol·K)；β为升温速率，℃/min；T为热力学温度，K。

 

结合图4(a)～(c)不同样品的热重及其微分曲线并计算得到不同粒径谷壳糠粉的燃烧特性参数如表6所示。从表6中可见：随着谷壳糠粉粒径的减小，着火温度 Ti、最大燃烧温度 Tp和燃尽温度 Th均略微增大，最大燃烧速率 DTGmax显著增大，其中当粒径为180目样品的最大燃烧速率达16.63%/min，较80目样品提高了76%；另一方面燃烧特性指数SN也随样品粒径减小而增大，其中180目样品的SN达到 3.82×10-7，较80目样品提高了 57.2%，表明随着粒径减小，燃烧特性显著增加，样品在有焰燃烧阶段变得更加剧烈。

结合粉尘层实验时的观察结果：厚度较薄粉饼燃烧时易开裂、翘起，厚度较大粉饼燃烧时易产生烟尘，但当除去表层燃尽灰壳时，产生瞬间复燃现象。由于谷壳糠灰中含有一定的碱金属元素[22]，可降低灰渣的熔点[23-25]；还有一定量的SiO2[26]，其在燃烧时易于结渣，促使形成一定厚度的灰壳层，可阻碍火焰的扩散，导致点燃时间延长，MIT升高。

有限元模拟分析得到的应力云图和试验时试件破坏的图片如图5所示，在与试验加载工况相同的情况下，由应力云图和等效塑性应变(PEEQ)对比试验结束后试件的破坏情况可知，在靠近翼缘处应力集中现象较明显，底部翼缘部位预埋件处的应力较集中，有限元模拟与试验结果基本相吻合。综上所述，所建立的新型消能减震复合墙板有限元模型具有较高的准确性，能够较好地反映构件的力学性能，进一步论证了有限元模型的正确性。

根据以上的实验数据分析可知：在持续受热及空气环境下，谷壳糠粉颗粒表面可燃挥发物首先吸热并分解，形成的可燃性气体与氧气发生放热反应产生火焰；随后表层部分碳处于过渡燃烧态，可形成部分薄灰壳；谷壳糠粉颗粒内部热解的可燃气体向表面扩散并发生燃烧，颗粒被火焰包围；但随着颗粒粒径的不断减小，增大颗粒与O2接触面积的同时，其蓄热性能也增强，导致样品的点燃时间延长，当积蓄的能量达到一定程度后易发生“轰燃”，此时由于颗粒表层易燃的半纤维素等与O2燃烧，温度迅速升高，导致颗粒内层有机物迅速炭化形成灰壳，致使反应的活化能增大，谷壳粉尘不能发生完全氧化，热重损失减小。LI等[27]发现随着金属铝粉颗粒粒径的减小，其燃烧过程由扩散控制逐渐转变为反应控制模型。本研究发现随着谷壳糠粉粒径的减小，其燃烧过程先随氧气浓度的增大而加速，随着粒径进一步减小，其燃烧因炭化灰壳的阻滞作用而减缓，表明其燃烧过程由扩散控制转变为反应动力学控制。

3 结论

（1）粉尘层和粉尘云实验结果表明：谷壳糠粉的MIT先随粒径的减小而降低，随后保持不变，同种样品随堆积厚度的增大而降低。

（2）1.2L哈特曼管实验结果表明：随着粉尘粒径增大，其LEL浓度随之增大，但最大爆炸压力P和爆炸压力上升速率均随之降低。

4.2 水分：播前底水必须用开水，而且要足，播种后覆土，马上盖地膜，以保持较高土温。苗床用水应事先予热或对热水，使水温达到25℃左右后使用，而且浇水也应选在晴天的早晨，苗床浇水切忌大水漫灌，浇水后更要加强通风，使苗床空气湿度控制在70%以下，保证幼苗的叶片早晨不结露。在灌水后遇阴雨天，可在露水下去后向苗床撒些草木灰或细干土吸湿，同时也防止土壤排湿。

（3）CONE结果表明：随着颗粒粒径的减小，谷粉的蓄热能力增强，导致最大释热峰的出现时间由169s（80目）延滞为307s（180目），热释放速率和烟气温度变化曲线均呈现一个峰。

（4）TG/DSC结果表明：随着谷壳糠粉粒径的减小，其燃烧过程先随氧气浓度的增大而加速，随着粒径进一步减小，其燃烧因炭化灰壳的阻滞作用而减慢，热解过程的反应活化能由 35.35kJ/mol（180～1250µm样品）增大至51.15kJ/mol，其燃烧过程由扩散控制转变为反应动力学控制。

参考文献

[1]WANG M，ZHAO R，SHAN Q，et al.Study on combustion characteristics of young lignite in mixed O2/CO2，atmosphere[J].Applied Thermal Engineering，2017，110：1240-1246.

[2]ZHANG Q，XING.Pyrolysis characteristic and kinetic analysis of different deep dewatered sludges[J].Environmental Chemistry，2014，170：325-330.

[3]SUNDARAM D S，PURI P，YANG V.A general theory of ignition and combustion of nano- and micron-sized aluminum particles[J].Combustion & Flame，2016，169：94-109.

[4]GANI A，NARUSE I.Effect of cellulose and lignin content on pyrolysis and combustion characteristics for several types of biomass[J].Renewable Energy，2007，32（4）：649-661.

[5]HAYKIRI H.Combustion characteristics of different biomass materials[J].Energy Conversion & Management，2003，44（1）：155-162.

[6]MI B B，LIU Z J，HU W H，et al.Investigating pyrolysis and combustion characteristics of torrefied bamboo，torrefied wood and their blends[J].Bioresource Technology，2016，209：50-55.

[7]LÓPEZ D，AVALOS A，GIROIR A，et al.Combustion kinetic study of woody and herbaceous crops by thermal analysis coupled to mass spectrometry[J].Energy，2015，90：1626-1635.

[8]WU D，VERPLAETSEN F，BULCK E V D.Experimental analysis of minimum ignition temperature of coal dust clouds in oxy-fuel atmospheres[J].Clean Coal Technologies，2014，84：330-339.

[9]POLKA M，SALAMONOWICZ Z，WOLINSKI M，et al.Experimental analysis of minimal ignition temperatures of a dust layer and clouds on a heated surface of selected flammable dusts[J].Procedia Engineering，2012，45（3）：414-423.

[10]BIDABADI M，HAJILOU M，POORFAR A K，et al.Modeling flame propagation of micron-sized iron dust particles in media with spatially discrete sources[J].Fire Safety Journal，2014，69：111-116.

[11]BIDABADI M，MOHEBBI M，POORFAR A K，et al.Modeling quenching distance and flame propagation speed through an iron dust cloud with spatially random distribution of particles[J].Journal of Loss Prevention in the Process Industries，2016，43：138-146。

[12]GAO W，MAGI T，SUN J，et al.Effects of particle size distributions on flame propagation mechanism during octadecanol dust explosions[J].Powder Technology，2013，249：168-174.

[13]GEGENFURTNER A.Investigation of combustion characteristics and kinetics of coal gangue with different feedstock properties by thermogravimetric analysis[J].Thermochimica Acta，2015，614（8）：137-148.

[14]ZHAO D，ZHANG J，WANG G，et al.Structure characteristics and combustibility of carbonaceous materials from blast furnace flue dust[J].Applied Thermal Engineering，2016，108：1168-1177.

[15]CHEN Z H，HU M，ZHU X L，et al.Characteristics and kinetic study on pyrolysis of five lignoc ellulosic biomass via thermogravimetric analysis[J].Bioresource Technology，2015，192：441-450.

[16]MA Z Q，CHEN D Y，GU J，et al.Determination of pyrolysis characteristics and kinetics of palm kernel shell using TGA-FTIR and model-free integral methods[J].Energy Conversion & Management，2015，89：251-259.

[17]ELSAYED S A，Mostafa M E.Pyrolysis characteristics and kinetic parameters determination of biomass fuel powders by differential thermal gravimetric analysis（TGA/DTG）[J].Energy Conversion &Management，2014，85：165-172.

[18]SUN Y，CHENG J.Hydrolysis of lignocellulosic materials for ethanol production：a review[J].Bioresource Technology，2002，83（1）：1-11.

[19]BRANCA C，BLASI C D.A unified mechanism of the combustion reactions of lignocellulosic fuels[J].Thermochimica Acta，2013，565：58-64.

[20]ZHANG Y，NIU Y，ZOU H，et al.Characteristics of biomass fast pyrolysis in a wire-mesh reactor[J].Fuel，2017，200：225-235.

[21]CAI J，WANG S，KUANG C，et al.Insight into the kinetic analysis of catalytic combustion for biomass after alkaline metals loaded pretreatment[J].Fuel，2017，203：501-513.

[22]BAZARGAN A，GEBREEGZIABHER T，HUI C W，et al.The effect of alkali treatment on rice husk moisture content and drying kinetics[J].Biomass & Bioenergy，2014，70：468-475.

[23]SAHU S G，CHAKRABORTY N，SARKAR P.Coal-biomass co-combustion：an overview [J].Renewable & Sustainable Energy Reviews，2014，39（6）：575-586.

[24]DU S，YANG H，QIAN K，et al.Fusion and transformation properties of the inorganic components in biomass ash[J].Fuel，2014，117（1）：1281-1287.

[25]SUN Z，SHEN J Z，JIN B S，et al.Combustion characteristics of cotton stalk in FBC [J].Biomass & Bioenergy，2010，34（5）：761-770.

[26]BAZARGAN A，BAZARGAN M，MCKAY G.Optimization of rice husk pretreatment for energy production [J].Renewable Energy，2015，77：512-520.

[27]LI Q，WANG K，ZHENG Y，et al.Explosion severity of micro-sized aluminum dust and its flame propagation properties in 20L spherical vessel [J].Powder Technology，2016，301：1299-1308.