李珊，閆勇，吳佳麗，錢相臣

(1.華北電力大學(xué)控制與計(jì)算機(jī)工程學(xué)院，北京，102206；2.肯特大學(xué)工程與數(shù)字藝術(shù)學(xué)院，坎特伯雷，肯特CT2 7NT，英國(guó))

能源的低碳化和可再生化是全球能源結(jié)構(gòu)轉(zhuǎn)型的發(fā)展趨勢(shì)。由于沼氣、生物質(zhì)等資源在碳排放和可再生性方面的優(yōu)勢(shì)以及相關(guān)技術(shù)突破和升級(jí)，它們替代高污染的煤炭等傳統(tǒng)化石能源已經(jīng)成為不可逆轉(zhuǎn)的趨勢(shì)[1]。甲烷是沼氣的主要成分和生物質(zhì)裂解的重要產(chǎn)物，相對(duì)于煤炭和石油，直接燃燒甲烷可以大大降低煙氣和焦油等雜質(zhì)的排放。然而，清潔能源的高效燃燒利用在技術(shù)和應(yīng)用層面仍存在諸多問(wèn)題亟待解決，需要研究甲烷及生物質(zhì)的燃燒過(guò)程，以優(yōu)化燃燒系統(tǒng)設(shè)計(jì)，達(dá)到提高燃料燃燒效率和降低污染物排放的目的?；鹧媸侨紵^(guò)程最直觀的表現(xiàn)形式，針對(duì)火焰的不同特性可以使用相應(yīng)類型的火焰檢測(cè)方法來(lái)進(jìn)行測(cè)量。目前，工業(yè)上廣泛使用的紅外、紫外和可見(jiàn)光檢測(cè)裝置利用了火焰的輻射特性[2?3]，主要用于檢測(cè)火焰的有無(wú)，而且容易被惡劣的使用環(huán)境污染，所以無(wú)法長(zhǎng)期、實(shí)時(shí)監(jiān)測(cè)火焰參數(shù)的動(dòng)態(tài)變化。圖像火檢技術(shù)通過(guò)高速相機(jī)等成像設(shè)備實(shí)時(shí)采集火焰圖像并通過(guò)圖像處理得到特征參數(shù)反映燃燒狀況[4?5]。這類方法實(shí)時(shí)性強(qiáng)準(zhǔn)確度高，但是運(yùn)行維護(hù)成本高和可視范圍小等缺點(diǎn)仍阻礙了該方法的大規(guī)模應(yīng)用。此外，利用光譜儀進(jìn)行火焰的燃燒狀態(tài)檢測(cè)主要在實(shí)驗(yàn)室規(guī)模下進(jìn)行[6]，并不適用于工業(yè)條件下的長(zhǎng)期在線監(jiān)測(cè)。燃燒過(guò)程包含一系列復(fù)雜的化學(xué)反應(yīng)，期間會(huì)產(chǎn)生大量帶電粒子，包括陰陽(yáng)離子、自由電子以及碳煙顆粒，因此火焰具有微弱的帶電特性[7]。非侵入式靜電傳感技術(shù)通過(guò)對(duì)火焰的電學(xué)特性進(jìn)行分析從而對(duì)火焰燃燒狀態(tài)進(jìn)行檢測(cè)，具有靈敏度高、環(huán)境適應(yīng)性強(qiáng)、成本低等優(yōu)勢(shì)[8?9]。而獲取的火焰靜電信號(hào)通常較為微弱，其時(shí)域頻域包含大量有用信息，但容易受到環(huán)境背景噪聲的干擾。經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition，EMD)法是HUANG等[10]于1998年提出的一種自適應(yīng)的信號(hào)分析方法。該方法能根據(jù)信號(hào)自身時(shí)間尺度特征(不需要任何先決條件)自適應(yīng)的將復(fù)雜信號(hào)分解為有限個(gè)頻率由高到低順序排列的本征模函數(shù)(intrinsic mode function，IMF)和殘余分量，各IMF 分量包含了原信號(hào)在不同時(shí)間尺度的局部特征信號(hào)[10]。EMD避免了小波變換等方法中基函數(shù)選擇的問(wèn)題，常被應(yīng)用于非平穩(wěn)隨機(jī)信號(hào)的分析，具有很高的信噪比[10]?；鹧嬷徐o電信號(hào)為類隨機(jī)信號(hào)，其具體特征目前還在探索階段，額外引入基函數(shù)可能對(duì)信號(hào)本身造成干擾，因此，EMD 分解是一種較好的適用于靜電信號(hào)的去噪方法。信息熵是信息論中用于描述信息不確定度的重要概念，近年來(lái)被廣泛應(yīng)用于機(jī)械軸承、電力電纜的故障診斷等領(lǐng)域[11]。將靜電信號(hào)作為信息載體同時(shí)結(jié)合EMD 分解與信息熵得到火焰靜電信號(hào)的能量熵，實(shí)現(xiàn)對(duì)火焰靜電信號(hào)頻率和能量分布的深入了解?；鹧娴奶卣魈崛∈茄芯炕鹧嫒紵隣顟B(tài)和穩(wěn)定性的基礎(chǔ)，本文作者利用非侵入式靜電傳感器陣列實(shí)現(xiàn)對(duì)甲烷及生物質(zhì)火焰靜電信號(hào)的檢測(cè)，然后結(jié)合經(jīng)驗(yàn)?zāi)B(tài)分解和信息熵理論提取火焰靜電信號(hào)的能量熵等特征參數(shù)，分析火焰靜電信號(hào)能量熵與燃料燃燒狀態(tài)間的聯(lián)系，提供了一種適用于沼氣和生物質(zhì)燃燒狀態(tài)分析、特征提取和穩(wěn)定性檢測(cè)的新方法。

