李偉峰，劉忠長(zhǎng)，王忠恕，許 允，王建堃

（吉林大學(xué) 汽車(chē)仿真與控制國(guó)家重點(diǎn)實(shí)驗(yàn)室，長(zhǎng)春130022）

0 引 言

天然氣因其具有儲(chǔ)量大、辛烷值高、排放低和價(jià)格便宜等優(yōu)點(diǎn)［1］，已成為一種非常有前途的發(fā)動(dòng)機(jī)替代燃料。雖然排氣再循環(huán)（EGR）技術(shù)可以有效降低發(fā)動(dòng)機(jī)的NOx排放，但是天然氣燃燒初期的火焰?zhèn)鞑ニ俾瘦^低，而EGR 使缸內(nèi)氧濃度降低且混合氣熱容增加［2－4］，火焰?zhèn)鞑ニ俾蔬M(jìn)一步降低，燃料經(jīng)濟(jì)性下降。因此，使用EGR 時(shí)be－NOx往往存在“此消彼長(zhǎng)”的變化關(guān)系。

N2和CO2都是EGR 的主要成分，分別研究N2及CO2稀釋對(duì)天然氣發(fā)動(dòng)機(jī)燃燒過(guò)程的影響，有助于深入分析EGR 成分對(duì)be－NOx“此消彼長(zhǎng)”的 變 化 關(guān) 系 的 貢 獻(xiàn) 率［5－9］。Yasufumi等［5］發(fā) 現(xiàn)N2和CO2稀釋都可以明顯改善柴油／天然氣雙燃料發(fā)動(dòng)機(jī)的NOx－PM“此消彼長(zhǎng)”的變化關(guān)系。Tie等［6］在一臺(tái)柴油機(jī)上開(kāi)展了類似的研究。Huang等［7］發(fā)現(xiàn)CO2稀釋可以降低天然氣發(fā)動(dòng)機(jī)的NOx排放，但是發(fā)動(dòng)機(jī)熱效率降低，HC 增加。Cheolwoong等［8］得出了類似的結(jié)論。Brecq等［9］指出N2和CO2稀釋都可以降低天然氣發(fā)動(dòng)機(jī)的爆燃傾向，但CO2的作用是N2的兩倍。然而，目前還鮮見(jiàn)對(duì)比研究N2和CO2稀釋對(duì)天然氣發(fā)動(dòng)機(jī)燃燒和NOx排放影響方面的報(bào)道。因此，為改善天然氣發(fā)動(dòng)機(jī)be－NOx“此消彼長(zhǎng)”的變化關(guān)系，開(kāi)展此方面的研究十分必要。

工業(yè)制取N2和CO2的技術(shù)已經(jīng)非常成熟［10］，這為N2和CO2稀釋在發(fā)動(dòng)機(jī)上的應(yīng)用奠定了基礎(chǔ)。另外，以空氣為原料，借助發(fā)動(dòng)機(jī)氣缸內(nèi)活塞下行或者廢氣渦輪增壓器產(chǎn)生的負(fù)壓，采用膜分離法可以制取N2并可將其直接引入發(fā)動(dòng)機(jī)的進(jìn)氣系統(tǒng)中。

作者在一臺(tái)由CA6DE柴油機(jī)改裝的渦輪增壓、中冷、重型、SI天然氣發(fā)動(dòng)機(jī)上開(kāi)展了N2和CO2稀釋對(duì)降低NOx排放的研究。

1 試驗(yàn)方案

1.1 試驗(yàn)發(fā)動(dòng)機(jī)

試驗(yàn)在一臺(tái)電控、多點(diǎn)噴射、6 缸、重型天然氣發(fā)動(dòng)機(jī)（CA6SE）上進(jìn)行，發(fā)動(dòng)機(jī)的主要參數(shù)如下：缸徑×行程為106 mm×125 mm；總排量為6.618L；壓縮比為12；標(biāo)定功率為155kW（2300 r·min－1）；最大轉(zhuǎn)矩為700（N·m）（1400／r·min－1）；燃料噴射方式為進(jìn)氣道多點(diǎn)順序噴射；發(fā)動(dòng)機(jī)型式為直列六缸四沖程增壓。該發(fā)動(dòng)機(jī)是在CA6DE柴油機(jī)的基礎(chǔ)上開(kāi)發(fā)的。在開(kāi)發(fā)過(guò)程中，保持原機(jī)基本結(jié)構(gòu)不變，并與原機(jī)保持最大限度的零部件通用性。

