卓增泗 胡立生

(上海交通大學電子信息與電氣工程學院)

基于小波變換的燃燒不穩(wěn)定性檢測方法研究

卓增泗 胡立生

(上海交通大學電子信息與電氣工程學院)

利用小波變換多尺度分析分解和重構(gòu)燃氣輪機燃燒室的高頻、中頻和低頻壓力信號，并在時域上定位燃燒不穩(wěn)定性發(fā)生的時刻。對于具有不穩(wěn)定成分的壓力信號，該方法能夠正確地檢測出不穩(wěn)定性，結(jié)果與電廠實際運行情況一致。

燃燒不穩(wěn)定性 小波變換 壓力振動 多尺度分析

我國早期引入的E級燃氣輪機在擴散燃燒方式下，采用輕油所排放的NOx的質(zhì)量濃度值為410mg/m3(標準狀態(tài), 15% O2，下同)，而采用天然氣排放的NOx的質(zhì)量濃度值則是160～220mg/m3[1]。我國在《火電廠大氣污染物排放標準》(GB 13223-2011)中規(guī)定[2]，采用燃油的燃氣輪機發(fā)電廠排放的NOx質(zhì)量濃度值不高于120mg/m3，而采用天然氣的燃氣輪機發(fā)電廠排放的NOx質(zhì)量濃度值不高于50mg/m3[2]，所以早期引入的E級燃氣輪機的NOx排放濃度不能滿足GB 13223-2011中的規(guī)定。

為了降低燃氣輪機排放污染物的濃度，我國近年來引入的F級燃氣輪機采用的是干低式NOx(Dry Low NOx， DLN)燃燒方式[3,4]，包含有預(yù)混式燃燒器[5]，燃料在進入燃燒區(qū)域之前先和過量的空氣混合均勻，再進入燃燒區(qū)域燃燒，以此降低燃燒溫度，減少燃燒過程中產(chǎn)生熱力型NOx，目前最低可以達到9μL/L，遠低于國家標準。

但是，對于使用干低式NOx技術(shù)的燃氣輪機，由于燃燒溫度相對較低，而且燃料空氣比接近貧燃料熄火極限，導致燃燒室產(chǎn)生燃燒不穩(wěn)定現(xiàn)象。燃燒不穩(wěn)定現(xiàn)象是目前采用DLN技術(shù)的大型燃機所面臨的共同問題，而且燃燒不穩(wěn)定可能會引起回火、噴嘴熄火，從而導致機組跳機，嚴重時則會損壞燃燒器、火焰筒及透平等昂貴設(shè)備，增加設(shè)備的故障率。要減少燃燒不穩(wěn)定現(xiàn)象，首先需要了解燃燒不穩(wěn)定現(xiàn)象的產(chǎn)生機理，進而檢測燃燒不穩(wěn)定性，根據(jù)檢測結(jié)果對運行工況進行相應(yīng)的調(diào)整，減弱燃燒不穩(wěn)定性帶來的影響，所以研究檢測燃燒不穩(wěn)定現(xiàn)象的方法具有一定的研究意義和價值。

1 燃燒不穩(wěn)定的產(chǎn)生機理

不穩(wěn)定的燃燒過程會發(fā)出聲音，如果聲音在一個自由開放的空間中，則聲音就會輻射殆盡。但是，如果在燃燒室這樣封閉且狹窄的空間中，那么不穩(wěn)定燃燒過程中產(chǎn)生的聲音將會在邊界(燃燒室壁)反射，反射的聲波和燃燒過程相互作用，而燃燒過程對流場比較敏感，則可能產(chǎn)生不穩(wěn)定性現(xiàn)象。

如圖1所示，若在某個時間點由于燃燒過程中氣體的流動和混合過程受到某些因素的影響而產(chǎn)生波動，則會引起熱釋放率的波動，熱釋放率的波動則產(chǎn)生了燃燒不穩(wěn)定現(xiàn)象，在封閉的環(huán)境中生成聲學振動，聲學振動則會使得流動和混合過程繼續(xù)產(chǎn)生波動，以此形成閉環(huán)反饋。

圖1 燃燒不穩(wěn)定性反饋簡圖

國內(nèi)外很多學者對燃燒不穩(wěn)定產(chǎn)生機理的研究表明，燃燒不穩(wěn)定性是由于周期性流體和燃燒過程中周期性地熱釋放率之間的復(fù)雜的相互反饋過程，產(chǎn)生低頻高幅的壓力波動所引起的[6]。

在燃燒不穩(wěn)定的產(chǎn)生機理中，燃燒過程的一個重要特征是壓力振動，而事實上燃燒不穩(wěn)定的外在表現(xiàn)是燃燒室在不同頻率下的壓力振動。通過壓電傳感器測量出燃燒室中的壓力，然后利用信號處理的方法分析信號中頻率的成分。小波變換可以在時域和頻域上同時定位信號，利用小波變換多尺度分析可以提取壓力信號的頻率分量，進而檢測燃燒的不穩(wěn)定性。

2 小波變換簡介