1 靜電法火焰檢測(cè)原理

1.1 火焰的帶電特性

燃料燃燒過(guò)程中與離子相關(guān)的反應(yīng)包括火焰的熱電離、化學(xué)電離以及電離之后的離子交換和離子重組。在化學(xué)電離中，基本放熱化學(xué)反應(yīng)的能量使產(chǎn)物之一處于電離狀態(tài)，一旦通過(guò)化學(xué)電離產(chǎn)生1個(gè)離子，更多離子的產(chǎn)生可以通過(guò)電荷交換實(shí)現(xiàn)[12?13]，比如烴類火焰由化學(xué)電離產(chǎn)生初始離子CHO+：

火焰中常見(jiàn)的H3O+則可以通過(guò)由化學(xué)電離產(chǎn)生的初始離子由下列過(guò)程生成：

進(jìn)而在含有過(guò)量烴基的富火焰中，烴類離子的形成很容易通過(guò)諸如此類的反應(yīng)進(jìn)行：

由于燃料的不完全燃燒，火焰中存在著大量的碳煙顆粒。SUGDEN 等[14]首次展示了碳煙顆粒荷電的實(shí)驗(yàn)結(jié)果，隨后被證明火焰中每個(gè)顆粒最多荷有1～2個(gè)電子的電量[15]。較早的研究指出火焰中荷正負(fù)電的碳煙顆粒各占其總數(shù)的30%左右，并且隨著火焰高度的變化，碳煙顆粒的荷電情況也有所改變[16]。但是，王宇等[17]則通過(guò)正負(fù)電場(chǎng)不對(duì)稱實(shí)驗(yàn)證明火焰中有大于50%的碳煙顆粒荷正電，進(jìn)一步明確了火焰的帶電特性。

1.2 火焰靜電信號(hào)測(cè)量

目前，靜電檢測(cè)技術(shù)已被成熟應(yīng)用于氣固兩相流參數(shù)檢測(cè)領(lǐng)域，實(shí)現(xiàn)氣力輸送管道中固體顆粒的速度、濃度、質(zhì)量流量等參數(shù)的測(cè)量[8,18?19]。類似的，火焰中存在復(fù)雜的氣、固和等離子多相運(yùn)動(dòng)，燃料燃燒及顆粒碰撞摩擦產(chǎn)生的靜電電荷能夠反映火焰中顆粒的流動(dòng)特性。由于火焰的波動(dòng)帶電顆粒接觸到裸露靜電電極后在電極表面產(chǎn)生傳導(dǎo)電荷，從而產(chǎn)生靜電信號(hào)，因此可以在火焰周圍布置靜電傳感裝置進(jìn)行燃料燃燒狀態(tài)檢測(cè)。

