浙江大学学报(工学版), 2026, 60(3): 670-678 doi: 10.3785/j.issn.1008-973X.2026.03.023

环境工程

阴燃治理含油污泥反应特性的实验及模拟分析

张烨,, 籍龙杰, 李鸿炫, 刘鹏, 李书鹏, 冯时, 王进卿, 徐旭,, 詹明秀

1. 中国计量大学 能源环境与安全工程学院,浙江 杭州 310018

2. 北京建工环境修复股份有限公司,北京 100015

3. 污染场地安全修复技术国家工程实验室,北京 100015

Experimental and simulation analysis on reaction characteristic of smoldering treatment of oily sludge

ZHANG Ye,, JI Longjie, LI Hongxuan, LIU Peng, LI Shupeng, FENG Shi, WANG Jinqing, XU Xu,, ZHAN Mingxiu

1. College of Energy Environment and Safety Engineering, China Jiliang University, Hangzhou 310018, China

2. Beijing Construction Engineering Group Environmental Remediation Limited Company, Beijing 100015, China

3. National Engineering Laboratory for Site Remediation Technologies, Beijing 100015, China

通讯作者: 徐旭,女,教授. orcid.org/0000-0003-3697-2787. E-mail:xuxu@cjlu.edu.cn

收稿日期: 2025-04-4  

基金资助: 浙江省自然科学基金资助项目(LY23E060002).

Received: 2025-04-4  

Fund supported: 浙江省自然科学基金资助项目(LY23E060002).

作者简介 About authors

张烨(1998—),男,硕士生,从事固体废弃物处置的研究.orcid.org/0009-0002-7499-4944.E-mail:1287509205@qq.com , E-mail:1287509205@qq.com

摘要

针对现有含油污泥处置技术存在的成本高、难以兼顾无害化和资源化的问题,采用阴燃法处置含油污泥,探究处置效果以及在不同达西空气流速、热值条件下的阴燃特性差异,建立一维的含油污泥阴燃数值模型. 结果表明,阴燃处理后污泥中的污染物去除率大于99.7%,油回收率可达57.8%. 明确了污泥阴燃的平均峰值温度和传播速度与水的质量分数、达西空气流速、孔隙率的关系. 对比阴燃前、后的物质种类和质量分数变化并结合分子动力学模拟,发现阴燃过程的反应不仅涉及长链石油烃的裂化和短链石油烃的氧化,还包含石油烃的挥发、环化和芳构化、异构化等. 数值模拟结果表明,利用该模型,能够较好地模拟含油污泥阴燃过程的温度和速度,平均误差约为11%. 当达西空气流速和热值改变时,模拟结果与实验结果的阴燃速率变化趋势一致,模拟结果与实验结果在不同条件下符合较好.

关键词: 含油污泥 ; 阴燃特性 ; 降解路径 ; 分子动力学模拟

Abstract

Smoldering treatment was applied to oily sludge to evaluate the disposal effectiveness and analyze the difference in smoldering characteristic under different Darcy air flow rate and heating value aiming at the problem of the high cost and the difficulty of simultaneously achieving harmlessness and resource recovery of existing oily sludge treatment technology. A one-dimensional numerical model of oily sludge smoldering was established. Results showed that the pollutant removal rate in the sludge exceeded 99.7% and the oil recovery rate reached 57.8% after smoldering treatment. The relationship of the average peak temperature and propagation velocity with the mass fraction of water, Darcy air flow rate, and porosity was analyzed. The types and mass fractions of substances before and after smoldering were compared by combining with molecular dynamics simulations. The reactions involved in the smoldering process include not only the cracking of long-chain petroleum hydrocarbons and the oxidation of short-chain hydrocarbons, but also volatilization, cyclization and aromatization, and isomerization of petroleum hydrocarbons. The numerical simulation results showed that the model could reasonably simulate the temperature and propagation velocity during oily sludge smoldering, with an average error of about 11%. The simulated trends in smoldering rate accorded with the experimental results when the Darcy air flow rate and heating value varied, and good agreement was obtained under different conditions.

Keywords: oily sludge ; smoldering characteristic ; degradation pathway ; molecular dynamic simulation

