曾 文，劉 靖，劉 宇，劉愛虢

（沈陽航空航天大學(xué) 航空航天工程學(xué)院，沈陽 110136）

0 引 言

隨著中國社會經(jīng)濟的發(fā)展，城市數(shù)量、城市人口的增加以及人民生活水平的提高，城市生活垃圾的產(chǎn)量也迅速地增加。預(yù)計在未來的10～15 a內(nèi)，中國城市生活垃圾將以每年 2.3%～3.3%的速度遞增，城市生活垃圾污染問題已成為一個十分嚴峻的社會及環(huán)境問題。如何實現(xiàn)城市生活垃圾的減量化、資源化和無害化，有效控制其對環(huán)境的污染，已成為中國城市化進程中迫切需要解決的一個重大課題。垃圾填埋氣體是垃圾填埋后在缺氧條件下，垃圾中有機物厭氧降解所產(chǎn)生的混合氣體（沼氣），其熱值約為19 256 kJ/m3，屬于低熱值氣體。隨著世界范圍日益緊張的能源趨勢，且由于其具有可循環(huán)性、原料多樣性、環(huán)保性等優(yōu)勢特點，生物質(zhì)氣的資源化利用已成為能源節(jié)約與生態(tài)環(huán)保的重要研究方向。

沼氣中可燃成分主要為CO、H2、CH4等，惰性氣體成分主要為N2、CO2等。但是，由于受諸多因素的影響，沼氣中各組分比例的含量變化明顯。由于各組分比例含量的變化，將對其熱值、著火特性與燃燒特性等產(chǎn)生較為顯著的影響。目前，燃氣輪機的設(shè)計燃料基本為天然氣或輕柴油。以沼氣為燃料的燃氣輪機，在將燃料由天然氣改為沼氣時，燃料成分發(fā)生變化，熱值降低（甲烷熱值為35 916 kJ/m3、柴油熱值為46 000 kJ/kg），且需要滿足環(huán)境保護和燃氣輪機可靠性、經(jīng)濟性要求。因此，需在充分了解沼氣燃燒特性及沼氣組分比例的變化對其燃燒特性的影響規(guī)律[1-4]（燃燒穩(wěn)定性及層流燃燒特性）的基礎(chǔ)上對燃氣輪機的燃燒室結(jié)構(gòu)和運行參數(shù)進行相應(yīng)調(diào)整與改進[5-9]。

Yan等[10]對 4種生物質(zhì)氣在平板火焰燃燒器中的層流燃燒速度進行了試驗測試與數(shù)值計算，同時分析了當量比與氫氣含量對其層流燃燒速度的影響規(guī)律。Vu等[11]在定容彈中試驗測試了 3種生物質(zhì)氣的火焰發(fā)展特性及層流燃燒特性，并對其層流燃燒速度進行了數(shù)值計算。Amell等[12]對高海拔下某種生物質(zhì)氣的層流燃燒特性進行了試驗測試與數(shù)值計算，并分析了初始壓力與當量比等對其層流燃燒特性的影響規(guī)律。Park等[13]對某種生物質(zhì)氣的對沖非預(yù)混燃燒時的燃燒特性進行了試驗測試與數(shù)值計算，并分析了生物質(zhì)氣的組分比例及氮氣稀釋等對其火焰結(jié)構(gòu)及 NOx生成的影響規(guī)律。但是，到目前為止，還無沼氣層流燃燒特性的相關(guān)報道。

1 試驗裝置與原理

在定容燃燒彈中，對 3種不同組分比例比例的沼氣的火焰發(fā)展特性進行了試驗測試，并根據(jù)火焰發(fā)展特性，獲得了 3種不同組成比例沼氣的燃燒穩(wěn)定性與層流燃燒速度。Hu等[14]對本文所采用的定容燃燒彈測試系統(tǒng)及相應(yīng)的試驗原理[15-16]進行了詳細闡述，在此不再贅述。

對沼氣的主要組分比例進行了分析，并根據(jù)組分比例及所占比例范圍配置了3種試驗配氣，如表1所示。從表 1可以看出，在試驗配氣過程中，主要通過改變沼氣中甲烷體積分數(shù)（2種極限情況及中間情況）來改變沼氣中組分比例的含量，從而分析不同組分比例比例對沼氣層流燃燒特性的影響。

表1 沼氣的主要組分比例(體積分數(shù))Table 1 Main compositions of biogas(volume fraction) %

2 試驗結(jié)果與分析

2.1 組分比例對火焰發(fā)展特性的影響