由于火焰中靜電信號(hào)通常較為微弱，并且原始靜電信號(hào)中包含了其工作環(huán)境中的背景噪聲，為提高特征提取的準(zhǔn)確性需要通過(guò)信號(hào)處理電路對(duì)原始信號(hào)進(jìn)行電流/電壓轉(zhuǎn)換、電壓放大和基礎(chǔ)濾波，隨后通過(guò)程序?qū)鬏斨劣?jì)算機(jī)的電壓信號(hào)進(jìn)行降噪和特征提取?；鹧嬷袔щ婎w粒濃度、運(yùn)動(dòng)狀態(tài)的變化都會(huì)導(dǎo)致靜電信號(hào)的變化，因此可以對(duì)獲取的靜電信號(hào)進(jìn)行分析進(jìn)一步了解燃料的燃燒狀態(tài)。

火焰中帶電顆粒的濃度可以通過(guò)靜電傳感器電極所獲靜電信號(hào)的強(qiáng)度體現(xiàn)。均方根(root mean square，RMS)是指一組數(shù)據(jù)的有效值，常用于表征采樣系統(tǒng)信號(hào)的平均水平。氣固兩相流測(cè)量領(lǐng)域使用靜電信號(hào)的RMS 表征管道內(nèi)固體顆粒相對(duì)體積濃度的高低[18?19]，本文同樣利用靜電信號(hào)的RMS 評(píng)估靜電電極對(duì)應(yīng)區(qū)域火焰中帶電顆粒的相對(duì)體積濃度[20]，靜電信號(hào)的RMS計(jì)算公式如下：

式中：N為采樣點(diǎn)數(shù)；xi為第i個(gè)采樣點(diǎn)值；為采樣信號(hào)的均值。

2 基于EMD 分解的靜電信號(hào)能量熵

2.1 基于EMD分解的靜電信號(hào)去噪方法

火焰靜電信號(hào)x(t)的EMD 分解流程如圖1所示，其中SD為h(t)的標(biāo)準(zhǔn)偏差，計(jì)算公式為

SDT為篩分門限值，一般取值0.2～0.3[10]。原始信號(hào)x(t)與所有符合條件的h(t)之和相減即為余項(xiàng)r(t)。

圖1 EMD分解流程圖Fig.1 Flowchart of EMD decomposition

火焰靜電信號(hào)經(jīng)過(guò)EMD 分解后得到的IMF 分量，隨著分解階數(shù)的增加，其頻率逐漸減小。而靜電信號(hào)包含的噪聲通常集中在相對(duì)高頻的區(qū)域，因此可以將高頻IMF 去除后將剩余分量進(jìn)行重構(gòu)達(dá)到去噪目的，但這樣有可能丟失高頻信號(hào)中的有用信號(hào)。為避免直接剔除高頻分量造成的信息損失，本文將對(duì)IMF進(jìn)行如下處理[21]：

1)對(duì)靜電信號(hào)所有的IMF 分量進(jìn)行自相關(guān)計(jì)算，由于隨機(jī)噪聲具有較弱的時(shí)間相關(guān)性，因此其自相關(guān)函數(shù)曲線通常接近沖擊函數(shù)，即在零點(diǎn)處取得最大值后，在零點(diǎn)兩旁迅速衰減至0。具有此類特征的IMF分量，即為噪聲主導(dǎo)的IMF分量；

2)采用固定閾值規(guī)則和軟閾值函數(shù)對(duì)噪聲主導(dǎo)的IMF分量進(jìn)行小波降噪處理；

