CN104315535A - 一种火焰燃烧状态检测装置及检测方法 - Google Patents

Abstract

本发明公开了一种火焰燃烧状态检测装置及检测方法，包括：光纤的一端安装于炉膛内，另一端与光电转换电路连接；光电转换电路通过光敏电阻将红外辐射强度转换为电流信号以实现远距离传输；所述光电转换电路与智能信号处理器连接，智能信号处理器通过网络接口与后台监控主机连接；智能信号处理器用于实现信号的采样、A/D转换以及数字滤波，所述智能信号处理器的数量根据炉膛中喷燃器的数量确定。采用本发明提出的火焰检测方法解决了目前火检装置领域长期存在的技术难题，为火焰的准确检测提供了理论依据。

Description

一种火焰燃烧状态检测装置及检测方法

技术领域

本发明涉及一种火焰燃烧状态检测装置及检测方法。

背景技术

在火力发电领域，煤碳发电有两种燃烧方式，一种为煤粉燃烧方式，另一种为循环流化床燃烧方式。前者，煤炭被磨煤机磨成粉状后通过送风机的强大风力被送进炉膛燃烧，目前绝大多数的火力发电机组为煤粉方式。后者主要是近十几年来为了环保的需要，对一些含有大量煤矸石的劣质煤炭进行处理而出现的一种非主力发电方式。对于运行中的煤粉燃烧锅炉来说，每时每刻都有大量的煤粉通过布置在多层、多角(墙)的喷燃器将大量的煤粉吹向炉膛。由于煤粉是一种高温易爆物质，煤粉如在炉膛内不能及时燃烧，将带来爆燃或爆炸的危险。因此，需要对每个喷燃器的燃烧状况进行准确检测，以便在喷燃器火焰熄灭后及时切断燃料，以保障炉膛安全。这种检测火焰燃烧状况的设备在行业内简称为火检。

但纵观火检设备发展历史，无论是300MW以下的小机组还是目前的1000MW机组，火检装置一直以来都由国外公司垄断。目前国内火检市场主要由ABB和FORNEY两家公司垄断。虽然也有哈尔滨、大连、上海、徐州等几家国内厂家，但由于技术和品牌无法与ABB和FORNEY相抗衡，因此，只在一些小机组中占有一定份额。但从技术的角度，包括ABB和FORNEY的产品也并非无懈可击。目前火检装置普遍存在一个技术难题，就是，行业中习惯地称之为“偷看”问题。所谓“偷看”，就是，随着煤种或燃烧状况的变化，虽然被观测的喷燃器的火焰强度已经很低了，甚至已经灭火，但此时火检检测的信号因未低于设定好的灭火门槛值，因此给出的信号依然是“有火”。但抬高门槛值往往又意味着频繁发生“灭火”误报警。因此，门槛值的整定很伤脑筋，解决的办法是引进给粉机运转状态等外部信号。这事实上已说明“偷看”确实已是一个老大难问题。“偷看”具有很大的危害性和危险性，必须彻底解决。

从火焰检测原理上看，由于炉膛中红外光约占90％，可见光约10％，因此，目前火检以检测红外光为主。红外光的检测一般采用硫化铅光敏电阻作为传感器。红外光辐射强度的变化通过光敏电阻再转换为电信号，最后由硬件电路或计算机检测处理。从信号特征上看，燃烧中的火焰既有直流成分，也有交流成分。交流成分实际反映了火焰的闪烁状态，因此，能否获取闪烁特征是火焰状态检测成败的关键因素。FORNEY对火焰的闪烁特征采用了的硬件带通滤波器，带通范围可设定，之后以带通滤波器之后的信号总强度与设定的门槛值相比较作为灭火条件。FORNEY的这一方法，虽然获得了闪烁特征，但阻断了直流成分，无法全面掌握火焰燃烧状况。ABB的主流产品810系列采用了模拟信号直接采样，配合FPGA数字滤波技术和傅里叶滤波技术，既可以获得直流成分，也可得到交流成分，但根据ABB燃烧品质计算公式：FQ＝20*Log(FS)–Background(其中FQ表示火焰质量，FS表示火焰信号，Background表示背景值B)可知，ABB采用了对数放大器对信号进行了转换，并从中减去背景值B。背景值随煤种和工况而变，因此，增大了灭火保护值的整定难度。这也是目前火检存在“偷看”难题的主要原因。比较FORNEY和ABB的火焰特征检测原理，FORNEY是用硬件电路滤除了背景成分，ABB则从信号中减去背景成分。实际上，背景成分还是非常重要的，是可以加以利用的。