圖1、圖2與圖3分別顯示了初始溫度為320 K，當量比為1.0，初始壓力分別為0.1、0.2、0.3 MPa時，沼氣組分比例對火焰形狀隨時間變化的影響。

碳氫燃料在層流燃燒過程中，會存在優(yōu)先擴散、流體動力學(xué)及浮力不穩(wěn)定等 3種火焰不穩(wěn)定性[17-20]。當馬克斯坦長度的值較小或為負值、混合氣當量比較小時將出現(xiàn)優(yōu)先擴散不穩(wěn)定性，火焰前鋒面呈現(xiàn)不規(guī)則的胞狀結(jié)構(gòu)和火焰突起。在火焰半徑大于30 mm時能觀察到流體動力學(xué)不穩(wěn)定性，火焰前鋒面呈現(xiàn)規(guī)則的胞狀結(jié)構(gòu)。由于本文試驗只考慮火焰半徑在6～25 mm之間的部分，不會顯示該種不穩(wěn)定性。當層流燃燒速度小于 0.15 m/s時將出現(xiàn)浮力不穩(wěn)定性，火焰受到浮力的影響從火焰中心向上飄起，隨著火焰的進一步發(fā)展會出現(xiàn)上半球大下半球小的圖像。本文試驗中，當初始溫度為320 K、初始壓力為0.1與0.2 MPa時，絕大部分工況下沼氣的層流燃燒速度都大于0.15 m/s；當初始壓力為0.3 MPa時，絕大部分工況下沼氣的層流燃燒速度都小于0.15 m/s。因此，本文試驗中在某些工況下將顯現(xiàn)浮力不穩(wěn)定性現(xiàn)象。

圖1 沼氣組分比例對火焰形狀隨時間變化的影響(P=0.1 MPa)Fig.1 Effects of biogas compositions on variations of flame shapes versus time (P=0.1 MPa)

圖2 沼氣組分比例對火焰形狀隨時間變化的影響(P=0.2 MPa)Fig.2 Effects of biogas compositions on variations of flame shapes versus time (P=0.2 MPa)

圖3 沼氣組分比例對火焰形狀隨時間變化的影響(P=0.3 MPa)Fig.3 Effects of biogas compositions on variations of flame shapes versus time (P=0.3 MPa)

如圖1、圖2與圖3所示，當初始溫度為320 K、當量比為1.0、初始壓力為0.1 MPa時，在沼氣中不同甲烷體積分數(shù)下火焰前鋒面均較為光滑，并呈準球形向外發(fā)展。在這些工況下，馬克斯坦長度值均較大，且沼氣的層流燃燒速度均大于0.15 m/s，因此，優(yōu)先擴散與浮力不穩(wěn)定性均不會出現(xiàn)。當初始壓力升高至0.2 MPa時，大部分工況下馬克斯坦長度值低于0.5 mm，火焰前鋒面出現(xiàn)了裂紋與火焰突起，優(yōu)先擴散不穩(wěn)定性初步顯示。當初始壓力繼續(xù)升高至0.3 MPa時，在沼氣中不同甲烷體積分數(shù)下馬克斯坦長度值繼續(xù)降低，優(yōu)先擴散不穩(wěn)定性逐步明顯；同時，由于在此工況下，沼氣中甲烷體積分數(shù)為47%與55.5%時，沼氣的層流燃燒速度均小于0.15 m/s，浮力不穩(wěn)定性將初步顯現(xiàn)。

圖4顯示了初始溫度為320 K，當量比為1.0，初始壓力分別為0.1、0.2、0.3 MPa時，沼氣組分比例對火焰半徑隨時間變化趨勢的影響規(guī)律。

如圖 4所示，在不同初始壓力與不同甲烷體積分數(shù)下火焰半徑隨時間的延長均呈現(xiàn)線性增長趨勢。當初始壓力為0.1與0.3 MPa時，隨著沼氣中甲烷體積分數(shù)的增加，火焰發(fā)展速度明顯加快；但是，當初始壓力為0.2 MPa時，沼氣中甲烷體積分數(shù)由47%升高至55.5%時，火焰發(fā)展速度增加明顯，但當甲烷體積分數(shù)繼續(xù)升高至59%時，火焰發(fā)展速度增加幅度較小。同時，隨著初始壓力的升高，在沼氣中甲烷體積分數(shù)相同的工況下，火焰發(fā)展速度呈現(xiàn)降低趨勢。

圖4 組分比例對火焰半徑r隨時間t變化趨勢的影響(320=T K，0.1=φ)Fig.4 Effects of compositions on variations of flame radius r versus time t (320=TK，0.1=φ)