3)將經(jīng)過(guò)小波閾值降噪處理的IMF 分量和信號(hào)主導(dǎo)的IMF分量與余項(xiàng)r(t)進(jìn)行信號(hào)重構(gòu)，從而得到去噪后的信號(hào)。

這種方法克服了傳統(tǒng)的降噪算法無(wú)法確定噪聲成分所處頻段的缺陷。為驗(yàn)證本方法的降噪效果，使用降噪后甲烷體積流量為0.8 L/min 擴(kuò)散燃燒火焰靜電信號(hào)模擬純火焰靜電信號(hào)，加入信噪比為8 dB 的高斯白噪聲，隨后使用上述方法進(jìn)行降噪，處理后信號(hào)的信噪比為16.986 dB，得到了明顯提升。

2.2 火焰靜電信號(hào)的信息熵和能量熵

SHANNON[22]于1948年提出信息熵的概念，用于量化隨機(jī)變量的不確定程度，作為系統(tǒng)有序化程度的度量，可以理解為某種特定信息的出現(xiàn)概率。假設(shè)某系統(tǒng)X可能處于n種不同的狀態(tài)x1，x2，…，xn，P(xi)代表狀態(tài)xi(i=1，2，…，n)出現(xiàn)的概率，則該系統(tǒng)的信息熵H(X)定義為

由于信息熵的非負(fù)性和可加性，0≤P(xi)≤1 且當(dāng)P(xi)=0時(shí)規(guī)定0lg0=0。

根據(jù)EMD 分解，可以提取出信號(hào)各個(gè)頻段的成分，假設(shè)一個(gè)火焰的靜電信號(hào)可以被分解為n個(gè)IMF 分量，則可以計(jì)算出各分量的能量E1，E2，…，En：

由于EMD 分解具有正交性，因此這n個(gè)IMF分量的能量之和E0與原靜電信號(hào)總能量相等：

根據(jù)信息熵理論，EMD分解的能量熵表示為

其中：pi表示第i個(gè)IMF 分量的能量占全部能量的比例，即pi=Ei/E0。

EMD 分解中能量熵表征的是信號(hào)能量與頻率的變化關(guān)系，能夠反映頻率分布和能量分布的等細(xì)節(jié)信息[11]。根據(jù)能量熵的極值性和最大熵定理，當(dāng)各個(gè)頻段能量分布出現(xiàn)的概率相等時(shí)，能量熵最大。因此，能量熵越小，代表靜電信號(hào)的能量集中在個(gè)別的IMF 分量當(dāng)中，火焰中帶電顆粒波動(dòng)頻率單一，燃燒較為穩(wěn)定。反之，則意味著靜電信號(hào)能量分散在各IMF 分量中，帶電顆粒波動(dòng)頻率范圍寬，顆粒運(yùn)動(dòng)復(fù)雜燃燒穩(wěn)定性差。實(shí)驗(yàn)中對(duì)重構(gòu)后的火焰靜電信號(hào)進(jìn)行EMD 分解并計(jì)算各IMF分量的能量，從而得到各電極信號(hào)的能量熵。

3 實(shí)驗(yàn)結(jié)果分析

3.1 實(shí)驗(yàn)裝置與實(shí)驗(yàn)條件

本研究使用非侵入式靜電傳感技術(shù)通過(guò)測(cè)量甲烷及生物質(zhì)火焰燃燒時(shí)的電荷波動(dòng)實(shí)現(xiàn)燃料燃燒過(guò)程的特征提取。圖2所示為火焰檢測(cè)系統(tǒng)示意圖，該系統(tǒng)主要包括燃燒器、非侵入式靜電傳感器陣列、信號(hào)調(diào)理電路和信號(hào)接收及處理單元。實(shí)驗(yàn)時(shí)，配氣系統(tǒng)按設(shè)定流量向燃燒器提供甲烷和空氣，生物質(zhì)則通過(guò)給料機(jī)以一定的質(zhì)量流量送至管道后由給料空氣吹至燃燒器。實(shí)驗(yàn)所使用的靜電傳感器參數(shù)如圖3所示，該傳感器包含10個(gè)靜電電極，其高度和寬度能夠覆蓋實(shí)驗(yàn)中燃料火焰，測(cè)量火焰各個(gè)區(qū)域帶電顆粒的動(dòng)態(tài)特性，每個(gè)電極配有相同參數(shù)的獨(dú)立信號(hào)處理電路[8?9]。信號(hào)調(diào)理電路包含電流?電壓轉(zhuǎn)換和信號(hào)放大及濾波預(yù)處理的功能。傳感器通過(guò)底部支撐柱固定于圓盤狀支撐架的卡槽內(nèi)，通過(guò)移動(dòng)支撐柱在卡槽中的位置調(diào)整傳感器與火焰間的距離。采集到的靜電信號(hào)由數(shù)據(jù)采集卡(DAQ)接收后傳至計(jì)算機(jī)進(jìn)行處理，提取火焰靜電信號(hào)的特征參數(shù)。

