CN104849260A - 一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法 - Google Patents

Abstract

本发明公开了一种直接采用火焰发射光谱分析来在线检测锅炉火焰中气相碱金属浓度的方法，首先选择要测量碱金属的特征谱线，然后构建特征谱线辐射强度关于火焰温度、气相碱金属浓度的拟合模型，再次通过实验标定确定拟合系数，最后由便携式光纤光谱仪测量火焰对象的辐射光谱，得到火焰温度和碱金属特征谱线辐射强度，代入建立的已知拟合系数的拟合模型计算得到火焰中气相碱金属的浓度。本发明能够简便、准确的在线检测炉内燃烧火焰中气相碱金属的浓度，并可以同时检测多种碱金属含量，设备简单、成本低廉，便于工业现场测量。

Description

一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法

技术领域

本发明属于火焰发射光谱检测技术领域，更具体地，涉及一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法。

背景技术

随着国家发展，用电量逐渐增加，能源消耗大幅增长，同时环境问题日益突出，使得电厂不得不燃用新煤种，并开发新型燃烧方式来应对能源与环境的双重压力。近年发现的准东煤储量大、煤质较好、灰分含量低，是良好的动力煤，但因其碱金属含量较高，导致炉膛结渣及受热面积灰等问题，成为研究利用的难点。林木废弃物、农业废弃物、水生植物、油料植物、有机物加工废料、人畜粪便及城市生活垃圾等生物质能被认为是一种具有应用前景的可再生能源，然而与常规化石燃料相比，生物质燃料中碱金属含量普遍较高，燃烧过程中极易进入气相，引起高温腐蚀、受热面结垢、炉内结渣等问题。因此，准确测量燃烧火焰中气相碱金属浓度对燃高碱煤和生物质锅炉的安全运行有着十分重要的意义，同时对研究碱金属释放机理及模型建立也有重要的指导作用。

以CCD阵列作为光谱测量元件的光纤光谱仪可以检测到燃烧火焰在不同波长下的光谱辐射信息，根据National Institute of Standards andTechnology(NIST)Atomic Spectra Database可知，主要碱金属的原子特征谱线为：锂Li(670.776nm、670.791nm)，钠Na(588.995nm、589.592nm)，钾K(766.490nm、769.896nm)和铷Rb(780.027nm、794.760nm)，这为碱金属的定性和定量分析提供了条件。传统碱金属检测方法主要是离线检测，如火焰原子吸收光谱法、火焰光度法、电感耦合等离子体原子发射光谱法、石墨炉原子吸收光谱法等，但这些方法都需要对采集的样品进行前处理，以达到实验条件，操作复杂，而且这种离线检测方法无法获取燃烧过程中炉内气相碱金属分布规律。随着激光技术的发展，研究者利用激光诱导击穿光谱技术Laser-Induced Breakdown Spectroscopy(LIBS)成功测量了褐煤和松木燃烧时气体中碱金属Na、K含量。但对于大型工业炉来说，LIBS在线检测技术也有明显的缺点：设备结构复杂，操作繁琐，成本高昂，并且LIBS检测没有考虑温度变化，要维持温度恒定，测量条件要求苛刻，不便于工业推广。

发明内容

针对现有技术的上述缺点和/或改进需求，本发明提供了一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法，其中根据锅炉炉膛中燃烧火焰的特点，直接采用火焰发射光谱分析在线检测锅炉火焰中气相碱金属的浓度，相应的可有效解决现有检测方法难以直接应用于大型工业锅炉中气相碱金属浓度检测的问题，具有检测简便、准确、设备简单、成本低廉等优点，因而尤其适用于大型工业锅炉中气相碱金属浓度的在线检测。

为实现上述目的，本发明提出了一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法，其特征在于，包括如下步骤：