通過計算式 Sn= d rdt可獲得拉伸火焰?zhèn)鞑ニ俣萐n，其中，r為火焰半徑，t為時間。圖5顯示了初始溫度為320 K，當量比為1.0，初始壓力分別為0.1、0.2、0.3 MPa時，沼氣組分比例對拉伸火焰?zhèn)鞑ニ俣入S火焰半徑變化趨勢的影響規(guī)律。從圖 5可以看出，在不同初始壓力與不同甲烷體積分數(shù)下，拉伸火焰?zhèn)鞑ニ俣入S火焰半徑的增加均呈現(xiàn)升高的趨勢；同時，隨著沼氣中甲烷體積分數(shù)的降低或初始壓力的升高，拉伸火焰?zhèn)鞑ニ俣戎饾u降低。

根據(jù)計算式 Sl- Sn=Lbα，可獲得表怔火焰穩(wěn)定性的馬克斯坦長度（ Lb，在圖6中為直線斜率的相反值）以及無拉伸火焰?zhèn)鞑ニ俣龋?Sl，在圖6中為直線在y軸上的截距），其中，α為火焰拉伸率。圖6顯示了初始溫度為320 K，當量比為1.0，初始壓力分別為0.1、0.2、0.3 MPa時，沼氣組成對拉伸火焰?zhèn)鞑ニ俣入S拉伸率變化趨勢的影響規(guī)律。

從圖 6可以看出，在不同初始壓力與不同甲烷體積分數(shù)下，隨火焰拉伸率的增加，拉伸火焰?zhèn)鞑ニ俣染尸F(xiàn)下降趨勢，說明馬克斯坦長度均為正值；同時，隨著沼氣中甲烷體積分數(shù)的降低或初始壓力的升高，無拉伸火焰?zhèn)鞑ニ俣戎饾u降低。

圖5 沼氣組成成分對拉伸火焰?zhèn)鞑ニ俣萐n隨火焰半徑r變化趨勢的影響(320=TK,0.1=φ)Fig.5 Effects of biogas compositions on variations of stretched flame speed versus flame radius(320=TK,0.1=φ)

圖6 沼氣組成成分對拉伸火焰?zhèn)鞑ニ俣萐n隨拉伸率α變化趨勢的影響(320=TK，0.1=φ)Fig.6 Effects of biogas compositions on variations of stretched flame speed Sn versus flame stretch rate α (320=TK，0.1=φ)

2.2 組分比例對火焰穩(wěn)定性的影響

火焰在向外發(fā)展過程中火焰前鋒面是否穩(wěn)定，主要由馬克斯坦長度的值來反映。馬克斯坦長度為正值時，反映出火焰發(fā)展速度隨拉伸率的增加而減小，當火焰前鋒面出現(xiàn)凸起時（拉伸增加），凸起部分的火焰發(fā)展速度將得到抑制，使火焰趨于穩(wěn)定；反之，馬克斯坦長度為負值時，火焰發(fā)展速度隨拉伸率的增加而增加，當火焰前鋒面出現(xiàn)凸起時，凸起部分的火焰發(fā)展速度將進一步增加，火焰的不穩(wěn)定性增加。

圖7顯示了初始溫度為320 K，初始壓力分別為0.1、0.2、0.3 MPa時，沼氣組分比例對馬克斯坦長度隨當量比變化趨勢的影響。由圖 7可見，隨著當量比的增加、初始壓力的降低或沼氣中甲烷體積分數(shù)的增加，馬克斯坦長度逐漸增大。

圖7 馬克斯坦長度Lb隨當量比的變化趨勢( 320=T K)Fig. 7 Variations of Markstein length Lb with equivalence ratio( 320=T K)

當初始壓力分別為0.2 MPa，當量比為0.8，沼氣中甲烷體積分數(shù)為 47%時，馬克斯坦長度為負值（初始壓力升高至0.3 MPa時，在相同當量比下，沼氣中甲烷體積分數(shù)為55.5%時馬克斯坦長度也為負值），說明此時火焰處于不穩(wěn)定狀態(tài)；其他工況下馬克斯坦長度均為正值，說明這些工況下火焰均處于穩(wěn)定狀態(tài)。

2.3 組分比例對層流燃燒速度的影響

圖8顯示了初始溫度為320 K，初始壓力分別為0.1、0.2、0.3 MPa時，沼氣組分比例對無拉伸火焰?zhèn)鞑ニ俣入S當量比變化趨勢的影響。同時，根據(jù)火焰前鋒面上的質(zhì)量守恒，有 A ρuUl= AρbSl，從而可獲得各工況下層流燃燒速度Ul。式中A為火焰前鋒面面積，ρu、ρb分別為未燃與已燃區(qū)混合氣的密度。