圖2 火焰檢測(cè)系統(tǒng)示意圖Fig.2 Schematic diagram of flame monitoring system

圖3 靜電傳感器結(jié)構(gòu)及幾何參數(shù)Fig.3 Structure and dimensions of electrostatic sensor

表1所示為非侵入式靜電傳感器陣列測(cè)量不同燃料火焰實(shí)驗(yàn)工況設(shè)置，I 組為甲烷擴(kuò)散燃燒工況，共進(jìn)行3組實(shí)驗(yàn)。II組為甲烷與空氣預(yù)混燃燒工況，在3種甲烷體積流量下按照表1所示分別改變空氣流量，共進(jìn)行12 組實(shí)驗(yàn)。III 組為生物質(zhì)預(yù)混燃燒實(shí)驗(yàn)，燃料選擇粒徑范圍為210～230 μm 的玉米芯顆粒，除表1所示的助燃空氣外，給料空氣流量為6 L/min。實(shí)驗(yàn)時(shí)綜合考慮信號(hào)的幅值，將靜電傳感器陣列距離燃燒器中心水平距離固定為25 mm，靜電信號(hào)采樣頻率為250 Hz，每組實(shí)驗(yàn)采樣時(shí)間均大于20 s。

表1 甲烷及生物質(zhì)燃燒條件Table 1 Test Matrix of methane and biomass combustion

3.2 火焰靜電信號(hào)的EMD去噪

本研究采用基于EMD 的降噪方法對(duì)純甲烷擴(kuò)散燃燒、甲烷與空氣預(yù)混燃燒和甲烷助燃的生物質(zhì)燃燒火焰典型工況下的靜電信號(hào)進(jìn)行去噪處理。對(duì)甲烷空氣體積流量比為1.0:3.0(燃空當(dāng)量比為3.18)的預(yù)混火焰靜電信號(hào)進(jìn)行EMD分解，共獲得7 個(gè)IMF 分量及余項(xiàng)，圖4(a)所示為分量IMF 1～I(xiàn)MF4，圖4(b)所示為這4 個(gè)分量自相關(guān)運(yùn)算結(jié)果，IMF5～I(xiàn)MF7 直接參與重構(gòu)故不再列出。根據(jù)前文所述原則IMF1 和IMF2 為噪聲主導(dǎo)的分量，其余分量不具備零點(diǎn)衰減特性，可見(jiàn)噪聲主要分布在信號(hào)的高頻段，為最大程度保留有用信號(hào)，對(duì)此部分分量進(jìn)行小波閾值處理后，與剩余的IMF分量一同進(jìn)行重構(gòu)完成降噪，結(jié)果如圖5(b)所示。由圖5可以看出：甲烷擴(kuò)散火焰靜電信號(hào)幅值比預(yù)混火焰和生物質(zhì)火焰信號(hào)幅值小。此外，由于給料機(jī)給料不均勻，玉米芯燃燒顆粒成簇的噴出燃燒口，并且火焰中時(shí)常噴濺出火星，導(dǎo)致傳導(dǎo)電荷在某個(gè)瞬間突然增大，因此，甲烷助燃的玉米芯火焰靜電信號(hào)具有規(guī)律的類尖峰信號(hào)。

