柴油機(jī)顆粒結(jié)構(gòu)特征與氧化特性分析

沈穎，趙洋，劉帥，劉理凡

(1.江蘇城鄉(xiāng)建設(shè)職業(yè)學(xué)院公用事業(yè)學(xué)院，江蘇 常州 213147；2.常熟理工學(xué)院汽車工程學(xué)院，江蘇 蘇州 215500；3.江蘇大學(xué)汽車與交通工程學(xué)院，江蘇 鎮(zhèn)江 212013)

柴油機(jī)具有熱效率高、動(dòng)力性好等優(yōu)點(diǎn)，在商用車、工程機(jī)械以及船舶等領(lǐng)域得到廣泛的應(yīng)用。柴油機(jī)的顆粒物排放是大氣中PM2.5的主要來源之一，對人體健康具有較大危害[1-2]。針對柴油機(jī)的顆粒物排放，各國相繼出臺一系列排放法規(guī)，且排放法規(guī)日益嚴(yán)格，我國第六階段機(jī)動(dòng)車排放法規(guī)同時(shí)對機(jī)動(dòng)車顆粒的數(shù)量和質(zhì)量排放提出了嚴(yán)格要求[3]。目前，顆粒物捕集器(DPF)是降低柴油機(jī)顆粒物排放的主要措施之一。當(dāng)DPF中顆粒的捕集量達(dá)到一定程度時(shí)，需對DPF中捕集的顆粒進(jìn)行氧化燃燒，實(shí)現(xiàn)DPF的再生。顆粒的幾何尺度、微觀結(jié)構(gòu)對顆粒在DPF中的氧化燃燒具有重要的影響[4-5]。因此，研究顆粒物的結(jié)構(gòu)特征以及氧化特性具有一定的意義。

針對柴油機(jī)排放顆粒的結(jié)構(gòu)特征和氧化特性，國內(nèi)外的學(xué)者開展了許多相關(guān)的研究。D.B.Kittelson等[6]采用N2吸附的方法對輕型柴油機(jī)排放的顆粒進(jìn)行了比表面積分析，結(jié)果表明，柴油機(jī)顆粒的比表面積與碳黑比表面積相當(dāng)，比表面積值約120 m2/g。丁焰等[7]對柴油機(jī)燃用不同燃料時(shí)顆粒的比表面積進(jìn)行了測量與分析，結(jié)果表明，柴油機(jī)燃用生物柴油時(shí)顆粒的比表面積比燃用柴油時(shí)的大。Rockne等[8]采用氣體吸附的方法，對車用和船用柴油機(jī)顆粒的孔隙結(jié)構(gòu)進(jìn)行了測量與分析，研究表明，柴油機(jī)顆粒的孔隙尺寸分布于微孔(孔徑小于2 nm)、中孔(孔徑大于2 nm小于50 nm)和大孔(孔徑大于50 nm)范圍，且船用柴油機(jī)較車用柴油機(jī)顆粒的孔隙數(shù)量和孔容積更多。趙洋等[9]采用小角散射的方法對顆粒的孔隙半徑進(jìn)行了測量與分析，研究表明，顆?？紫栋霃椒植加?～11 nm之間。Alfe等[10]采集了不同條件下預(yù)混火焰中生成的顆粒，并對顆粒的納觀結(jié)構(gòu)與氧化特性進(jìn)行了分析，研究表明，顆粒的納觀結(jié)構(gòu)與氧化特性存在一定的關(guān)系。Li等[11]對柴油機(jī)燃用生物柴油排放的顆粒進(jìn)行了氧化特性分析，研究表明，顆粒在氧化初期結(jié)構(gòu)發(fā)生變化，加劇了顆粒的氧化。馬志豪等[12]對柴油機(jī)顆粒的微觀結(jié)構(gòu)和氧化特性進(jìn)行了研究，結(jié)果表明，顆粒比表面積與其氧化特性有關(guān)，隨著顆粒比表面積增大，顆粒的吸附能力增強(qiáng)，更容易氧化。這些研究工作對柴油機(jī)顆粒的結(jié)構(gòu)特征參數(shù)進(jìn)行了測量，并對顆粒的氧化特性進(jìn)行了分析，定性分析了顆粒結(jié)構(gòu)參數(shù)與氧化特性之間的聯(lián)系，加深了對顆粒氧化特性變化規(guī)律的認(rèn)識，但未更進(jìn)一步地建立顆粒結(jié)構(gòu)特征參數(shù)與氧化特性之間的關(guān)聯(lián)。