小波變換以傅里葉變換為基礎(chǔ)，其基本思想是在平方可積的函數(shù)空間L2中尋找一組正交基，用這組正交基表示L2中任意的函數(shù)，基函數(shù)可以由預(yù)先選定的基小波(又稱為母小波，Mother Wavelet)通過伸縮平移得到。設(shè)ψ(t)是平方可積函數(shù)，且ψ(t)的傅里葉變換Ψ(ω)滿足以下容許條件：

(1)

則稱ψ(t)為一個基本小波或者母小波。將母小波通過伸縮和平移即可得到小波序列，在連續(xù)的情況下，小波序列為：

(2)

其中a是伸縮因子，τ是平移因子。

對于任意f(t)∈L2(R)，其連續(xù)小波變換為：

(3)

小波變換的時頻窗大小是固定的。也就是說，對于固定的τ，隨著a的增大,小波的時窗就變寬，頻窗就變窄；相反，隨著a的減小，小波變換的時窗就變窄，頻寬就變寬。所以當a較小時，反映的是t = τ附近的高頻成分，當a較大時，反映的是t = τ附近的低頻成分，這樣通過調(diào)節(jié)伸縮因子就可以分析信號的高頻和低頻成分了。

對于離散情況下，采用二進小波序列：

ψj,k(t)=2-j/2ψ(2-jt-k),j,k∈Z

(4)

所以f(t)的離散小波變換為：

(5)

j尺度下重構(gòu)高頻信號為：

(6)

故信號f(t)的高頻分量可以通過式(6)構(gòu)建出來。

3 檢測燃燒不穩(wěn)定性

當燃燒室中低頻、中頻和高頻的壓力振動超過限定幅值時則可能表現(xiàn)出燃燒不穩(wěn)定性。低頻動態(tài)(LowFrequencyDynamics,LFD)壓力振動的頻率分布一般在5～50Hz，常常發(fā)出轟轟低沉的聲音。LFD的幅值會隨著燃燒溫度的降低而增大，所以也稱為“冷振蕩峰值(ColdTone)”[7]，在接近貧燃料熄火極限狀態(tài)容易觀察到。冷振蕩峰值過高容易導致燃燒熱部件和其他重要部件損壞。中頻動態(tài)(MediumFrequencyDynamics,MFD)壓力的振動頻率一般在100～250Hz。因為壓力幅值會隨著燃燒溫度的升高而增大，所以也稱之為“熱振蕩峰值(HotTone)”，MFD對熱部件產(chǎn)生的影響沒有LFD的大。高頻動態(tài)(HighFrequencyDynamics,HFD)壓力振動的振動頻率在250Hz以上，常常伴有“尖叫”。高幅值的HFD嚴重損害機組部件，會在短時間內(nèi)讓部件產(chǎn)生周期疲勞破壞。

用小波分析可以提取壓力脈動信號中的高頻、中頻和低頻分量[8]。圖2為壓力信號的三尺度小波分解。

圖2 壓力信號三尺度分解

壓力信號為S，第1層小波分解將信號大致分解為高頻分量D1和低頻分量A1，A1是信號S大致形狀的近似，而D1則體現(xiàn)了S的細節(jié)部分。將A1繼續(xù)進行小波分解，分解為次高頻分量D2和低頻分量A2，同理，A2反映了A1的大致形狀，D2則體現(xiàn)了A1的細節(jié)部分。將A2繼續(xù)分解為D3和A3。一般對于j尺度的小波分解，S可以分解為：

(7)

j-1尺度的近似分量Aj-1分解為Aj和Dj，Aj是Aj-1的低頻近似分量，而Dj則是Aj-1的高頻細節(jié)分量。

4 結(jié)果分析

4.1 采集系統(tǒng)簡介

某9F型燃氣輪機的燃燒室中采用的是Vibro-MeterCP系列壓電傳感器，其采樣頻率為2 560Hz，安裝在燃燒室法蘭的壓電傳感器將壓力信號轉(zhuǎn)換為電荷信號，電荷信號經(jīng)過解調(diào)器轉(zhuǎn)換為電壓或電流信號，隔離柵起到隔離和保護測量系統(tǒng)的作用，最后信號由分析系統(tǒng)進行處理，流程如圖3所示。

圖3 信號采集和分析流程

4.2 壓力信號分量提取

用三尺度小波分析將壓力信號分解為D1(高頻分量)、D2(中頻分量)、D3(低頻分量)和A3(近似形狀)。由于采用不同的小波在重構(gòu)信號時產(chǎn)生不同的分解重構(gòu)誤差，因此選取精度最高的小波。Db1是小波重構(gòu)信號誤差最小的，故選取Db1小波作為小波分析的母小波[9]。

圖4為調(diào)試電廠燃燒室穩(wěn)定燃燒時的壓力波動曲線，大部分的時間壓力幅值在-40～40mbar(1bar=101kPa)。低頻的壓力脈動主要由燃燒室中熱流體場的波動引起，而高頻的壓力脈動主要由火焰峰面快速脈動引起。

圖4 穩(wěn)定燃燒時的壓力波動曲線