3.3 火焰靜電信號(hào)能量分析

圖4 靜電信號(hào)IMF分量及自相關(guān)函數(shù)Fig.4 IMF components and autocorrelation function of senior signals

圖5 不同工況下火焰靜電信號(hào)去噪前后對(duì)比Fig.5 Comparison of denoising effect of flame signals under different combustion conditions

圖6 甲烷在不同流量條件下火焰靜電信號(hào)能量熵分布Fig.6 Energy entropy distribution of the sensor signals under different methane flow rate conditions

由于實(shí)驗(yàn)中火焰高度最多到達(dá)第10 個(gè)電極，因此將10 個(gè)電極計(jì)算各自能量熵后取平均值作為各工況下整體燃料燃燒靜電信號(hào)能量熵，實(shí)驗(yàn)結(jié)果如圖6所示，其中橙色圓圈標(biāo)注的為甲烷在不同體積流量下擴(kuò)散燃燒的結(jié)果，其余為甲烷體積流量分別為0.6，0.8 和1.0 L/min 時(shí)改變空氣量從而形成的預(yù)混火焰的計(jì)算結(jié)果。此外，甲烷助燃的玉米芯火焰靜電信號(hào)能量熵為1.518。由圖6可以看出：純甲烷擴(kuò)散火焰靜電信號(hào)的能量熵普遍低于相同甲烷流量下預(yù)混火焰靜電信號(hào)的能量熵，并且在甲烷流量相同的情況下，隨著空氣量增大，相應(yīng)靜電信號(hào)的能量熵呈增大趨勢(shì)。與氣態(tài)燃料相比，玉米芯火焰靜電信號(hào)的能量熵明顯較大。根據(jù)能量熵的性質(zhì)，熵值越小能量分布越集中，反之則能量分布較為分散。如圖7所示，對(duì)比火焰中部同一電極在3種典型燃燒工況下靜電信號(hào)各IMF 分量間的能量分布，擴(kuò)散火焰能量集中在較低頻段的IMF 分量中，與其能量熵最小的結(jié)論相符。預(yù)混火焰和生物質(zhì)火焰的能量則分散分布在高頻IMF 分量中，生物質(zhì)火焰能量分布更為分散，各IMF能量占比均在0.5以下，對(duì)應(yīng)相對(duì)較大的能量熵。結(jié)合3種燃燒工況下火焰靜電信號(hào)IMF能量分布圖和能量熵計(jì)算結(jié)果，由于顆粒微觀運(yùn)動(dòng)復(fù)雜度和顆粒運(yùn)動(dòng)的能量間存在相關(guān)關(guān)系[23]，因此靜電信號(hào)能量熵能夠反映靜電傳感器電極敏感區(qū)域內(nèi)顆粒運(yùn)動(dòng)的復(fù)雜程度，熵值越大帶電顆粒運(yùn)動(dòng)復(fù)雜度越高。