在試驗(yàn)中所使用的電控系統(tǒng)具有實(shí)現(xiàn)多點(diǎn)電子控制噴射、稀燃、節(jié)氣門(mén)控制、錯(cuò)誤診斷及通信等功能。發(fā)動(dòng)機(jī)以壓縮天然氣作為燃料。天然氣在壓縮瓶中的壓力為20MPa，通過(guò)減壓閥將壓力減小到0.8 MPa，之后通過(guò)一個(gè)計(jì)量裝置進(jìn)入到進(jìn)氣歧管內(nèi)。為了防止天然氣在減壓過(guò)程中由于體積膨脹而引起“結(jié)冰”現(xiàn)象，在減壓閥上安裝了一個(gè)熱交換器。天然氣的噴射時(shí)刻及火花塞點(diǎn)火時(shí)刻都通過(guò)控制發(fā)動(dòng)機(jī)運(yùn)行狀態(tài)的電腦控制。天然氣的流量通過(guò)一個(gè)質(zhì)量流量計(jì)（MFM）測(cè)量。圖1給出了發(fā)動(dòng)機(jī)的本體、電控系統(tǒng)以及各測(cè)試儀器。

圖1 試驗(yàn)平臺(tái)示意圖Fig.1 Schematic diagram of engine testing bench

試驗(yàn)中使用的N2和CO2的純度為99.9%。N2和CO2的壓力通過(guò)穩(wěn)壓裝置控制并穩(wěn)定在0.2 MPa左右，之后稀釋氣體在進(jìn)氣總管中與空氣混合，并隨空氣進(jìn)入缸內(nèi)。稀釋氣體的流量直接通過(guò)質(zhì)量流量計(jì)（MFM）測(cè)量。表1 給出了N2、CO2、空氣和CH4主要理化性質(zhì)的比較［11］。

表1 稀釋氣體、空氣及燃料的性質(zhì)（0.1 MPa，25 ℃）Table 1 Properties of dilution gases，air and fuel（0.1 MPa，25 ℃）

1.2 試驗(yàn)方案

試驗(yàn)在轉(zhuǎn)速為1450r·min－1、節(jié)氣門(mén)開(kāi)度為23.7%下進(jìn)行，天然氣的流量為12.44kg·h－1，天然氣的噴射時(shí)刻為315°CA BTDC，火花塞點(diǎn)火時(shí)刻為20°CA BTDC。在試驗(yàn)中，以A50（1450r／min，50%負(fù)荷）工況作為基準(zhǔn)點(diǎn)，保持邊界條件（轉(zhuǎn)速、燃料噴射時(shí)刻及脈寬、點(diǎn)火時(shí)刻、節(jié)氣門(mén)開(kāi)度等）不變，在進(jìn)氣中分別添加CO2和N2，直到發(fā)動(dòng)機(jī)轉(zhuǎn)矩下降超過(guò)5%為止，如表2所示。N2和CO2的質(zhì)量流量分別通過(guò)質(zhì)量流量控制裝置進(jìn)行調(diào)節(jié)。

表2 試驗(yàn)工況點(diǎn)Table 2 Test cases examined

在發(fā)動(dòng)機(jī)運(yùn)行一段時(shí)間并達(dá)到穩(wěn)定狀態(tài)后展開(kāi)試驗(yàn)。在試驗(yàn)中，進(jìn)氣中冷后溫度保持在（40±1.5）°C，發(fā)動(dòng)機(jī)冷卻水溫度保持在（80±2）°C，機(jī)油溫度保持在（95±2）°C。另外，由圖1可知，稀釋氣體在中冷器前進(jìn)入進(jìn)氣管內(nèi)并與空氣混合，因此其溫度保持在（40±1.5）°C。

在本研究中用稀釋比表征稀釋氣體對(duì)混合氣的稀釋程度，其定義如下：

2 試驗(yàn)結(jié)果及分析