PDF (1331KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张烨, 籍龙杰, 李鸿炫, 刘鹏, 李书鹏, 冯时, 王进卿, 徐旭, 詹明秀. 阴燃治理含油污泥反应特性的实验及模拟分析. 浙江大学学报(工学版)[J], 2026, 60(3): 670-678 doi:10.3785/j.issn.1008-973X.2026.03.023

ZHANG Ye, JI Longjie, LI Hongxuan, LIU Peng, LI Shupeng, FENG Shi, WANG Jinqing, XU Xu, ZHAN Mingxiu. Experimental and simulation analysis on reaction characteristic of smoldering treatment of oily sludge. Journal of Zhejiang University(Engineering Science)[J], 2026, 60(3): 670-678 doi:10.3785/j.issn.1008-973X.2026.03.023

含油污泥是原油开采、炼油厂污水处理过程中产生的固体废弃物,产量高,全球范围内因缺乏适当处理而积累的含油污泥超过10亿吨[1]. 含油污泥的成分复杂,具有含水率高、流动性差、体积大等特点,处理难度较高[2-3]. 目前,含油污泥的处理方法有热解、焚烧和溶剂萃取等方法[4-5],热解处理的设备制造、运行和维护成本高,焚烧法无法回收含油污泥中的油分,亟需处理效果好、成本低的新型含油污泥治理技术[6-7].

阴燃是缓慢、低温、无焰的燃烧形式,常发生于比表面积较大的多孔有机固体[8]. 在适当条件下,阴燃燃烧可以自我维持[9]. Zhao等[10]验证了阴燃处置含油污泥的可行性,明确回收油可资源化利用,探讨水的质量分数和物料配比对阴燃性能的影响. 杨高玄等[11-12]研究含油率、水的质量分数及介质粒径等参数对含油污泥阴燃的影响和石油烃降解效率. Gan等[13-14]探究水的质量分数、介质粒径、空气流速等对含油污泥阴燃的影响,分析回收油的成分. 雷大鹏等[15]的中试研究表明,阴燃处置后含油污泥中总石油烃的去除率可以大于98%. 目前,含油污泥阴燃的研究主要侧重污染物去除效果、影响因素以及热解过程元素迁移转化特性等[16],较少涉及含油污泥阴燃过程的反应机制和污染物降解路径. 此外,含油污泥阴燃的数值模拟工作开展较少,阴燃速率的预测结果存在较大的误差[14]. 开展含油污泥阴燃过程的反应机制、污染物迁移路径和含油污泥阴燃预测模型研究,可为其无害化、资源化处置提供理论支撑.

本研究旨在探究含油污泥处置的效果,探讨达西空气流速、热值、物料配比等关键参数对阴燃特性的影响,分析阴燃前、后物料的无害化和资源化程度,阐明含油污泥的阴燃反应机制. 构建优化的含油污泥阴燃数值模型,为阴燃温度和速率的精准预测提供理论依据.

1. 材料与方法

1.1. 实验原料

以某落地含油污泥为研究对象,利用元素分析仪和工业分析仪分析样品的成分,如表1所示. 含油污泥中固定碳的质量分数较小,占比为0.59%;灰分的质量分数较高,占比为55.48%;水分和挥发分的质量分数相近,质量占比分别为20.43%和23.50%;含油污泥中碳元素的质量分数为12.74%,与石油烃的主要成分特征相符[17-18]. 利用氧弹热量仪测定含油污泥的高位热值为11 178 kJ/kg. 利用气相色谱仪测定原始含油污泥中的石油烃组分,其烃类组成以正构烷烃为主,碳数主要分布在C17~C25. 其中二十烷(C20H42)的质量分数最高,达到15.73%;其次为十七烷(C17H36,8.12%)和二十二烷(C22H46,6.32%),同时检测到十八烷(5.60%)、二十四烷(5.65%)组分,如表2所示.

表 1   落地含油污泥的工业分析和元素分析

Tab.1  Proximate analysis and ultimate analysis of landed oily sludge

成分元素分析 wad / %工业分析 wad / %
CHO1)NS水分灰分挥发分固定碳
注:1)由差减法测出.
含油污泥12.741.648.880.310.5220.4355.4823.500.59

新窗口打开| 下载CSV


表 2   原始含油污泥的石油烃种类及质量分数

Tab.2  Type and mass fraction of petroleum hydrocarbon in original oily sludge

成分wB / %成分wB / %成分wB / %成分wB / %
2-丁烯4.44十六烷3.361-十二烯0.54二十三烷5.52
1-戊烯1.422,6,10,14-四甲基十六烷4.462,7,10-三甲基十二烷0.45二十四烷5.65
1-己烯1.00十七烷8.121-十三烯0.61二十五烷5.36
1-庚烯0.541-十七烯0.87十四烷0.80二十六烷5.00
1-辛烯0.44十八烷5.602,6,10-三甲基十五烷1.52二十八烷4.07
1-壬烯0.48十九烷5.06十五烷2.56二十九烷2.89
1-癸烯0.65二十烷15.731-十五烯1.00四十四烷4.36
1-十一烯0.60二十二烷6.323-十六碳烯0.55

新窗口打开| 下载CSV


1.2. 实验装置

阴燃实验装置如图1所示,由通气装置、加热装置、阴燃反应室、尾气冷凝和净化装置组成.

图 1

图 1   含油污泥阴燃装置的示意图

Fig.1   Schematic diagram of oily sludge smoldering device


通气装置在反应器底部,进风口与空气压缩机管路连接,气流通过带有气孔的圆板进入阴燃反应室,确保反应气流的均匀性. 加热装置用于底部物料的初步加热和点火. 阴燃反应室用于装填含油污泥,在侧面每隔3 cm共布置9个K型热电偶TC1~TC9. 尾气从阴燃反应室排出后进入冷凝和净化装置,尾气先冷凝再回收. 通过硅胶吸收水分,用活性炭吸附异味,用碱性溶液吸收SO2和H2S酸性气体,排放到大气中[19].

1.3. 实验方法

称取实验所需的含油污泥和石英砂,混合均匀后填入阴燃反应室,连接各部分装置后,打开电阻加热器持续加热. 待TC1测量值达到200 ℃后打开通气装置,调节流量计向阴燃反应室供气. 继续观察TC1所测的温度,在TC1温度到达峰值并持续下降5 min后,关闭加热装置,收集所需样品待测. 使用Agilent 5977气相色谱质谱联用仪,参照HJ 1021-2019《土壤和沉积物石油烃(C10~C40)的测定气相色谱法》测试分析具体成分[20].

1.4. 数值模型

在COMSOL中建立含油污泥阴燃的模型. 阴燃过程涉及的过程很多,包括多孔介质内的流动、传热和传质以及物质的反应化学等[21]. 含油污泥的成分复杂,实际的动力学机制难以给出. 通常采用集总反应简化,该方法仅能对阴燃过程中各类化学过程的重要性进行量化,却无法准确描述这些化学过程的具体特征[22]. 采用简单的一步氧化反应机制描述含油污泥的阴燃过程,认为氧化反应速率符合一级Arrhenius形式.