針對以上顆粒結(jié)構(gòu)特征與氧化特性研究中存在的不足，本研究通過搭建186FA柴油機(jī)試驗(yàn)臺架，用濾芯式顆粒采集器，采集了柴油機(jī)不同工況下排放的顆粒；采用N2吸附方法，對顆粒的比表面積和孔隙進(jìn)行分析，得到顆粒的比表面積和孔隙結(jié)構(gòu)；采用熱重分析儀分析了顆粒的氧化特性；結(jié)合顆粒的比表面積、孔隙結(jié)構(gòu)與顆粒的氧化特性，探討了顆粒結(jié)構(gòu)特征參數(shù)與顆粒氧化特性的相關(guān)性。

1 顆粒結(jié)構(gòu)特征與氧化動(dòng)力學(xué)參數(shù)

1.1 顆粒比表面積

比表面積為單位質(zhì)量固體所具有的表面面積，單位為m2/g。對于無孔性固體，比表面積單指固體外表面面積；對于孔隙固體，比表面積包含了外表面和孔隙內(nèi)部壁面面積。柴油機(jī)排放的顆粒物為粉末狀多孔介質(zhì)，可采用N2吸附并通過BET方程對物質(zhì)的吸附數(shù)據(jù)進(jìn)行擬合獲得其比表面積[13-14]。

1.2 顆?？紫督Y(jié)構(gòu)

針對顆粒的孔隙結(jié)構(gòu)，可通過孔徑分布、孔徑大小、孔隙容積等參數(shù)對其進(jìn)行定量分析。通過N2吸附，結(jié)合DFT，BJH，DR等方法對顆粒的脫附數(shù)據(jù)進(jìn)行擬合可獲得其孔隙信息。由于DFT方法能夠同時(shí)對孔性固體的微孔(孔半徑小于2 nm)、中孔(孔半徑介于2～50 nm)進(jìn)行分析，因此本研究中采用DFT方法分析顆粒的孔隙[15]。

1.3 特征溫度與活化能

顆粒在氧化過程中的特征溫度在一定程度上可反映顆粒對氧化溫度的敏感性和氧化難易程度。定義顆粒在起始和結(jié)束時(shí)失重速率為-0.001%/℃對應(yīng)的溫度分別為顆粒的起燃溫度Ti和燃盡溫度Th，失重速率峰值對應(yīng)的溫度為峰值溫度Tp。

活化能表示顆粒被氧化時(shí)所需的最低能量。采用Coast-Redfern積分法對顆粒的氧化動(dòng)力學(xué)方程進(jìn)行求解便可得到顆粒被氧化的活化能，其表達(dá)式[16]為

(1)

(2)

(3)

式中：α為顆粒在氧化過程中的失重率；m0為顆粒的初始質(zhì)量；m為顆粒氧化過程中瞬時(shí)質(zhì)量；β為升溫速率；R為摩爾氣體常數(shù)，其值為8.314 J/(mol·K)；T為熱力學(xué)溫度；A1為總指前因子；Sa,0為顆粒的初始表面積；P(O2)為反應(yīng)氣中氧氣的分壓；E為活化能；n為反應(yīng)級數(shù)，其值介于0.5～1，本研究中n=0.7[17-18]。

2 試驗(yàn)方法與設(shè)備

為了分析柴油機(jī)實(shí)際排放顆粒的結(jié)構(gòu)特征和氧化特性，通過搭建的186FA柴油機(jī)試驗(yàn)臺架，采用濾芯式顆粒采集器對柴油機(jī)實(shí)際排放的顆粒進(jìn)行采集。186FA柴油機(jī)結(jié)構(gòu)與性能參數(shù)見表1。最大扭矩工況和標(biāo)定工況是考核柴油機(jī)性能最為典型的兩個(gè)工況點(diǎn)，同時(shí)顆粒排放較多，便于采集。因此，選取了柴油機(jī)最大扭矩工況(100%負(fù)荷，2 700 r/min)、標(biāo)定工況(100%負(fù)荷，3 600 r/min)以及低速滿負(fù)荷工況(100%負(fù)荷，1 500 r/min)3個(gè)工況點(diǎn)，分別代表柴油機(jī)低中高轉(zhuǎn)速工況。通過測功機(jī)控制系統(tǒng)調(diào)節(jié)柴油機(jī)的工況，在排氣管同一位置安裝濾芯式顆粒采集器，對柴油機(jī)上述3個(gè)工況排放的顆粒進(jìn)行采集，每完成一次采集便停機(jī)更換濾芯，采集時(shí)長約60 min，反復(fù)進(jìn)行多次采集。