(1)基于特征谱线构建拟合模型：选取要测量的碱金属特征谱线，基于上述选取的碱金属特征谱线，构建燃烧火焰中碱金属特征谱线辐射强度I_Alkali关于气相碱金属浓度C_Alkali与火焰温度T之间的拟合模型：其中a_mn是多项式拟合系数，阶数m、n为正整数，T表示作为测量对象的火焰的温度，其单位为K；

(2)执行标定处理：将检测光谱仪进行绝对辐射强度标定，然后用该标定后的检测光谱仪测量已知碱金属浓度的火焰对象，得到其碱金属特征谱线辐射强度和火焰温度，把测量的一系列不同碱金属浓度下对应的特征谱线辐射强度和火焰温度带入所述拟合模型中，得到模型的多项式拟合系数a_mn的具体数值；

(3)在线检测：采用上述标定后的检测光谱仪实时测量锅炉燃烧火焰辐射光谱，获得相对应的碱金属特征谱线辐射强度，同时获得火焰温度，将所述碱金属特征谱线辐射强度、所述火焰温度以及步骤(2)获得的多项式拟合系数a_mn的具体数值代入步骤(1)中所构建的拟合模型中，计算获得锅炉燃烧火焰中气相碱金属的实际浓度，以此实现锅炉燃烧火焰中气相碱金属浓度的在线检测。

作为进一步优选的，步骤(1)中所述拟合模型的构建步骤如下：

(1.1)根据火焰光谱分析，在热力学平衡时，元素的原子或离子吸收高温火焰中的能量由基态变为激发态，单位体积的基态原子数N₀与激发态原子数N_i之间的分布满足统计力学中麦克斯韦—玻尔兹曼分布定律，即：

$N_i=\frac{g_i}{g_0}\·N_0\·e^{-\frac{E_i}{\mathrm{kT}}}---\left(1\right);$

其中，g_i、g₀分别为激发态与基态的统计权重，E_i为激发能，k为玻尔兹曼常数，T为激发温度。

(1.2)两能级i、j之间跃迁所产生的谱线强度I_ij与激发态原子数目的关系表示如下：

I_ij＝N_i·A_ij·h·v_ij    (2)；

其中，A_ij为两个能级间的跃迁几率，h为普朗克常数，v_ij为发射谱线的频率；将上述式(1)带入该式(2)中得到下式：

$I_{\mathrm{ij}}=\frac{g_i}{g_0}\·A_{\mathrm{ij}}\·h\·v_{\mathrm{ij}}\·N_0\·e^{-\frac{E_i}{\mathrm{kT}}}---\left(3\right);$

(1.3)根据式(3)可知，原子特征谱线辐射强度与基态原子数N₀和激发温度T有关，而基态原子数N₀又与浓度有关，因此碱金属特征谱线辐射强度I_Alkali可以表示为碱金属浓度C_Alkali与火焰温度T的函数，采用最小二乘法拟合得到拟合模型：

$I_{\mathrm{Alkali}}=\underset{0}{\overset{m}{\mathrm{\Σ}}}\underset{0}{\overset{n}{\mathrm{\Σ}}}{a}_{\mathrm{mn}}{\left(C_{\mathrm{Alkali}}\right)}^m{T}^n.$

作为进一步优选的，所述检测光谱仪在黑体炉上进行标定。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，主要具备以下的技术优点：

1.本发明基于对含碱金属火焰在可见光范围内发射光谱的分析，构建了碱金属特征谱线辐射强度关于火焰中气相碱金属浓度、火焰温度三者之间的拟合模型，建立了一种简单有效的燃烧火焰中气相碱金属浓度检测方法，从而获得炉内气相碱金属分布规律，对于分析预防碱金属腐蚀和结渣具有重要的意义，且本发明属于非接触式测量方法，对测量火焰对象无干扰。

2.本发明的检测方法能够简便、准确的在线检测炉内燃烧火焰中气相碱金属的浓度，并可以同时在线检测多种碱金属的含量，设备简单、成本低廉，便于工业现场测量，可直接应用于大型工业锅炉中气相碱金属浓度的检测，实现炉内含碱金属燃料燃烧过程中气相碱金属浓度的在线测量。