為了分析不同稀釋氣體對(duì)發(fā)動(dòng)機(jī)燃燒和NOx排放的影響，本文采用稀釋比作為自變量對(duì)結(jié)果進(jìn)行分析［5－8］。另外，為了說(shuō)明N2和CO2稀釋對(duì)混合氣熱物理性質(zhì)、氣體濃度、燃料消耗率、燃燒過(guò)程等影響的差異以及兩者降低NOx排放的原因，本文還采用了NOx排放作為自變量對(duì)結(jié)果進(jìn)行分析。

2.1 NOx 排放特性

圖2給出了不同稀釋比時(shí)NOx排放的變化情況。無(wú)論是使用N2還是CO2，隨著稀釋比增加，NOx排放都明顯降低。由澤爾多維奇定理可知，影響NOx生成的因素主要有缸內(nèi)氧濃度、最高燃燒溫度和燃?xì)庠诟邷?、富氧環(huán)境中的停留時(shí)間［12］。一般而言，只有很少一部分NO 產(chǎn)生于很薄的火焰反應(yīng)帶中，大部分NO 在火焰離開(kāi)后的已燃?xì)怏w中生成［12］。試驗(yàn)中，發(fā)動(dòng)機(jī)節(jié)氣門(mén)開(kāi)度、燃料噴射時(shí)刻及脈寬不變，隨著稀釋氣體增加，由于其稀釋效應(yīng)，部分空氣被稀釋氣體取代，缸內(nèi)氧濃度降低，從而抑制了NOx的生成，如圖3（a）所示；稀釋氣體的熱效應(yīng)使缸內(nèi)混合氣熱容增加，當(dāng)上升相同的溫度時(shí)混合氣需要吸收更多的熱量，因此最高燃燒溫度降低［8］，如圖4（a）和圖5（a）所示。在稀釋效應(yīng)和熱效應(yīng)的共同作用下，NOx排放隨著稀釋比增加而急劇下降［5－8］。

圖2 NOx 排放隨著稀釋比的變化Fig.2 Variation of NOxemissions vs.DR

由圖3（a）和圖4（a）可知，當(dāng)稀釋比相同時(shí)，CO2和N2稀釋所對(duì)應(yīng)的缸內(nèi)氧濃度幾乎相同，而CO2稀釋對(duì)應(yīng)的缸內(nèi)混合氣熱容較高，這說(shuō)明就降低NOx排放而言，相比之下，兩者的稀釋效應(yīng)幾乎相同，但CO2的熱效應(yīng)明顯較大，所以其最高燃燒溫度較低，如圖5（a）所示。因此，當(dāng)稀釋比相同時(shí)，N2降低NOx的能力低于CO2，其所對(duì)應(yīng)的NOx排放相對(duì)更高。相同條件下，為達(dá)到相同NOx排放水平，與CO2稀釋相比，N2的稀釋比更高。

圖3 不同稀釋比和不同NOx排放時(shí)點(diǎn)火前缸內(nèi)平均氧濃度的變化Fig.3 Variation of in－cylinder O2concentration before sparking timing vs.DR and NOxemissions

圖4 不同稀釋比和不同NOx 排放時(shí)混合氣熱容的變化Fig.4 Variation of heat capacity of mixture in cylinder vs.DR and NOxemissions

圖5 不同稀釋比和不同NOx 排放時(shí)缸內(nèi)最高平均燃燒溫度的變化Fig.5 Variation of maximum mean combustion temperature in cylinder vs.DR and NOxemissions

由圖3（b）、圖4（b）和圖5（b）可知，當(dāng)NOx排放降低到相同水平時(shí)，與CO2稀釋相比，N2稀釋對(duì)應(yīng)的氧濃度較低，而混合氣熱容較高，這說(shuō)明相比之下，就降低NOx排放而言，雖然稀釋氣體的稀釋效應(yīng)和熱效應(yīng)都會(huì)降低NOx排放，但N2稀釋抑制NOx生成的主要因素是較低的缸內(nèi)氧濃度，即N2的稀釋效應(yīng)；而CO2稀釋抑制NOx生成的主要因素是較低的缸內(nèi)最高燃燒溫度，即CO2的熱效應(yīng)。由于兩者降低NOx排放的原因不同，因此其對(duì)發(fā)動(dòng)機(jī)燃燒過(guò)程的影響也有所不同。

