伴随经济发展,社会生产和生活过程中污水产量急剧增加,据统计,我国城市污水处理规模超过2亿m3/d,位居世界第一,由此产生的污泥量突破6 000万m3/a(以含水率80%计)[1]。市政污泥作为生活污水处理的主要产物,富含有机质,资源化利用是实现碳中和的重要选择。传统的污泥处置方式如发酵堆肥、直接填埋和掺烧发电等技术在大气污染物控制、重金属污染等方面相对落后,成为制约我国污泥处理行业的主要问题[2-3]。基于化学链技术固有的载氧体反应热效应和低污染排放的优势,污泥的化学链燃烧备受关注[4],与传统掺烧发电需污泥预干燥耗能巨大不同,化学链燃烧技术利用载氧体在2个串联反应器中(燃料反应器(FR)和空气反应器(AR))中不断还原、氧化再生的特性实现氧传递,并通过调节载氧体配比使得高含湿污泥可流化燃烧[5-6]。但污泥中灰分高,可能会造成灰沉积,载氧体的选取需保证活性和控制成本,因此更多学者选用活性较高、成本低廉的铁基载氧体作为污泥气化载氧体[7-8]。
目前污泥燃烧技术的研究主要是宏观层面的两步法(气化-燃烧)[9],郭晋荣等[10]研究了污泥热解气化阶段的特性,分析了污泥理化特性对燃烧性能的影响;聂鑫等[11]通过模拟研究了污泥与煤掺烧的燃烧排放;LYNGFELT等[12]模拟了大容量锅炉中污泥化学链燃烧的效率和排放等;李莹莹等[13]通过分子动力学模拟研究了煤热解阶段的微观机理,为理解煤中组分的迁移演变提供了思路。
上述工作仅考察了燃料的独立热解,对于载氧体的存在,如铁铝载氧体的晶格表面对污泥化学结构热解的影响鲜见报道,为了解污泥化学链燃烧过程中,载氧体促进热解的微观机理,笔者构建低阶污泥的化学结构,根据试验铁基载氧体的成分[14],选用AlFeO3为载氧体,构建污泥独立热解和在载氧体表面热解的模型,选取不同的升温速率和热解温度,利用LAMMPS软件进行分子动力学模拟,探究污泥分子团在化学链燃烧初始阶段的热解过程,获得典型产物生成机理及污染性元素N、S的迁移行为。
研究所用市政污泥取自青岛市团岛某市政污水处理厂,具有北方泥质属性。为准确构建污泥的化学结构,对污泥进行表征分析。将压滤后污泥放入鼓风干燥箱中105 ℃恒温干燥至恒重,将干燥后的污泥充分研磨,过筛至75 μm以下待用。采用元素分析仪(型号VARIO EL Ⅲ)测定污泥碳、氢、氮元素含量,采用全自动定硫仪(型号ZDL-9)测定污泥硫含量,氧元素由差减法计算所得。污泥的工业分析参考GB/T 212—2008《煤的工业分析方法》进行。市政污泥的工业分析和元素分析结果见表1。
表1 市政污泥的燃料分析
Table 1 Composition analysis of municipal sludge
工业分析/%MarFCarVarAarMdFCdVdAd78.951.1813.276.6005.5963.0331.37元素分析/%CdHdOdNdSdCdafHdafOdafNdafSdaf35.953.4223.315.350.6052.384.9833.967.790.87
为了更直观比较各元素含量,利用归一化原理计算得出H/C、O/C、N/C、S/C原子比分别为1.141 5、0.486 3、0.127 5、0.006 2。可知市政污泥的硫元素含量较低,为保证模型构建过程中至少有1个硫原子,碳原子数量应为160个,计算所构建的市政污泥干燥无灰基化学结构大致为C160H182O77N19S,分子量约3 633.03 g/mol。
为探究市政污泥原子间的具体成键方式与化学态,利用Thermo Scientific K-Alpha X射线电子光谱仪对干燥后的市政污泥进行XPS全谱与精细谱分析测试,分别如图1、2所示。
图1 市政污泥干燥基XPS全谱
Fig.1 XPS full spectra of municipal sludge drying base
图2 市政污泥干燥基XPS精细谱图
Fig.2 XPS fine spectra of municipal sludge sludge drying base
由图1可知,市政污泥中除C、H、O、N、S、P、Cl等元素外,还含有Ca、Na、Fe等金属元素,与Si共同构成污泥的灰分,符合污泥富含N、P等元素且灰分较高的特点。
通过Avantage软件对各光子能谱进行分峰拟合,总结出C、N、O、S等元素的化学态、成键方式,具体见表2,峰位信息参考Avantage自带数据库。
由表2可知,市政污泥中含C的化学基团主要有4种:分别为284.80 eV处的碳氢化合物中C—C键和C—H键;286.14 eV处包含醚、苯酚和羟基等中的C—O键;288.07 eV处羧基中的—COOH;289.03 eV处构成羰基中的CO。O的存在方式经检测认定为3种,主要为532.63 eV的碳氧双键(羰基氧、酯、醛、酮、羧酸),且531.50 eV处存在部分碳氧单键(醚键、羟基),530.33 eV处检测出少量金属化合物。N的存在形式主要是吡啶(399.37 eV)和吡咯(400.11 eV),此外还有季氮(401.03 eV)和少量氮氧化物(402.37 eV)。含S的化学基团主要为R—SH键(163.27 eV),即硫醇、硫酚,还有164.37 eV处的噻吩硫,168.60 eV处的砜型硫,以及169.70 eV处的少量无机硫。考虑到所构建结构中的S原子个数,建模时将巯基作为S的主要连接方式。
表2 污泥样品 C 1s、O 1s、N 1s、S 2p XPS图谱分析
Table 2 XPS C 1s、O 1s、N 1s、S 2p data of municipal sludge
元素峰官能团结合能/eV物质的量分数/% C 1sC—C、C—H284.8061.95C—O286.1427.98—COOH288.076.83CO289.033.27O 1s金属氧化物530.3317.11C—O531.5029.04CO532.6353.85N 1s吡啶氧399.3737.82吡咯氧400.1139.01季氮401.0315.29氮氧化物402.377.88S 2pR—SH163.2740.54噻吩硫164.3720.22砜型硫168.6023.65无机硫169.7015.59
为建立污泥样品的平均碳骨架结构,利用德国布鲁克公司Bruker 400M型核磁共振谱仪对市政污泥样品进行13C-NMR测试分析,结果如图3所示。乘用Origin软件对图谱进行分峰拟合,其中化学位移0~80×10-6为脂肪碳峰,100×10-6~165×10-6为芳香碳峰,175×10-6~200×10-6为羰基、羧基峰。峰值详细信息及物质的量分数见表3。
图3 市政污泥13C-NMR图谱
Fig.3 13C-NMR spectra of municipal sludge
表3 市政污泥结构参数
Table 3 Percentage of structural parameters of municipal sludge
参数fafCaf a′fNafHafPafSafBafalf∗alfHalfOal物质的量分数/%38.6727.8610.814.626.192.420.162.0460.0435.5717.656.82
注:fa为总sp2杂交碳;为羰基或羧基;f a′为芳香碳;
为非芳香碳;
为质子化芳香碳;
为键合羟基或醚基芳香碳;
为烷基化芳香碳;
为桥头芳香碳;fal为总sp3杂化碳;
为
为CH或
为键合氧的脂肪碳。
根据表3的12个峰值数据,计算桥碳与周碳比XBP,可以确定市政污泥结构中芳环的平均缩聚程度为
(1)
通过代入,可得计算结果为0.23,更接近萘的XBP(0.25),可知样品中芳香结构以萘环和苯环为主,其他芳香结构含量较少。
选取1个苯环、1个萘环、6个吡咯型氮、8个吡啶型氮、3个质子化氮、2个氮氧化物作为基础结构单元,在ChemDraw中建立整体市政污泥化学结构(图4)。其中各结构单元和官能团数量见表4,S原子以硫醇形式存在于碳链上,该结构首先于Chem3D中生成初始三维结构查看原子是否重叠,然后进行初步结构优化。
通常情况下,搭建的初始模型与平衡状态下的模型有较大差距。为减少后续模型搭建的计算过程,需对初始模型优化以使局部能量尽可能小。优化过程为在Material Studio软件中利用Forcite模块对化学结构模型进行循环10次的分子动力学退火模拟,退火温度设置为500~1 000 K,使势能面上分子结构优化且能量最小;最后进行几何结构优化,获得稳态几何构型与初始模型进行对比如图5所示。
利用Material Studio软件Amorphous Cell模块进行热解模型的结构优化。首先建立低密度(0.5 g/cm3)周期性边界条件的热解模型box1,每个盒子中包括2个优化后的化学结构,再利用VMD软件以charge原子类型输出LAMMPS可识别的data文件。在LAMMPS中调用ReaxFF反应力场于低温下进行100 ps的初始结构优化,采用NPT系综,温度设置为273 K,采用Berendsen控温器控温,温度阻尼系数500;压力为101 325 Pa,压力阻尼系数为500,时间步长为0.1 fs,计算步数则为1 000 000步,最终获得了密度约0.82 g/cm3的初始热解模型。污泥热解模型如图6所示,优化后box内部原子变得更紧密,污泥化学结构出现相互聚合缩聚,且无新物质产生,最终获得稳定的结构。
图4 市政污泥化学结构模型
Fig.4 Chemical structure model of municipal sludge
表4 市政污泥化学结构单元及官能团数量
Table 4 Number of chemical structural unit and functional groups of excess activated sludge
类型数量类型数量1—OH431—O—3吡咯型氮(N-5)6CO5吡啶型氮(N-6)8—COOH12质子化氮(N-Q)3—SH1氮氧化物(N-X)2
图5 市政污泥化学结构3D模型
Fig.5 Chemical structure 3D model diagram of excess activated sludge
图6 污泥热解模型
Fig.6 Sludge pyrolysis model
污泥在载氧体表面热解的模型选取了具有Pna21空间群的AlFeO3分子原胞,模型晶胞3组棱长为a=0.498 39 nm,b=0.85 544 nm,c=0.924 13 nm,晶胞3组棱之间夹角为α=β=γ=90°,包含8个Fe原子、8个Al原子和24个O原子,其晶胞的稳定结构如图7所示。
图7 AlFeO3晶胞的稳定结构
Fig.7 Stable structure of AlFeO3 unit cell
AlFeO3(001)为自然状态下最常见的暴露表面,建立8×5×1周期性平板模型,为减少原子间相互作用,在模型z方向上设置5 nm真空层,并将底部原子抬高0.2 nm避免计算时周期性影响,将污泥独立热解模型导入,得到污泥在载氧体表面热解模型如图8所示。
图8 污泥在载氧体表面热解模型
Fig.8 Pyrolysis model of sludge on the surface of oxygen carrier
因本文只关注污泥化学结构模型的热解,固定铁铝载氧体物质层,将污泥化学结构的分解和气体分子(团)向铁铝载氧体表面的吸附视为气固间相对运动的合成。因此,不再表述载氧体的晶格氧传递行为和载氧作用。
研究污泥独立热解时,针对不同升温速率对热解产物的影响,选择合适的热解温度模拟污泥在载氧体表面的恒温热解过程。具体模拟过程与模拟试验条件设置如下:① 在温升区间1 000~3 000 K内,分别采用升温速率8、16、32、64 K/ps,对污泥独立热解过程进行模拟,升温速率选取的主要依据是方便模拟步数的取整。控温方式为Berendsen温控器,温度阻尼系数为100;模拟过程产物输出截止值为0.3;时间步长为0.25 fs。② 选取2 500、3 000 K,对污泥独立热解与污泥在载氧体表面的热解模型分别进行300 ps恒温热解过程模拟。控温方式为Berendsen温控器,温度阻尼为100;模拟过程产物输出截止值为0.3;时间步长为0.25 fs,模拟步数为1 200 000。
截止值参考LAMMPS使用手册对ReaxFF反应力场分析的建议值,取值0.3时通常可获得良好且准确的结果。模拟过程发生在较短时间内(通常为皮秒级),远低于试验时间尺度(通常为秒级),因此采用较高的模拟温度和升温速率,增加原子间的碰撞使反应能在极短时间内发生,缩短化学反应的发生周期[15]。研究表明,温度评估策略能成功反映试验温度下的反应机理,参考煤热解ReaxFF MD在2 500和3 000 K的良好结果[16-17],选取该温度作为污泥恒温热解的模拟温度。模拟时间参考煤热解选取200~250 ps,为使污泥热解过程更彻底,模拟时间设置为300 ps。
通过统计LAMMPS输出的热解产物,总结了4种升温速率下最终输出产物见表5。
表5 1 000~3 000 K不同升温速率的热解产物
Table 5 Pyrolysis products with different heating rates in 1 000-3 000 K
8 K/ps 16 K/ps 32 K/ps 64 K/ps分子式数量分子式数量分子式数量分子式数量C12H9O3N21C11H10O3N21C137H110O35N181C29H22O10N21C114H94O33N131C53H45O19N81C11H12O6N21C77H57O22N121C17H18O7N21C5H2O4N1C16H16O21C6H6O4N1C25H18O8N41C5H7O1C2H21C8H12ON1C11H8O4N21C24H22O7N1CHN2C6H4ON2S1CH2O1C8H10O5N21CO9H2O31CO210C6H3O3N1C10H8O2N21CH2O1H2O29C8H9O5N1H2O26CO6C20H16O8N31CHO21C2H2O2C8H15O41CO4H3O1C2H61C2H2O1C9H11O3N1C18H12O9N31CH4O4C5H4O4N1C10H8O3N3S1H2O15H2N1C3H3O2N1C18H22O7N1CON1C2H1CO27N1C13H9O5N21C19H11O8N31C5H6O1C5H8O1C71H53O22N91CO29C2H4O22C3H6O1CO25C5HON1CH5O1C42H42O10N51C2H2O2C12H13O5N21C3H6O1CHO1C6H10O1C16H15O6N41C2H4O3
续表
8 K/ps 16 K/ps 32 K/ps 64 K/ps分子式数量分子式数量分子式数量分子式数量H4C2H4O4C13H12O4N1C10H9ON3S1H2S1CH2O2C3H6O1C3H61C3H4O21C4H8O21CH2O22C3H6O31CH2O21C2H31C6H8O21C36H25O6N41H3O21CH41C2H4O4C18H13O8N31C8H8O2N1C2H4O22C4H7O21C25H23O3N31CH4O1C10H11ON3S1C5H2O3N1C3H6O21CH31H7C8H7O4N1C20H14O10N31C9H11O5N1H2S2C12H9O3N21C12H12O3N21CH2O1C2H3O2CH2O21HO2C2HO1C8H16O31C2H41CHO1C8H7O4N1CHO1CH4O1H2S1C3H61CH31CO2O2HO2CHO1C2H21C3H61
热解产物分布大致可分为3类:活性自由基(·H、·OH、·O、·N等)、无机气体(CO、CO2等)和有机物分子。为描述有机物分子的类别,将其划分为有机气体(C1~C4)、轻质焦油产物(C5~C15)、重质焦油产物(C16~C40)及焦炭产物(C40+)[18],根据热解产物统计了不同升温速率下有机物的生成情况见表6。
表6 不同升温速率下各种有机物生成数量
Table 6 Quantities of various organic compounds generated at different heating rates
组分数量8 K/ps16 K/ps32 K/ps64 K/psC1~C47192319C5~C1561388C16~C404235C40+2211
由表6可知,不同升温速率下,污泥热解主要有机产物均为有机气体,且考虑热解过程生成的无机气体,随升温速率增大,污泥热解程度降低,重质焦油产物在16 K/ps后数量开始上升,说明过高的升温速率不利于污泥热解。为了粗略定量分析,将气体产物及轻质焦油看作污泥挥发分,计算挥发分各组分分子量总和并与模拟体系分子量相比,升温速率分别为8、16、32、64 K/ps时,挥发分占比分别为35.12%、54.31%、50.93%、48.26%,与污泥实际干燥基中挥发分63.03%相比,16 K/ps升温速率更接近实际值。
通过分析总产物发现热解过程中N元素主要迁移到焦油中,为了研究N、S元素在不同升温速率下的迁移情况同时避免对焦油产物做过多赘述,汇总不同升温速率下除焦油外的含N产物和含S产物,具体见表7。
表7 不同升温速率下除焦油外的含N产物和含S产物
Table 7 N-containing products except tar and all S-containing products at different heating rates
产物8 K/ps16 K/ps32 K/ps64 K/ps分子式数量分子式数量分子式数量分子式数量含N元素组分C114H94O33N131C53H45O19N81C137H110O35N181C77H57O22N121N1CON1CHN2C3H3O2N1C42H42O10N51C71H53O22N91H2N1
续表
产物8 K/ps16 K/ps32 K/ps64 K/ps分子式数量分子式数量分子式数量分子式数量含S元素组分 C10H8O3N3S1C10H11ON3S1H2S2C6H4ON2S1H2S1H2S1C10H9ON3S1
由表7可知,不同升温速率下N元素的最终产物除部分进入焦油组分外,均主要存在于焦炭和活性自由基中,这些可燃组分在后续燃烧阶段释放NOx气体。而不同升温速率对S元素的最终形态影响不大,S元素会转移到轻质焦油中或变成H2S气体,而这些组分在富氧气氛中易燃烧,H2S生成SO2。
在2 500、3 000 K,对污泥独立热解和载氧体表面的热解进行恒温热解模拟,通过比较产物分布探究载氧体在热解过程中的作用,最终产物的输出见表8。
表8 不同热解温度下污泥独立热解与载氧体表面热解的产物分布
Table 8 Product distribution of sludge pyrolysis alone and with oxygen carrier at different temperatures
污泥独立热解在载氧体表面热解2 500 K3 000 K2 500 K3 000 K分子式数量分子式数量分子式数量分子式数量C154H120O39N181C19H14O7N41C161H168O78N91C168H176O58N81C5H3O3N1C32H24O10N21H2O9CO41C9H13ON1C4H4O2N1CO30H2O20C2H3ON2C5H4O31C3H5O2CH2O3C7H5O31CH2N1H2ON1H2ON1C19H11O5N31C118H83O25N151C9H3O2N1HON3CO27C2HO21N22HO2H2O24CHON2HN21CO22N1C5H6O3N1C21H12O6N1C2H2O3C6H7O2N1CH3ON2HN4CH3O1C3H6O21C3H6O1H2N2CHO5C21H18O6N51C2H2N21C4H2N1C2H2O21CO10C2H4O21C4H1N3C2H4O22CH3O2N1CN2CN4C2H6O2C2H2O5HON2ON1CHO1CH42C2H4O1C3HN1CH4O21C4H4O21C2H42C3H42C2H7O31C4H7O1HO3C4H41C2H2O21C4H3O2N1CH2O2C2H2C3H6O31C7H9O6N1C2O1CH33C18H11O10N31CO14N2CON1C7H12O21C2H23CH3N1N24C3H6O1C4H4O1C2H21CHN1CH4O2CO22CO21H2N1C2H4O1H2O29C2H3O1C4H51HO2C5H6O1C2HO2C2HO2H2S2C2H4O1C2H6O1C4H4O1C13H10O2N21C3H4O1CH35C6H51CHO21C5H3O31C3H62C25
续表
污泥独立热解在载氧体表面热解2 500 K3 000 K2 500 K3 000 K分子式数量分子式数量分子式数量分子式数量CHN1C11H12ON1CON1C2H24C4H8O21CHO21C8H8O2N1HN3C2H41C8H6O3N1CHN2H2N21CH31C3H4O21H2S2C2H4O1C3H61H23CHNS1C2H41C2H1C2H51CH2O21C3H51C3H6O1C2H2O21CH41CH23CH2ON2C5H4O1C2O2C2H3O21C6H9N1C2H41C9H10O2N1C8H31H2O1ON1CH42CH2O1C2H31C2H32CH4O1C3H7O1C3H51H2S2C3H3O1H2S1O1C2H2O1H4
为直观对比不同温度下2种热解模型的产物生成情况,统计了热解温度2 500 K和3 000 K时产物分布见表9。
表9 不同热解温度下热解产物分布
Table 9 Distribution of pyrolysis products at different pyrolysis temperatures
组分数量(无载氧体)数量(有载氧体)2 500 K3 000 K2 500 K3 000 K无机气体和活性基团33477095C1~C421382944C5~C156751C16~C403210C40+1111
由表9可知,热解温度由2 500 K升高至3 000 K时,产物中无机气体、活性基团和轻质焦油数量明显上升,对应的重质焦油和焦炭数量减少,这是因为升高温度会促进如C—O键和桥键的断裂,释放出更多的有机分子和活性自由基,从而提高污泥热解程度。
综上所述,相同热解温度时,载氧体的存在有效促进了污泥的热解程度,其中无机气体与活性基团数量的提升程度超过了100%,主要原因为载氧体的存在催化了污泥裂解,生成了多种活性自由基,而生成的活性自由基与热解产物之间存在二次反应,促进了热解过程,加快污泥热解且热解程度变高;载氧体的存在能够有效减少焦油产生,这与实际污泥化学链燃烧过程相同,再次证明了热解模型的合理性,模拟结果具有参考价值。
由产物分布情况可知,2种结构体系中的S原子最终基本都转移到H2S气体中,仅在3 000 K有载氧体的情况下S原子还有部分存在于活性基团中,说明热解过程中以巯基(R—SH)存在的S原子最终以H2S方式从污泥中释放,然后参与后续燃烧反应。
无载氧体作用时,污泥热解过程中N元素主要存在于各种有机组分产物中,基本无含氮活性基团生成;而污泥在载氧体表面的热解过程出现了大量含氮活性基团及·H、·OH等自由基。污泥化学链转化过程中,污泥中燃料氮会转化为中间产物NH3和HCN,中间产物与自由基(·H、·O、·OH等)发生一系列基元反应[19-20]。NH3生成NO的反应机理为:NH3可与自由基·H、·OH作用生成NHi,生成的NHi继续与自由基生成NO,具体过程为
(2)
HCN生成NO的反应机理为:HCN则可被自由基·O、·OH等氧化成NCO,生成的NCO与自由基·O、·OH、·H反应进一步生成NO,具体为
(3)
HCN+·ONCO+·H,
(4)
NCO+NON2O+CO,
(5)
NCO+·H·NH+CO,
(6)
NCO+·ONO+CO。
(7)
载氧体的存在加剧了污泥结构中直链与侧链断裂,产生小分子自由基,研究表明[21],·H、·OH等自由基对NH3及HCN的形成有促进作用,载氧体会促进污泥热解过程自由基的释放,从而加速NH3与HCN生成。但产物中存在HCN与NO,但无NH3生成,可能是由于铁基载氧体促使NH3转化为HCN。
污泥在载氧体表面的热解过程中开始出现NO及部分N2。N2生成主要是NO与中间产物NH2的还原反应(式(8)),且铁基载氧体的存在促进了NO分解,有利于NO还原[22],可将NO生成看作NH3氧化生成NO与NO分解成N2的竞争反应。因此载氧体的存在能有效减少NOx排放。
NH2+NON2+H2O。
(8)
1)不同升温速率下,污泥独立热解的主要热解产物均为有机气体,过高的升温速率不利于污泥热解。N元素主要迁移至重质焦油中,此外还存在于焦炭和活性自由基团中,在后续燃烧中,N进入燃料型NOx。不同升温速率对S元素的最终形态影响不大,S元素主要转移到轻质焦油中或进入H2S气体,参与后续燃烧。
2)热解温度升高和载氧体存在会促进污泥热解,使产物中无机气体、活性基团和轻质焦油的生成数量增加,说明载氧体能有效催化重质焦油裂解。
3)载氧体的存在促使N元素迁移至活性基团中,活性基团与热解的中间产物反应生成NO,但铁基载氧体同时会促进NO分解,部分还原为N2,降低了后续燃烧中NOx的生成。
[1] 戴晓虎,张辰,章林伟,等. 碳中和背景下污泥处理处置与资源化发展方向思考[J]. 给水排水,2021,57(3):1-5.
DAI Xiaohu,ZHANG Chen,ZHANG Linwei,et al. Thoughts on the development direction of sludge treatment and resource recovery under the background of carbon neutrality[J]. Water & Wastewater Engineering,2021,57(3):1-5.
[2] FYTILI Despoina,ZABANIOTOU Anastasia. Utilization of sewage sludge in EU application of old and new methods:A review,renew[J]. Renewable and Sustainable Energy Reviews,2008,12:116-140.
[3] 杨裕起. 城市污泥处理处置技术研究进展[J]. 化工设计通讯,2020,46(2):223-231.
YANG Yuqi. Research progress of municipal sludge treatment and disposal technology[J]. Chemical Engineering Design Communications,2020,46(2):223-231.
[4] 吴健,沈来宏,蒋守席,等.基于双级燃料反应器的污泥化学链燃烧特性[J].化工进展,2018,37(7):2830-2836.
WU Jian,SHEN Laihong,JIANG Shouxi,et al. Chemical looping combustion of sewage sludge based on two-stage fuel reactor[J]. Chemical Industry and Engineering Progress,2018,37(7):2830-2836.
[5] 牛欣,肖军.污泥化学链燃烧过程中氮迁移转化特性研究[J].燃料化学学报,2017,45(4):505-512.
NIU Xin,XIAO Jun. Nitrogen transformation in chemical looping combustion of sewage sludge[J]. Journal of Fuel Chemistry and Technology,2017,45(4):505-512.
[6] 刘帅,黄振,何方,等. NiFe2O4为载氧体的生物质半焦化学链燃烧热力学模拟研究[J].新能源进展,2016,4(3):172-178.
LIU Shuai,HUANG Zhen,HE Fang,et al. Thermodynamic analysis of biomass char chemical looping combustion with NiFe2O4 as oxygen carrier[J]. Advances in New and Renewable Energy,2016,4(3):172-178.
[7] SOZEN Huseyin,GUO Qiangwei,LI Haibin,et al. Chemical-looping gasification of biomass in a 10 kWth interconnected fluidized bed reactor using Fe2O3/Al2O3 oxygen carrier[J]. Journal of Fuel Chemistry and Technology,2014,42(8):922-931.
[8] 魏泽华,刘道诚,荆洁颖,等.化学链燃烧中铁基载氧体研究进展[J].洁净煤技术,2019,25(3):19-27.
WEI Zehua,LIU Daocheng,JING Jieying,et al.Research progress on Fe-based oxygen carrier in chemical looping combustion[J].Clean Coal Technology,2019,25(3):19-27.
[9] 胡东海,黄戒介,李春玉,等. 固体化学链燃烧技术及污染物释放研究进展[J].洁净煤技术,2020,26(4) :1-10.
HU Donghai,HUANG Jiejie,LI Chunyu,et al.Research progress on chemical-looping combustion and pollutant release for solid fuels[J]. Clean Coal Technology,2020,26(4):1-10.
[10] 郭晋荣,贾里,王彦霖,等 . 城市污泥热解特性及理化性能研究[J]. 中南大学学报(自然科学版),2021,52(6):2023-2031.
GUO Jinrong,JIA Li,WANG Yanlin,et al. Study on municipal sludge physical-chemical properties and pyrolysis characteristics[J]. Journal of Central South University(Science and Technology),2021,52(6):2023-2031.
[11] 聂鑫,王一坤,谢海燕,等. 燃煤电站耦合污泥加压热水解系统发电的经济性分析[J]. 洁净煤技术,2021,27(S2):284-289.
NIE Xin,WANG Yikun,XIE Haiyan,et al. Economic analysis of the coal-fired power station coupling with sludge pressurized thermos-hydrolysis system[J]. Clean Coal Technology,2021,27(S2):284-289.
[12] ANDERS Lyngfelt,BO Leckner. A 1 000 MWth boiler for chemical-looping combustion of solid fuels-discussion of design and
costs[J]. Applied Energy,2015,157:475-487.
[13] 李莹莹,王杰平,李光跃,等. ReaxFF动力学分析煤热反应机理[J]. 洁净煤技术,2017,23(1):48-51,56.
LI Yingying,WANG Jieping,LI Guangyue,et al. Application of ReaxFF molecular dynamics on coal thermal reaction mechanism[J]. Clean Coal Technology,2017,23(1):48-51,56.
[14] WANG Kun,AN Zewen,WANG Fengyin,et al. Effect of ash on the performance of iron-based oxygen carrier in the chemical looping gasification of municipal sludge[J]. Energy,2021,231:120939.
[15] HONG Dikun,GUO Xin. Molecular dynamics simulations of Zhundong coal pyrolysis using reactive force field[J]. Fuel,2017,210:58-66.
[16] ZHENG Mo,LI Xiaoxia,NIE Fengguang,et al. Investigation of overall pyrolysis stages for liulin bituminous coal by large-scale reaxff molecular dynamics[J]. Energy & Fuels,2017,31(4):3675-3683.
[17] 冯炜,高红凤,王贵,等.枣泉煤分子模型构建及热解的分子模拟[J].化工学报,2019,70(4):1522-1531.
FENG Wei,GAO Hongfeng,WANG Gui,et al. Molecular model and pyrolysis simulation of zaoquan coal[J]. CIESC Journal,2019,70(4):1522-1532.
[18] FIDEL Castro-Marcano,MICHAEL F Russo Jr,ADRI C T van Duin,et al. Pyrolysis of a large-scale molecular model for illinois no. 6 coal using the reaxff reactive force field[J]. Journal of Analytical and Applied Pyrolysis,2014,109:79-89.
[19] ISABEL Vermeulen,CHANTAL Block,CARLO Vandecasteele.Estimation of fuel-nitrogen oxide emissions from the element composition of the solid or waste fuel[J]. Fuel,2012,94:75-80.
[20] CAO Jingpei,LI Liuyun,KAYOKO Morishit,et al. Nitrogen transformations during fast pyrolysis of sewage sludge[J]. Fuel,2013,104:1-6.
[21] 李彦坤,马晶晶,常国璋,等. 基于铁基载氧体的煤化学链气化还原过程中氮元素迁移行为[J]. 化工进展,2019,38(8):3631-3639.
LI Yankun,MA Jingjing,CHANG Guozhang,et al. Migration behavior of nitrogen in chemical-looping gasification reduction process of coal with an iron-based oxygen carrier[J]. Chemical Industry and Engineering Progress,2019,38(8):3631-3639.
[22] PELS J R,WOITOWICZ M A,KAPTEIJN F,et al. Trade-off between NOx and N2O in fluidized-bed combustion of coals[J]. Energy & Fuels,1995,9(5):743-752.
移动阅读
WANG Cuiping,LIANG Wenzheng,WANG Kun,et al.Molecular dynamics simulation of chemical structure pyrolysis of sludge based on iron-based oxygen carrier[J].Clean Coal Technology,2022,28(3):139-149.