圖8所示為各工況下火焰靜電信號(hào)RMS分布。擴(kuò)散燃燒沒(méi)有額外助燃空氣參與，隨著甲烷在環(huán)境中的擴(kuò)散，逐漸與氧氣混合，由于燃料的不完全燃燒而產(chǎn)生的碳煙顆粒致使火焰帶有微弱電量。擴(kuò)散燃燒過(guò)程緩慢，火焰較為無(wú)力容易發(fā)生幾何波動(dòng)[5]，帶電顆粒隨火焰的波動(dòng)而運(yùn)動(dòng)，從而造成靜電信號(hào)能量在IMF 低頻段的集中。擴(kuò)散燃燒下帶電顆粒的低濃度和氣流的緩慢波動(dòng)造成較小的靜電信號(hào)能量熵，帶電顆粒運(yùn)動(dòng)相對(duì)簡(jiǎn)單規(guī)律，燃燒穩(wěn)定。預(yù)混燃燒條件下，由于甲烷和空氣充分混合后再點(diǎn)燃，其火焰面(燃燒產(chǎn)物與預(yù)混氣體的分界面)附近的燃燒化學(xué)反應(yīng)復(fù)雜而劇烈。同時(shí)伴隨著光能和熱能等能量的釋放，產(chǎn)生大量的陰陽(yáng)離子、自由電子以及部分燃料的不完全燃燒產(chǎn)生的碳煙顆粒，因此帶電顆粒濃度大于擴(kuò)散火焰。這個(gè)過(guò)程中火焰中間自由基間的能量轉(zhuǎn)換和反應(yīng)物能量釋放速率的轉(zhuǎn)變是造成能量集中在IMF 高頻段的主要原因。此外，由于預(yù)混火焰中的帶電顆粒濃度較大，而且甲烷氣體和空氣流速之間的差異引起火焰面附近強(qiáng)烈的熱量和質(zhì)量交換，所以帶電顆粒在較為復(fù)雜的氣流條件下發(fā)生摩擦碰撞的概率增大使得帶電顆粒運(yùn)動(dòng)復(fù)雜度增大，這就降低了燃燒的穩(wěn)定性，嚴(yán)重時(shí)過(guò)大的空氣流量會(huì)引發(fā)火焰劇烈抖動(dòng)甚至脫火。

圖7 甲烷擴(kuò)散燃燒、預(yù)混燃燒及生物質(zhì)顆粒燃燒火焰靜電信號(hào)IMF能量分布Fig.7 IMF energy distribution of sensor signals of methane and methane-biomass co-firing flames

甲烷助燃的生物質(zhì)火焰中同時(shí)存在固體和氣體燃料，顆粒在高溫氣體的裹挾下進(jìn)行揮發(fā)分和固定碳的逐級(jí)燃燒。由于生物質(zhì)固體顆粒形狀不規(guī)則，燃燒過(guò)程中顆粒間較為分散，難以完全氧化，并且固定碳完全燃燒需要較長(zhǎng)的燃燒時(shí)間和較高的反應(yīng)溫度，因此生物質(zhì)火焰中長(zhǎng)時(shí)間存在未完全燃燒的固體碳顆粒。除燃燒產(chǎn)生的帶電顆粒外，生物質(zhì)顆粒和未完全燃燒的碳顆粒隨氣流發(fā)生不規(guī)則運(yùn)動(dòng)，從而發(fā)生碰撞摩擦而荷電[16]，造成較大的帶電顆粒濃度。期間火焰中進(jìn)行的氣體間以及氣固間熱質(zhì)交換造成能量在IMF 高頻段的集中分布。此外，生物質(zhì)燃燒過(guò)程中給料空氣與燃料間流速差異更大氣流條件更為復(fù)雜，造成相對(duì)較大的能量熵值，整體火焰中的帶電顆粒反映出更大的運(yùn)動(dòng)復(fù)雜程度，對(duì)應(yīng)較為復(fù)雜且不穩(wěn)定的燃燒狀況。

圖8 甲烷擴(kuò)散燃燒、預(yù)混燃燒及生物質(zhì)顆粒燃燒火焰靜電信號(hào)RMS分布Fig.8 RMS distribution of sensor signals of methane and methane-biomass co-firing flames

圖9所示為3種典型燃燒工況下不同電極的能量熵，電極自下向上進(jìn)行編號(hào)分別對(duì)應(yīng)火焰自根部至頂部的燃燒區(qū)域。甲烷擴(kuò)散火焰和預(yù)混火焰能量熵呈現(xiàn)根部、頂部小，中上部大的分布，預(yù)混火焰分布得更為規(guī)律。這是因?yàn)闅鈶B(tài)火焰不同區(qū)域呈現(xiàn)不同的燃燒狀態(tài)[16]，火焰自底部逐漸與空氣混合后在火焰中部進(jìn)行劇烈、復(fù)雜的化學(xué)反應(yīng)導(dǎo)致這一區(qū)域氣流紊亂。大量的陰陽(yáng)離子和碳煙顆粒在氣流的影響下劇烈運(yùn)動(dòng)，因而火焰中部帶電顆粒的運(yùn)動(dòng)復(fù)雜度較大，相比較下，火焰根部燃燒則更加穩(wěn)定。擴(kuò)散火焰整體氣流平穩(wěn)，帶電顆粒容易聚集在化學(xué)反應(yīng)劇烈的火焰中部，缺少空氣氣流的影響而向兩邊擴(kuò)散較慢。中部較大的帶電顆粒濃度使得顆粒間碰撞摩擦概率增大，火焰根部雖有電荷聚集但較難發(fā)生明顯波動(dòng)呈穩(wěn)定燃燒狀態(tài)，因此，擴(kuò)散火焰中部能量熵明顯大于火焰根部和頂部能量熵。