2.2 CO 和THC排放

圖6給出了不同稀釋比時(shí)CO 和THC 排放的變化情況。CO 是碳?xì)淙剂显谌紵^(guò)程中生成的中間產(chǎn)物。當(dāng)反應(yīng)氣的氧濃度及溫度足夠高，化學(xué)反應(yīng)所占有的時(shí)間足夠長(zhǎng)時(shí)，CO 會(huì)氧化成CO2［12］。由表2可知，隨著N2稀釋比增加，部分空氣被N2取代，過(guò)量空氣系數(shù)減小，缺氧引起燃料燃燒不完全，因此CO 排放增加；另外缸內(nèi)溫度隨著稀釋比增加而下降，局部區(qū)域溫度過(guò)低，使CO 失去溫度條件而不能繼續(xù)氧化為CO2，也會(huì)使CO 排放增加。隨著CO2稀釋比增加，由于CO2的密度較大，其所替代的空氣量較少，天然氣－空氣混合氣略微變濃，但由于缸內(nèi)溫度明顯下降，燃燒持續(xù)期延長(zhǎng)，這就增加了CO 氧化成CO2的化學(xué)反應(yīng)時(shí)間，因此隨著CO2稀釋比增加，CO先略有上升，后呈下降趨勢(shì)。另外，對(duì)于這兩種稀釋氣體，隨著稀釋比增加，混合氣變濃且缸內(nèi)溫度下降，大量CH4未經(jīng)燃燒便排出缸外，因而THC排放明顯上升。與CO2稀釋相比，當(dāng)稀釋比相同時(shí)，N2稀釋引起的CO 和THC 排放較高，這是由于其混合氣較濃的緣故，如表2所示。

圖6 CO 和THC排放隨著稀釋比的變化Fig.6 Variation of CO and THC emissions vs.DR

2.3 比燃料消耗率

圖7 給出了be的變化情況。由圖7（a）可知，隨著CO2稀釋比增加，be明顯升高。在試驗(yàn)中節(jié)氣門(mén)位置不變，因此影響be的因素主要有2個(gè)：混合氣的比熱比和燃燒過(guò)程。隨著CO2流量的增加，一部分空氣被CO2取代，由于CO2的比熱比較低，因此缸內(nèi)混合氣的比熱比下降，如圖8（a）所示。由內(nèi)燃機(jī)理論循環(huán)可知，混合氣的比熱比越小，理論熱效率越低，則be越高［12］。由圖9（a）可知，隨著稀釋比的增加，CA10（燃料化學(xué)能釋放10%所對(duì)應(yīng)的曲軸轉(zhuǎn)角）和CA90（燃料化學(xué)能釋放90%所對(duì)應(yīng)的曲軸轉(zhuǎn)角）后移，混合氣的滯燃期（火花點(diǎn)火時(shí)刻與CA10之間的曲軸轉(zhuǎn)角）和燃燒持續(xù)期（CA10 和CA90 之間的曲軸轉(zhuǎn)角）延長(zhǎng)，燃燒定容度降低，混合氣在膨脹行程做功能力下降，熱效率降低，be升高。

圖7 不同稀釋比和不同NOx 排放時(shí)be的變化Fig.7 Variation of bevs.DR and NOxemissions

圖8 不同稀釋比和不同NOx 排放時(shí)混合氣比熱比的變化Fig.8 Variation of specific heat ratio of charge mixture vs.DR and NOxemissions

由圖7（a）、圖8（a）和圖9（a）可知，當(dāng)稀釋比相同時(shí)，與CO2稀釋相比，N2稀釋時(shí)可獲得較低的be，這是由混合氣較高的比熱比及較快的燃燒速率（燃燒定容度較高）決定的。

由圖7（b）可知，隨著NOx排放降低，be明顯升高，be－NOx呈現(xiàn)出“此消彼長(zhǎng)”的變化關(guān)系。由圖7（b）、圖8（b）和圖9（b）可知，當(dāng)NOx排放降低到相同水平時(shí)，與CO2稀釋相比，N2稀釋時(shí)混合氣的比熱比較高，滯燃期和燃燒持續(xù)期較短，be明顯較低。隨著NOx排放減少，與CO2相比，N2作為稀釋氣體時(shí)，be上升的幅度明顯較低，因此N2稀釋改善了be－NOx“此消彼長(zhǎng)”的變化關(guān)系。

