任昕

摘? ? ? 要： 面對CO2減排的迫切需要和新能源發(fā)展受限的挑戰(zhàn)，集于節(jié)能減排于一體的富氧燃燒技術(shù)應(yīng)運而生。該前沿技術(shù)可通過對溫室氣體進(jìn)行系統(tǒng)控制，從而創(chuàng)造低碳環(huán)保及增產(chǎn)創(chuàng)收的雙贏的局面。利用FLUENT軟件對O2/CO2氛圍下天然氣富氧燃燒進(jìn)行數(shù)值模擬，在最優(yōu)氧氣濃度30%，助燃?xì)怏w溫度為300 K的條件下，選取不同的O2/CO2配比，即設(shè)置CO2的濃度分別為5%、10%、15%、20%、25%、30%、35%、40%，對比分析這八組實驗的燃燒溫度分布特性、組分濃度分布特性、出口平均速度分布特性、污染物NOx的排放水平，從而確定不同的燃燒條件對燃燒特性的影響規(guī)律，找出最佳O2/CO2配比濃度范圍，并對結(jié)果進(jìn)行擬合總結(jié)分析，最終得到天然氣最佳的燃燒條件，使天然氣燃燒效率最高，達(dá)到減少污染物的排放、節(jié)約能源目的。

關(guān)? 鍵? 詞：天然氣;富氧燃燒;燃燒特性

中圖分類號：TF 055? ? ? ?文獻(xiàn)標(biāo)識碼： A? ? ? ?文章編號： 1671-0460（2019）09-2083-04

Abstract： In response to the urgent need for carbon dioxide emission reduction and the challenge of limited development of new energy， oxy-fuel combustion technology emerges as the times require， which integrates energy-saving and emission-reduction. This cutting-edge technology can be used to systematically control greenhouse gases， thereby creating a win-win situation for low-carbon environmental protection and increasing production and income. In this paper， FLUENT software was used to simulate the oxy-fuel combustion of natural gas in O2/CO2 atmosphere. Under the conditions of optimum oxygen concentration of 30% and combustion-supporting gas temperature of 300 K， different O2/CO2 ratios were selected， i.e. the concentration of CO2 was respectively 5%， 10%， 15%， 20%， 25%， 30%， 35%， 40%. The combustion temperature distribution characteristics， component concentration distribution characteristics， outlet average velocity distribution and pollutant NOx emission level of these eight groups of experiments were compared and analyzed. The effect of different combustion conditions on the combustion characteristics was investigated， and the optimum range of O2/CO2 concentration ratio was determined， and a fitting summary and analysis of the results was carried out to obtain the optimum combustion conditions of natural gas， which can make the combustion efficiency of natural gas be the highest， and the emission of pollutants be lowest.

Key words： Oxygen; Enriched combustion; Combustion characteristics

富氧燃燒又稱為空氣分離/煙氣再循環(huán)技術(shù)，該技術(shù)采用比常規(guī)空氣含氧濃度高的富氧空氣進(jìn)行燃燒，用燃燒生成的CO2代替空氣中的N2反復(fù)循環(huán)使用，通過調(diào)整助燃空氣與循環(huán)煙氣的比例控制O2/CO2配比以適應(yīng)不同的燃燒要求[1，2]。富氧燃燒排煙中富含高濃度CO2，便于后續(xù)碳捕集的順利進(jìn)行。相較于傳統(tǒng)燃燒，其優(yōu)越性體現(xiàn)在：

（1）降低NOx的排放量[3]：Sun[4]等分析富氧燃燒中溫度對NOx排放水平的影響，隨著溫度升高及CO2濃度增加，NOx生成量顯著降低;

（2）提高火焰溫度、加快燃燒速度、促進(jìn)燃燒反應(yīng)完全[5]：鐘孝蛟[6]等利用Chemkin軟件模擬O2/CO2，O2/N2，O2/CO2/Ar三種氣氛下直鏈烷烴燃燒的火焰?zhèn)鞑ニ俣龋Y(jié)果表明在相同條件下，甲烷在O2/CO2氣氛中燃燒火焰?zhèn)鞑ニ俣燃s為O2/N2氣氛中燃燒火焰?zhèn)鞑ニ俣鹊?/18，乙烷約為1/7;

（3）增強煙氣的傳熱能力及降低煙黑污染物排放水平：梁軍輝[7]分析富氧程度的提高對擴(kuò)散火焰中煙黑生成規(guī)律的影響，得出氧濃度升高會促進(jìn)燃燒溫度的升高從而抑制煙黑生成量的結(jié)論。

本文利用FLUENT軟件進(jìn)行O2/CO2氛圍下天然氣燃燒的數(shù)值模擬，在最優(yōu)氧氣濃度30%，助燃?xì)怏w溫度為300K的條件下，設(shè)置CO2的濃度分別為5%、10%、15%、20%、25%、30%、35%、40%，選取這八組實驗的燃燒溫度分布特性、組分濃度分布特性、出口平均速度分布特性、污染物NOx的排放水平等進(jìn)行對比分析，確定不同的燃燒條件對燃燒特性的影響規(guī)律，找出最佳O2/CO2配比濃度范圍，對結(jié)果進(jìn)行擬合總結(jié)分析，最終得到天然氣最佳的燃燒條件，使天然氣燃燒效率最高，達(dá)到減少污染物的排放、節(jié)約能源目的。