表1 186FA柴油機(jī)結(jié)構(gòu)與性能參數(shù)

N2吸附法廣泛應(yīng)用于粉末固體比表面積和孔隙的測量與分析。采用N2吸附的方法，利用NOVA3000e比表面積和孔隙分析儀對顆粒進(jìn)行N2等溫吸附和脫附測量，結(jié)合多點(diǎn)BET和DFT方法分析顆粒的N2等溫吸附數(shù)據(jù)，分別獲得顆粒的比表面積和孔隙結(jié)構(gòu)參數(shù)。為保證測量結(jié)果的準(zhǔn)確性，待測顆粒樣品質(zhì)量應(yīng)大于100 mg[19]，用于N2等溫吸附測量的顆粒樣品質(zhì)量分別為198 mg(100%負(fù)荷，1 500 r/min)、136 mg(100%負(fù)荷，2 700 r/min)、152 mg(100%負(fù)荷，3 600 r/min)。為避免顆粒在缸內(nèi)燃燒后期以及排氣過程中吸附的水蒸氣和可溶性有機(jī)物對測量結(jié)果產(chǎn)生影響，進(jìn)行顆粒N2等溫吸附與脫附測量前，通過真空恒溫加熱去除顆?？紫吨幸盐降乃魵?、可溶性有機(jī)物等物質(zhì)，脫氣溫度為300 ℃，脫氣時(shí)間為3 h。

熱重分析法在研究顆粒的氧化特性中得到廣泛的應(yīng)用，可得到顆粒質(zhì)量隨溫度或時(shí)間的變化關(guān)系。采用METTLER TGA/DSCI熱重分析儀對采集的顆粒進(jìn)行氧化特性測量。進(jìn)行顆粒熱重試驗(yàn)時(shí)，顆粒樣品質(zhì)量約3 mg，升溫區(qū)間為40～900 ℃，升溫速率為20 ℃/min；反應(yīng)氣氛圍為10%O2和 90%N2組成的混合氣，反應(yīng)氣流量為100 mL/min。進(jìn)行熱重試驗(yàn)前，按照上述條件首先進(jìn)行一組無顆粒的空白試驗(yàn)，作為對照試驗(yàn)組。

3 結(jié)果分析

3.1 顆粒結(jié)構(gòu)特征分析

表2列出采用BET方程和DFT法分別對顆粒的N2吸附和脫附數(shù)據(jù)擬合得到的顆粒比表面積和孔隙結(jié)構(gòu)。從表2可以看出：與比表面積分析中常用的標(biāo)準(zhǔn)粉末狀碳黑相比，柴油機(jī)排氣顆粒的比表面積與其相當(dāng)；柴油機(jī)排氣顆粒屬于多孔介質(zhì)，從平均孔隙半徑大小來看，顆粒的孔隙屬于中孔范疇(大于2 nm)；在柴油機(jī) 100%負(fù)荷條件下，隨著柴油機(jī)轉(zhuǎn)速的增加，顆粒的比表面積和孔容積均呈增大的趨勢，顆粒的孔隙平均半徑減小。

表2 顆粒的比表面積和孔隙結(jié)構(gòu)

3.2 顆粒氧化特性分析

圖1示出顆粒在氧化過程中的失重(TG)和失重速率(DTG)曲線。從圖1a可以看出，柴油機(jī)在3個(gè)工況下排放的顆粒在低溫階段的質(zhì)量基本保持不變，隨著溫度的升高，顆粒逐漸發(fā)生失重，在溫度超過一定值時(shí)，顆粒的質(zhì)量基本不再發(fā)生變化。柴油機(jī)轉(zhuǎn)速對顆粒的失重特性具有顯著的影響，隨著柴油機(jī)轉(zhuǎn)速的增加，顆粒失重溫度提前。