发明内容

为了解决上述问题，本发明提出了一种火焰燃烧状态检测装置及检测方法，该装置及方法充分利用火焰的背景成分(直流成分)和闪烁成分(交流成分)，并以其相对大小作为灭火判断依据。

为了实现上述目的，本发明采用如下技术方案：

一种火焰燃烧状态检测装置，包括：光纤、光电转换电路、智能信号处理器以及后台监控主机；

光纤的一端安装于炉膛内，另一端与光电转换电路连接；光电转换电路通过光敏电阻将红外辐射强度转换为电流信号以实现远距离传输；所述光电转换电路与智能信号处理器连接，智能信号处理器通过网络接口与后台监控主机连接；

所述智能信号处理器用于实现信号的采样、A/D转换以及数字滤波，所述智能信号处理器的数量根据炉膛中喷燃器的数量确定。

所述光电转换电路包括：

反相运算放大器U1的正向输入端通过电阻R2接地，反相运算放大器U1的反向输入端经过反相运放输入电阻R1后接入负的基准电压-Vref，红外光敏电阻Rt串接在反相运算放大器U1的反向输入端和输出端之间；反相运算放大器U1的输出端连接电压-电流转换器后输出。

所述智能信号处理器包括主控单元MCU，所述主控单元MCU内部集成A/D转换器、数据存储器RAM、程序存储器和通信接口UART；

光电转换电路输出的电流接入同相输入跟随方式运算放大器U2的正向输入端，所述同相输入跟随方式运算放大器U2的正向输入端经电阻R3后接地，同相输入跟随方式运算放大器U2的反向输入端与输出端连接，同相输入跟随方式运算放大器U2的输出端连接主控单元MCU的A/D转换器；所述主控单元MCU的通信接口UART经过差分发送接收适配器TxRxU与网络接口连接。

一种火焰燃烧状态检测方法，包括以下步骤：

1)光电转换：炉膛的红外光通过光纤照射到光敏电阻，光电转换电路将电阻信号转换为电流信号输出，输出电流的大小与辐射强度呈反比；

2)信号处理：智能信号处理器对来自光电转换电路的电流信号转换为电压信号后依次进行快速采样、A/D转换和数字滤波；

3)频谱分析：对采样信号u(t)进行数字频谱分析，得到其中的直流成分，以及若干频率点的交流成分的大小；

4)指数计算：计算所述交流成分之和，并除以直流成分，得到火焰的闪烁成分与背景成分之比，所述比值称之为燃烧指数；

5)逻辑判断：将所得燃烧指数与设定的门槛值相比较，如果燃烧指数小于门槛值且交流变化幅度下降到门槛值以下，则认为火焰灭火，并发出报警信号，否则，认为火焰燃烧正常。

所述步骤1)的具体方法为：

用精密基准源产生一个负参考电压-Vref，并接入运算放大器的反相端，将光敏电阻Rt作为反馈电阻，放大器的输出电压Vo与输入电压Vref以及光敏电阻Rt呈现如下关系：

$\mathrm{Vo}=\frac{R_t}{R_1}{V}_{\mathrm{ref}}=k\·R_t---\left(0\right)$

由于光敏电阻的阻值随辐射强度增大而减小，因此，输出电压与辐射强度呈反比关系；由于电压-电流变换是线性变换，因此，光电转换电路的电流也与辐射强度呈反比关系。

所述步骤2)的具体方法为：

将光电转换后的信号u(t)表示为一个基频为ω的傅立叶函数：