圖9顯示了初始溫度為320 K，初始壓力分別為0.1、0.2、0.3 MPa時，沼氣組分比例對層流燃燒速度隨當量比變化趨勢的影響。

從圖8與圖9可以看出，在不同初始壓力與沼氣中不同甲烷體積分數(shù)下，無拉伸火焰?zhèn)鞑ニ俣扰c層流燃燒速度隨當量比的升高呈現(xiàn)先增加后降低的趨勢，在當量比為1.1時達到最大。同時，在沼氣不同甲烷體積分數(shù)下，無拉伸火焰?zhèn)鞑ニ俣扰c層流燃燒速度隨初始壓力的升高而逐漸降低；在不同初始壓力下，隨著沼氣中甲烷體積分數(shù)的增加，無拉伸火焰?zhèn)鞑ニ俣扰c層流燃燒速度逐漸升高。

表2顯示了初始溫度為320 K，初始壓力分別為0.1、0.2、0.3 MPa、不同當量比與甲烷體積分數(shù)下沼氣的層流燃燒速度。從表2可看出，當沼氣中甲烷體積分數(shù)由47%升高至59%時，沼氣層流燃燒速度均增大20%以上，某些工況（初始壓力為0.3 MPa、當量比為1.3時）下沼氣的層流燃燒速度甚至增大 46%。由此可見，沼氣組分及其所占比例對沼氣層流燃燒速度的影響甚為顯著。因此，在設(shè)計以沼氣為燃料的燃氣輪機時，燃燒室長度的設(shè)計需根據(jù)具體的沼氣成分及其所占比例進行設(shè)計，這樣才能保證沼氣在燃燒室內(nèi)進行充分燃燒，達到高燃燒效率、低污染物排放的效果。

圖8 無拉伸層流火焰?zhèn)鞑ニ俣萐l隨當量比的變化趨勢( 320=T K)Fig.8 Variations of unstretched flame speed Sl with equivalence ratio ( 320=T K)

圖9 層流燃燒速度Ul隨當量比的變化趨勢φ( 320=T K)Fig.9 Variations of laminar burning velocity Ul with equivalence ratio φ ( 320=T K)

表2 不同工況下沼氣的層流燃燒速度Table 2 Laminar burning velocity of biogas at different conditions m·s?1

3 結(jié) 論

試驗測量了當量比范圍為 0.8～1.3，初始壓力范圍為0.1～0.3 MPa，初始溫度為320 K，甲烷體積分數(shù)分別為47%、55.5%、59%的沼氣火焰發(fā)展特性，并分析了組分比例比例對其火焰發(fā)展特性、火焰穩(wěn)定性與層流燃燒速度的影響規(guī)律。得到以下結(jié)論：

1) 當初始溫度為320 K，當量比為1.0，初始壓力為0.1 MPa時，在沼氣中不同甲烷體積分數(shù)下火焰前鋒面均較為光滑，并呈準球形向外發(fā)展。當初始壓力升高至0.2 MPa時，在沼氣中不同甲烷體積分數(shù)下火焰前鋒面均出現(xiàn)了裂紋與火焰突起，優(yōu)先擴散不穩(wěn)定性初步顯示。當初始壓力繼續(xù)升高至0.3 MPa時，優(yōu)先擴散不穩(wěn)定性逐步明顯；同時浮力不穩(wěn)定性將初步顯現(xiàn)。

2) 隨著當量比的增加、初始壓力的降低或沼氣中甲烷體積分數(shù)的增加，馬克斯坦長度逐漸增大；當初始壓力為0.2 MPa，當量比為0.8，沼氣中甲烷體積分數(shù)為47%時，馬克斯坦長度為負值（初始壓力升高至0.3 MPa時，在相同當量比下，沼氣中甲烷體積分數(shù)為55.5%時馬克斯坦長度也為負值），說明此時火焰處于不穩(wěn)定狀態(tài)；其他工況下馬克斯坦長度均為正值，說明這些工況下火焰均處于穩(wěn)定狀態(tài)。

3) 在不同初始壓力與沼氣中不同甲烷體積分數(shù)下，無拉伸火焰?zhèn)鞑ニ俣扰c層流燃燒速度隨當量比的升高呈現(xiàn)先增加后降低的趨勢，在當量比為1.1時達到最大。同時，在不同初始壓力下，隨著沼氣中甲烷體積分數(shù)的增加，無拉伸火焰?zhèn)鞑ニ俣扰c層流燃燒速度逐漸升高。