圖9 甲烷擴(kuò)散燃燒、預(yù)混燃燒及生物質(zhì)顆粒燃燒火焰不同區(qū)域靜電信號(hào)能量熵分布Fig.9 Energy entropy distribution in different regions of methane and methane-biomass co-firing flames

預(yù)混火焰中部火焰面附近燃燒劇烈，火焰面外側(cè)燃燒產(chǎn)物聚集，使得帶電顆粒在火焰中上部濃度較大，帶電顆粒在空氣氣流和重力的共同影響下造成火焰中部較大的能量熵和復(fù)雜的燃燒狀況。由于復(fù)雜的氣固運(yùn)動(dòng)和氣流影響，生物質(zhì)火焰靜電信號(hào)能量熵在整個(gè)火焰高度范圍內(nèi)差別不大。生物質(zhì)燃燒火焰根部未燃燒顆粒在給料空氣的氣流中劇烈運(yùn)動(dòng)，碰撞摩擦而帶電使得靜電信號(hào)能量熵增大，火焰中部與氣態(tài)火焰相似同樣為劇烈燃燒區(qū)，燃燒產(chǎn)生的大量帶電顆粒與生物質(zhì)顆粒一同在復(fù)雜的氣流下劇烈運(yùn)動(dòng)。由于空氣流速較大，火焰根部和中部的部分帶電粒子運(yùn)動(dòng)到火焰頂部，火焰頂部發(fā)生劇烈抖動(dòng)。此外，生物質(zhì)火焰尖端還會(huì)有未完全燃燒的碳顆粒造成的火星飛濺與靜電電極進(jìn)行接觸，因此，生物質(zhì)燃燒火焰內(nèi)帶電粒子整體都呈現(xiàn)較大的運(yùn)動(dòng)復(fù)雜度

，各部位能量熵差別不大，燃燒穩(wěn)定性較差。

4 結(jié)論

1)分析了非侵入式靜電傳感器陣列測(cè)得的純甲烷擴(kuò)散火焰、甲烷空氣預(yù)混火焰和甲烷助燃的生物質(zhì)火焰靜電信號(hào)，并針對(duì)火焰靜電信號(hào)的類隨機(jī)特征采用基于EMD 分解的降噪方法進(jìn)行降噪處理去除了高頻隨機(jī)噪聲，達(dá)到了良好的去噪效果。

2)擴(kuò)散火焰靜電信號(hào)的平均能量熵隨燃料流量的增大而增大，且相同甲烷流量下空氣流量越大預(yù)混火焰能量熵越高，甲烷助燃的生物質(zhì)火焰能量熵較甲烷氣態(tài)燃燒火焰大。擴(kuò)散火焰和預(yù)混火焰靜電信號(hào)能量熵均呈現(xiàn)火焰中部大，根部和頂部小的趨勢(shì)，生物質(zhì)火焰各部位能量熵差別較小。

3)火焰靜電信號(hào)能量熵越大帶電顆粒運(yùn)動(dòng)復(fù)雜度越高，燃料燃燒的穩(wěn)定性越差。從能量熵出發(fā)聯(lián)系火焰帶電顆粒運(yùn)動(dòng)情況與燃燒狀態(tài)，為火焰燃燒狀態(tài)分析、特征提取和穩(wěn)定性檢測(cè)提供一種新途徑，為可再生低碳排放燃料的高效利用提供新思路。