1? 數(shù)值模擬方法

1.1? 幾何模型及網(wǎng)格劃分

本模擬采用平面圖如圖1所示的軸對稱圓筒形燃燒器：燃燒器壁面溫度保持在300 K，中心處有一長為10 mm的噴嘴，其外徑為12 mm。天然氣以60 m/s的速度經(jīng)噴嘴進(jìn)入燃燒器，助燃?xì)怏w則以0.6 m/s的速度從噴嘴周圍進(jìn)入燃燒器，二者邊混合邊燃燒，產(chǎn)生均勻的湍流擴(kuò)散火焰。燃燒器每間隔0.5 m設(shè)置一個厚度為5 mm、高為0.1 m的擋環(huán)，可有效加速燃燒進(jìn)程。從噴嘴進(jìn)入的天然氣在燃燒器中高速流動，并與低速流動的空氣進(jìn)行混合。

燃燒器網(wǎng)格的精確度和質(zhì)量直接關(guān)系到模擬計算精確與否，利用GAMBIT軟件對燃燒器劃分四邊形單元結(jié)構(gòu)網(wǎng)格，共生成13 882個節(jié)點，27 436個網(wǎng)格。燃燒器噴嘴及擋板附近的氣體流速較快、反應(yīng)進(jìn)行較激烈，故采用細(xì)密網(wǎng)格。在燃燒器的其他區(qū)域，氣體流速較低、反應(yīng)進(jìn)行相對平緩，故采用稀疏網(wǎng)格（圖2）。

1.2? 控制方程及邊界條件

本文選用的天然氣氣質(zhì)條件：CH4 94.48%、C2H6 3.87%、C3H8 0.58%、CO2 0.52%及N2 0.55%。采用流體動力學(xué)基本控制方程[9]，其普遍形式如下：

定義燃燒模型為軸對稱模型，采用標(biāo)準(zhǔn)k-epsilon模型[10]模擬湍流擴(kuò)散運動，選擇Volumetric反應(yīng)、甲烷-空氣組合Mixture Material和Eddy-dissipation（渦流—耗散模型）進(jìn)行化學(xué)反應(yīng)動力學(xué)數(shù)值模擬，并對邊界條件進(jìn)行如表1所示的設(shè)置。

2? 數(shù)值模擬結(jié)果對比分析

2.1? 燃燒溫度分布特性

在O2濃度為30%，助燃?xì)怏w溫度為300 K的富氧燃燒條件下，以每5%為一個區(qū)間，將CO2的濃度范圍從5%調(diào)整至40%，燃燒器的噴嘴處的平均溫度基本上都低于600 K。由于本模型為非預(yù)混燃燒，在噴嘴處助燃?xì)怏w與天然氣剛開始混合，燃燒還未快速放熱。在距離噴嘴約為10 mm的位置溫度開始急劇上升，其原因在于助燃?xì)怏w與天然氣開始快速混合并發(fā)生燃燒反應(yīng)，反應(yīng)持續(xù)放熱（圖3）。

圖3（a）至圖3（e）為變CO2濃度條件下的主要溫度等值線，可以看出，隨著CO2的濃度升高，溫度等值線圖中的高溫區(qū)的形成部位略微縮小，溫度梯度總體變化不大。

圖4為CO2濃度變化對燃燒器平均出口溫度的影響規(guī)律分布圖，在O2濃度一定的情況下，隨著CO2的濃度升高，天然氣的平均出口溫度呈下降趨勢。

產(chǎn)生這一現(xiàn)象的原因在于：

（1）助燃?xì)怏w中的CO2濃度增加，O2濃度不足，使天然氣燃燒不完全，導(dǎo)致燃燒器出口處溫度下降;

（2）CO2的比熱大，吸收的熱量多，因此隨著CO2濃度的升高，散熱逐漸變緩，出口處的溫度下降。

2.2? 組分濃度分布特性

燃燒器出口處的CO2質(zhì)量分?jǐn)?shù)隨著助燃?xì)怏w中的CO2濃度的升高而升高，并且升高趨勢十分明顯。在助燃?xì)怏w中CO2濃度為40%時，燃燒器的出口處CO2質(zhì)量分?jǐn)?shù)接近于50%。

說明助燃空氣中CO2的配比越高，在燃燒反應(yīng)中生成的CO2濃度也越高便于富集和封存，故富氧燃燒采用液化回收大部分的煙氣及少部分煙氣與富氧空氣成比例混合后循環(huán)利用的方法，以實現(xiàn)CO2的資源化利用（圖5）。

2.3? 速度分布特性

選取O2濃度為30%時候的天然氣燃燒的速度矢量圖來研究天然氣的速度分布特性（圖6）。