從圖1b可以看出，柴油機(jī)3個(gè)工況下排放的顆粒在失重過程中呈現(xiàn)單峰失重過程，并未表現(xiàn)出典型的“雙峰”失重過程[20]，即在低溫階段基本顆粒中基本沒有物質(zhì)揮發(fā)或參與氧化，主要是因?yàn)椴蓸游恢脺囟冗^高(顆粒采集于排氣道出口處，溫度在300 ℃以上)，使得廢氣中的HC難以吸附于顆粒表面和孔隙中；隨著氧化溫度的升高，顆粒的失重速率呈先增加后減小的趨勢；隨著柴油機(jī)轉(zhuǎn)速的增加，顆粒的失重速率峰值溫度明顯降低，說明顆粒更容易被氧化。

圖1 顆粒的失重與失重速率曲線

為進(jìn)一步對顆粒的氧化特性進(jìn)行分析，通過顆粒的DTG曲線提取了顆粒在氧化過程中的起燃溫度Ti、峰值氧化溫度Tp以及燃盡溫度Th，結(jié)果見表3。從表3可以看出，隨著柴油機(jī)轉(zhuǎn)速增加，顆粒的起燃溫度、峰值溫度和燃盡溫度均降低；相比于1 500 r/min工況下的顆粒，2 700 r/min和3 600 r/min工況下的顆粒起燃溫度降低了15.3%和10.2%，峰值溫度降低了9.2%和33.6%，燃盡溫度分別降低了17.9%和17.7%。顆粒在缸內(nèi)產(chǎn)生后，在燃燒后期經(jīng)歷了部分的氧化過程，郝斌[21]分析認(rèn)為，隨著柴油機(jī)轉(zhuǎn)速的增加，顆粒在缸內(nèi)停留的絕對時(shí)間變短，導(dǎo)致了顆粒在缸內(nèi)被氧化程度減小，使得顆粒內(nèi)部無序程度增加，因此氧化活性增加，進(jìn)而使顆粒更容易被氧化。高轉(zhuǎn)速下顆粒的比表面積和孔隙容積增大，在氧化過程中氧氣分子與顆粒之間發(fā)生接觸的概率增大；顆?？讖綔p小，有利于顆粒對氧氣分子的吸附，使氧氣分子附著在顆粒內(nèi)部孔隙壁面；因此，高轉(zhuǎn)速下顆粒更容易與氧氣發(fā)生反應(yīng)，最終表現(xiàn)為柴油機(jī)高轉(zhuǎn)速下排放的顆粒氧化活性更高。

表3 顆粒氧化特征溫度

活化能能夠定量地分析顆粒在氧化過程中所需的能量。依據(jù)顆粒起燃溫度至燃盡溫度的熱重試驗(yàn)數(shù)據(jù)，采用式(2)對顆粒失重?cái)?shù)進(jìn)行處理，將ln[(1-(1-α)1-n)/((1-α)T2)]對1/T作圖并進(jìn)行線性擬合，通過擬合直線的斜率得到顆粒的活化能。圖2示出柴油機(jī)在轉(zhuǎn)速1 500 r/min條件下顆粒的氧化動(dòng)力學(xué)參數(shù)擬合曲線，其他工況顆粒擬合情況類似。從圖2可以看出，顆粒在起燃至燃盡反應(yīng)過程中，大致分為兩個(gè)階段。顆粒主要由可溶性有機(jī)物(SOF)、干炭煙(Dry soot)和灰分(Ash)構(gòu)成[22]，從前文的顆粒失重和失重速率曲線可知試驗(yàn)顆粒樣品僅包含了干炭煙和灰分物質(zhì)。其中，灰分為顆粒氧化后剩余的物質(zhì)，所以在起燃溫度至燃盡溫度區(qū)間主要為干炭煙的氧化。干炭煙由碳組成，第一階段反應(yīng)為干炭煙中活性較高的碳與氧氣發(fā)生反應(yīng)，由于活性高，所需的能量較少；第二反應(yīng)階段為活性較低的石墨化碳發(fā)生氧化燃燒，由于石墨化碳有序性高，反應(yīng)所需的能量高。

圖2 顆粒動(dòng)力學(xué)參數(shù)擬合曲線