$u\left(t\right)=U_0{e}^{-\mathrm{\λt}}+\underset{k=1}{\overset{M}{\mathrm{\Σ}}}{U}_k\sin (\mathrm{k\ωt}+{\mathrm{\θ}}_k)---\left(1\right)$

其中，U₀为直流分量；λ为直流分量衰减系数；U_k为第k次谐波的幅值；ω为基波电流或电压的角频率；θk为第k次谐波的初始相角，k＝1,2,…,M；

将(1)式中的U₀e^-λ^t和sin(kωt+θ_k)分别按泰勒级数和三角函数展开，得到：

U₀e^-λt≈U₀-U₀λt                          (2)

U_ksin(kωt+θ_k)＝sin(kωt)U_kcos(θ_k)+cos(kωt)U_ksin(θ_k)  (3)

上式中，U_kcos(θ_k)和U_ksin(θ_k)分别是需要测量的第k次谐波相量的实部和虚部；

根据最小二乘方判据，每次采样后将得到如下一个方程：

$U_0-U_0\mathrm{\λ}{t}_i+\underset{k=1}{\overset{M}{\mathrm{\Σ}}}\left\{\sin \left(\mathrm{k\ω}{t}_i\right)\right[U_k\cos \left({\mathrm{\θ}}_k\right)]+\cos \left(\mathrm{k\ω}{t}_i\right)[U_k\sin \left({\mathrm{\θ}}_k\right)\left]\right\}=u\left(t_i\right)---\left(4\right)$

t_i为第i次采样时刻，经过连续N次采样后，将得到N个方程；

将U₀、U₀λ，以及所有谐波的实部和虚部作为未知量，且用矩阵表示，则N次采样结果可用如下采样矩阵方程表示：

$\left(\begin{array}{ccccccc}1& t_1& \sin \left(\mathrm{\ω}{t}_1\right)& \cos \left(\mathrm{\ω}{t}_1\right)& ...& \sin \left(\mathrm{M\ω}{t}_1\right)& \cos \left(M{\mathrm{\ωt}}_1\right)\\ 1& t_2& \sin \left(\mathrm{\ω}{t}_2\right)& \cos \left(\mathrm{\ω}{t}_2\right)& ...& \sin \left(\mathrm{M\ω}{t}_2\right)& \cos \left(\mathrm{M\ω}{t}_2\right)\\ & .& .& .& ...& .& .\\ 1& t_N& \sin \left(\mathrm{\ω}{t}_N\right)& \cos \left(\mathrm{\ω}{t}_N\right)& ...& \sin \left(\mathrm{M\ω}{t}_N\right)& \cos \left(\mathrm{M\ω}{t}_N\right)\end{array}\right)\left(\begin{array}{c}{U}_0\\ -\mathrm{\λ}{U}_0\\ U_1\cos \left({\mathrm{\θ}}_1\right)\\ U_1\sin \left({\mathrm{\θ}}_1\right)\\ .\\ U_M\cos \left({\mathrm{\θ}}_M\right)\\ U_M\sin \left({\mathrm{\θ}}_M\right)\end{array}\right)=\left(\begin{array}{c}u\left(1\right)\\ u\left(2\right)\\ .\\ u\left(N\right)\end{array}\right)$

采用等间隔采样，则且

用A表示N行，2(M+1)列的系数矩阵，用X表示单列2(M+1)行的待测变量矩阵，用U表示单列采样矩阵，则上述采样矩阵方程可表示为：

A·X＝U               (5)

如果系数矩阵A存在逆矩阵A^-1，则：

X＝A^-1·U               (6)

其中，逆矩阵A^-1的维数为2(M+1)行，N列；

A^-1即是最小二乘滤波器，如果用a_ij表示其第i行第j列的元素，则，i＝1～2(M+1)，j＝1～N；

根据矩阵X的定义和方程(6)，则直流成分为：

$U_0=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{1j}u\left(j\right)---\left(7\right)$

第k次谐波(复数)的实部为：

$U_k\cos \left({\mathrm{\θ}}_k\right)=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{(2k+1)j}u\left(j\right)---\left(8\right)$

第k次谐波(复数)的虚部为：