附图说明

图1是本发明的流程框图；

图2是实验标定处理系统图；

图3是垃圾燃烧火焰在可见光波段的发射谱线。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

图1给出了本发明的流程框图，所述燃烧火焰中气相碱金属浓度的在线检测方法包括如下步骤：

(1)选取特征谱线并构建拟合模型步骤

正常状态下，元素处于基态，当元素的原子或离子吸收高温火焰中的能量时，被激活成激发态(高能态)，由于激发态极不稳定，会跃迁到低能态或基态，同时释放多余能量，从而发射出一系列特征光谱线，由最低激发态返回到基态所发射的谱线称为第一共振线，它是众多光谱线中最强的，也是元素分析中最灵敏的谱线。根据NIST Atomic Spectra Database数据库可知，每种碱金属都有自己的原子特征谱线，其第一共振线中强度较强的特征谱线分别为：Li(670.776nm、670.791nm)，Na(588.995nm、589.592nm)，K(766.490nm、769.896nm)和Rb(780.027nm、794.760nm)，这就是理论中要选取的碱金属特征谱线。实际测量时，由于光谱仪的分辨率不同，选择特征谱线时还要结合所用光谱仪能分辨的最小波长选取。

本发明用AvaSpec-2048-USB2光谱仪去测量垃圾焚烧炉中火焰的辐射光谱，由于该光谱仪的分辨率为0.8nm，实际检测到的碱金属原子特征谱线分别为锂Li(670.451nm)，钠Na(589.139nm)，钾K(766.538nm、769.886nm)和铷Rb(779.920nm、794.388nm)，其中Na(589.139nm)和K(766.538nm、769.886nm)的谱线最明显。

理论上上述检测到的特征谱线都可以用于定量测量，但考虑到会有其他元素的特征谱线干扰时，如氩Ar会在772.376nm和772.421nm出现特征谱线，因此K、Rb元素特征谱线的选取应尽量避开其干扰。因此，针对本发明碱金属特征谱线选取为：锂Li(670.451nm)、钠Na(589.139nm)、钾K(766.538nm)和铷Rb(794.388nm)。

基于所选择的特征谱线构建火焰中碱金属特征谱线辐射强度I_Alkali关于气相碱金属浓度C_Alkali与火焰温度T之间的拟合模型：

(1.1)根据火焰光谱分析，在热力学平衡时，元素的原子或离子吸收高温火焰中的能量，由基态变为激发态，单位体积的基态原子数N₀与激发态原子数N_i之间的分布满足统计力学中麦克斯韦—玻尔兹曼分布定律，即：

$N_i=\frac{g_i}{g_0}\·N_0\·e^{-\frac{E_i}{\mathrm{kT}}}---\left(1\right)$

其中，g_i、g₀分别为激发态与基态的统计权重；E_i为激发能；k为玻尔兹曼常数；T为激发温度。

(1.2)而i、j两能级之间跃迁所产生的谱线强度I_ij与激发态原子数目成正比，可以表示为：

I_ij＝N_i·A_ij·h·v_ij    (2)

其中，h为普朗克常数；A_ij为两个能级间的跃迁几率；v_ij为发射谱线的频率；式(1)带入式(2)得：

$I_{\mathrm{ij}}=\frac{g_i}{g_0}\·A_{\mathrm{ij}}\·h\·v_{\mathrm{ij}}\·N_0\·e^{-\frac{E_i}{\mathrm{kT}}}---\left(3\right)$

(1.3)根据式(3)可知，原子特征谱线辐射强度与基态原子数N₀和激发温度T有关，而基态原子数N₀又与浓度有关，因此碱金属特征谱线辐射强度I_Alkali可以表示为碱金属浓度C_Alkali与火焰温度T的函数，采用最小二乘法拟合得到形如式(4)的二元多项式函数：