表4列出對圖2中顆粒失重?cái)?shù)據(jù)處理的線性擬合結(jié)果以及對應(yīng)的活化能。從表4可以看出，顆粒第一反應(yīng)階段的活化能介于40～50 kJ/mol，與孟忠偉等[18]得到的結(jié)果相似，隨著柴油機(jī)轉(zhuǎn)速的增加，顆粒第一階段氧化活化能E1呈增加的趨勢，可能是因?yàn)殡S著柴油機(jī)轉(zhuǎn)速的增加，顆粒在缸內(nèi)停留時(shí)間短，干炭煙中活性較高的碳含量較多，導(dǎo)致這部分物質(zhì)氧化所需的活化能增加。顆粒第二階段反應(yīng)所需的活化能E2變化較大，相比于1 500 r/min工況下排放的顆粒，2 700 r/min和3 600 r/min工況下顆粒的第二階段反應(yīng)活化能減少了21.1%和48.6%，主要是由于顆粒比表面積和孔隙容積增大、孔隙平均半徑減小，使得顆粒更容易與氧氣分子接觸，在高溫的條件下發(fā)生氧化反應(yīng)。

表4 顆?；罨軘M合結(jié)果

3.3 結(jié)構(gòu)特征和氧化特性相關(guān)性分析

偏最小二乘回歸分析方法(PLS)能夠解決變量之間存在多重相關(guān)性的問題，且適用于樣本容量小于變量個(gè)數(shù)的情況[23]。通過PLS方法，以顆粒比表面積、孔隙平均半徑和孔隙容積為自變量，以顆粒第二階段反應(yīng)所需的活化能為因變量，通過oringin 2018軟件數(shù)據(jù)統(tǒng)計(jì)模塊分析3個(gè)自變量對活化能的影響權(quán)重。由于各自變量單位不同，進(jìn)行分析時(shí)需對各變量進(jìn)行無量綱化，其表達(dá)式為

(4)

式中：xi為第i個(gè)變量的樣本值；ui(1)為第i個(gè)變量的樣本均值；si(1)為第i個(gè)變量的樣本標(biāo)準(zhǔn)差。進(jìn)行變量的標(biāo)準(zhǔn)化后，通過PLS統(tǒng)計(jì)分析便可得到顆粒比表面積、孔隙平均半徑以及孔隙容積對顆粒第二反應(yīng)階段活化能的回歸系數(shù)，結(jié)果見圖3。

圖3 標(biāo)準(zhǔn)化PLS模型的回歸系數(shù)

從圖3可以看出，顆粒第二反應(yīng)階段活化能與顆粒比表面積、孔隙容積呈負(fù)相關(guān)，與顆粒的孔隙平均半徑呈正相關(guān)的關(guān)系；顆粒的比表面積、孔隙平均半徑和孔隙容積的回歸系數(shù)絕對值基本相同，均在0.33左右，因此可以認(rèn)為三者對顆粒第二反應(yīng)階段活化能的具有同等的影響權(quán)重，說明顆粒的比表面積、孔隙平均半徑和孔隙容積對顆粒在氧化過程中與氧氣分子的接觸均有重要的作用。為檢驗(yàn)PLS模型的準(zhǔn)確性，以得到的原始數(shù)據(jù)PLS模型對顆粒第二反應(yīng)階段的活化能進(jìn)行預(yù)測，結(jié)果見圖4。從圖4可以看出，不同顆粒樣本的第二反應(yīng)階段的活化能試驗(yàn)測量值與PLS模型預(yù)測值吻合度較高，兩者絕對誤差最大為3.4%，證明PLS模型的分析結(jié)果具有較高的可靠性。

圖4 活化能預(yù)測值與試驗(yàn)測量值的對比

4 結(jié)論

a)柴油機(jī)排氣顆粒為多孔介質(zhì)，隨著柴油機(jī)轉(zhuǎn)速增加，顆粒的比表面積和孔隙容積逐漸增大，顆粒的孔隙平均半徑減小；

b)柴油機(jī)排氣顆粒在低溫階段質(zhì)量基本保持不變，在高溫氛圍下快速失重；隨著柴油機(jī)轉(zhuǎn)速增加，顆粒的氧化溫度向低溫區(qū)間偏移，顆粒的Ti、Tp和Th明顯降低，顆粒更容易被氧化；顆粒在起燃溫度至燃盡溫度區(qū)間的氧化主要分為了兩個(gè)反應(yīng)階段，E1較低，E2較大，隨著柴油機(jī)轉(zhuǎn)速的增加，E1略有增加，E2明顯減?。?/p>

c)根據(jù)PLS分析結(jié)果，顆粒的比表面積、孔隙容積、孔隙平均半徑對顆粒E2具有近似相同的影響權(quán)重，約為33.3%。