JP4333180B2 - Exhaust gas purification system - Google Patents
Exhaust gas purification system Download PDFInfo
- Publication number
- JP4333180B2 JP4333180B2 JP2003086921A JP2003086921A JP4333180B2 JP 4333180 B2 JP4333180 B2 JP 4333180B2 JP 2003086921 A JP2003086921 A JP 2003086921A JP 2003086921 A JP2003086921 A JP 2003086921A JP 4333180 B2 JP4333180 B2 JP 4333180B2
- Authority
- JP
- Japan
- Prior art keywords
- rate
- differential pressure
- dpf
- amount
- purification
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000746 purification Methods 0.000 title claims description 103
- 230000008929 regeneration Effects 0.000 claims description 68
- 238000011069 regeneration method Methods 0.000 claims description 68
- 239000000446 fuel Substances 0.000 claims description 51
- 238000013178 mathematical model Methods 0.000 claims description 37
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000009825 accumulation Methods 0.000 claims description 24
- 238000012360 testing method Methods 0.000 claims description 19
- 238000012937 correction Methods 0.000 claims description 7
- 239000013618 particulate matter Substances 0.000 description 206
- 239000007789 gas Substances 0.000 description 29
- 239000003054 catalyst Substances 0.000 description 23
- 238000005259 measurement Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 13
- 238000000034 method Methods 0.000 description 13
- 238000007254 oxidation reaction Methods 0.000 description 11
- 230000003647 oxidation Effects 0.000 description 10
- 238000002485 combustion reaction Methods 0.000 description 9
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 239000000919 ceramic Substances 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000011144 upstream manufacturing Methods 0.000 description 5
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 4
- 230000001172 regenerating effect Effects 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 3
- 229910000420 cerium oxide Inorganic materials 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 description 3
- BMMGVYCKOGBVEV-UHFFFAOYSA-N oxo(oxoceriooxy)cerium Chemical compound [Ce]=O.O=[Ce]=O BMMGVYCKOGBVEV-UHFFFAOYSA-N 0.000 description 3
- 229910002089 NOx Inorganic materials 0.000 description 2
- MWUXSHHQAYIFBG-UHFFFAOYSA-N Nitric oxide Chemical compound O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 239000012784 inorganic fiber Substances 0.000 description 2
- 238000002844 melting Methods 0.000 description 2
- 230000008018 melting Effects 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052697 platinum Inorganic materials 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- MGWGWNFMUOTEHG-UHFFFAOYSA-N 4-(3,5-dimethylphenyl)-1,3-thiazol-2-amine Chemical compound CC1=CC(C)=CC(C=2N=C(N)SC=2)=C1 MGWGWNFMUOTEHG-UHFFFAOYSA-N 0.000 description 1
- PNEYBMLMFCGWSK-UHFFFAOYSA-N aluminium oxide Inorganic materials [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000000498 cooling water Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000012447 hatching Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910003446 platinum oxide Inorganic materials 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004071 soot Substances 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Landscapes
- Processes For Solid Components From Exhaust (AREA)
- Filtering Of Dispersed Particles In Gases (AREA)
- Exhaust Gas After Treatment (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、連続再生型DPFを備えて、エンジンの排気ガスを浄化する排気ガス浄化システムに関する。
【0002】
【従来の技術】
ディーゼル内燃機関から排出されるパティキュレート(パティキュレート・マター:粒子状物質:以下PMという)の排出量は、NOx,COそしてHC等と共に年々規制が強化されてきており、このPMをDPF(ディーゼルパティキュレートフィルタ:Diesel Particulate Filter )と呼ばれるフィルタで捕集して、外部へ排出されるPMの量を低減する技術が開発されている。
【0003】
このPMを捕集するDPFにはセラミック製のモノリスハニカム型ウオールフロータイプのフィルタや、セラミックや金属を繊維状にした繊維型タイプのフィルタ等があり、これらのDPFを用いた排気ガス浄化システムは、他の排気ガス浄化システムと同様に、内燃機関の排気通路の途中に設置され、内燃機関で発生する排気ガスを浄化して排出している。
【0004】
これらのDPF装置に、CRT(Continuously Regenerating Trap)方式と呼ばれる、DPFの上流側に酸化触媒(Diesel Oxdation Catalyst:DOC)を設けた連続再生型DPFや、CSF(Catalyzed Soot Filter )方式と呼ばれる、フィルタに担持させた触媒の作用によってPMの燃焼温度を低下させ、排気ガスによってPMを焼却する連続再生型DPFがある。
【0005】
このCRT方式の連続再生型DPFは、NO2 (二酸化窒素)によるPMの酸化が、排気ガス中の酸素によりPMを酸化することより低温で行われることを利用したもので、酸化触媒とフィルタとから構成され、この上流側の白金等を担持した酸化触媒により、排気ガス中のNO(一酸化窒素)を酸化してNO2 にして、このNO2 で、下流側のフィルタに捕集されたPMを酸化してCO2 (二酸化炭素)とし、PMを除去している。
【0006】
また、CSF方式の連続再生型DPFは、酸化セリウム(CeO2 )等の触媒を有する触媒付きフィルタで構成され、低温域(300℃〜600℃程度)では、触媒付きフィルタにおける排気ガス中のO2 (酸素)を使用した反応(4CeO2 +C→2Ce2 O3 +CO2 ,2Ce2 O3 +O2 →4CeO2 等)によりPMを酸化し、高温域(600℃程度以上)では、排気ガス中のO2 により直接PMを酸化している。
【0007】
しかしながら、これらの連続再生型DPFにおいても、排気温度が低い場合やNOの排出が少ない内燃機関の運転状態、例えば、内燃機関のアイドル運転や低負荷・低速度運転等の低排気温度状態が継続した場合においては、排気ガス温度が低く触媒の温度が低下して活性化しないため、酸化反応が促進されず、また、NOが不足するので、上記の反応が生ぜず、PMを酸化してフィルタを再生できない。
【0008】
そのため、PMのフィルタへの堆積が継続されて、フィルタの目詰まりが進行すると、排圧が上昇し、燃費の悪化を招き、また、エンジン停止等のトラブルが生じるおそれがある。
【0009】
そこで、これらの連続再生型DPFでは、このPMの蓄積量が予め設定されたPMの蓄積限界値に到達した時に、エンジンの運転状態を再生モード運転に変更して、排気温度を強制的に上昇させたり、NOxの量を増加させたりして、フィルタに捕集されたPMを酸化して除去し、連続再生型DPFを強制再生している。
【0010】
そして、このDPFの再生開始の判断には、DPFの前後の差圧(DPFの圧力損失)が利用される場合もあるが、エンジン回転速度とエンジン負荷から予め記憶したこれらに対応したパティキュレート排出特性とからPM排出量を求め、これを積算して捕集量(PM蓄積量)を求めたり(例えば、特許文献1参照。)、エンジン負荷、回転数、走行距離、走行時間、または燃料消費量の少なくともいずれか一つに基づいて、より具体的にはエンジン負荷と回転数から求めたPM捕集量を走行距離で補正したりして求めた単位時間当たりのPM捕集量を積算して捕集量(PM蓄積量)を求めている(例えば、特許文献2参照。)。
【0011】
また、差圧と差圧を検出した時点までのエンジン回転数の積算値と差圧に対する基準積算値との比に基づいて捕集量(PM蓄積量)を求め、更に、この捕集量にエンジン負荷に応じた補正係数を乗じて補正して捕集量を求めたり(例えば、特許文献3参照。)、差圧と回転数をニュートルネットを応用して捕集量(PM蓄積量)を求めたりしている(例えば、特許文献4参照。)。
【0012】
【特許文献1】
特公平01−23647号公報 (第2頁)
【特許文献2】
特許第2626111号公報 (第3頁、第4頁)
【特許文献3】
特開平08−210122号公報 (第3頁、第4頁)
【特許文献4】
特開平08−284638号公報 (第6頁、第7頁)
【0013】
【発明が解決しようとする課題】
しかしながら、これらのPM蓄積量の算定においては、予めベンチテスト等で求めたマップデータを使用したり、ニュートラルを使用したりするため、データの蓄積やプログラム等のソフト作成に手間が掛かるという問題がある。
【0014】
また、DPFの前後の差圧から、PM捕集量、PM浄化量を数学モデルで推定する方法を使用する場合は、DPFの前後の差圧がDPF周囲の気圧及びPM以外の成分の堆積の影響によっても変化するため、DPFの前後の差圧からDPF内のPM蓄積量を正確に推定することは難しいという問題がある。
【0015】
本発明は、上述の問題を解決するべくなされたものであり、その目的は、エンジンの主要パラメータであるエンジン回転数と燃料流量と、DPF差圧変化率を用いて、エンジンで生成されるPM生成量と化学反応によりDPF内から浄化されるPM浄化量を推測する数学モデルを使用して、簡便で精度良くPM捕集量を推定し、適切な時期に連続再生型DPFの再生を行うことができる排気ガス浄化システムを提供することにある。
【0016】
【課題を解決するための手段】
以上のような目的を達成するための排気ガス浄化システムは、ディーゼルエンジンの排気通路に設けられたパティキュレートを捕集する連続再生型DPFと該連続再生型DPFの前後の差圧を検出する差圧センサと前記連続再生型DPFの再生制御を行う再生制御装置を備えた排気ガス浄化システムにおいて、該再生制御手段は、エンジン回転数及び燃料流量におけるDPF差圧率及びPM生成率の試験データを複数記憶する記憶手段と、前記複数の試験データの中から選択された試験データを用いて数学モデルの係数の決定を行う手段と、前記差圧とエンジン回転数とエンジンの燃料流量を入力するデータ入力手段と、前記差圧の時間に関する変化率であるDPF差圧変化率を算定する差圧変化率算定手段と、前記回転数と前記燃料流量から数学モデルに基づいてPM生成率を算定するPM生成率算定手段と、前記回転数と前記燃料流量から数学モデルに基づいてPM浄化率を算定するPM浄化率算定手段と、前記PM生成率から所定時間の間のPM生成量を算出するPM生成量算定手段と、前記PM浄化率から所定時間の間のPM浄化量を算出するPM浄化量算定手段と、前記PM生成量と前記PM浄化量とから所定時間後のPM蓄積量を算定するPM蓄積量算定手段を備え、前記再生制御装置が、前記DPF差圧変化率が所定の設定値を超える場合は、エンジン回転数、燃料流量、DPF差圧変化率を記録し、このDPF差圧変化率が所定の設定値を超えた回数が前記所定の設定回数を超える場合は、このDPF差圧変化率に応じた範囲の試験データを用いて、前記数学モデルの係数を決定することにより、前記数学モデルの係数を修正する係数修正手段を備えて構成される。
【0017】
この排気ガス浄化システムによれば、連続再生型DPFの再生開始等の判断に使用するPMの堆積量の推定において、エンジン回転数、燃料流量とDPF前後差圧の時間変化率であるDPF差圧変化率とから、PM生成量を算出する数学モデルのPM生成量予測モデル式を用いてPM生成量の時間変化量を求め、また、PM浄化量を算出する数学モデルのPM浄化量予測モデル式を用いてPM浄化量の時間変化量を求めて、これらのPM生成量の時間変化量とPM浄化量の時間変化量と連続再生型DPFの捕集率とから、連続再生型DPFにおけるPM蓄積量を推定する。
【0018】
この数学モデルのPM生成量予測モデル式とPM浄化量予測モデル式の係数は、予め、エンジンの運転条件を表す回転数と燃料流量とDPF差圧変化率とPM生成量の関係を予め測定して、最小二乗法等により決定しておく。
【0019】
つまり、予め、複数のエンジン運転条件(エンジン回転数、燃料流量)に対応したPM生成量及び差圧を測定して係数を決定した数学モデルの多項式を用いて、走行中のDPF差圧変化率からPM生成量、PM浄化量の時間変化量を算定し、更に、これらから、連続再生型DPF内のPM蓄積量を算定する。
【0020】
そして、本発明の排気ガス浄化システムによれば、数学モデルに基づいて、エンジン回転数、燃料流量とDPF差圧変化率とから、PM蓄積量を正確に推定することができるので、適切な時期にDPF強制再生処理制御に入ることができ、DPFの溶損、エンジンの燃費悪化等を防止できる。
【0022】
この係数修正手段を備えることにより、DPF差圧変化率の範囲に応じて係数を設定した数学モデルで、精度良く、PM生成率とPM浄化率を算定できる。
【0023】
また、PM生成率算定手段で、PM生成率を算定する数学モデルを、回転数と燃料流量とに関する二次式とし、PM浄化率算定手段で、PM浄化率を算定する数学モデルを、回転数と燃料流量の積に関する二次式とすることにより、比較的簡単な近似式で、PM生成率とPM浄化率を精度よく算定できる。
【0024】
【発明の実施の形態】
以下、本発明に係る実施の形態の排気ガス浄化システムについて、酸化触媒(DOC)と触媒付きフィルタ(CSF)の組合せで構成される連続再生型DPFを備えた排気ガス浄化システムを例にして図面を参照しながら説明する。
【0025】
図1に、この実施の形態の排気ガス浄化システム1の構成を示す。この内燃機関の排気ガス浄化システム1では、ディーゼルエンジン(内燃機関)10の排気マニホールド11に接続する排気通路12に連続再生型DPF13が設けられている。この連続再生型DPF13は、上流側に酸化触媒13aを下流側に触媒付きフィルタ13bを有して構成される。
【0026】
この酸化触媒13aは、多孔質のセラミックのハニカム構造等の担持体に、白金(Pt)等の酸化触媒を担持させて形成され、触媒付きフィルタ13bは、多孔質のセラミックのハニカムのチャンネルの入口と出口を交互に目封じしたモノリスハニカム型ウオールフロータイプのフィルタや、アルミナ等の無機繊維をランダムに積層したフェルト状のフィルタ等で形成される。このフィルタの部分に白金や酸化セリウム等の触媒を担持する。
【0027】
そして、触媒付きフィルタ13bに、モノリスハニカム型ウオールフロータイプのフィルタを採用した場合には、排気ガスG中のPM(パティキュレートマター:粒子状物質)は多孔質のセラミックの壁で捕集(トラップ)され、繊維型フィルタタイプを採用した場合には、フィルタの無機繊維でPMを捕集する。
【0028】
そして、触媒付きフィルタ13bの再生制御用に、連続再生型DPF13の前後の差圧を検出する差圧センサ21が設けられ、更に、酸化触媒13aと触媒付きフィルタ31bの上流側、中間及び下流側に、それぞれ、上流側温度センサ22、中間温度センサ23及び下流側温度センサ24が設けられる。
【0029】
また、エンジン10の制御用に、エンジン回転数を検出する回転数センサ(図示しない)と、エンジンの負荷を検出するアクセル開度センサ等の負荷センサ(図示しない)が設けられる。
【0030】
これらのセンサの出力値は、エンジン10の運転の全般的な制御を行うと共に、触媒付きフィルタ13bの再生制御も行う制御装置(ECU:エンジンコントロールユニット)30に入力され、この制御装置30から出力される制御信号により、エンジン10の燃料噴射装置14や、吸気マニホールド15への吸気量を調整する吸気弁16や、EGR通路17にEGRクーラ18と共に設けられたEGR量を調整するEGRバルブ19等が制御される。
【0031】
この燃料噴射装置14は燃料ポンプ(図示しない)で昇圧された高圧の燃料を一時的に貯えるコモンレール(図示しない)に接続されており、制御装置30には、エンジンの運転のために、PTOのスイッチのON/OFF,ニュートラルスイッチのON/OFF,車両速度,冷却水温度Tw,エンジン回転数N,アクセル開度Acc等の情報も入力される。
【0032】
そして、排気ガス浄化システム1では、図2に示すように、連続再生型DPF13の再生処理等を制御するための再生制御装置30aが設けられる。この再生制御装置30aは、エンジン10の制御装置30の一部に組み込まれ、連続再生型DPF13の再生制御はエンジン10の制御と並行して行われる。
【0033】
そして、再生制御装置30aには、エンジン10からの測定データとしてエンジン回転数N、燃料流量Qが入力され、連続再生型DPF13からのデータとして、差圧ΔPが入力され、データ処理した後、PM蓄積量PMsが算定され、このPM蓄積量PMsに基づいて、連続再生型DPF13の再生処理開始が判断され、この判断に基づいて再生処理が行われる。
【0034】
この再生制御装置30aは、連続再生型DPF13の再生処理を制御する再生制御手段C1を備えて構成され、図3に示すように、この再生制御手段C1は、再生開始判定手段C10、通常制御運転手段C20、再生処理制御運転手段C30等を備えて構成され、再生開始判定手段C10は、データ入力手段C11、差圧変化率算定手段C12、PM生成率算定手段C13、PM浄化率算定手段C14、PM生成量算定手段C15、PM浄化量算定手段C16、PM蓄積量算定手段C17、及び、係数修正手段C18を備えて構成される。
【0035】
そして、データ入力手段C11は、連続再生型DPF装置13の前後の差圧ΔPとエンジン回転数Nとエンジンの燃料流量Qを入力し、差圧変化率算定手段C12は、差圧ΔPから差圧の時間に関しての変化率であるDPF差圧変化率Rdpを算定し、PM生成率算定手段C13は、エンジン回転数Nと燃料流量Qから数学モデルに基づいてPM生成率Rm1を算定し、PM浄化率算定手段C14は、エンジン回転数Nと燃料流量Qから数学モデルに基づいてPM浄化率Rm2を算定するように構成される。
【0036】
また、PM生成量算定手段C15は、PM生成率Rm1から所定の時刻tk から所定の時間間隔Δtの間のPM生成量ΔM1を算定し、PM浄化量算定手段C16は、PM浄化率Rm2から所定の時刻tk から所定の時間間隔Δtの間のPM浄化量ΔM2を算定し、PM捕集量積算手段C17は、PM生成量ΔM1とPM浄化量ΔM2とから所定の時刻tk から所定の時間間隔Δt後の時刻tk+1 のPM蓄積量PMsを算定するように構成される。
【0037】
そして、係数修正手段C18は、DPF差圧変化率Rdpが所定の設定値Rdp0を超える場合は、エンジン回転数N、燃料流量Q、DPF差圧変化率Rdpを記録し、このDPF差圧変化率Rdpが所定の設定値Rdp0を超えた回数Lが所定の設定回数L0を超える場合は、PM生成率算定手段C13とPM浄化率算定手段C14の数学モデルの係数を修正するように構成される。
【0038】
次に、PM生成率算定手段C13とPM浄化率算定手段C14で使用される数学モデルと、実験データからこれらの係数を推定するための数学モデルについて説明する。
【0039】
この数学モデルは、DPF差圧ΔP、PM生成量M1、PM浄化量M2、PM捕集量M3についての時間変化量Rdp(=dΔP/dt)、Rm1(=dM1/dt)、Rm2(=dM2/dt)、Rm3(=dM3/dt)における、Rm1を推定するPM生成量予測モデル式、及び、Rm2を推定するPM浄化量予測モデル式、及び、実験データからこれらのモデル式の係数を決定するために用いるRdp、Rm1、Rm2、Rm3相互間の関係を示す数式である。
【0040】
先ず、PM生成量予測モデル式に関しては、PM生成量M1の時間変化率であるPM生成率Rm1を、エンジン回転数Nと燃料流量Qをパラメータとし、係数a、b、c、d、eを定数とする二次多項式である(1)式で表す。
【0041】
次に、PM浄化量予測モデル式に関しては、PM浄化量M2の時間変化率であるPM浄化率Rm2を、エンジン回転数Nと燃料流量Qの積(N*Q)をパラメータとし、係数A、B、Cを定数とする二次多項式である(2)式で表す。
【0042】
【数1】
【0043】
一方、特定のエンジン回転数N、燃料流量Qのエンジン運転条件におけるDPF差圧ΔPの時間変化Rdpは、Kを比例定数とする(3)式で表わす。つまり、ここで、α(0<α<1)をDPFの捕集率とし、DPF差圧ΔPの時間変化Rdpは、DPF内でのPM捕集量M3(=M1*α)の時間変化率であるPM捕集率Rm3(=Rm1*α)とPM浄化率Rm2の差に比例するものとする。
【0044】
ここで、PMの浄化が発生しないエンジンの運転条件ではRm2=0であるから(4)式となり、また、(3)式を変形すると(5)式となる。
【0045】
【数2】
【0046】
従って、この(4)式からエンジンの運転条件(N,Q)におけるRdp、Rm1の試験データにより比例定数Kを決定でき、5点以上の異なるエンジン運転条件(N、Q)におけるRdp、Rm1の試験データがあれば、PM生成量予測モデル式(1)の係数a、b、c、d、eを決定できる。なお、DPFの捕集率αは、予め実験等で測定し、求めておく。
【0047】
また、(5)式からPM浄化率Rm2を算定できるので、3点以上の異なる、PM浄化が起こると推定されるエンジン運転条件(N、Q)についてのRdp、Rm1のデータがあれば、PM浄化量予測モデル式(2)の係数A、B、Cを決定できる。
【0048】
そして、この再生制御装置30aの再生制御手段C1で使用される(1)式と(2)式の係数の決定は、図4及び図5の制御フローによって行われ、連続再生型DPF13の再生処理の制御は、図6の制御フローに従って行われる。
【0049】
最初に、数学モデルの(1)式と(2)式の係数決定について説明する。
【0050】
図4のPM生成量予測モデル式(1)の係数決定の制御フローでは、PM生成量予測モデル式である(1)式の係数a,b,c,d,eを、最小二乗法により五元連立一次方程式を立ててこれを解いて求める。なお、連立一次方程式の解法は公知であるので、図4には示さない。
【0051】
この図4の制御フローがスタートすると、ステップS11で、エンジン回転数Ni,燃料流量Qi ,PM生成率Rm1i ,i=1〜nの試験データを入力する。但し、nはデータ数であり、5以上の整数である。
【0052】
次のステップS12では、係数a,b,c,d,eを求めるための五元連立一次方程式(ステップS13の式)の係数等R,S,T,U,V,F,G,X,Y,H,I,J,K,M,W,Zを算定する。但し、式中の記号「Σi 」はi=1〜nの総和を示す。
【0053】
そして、ステップS13では、この五元連立一次方程式を解き、係数a,b,c,d,eを算定する。これにより、PM生成量予測モデル式(1)が決定できるので、この制御フローを終了しリターンする。
【0054】
また、図5のPM浄化量予測モデル式(2)の係数決定の制御フローでは、(5)式を変形した(6)式からPM浄化率とPM生成率の比R12i =Rm 2i /Rm 1i が1以上となった場合をPM浄化が起こるエンジン運転条件と仮定し、このエンジン運転条件について、(5)式からPM浄化率Rm2iを算定し、PM浄化量予測モデル式である(2)式の係数A,B,Cを、最小二乗法により三元連立一次方程式を立ててこれを解いて求める。なろ、連立一次方程式の解法は公知であるので、図5には示さない。
【0055】
【数3】
【0056】
この図5の制御フローがスタートすると、ステップS21で、エンジン回転数Ni ,燃料流量Qi ,PM生成率Rm1i ,i=1〜nの試験データを入力する。但し、nはデータ数であり、5以上の整数である。
【0057】
ステップS22で、n個のデータについて、PM生成率Rm1i とDPF差圧変化率Rdpi と予め試験データにより決定した比例定数Kとから、(6)式によりPM浄化率とPM生成率の比R12i を算定する。
【0058】
ステップS23では、PM浄化率とPM生成率の比R12i が1以上となった場合のデータ、エンジン回転数Nj ,燃料流量Qj ,PM生成率Rm1j ,DPF差圧変化率Rdpj ,j=1〜mを選定する。但し、mは選択されたデータの数であり、3以上でn以下の整数である。
【0059】
ステップS24では、ステップS23で選び出された測定データNj ,Qj ,Rm1j ,Rdpj ,j=1〜mから、(5)式よりPM浄化率Rm2jを計算する。
【0060】
次のステップS25では、エンジン回転数Nj ,燃料流量Qj とステップS24で算定された、PM浄化率Rm2j ,j=1〜mから、PM浄化率予測モデル式(2)の係数A,B,Cを求めるための三元連立一次方程式(ステップS26の式)の係数等R,S,T,F,G,H,Jを算定する。但し、式中の記号「Σj 」はj=1〜mの総和を示す。
【0061】
そして、ステップS26では、この三元連立一次方程式を解き、係数A,B,Cを算定する。これにより、PM浄化量予測モデル式(2)が決定できるので、この制御フローを終了してリターンする。
【0062】
次に、測定データにより各係数を決定されたPM生成量予測モデル式(1)とPM浄化量予測モデル式(2)に基づいて、PM蓄積量を算定する制御について、図6の制御フローに基づいて説明する。
【0063】
この制御フローが再生制御フローから呼ばれてスタートすると、ステップS31で時刻tk をカウントし、時刻tk 及び時間間隔Δtを記録する。なお、この時間間隔Δtはこの制御用の時間間隔であり、計測データNi ,Qi ,Rm1i ,Rdpi のサンプリング間隔とは異なる。
【0064】
ステップS32で時刻tk におけるエンジン回転数Nk ,燃料流量Qk ,DPF差圧ΔPk を検出し、ステップS33で、これらのデータを用いて、時刻tk におけるDPF差圧変化率Rdpk (=ΔPk /Δt)を算定し、また、数学モデル(1)式と(2)式からPM浄化率Rm1k とPM生成率Rm2k を算定する。
【0065】
そして、次のステップS34で、時刻tk から時刻tk+1 (=tk +Δt)のΔt時間内におけるPM生成量ΔM1k (=Rm1k *Δt)とPM浄化量ΔM2k (=Rm2k *Δt)を算定し、ステップS35で、時刻tk+1 におけるPM蓄積量(PMk+1 =PMk +ΔM1k *α−ΔM2k )を算定する。ここで、α(0<α<1)はDPFの捕集率である。
【0066】
また、ステップS36で、DPF差圧変化率Rdpk が、現在の数学モデルの式(1)(2)の係数を決定した設定範囲内(Rdpk ≦Rdp0)にあるかを判定し、設定範囲内であった場合には、そのままステップS38に行き、設定範囲外であった場合には、ステップS37で、設定範囲外の回数Lをカウント(L=L+1)し、その時の測定データであるエンジン回転数Nk ,燃料流量Qk ,DPF差圧変化率Rdpk を記録し、ステップS38に行く。
【0067】
ステップS38では、設定範囲外の回数Lが所定の設定回数L0を超えたか否かを判定し、超えない場合は、そのままリターンし、超えた場合には、ステップS39で、PM生成量予測モデル式(1)の係数決定制御フロー(図4)とPM浄化量予測モデル式(2)の係数決定の制御フロー(図5)を呼んで、各数式モデル(1)式及び(2)式の係数の修正を行ってから、リターンする。
【0068】
以上のPM蓄積量を算定する制御によれば、本数学モデルでは、連続再生型DPF13内のPM生成量M1、PM浄化量M2及びPM捕集量M3の時間変化率Rm1,Rm2,Rm3(=Rm1*α)が、エンジン回転数Nと燃料流量Qを変数とする(1)式、(2)式で算定される。
【0069】
そのため、エンジン回転数N、燃料流量Q、DPF差圧変化率Rdp及びその間の経過時間tを記録した運転履歴情報に基づいて、この数学モデルの(1)式、(2)式を用いることにより、PM生成量M1、PM浄化量M2及びPM捕集量M3の、所定の時間Δt内の各変化量ΔM1、PM浄化量ΔM2及びPM捕集量ΔM3を推定できる。
【0070】
またPM蓄積量(PMの重量)の初期値PMs1 が既知であれば、PM蓄積量PMsの変化量ΔPMs(=ΔM3−ΔM2=ΔM1*α−ΔM2)を積算することにより、その時点tk+1 におけるPM蓄積量PMsK+1 を推定できる。
【0071】
従って、本発明の排気ガス浄化システム1によれば、PM蓄積量PMsを精度良く検知でき、再生開示時期を適切に判定でき、より適切な時期にDPF再生処理制御に入ることができるので、連続再生型DPF13の溶損や燃費の悪化等を回避できる。
【0072】
次に、上記の制御により算定したPM捕集量、PM浄化量の算出例を示す。
【0073】
図7に示したような試験状態において、PM発生量M1及びDPF差圧変化率Rdpの測定試験を実施した。この測定試験では、エンジン回転数Nを一定(2000rpm)とし、DPF入口温度Tinの時間変化が階段状に上昇するように、燃料流量Qを一定の時間間隔で増加し、この時の排気温度、排気流量及びPM発生量M1に応じてDPF差圧ΔPを階段状に上昇させた。この時のDPF入口温度Tin及びDPF差圧ΔPの測定値の一例を図8に示す。
【0074】
このDPF入口温度Tinを一定とした時間帯では、温度レベルがPM浄化反応の起こる温度よりも低ければ、PMがDPFに捕集されるため、DPF差圧ΔPは一定又は増加となる。一方、温度レベルがPM浄化反応が起きる温度以上ならば、DPF差圧ΔPは捕集されるPM量と浄化されるPM量に応じて、DPF内のPM量が減少、平衡、または増加するため、DPF差圧ΔPは減少、平衡、または増加のいずれかとなる。
【0075】
そして、DPF入口温度Tinを一定とした時間帯(10分間)の中間点付近の90秒間におけるDPF差圧変化率Rdpの平均値を、予測式モデル(1)の係数決定用測定データRdpi とした。この試験により得られたPM生成量M1及びDPF差圧変化率Rdpの測定データM1m ,Rdpm を、エンジン回転数N、燃料流量Qをパラメータとするマップ表示で、図9及び図10に示す。
【0076】
これらの測定データM1m ,Rdpm とDPF差圧モデル(4)式より比例定数Kを0.004(mmAq/s)と推定し、この比例定数Kと測定データM1m ,Rdpm からPM生成量予測モデル式(1)の係数を決定した。このPM生成量予測モデル式(1)で算定されたPM生成量M1の計算値M1c をエンジン回転数N、燃料流量Qをパラメータとするマップ表示で図11に示す。
【0077】
また、比例定数Kと測定データM1,Rdpと(5)式からPM浄化量M2の推定M2e を算定した。さらに、PM浄化量予測モデル式(2)の係数を決定し、このPM浄化量予測モデル式(2)で算定された計算値M2c で推定値M2e を補足した結果をエンジン回転数N、燃料流量Qをパラメータとするマップ表示で図12に示す。
【0078】
また、図13に、PM浄化率とPM生成率の比R12をエンジン回転数N、燃料流量Qをパラメータとするマップ表示で示す。この比R12が1よりも大きい場合には、PM生成率Rm1よりPM浄化率Rm2が上回っていると推定できるので、(6)式で算定される比R12が1より大きくなる測定データを選定し、この選定された測定データからPM浄化量予測モデル式(2)の係数を決定している。この比R12が1よりも大きい領域を斜線で示す。
【0079】
【発明の効果】
以上に説明をしたように、本発明の排気ガス浄化システムによれば、数学モデルに基づいて、エンジン回転数、燃料流量とDPF前後差圧の時間変化率であるDPF差圧変化率とから、PM積算量を正確に推定することができるので、連続再生型DPFに蓄積されたPMの量が再生開始のPM蓄積量に達しているか否かを正確に検知でき、適切な時期にDPFの再生処理制御に入ることができるので、DPFの溶損、燃費の悪化等を防止できる。
【0080】
また、DPF差圧変化率が所定の設定値を超える回数が、所定の設定回数を超える場合は、数学モデルの係数を修正する係数修正手段を備えることにより、差圧の変化率の範囲に応じて係数を設定した数学モデルで、精度良く、PM生成率とPM浄化率を算定できる。
【図面の簡単な説明】
【図1】本発明に係る実施の形態の排気ガス浄化システムのシステム構成図である。
【図2】本発明に係る実施の形態のエンジン、連続再生型DPF、計測系及びデータ処理系の構成を示す図である。
【図3】本発明に係る排気ガス浄化システムの制御手段の構成を示す図である。
【図4】PM生成量予測モデル式の係数決定の制御フローの一例を示す図である。
【図5】PM浄化量予測モデル式の係数決定の制御フローの一例を示す図である。
【図6】PM蓄積量の算定の制御フローの一例を示す図である。
【図7】測定試験時におけるエンジン、連続再生型DPF、計測系及びデータ処理系の構成を示す図である。
【図8】DPF入口温度とDPF差圧の実測例を示す図である。
【図9】PM生成量の実測値をマップ表示で示す図である。
【図10】DPF差圧変化率の実測値をマップ表示で示す図である。
【図11】PM生成量の算定値をマップ表示で示す図である。
【図12】PM浄化量の算定値をマップ表示で示す図である。
【図13】PM浄化量とPM生成量の比をマップ表示で示す図である。
【符号の説明】
1 排気ガス浄化システム
10 ディーゼルエンジン(内燃機関)
12 排気通路
13 連続再生型DPF装置
13a 酸化触媒
13b 触媒付きフィルタ(DPF)
21 差圧センサ
30 制御装置(ECU)
30a 再生制御装置
C1 再生制御手段
C10 再生開始判定手段
C11 データ入力手段
C12 差圧変化率算定手段
C13 PM生成率算定手段
C14 PM浄化率算定手段
C15 PM生成量算定手段
C16 PM浄化量算定手段
C17 PM蓄積量算定手段
C18 係数修正手段
C20 通常制御運転手段
C30 再生処理制御運転手段
L 回数
L0 所定の設定回数
N エンジン回転数
PMs PM蓄積量
Q 燃料流量
Rdp DPF差圧変化率
Rm1 PM生成率
Rm2 PM浄化率
Rdp0 所定の設定値
tk 所定の時刻
ΔP 差圧
Δt 所定の時間間隔
ΔM1 PM生成量
ΔM2 PM浄化量[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an exhaust gas purification system that includes a continuous regeneration type DPF and purifies engine exhaust gas.
[0002]
[Prior art]
The amount of particulate matter (particulate matter: hereinafter referred to as PM) emitted from a diesel internal combustion engine is being regulated more and more with NOx, CO, HC, etc., and this PM is converted into DPF (diesel). A technology for reducing the amount of PM collected by a filter called a particulate filter (Diesel Particulate Filter) has been developed.
[0003]
DPFs that collect PM include ceramic monolith honeycomb wall flow type filters, fiber type filters made of ceramic or metal fibers, and exhaust gas purification systems using these DPFs. Like other exhaust gas purification systems, the exhaust gas is installed in the middle of the exhaust passage of the internal combustion engine, and exhaust gas generated in the internal combustion engine is purified and discharged.
[0004]
In these DPF devices, a continuous regeneration type DPF in which an oxidation catalyst (Diesel Oxdation Catalyst: DOC) is provided upstream of the DPF, a filter called a CRT (Continuously Regenerating Trap) system, and a filter called a CSF (Catalyzed Soot Filter) system There is a continuous regeneration type DPF in which the combustion temperature of PM is lowered by the action of a catalyst carried on the catalyst, and PM is incinerated with exhaust gas.
[0005]
This CRT type continuous regeneration type DPF is NO2Utilizes the fact that the oxidation of PM by (nitrogen dioxide) is performed at a lower temperature than the oxidation of PM with oxygen in the exhaust gas. It consists of an oxidation catalyst and a filter, and supports this upstream platinum etc. NO (nitrogen monoxide) in the exhaust gas is oxidized by the oxidized catalyst2And this NO2Then, the PM collected in the downstream filter is oxidized to CO.2(Carbon dioxide) and PM is removed.
[0006]
Also, the CSF type continuous regeneration type DPF is cerium oxide (CeO).2) In the exhaust gas in the filter with a catalyst in a low temperature range (about 300 ° C. to 600 ° C.).2Reaction using (oxygen) (4CeO2+ C → 2Ce2OThree+ CO2, 2Ce2OThree+ O2→ 4CeO2Etc.) to oxidize PM, and in a high temperature range (about 600 ° C. or higher), O in exhaust gas2This directly oxidizes PM.
[0007]
However, even in these continuous regeneration type DPFs, when the exhaust temperature is low or the operating state of the internal combustion engine with low NO emission, for example, the low exhaust temperature state such as idle operation or low load / low speed operation of the internal combustion engine continues. In this case, since the exhaust gas temperature is low and the temperature of the catalyst is not lowered and activated, the oxidation reaction is not promoted, and NO is insufficient, so that the above reaction does not occur and the PM is oxidized and filtered. Cannot be played.
[0008]
Therefore, when PM is continuously deposited on the filter and the filter is clogged, the exhaust pressure increases, fuel consumption is deteriorated, and troubles such as engine stop may occur.
[0009]
Therefore, in these continuous regeneration type DPFs, when the accumulated amount of PM reaches a preset accumulation limit value of PM, the engine operating state is changed to regeneration mode operation, and the exhaust temperature is forcibly increased. Or increasing the amount of NOx to oxidize and remove the PM trapped in the filter to forcibly regenerate the continuously regenerating DPF.
[0010]
In some cases, the differential pressure before and after the DPF (pressure loss of the DPF) is used to determine the start of regeneration of the DPF, but the particulate discharge corresponding to these stored in advance from the engine rotational speed and the engine load. The PM emission amount is obtained from the characteristics, and this is integrated to obtain the collection amount (PM accumulation amount) (see, for example, Patent Document 1), engine load, rotation speed, travel distance, travel time, or fuel consumption. Based on at least one of the quantities, more specifically, the PM collection amount per unit time obtained by correcting the PM collection amount obtained from the engine load and the rotational speed with the travel distance is integrated. Thus, the amount collected (PM accumulated amount) is obtained (see, for example, Patent Document 2).
[0011]
In addition, the amount of collected (PM accumulated amount) is obtained based on the ratio between the differential pressure and the integrated value of the engine speed up to the time when the differential pressure is detected and the reference integrated value with respect to the differential pressure. The correction amount is multiplied by a correction coefficient corresponding to the engine load to obtain a collection amount (see, for example, Patent Document 3), or the differential pressure and the rotational speed are applied using a neutron net to collect the amount (PM accumulation amount). (For example, refer to Patent Document 4).
[0012]
[Patent Document 1]
Japanese Patent Publication No. 01-23647 (Page 2)
[Patent Document 2]
Japanese Patent No. 2626111 (
[Patent Document 3]
JP 08-210122 A (
[Patent Document 4]
JP 08-284638 A (
[0013]
[Problems to be solved by the invention]
However, in the calculation of the amount of accumulated PM, the map data previously obtained by a bench test or the like is used, or the neutral is used. is there.
[0014]
In addition, when using a method for estimating the amount of PM collected and the amount of PM purification using a mathematical model from the differential pressure before and after the DPF, the differential pressure before and after the DPF is caused by the pressure around the DPF and the accumulation of components other than PM. Since it also changes due to the influence, there is a problem that it is difficult to accurately estimate the PM accumulation amount in the DPF from the differential pressure before and after the DPF.
[0015]
The present invention has been made to solve the above-described problems, and its object is to generate PM generated by an engine using engine revolution speed, fuel flow rate, and DPF differential pressure change rate, which are main engine parameters. Use a mathematical model to estimate the amount of PM purified from the DPF by the amount of product and chemical reaction, estimate the amount of PM trapped easily and accurately, and regenerate the continuously regenerating DPF at an appropriate time It is an object of the present invention to provide an exhaust gas purification system capable of achieving the above.
[0016]
[Means for Solving the Problems]
An exhaust gas purification system for achieving the above object is a difference between a continuous regeneration type DPF that collects particulates provided in an exhaust passage of a diesel engine and a differential pressure that detects a differential pressure before and after the continuous regeneration type DPF. In the exhaust gas purification system including a pressure sensor and a regeneration control device that performs regeneration control of the continuous regeneration type DPF, the regeneration control means includes:A storage means for storing a plurality of DPF differential pressure rate and PM generation rate test data at an engine speed and a fuel flow rate, and a coefficient of a mathematical model are determined using test data selected from the plurality of test data. Means,Data input means for inputting the differential pressure, engine speed and engine fuel flow rate, differential pressure change rate calculating means for calculating a DPF differential pressure change rate which is a change rate with respect to time of the differential pressure, and the rotational speed PM generation rate calculation means for calculating a PM generation rate from the fuel flow rate based on a mathematical model, and based on the mathematical model from the rotation speed and the fuel flow ratePM purification rateA PM purification rate calculating means for calculating a PM generation amount calculating means for calculating a PM generation amount for a predetermined time from the PM generation rate, and a PM for calculating a PM purification amount for a predetermined time from the PM purification rate Purification amount calculation means, and PM accumulation amount calculation means for calculating a PM accumulation amount after a predetermined time from the PM generation amount and the PM purification amount, and the regeneration control device has a predetermined rate of change in DPF differential pressure. When the set value is exceeded, the engine speed, the fuel flow rate, and the DPF differential pressure change rate are recorded, and when the number of times that the DPF differential pressure change rate exceeds the predetermined set value exceeds the predetermined set number of times,By determining the coefficient of the mathematical model using test data in a range corresponding to the DPF differential pressure change rate,Coefficient correcting means for correcting the coefficient of the mathematical model is provided.
[0017]
According to this exhaust gas purification system, in estimating the amount of accumulated PM used to determine the start of regeneration of a continuous regeneration type DPF, the DPF differential pressure that is the rate of change over time of the engine speed, the fuel flow rate, and the differential pressure across the DPF Based on the rate of change, the PM generation amount prediction model formula of the mathematical model for calculating the PM generation amount is used to obtain the time change amount of the PM generation amount, and the PM purification amount prediction model formula of the mathematical model for calculating the PM purification amount Is used to determine the amount of time change in the PM purification amount, and the PM accumulation in the continuous regeneration type DPF is obtained from the time variation amount of the PM generation amount, the time variation amount of the PM purification amount, and the collection rate of the continuous regeneration type DPF. Estimate the amount.
[0018]
The coefficients of the mathematical model PM generation amount prediction model equation and the PM purification amount prediction model equation are obtained by measuring in advance the relationship among the engine speed, the fuel flow rate, the DPF differential pressure change rate, and the PM generation amount. Thus, it is determined by the least square method or the like.
[0019]
That is, the DPF differential pressure change rate during traveling is determined using a mathematical model polynomial in which the coefficient is determined by measuring the PM generation amount and the differential pressure corresponding to a plurality of engine operating conditions (engine speed, fuel flow rate) in advance. From this, the amount of time change of the PM generation amount and the PM purification amount is calculated, and further, the PM accumulation amount in the continuous regeneration type DPF is calculated from these.
[0020]
According to the exhaust gas purification system of the present invention, the PM accumulation amount can be accurately estimated from the engine speed, the fuel flow rate, and the DPF differential pressure change rate based on a mathematical model. Thus, the DPF forced regeneration processing control can be entered, and it is possible to prevent melting of the DPF, deterioration of fuel consumption of the engine, and the like.
[0022]
By providing this coefficient correction means, the PM generation rate and the PM purification rate can be accurately calculated with a mathematical model in which a coefficient is set according to the range of the DPF differential pressure change rate.
[0023]
Further, the mathematical model for calculating the PM generation rate by the PM generation rate calculating means is a quadratic expression relating to the rotational speed and the fuel flow rate, and the mathematical model for calculating the PM purification rate by the PM purification rate calculating means is the rotational speed. By using a quadratic expression related to the product of the fuel flow rate, the PM generation rate and the PM purification rate can be calculated with high accuracy by a relatively simple approximate expression.
[0024]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, an exhaust gas purification system according to an embodiment of the present invention will be described with reference to an example of an exhaust gas purification system provided with a continuous regeneration type DPF constituted by a combination of an oxidation catalyst (DOC) and a filter with catalyst (CSF). Will be described with reference to FIG.
[0025]
FIG. 1 shows a configuration of an exhaust
[0026]
This
[0027]
When a monolith honeycomb wall flow type filter is adopted as the filter with
[0028]
For the regeneration control of the filter with
[0029]
For controlling the
[0030]
The output values of these sensors are input to a control device (ECU: engine control unit) 30 that performs overall control of the operation of the
[0031]
The
[0032]
In the exhaust
[0033]
Then, the engine speed N and the fuel flow rate Q are input to the
[0034]
The
[0035]
The data input means C11 inputs the differential pressure ΔP before and after the continuous regeneration
[0036]
Further, the PM generation amount calculation means C15 determines the predetermined time t from the PM generation rate Rm1.kTo calculate a PM generation amount ΔM1 during a predetermined time interval Δt, and the PM purification amount calculation means C16 determines a predetermined time t from the PM purification rate Rm2.kThe PM purification amount ΔM2 during the predetermined time interval Δt is calculated, and the PM trapping amount integrating means C17 calculates the PM time from the PM generation amount ΔM1 and the PM purification amount ΔM2.kT after a predetermined time interval Δt fromk + 1The PM accumulation amount PMs is calculated.
[0037]
When the DPF differential pressure change rate Rdp exceeds a predetermined set value Rdp0, the coefficient correcting means C18 records the engine speed N, the fuel flow rate Q, the DPF differential pressure change rate Rdp, and this DPF differential pressure change rate. When the number of times Ldp exceeds the predetermined set value Rdp0 exceeds the predetermined number of times L0, the coefficient of the mathematical model of the PM generation rate calculation means C13 and the PM purification rate calculation means C14 is modified.
[0038]
Next, a mathematical model used in the PM generation rate calculation unit C13 and the PM purification rate calculation unit C14 and a mathematical model for estimating these coefficients from experimental data will be described.
[0039]
This mathematical model is expressed as follows: DPF differential pressure ΔP, PM generation amount M1, PM purification amount M2, PM trapping amount M3, time change amounts Rdp (= dΔP / dt), Rm1 (= dM1 / dt), Rm2 (= dM2) / Dt), Rm3 (= dM3 / dt), a PM generation amount prediction model equation for estimating Rm1, a PM purification amount prediction model equation for estimating Rm2, and coefficients of these model equations are determined from experimental data. It is a mathematical formula showing the relationship among Rdp, Rm1, Rm2, and Rm3 used for the purpose of
[0040]
First, regarding the PM generation amount prediction model formula, the PM generation rate Rm1 which is the time change rate of the PM generation amount M1 is used, the engine speed N and the fuel flow rate Q are parameters, and the coefficients a, b, c, d and e are set. It is expressed by equation (1), which is a quadratic polynomial that is a constant.
[0041]
Next, regarding the PM purification amount prediction model equation, the PM purification rate Rm2 that is the time change rate of the PM purification amount M2 is used as a parameter with the product (N * Q) of the engine speed N and the fuel flow rate Q as a parameter. This is expressed by equation (2), which is a quadratic polynomial with B and C as constants.
[0042]
[Expression 1]
[0043]
On the other hand, the time change Rdp of the DPF differential pressure ΔP under the engine operating conditions of a specific engine speed N and fuel flow rate Q is expressed by the equation (3), where K is a proportional constant. That is, here, α (0 <α <1) is the DPF collection rate, and the time change Rdp of the DPF differential pressure ΔP is the time change rate of the PM collection amount M3 (= M1 * α) in the DPF. It is assumed that it is proportional to the difference between the PM collection rate Rm3 (= Rm1 * α) and the PM purification rate Rm2.
[0044]
Here, since Rm2 = 0 under the operating condition of the engine in which PM purification does not occur, Equation (4) is obtained, and Equation (5) is obtained by transforming Equation (3).
[0045]
[Expression 2]
[0046]
Therefore, from this equation (4), the proportionality constant K can be determined from the test data of Rdp and Rm1 under the engine operating conditions (N, Q), and Rdp and Rm1 under five different engine operating conditions (N, Q) can be determined. If there is test data, the coefficients a, b, c, d, and e of the PM generation amount prediction model equation (1) can be determined. The DPF collection rate α is determined in advance by experiments and the like.
[0047]
In addition, since the PM purification rate Rm2 can be calculated from the equation (5), if there is data of Rdp and Rm1 for three or more different engine operating conditions (N, Q) where PM purification is estimated to occur, PM The coefficients A, B, and C of the purification amount prediction model equation (2) can be determined.
[0048]
The coefficients of the equations (1) and (2) used in the regeneration control means C1 of the
[0049]
First, the coefficient determination of the mathematical models (1) and (2) will be described.
[0050]
In the control flow for determining the coefficient of the PM generation amount prediction model equation (1) in FIG. 4, the coefficients a, b, c, d, and e of the equation (1), which is the PM generation amount prediction model equation, are calculated by the least square method. Establish the original simultaneous linear equation and solve it. In addition, since the solution method of simultaneous linear equations is well-known, it is not shown in FIG.
[0051]
When the control flow of FIG. 4 starts, in step S11, the engine speed Ni, the fuel flow rate Qi, PM generation rate Rm1i, I = 1 to n test data is input. However, n is the number of data and is an integer of 5 or more.
[0052]
In the next step S12, the coefficients, etc., of the five-way simultaneous linear equation (the expression of step S13) for obtaining the coefficients a, b, c, d, e, R, S, T, U, V, F, G, X, Y, H, I, J, K, M, W, and Z are calculated. However, the symbol “Σi"Indicates the sum of i = 1 to n.
[0053]
In step S13, the five-way simultaneous linear equations are solved, and coefficients a, b, c, d, and e are calculated. Thereby, since the PM generation amount prediction model equation (1) can be determined, this control flow is ended and the process returns.
[0054]
In the control flow for determining the coefficient of the PM purification amount prediction model equation (2) in FIG. 5, the ratio R12 of the PM purification rate and the PM generation rate from the equation (6) obtained by modifying the equation (5).i = Rm 2i / Rm 1i Is assumed to be an engine operating condition in which PM purification occurs, and for this engine operating condition, the PM purification rate Rm from the equation (5)2iAnd the coefficients A, B, and C of the formula (2) that is a PM purification amount prediction model formula are obtained by solving a three-dimensional simultaneous linear equation by the least square method. Of course, the solution of simultaneous linear equations is well known and is not shown in FIG.
[0055]
[Equation 3]
[0056]
When the control flow of FIG. 5 starts, in step S21, the engine speed Ni, Fuel flow Qi, PM generation rate Rm1i, I = 1 to n test data is input. However, n is the number of data and is an integer of 5 or more.
[0057]
In step S22, PM generation rate Rm1 for n data.iAnd DPF differential pressure change rate RdpiAnd the proportionality constant K previously determined from the test data, the ratio R12 of the PM purification rate and the PM generation rate according to the equation (6)iIs calculated.
[0058]
In step S23, the ratio R12 of the PM purification rate and the PM generation rateiData when engine is 1 or more, engine speed Nj, Fuel flow Qj, PM generation rate Rm1j, DPF differential pressure change rate Rdpj, J = 1 to m. However, m is the number of selected data, and is an integer of 3 or more and n or less.
[0059]
In step S24, the measurement data N selected in step S23.j, Qj, Rm1j, Rdpj, J = 1 to m, the PM purification rate Rm2j is calculated from the equation (5).
[0060]
In the next step S25, the engine speed Nj, Fuel flow QjAnd PM purification rate Rm2 calculated in step S24j, J = 1 to m, the coefficients of ternary simultaneous linear equations (formula of step S26) for obtaining the coefficients A, B, and C of the PM purification rate prediction model formula (2), etc. R, S, T, F, Calculate G, H, and J. However, the symbol “Σj"Indicates the sum of j = 1 to m.
[0061]
In step S26, the ternary simultaneous linear equations are solved and the coefficients A, B, and C are calculated. Thereby, since the PM purification amount prediction model equation (2) can be determined, this control flow is ended and the process returns.
[0062]
Next, regarding the control for calculating the PM accumulation amount based on the PM generation amount prediction model equation (1) and the PM purification amount prediction model equation (2) in which each coefficient is determined by the measurement data, the control flow of FIG. This will be explained based on.
[0063]
When this control flow is called from the regeneration control flow and starts, in step S31, the time tkAnd the time tkAnd record the time interval Δt. This time interval Δt is a time interval for this control, and the measurement data Ni, Qi, Rm1i, RdpiDifferent from the sampling interval.
[0064]
At time t32, time tkEngine speed Nk, Fuel flow Qk, DPF differential pressure ΔPkIn step S33, using these data, the time tkDPF differential pressure change rate Rdpk(= ΔPk/ Δt), and the PM purification rate Rm1 from the mathematical models (1) and (2)kAnd PM production rate Rm2kIs calculated.
[0065]
In the next step S34, the time tkTo time tk + 1(= TkPM generation amount ΔM1 within Δt time of + Δt)k(= Rm1k* Δt) and PM purification amount ΔM2k(= Rm2k* Δt) is calculated, and in step S35, the time tk + 1PM accumulation amount (PMk + 1= PMk+ ΔM1k* Α-ΔM2k) Is calculated. Here, α (0 <α <1) is the DPF collection rate.
[0066]
In step S36, the DPF differential pressure change rate RdpkWithin the set range (Rdp) where the coefficients of the equations (1) and (2) of the current mathematical model are determined.k≦ Rdp0), if it is within the set range, the process goes to step S38 as it is, and if it is outside the set range, the number of times L outside the set range is counted (L = L + 1), and the engine speed N as the measurement data at that timek, Fuel flow Qk, DPF differential pressure change rate RdpkAnd go to step S38.
[0067]
In step S38, it is determined whether or not the number of times L outside the set range has exceeded the predetermined set number of times L0.In step S39,Coefficient determination control flow of PM generation amount prediction model equation (1)(Fig. 4)And control flow for determining the coefficient of PM purification amount prediction model equation (2)(Fig. 5), And after correcting the coefficients of each mathematical model (1) and (2), the process returns.
[0068]
According to the control for calculating the PM accumulation amount described above, in this mathematical model, the time change rates Rm1, Rm2, and Rm3 of the PM generation amount M1, the PM purification amount M2, and the PM collection amount M3 in the continuous regeneration type DPF 13 (= Rm1 * α) is calculated by the equations (1) and (2) using the engine speed N and the fuel flow rate Q as variables.
[0069]
Therefore, based on the operation history information in which the engine speed N, the fuel flow rate Q, the DPF differential pressure change rate Rdp, and the elapsed time t between them are recorded, the mathematical models (1) and (2) are used. The amount of change ΔM1, PM purification amount ΔM2, and PM collection amount ΔM3 within a predetermined time Δt of the PM generation amount M1, the PM purification amount M2, and the PM collection amount M3 can be estimated.
[0070]
In addition, the initial value PMs of the PM accumulation amount (PM weight)1Is known, by accumulating the change amount ΔPMs (= ΔM3−ΔM2 = ΔM1 * α−ΔM2) of the PM accumulation amount PMs, the time tk + 1PM accumulation amount PMsK + 1Can be estimated.
[0071]
Therefore, according to the exhaust
[0072]
Next, a calculation example of the PM trapping amount and the PM purification amount calculated by the above control will be shown.
[0073]
In the test state as shown in FIG. 7, a measurement test of the PM generation amount M1 and the DPF differential pressure change rate Rdp was performed. In this measurement test, the engine speed N is constant (2000 rpm), and the fuel flow rate Q is increased at regular time intervals so that the time change of the DPF inlet temperature Tin rises stepwise. The DPF differential pressure ΔP was raised stepwise according to the exhaust gas flow rate and the PM generation amount M1. An example of measured values of the DPF inlet temperature Tin and the DPF differential pressure ΔP at this time is shown in FIG.
[0074]
In a time zone in which the DPF inlet temperature Tin is constant, if the temperature level is lower than the temperature at which the PM purification reaction occurs, PM is collected in the DPF, so the DPF differential pressure ΔP becomes constant or increases. On the other hand, if the temperature level is equal to or higher than the temperature at which the PM purification reaction occurs, the DPF differential pressure ΔP decreases, balances, or increases depending on the amount of collected PM and the amount of PM to be purified. The DPF differential pressure ΔP is either decreased, balanced, or increased.
[0075]
Then, the average value of the DPF differential pressure change rate Rdp in the vicinity of the middle point of the time zone (10 minutes) where the DPF inlet temperature Tin is constant is used as the coefficient determination measurement data Rdp of the prediction model (1).iIt was. The measurement data M1m and Rdpm of the PM generation amount M1 and the DPF differential pressure change rate Rdp obtained by this test are shown in a map display using the engine speed N and the fuel flow rate Q as parameters, as shown in FIGS.
[0076]
The proportional constant K is estimated to be 0.004 (mmAq / s) from these measured data M1m and Rdpm and the DPF differential pressure model (4), and the PM generation amount prediction model formula is calculated from the proportional constant K and the measured data M1m and Rdpm. The coefficient of (1) was determined. FIG. 11 shows a map display in which the calculated value M1c of the PM generation amount M1 calculated by the PM generation amount prediction model equation (1) uses the engine speed N and the fuel flow rate Q as parameters.
[0077]
Further, the estimated M2e of the PM purification amount M2 was calculated from the proportionality constant K, the measurement data M1, Rdp and the equation (5). Further, the coefficient of the PM purification amount prediction model equation (2) is determined, and the result obtained by supplementing the estimated value M2e with the calculated value M2c calculated by this PM purification amount prediction model equation (2) is the engine speed N, fuel flow rate. A map display with Q as a parameter is shown in FIG.
[0078]
Further, FIG. 13 shows a ratio R12 between the PM purification rate and the PM generation rate in a map display using the engine speed N and the fuel flow rate Q as parameters. If this ratio R12 is greater than 1, it can be estimated that the PM purification rate Rm2 is greater than the PM generation rate Rm1, so select measurement data that will cause the ratio R12 calculated by equation (6) to be greater than 1. The coefficient of the PM purification amount prediction model equation (2) is determined from the selected measurement data. A region where the ratio R12 is larger than 1 is indicated by hatching.
[0079]
【The invention's effect】
As described above, according to the exhaust gas purification system of the present invention, based on the mathematical model, from the engine speed, the fuel flow rate, and the DPF differential pressure change rate that is the time change rate of the differential pressure across the DPF, Since the PM accumulated amount can be accurately estimated, it is possible to accurately detect whether the amount of PM accumulated in the continuous regeneration type DPF has reached the PM accumulation amount at the start of regeneration, and regeneration of the DPF at an appropriate time. Since the process control can be entered, it is possible to prevent the DPF from being melted and the fuel consumption from being deteriorated.
[0080]
In addition, when the number of times that the DPF differential pressure change rate exceeds the predetermined set value exceeds the predetermined set number, by providing coefficient correction means for correcting the coefficient of the mathematical model, it is possible to respond to the range of the differential pressure change rate. The PM generation rate and the PM purification rate can be calculated accurately with a mathematical model in which coefficients are set.
[Brief description of the drawings]
FIG. 1 is a system configuration diagram of an exhaust gas purification system according to an embodiment of the present invention.
FIG. 2 is a diagram showing a configuration of an engine, a continuous regeneration type DPF, a measurement system, and a data processing system according to an embodiment of the present invention.
FIG. 3 is a diagram showing a configuration of control means of the exhaust gas purification system according to the present invention.
FIG. 4 is a diagram illustrating an example of a control flow for determining a coefficient of a PM generation amount prediction model equation;
FIG. 5 is a diagram illustrating an example of a control flow for determining a coefficient of a PM purification amount prediction model equation;
FIG. 6 is a diagram illustrating an example of a control flow for calculating a PM accumulation amount.
FIG. 7 is a diagram showing a configuration of an engine, a continuous regeneration type DPF, a measurement system, and a data processing system during a measurement test.
FIG. 8 is a diagram showing an example of actual measurement of DPF inlet temperature and DPF differential pressure.
FIG. 9 is a diagram showing a measured value of the PM generation amount in a map display.
FIG. 10 is a diagram showing a measured value of a DPF differential pressure change rate in a map display.
FIG. 11 is a diagram showing a calculated value of the PM generation amount in a map display.
FIG. 12 is a diagram showing a calculated value of the PM purification amount in a map display.
FIG. 13 is a diagram showing the ratio between the PM purification amount and the PM generation amount in a map display.
[Explanation of symbols]
1 Exhaust gas purification system
10 Diesel engine (internal combustion engine)
12 Exhaust passage
13 Continuous regeneration type DPF device
13a Oxidation catalyst
13b Filter with catalyst (DPF)
21 Differential pressure sensor
30 control unit (ECU)
30a Reproduction control device
C1 Reproduction control means
C10 reproduction start judging means
C11 data input means
C12 Differential pressure change rate calculation means
C13 PM generation rate calculation means
C14 PM purification rate calculation means
C15 PM generation amount calculation means
C16 PM purification amount calculation means
C17 PM accumulated amount calculation means
C18 coefficient correction means
C20 Normal control operation means
C30 regeneration processing control operation means
L times
L0 Preset number of times
N engine speed
PMs PM accumulation amount
Q Fuel flow rate
Rdp DPF differential pressure change rate
Rm1 PM production rate
Rm2 PM purification rate
Rdp0 Preset value
tk Predetermined time
ΔP differential pressure
Δt Predetermined time interval
ΔM1 PM production
ΔM2 PM purification amount
Claims (3)
該再生制御手段は、
エンジン回転数及び燃料流量におけるDPF差圧率及びPM生成率の試験データを複数記憶する記憶手段と、
前記複数の試験データの中から選択された試験データを用いて数学モデルの係数の決定を行う手段と、
前記差圧とエンジン回転数とエンジンの燃料流量を入力するデータ入力手段と、
前記差圧の時間に関する変化率であるDPF差圧変化率を算定する差圧変化率算定手段と、
前記回転数と前記燃料流量から数学モデルに基づいてPM生成率を算定するPM生成率算定手段と、
前記回転数と前記燃料流量から数学モデルに基づいてPM浄化率を算定するPM浄化率算定手段と、
前記PM生成率から所定時間の間のPM生成量を算出するPM生成量算定手段と、
前記PM浄化率から所定時間の間のPM浄化量を算出するPM浄化量算定手段と、
前記PM生成量と前記PM浄化量とから所定時間後のPM蓄積量を算定するPM蓄積量算定手段を備え、
前記再生制御装置が、前記DPF差圧変化率が所定の設定値を超える場合は、エンジン回転数、燃料流量、DPF差圧変化率を記録し、このDPF差圧変化率が所定の設定値を超えた回数が前記所定の設定回数を超える場合は、このDPF差圧変化率に応じた範囲の試験データを用いて、前記数学モデルの係数を決定することにより、前記数学モデルの係数を修正する係数修正手段を備えることを特徴とする排気ガス浄化システム。A continuous regeneration type DPF that collects particulates provided in an exhaust passage of a diesel engine, a differential pressure sensor that detects a differential pressure before and after the continuous regeneration type DPF, and a regeneration control device that performs regeneration control of the continuous regeneration type DPF In an exhaust gas purification system equipped with
The reproduction control means includes
Storage means for storing a plurality of DPF differential pressure rate and PM generation rate test data at engine speed and fuel flow rate;
Means for determining coefficients of a mathematical model using test data selected from the plurality of test data;
Data input means for inputting the differential pressure, the engine speed and the fuel flow rate of the engine;
A differential pressure change rate calculating means for calculating a DPF differential pressure change rate which is a change rate with respect to time of the differential pressure;
PM generation rate calculation means for calculating a PM generation rate based on a mathematical model from the rotational speed and the fuel flow rate;
And PM purification ratio calculating means for calculating a PM purification ratio on the basis of the mathematical model from the fuel flow rate and the rotational speed,
PM generation amount calculation means for calculating a PM generation amount during a predetermined time from the PM generation rate;
PM purification amount calculation means for calculating a PM purification amount for a predetermined time from the PM purification rate;
PM accumulation amount calculation means for calculating the PM accumulation amount after a predetermined time from the PM generation amount and the PM purification amount,
If the DPF differential pressure change rate exceeds a predetermined set value, the regeneration control device records the engine speed, the fuel flow rate, the DPF differential pressure change rate, and the DPF differential pressure change rate is set to a predetermined set value. If the number of times exceeds the predetermined number of times, the coefficient of the mathematical model is corrected by determining the coefficient of the mathematical model using test data in a range corresponding to the DPF differential pressure change rate. An exhaust gas purification system comprising coefficient correction means.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003086921A JP4333180B2 (en) | 2003-03-27 | 2003-03-27 | Exhaust gas purification system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003086921A JP4333180B2 (en) | 2003-03-27 | 2003-03-27 | Exhaust gas purification system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004293413A JP2004293413A (en) | 2004-10-21 |
JP4333180B2 true JP4333180B2 (en) | 2009-09-16 |
Family
ID=33401410
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003086921A Expired - Fee Related JP4333180B2 (en) | 2003-03-27 | 2003-03-27 | Exhaust gas purification system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4333180B2 (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1788210A3 (en) * | 2004-01-13 | 2012-01-11 | EMCON Technologies LLC | Method and apparatus for directing exhaust gas through a fuel-fired burner of an emission abatement assembly |
JP5325249B2 (en) | 2011-03-18 | 2013-10-23 | 株式会社小松製作所 | Particulate matter accumulation amount estimation device, exhaust gas purification system, and particulate matter accumulation amount estimation method |
KR101317410B1 (en) | 2011-11-22 | 2013-10-10 | 서울대학교산학협력단 | Nox mass prediction method |
CN104411928B (en) * | 2012-11-13 | 2017-03-08 | 丰田自动车株式会社 | The emission-control equipment of internal combustion engine |
WO2017179277A1 (en) * | 2016-04-11 | 2017-10-19 | ボッシュ株式会社 | Filter reproduction system, ecu and filter reproduction program |
CN107002528B (en) | 2016-12-14 | 2019-10-25 | 株式会社小松制作所 | Waste gas cleaning system and exhaust gas purifying method |
JP7113598B2 (en) | 2017-07-14 | 2022-08-05 | 株式会社Soken | Failure detection device and failure detection method for particulate filter |
JP6780763B2 (en) * | 2018-12-25 | 2020-11-04 | トヨタ自動車株式会社 | Internal combustion engine control device |
CN112983608A (en) * | 2021-02-09 | 2021-06-18 | 联合汽车电子有限公司 | Particle trap carbon load calculation method and system, controller and storage medium |
CN113719366B (en) * | 2021-09-22 | 2022-08-23 | 潍柴动力股份有限公司 | DPF parking regeneration control method and device for vehicle |
CN115779589B (en) * | 2023-02-06 | 2023-05-05 | 山东国瀚机电科技有限公司 | Intelligent control method for electric control system of mining dry dust collector |
-
2003
- 2003-03-27 JP JP2003086921A patent/JP4333180B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2004293413A (en) | 2004-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4042476B2 (en) | Exhaust gas purification device for internal combustion engine | |
US7797927B2 (en) | Method for control of exhaust gas purification system, and exhaust gas purification system | |
JP4513593B2 (en) | Exhaust gas purification device for internal combustion engine | |
US8181449B2 (en) | Control method of exhaust gas purification system and exhaust gas purification system | |
JP4389606B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP4403961B2 (en) | Exhaust gas purification device for internal combustion engine | |
US7628010B2 (en) | Exhaust purification system | |
US20070209333A1 (en) | Particulate matter accumulation amount detection apparatus and method | |
US7210285B2 (en) | Internal combustion engine exhaust gas purifying system | |
JP2004340138A (en) | System and method for pressure sensor diagnosis by computer | |
US10273857B2 (en) | Method and apparatus for monitoring a particulate filter | |
JP6213118B2 (en) | Diagnostic equipment | |
JP4333180B2 (en) | Exhaust gas purification system | |
JP4506060B2 (en) | Particulate filter regeneration control device | |
JP4308702B2 (en) | Exhaust purification equipment | |
JP5206870B2 (en) | Exhaust gas purification device for internal combustion engine | |
US8161737B2 (en) | Apparatus for and method of purifying exhaust gas | |
JP2004132358A (en) | Filter control device | |
JP4008867B2 (en) | Exhaust purification equipment | |
JP4424071B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP5912494B2 (en) | Diesel engine exhaust purification system | |
JP2006063970A (en) | Exhaust emission control device of internal combustion engine | |
JP5493268B2 (en) | Exhaust gas purification system control method and exhaust gas purification system | |
JP7384114B2 (en) | Filter condition detection device | |
JP2012036859A (en) | Exhaust emission control device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051215 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081118 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090114 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090303 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090428 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090602 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090615 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120703 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120703 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130703 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |