伏軍 張續(xù)成

摘要：為了加快汽油機微粒捕集器（GPF，gasoline particulate filter）在實際工程中的應(yīng)用，本文以某國產(chǎn)1.5L發(fā)動機匹配GPF后處理系統(tǒng)為例，對微粒捕集器的壓降模型展開研究，提出一種GPF碳載量估算離線標定方法，闡述了GPF碳載量模型的標定內(nèi)容和相關(guān)參數(shù)的標定方法。結(jié)果表明碳載量整體誤差在±0.3g/L以內(nèi)，滿足工程使用要求。

Abstract： In order to speed up the application of GPF （gasoline particulate filter） in practical engineering， this paper takes a domestic 1.5L engine matching GPF after-treatment system as an example， studies the pressure drop model of GPF， puts forward an off-line calibration method of GPF carbon load estimation， and expounds the calibration content of GPF carbon load model and the calibration method of relevant parameters. The results show that the overall error of carbon load is within ±0.3g/L， which meets the engineering requirements.

關(guān)鍵詞：汽油機;GPF;碳載量;標定方法

Key words： gasoline engine;GPF;soot load;calibration method

0? 引言

隨著中國汽車保有量的不斷增加，以及環(huán)境面臨的嚴峻挑戰(zhàn)，環(huán)境保護部在2016年12月發(fā)布了《輕型汽車污染物排放限值及其測量方法（第六階段）》[1]。其中，要求輕型汽車汽油機的微粒排放標準第一階段（2020年7月1日實行）要低于4.5mg/km及6.0×1011PN/kg，第二階段（2023年7月1日開始實行）要低于3mg/km。汽油機加裝微粒捕集器是目前公認的對顆粒物質(zhì)最有效的過濾方法。目前國內(nèi)外對汽油機微粒捕集器的研究只限于對GPF物理特性和碳載量的估算方法展開研究。國外的康寧公司的Per Nicolin，Dominik Rose[2]等人建立了GPF的氧化動力學(xué)模型，并進行了整車試驗驗證;福特公司的Michiel Van Nieuwstadt and Joseph Ulrey，建立了一種開環(huán)碳載量預(yù)估模型，并對GPF相關(guān)的控制策略進行了簡單的闡述[3]。國內(nèi)的研究針對于汽油機微粒捕集器的物理特性，文獻[4-5]分別對GPF的物理模型進行建模與優(yōu)化以提高GPF的捕集效率，并通過試驗數(shù)據(jù)分析了GPF不同涂層、不同的布置方式對發(fā)動機排氣背壓的影響。對于GPF碳載量模型的標定方法國內(nèi)外幾乎沒有研究，而GPF碳載量模型的標定是GPF能否投入實際工程應(yīng)用的關(guān)鍵一步?；诖?，本文研究了微粒捕集器的壓降模型，提出一種GPF碳載量的標定方法，闡述標定內(nèi)容和相關(guān)參數(shù)的標定方法，并設(shè)計試驗方法與評價指標，為GPF快速實現(xiàn)商用化提供一種可行性方案。

1? GPF的碳載量估算方法

汽油機微粒捕集器在實際工作中是安裝在排氣系統(tǒng)中，對從發(fā)動機排出廢氣中的顆粒物、碳煙以及一些難以氧化的灰分進行捕集，隨著碳煙加載量的增加，會導(dǎo)致排氣背壓升高，使發(fā)動機的動力性能下降以及燃油經(jīng)濟性能變差。因此，準確的估算GPF碳載量并選擇合適的時機對碳煙進行再生，不僅可以保證發(fā)動機性能還有利于GPF的使用壽命[6]。綜上可知，對于加裝GPF的汽油發(fā)動機，精確的碳載量估算十分必要。

關(guān)于微粒捕集器的碳載量計算，最直接方法是對GPF進行稱重，但在車輛實際運行中，不可能隨時對GPF拆卸稱重，因此該方法在使用過程中存在諸多不便，在實際中應(yīng)用較少。目前國內(nèi)外對微粒捕集器間接估算方法進行了大量的研究，文獻[7-8]研究了柴油機微粒捕集器（DPF， diesel particulate filter）前后壓差與碳載量的關(guān)系，通過不同流阻下的壓差計算DPF內(nèi)的碳載量，文獻[9-10]研究了基于模型估算碳載量的方法，對DPF的壓降模型進行研究，將簡化后的壓降模型應(yīng)用到GPF碳載量標定中。

由文獻[11]可知，干凈狀態(tài)下的GPF壓降可表示為：

式中：ΔPclean為干凈狀態(tài)下GPF壓降（kPa）;μ為排氣動力粘度（（N·s）/m2）;Qv為排氣體積流量（m3/s）;Vtrap為過濾體體積（m3）;a為過濾體孔道寬度（m）;ωs為過濾體通道壁厚（m）;F為摩擦系數(shù)，其取值28.454;L為過濾體孔道長度（m）;k0為干凈狀態(tài)下壁面滲透率（m2），由試驗得到;ρg為排氣密度（kg/m3）;D為過濾體直徑（m）;ξ為進出口局部壓力損失系數(shù)之和，由試驗得到。

考慮k0與ξ對GPF干凈狀態(tài)下壓降的影響，引入兩個標定系數(shù)C0與C1，并使：

通過公式（2）、（3）分析，k0為干凈狀態(tài)下的壁面滲透率，GPF選定其為定值，因此可知C0為定值;ξ為進出口局部壓力損失系數(shù)之和，與排氣體積流量有關(guān)，為了更加準確的計算干凈狀態(tài)下GPF壓降，因此C1可以表達為關(guān)于排氣體積流量的一維表格，后續(xù)將對C0與C1標定方法進行詳細闡述。

綜上，公式（1）可以簡化為：

公式（4）結(jié)合ρg、Qv可以通過試驗獲取，μ為排氣動力粘度與溫度有關(guān)，本文取GPF入口溫度為Tin（K），則μ可通過下式計算[12]：

本文排氣動力粘度與GPF入口溫度關(guān)系如表1所示。

通過公式（4）可以得到不同發(fā)動機排氣體積流量下GPF干凈狀態(tài)下的壓降，當GPF內(nèi)累積一定量的碳煙，可以由壓差傳感器直接讀出，設(shè)GPF壓差傳感器讀出的壓降為，則：

式中，ΔPsoot為碳煙顆粒產(chǎn)生的壓降（kPa），理論情況下由發(fā)動機讀取的GPF壓差傳感器讀數(shù)會包含由灰分產(chǎn)生的壓降，通過文獻[13]可知，灰分對GPF壓降影響較小，在4000r/min和75%負荷下的前240小時內(nèi)，沒有發(fā)現(xiàn)明顯的排氣背壓升高。在隨后的200小時內(nèi)，在4000r/min的滿載工況下，測量到GPF排氣背壓只增加了30mbar。因此根據(jù)本次試驗周期可以忽略灰分產(chǎn)生的壓降，并假設(shè)GPF壓差傳感器所測量的壓降全部由碳煙顆粒提供。

由上述的ΔPsoot與實際碳載量的關(guān)系可以得出碳載量估算值，但考慮在實際發(fā)動機運行過程中，當GPF入口溫度超過500℃時，碳煙顆粒開始燃燒[14]，為了更加準確的預(yù)估GPF內(nèi)的碳載量引入碳煙指數(shù)（Soot Index），Soot Index可以對GPF入口溫度影響的碳載量進行修正。碳載量指數(shù)表達式如下：

通過公式（4）和公式（7）可知，ΔPclean和Soot Index都是關(guān)于排氣體積流量的函數(shù)，發(fā)動機在不同轉(zhuǎn)速、不同負載下會有不同的排氣體積流量，將多參數(shù)簡化為單一參數(shù)有利于提高標定效率，因此本文以排氣體積流量表征不同的發(fā)動機工況。

對于干凈狀態(tài)下GPF進行碳煙加載，并對加載試驗前后的GPF進行稱重，以得到實際的碳載量，通過試驗采集的數(shù)據(jù)以及壓差傳感器的讀數(shù)可以得到當前碳載量下的Soot Index;重復(fù)上述操作對GPF逐漸進行碳煙加載，直至最大碳載量為止，可以得到不同工況下，碳載量與Soot Index的關(guān)系。

2? GPF碳載量標定試驗

2.1 標定試驗準備

本次試驗是基于某國產(chǎn)四缸汽油發(fā)動機以及與其匹配的GPF后處理系統(tǒng)，試驗發(fā)動機如表2所示，表3列出了試驗所用GPF主要特征參數(shù)。