$ R=AY_{\mathrm{s}}Y_{{\mathrm{O}}_2}\exp\left(-\frac{E}{R_{\mathrm{g}}T_{\mathrm{s}}}\right). $

式中:R为反应速率,A为指前因子,E为活化能,Rg为通用气体常数,Ts为固体热力学温度,Ys为固体质量分数,$ Y_{{\mathrm{O}}_2} $为氧气质量分数。

控制方程组主要参考Zanoni等[23-24]的模型,包括固体和气体的质量守恒方程:

$ \frac{\partial Y_{\mathrm{s}}}{\partial t}=-R, $

$ \frac{\partial(\rho_{\mathrm{g}}\phi_{\mathrm{g}})}{\partial t}+\frac{\partial(\rho_{\mathrm{g}}u_{\mathrm{g}})}{\partial x}=\rho_{\mathrm{s}}R. $

式中:$ \rho\mathrm{_g} $为空气密度,$ \phi_{\mathrm{g}} $为气体孔隙率,$ u\mathrm{_g} $为入口达西空气流速,$ \rho\mathrm{_s} $为固体密度.

氧气的运输方程为

$ \phi_{\mathrm{g}}\frac{\partial(\rho_{\mathrm{g}}Y_{\mathrm{O}_2})}{\partial t}+\frac{\partial(\rho_{\mathrm{g}}u_{\mathrm{g}}Y_{\mathrm{O}_2})}{\partial x}=\phi_{\mathrm{g}}\frac{\partial}{\partial x}\left(\rho_{\mathrm{g}}D_{\mathrm{g}}\frac{\partial Y_{\mathrm{O}_2}}{\partial x}\right)-\rho_{\mathrm{s}}v_{\mathrm{O}_2}R. $

式中:Dg为氧气的扩散系数,$ v_{{\rm{O}}_2} $为氧气消耗系数.

能量守恒有局部热平衡和局部热非平衡2种形式. 当假设局部热平衡时,气体和固体的局部温度相同,此时只需要一个能量方程[25]. 对比局部热平衡和局部热非平衡的计算结果可知,虽然在峰值温度上差异很小,但在反应前沿的升温过程中,局部热非平衡模型更接近实验结果[23]. 采用局部热非平衡模型描述系统的能量守恒方程:

$ \begin{split} {(\rho {{c}_{p}})}_{{\rm{eff}}}\frac{\partial {T}_{{\rm{s}}}}{\partial t}=&\frac{\partial }{\partial x}\left({\lambda}_{{\rm{eff}}}\frac{\partial {T}_{{\rm{s}}}}{\partial x}\right)-U\left(\frac{{A}_{{\rm{s,cl}}}}{{V}_{{\rm{cl}}}}\right)({T}_{{\rm{s}}}-{T}_{0})+\\&{h}_{{\rm{sg}}}\left(\frac{{A}_{{\rm{s,sp}}}}{{V}_{{\rm{sp}}}}\right)({T}_{{\rm{g}}}-{T}_{{\rm{s}}})-{\rho }_{1}\Delta {Q}_{1}R , \end{split} $

$ \begin{split} {\phi }_{{\rm{g}}}({\rho }_{{\rm{g}}}{c}_{p{\rm{g}}})\frac{\partial {T}_{{\rm{g}}}}{\partial t}+&{\rho }_{{\rm{g}}}{c}_{p{\rm{g}}}{u}_{{\rm{g}}}\frac{\partial {T}_{{\rm{g}}}}{\partial x}={\phi }_{{\rm{g}}}\frac{\partial }{\partial x}\left({\lambda}_{{\rm{g}}}\frac{\partial {T}_{{\rm{g}}}}{\partial x}\right)+\\& {h}_{{\rm{sg}}}\left(\frac{{A}_{{\rm{s,sp}}}}{{V}_{{\rm{sp}}}}\right)({T}_{{\rm{s}}}-{T}_{{\rm{g}}}) .\end{split} $

式中:$ \left(\rho c_p\right)_{{\mathrm{e f f}}} $为有效热容量,$ \lambda_{{\mathrm{e f f}}} $为有效导热系数,$ U $为自然对换热系数,$ A_{{\mathrm{s, c l}}} $为反应容器换热表面积,$ A_{{\mathrm{s, sp}}} $为固体表面积,$ V_{{\mathrm{c l}}} $为反应容器体积,$ V_{{\mathrm{sp}}} $为固体体积,$ T_{0} $为环境温度,Tg为气体热力学温度, $ \rho_1 $为含油污泥密度,$ \Delta Q_1 $为燃烧热,$ c_{p \mathrm{g}} $为气体比定压热容,$ \lambda_{\mathrm{g}} $为气体导热系数,$ h_{{\mathrm{s g}}} $ 为气固界面传热系数。由于考虑局部热非平衡,$ {h}_{{\mathrm{sg}}} $使用文献[24]的方法计算.

模型的初始条件如下:初始时刻t0反应器内的含油污泥质量分数为1,反应器内各处温度为室温,压力为大气压. 在底部入口边界x = 0 m处,氧气质量分数恒定为23%,在t1时刻前,不通入空气,t1时刻后,开始通入空气;在t2时刻前,入口边界处持续加热,t2时刻后停止加热. 出口处的空气压力恒定为大气压.

1.5. 分子动力学计算方法