$U_k\sin \left({\mathrm{\θ}}_k\right)=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{(2k+2)j}u\left(j\right)---\left(9\right)$

复数第k次谐波的峰值为：

$U_k=\sqrt{{\left[U_k\cos \left({\mathrm{\θ}}_k\right)\right]}^2+{\left[U_k\sin \left({\mathrm{\θ}}_k\right)\right]}^2}---\left(10\right)$

所述步骤3)中的燃烧指数计算公式为：

$P=\frac{\mathrm{\Σ}{U}_k}{U_0}---\left(11\right)$

其中，U_k为频率点交流成分之和，k为频率点的个数，U₀为背景的直流成分强度。

所述步骤3)中若干频率点的交流成分分别为低频9Hz、中频36Hz和高频81Hz三个频率点的交流成分。

本发明的有益效果：

1.采用本发明提出的火焰检测方法解决了目前火检装置领域长期存在的技术难题，为火焰的准确检测提供了理论依据；

2.采用本发明提出的火焰检测方法开发的产品在技术指标上已超越了目前ABB和FORNEY等外国公司的产品，大大提高了民族品牌产品的竞争力。

3.火检装置作为火力发电中锅炉保护的主要设备之一，且必须定期更换，因此，具有巨大的经济效益和社会效益。

附图说明

图1为本发明光电转换电路示意图。

图2为本发明智能信号处理器原理结构示意图。

图3为本发明实施例中火焰燃烧状态检测装置IFS10结构示意图。

具体实施方式

下面结合附图与实施例对本发明作进一步说明。

一种火焰燃烧状态检测装置，具体包括：光纤的一端安装于炉膛内，另一端安装于光电转换电路内，其作用为将观测点的红外光引出。光电转换电路通过硫化铅(PbS)光敏电阻将红外辐射强度转换为电流信号以实现远距离传输。智能信号处理器实现信号的采样、A/D转换，以及数字滤波和逻辑判断，并通过网络接口与后台主机相联。一个智能信号处理器检测一个喷燃器的燃烧状态，一个炉膛有多少喷燃器就需要多少个智能信号处理器。后台主机可将这些智能信号处理器联成一个网络，并集中监控和记录历史数据。

光电转换电路如图1所示，包括：

反相运算放大器U1的正向输入端通过电阻R2接地，反相运算放大器U1的反向输入端经过反相运放输入电阻R1后接入负的基准电压-Vref，红外光敏电阻Rt串接在反相运算放大器U1的反向输入端和输出端之间；反相运算放大器U1的输出端连接电压-电流转换器后输出。

智能信号处理器如图2所示，包括主控单元MCU，所述主控单元MCU内部集成A/D转换器、数据存储器RAM、程序存储器和通信接口UART；

光电转换电路输出的电流接入同相输入跟随方式运算放大器U2的正向输入端，所述同相输入跟随方式运算放大器U2的正向输入端经电阻R3后接地，同相输入跟随方式运算放大器U2的反向输入端与输出端连接，同相输入跟随方式运算放大器U2的输出端连接主控单元MCU的A/D转换器；所述主控单元MCU的通信接口UART经过差分发送接收适配器TxRxU与网络接口连接。

目前，采用该检测方法的我公司第一代智能火焰检测分析装置IFS10已在山东某电厂000MW机组投用。该机组共有16个煤粉燃烧器12个油燃烧器，共28个火焰检测装置。在IFS10中，每个火检装置包括一套就地安装的光电转换电路，一个智能信号处理器。智能信号处理器具有网络接口，可连成一个RS485网络，由一台监控主机监视和操作。每个智能模块内各有一个高性能数字处理器(DSP)。因此，一个300MW机组的28个智能模块通过RS-485总线与监控主机构成一套完整装置。系统结构如图3所示。光电转换器安装于就地；信号处理器安装于电子间；监控主机安装于工程师站或运行操作员站。因三者相互距离都较远，所以光电转换器与信号处理器之间的信号采用电流传输，信号处理器与主机采用RS485网络。

一种火焰燃烧状态检测方法，包括以下步骤：