圖1是GPF在發(fā)動機后處理系統(tǒng)中安裝位置以及傳感器安裝位置示意圖。GPF安裝方式為后置式，即安裝在三元催化器下游750mm處。為了標定還需要在GPF進口和入口分別安裝取氣口以測量GPF實際壓降，最后在GPF入口處安裝溫度傳感器用來采集GPF的入口溫度。

2.2 標定試驗內(nèi)容

前文對GPF碳載量估算方法進行了詳細的闡述，現(xiàn)對相關(guān)參數(shù)的標定試驗內(nèi)容進行介紹。

2.2.1 干凈狀態(tài)下GPF壓降的標定

開始試驗之前，進行發(fā)動機預(yù)熱，當發(fā)動機水溫達到80℃且機油溫度上升至85℃后，為確保GPF為干凈狀態(tài)，設(shè)定發(fā)動機運行轉(zhuǎn)速為3200r/min并減稀空燃比至16，運行30min以使GPF內(nèi)的碳煙完全氧化。完成上述操作后，將發(fā)動機工況設(shè)定為1200r/min、20Nm，待發(fā)動機排氣溫度穩(wěn)定后，設(shè)定發(fā)動機扭矩以20Nm步長至外特性工況點，轉(zhuǎn)速步長為400r/min至最大轉(zhuǎn)速點進行數(shù)據(jù)采集。采集的變量主要有轉(zhuǎn)速、扭矩、排氣體積流量、排氣密度、GPG入口溫度以及GPF壓降等。

通過上述的試驗數(shù)據(jù)可以得到GPF壓差傳感器實測壓降ΔPs與排氣體積流量的關(guān)系如圖2所示，且得到擬合方程為：

在相同的排氣體積流量下，由公式（4）計算得到的GPF壓降等于壓差傳感器的讀數(shù)，因此實測壓降與排氣體積流量的關(guān)系曲線應(yīng)該與干凈狀態(tài)下計算得到的壓降與排氣體積流量的關(guān)系曲線重合。先給C0與C1賦以初始值，再將計算得到干凈狀態(tài)下的壓降與排氣體積流量的關(guān)系也進行擬合，然后調(diào)整標定系數(shù)C0與C1，使計算得到干凈狀態(tài)下壓降與排氣體積流量曲線盡可能與壓差傳感器讀數(shù)與排氣體積流量的關(guān)系曲線貼近，可以得到標定系數(shù)C0為1500，C1與排氣體積流量的關(guān)系如表4所示。

為了驗證C0與C1的標定精度，用Soot Index進行評價，由公式（7）可知，碳載量為零時Soot Index也為零，通過干凈狀態(tài)下的試驗數(shù)據(jù)得到干凈狀態(tài)下的Soot Index與排氣體積流量的關(guān)系如圖3所示，經(jīng)過排氣體積流量與GPF入口溫度的修正后，得到發(fā)動機萬有特性試驗后計算得到的Soot Index基本上都接近于零，結(jié)果表明C0與C1標定精度滿足要求。

2.2.2 GPF碳載量估算及評價指標

通過上文得到的標定系數(shù)C0與C1，當GPF內(nèi)累積一定質(zhì)量的碳煙后，在任意工況下可以通過試驗采集的數(shù)據(jù)和對C0與C1查表可以計算當前工況下干凈狀態(tài)下的壓降，結(jié)合GPF壓差傳感器實測壓降根據(jù)公式（7）可以得到當前的Soot Index。然后對GPF進行碳煙加載，得到不同碳載量與Soot Index的關(guān)系如表5所示。

由表5可知，根據(jù)發(fā)動機的運行工況采集相關(guān)數(shù)據(jù)計算出當前工況下的Soot Index，并根據(jù)Soot Index與實際碳載量的關(guān)系進行查表可以得到當前工況下的碳載量。本文分別以1.322和2.299g/L的實際碳載量對這種碳載量標定方法的準確性進行評估及分析。對碳煙加載后的GPF進行發(fā)動機萬有特性試驗，將采集的數(shù)據(jù)代入公式（7）得到當前碳載量下的Soot Index，圖4表達出了發(fā)動機全工況下的Soot Index與排氣體積流量關(guān)系。