反应力场分子动力学方法(reactive force field molecular dynamics, ReaxFF MD)是结合分子力场与分子动力学模拟的计算化学方法,旨在模拟和研究化学反应的高维动力学行为[26]. 该方法的核心优势是能够精确描述分子系统中化学键的形成与断裂过程,通过优化力场参数和反应条件对反应路径进行系统性探索[27]. 为了对比实验结果,对二十烷模型在3 000、3 500和4 000 K温度下进行等温热解模拟,对分解程度进行统计,二十烷的分子结构如图2所示.

图 2

图 2   二十烷的分子结构图

Fig.2   Molecular structure diagram of eicosane


2. 结果与讨论

2.1. 含油污泥治理的效果

含油污泥阴燃后的固体残渣占产物的60%,主要为含油污泥阴燃后剩余的灰分. 液体产物包含油分和水分,分别占6.5%和17.05%,源于含油污泥中的水分和有机物反应后的产物,其余16.45%为反应生成的不凝性气体,如CO、CO2.

含油污泥的主要污染物为石油烃[28],以阴燃前、后含油污泥中石油烃总质量的变化作为无害化程度评估依据. 原始含油污泥中的石油烃质量分数为11.25%,固体残渣的石油烃质量分数为0.01%~0.13%,均低于0.3%,实现了无害化处理. 根据尾气冷凝液中回收油的质量与原始油质量分数,计算得到油回收率为36%~57.8%,资源化效果良好. 热解技术需要外部持续加热注入能量,阴燃依靠燃料的自身反应放出热量,维持反应进行,部分油分在阴燃过程中消耗,导致无法回收更多的油分.

2.2. 含油污泥阴燃特性的影响因素

2.2.1. 水的质量分数

图3可知,在达西空气流速为3 cm/s,热值为11178 kJ/kg的条件下,随着水质量分数ww的增加,含油污泥的点火时间t不断增加,近似呈线性关系,点火时间和水质量分数关系方程为$ {t}=3.27{w_{\mathrm{w}}}-4.5 $($ {R}^{2}=0.970~6 $),其中t的单位为min. 含油污泥的着火温度高于水的沸点,在加热过程中,水分蒸发的汽化潜热会消耗较多的热量,从而减缓升温速率,且水分越多消耗越大,不利于阴燃的初步点火. 虽然水质量分数的增加显著延缓了点火时间,但水质量分数为10%~50%的含油污泥阴燃均实现了完全的自维持阴燃,自维持阴燃是指含油污泥在点火后能否依靠自身反应放热在反应器内传播完全. 水分的增加必然会降低反应体系的温度,对比不同水质量分数的含油污泥的阴燃平均峰值温度θ和阴燃速率v的差异,如图4所示. 水质量分数从10%增加到50%,含油污泥阴燃平均峰值温度从707 ℃降低到了493 ℃,阴燃速率从3.69 mm/min降低到了2.60 mm/min,这与其他阴燃研究[21, 29]的结论类似.

图 3

图 3   阴燃点火时间随水质量分数的变化曲线

Fig.3   Curve of ignition time of smoldering with moisture mass fraction


图 4

图 4   平均峰值温度和阴燃速率随水质量分数的变化曲线

Fig.4   Variation curve of average peak temperature and propagation speed of smoldering with moisture mass fraction


2.2.2. 达西空气流速

不同达西空气流速下的阴燃峰值温度和阴燃速率如图5所示. 在水质量分数为20%,热值为11 178 kJ/kg的条件下,达西空气流速增加,平均阴燃峰值温度和阴燃速率与达西空气流速呈现出明显的线性关系.

图 5

图 5   平均峰值温度和阴燃速度随达西空气流速的变化曲线

Fig.5   Curve of average peak smoldering temperature and smoldering speed with Darcy air velocity variation


当达西空气流速从1 cm/s增加到5 cm/s时,平均阴燃峰值温度由630 K提升至825 K,阴燃速率从2.8 mm/min增长至3.8 mm/min. 这种正相关关系表明流速参数对反应动力学具有显著的调控作用. 达西空气流速的增加不仅通过增强氧气扩散系数加速了挥发分与固定碳的氧化反应,而且改变了多孔介质内的温度场分布特征. 高温反应区厚度随着流速的增大而扩展,这显著提高了反应前沿与未反应区域间的热量传递效率.

2.2.3. 含油污泥热值

热值直观反映了燃料燃烧时的放热能力,由于阴燃在点火后仅依靠燃料自身的反应放热维持阴燃推进,在阴燃的自维持研究中该指标有重要作用. 在水质量分数为20%,达西空气流速为3 cm/s的条件下,不同热值q的含油污泥阴燃平均峰值温度如图6所示. 随着热值的升高,平均峰值温度随之升高.

图 6

图 6   阴燃平均峰值温度与热值的关系曲线

Fig.6   Curve of relationship between smoldering average peak temperature and calorific value


热值低于2 235 kJ/kg的含油污泥均不能完全阴燃,以上述含油污泥为例,点火后仅有TC1和TC2的峰值温度超过400 ℃,其余测点温度均低于200 ℃且在点火的30 min后开始下降,阴燃逐渐熄灭. 这是由于该工况下含油污泥的热值较低,反应放热量持续小于系统的热损失,系统温度持续下降,进一步导致反应强度减弱,阴燃过程最终停止. 其余热值更高的含油污泥在点火后阴燃稳定,可以在反应器内传播完全,这些点火后可完全阴燃的污泥平均峰值温度与热值之间的关系为$ \theta =0.039\;3q+425.6 $ ($ {R}^{2}=0.935~6 $).

2.2.4. 孔隙率对阴燃特性的影响