$I_{\mathrm{Alkali}}=\underset{0}{\overset{m}{\mathrm{\Σ}}}\underset{0}{\overset{n}{\mathrm{\Σ}}}{a}_{\mathrm{mn}}{\left(C_{\mathrm{Alkali}}\right)}^m{T}^n---\left(4\right)$

其中，a_mn是多项式拟合系数，阶数m和n根据需要选取，使得拟合误差足够小(1％-5％)即可，此二元多项式函数即为所构建的拟合模型。

在本实施例中，所述阶数m、n最大为4，将上述数值代入式(4)中，得到拟合关系式为：

I_Alkali＝a₀₀+a₀₁*T+a₀₂*T²+a₀₃*T³+a₀₄*T⁴+a₁₀*C_Alkali+a₁₁*C_Alkali*T

+a₁₂*C_Alkali*T²+a₁₃*C_Alkali*T³+a₂₀*C_Alkali ²+a₂₁*C_Alkali ²*T      (5)

+a₂₂*C_Alkali ²*T²+a₃₀*C_Alkali ³+a₃₁*C_Alkali ³*T+a₄₀*C_Alkali ⁴

其中，碱金属特征谱线辐射强度I_Alkali的单位为W/m³/sr，温度T的单位为K，气相碱金属浓度C_Alkali的单位为ppm。

(2)标定处理步骤

为了保证进入光纤，并被光谱仪接收的信号全部是来自测量位置视线累积方向上的火焰辐射信号，在光纤探测一端安装准直透镜；光谱仪的原始输出信号是一个通过光电转换得到的电压值，不能完全反映相应辐射能量的大小，为了得到光谱绝对辐射强度，光谱仪系统在黑体炉上进行标定，用进行过绝对辐射强度标定的光谱仪系统测量已知碱金属浓度的火焰对象，得到其碱金属特征谱线辐射强度和火焰温度，通过实验标定，把测量的一系列不同碱金属浓度下对应的特征谱线辐射强度和火焰温度带入式(5)中，得到模型的多项式拟合系数a_mn。

模型多项式拟合系数a_mn的确定是在图2所示的实验标定处理系统上进行的，主要设备包括层流预混火焰多射流燃烧器、燃料及氧化剂供给装置、空气供给装置、质量流量控制器、光谱检测系统等。乙烯、氧气、氮气分别从钢瓶中流出，通过调压阀调整压力，然后进入质量流量控制器1，通过四通阀2进入气体混合器3；氩气从钢瓶中流出后依次经过调压阀和质量流量计，并且携带超声波雾化器4中生成的碱金属盐液雾，进入气体混合器3，四种气体与碱金属液雾混合均匀后一并进入燃烧器5的射流管路，伴随气流由空气压缩机6供给，进入质量流量控制器，最后进入射流燃烧器5的伴随气流管路。光谱检测系统主要有准直透镜7、光纤8、光谱仪9、USB连接线10及笔记本电脑11组成。表1给出了实验标定的具体工况：

表1实验工况设定

实验过程具体如下：

(2.1)计算气相碱金属浓度C_Alkali和碱金属特征谱线辐射强度I_Alkali：

在一个工况下，保持燃料与氧化剂流量恒定即维持一个恒定的火焰温度，改变碱金属溶液的浓度C_s，得到一系列火焰中气相碱金属浓度C_Alkali(0ppm，2ppm，4ppm，6ppm，8ppm，10ppm)以及其对应的碱金属特征谱线辐射强度值I_Alkali；

对于碱金属在气相中浓度的计算，由于雾化碱金属溶液是与燃气充分混合后再进入射流管，因此可以假设火焰中碱金属分布是均匀的。碱金属溶液的溶质质量流量ν_m(g/s)通过式(6)计算得到：

ν_m＝C_s·ν_s    (6)

其中，C_s为碱金属溶液的浓度(％)，ν_s为碱金属溶液的消耗速率(g/s)，可以通过让雾化器长时间工作，统计平均得到；