2.4 燃燒過(guò)程

圖9 不同稀釋比和不同NOx 排放時(shí)CA10和CA90的變化Fig.9 Variation of CA10and CA90vs.DR and NOxemissions

圖9 （a）給出CA10和CA90隨著稀釋比的變化情況。隨著CO2稀釋比增加，CA10 和CA90明顯后移，滯燃期和燃燒持續(xù)期延長(zhǎng)。阿雷尼烏斯定律指出，在影響化學(xué)反應(yīng)速率的諸多因素中，溫度對(duì)反應(yīng)速率的影響最為顯著［13］。根據(jù)范特荷夫由試驗(yàn)數(shù)據(jù)歸納的反應(yīng)速率與溫度的近似關(guān)系，當(dāng)溫度升高10K 且其他條件不變的情況下，化學(xué)反應(yīng)速率將增至2～4倍［14］。由圖10（a）可知，隨著CO2稀釋比的增加，火花塞點(diǎn)火時(shí)（20°CA BTDC）的缸內(nèi)平均溫度明顯下降，這是因?yàn)榛旌蠚獾臒崛蓦S著稀釋比增加而增加，如圖4（a）所示。點(diǎn)火時(shí)刻較低的缸內(nèi)溫度不利于天然氣早期火焰發(fā)展，使火焰發(fā)展初期所需時(shí)間延長(zhǎng)，火焰?zhèn)鞑ニ俾式档?，因此CA10和CA90后移，滯燃期和燃燒持續(xù)期延長(zhǎng)，燃燒定容度降低。

由圖9和圖10可知，當(dāng)稀釋比或NOx排放相同時(shí)，與CO2稀釋相比，N2稀釋所對(duì)應(yīng)的CA10和CA90前移，燃燒定容度增加，這是因?yàn)镹2稀釋對(duì)應(yīng)的點(diǎn)火時(shí)刻缸內(nèi)溫度較高。

另外，導(dǎo)溫系數(shù)對(duì)燃燒火焰的傳播有非常大的影響。N2的導(dǎo)溫系數(shù)明顯高于CO2，如圖11所示。由燃燒學(xué)可知，導(dǎo)溫系數(shù)高的氣體有利于燃燒火焰的傳播［15］。因此，當(dāng)NOx排放相同時(shí)，與CO2稀釋相比，N2稀釋時(shí)天然氣燃燒速率較高的另外一個(gè)原因是N2具有較高的導(dǎo)溫系數(shù)。

圖10 不同稀釋比和不同NOx 排放時(shí)點(diǎn)火時(shí)刻（20°CA BTDC）缸內(nèi)平均溫度的變化Fig.10 Variation of mean temperature in cylinder at spark timing（20°CA BTDC）vs.DR and NOxemissions

圖11 N2和CO2的導(dǎo)溫系數(shù)隨著溫度的變化Fig.11 Variation of thermal diffusivities of N2 and CO2vs.temperatures

圖12 給出了缸內(nèi)壓力、放熱率和缸內(nèi)溫度隨著曲軸轉(zhuǎn)角的變化情況。隨著稀釋比增加，無(wú)論是N2稀釋還是CO2稀釋，缸內(nèi)壓力峰值、放熱率峰值和最高平均燃燒溫度都明顯降低，放熱率重心后移，燃燒定容度降低，這是由天然氣的滯燃期和燃燒持續(xù)期延長(zhǎng)引起的。試驗(yàn)結(jié)果說(shuō)明：進(jìn)氣中添加N2或者CO2都可以降低缸內(nèi)的熱負(fù)荷和壓力負(fù)荷，有利于避免早燃、爆燃及表面點(diǎn)火等不正常燃燒現(xiàn)象。

圖12 缸內(nèi)壓力、放熱率及溫度隨著曲軸轉(zhuǎn)角的變化Fig.12 Variation of pressure，heat release rate and temperature in cylinder vs.crank angle