改变含油污泥与石英砂的质量比以改变孔隙率,如图7所示为不同孔隙率下的阴燃特性差异图. 可以看出,阴燃平均峰值温度随着孔隙率的减小而升高. 这是因为低孔隙率意味着单位体积内的含油污泥更多,可燃物质量分数更高,放热量更大. 阴燃速率随着孔隙率的增大而增大,但孔隙率为25%和35%时的阴燃速率仅差0.13 mm/min,说明石英砂对含油污泥的传热和透气性能的改善作用较明显. 当孔隙率≤20%时,含油污泥占据主体,此时含油污泥几乎完全包裹石英砂,对阴燃速率的影响起主导作用,石英砂的作用较小.

图 7

图 7   阴燃平均峰值温度和传播速度随孔隙率的变化曲线

Fig.7   Curve of average peak temperature and propagation speed of smoldering with porosity variation


2.3. 含油污泥阴燃的污染物降解路径分析

含油污泥阴燃是动态过程,在不同区域可能同时存在多个过程. 例如在阴燃反应区,是富氧的,会发生氧化反应,在阴燃反应区前部,氧气量较少,可能会发生热解. 随着反应的推进,发生热解的区域在一段时间后可能会发生氧化. 除了热解和氧化,还可能有物质在反应区挥发而后在未反应区冷凝,或通入的空气带走部分未发生进一步反应的有机物. 为了深入探究阴燃过程,以原始含油污泥及阴燃产物成分和分子动力学模拟结果为依据,推断含油污泥阴燃过程的反应类型和污染物的降解路径,如图8所示.

图 8

图 8   含油污泥阴燃过程可能的污染物降解路径示意图

Fig.8   Schematic diagram of possible pollutant degradation pathway during smoldering process of oily sludge


表3所示,阴燃初期主要表现为石油烃的物理挥发:对比原始含油污泥与尾气冷凝液的GC-MS数据,检测到1-十三烯(232 ℃)、十四烷(252 ℃)及十五烷(268 ℃)等共有组分. 这些物质的沸点对应阴燃反应区200 ℃以上的高温带,处于点火完成后的快速升温阶段[14]. 空间分布分析表明,此类物质可能富集于物料表层,在阴燃锋面推进过程中受热挥发后被载气夹带逸出,未参与后续的化学反应.

表 3   含油污泥阴燃尾气冷凝液主要物质的质量分数

Tab.3  Mass fraction of major component in condensate from oil sludge smoldering exhaust

成分wB/%成分wB/%
丙烷1.50十三烷3.21
辛烷2.451-十三烯0.19
壬烷7.20十四烷1.66
癸烷10.43十五烷0.55
十一烷11.212,6-二甲基十一烷1.34
4-甲基癸烷3.56邻二甲苯2.28
十二烷6.50

新窗口打开| 下载CSV


碳数分布变化揭示显著的裂化现象:原始物料中C1~C15组分占比为15.6%,经阴燃处理后尾气中占比升至78.44%;C16+组分占比从84.4%降至21.55%. 这种碳数分布逆转(尾气冷凝液呈C11中心对称分布,原始物料以C20为主)表明,长链烷烃发生选择性裂解. 实验检测到以丙烷为主的新生短链烃证实了该推断,推测C20烷烃在本实验工况下优先裂解至C11附近即发生迁移终止.

组分结构分析显示异构化特征:尾气中检出2-甲基壬烷、3,4-二甲基庚烷及其支链异构体以及邻/对二甲苯和芳构异构体. 对比原始物料以直链烷烃为主的特征,推断高温条件(300~500 ℃)促使分子结构重排. 这类转化仅改变分子结构而未改变相对分子质量.

尾气中的环状化合物包括环烷烃(如乙基环己烷)、单环芳烃(甲苯、二甲苯)及多环芳烃(茚、萘),生成过程符合自由基反应路径:长链裂解产生的·C2H3自由基在高于500 ℃的高温区发生环化缩合,形成热力学更稳定的环状结构,该过程与文献[30]的多环芳烃生成机理一致.

质量平衡计算显示,42.2%~64%的石油烃经完全氧化生成CO2/H2O. 残渣油分检测与冷凝液完全回收的实验结果证实,氧化是阴燃过程的主要转化途径. 氧化反应贯穿于阴燃的各个阶段,与挥发、裂化过程形成竞争关系[21].

图9中,tr为反应时间,Nt为产物总数,$N_{{\mathrm{H}}_2} $为H2分子数. 如图9所示,通过对二十烷模型分子动力学模拟结果的分析,揭示了不同温度(3 000、3 500、4 000 K)下碳链产物分布及反应路径的演变规律. 如图9(a)~(c)所示分别为C1~C4、C5~C15及C16+组分的质量分数动态变化,如图9(d)所示为4 000 K温度下H2生成与产物总数的关联性.

图 9

图 9   不同温度下含碳组分质量分数以及4 000 K温度下H2数量和产物总数等温模拟变化的统计

Fig.9   Statistics of carbon-containing mass fraction at different temperature and isothermal evolution of number of H2 and total amount of product at 4 000 K


在3 000 K条件下(见图9(a)),体系在70 ps时发生初始分解,120ps时C16+与C5~C15的质量分数达到平衡. 此后C5~C15成为主导产物(>50%),显著高于C1~C4(约20%)及C16+(终值>30%). 该温度特征表明:分解反应不彻底,残留C16+组分;中长链产物(C5~C15)具有较高的热稳定性,抑制了进一步裂解为C1~C4小分子.