火焰中碱金属的摩尔流量ν_ms(mol/s)通过式(7)计算得到：

${\mathrm{\ν}}_{\mathrm{ms}}=a_m\·v_m/M=\frac{a_m\·C_s\·{\mathrm{\ν}}_s}{M}---\left(7\right)$

其中，a_m为碱金属盐的摩尔分数(％)；M为碱金属盐的摩尔质量(g/mol)；

火焰区域总气体流量v_f(L/s)可以表示为：

$v_f=(v_{o_2}+v_{N_2}+v_{C_2{H}_4})\·T_f/T_r---\left(8\right)$

其中，和分别为氧气、氮气和乙烯的体积流量，通过流量控制器控制得到；T_f为火焰平均温度K；T_r为实验室房室温K；

根据理想气体的状态方程，火焰区域总气体的摩尔流量ν_mf(mol/s)可以表示为：

${\mathrm{\ν}}_{\mathrm{mf}}=\frac{P\·v_f}{R\·T_f}---\left(9\right)$

其中，R为理想气体常数，P为压力；

将式(8)带入式(9)整理得：

${\mathrm{\ν}}_{\mathrm{mf}}=\frac{P\·(v_{o_2}+v_{N_2}+v_{C_2{H}_4})}{R\·T_r}---\left(10\right)$

将式(7)与式(10)做比值，即可得到火焰区域气相碱金属的体积浓度C_Alkali(ppm)：

$C_{\mathrm{Alkali}}=v_{\mathrm{ms}}/v_{\mathrm{mf}}=\frac{a_m\·C_s\·{\mathrm{\ν}}_s}{M}\·\frac{R\·T_r}{P\·(v_{o_2}+v_{N_2}+v_{C_2{H}_4})}---\left(11\right)$

至此，在已知碱金属溶液浓度、设定气体流量、已知室温的前提下，可以通过式(11)计算得到火焰中气相碱金属的体积浓度；

然后，改变工况，即营造不同温度环境，重复上述过程，从标定过的光谱仪中测量得到不同火焰温度下的碱金属特征谱线辐射强度值I_Alkali；

(2.2)计算火焰温度：

对于火焰温度的计算，利用热辐射定律，对经过黑体炉标定的绝对辐射强度I_d(λ)，选取间隔为Δλ的两波长λ、λ+Δλ，计算得到火焰温度T_λ：

$T_{\mathrm{\λ}}=-C_2(\frac{1}{\mathrm{\λ}}-\frac{1}{\mathrm{\λ}+\mathrm{\Δ\λ}})/\ln \left(\frac{I_d\left(\mathrm{\λ}\right)}{I_d(\mathrm{\λ}+\mathrm{\Δ\λ})}\frac{{\mathrm{\λ}}^5}{\mathrm{\λ}+\mathrm{\Δ\λ}}\right)---\left(12\right)$

其中，Planck常数C₂＝1.4388E-2；I_d(λ)、I_d(λ+Δλ)分别为波长λ、λ+Δλ下的光谱辐射强度。

(2.3)计算拟合系数a_mn：

把表1中五个实验工况下的不同火焰温度和不同气相碱金属浓度及对应的碱金属特征谱线辐射强度值，代入式(5)中建立的拟合关系式，计算得到拟合系数a_mn，至此，确立碱金属特征谱线辐射强度关于气相碱金属浓度和火焰温度之间的函数模型。

(3)在线检测步骤

用标定过的光谱仪系统去测量炉内燃烧火焰辐射光谱，利用热辐射定律，计算得到火焰温度，提取火焰辐射光谱中要测量的碱金属特征谱线辐射强度值，将计算的火焰温度和特征谱线辐射强度值代入已知拟合系数的式(5)中，求解出燃烧火焰中气相碱金属浓度。