圖13 給出了分別使用N2和CO2作為稀釋氣體將NOx排放降低到相同水平時(shí)缸內(nèi)壓力、放熱率及溫度隨曲軸轉(zhuǎn)角的變化情況。由圖13可知，與CO2稀釋相比，N2稀釋時(shí)缸壓峰值較高，著火延遲期較短，放熱率重心前移，缸內(nèi)溫度峰值較高，這說(shuō)明當(dāng)NOx排放降低到相同水平時(shí)，與CO2稀釋相比，N2稀釋縮短了混合氣的滯燃期和燃燒持續(xù)期，增加了天然氣的燃燒速率，提高了燃燒定容度，從而優(yōu)化了天然氣－空氣混合氣的燃燒過(guò)程。

圖14給出了不同稀釋比和不同NOx排放時(shí)缸內(nèi)壓力峰值及放熱率峰值的變化情況。由圖5和圖14 可知，當(dāng)稀釋比或NOx排放相同時(shí)，與CO2稀釋相比，N2稀釋時(shí)缸內(nèi)壓力峰值、放熱率峰值及溫度峰值較高，這是因?yàn)镹2稀釋時(shí)混合氣的滯燃期和燃燒持續(xù)期較短，這是由混合氣的較低的熱容和較高的導(dǎo)溫系數(shù)引起的。雖然稀釋氣體的稀釋效應(yīng)和熱效應(yīng)都會(huì)降低NOx排放，但由2.1節(jié)可知，N2稀釋抑制NOx生成的主要原因是較低的缸內(nèi)氧濃度（稀釋效應(yīng)），CO2稀釋抑制NOx生成的主要原因是較低的缸內(nèi)溫度（熱效應(yīng)）。已有研究表明：與氧濃度相比，溫度對(duì)混合氣燃燒速率的影響更大［14］。因此，當(dāng)NOx排放降低到相同水平時(shí)，與CO2稀釋相比，N2稀釋對(duì)天然氣－空氣混合氣燃燒速率的影響較小。所以，當(dāng)NOx排放降低到相同水平時(shí)，與CO2稀釋相比，N2稀釋可以優(yōu)化混合氣的燃燒過(guò)程。在選擇降低NOx排放的策略時(shí)，同時(shí)降低缸內(nèi)最高燃燒溫度和氧濃度的策略（如N2稀釋）優(yōu)于過(guò)多降低最高燃燒溫度的策略（如CO2稀釋）。

圖13 NOx 排放降低到相同水平時(shí)缸內(nèi)壓力、放熱率和溫度隨曲軸轉(zhuǎn)角的變化Fig.13 Variation of pressure，heat release rate and temperature in cylinder vs.crank angle at the same NOxemissions

圖14 不同稀釋比和不同NOx 排放時(shí)缸內(nèi)壓力峰值及放熱率峰值的變化Fig.14 Variation of pmaxand ROHRmaxvs.DR and NOxemissions

3 結(jié) 論

（1）無(wú)論是使用N2還是CO2作為稀釋氣體，隨著稀釋比增加，NOx排放都明顯降低。當(dāng)稀釋比相同時(shí)，與使用N2相比，使用CO2能更有效地降低NOx排放。

（2）隨著N2稀釋比的增加，CO 和THC 排放明顯上升；隨著CO2稀釋比的增加，THC 排放明顯上升，但CO 排放先上升后下降。

（3）與CO2稀釋相比，當(dāng)NOx排放降低到相同水平時(shí)，N2稀釋緩和了be－NOx“此消彼長(zhǎng)”的變化關(guān)系。

（4）與CO2稀釋相比，當(dāng)NOx排放相同時(shí)，N2稀釋使CA10和CA90前移，燃燒定容度增加，從而優(yōu)化了混合氣的燃燒過(guò)程。

（5）雖然稀釋氣體的稀釋效應(yīng)和熱效應(yīng)都可以降低NOx排放，但相比之下，N2稀釋降低NOx排放的主要因素是缸內(nèi)較低的氧濃度，即N2的稀釋效應(yīng)；CO2稀釋降低NOx排放的主要因素是較低的缸內(nèi)最高燃燒溫度，即CO2的熱效應(yīng)。