当温度升至3 500 K(见图9(b))时,分解起始时间提前至25 ps,C16+在80 ps后迅速耗尽(终值≈0). 在75、90、130 ps等时间节点出现C1~C4质量分数超过C5~C15的现象,表明高温显著促进C—C键断裂,但动态平衡下中长链产物保持相对优势.

4 000 K条件下的演化(见图9(c))呈现加速特征:C16+在70 ps后完全消失,C5~C15与C1~C4在75 ps后趋近平衡(质量分数均为40%). 该温度下系统产生大量的以甲烷、乙烯为主的小分子产物,其高反应活性暗示在氧化环境中可能引发CO2的快速生成. 温度从3 500 K升至4 000 K对产物分布的影响弱化,提示该温区可能接近热解反应的特征温度.

分析4 000 K体系(见图9(d))可知,产物总数随反应时间单调递增至稳定值,证实裂解反应主导于交联过程,同时H2生成量持续增长. 结合产物分布特征,表明自由基浓度的上升促进了环烃/芳烃的生成. 此反应路径可能显著提升体系在氧化条件下的H2O生成倾向.

2.4. 阴燃模拟结果
2.4.1. 达西空气流速为2 cm/s的模拟结果

当达西空气流速为2 cm/s时,含油污泥的阴燃模拟结果与实验结果的对比如图10所示. 其中,实线表示实验中记录的温度曲线,虚线表示相同工况下模拟结果的温度曲线,颜色相同的线条表示反应器内同一位置的温度曲线.

图 10

图 10   含油污泥阴燃模拟与实验结果的对比

Fig.10   Comparison between simulation and experimental result of oily sludge smoldering


对比模拟结果和实验结果可知,模拟结果基本复现了实验中的升温过程,峰值温度和阴燃速率与实验结果相近,符合较好,其中阴燃速率的结果相较于已有的含油污泥阴燃模型更优[14],误差更小.

模拟结果的峰值温度大多略高于实验结果,可能的原因是模型中仅考虑了一步简单氧化反应. 由于含油污泥阴燃过程的反应复杂,除了氧化放热,还有热解及水分蒸发的吸热过程,此类过程会降低平均峰值温度. 文献中简化为一步氧化反应的活性炭阴燃模拟温度与实验符合程度高[24],其原因是颗粒活性炭90%以上为固定碳,与纯碳发生的氧化反应更接近.

为了量化模拟结果与实验结果的误差,确定数值模型的可靠性,以相同的权重从平均峰值温度、平均阴燃速率、温度-时间曲线、温度-距离曲线4个方面计算模拟结果与实验的误差,计算公式如下:

$ \begin{split} e=&\frac{1}{4}\left| \frac{{T}_{{{\exp }}}-{T}_{{{{\mathrm{num}}}}}}{{T}_{{{{\mathrm{num}}}}}}\right|+\frac{1}{4}\left| \frac{{v}_{{{\exp }}}-{v}_{{{{\mathrm{num}}}}}}{{v}_{{{{\mathrm{num}}}}}}\right|+\\&\frac{1}{4}{\rm{NRMSD}}_{T(t)}+\frac{1}{4}{\rm{NRMSD}}_{T(x)} .\end{split} $

式中:Texp为实验平均峰值温度,Tnum为数值模拟平均峰值温度,vexp为实验平均阴燃速率,vnum为数值模拟平均阴燃速率,NRMSDT(t)为温度–时间曲线的归一化均方根偏差,NRMSDT(x)为温度–空间曲线的归一化均方根偏差. 平均归一化均方根偏差使用文献[23]的方法计算,得到平均误差约为11%.

2.4.2. 达西空气流速对模拟结果的影响

当其他条件不变,达西空气流速改变时,阴燃平均峰值温度和阴燃传播速率的模拟结果与实验结果的对比如图11所示.

图 11

图 11   达西空气流速改变时平均峰值温度和阴燃速率模拟结果和实验结果的差异

Fig.11   Difference of simulation and experimental result in average peak temperature and smoldering rate with variation of Darcy air velocity


随着达西空气流速的增加,阴燃平均峰值温度和阴燃传播速率的模拟结果随达西空气流速的变化趋势和实验结果是一致的,平均峰值温度的模拟结果的误差小于10%. 阴燃速率的模拟结果随着达西空气流速的增大而近似线性增加,但实验结果对该条件的变化更敏感,当达西空气流速增大时,实验中的阴燃速率变化更大. 其原因是模拟中的化学反应简化较多,而实际上非常复杂,在平均峰值温度较高的情况下,模拟反应速率与实际反应速率的差距较大,阴燃速率的差异稍大.

2.4.3. 热值对模拟结果的影响

图12所示为热值改变时阴燃平均峰值温度的模拟结果与实验结果的差异. 可以看出,阴燃平均峰值温度和阴燃传播速率的模拟结果与实验结果在污泥热值改变时均符合较好.

图 12

图 12   不同热值的平均峰值温度和阴燃速率模拟值和实验值的差异

Fig.12   Difference of simulation and experimental result in average peak temperature and smoldering rate of different calorific value


对于较低的热值,如2 235 kJ/kg的污泥,利用该模型可以模拟点火后不能完全传播的现象,表明该模型对预测污泥能否阴燃有一定的作用.

3. 结 论

(1)阴燃后的含油污泥中石油烃质量分数低于0.3%,污染物去除率大于99.7%,无害化程度高,油回收率可达57.8%,资源化程度良好. 水的质量分数对点火时间的影响较大. 达西空气流速的提高对阴燃平均峰值温度的影响显著. 含油污泥热值对阴燃点火的影响较小,对阴燃能否传播完全的影响较大.