由于擋板的存在，形成了明顯的回流區(qū)，整個圓筒燃燒器內(nèi)的流動可看成是二維不可壓縮流動，其流動規(guī)律可由伯努利方程和二維連續(xù)性方程進(jìn)行描述。當(dāng)氣體流經(jīng)擋板時，根據(jù)連續(xù)性方程可以推出，空氣的流速增加，同時由伯努利方程得到此時靜壓力將會減少，壓力差由此產(chǎn)生，從而導(dǎo)致空氣加速流動，故產(chǎn)生了回流區(qū)，加速了氣體流動，增加氧氣分子與燃料分子的碰撞概率，從而加快燃料氣體與助燃?xì)怏w的混合速度，更好的促進(jìn)燃燒反應(yīng)的進(jìn)行。

圖7為不同CO2濃度下燃燒器出口平均速度分布圖。

可此看出，隨著CO2的濃度升高，燃燒器面積平均出口速度呈下降趨勢。CO2的摩爾質(zhì)量比O2大，因此CO2的濃度上升使氣體流動速度變緩，導(dǎo)致平均出口速度減小。

2.4? 污染物NOx的排放水平

燃燒排放的污染物NOx中主要成分為NO（至少占比95%），因此NOx排放量可通過燃燒器出口處NO的濃度來替代表示。隨著助燃?xì)怏w中CO2濃度增加，燃燒器出口NO質(zhì)量分?jǐn)?shù)逐漸降低，且下降梯度較為均勻。為降低NOx的生成量，應(yīng)該盡量使CO2與助燃?xì)怏w混合，且在富氧程度一定時，助燃?xì)怏w所占百分比越大，NOx的生成量越少，對環(huán)境越有益[11]，同時驗證了富氧燃燒技術(shù)是一種能綜合控制污染物排放的新型潔凈節(jié)能燃燒技術(shù)（圖8）。

3? 結(jié) 論

本文在閱讀了大量現(xiàn)階段國內(nèi)外文獻(xiàn)的前提下，以fluent模擬燃燒為基礎(chǔ)，對O2/CO2氛圍下天然氣燃燒進(jìn)行數(shù)值模擬分析，在最優(yōu)氧氣濃度30%，助燃?xì)怏w溫度為300 K的富氧燃燒條件下，改變CO2的濃度分別為5%、10%、15%、20%、25%、30%、35%、40%，對天然氣進(jìn)行燃燒模擬，進(jìn)而分析其燃燒特性。主要分析了對天然氣燃燒產(chǎn)生的CO2及NOx的生成量、燃燒平均溫度和出口速度的影響，得出以下結(jié)論：

（1）在O2濃度一定的情況下，隨著CO2的濃度升高，天然氣的平均出口溫度及出口處的速度都呈下降趨勢。

（2）出口處的CO2質(zhì)量分?jǐn)?shù)隨著助燃?xì)怏w中的CO2濃度的升高而升高，并且升高趨勢十分明顯。

（3）NOx的生成量隨著助燃?xì)怏w中CO2與氧氣配比的增加而降低，且下降梯度較為均勻。

參考文獻(xiàn)：

[1]鄭楚光，趙永椿，郭欣.中國富氧燃燒技術(shù)研發(fā)進(jìn)展[J].中國電機工程學(xué)報，2014，34（23）：3856-3864.

[2]李姍. O2/CO2氛圍下天然氣燃燒數(shù)值模擬及熱物性檢測研究[D].長江大學(xué)，2016：3-4.

[3]游卓.富氧燃燒過程中的NOx控制及其系統(tǒng)效率研究[D].浙江大學(xué)， 2013：12-14.

[4]劉宏衛(wèi).氣氛下爐內(nèi)燃燒過程數(shù)值模擬研究[D].保定華北電力大學(xué)，2009：35-37.

[5]Chen R， Axelbaum R. L.Scalar. Dissipation rate at extinction and the effects of oxygen-enriched combustion[J]. Combust flame， 2005， 142： 62-71.

[6]鐘孝蛟，劉豪，趙然，王子劍，胡翰，閆志強，邱建榮.O2/CO2氣氛下火焰?zhèn)鞑ニ俣扔绊懸蛩胤治鯷J].2011.

[7]梁軍輝.氧濃度對擴(kuò)散火焰中煙黑的生成及分布特性的影響研究[D]. 杭州：浙江大學(xué)，2012：47-48.

[8]韓占忠. fluent流體工程仿真計算實例與應(yīng)用[M]. 北京：北京理工大學(xué)出版社， 2004： 233-236.

[9]次英. 流體力學(xué)與傳熱學(xué)基礎(chǔ)[M]. 北京：科學(xué)出版社，2016：16-20.

[10]劉方，翁廟成，龍?zhí)煊? CFD基礎(chǔ)及應(yīng)用[M].重慶：重慶大學(xué)出版社.2015：36-40.

[11]劉暢.天然氣富氧燃燒特性分析[J].當(dāng)代化工，2018，47（04）：854-857.