（6）除了控制混合氣氧濃度之外，溫度對(duì)混合氣燃燒過(guò)程的影響不可忽視。在選擇降低NOx排放的策略時(shí)，同時(shí)降低缸內(nèi)燃燒溫度和氧濃度的策略（如N2稀釋）優(yōu)于過(guò)多降低缸內(nèi)燃燒溫度的策略（如CO2稀釋）。

［1］楊曉萍，洪偉，趙梅.天然氣發(fā)動(dòng)機(jī)自適應(yīng)控制軟件抗干擾設(shè)計(jì)［J］.吉林大學(xué)學(xué)報(bào)：工學(xué)版，2004，34（3）：348－351.Yang Xiao－ping，Hong Wei，Zhao Mei.Anti－disturbing design of self－adaptive control softwares for CNG engine［J］.Journal of Jilin University（Engineering and Technology Edition），2004，34（3）：348－351.

［2］Korakianitis T，Namasivayam A M，Crookes R J.Natural－gas fueled spark－ignition（SI）and compression－ignition（CI）engine performance and emissions［J］.Progress in Energy and Combustion Science，2011，37：89－112.

［3］余小草，劉忠長(zhǎng)，王忠恕，等.不同稀釋方式降低天然氣發(fā)動(dòng)機(jī)NOx排放的研究［J］.內(nèi)燃機(jī)工程，2012，33（6）：22－27.Yu Xiao－cao，Liu Zhong－chang，Wang Zhong－shu，et al.Research on NOxemissions reduction of CNG engine by various dilution methods［J］.Chinese Internal Combustion Engine Engineering，2012，33（6）：22－27.

［4］Hu E，Huang Z，Liu B，et al.Experimental study on combustion characteristics of a spark－ignition engine fueled with natural gas－h(huán)ydrogen blends combining with EGR［J］.International Journal of Hydrogen Energy，2009，34（2）：1035－1044.

［5］Yasufumi Y，Eiji K.Influence of intake air dilution with N2or CO2gases on the combustion characteristics of a dual fuel diesel engine with natural gas［C］∥SAE Paper 2013－01－2691.

［6］Tie L，Hiroyuki I，Toshio S，et al.Characterization of low temperature diesel combustion with various dilution gases［C］∥SAE Paper 2007－01－0126.

［7］Huang J，Crookes R J.Assessment of simulated biogas as a fuel for the spark ignition engine［J］.Fuel，1998，77（15）：1793－1801.

［8］Cheolwoong P，Sangyeon W，Changgi K，et al.Effect of mixing CO2with natural gas－h(huán)ydrogen blends on combustion in heavy－duty spark ignition engine［J］.Fuel，2012，102：299－304.

［9］Brecq G，Bellettre J，Tazerout M，et al.Knock prevention of CHP engines by addition of N2and CO2to the natural gas fuel［J］.Apply Thermal Engineering，2003，23：1359－1371.

［10］Reich M，Kapenekas H.Nitrogen purification－pilot plant removal of oxygen［J］.Industrial and Engineering Chemistry，1957，49（5）：869－873.

［11］陳貴堂，王永珍.工程熱力學(xué)［M］.北京：北京理工大學(xué)出版社，2008.

［12］周龍保，劉忠長(zhǎng)，高宗英.內(nèi)燃機(jī)學(xué)［M］.北京：機(jī)械工業(yè)出版社，2010.

［13］嚴(yán)傳俊，范瑋.燃燒學(xué)［M］.西安：西北工業(yè)大學(xué)出版社，2008.

［14］李向容，魏镕，孫柏剛，等.內(nèi)燃機(jī)燃燒科學(xué)與技術(shù)［M］.北京：北京航空航天大學(xué)出版社，2012.

［15］姚春德，居鈺生，劉增勇，等.惰性氣體對(duì)準(zhǔn)均質(zhì)壓燃燃燒過(guò)程的影響［J］.內(nèi)燃機(jī)學(xué)報(bào)，2003，21（4）：205－210.Yao Chun－de，Ju Yu－sheng，Liu Zeng－yong，et al.Study on effect of inertia gases on combustion of quas－i HCCI［J］.Transactions of CSICE，2003，21（4）：205－210.