(2)对比阴燃前、后产物种类和质量分数变化并结合分子动力学模拟可见,含油污泥阴燃过程发生的反应众多并且产物复杂. 主要涉及的反应有长链石油烃的裂化和短链石油烃的氧化,还包含石油烃的挥发、环化和芳构化过程,并给出可能的反应路径.

(3)建立的含油污泥阴燃数值模型模拟结果与实验结果较接近,模型的平均误差约为11%. 当达西空气流速和含油污泥热值改变时,模拟结果的阴燃平均峰值温度和阴燃传播速率变化趋势与实验结果一致,结果符合较好. 含油污泥的阴燃传播速率对达西空气流速的敏感性较高,当达西空气流速增大时,模拟结果的误差稍大.

参考文献

LI H, SONG S F, QU C T, et al. Study of treatment process and development of oily sludge [J]. Advanced Materials Research, 2013, 838/839/840/841: 2667–2672.

[本文引用: 1]

王玉华, 陈传帅, 孟娟, 等

含油污泥处置技术的新发展及其应用现状

[J]. 安全与环境工程, 2018, 25 (3): 103- 110

DOI:10.13578/j.cnki.issn.1671-1556.2018.03.017      [本文引用: 1]

SILVA L J, ALVES F C, DE FRANÇA F P

A review of the technological solutions for the treatment of oily sludges from petroleum refineries

[J]. Waste Management and Research, 2012, 30 (10): 1016- 1030

DOI:10.1177/0734242X12448517      [本文引用: 1]

HU G, LI J, ZENG G

Recent development in the treatment of oily sludge from petroleum industry: a review

[J]. Journal of Hazardous Materials, 2013, 261: 470- 490

DOI:10.1016/j.jhazmat.2013.07.069      [本文引用: 1]

WANG Yuhua, CHEN Chuanshuai, MENG Juan, et al

Development and application of disposal techniques on oil sludge

[J]. Safety and Environmental Engineering, 2018, 25 (3): 103- 110

DOI:10.13578/j.cnki.issn.1671-1556.2018.03.017      [本文引用: 1]

ZORIN A D, KARATAEV E N, ZANOZINA V F, et al

Thin-film cracking of acid tars

[J]. Petroleum Chemistry, 2012, 52 (4): 245- 252

DOI:10.1134/S0965544112040135      [本文引用: 1]

陈晓丽, 李龙, 雷大鹏, 等

阴燃技术在土壤修复中应用的中试

[J]. 化工设计通讯, 2019, 45 (10): 37- 38

DOI:10.3969/j.issn.1003-6490.2019.10.027      [本文引用: 1]

CHEN Xiaoli, LI Long, LEI Dapeng, et al

Pilot test of smoldering technology applied in soil remediation

[J]. Chemical Engineering Design Communications, 2019, 45 (10): 37- 38

DOI:10.3969/j.issn.1003-6490.2019.10.027      [本文引用: 1]

葛传芹, 雷大鹏, 刘杉, 等

采用异位阴燃修复技术处理含油固废中试实验

[J]. 环境工程学报, 2022, 16 (2): 601- 611

DOI:10.12030/j.cjee.202108168      [本文引用: 1]

GE Chuanqin, LEI Dapeng, LIU Shan, et al

Pilot studies on treatment of oily solid waste with ex situ smoldering remediation technology

[J]. Chinese Journal of Environmental Engineering, 2022, 16 (2): 601- 611

DOI:10.12030/j.cjee.202108168      [本文引用: 1]

王天宇, 蒋文明, 刘杨

含油污泥阴燃处理技术研究与进展

[J]. 化工学报, 2020, 71 (4): 1411- 1423

DOI:10.11949/0438-1157.20191085      [本文引用: 1]

WANG Tianyu, JIANG Wenming, LIU Yang

Research and progress of smoldering combustion technology for oily sludge

[J]. Journal of Chemical Industry and Engineering(China), 2020, 71 (4): 1411- 1423

DOI:10.11949/0438-1157.20191085      [本文引用: 1]

SOLINGER R, GRANT G P, SCHOLES G C, et al

STARx Hottpad for smoldering treatment of waste oil sludge: proof of concept and sensitivity to key design parameters

[J]. Waste Management and Research, 2020, 38 (5): 554- 566

[本文引用: 1]

ZHAO C, LI Y, GAN Z, et al

Method of smoldering combustion for refinery oil sludge treatment

[J]. Journal of Hazardous Materials, 2021, 409: 124995

DOI:10.1016/j.jhazmat.2020.124995      [本文引用: 1]

杨高玄. 含油污泥的阴燃处理技术研究 [D]. 大连: 大连理工大学, 2021.

[本文引用: 1]

YANG Gaoxuan. Research on smoldering treatment technology of oil-contaminated sludge [D]. Dalian: Dalian University of Technology, 2021.

[本文引用: 1]

杨高玄, 占敬敬

阴燃处理华北油田含油污泥的研究

[J]. 应用化工, 2022, 51 (2): 317- 321

DOI:10.3969/j.issn.1671-3206.2022.02.005      [本文引用: 1]

YANG Gaoxuan, ZHAN Jingjing

Study on the treatment of oily sludge in North China oilfield by smoldering

[J]. Applied Chemical Industry, 2022, 51 (2): 317- 321

DOI:10.3969/j.issn.1671-3206.2022.02.005      [本文引用: 1]

GAN Z, ZHAO C, LI Y, et al

Experimental investigation on smoldering combustion for oil sludge treatment: influence of key parameters and product analysis