1)光电转换。通过光纤将炉膛的红外光照射到硫化铅光敏电阻，光电转换器将电阻信号转换为电压信号，之后再将电压信号变换为电流信号以便远距离传输，输出电流的大小与辐射强度呈反比；

用精密基准源产生一个负参考电压-Vref，并接入运算放大器的反相端，将硫化铅光敏电阻作为反馈电阻，放大器的输出电压Vo与输入电压Vref以及光敏电阻呈现如下关系：

$\mathrm{Vo}=\frac{R_t}{R_1}{V}_{\mathrm{ref}}=k\·R_t---\left(0\right)$

由于光敏电阻的阻值随辐射强度增大而减小，因此，输出电压便与辐射强度呈反比关系。由于电压-电流变换是线性变换，因此，光电转换器的电流也与辐射强度呈反比关系。

2)信号处理。信号处理器对来自光电转换器的电流信号转换为电压后进行快速采样、A/D转换和数字滤波；

3)频谱分析。信号处理器对采样信号进行数字频谱分析，得到其中的直流成分，以及低频9Hz、中频36Hz、高频81Hz三个频率点的交流成分的大小；

将光电转换后的信号u(t)表示为一个基频为ω的傅立叶函数：

$u\left(t\right)=U_0{e}^{-\mathrm{\λt}}+\underset{k=1}{\overset{M}{\mathrm{\Σ}}}{U}_k\sin (\mathrm{k\ωt}+{\mathrm{\θ}}_k)---\left(1\right)$

其中，U₀为直流分量；λ为直流分量衰减系数；U_k为第k次谐波的幅值；ω为基波电流或电压的角频率；θk为第k次谐波的初始相角，k＝1,2,…,M；

将(1)式中的U₀e^-λ^t和sin(kωt+θ_k)分别按泰勒级数和三角函数展开，得到：

U₀e^-λt≈U₀-U₀λt            (2)

U_ksin(kωt+θ_k)＝sin(kωt)U_kcos(θ_k)+cos(kωt)U_ksin(θ_k)  (3)

上式中，U_kcos(θ_k)和U_ksin(θ_k)分别是需要测量的第k次谐波相量的实部和虚部；

根据最小二乘方判据，每次采样后将得到如下一个方程：

$U_0-U_0\mathrm{\λ}{t}_i+\underset{k=1}{\overset{M}{\mathrm{\Σ}}}\left\{\sin \left(\mathrm{k\ω}{t}_i\right)\right[U_k\cos \left({\mathrm{\θ}}_k\right)]+\cos \left(\mathrm{k\ω}{t}_i\right)[U_k\sin \left({\mathrm{\θ}}_k\right)\left]\right\}=u\left(t_i\right)---\left(4\right)$

ti为第i次采样时刻，经过连续N次采样后，将得到N个方程；

将U₀、U₀λ，以及所有谐波的实部和虚部作为未知量，且用矩阵表示，则N次采样结果可用如下采样矩阵方程表示：

$\left(\begin{array}{ccccccc}1& t_1& \sin \left(\mathrm{\ω}{t}_1\right)& \cos \left(\mathrm{\ω}{t}_1\right)& ...& \sin \left(\mathrm{M\ω}{t}_1\right)& \cos \left(M{\mathrm{\ωt}}_1\right)\\ 1& t_2& \sin \left(\mathrm{\ω}{t}_2\right)& \cos \left(\mathrm{\ω}{t}_2\right)& ...& \sin \left(\mathrm{M\ω}{t}_2\right)& \cos \left(\mathrm{M\ω}{t}_2\right)\\ & .& .& .& ...& .& .\\ 1& t_N& \sin \left(\mathrm{\ω}{t}_N\right)& \cos \left(\mathrm{\ω}{t}_N\right)& ...& \sin \left(\mathrm{M\ω}{t}_N\right)& \cos \left(\mathrm{M\ω}{t}_N\right)\end{array}\right)\left(\begin{array}{c}{U}_0\\ -\mathrm{\λ}{U}_0\\ U_1\cos \left({\mathrm{\θ}}_1\right)\\ U_1\sin \left({\mathrm{\θ}}_1\right)\\ .\\ U_M\cos \left({\mathrm{\θ}}_M\right)\\ U_M\sin \left({\mathrm{\θ}}_M\right)\end{array}\right)=\left(\begin{array}{c}u\left(1\right)\\ u\left(2\right)\\ .\\ u\left(N\right)\end{array}\right)$