由圖4可知，在不同碳載量的情況下，由發(fā)動機全工況運行后采集的數(shù)據(jù)計算得出的Soot Index呈收斂趨勢，且通過Soot Index平均值查表得到的估算碳載量分別為1.257和2.274g/L，其與實際稱重的碳載量基本相等。此外，通過圖4也可以看出Soot Index存在一定幅度的波動，造成這一情況的原因是，在對加載碳煙后的GPF進行發(fā)動機萬有特性試驗時，由于發(fā)動機工況、空燃比以及排氣溫度等是不斷變化的，GPF內(nèi)的碳煙也在不斷地產(chǎn)生與氧化，所以Soot Index也相應(yīng)在一定范圍內(nèi)變化，通過下文對碳載量的估算偏差可知其總體結(jié)果滿足要求。

圖5為碳載量估算偏差。由圖5可見，在低排氣體積流量范圍內(nèi)，碳載量估算偏差較大，碳載量估算偏差在±0.6g/L以內(nèi)，在其它排氣體積流量范圍內(nèi)，碳載量估算偏差在±0.3g/L范圍內(nèi)。

原因分析：由圖3可知，干凈狀態(tài)下Soot Index在低排氣體積流量區(qū)域內(nèi)偏差較大，其主要是因為GPF壓差傳感器在較低的排氣體積流量范圍內(nèi)誤差較大引起的，考慮到本次試驗使用的GPF壓差傳感器量程較大（0-80kPa），后續(xù)試驗可以采用小量程、高精度的壓差傳感器進行試驗以提高碳載量估算精度。

3? 結(jié)論

綜上所述，本文對GPF碳載量的標定試驗進行研究，提出了一種GPF碳載量的標定方法并闡述了其評價指標，結(jié)果表明基于本文提出的GPF碳載量標定方法得出碳載量估算整體誤差在±0.3g/L以內(nèi)，能夠滿足工程需求，并通過Soot Index對碳載量估算進行了評價。為了進一步提高碳載量估算精度，在具體的試驗過程中可以根據(jù)GPF滿載狀態(tài)下的最大GPF壓降選擇合適的壓差傳感器量程范圍，以提高在低排氣體積流量內(nèi)的碳載量估算精度。

參考文獻：

[1]國家環(huán)保局，國家質(zhì)量監(jiān)督檢驗局.GB18352.6—2016，輕型汽車污染物排放限值及測量方法（中國第六階段）[S].北京：中國環(huán)境科學(xué)出版社，2016.

[2]Per Nicolin， Dominik Rose. Modeling of the Soot Oxidation in Gasoline Particulate Filters[C]. SAE Paper，2015.

[3]Michiel Van Nieuwstadt， Joseph Ulrey. Control Strategies for Gasoline Particulate Filters[J]. SAE Technical Paper 2017-01- 0931.

[4]尹乾熙.缸內(nèi)直噴汽油機微粒捕集器捕集特性研究[D].吉林：吉林大學(xué)，2015.

[5]李配楠.滿足國六標準的汽油機顆粒捕集器_GPF_的試驗研究[D].合肥：合肥工業(yè)大學(xué)，2017.

[6]伏軍.微粒捕集器噴油助燃再生過程熱工參數(shù)建模及控制系統(tǒng)研究[D].湖南大學(xué)，2011.

[7]Naohisa O，Tomohiko N，Shigeki D.? New concept catalyzed DPF for estimating soot loadings from pressure drop [C]// SAE Paper. Detroit，MI， USA，2008， 2008-01-0620.

[8]Dominik R，Thorsten B. Different approaches to soot estimation as key requirement for DPF applications[C]// SAE Paper. Detroit， MI， USA， 2009， 2009-01-1262.

[9]唐蛟，李國祥，王志堅，等.DPF 碳載量模型的建立及試驗[J].內(nèi)燃機學(xué)報，2015，33（1）：51-56.

[10]閆明星.基于模型的微粒捕集器碳煙加載量估計研究[D].長春：吉林大學(xué)，2015.

[11]Michael J，Klingmann R，Spiess? S.，et al. Application of catalyzed gasoline particulate filters to GDI vehicle [C]. SAE Paper 2012-01-1244.

[12]王丹.柴油機微粒捕集器及其再生技術(shù)研究[D].長春：吉林大學(xué)，2013.

[13]A Study of Ash Accumulation in the After-treatment System of a Gasoline Direct Injection Engine Equipped with a Gasoline Particulate Filter[C]. SAE Paper 2017-01-0879.

[14]Shimoda T，Ito Y，Yuuki K，et al. Potential of a low pressure drop filter concept for direct injection gasoline engines to reduce particulate number emission [C]. SAE Paper 2012-01-1241.