[J]. Fuel, 2022, 316: 123354

DOI:10.1016/j.fuel.2022.123354      [本文引用: 1]

甘宗伟. 阴燃法治理含油污泥关键特性实验及模拟研究 [D]. 济南: 山东大学, 2022.

[本文引用: 4]

GAN Zongwei. Experimental and simulation study on key characteristics of smoldering method for oily sludge treatment [D]. Jinan: Shandong University, 2022.

[本文引用: 4]

雷大鹏, 单晖峰, 杨登, 等

工程化阴燃技术治理含油污泥工程示范

[J]. 环境工程, 2022, 40 (10): 150- 155

DOI:10.13205/j.hjgc.202210020      [本文引用: 1]

LEI Dapeng, SHAN Huifeng, YANG Deng, et al

A demonstration project of oil sludge treatment with engineered smoldering technology

[J]. Environmental Engineering, 2022, 40 (10): 150- 155

DOI:10.13205/j.hjgc.202210020      [本文引用: 1]

苗芳. 含油污泥中重质组分在热解过程中的演化规律 [D]. 西安: 西安石油大学, 2023: 65.

[本文引用: 1]

MIAO Fang. Evolution patterns of heavy components in oily sludge during pyrolysis process [D]. Xi’an: Xi’an Shiyou University, 2023.

[本文引用: 1]

HEIDARZADEH N, GITIPOUR S, ALI ABDOLI M

Characterization of oily sludge from a Tehran oil refinery

[J]. Waste Management and Research, 2010, 28 (10): 921- 927

DOI:10.1177/0734242X09345794      [本文引用: 1]

CAUMETTE G, LIENEMANN C, MERDRIGNAC I, et al

Element speciation analysis of petroleum and related materials

[J]. Journal of Analytical Atomic Spectrometry, 2009, 24 (3): 263- 276

DOI:10.1039/b817888g      [本文引用: 1]

YU H, SHAN C, LI J, et al

Alkaline absorbents for SO2 and SO3 removal: a comprehensive review

[J]. Journal of Environmental Management, 2024, 366: 121532

DOI:10.1016/j.jenvman.2024.121532      [本文引用: 1]

LI T, CAO X, HE S, et al

An accelerated solvent extraction and gas chromatography-flame ionization detector method to rapidly determining and assessing total petroleum hydrocarbon contamination in soil from Fushan oilfield, China

[J]. Environmental Science and Pollution Research, 2020, 27 (30): 37444- 37454

DOI:10.1007/s11356-020-09418-4      [本文引用: 1]

YERMÁN L, WALL H, TORERO J L

Experimental investigation on the destruction rates of organic waste with high moisture content by means of self-sustained smoldering combustion

[J]. Proceedings of the Combustion Institute, 2017, 36 (3): 4419- 4426

DOI:10.1016/j.proci.2016.07.052      [本文引用: 3]

杨玖玲. 泥炭阴燃及阴燃气体生成规律的实验与机理研究 [D]. 合肥: 中国科学技术大学, 2017.

[本文引用: 1]

YANG Jiuling. Experimental and mechanistic study on peat smoldering and the generation law of smoldering gases [D]. Hefei: University of Science and Technology of China, 2017.

[本文引用: 1]

ZANONI M A B, TORERO J L, GERHARD J I

Determination of the interfacial heat transfer coefficient between forced air and sand at Reynold’s numbers relevant to smouldering combustion

[J]. International Journal of Heat and Mass Transfer, 2017, 114: 90- 104

DOI:10.1016/j.ijheatmasstransfer.2017.06.020      [本文引用: 3]

ZANONI M A B, WANG J, GERHARD J I

Understanding pressure changes in smouldering thermal porous media reactors

[J]. Chemical Engineering Journal, 2021, 412: 128642

DOI:10.1016/j.cej.2021.128642      [本文引用: 3]

TORERO J L, GERHARD J I, MARTINS M F, et al

Processes defining smouldering combustion: integrated review and synthesis

[J]. Progress in Energy and Combustion Science, 2020, 81: 100869

DOI:10.1016/j.pecs.2020.100869      [本文引用: 1]

LI X, ZHENG M, REN C, et al

ReaxFF molecular dynamics simulations of thermal reactivity of various fuels in pyrolysis and combustion

[J]. Energy and Fuels, 2021, 35 (15): 11707- 11739

[本文引用: 1]

LIU Q, LIU S, LV Y, et al

Atomic-scale insight into the pyrolysis of polycarbonate by ReaxFF-based reactive molecular dynamics simulation

[J]. Fuel, 2021, 287: 119484

DOI:10.1016/j.fuel.2020.119484      [本文引用: 1]

EGAZAR’YANTS S V, VINOKUROV V A, VUTOLKINA A V, et al

Oil sludge treatment processes

[J]. Chemistry and Technology of Fuels and Oils, 2015, 51 (5): 506- 515

DOI:10.1007/s10553-015-0632-7      [本文引用: 1]

YERMÁN L, HADDEN R M, CARRASCAL J, et al

Smouldering combustion as a treatment technology for faeces: exploring the parameter space

[J]. Fuel, 2015, 147: 108- 116

DOI:10.1016/j.fuel.2015.01.055      [本文引用: 1]

RAVINDRA K, SOKHI R, VAN GRIEKEN R

Atmospheric polycyclic aromatic hydrocarbons: source attribution, emission factors and regulation

[J]. Atmospheric Environment, 2008, 42 (13): 2895- 2921

DOI:10.1016/j.atmosenv.2007.12.010      [本文引用: 1]

/