采用等间隔采样，则且

用A表示N行，2(M+1)列的系数矩阵，用X表示单列2(M+1)行的待测变量矩阵，用U表示单列采样矩阵，则上述采样矩阵方程可表示为：

A·X＝U              (5)

如果系数矩阵A存在逆矩阵A^-1，则：

X＝A^-1·U               (6)

其中，逆矩阵A-1的维数为2(M+1)行，N列；

A^-1即是最小二乘滤波器，如果用a_ij表示其第i行第j列的元素，则，i＝1～2(M+1)，j＝1～N；

根据矩阵X的定义和方程(6)，则直流成分为：

$U_0=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{1j}u\left(j\right)---\left(7\right)$

复数第k次谐波的实部为：

$U_k\cos \left({\mathrm{\θ}}_k\right)=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{(2k+1)j}u\left(j\right)---\left(8\right)$

复数第k次谐波的虚部为：

$U_k\sin \left({\mathrm{\θ}}_k\right)=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{(2k+2)j}u\left(j\right)---\left(9\right)$

复数第k次谐波的峰峰值为：

$U_k=\sqrt{{\left[U_k\cos \left({\mathrm{\θ}}_k\right)\right]}^2+{\left[U_k\sin \left({\mathrm{\θ}}_k\right)\right]}^2}---\left(10\right)$

4)指数计算。信号处理器计算三个个交流成分之和，并除以直流成分，得到火焰的闪烁成分与背景成分之比，所述比值称之为燃烧指数；

燃烧指数计算公式为：

$P=\frac{\underset{k=\mathrm{9,36,81}}{\mathrm{\Σ}}{U}_k}{U_0}---\left(11\right)$

其中，分子中的交流和为9Hz、36Hz、81Hz谐波的大小之和，U₀为背景的直流成分强度。

5)逻辑判断。微处理器将所得燃烧指数与设定的门槛值相比较，如果燃烧指数小于门槛值且交流变化幅度也下降到门槛值以下，则认为火焰灭火，并发出报警信号，否则，认为火焰燃烧正常。门槛值由实验获得。

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (8)

1.一种火焰燃烧状态检测装置，其特征是，包括：光纤、光电转换电路、智能信号处理器以及后台监控主机； 

光纤的一端安装于炉膛内，另一端与光电转换电路连接；光电转换电路通过光敏电阻将红外辐射强度转换为电流信号以实现远距离传输；所述光电转换电路与智能信号处理器连接，智能信号处理器通过网络接口与后台监控主机连接； 

所述智能信号处理器用于实现信号的采样、A/D转换以及数字滤波，所述智能信号处理器的数量根据炉膛中喷燃器的数量确定。 

2.如权利要求1所述的一种火焰燃烧状态检测装置，其特征是，所述光电转换电路包括： 

反相运算放大器U1的正向输入端通过电阻R2接地，反相运算放大器U1的反向输入端经过反相运放输入电阻R1后接入负的基准电压-Vref，红外光敏电阻Rt串接在反相运算放大器U1的反向输入端和输出端之间；反相运算放大器U1的输出端连接电压-电流转换器后输出。 

3.如权利要求1所述的一种火焰燃烧状态检测装置，其特征是，所述智能信号处理器包括主控单元MCU，所述主控单元MCU内部集成A/D转换器、数据存储器RAM、程序存储器和通信接口UART； 

光电转换电路输出的电流接入同相输入跟随方式运算放大器U2的正向输入端，所述同相输入跟随方式运算放大器U2的正向输入端经电阻R3后接地，同相输入跟随方式运算放大器U2的反向输入端与输出端连接，同相输入跟随方式运算放大器U2的输出端连接主控单元MCU的A/D转换器；所述主控单元MCU的通信接口UART经过差分发送接收适配器TxRxU与网络接口连接。 