運用小波多尺度分析，將壓力信號分解為D1(高頻分量)、D2(中頻分量)、D3(低頻分量)，如圖5所示。其中高頻分量在平衡值為0處快速上下波動，其幅值主要分布在-20～20mbar；中頻分量也是在平衡值0處快速上下波動，不過其波動幅值比高頻略小，主要分布在-15～15mbar之間；D3反映了低頻分量的細節(jié)部分，其波動比高頻和中頻慢；A3是壓力信號的大致形狀，其波動的速度較慢。

圖5 穩(wěn)定燃燒壓力信號的小波分量

在電廠調(diào)試時，通過改變?nèi)紵鞑煌壷g的燃料分配或者調(diào)節(jié)燃氣輪機進口可調(diào)導葉(InletGuideVane,IGV)改變進氣量可以調(diào)節(jié)燃空比。某段時間在貧燃料極限附近降低燃空比，其壓力波動曲線和小波分量如圖6、7所示。燃燒室中壓力振動平衡值從第2s開始上升，維持大約6s后壓力振動平衡值恢復(fù)到正常水平。在分析結(jié)果中，壓力信號的高頻分量和中頻分量快速上下波動，但是其幅值大部分依舊在-20～ 20mbar。A3則反映了壓力信號的大致形狀，在2～10s之間有明顯的壓力幅值波動；對于低頻分量，能夠檢測出在2～10s之間的波動幅值遠大于其他時間段的壓力幅值，如果低頻幅值超過了規(guī)定限制，那么就會產(chǎn)生燃燒不穩(wěn)定現(xiàn)象。因為壓力低頻振動發(fā)生在貧燃料熄火極限附近，當調(diào)節(jié)IGV增加空氣量的時候，燃空比接近貧燃料熄火極限，燃燒溫度降低，所以低頻壓力幅值增大，運用多尺度小波分析檢測出來的低頻壓力振動符合實際情況。

圖6 受擾動時的壓力波動曲線

圖7 受擾動時的壓力信號小波分量

5 結(jié)束語

介紹了研究燃氣輪機燃燒室不穩(wěn)定性的必要性和燃燒不穩(wěn)定的產(chǎn)生機理，燃燒不穩(wěn)定性的一個重要特征是燃燒室的壓力振動，可以根據(jù)燃燒室壓力振動的特征頻率檢測燃燒不穩(wěn)定性。小波分析能夠同時在時域和頻域上定位壓力信號，并利用多尺度分析分解和重構(gòu)壓力信號，分別提取出壓力信號的高頻分量、中頻分量和低頻分量。利用小波變換提取和分析受擾動的壓力信號的分量，其低頻振動的分析結(jié)果和實際情況吻合。

[1]MooreMJ．NOxEmissionControlinGasTurbinesforCombinedCycleGasTurbinePlant[J]．ProceedingsoftheInstitutionofMechanicalEngineers，PartA：JournalofPowerandEnergy，1997，211(1)：43～52．

[2]GB13223-2011，火電廠大氣污染物排放標準[S].北京：中國環(huán)境科學出版社，2012.

[3]WasloJ，KuwataM，WashamRM．PremixedPilotNozzleforDryLowNOxCombustor[P].US：US4982570，1991-01-08．

[4] 黃素華，陳昱，陸靜，等．干式低NOx燃燒器火焰穩(wěn)定性監(jiān)測技術(shù)探討[J].華東電力，2011，(9)：1533～1537．

[5] 張建，李金科．裂解爐用低NOx燃燒器國內(nèi)技術(shù)進展[J]．化工機械，2012，39(6)：681～685．[6] 曹紅加．預(yù)混火焰燃燒不穩(wěn)定性及其主動控制[D]．北京：中國科學院研究生院(工程熱物理研究所)，2004．

[7] 劉聯(lián)勝，李志勇，佘凱，等．貧氧、預(yù)混火焰熱聲耦合振蕩的聲場分析[J]．燃燒科學與技術(shù)，2010，16(2)：175～180．

[8] 張永懷，白鵬，周金林，等．小波分析與虛擬儀器在紅外氣體分析器中的應(yīng)用[J]．化工自動化及儀表，2003，30(3)：66～68.

[9] 季俊杰，羅永浩．用壓力波動的小波高頻分量檢測燃燒不穩(wěn)定性[J].上海交通大學學報，2007，41(3)：342～346．

ResearchofCombustionInstabilityDetectionMethodBasedonWaveletTransform

ZHUO Zeng-si, HU Li-sheng

(School of Electronic Information and Electrical Engineering, Shanghai Jiaotong University)

Through making use of wavelet transform, the multiscale analysis, decomposition and reconstruction of the combustor’s high frequency, medium frequency and low frequency signals of the gas turbine were implemented, including determination of combustion instability’s location occurred in the time domain. As for the pressure signals with instability, this method can detect the combustion instability that agrees well with real plant operation.

combustion instability, wavelet transform, pressure vibration, multiscale analysis

TH865

A

1000-3932(2017)05-0439-05

卓增泗(1991-)，碩士研究生，從事燃氣輪機燃燒不穩(wěn)定性的研究，1534725079@qq.com。

2016-11-20，

2016-12-19)