图3是AvaSpec-2048-USB2便携式光谱仪在垃圾焚烧锅炉上检测得到的火焰辐射光谱图。

该光谱仪探测器是2048像素的CCD阵列，测量波长范围为200～1100nm，分辨率0.8nm，信噪比为200：1，根据炉内实际测量得到的火焰辐射光谱，计算得到火焰温度T,把要测量的碱金属特征谱线辐射强度I_Alkali和火焰温度T，代入已知拟合系数a_mn的式(5)中，即可得到该测量火焰中气相碱金属的具体浓度。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法，其特征在于，包括如下步骤：

(1)基于特征谱线构建拟合模型：选取要测量的碱金属特征谱线，基于上述选取的碱金属特征谱线，构建燃烧火焰中碱金属特征谱线辐射强度I_Alkali关于气相碱金属浓度C_Alkali与火焰温度T之间的拟合模型：其中a_mn是多项式拟合系数，阶数m、n为正整数，T表示作为测量对象的火焰的温度，其单位为K；

(2)执行标定处理：将检测光谱仪进行绝对辐射强度标定，然后用该标定后的检测光谱仪测量已知碱金属浓度的火焰对象，得到其碱金属特征谱线辐射强度和火焰温度，把测量的一系列不同碱金属浓度下对应的特征谱线辐射强度和火焰温度带入所述拟合模型中，得到模型的多项式拟合系数a_mn的具体数值；

(3)在线检测：采用上述标定后的检测光谱仪实时测量锅炉燃烧火焰辐射光谱，获得相对应的碱金属特征谱线辐射强度，同时获得火焰温度，将所述碱金属特征谱线辐射强度、所述火焰温度以及步骤(2)获得的多项式拟合系数a_mn的具体数值代入步骤(1)中所构建的拟合模型中，计算获得锅炉燃烧火焰中气相碱金属的实际浓度，以此实现锅炉燃烧火焰中气相碱金属浓度的在线检测。

2.如权利要求1所述的一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法，其特征在于，步骤(1)中所述拟合模型的构建步骤如下：

(1.1)根据火焰光谱分析，在热力学平衡时，元素的原子或离子吸收高温火焰中的能量由基态变为激发态，单位体积的基态原子数N₀与激发态原子数N_i之间的分布满足统计力学中麦克斯韦—玻尔兹曼分布定律，即：

$N_i=\frac{g_i}{g_0}\·N_0\·e^{-\frac{E_i}{\mathrm{kT}}}---\left(1\right);$

其中，g_i、g₀分别为激发态与基态的统计权重，E_i为激发能，k为玻尔兹曼常数，T为激发温度。

(1.2)两能级i、j之间跃迁所产生的谱线强度I_ij与激发态原子数目的关系表示如下：

I_ij＝N_i·A_ij·h·v_ij   (2)；

其中，A_ij为两个能级间的跃迁几率，h为普朗克常数，v_ij为发射谱线的频率；将上述式(1)带入该式(2)中得到下式：

$I_{\mathrm{ij}}=\frac{g_i}{g_0}\·A_{\mathrm{ij}}\·h\·v_{\mathrm{ij}}\·N_0\·e^{-\frac{E_i}{\mathrm{kT}}}---\left(3\right);$

(1.3)根据式(3)可知，原子特征谱线辐射强度与基态原子数N₀和激发温度T有关，而基态原子数N₀又与浓度有关，因此碱金属特征谱线辐射强度I_Alkali可以表示为碱金属浓度C_Alkali与火焰温度T的函数，采用最小二乘法拟合得到拟合模型：

$I_{\mathrm{Alkali}}=\underset{0}{\overset{m}{\mathrm{\Σ}}}\underset{0}{\overset{n}{\mathrm{\Σ}}}{a}_{\mathrm{mn}}{\left(C_{\mathrm{Alkali}}\right)}^m{T}^n.$

3.如权利要求1或2所述的一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法，其特征在于，所述检测光谱仪在黑体炉上进行标定。