4.一种火焰燃烧状态检测方法，其特征是，包括以下步骤： 

1)光电转换：炉膛的红外光通过光纤照射到光敏电阻，光电转换电路将电阻信号转换为电流信号输出，输出电流的大小与辐射强度呈反比； 

2)信号处理：智能信号处理器对来自光电转换电路的电流信号转换为电压信号后依次进行快速采样、A/D转换和数字滤波； 

3)频谱分析：对采样信号u(t)进行数字频谱分析，得到其中的直流成分，以及若干频率点的交流成分的大小； 

4)指数计算：计算所述交流成分之和，并除以直流成分，得到火焰的闪烁成分与背景成分之比，所述比值称之为燃烧指数； 

5)逻辑判断：将所得燃烧指数与设定的门槛值相比较，如果燃烧指数小于门槛值且交流 变化幅度下降到门槛值以下，则认为火焰灭火，并发出报警信号，否则，认为火焰燃烧正常。 

5.如权利要求4所述的一种火焰燃烧状态检测方法，其特征是，所述步骤1)的具体方法为： 

用精密基准源产生一个负参考电压-Vref，并接入运算放大器的反相端，将光敏电阻Rt作为反馈电阻，放大器的输出电压Vo与输入电压Vref以及光敏电阻Rt呈现如下关系： 

由于光敏电阻的阻值随辐射强度增大而减小，因此，输出电压与辐射强度呈反比关系；由于电压-电流变换是线性变换，因此，光电转换电路的电流也与辐射强度呈反比关系。 

6.如权利要求4所述的一种火焰燃烧状态检测方法，其特征是，所述步骤2)的具体方法为： 

将光电转换后的信号u(t)表示为一个基频为ω的傅立叶函数： 

其中，U₀为直流分量；λ为直流分量衰减系数；U_k为第k次谐波的幅值；ω为基波电流或电压的角频率；θk为第k次谐波的初始相角，k＝1,2,…,M； 

将(1)式中的U₀e^-λ^t和sin(kωt+θ_k)分别按泰勒级数和三角函数展开，得到： 

U₀e^-λt≈U₀-U₀λt                (2) 

U_ksin(kωt+θ_k)＝sin(kωt)U_kcos(θ_k)+cos(kωt)U_ksin(θ_k)     (3) 

上式中，U_kcos(θ_k)和U_ksin(θ_k)分别是需要测量的第k次谐波相量的实部和虚部； 

根据最小二乘方判据，每次采样后将得到如下一个方程： 

t_i为第i次采样时刻，经过连续N次采样后，将得到N个方程； 

将U₀、U₀λ，以及所有谐波的实部和虚部作为未知量，且用矩阵表示，则N次采样结果可用如下采样矩阵方程表示： 

采用等间隔采样，则且

用A表示N行，2(M+1)列的系数矩阵，用X表示单列2(M+1)行的待测变量矩阵，用U表示单列采样矩阵，则上述采样矩阵方程可表示为： 

A·X＝U                     (5) 

如果系数矩阵A存在逆矩阵A^-1，则： 

X＝A^-1·U                     (6) 

其中，逆矩阵A^-1的维数为2(M+1)行，N列； 

A^-1即是最小二乘滤波器，如果用a_ij表示其第i行第j列的元素，则，i＝1～2(M+1)，j＝1～N； 

根据矩阵X的定义和方程(6)，则直流成分为： 

复数第k次谐波的实部为： 

复数第k次谐波的虚部为： 

复数第k次谐波的峰值为： 

。

7.如权利要求4所述的一种火焰燃烧状态检测方法，其特征是，所述步骤3)中的燃烧 指数计算公式为： 

其中，U_k为频率点交流成分之和，k为频率点的个数，U₀为背景的直流成分强度。 

8.如权利要求4所述的一种火焰燃烧状态检测方法，其特征是，所述步骤3)中若干频率点的交流成分分别为低频9Hz、中频36Hz和高频81Hz三个频率点的交流成分。