JP2004293338A

JP2004293338A - Method for presuming amount of nox occlusion

Info

Publication number: JP2004293338A
Application number: JP2003083650A
Authority: JP
Inventors: Yoshiki Tanabe; 圭樹田邊; Susumu Koketsu; 晋纐纈; Shinji Nakayama; 真治中山; Minehiro Murata; 峰啓村田; Daisuke Haruhara; 大輔春原
Original assignee: Mitsubishi Fuso Truck and Bus Corp
Current assignee: Mitsubishi Fuso Truck and Bus Corp
Priority date: 2003-03-25
Filing date: 2003-03-25
Publication date: 2004-10-21
Also published as: KR100592414B1; US20040244361A1; DE102004013604B4; DE102004013604A1; KR20040084689A

Abstract

<P>PROBLEM TO BE SOLVED: To presume an amount of NOx occlusion at high accuracy in a method for presuming an amount of NOx occlusion occluded by a NOx occlusion type reducing catalyst interposed in an exhaust passage. <P>SOLUTION: In a method for presuming an amount of NOx occlusion of a NOx occlusion catalyst interposed in an exhaust passage 2 of an engine, the amount of NOx occlusion is presumed using a polynominal reflecting an NOx occlusion characteristic of the NOx occlusion catalyst 1, and a coefficient of the polynominal is sequentially corrected on the basis of a NOx purifying rate actually measured. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

【０００１】
【発明の属する技術分野】
本発明は、排気通路に設けられたＮＯｘ吸蔵型還元触媒に吸蔵されたＮＯｘの吸蔵量推定方法に関する。
【０００２】
【従来の技術】
一般に、ディーゼルエンジンでは、その構造上の特性から排ガス中の酸素が多く酸化雰囲気となる（つまり、排気空燃比がリーンとなる）。また、ガソリンエンジンであっても希薄燃焼可能なものであれば、希薄燃焼時にはやはり排気空燃比がリーンとなる。そこで、このような酸化雰囲気でもＮＯｘを浄化できるよう、従来よりリーンＮＯｘ触媒が開発されている。
【０００３】
このようなリーンＮＯｘ触媒としては、ＮＯｘ吸蔵型還元触媒（以下、ＮＯｘ吸蔵触媒、又は単に触媒という）がある。ＮＯｘ吸蔵触媒では、酸化雰囲気では排気中のＮＯｘが吸蔵され、低酸素濃度下においてＣＯが存在すると吸蔵されたＮＯｘが放出されるようになっており、ＮＯｘ吸蔵触媒から放出されたＮＯｘの殆どが、エンジンから排出された未燃ＨＣ，ＣＯ等の還元剤により無害なＮ_２へと還元されてから大気へと排出される。
【０００４】
そして、ＮＯｘ吸蔵触媒が装備された車両では一般的にいわゆるリッチスパイク制御が行なわれる。つまり、リーンな空燃比での運転が続くとＮＯｘ吸蔵触媒はこれ以上ＮＯｘを吸蔵できない飽和状態となるため、タイミングを見計らって空燃比を一時的に強制的にリッチにする（リッチスパイク制御を行なう）ことで排気中の酸素濃度を低下させるとともに、還元剤を供給することでＮＯｘ吸蔵触媒からＮＯｘを放出させ、これにより、ＮＯｘを吸蔵しうる状態にＮＯｘ吸蔵触媒を再生するのである。
【０００５】
ところで、上述のようなリッチスパイク制御を行なうには、ＮＯｘ吸蔵触媒におけるＮＯｘ吸蔵量を正確に推定（又は検出）する必要があり、このようなＮＯｘ吸蔵量の推定手法としては、例えば下記のような技術が知られている。
▲１▼第１の従来技術：ＮＯｘセンサを用いたＮＯｘ吸蔵量推定手法
この第１の従来技術では、図６に示すように、触媒１０１の上流と下流とにそれぞれＮＯｘの濃度を検出するＮＯｘセンサ１０２，１０３を設けるとともに、排気通路中に排ガス流量を検出するエアフローセンサ（図示省略）を設ける。
【０００６】
そして、ＥＣＵ１０４では、ＮＯｘセンサ１０２，１０３からの情報及びエアフローセンサからの情報に基づいて、下式により、ＮＯｘ吸蔵量を推定する。
ＮＯｘ吸蔵量＝∫排ガス流量×（触媒上流ＮＯｘ濃度−触媒下流ＮＯｘ濃度）
▲２▼第２の従来技術：数学的モデルによる推定手法（例えば、特許文献１）
この第２の従来技術では、吸蔵反応，酸化還元反応及び離脱反応等の触媒の化学現象及び物理現象に基づく数学的な触媒モデルを設け、この触媒モデル式によりＮＯｘ吸蔵量を推定する。
【０００７】
【特許文献１】
特開平９−７２２３５号公報
【０００８】
【発明が解決しようとする課題】
しかしながら、上記第１の従来技術においては、リーン運転中はＮＯｘ吸蔵量を算出することは可能だが、リッチ運転中のＮＯｘ放出量（吸蔵量の減少量）はＮＯｘセンサでは検出することができない。このためリッチ運転時間が不十分だった場合には、触媒に残留したＮＯｘが累積してＮＯｘ吸蔵量の演算誤差が生じるという課題がある。また、触媒の上流と下流とにそれぞれＮＯｘセンサが必要であり、コスト増を招くという課題がある。
【０００９】
また、第２の従来技術では、少なくとも触媒下流に１つのＮＯｘセンサを設けるだけでＮＯｘ吸蔵量を推定可能であるが、この技術は三元触媒の空燃比制御に適用することを主眼としているため、モデル構造がＮＯｘ吸蔵触媒とは異なり、ＮＯｘ吸蔵触媒への応用が困難であるという課題がある。また、その他の触媒モデルによりＮＯｘ吸蔵量を推定する場合においても触媒種毎に特性値が異なるためモデル式を変更する必要があり、加えて触媒の劣化などに対応するためには、より詳細なモデル化を必要とするという課題がある。
【００１０】
本発明は、このような課題に鑑み創案されたもので、触媒のモデルに常に最新の状態を反映させて高い精度でＮＯｘ吸蔵触媒のＮＯｘ吸蔵量を推定できるようにした、ＮＯｘ吸蔵量の推定方法を提供することを目的とする。
【００１１】
【課題を解決するための手段】
このため、請求項１記載の本発明のＮＯｘ吸蔵量の推定方法は、エンジンの排気通路に介装されたＮＯｘ吸蔵触媒のＮＯｘ吸蔵量を推定する方法であって、該ＮＯｘ吸蔵量を、該ＮＯｘ吸蔵触媒のＮＯｘ吸蔵特性を反映させた多項式を用いて推定するとともに、該各多項式の係数を、実測されたＮＯｘ浄化率に基づいて逐次補正することを特徴としている。
【００１２】
また、該ＮＯｘ吸蔵触媒におけるＮＯｘ放出量を下式により算出してもよい。ＮＯｘ放出量＝∫（触媒入口還元剤濃度×還元剤利用率
−０．５×触媒入口酸素濃度）×排ガス流量
なお、好ましくは、該還元剤利用率をＮＯｘ吸蔵触媒のＮＯｘ放出特性を反映させた多項式を用いて推定するとともに、該多項式の係数を還元剤の濃度に基づいて逐次補正する。
【００１３】
また、好ましくは、該エンジンが、排気空燃比がリーンとなるリーン運転と、該排気空燃比がリッチとなるリッチ運転とを切り替え可能に構成され、該リッチ運転中は該多項式の係数を保持するとともに、リーン運転開始時に該保持された係数を用いて得られるＮＯｘ浄化率と該実測されたＮＯｘ浄化率との差がしきい値以上であると該ＮＯｘ吸蔵量を補正する。
【００１４】
さらには、該各係数の平均値がいずれか一つでもそれぞれ設定された所定の範囲を逸脱すると、該触媒が異常であると判定するようにしてもよい。
【００１５】
【発明の実施の形態】
以下、図面により、本発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法について説明すると、図１は本発明が適用される車両の排気系の全体構成を示す模式図である。
図示するように、排気通路２にはＮＯｘ吸蔵触媒１が介装されており、このＮＯｘ吸蔵触媒１の下流側には触媒出口のＮＯｘ濃度を検出するＮＯｘセンサ３が設けられている。
【００１６】
また、図示はしないが、排気通路２には、Ｏ_２センサ，温度センサ及びエアフローセンサ（ＡＦＳ）等が接続されており、これらのセンサにより排気通路２内の酸素濃度，排ガス温度及び排ガス流量がそれぞれ検出されるようになっている。
また、各センサ３は制御手段としてのＥＣＵ４に接続されている。ここで、ＥＣＵ４は入出力装置，記憶装置（ＲＯＭ，ＲＡＭ，不揮発性ＲＡＭ等），演算装置（ＣＰＵ），タイマカウンタ等を備えて構成されており、このＥＣＵ４により、ＮＯｘ吸蔵触媒１におけるＮＯｘ吸蔵量の推定が行なわれるようになっている。
【００１７】
また、ＥＣＵ４には、種々のマップが設けられており、例えばアクセル開度やエンジン回転数をパラメータとして触媒入口のＮＯｘ，ＣＯ，ＨＣの各濃度（以下、それぞれ触媒入口［ＮＯｘ］，触媒入口［ＣＯ］，触媒入口［ＨＣ］で表す）がマップから読み出されるようになっている。なお、本実施形態では、これらの値をマップから求めているが、排気通路２に種々のセンサを設け、これらの値をセンサから直接検出してもよい。
【００１８】
ところで、上述した従来技術２ではＮＯｘ吸蔵触媒の特性そのものを数式モデルで記述してＮＯｘ吸蔵量を求めていたのに対して、本発明では、ＮＯｘ吸蔵量を４次元の線形多項式を用いて算出するようになっている。
すなわち、ＥＣＵ４には下式（１）に示すような多項式が記憶されており、この多項式に基づいてＮＯｘ浄化率及びＮＯｘ吸蔵触媒１におけＮＯｘ吸蔵割合が算出されるようになっている。

なお、式（１）において、ｒはＮＯｘ浄化率、ｘはＮＯｘ吸蔵割合、ｙは排ガス温度、ｚはＳＶ値（又は排ガス流速）、ｋ_０，ｋ_１，ｋ_２，・・・は係数である。
【００１９】
ここで、ＮＯｘ吸蔵触媒１のＮＯｘ吸蔵特性は、実験等により予め知ることができ、多項式で近似することができる。例えば、パラメータが３つ（３次元）であれば、触媒の特性は図２に示すような曲面で表すことができ、これを多項式で表すことができるのである。もちろん、式（１）は４次元であるので図２に示すような曲面とはならないが、考え方は３次元の場合と同様である。
【００２０】
また、各係数ｋ_ｉ（ｉ＝１，２，・・・）は、予め実験等で求めた特性から適当な値が初期値として入力されている。
そして、エンジン運転中に、触媒下流のＮＯｘセンサ３の検出値に基づいて逐次各係数ｋ_ｉを補正することにより常に、式（１）により、最新の触媒状態を表すことができるのである。
【００２１】
以下、ＮＯｘ推定の手法についてさらに詳しく説明する。まず、式（１）におけるＮＯｘ浄化率ｒは下式（２）でも求めることができる。
ｒ＝触媒出口ＮＯｘ濃度／触媒入口ＮＯｘ濃度・・・・・・・・・・・・（２）
ここで、触媒出口ＮＯｘ濃度はＮＯｘセンサ３で検出された値であり、触媒入口ＮＯｘ濃度はマップから得られる値である。したがって、式（２）で得られるＮＯｘ浄化率ｒはセンサ値（実測値）ということができ、式（１）で得られるＮＯｘ浄化率ｒは推定値（計算値）ということができる。
【００２２】
また、式（１）の排ガス温度［ｙ］及びＳＶ値［ｚ］は実測データ又はマップ値が適用される。これにより、式（１）における未知数はＮＯｘ吸蔵割合［ｘ］のみとなり、これ以外の値は既知となる。したがって、式（１）を変形した下式（３）により、ＮＯｘ吸蔵割合［ｘ］を求めることができる。
ｘ＝〔ｒ−（ｋ_０＋ｋ_２ｙ＋ｋ_３ｚ・・）〕／（ｋ_１＋ｋ_４ｙ＋・・）・・（３）
そして、この式（３）で求めたＮＯｘ吸蔵割合［ｘ］を再び式（１）に戻して次回の計算周期での計算に用いるようになっている。つまり、式（１）では前回求めたＮＯｘ吸蔵割合［ｘ］と、今回の計算周期で新たに検出された［ｙ］及び［ｚ］とを用いて、浄化率ｒを計算するようになっているのである。
【００２３】
ここで、式（１）で算出される浄化率ｒは推定値であるが、各係数ｋ_ｉが正確な値であれば、この浄化率ｒの推定値と式（２）で得られる実測値とが一致するはずである。
そこで、本発明では、浄化率ｒの推定値と実測値とを比較して、推定値と実測値との間に差があれば、推定値が実測値となるように各係数ｋ_ｉを最小二乗法を用いて逐次補正するようになっている。そして、このような計算を繰り返し行ない各係数ｋ_ｉを逐次補正することで、式（１）の多項式を、ＮＯｘ吸蔵触媒１の状態を正確に反映した式に随時変更することができるのである。
【００２４】
つまり、ＮＯｘ吸蔵触媒１は使用状況や経年劣化等によりその特性が変化するものであるが、何度も係数ｋ_ｉを繰り返し補正することにより、式（１）が正確に触媒１の状態を表す式となり、高い精度でＮＯｘ吸蔵量を推定することができるようになるのである。なお、式（３）で求めたＮＯｘ吸蔵割合［ｘ］は、初期はその計算周期における瞬間的なＮＯｘ吸蔵量であるが、この計算結果を繰り返し式（１）に戻して浄化率ｒを算出することで、触媒１に吸蔵されたＮＯｘの積算値を算出することができる。
【００２５】
このように、ＮＯｘを吸蔵するリーン運転時においては、触媒下流に設けられたＮＯｘセンサ３から触媒出口［ＮＯｘ］を取り込むとともに、式（２）による推定値とを比較し、その都度、最小二乗法により係数ｋ_ｉを補正することで正確なＮＯｘ吸蔵量を得ることができる。
ところで、リッチ運転中はＮＯｘが放出されるため、正確なＮＯｘ放出量を検出又は推定できないと、次回のリーン運転時におけるＮＯｘ吸蔵量の初期値を正確に算出できなくなる。
【００２６】
そこで、リッチ運転中はｋ_ｉを保持するとともに、下式（４），（５）によりＮＯｘ放出量を算出し、リッチ運転時開始直前のＮＯｘ吸蔵量から式（４），（５）により算出されたＮＯｘ放出量を差し引いてリッチ運転中のＮＯｘ吸蔵量を演算するようになっている。
ＮＯｘ放出量＝∫（触媒入口還元剤濃度×還元剤利用率［ｒ′］−０．５×触媒入口［Ｏ_２］）×排ガス流量・・・・・（４）
還元剤利用率［ｒ′］＝ｆ（ｙ，ｚ）・・・・・・・（５）
ここで、ＥＣＵ１には、図３に示すような還元剤利用率設定マップが設けられている。このマップには、排ガス温度［ｙ］及びＳＶ値［ｚ］に基づく還元剤利用率の特性が記憶されており、これらのパラメータ［ｙ］，［ｚ］に基づいて還元剤利用率が設定されるようになっている。
【００２７】
また、触媒入口還元剤濃度（触媒入口［ＣＯ］，触媒入口［ＨＣ］）は、上述したようにアクセル開度やエンジン回転数をパラメータとしてマップから求めることができ、他のパラメータも同様にマップ値又は検出値として得ることができる。
したがって、リーン運転からリッチ運転に切り替えられると、ＥＣＵ４ではリーン運転終了時の各係数ｋ_ｉを保持するとともにＮＯｘ吸蔵量を記憶しておき、再びリッチ運転からリーン運転に切り替えられると、上記ＮＯｘ吸蔵量からリッチ運転終了時のＮＯｘ放出量を引き算することによりＮＯｘの残存量が算出される。
【００２８】
また、ＥＣＵ４では、リーン運転開始時のＮＯｘ浄化率ｒの実測値と推定値との差がしきい値以上であると、ＮＯｘ吸蔵量を補正するようになっている。つまり、リッチ運転からリーン運転への切り替え直後において、式（１）で推定されるＮＯｘ浄化率と式（２）で得られるＮＯｘ浄化率との差が所定値以上であると、リーン運転開始時に算出されたＮＯｘ残存量（＝前回のリーン運転終了時のＮＯｘ吸蔵量−リッチ運転終了時のＮＯｘ放出量）が正確ではないと判断して、このＮＯｘ残存量［ｘ］（つまり、リーン運転開始時のＮＯｘ吸蔵量）を、式（３）に基づいて補正するようになっている。
【００２９】
さらに、ＥＣＵ４は突発的な触媒１の異常等により係数ｋ_ｉが正常ではない値になると、これを検出してドライバに知らせる機能も有している。具体的には、各係数ｋ_ｉについて、所定の連続した計算周期における移動平均を常時算出し、この移動平均の値が各係数毎に設定された所定の範囲から逸脱すると、触媒１が劣化又は破損した等の異常が生じたと判定するようになっている。なお、所定の範囲は、例えば各係数ｋ_ｉの初期値に±α（αは一定値）を加算して設定されるようになっている。
【００３０】
また、ＥＣＵ４により触媒１の異常が判定されると、例えばインストルメントパネルの警告灯を点灯させるとともにリッチ運転が禁止されるようなっている。そして、このように触媒１の異常時にリッチ運転を禁止することにより、ＣＯを排出するような事態を未然に回避することができる。
本発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法は、上述のように構成されているので、以下のようにしてＮＯｘ吸蔵量が推定される。
【００３１】
まず、リーン運転時において、式（２）から実ＮＯｘ浄化率ｒを算出するとともに、実測データ又はマップ値から排ガス温度［ｙ］及びＳＶ値［ｚ］を求める。また、式（１）の多項式において、各係数ｋ_ｉに実験等により求めた適当な値を初期値として入力する。そして、式（１）を変形した式（３）に各値を代入し、ＮＯｘ吸蔵割合［ｘ］を算出する。
【００３２】
次に、式（３）で得られたＮＯｘ吸蔵割合［ｘ］と、新たに求めた排ガス温度［ｙ］及びＳＶ値［ｚ］とを式（１）に代入してＮＯｘ浄化率ｒ（推定値）を求める。そして、この推定ＮＯｘ浄化率ｒと、式（２）で新たに算出した実ＮＯｘ浄化率ｒとを比較し、推定ＮＯｘ浄化率が、実ＮＯｘ浄化率と一致するように各係数ｋ_ｉを最小二乗法により補正する。
【００３３】
そして、このような計算を繰り返し実行し、各係数ｋ_ｉを随時更新することで式（１）が正確に触媒１の状態を表す式となり、高い精度でＮＯｘ吸蔵量を推定することができるようになる。また、式（３）を用いて繰り返しＮＯｘ吸蔵割合［ｘ］を算出することにより、触媒１に吸蔵されたＮＯｘの積算値を算出することができる。
【００３４】
また、リッチ運転中はＮＯｘが放出されるので、このときには式（１）によるＮＯｘ吸蔵量の算出を中断して各係数ｋ_ｉを記憶しておくとともに、上述した式（４），（５）によりＮＯｘ放出量を算出する。
そして、再びリッチ運転からリーン運転に切り替えられると、リーン運転終了時のＮＯｘ吸蔵量からリッチ運転終了時のＮＯｘ放出量を減じることにより、リーン運転開始時のＮＯｘ残存量が算出される。
【００３５】
また、リーン運転開始時のＮＯｘ浄化率ｒの実測値と推定値との差がしきい値以上であると、ＮＯｘ吸蔵量が補正される。ここで、ＮＯｘ吸蔵量の補正手法について図４のフローチャートに基づいて説明すると、まず、リーン運転時に多項式（１）の各係数ｋ_ｉを逐次補正して各係数ｋ_ｉを同定する（ステップＳ１）。
そして、リッチ運転が開始したと判定されると（ステップＳ２）、各係数ｋ_ｉを保持しておき（ステップＳ３）、その後、リッチ運転が終了すると（ステップＳ４）、ＮＯｘ浄化率の実測値と推定値とを比較し、これらの差がしきい値以上であるかが判定される（ステップＳ５）。
【００３６】
ＮＯｘ浄化率の実測値と推定値との差がしきい値未満であれば、再びステップ１に戻ってリーン運転時の各係数ｋ_ｉの更新を続行し、しきい値以上であればＮＯｘ吸蔵量を式（３）で算出される値に補正して（ステップＳ６）リターンする。
そして、このようにしてＮＯｘ吸蔵量を補正することにより、ＮＯｘ吸蔵量の推定精度をさらに高めることができる。
【００３７】
一方、触媒１の劣化や異常が判定されると、インストルメントパネルの警告灯を点灯させて、ドライバに知らせるとともに、リッチ運転が禁止される。これを図５のフローチャートに基づいて説明すると、まず、多項式（１）の各係数ｋ_ｉを逐次補正して各係数ｋ_ｉを同定し（ステップＳ１１）、次に各係数ｋ_ｉの一定期間の移動平均を算出する（ステップＳ１２）。
【００３８】
次に、算出された各係数ｋ_ｉの移動平均がそれぞれ所定の範囲内か否かを判定し（ステップＳ１３）、各係数ｋ_ｉの移動平均がいずれも所定の範囲内であれば、再びステップＳ２に戻り、各係数ｋ_ｉの移動平均のうちいずれか一つでも所定の範囲外であると、触媒１に異常が生じたと判定して警報を行なう（ステップＳ１４）とともに、リッチ運転を禁止するのである（ステップＳ１５）。
以上詳述したように、本願発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法によれば、ＮＯｘ吸蔵量を、ＮＯｘ吸蔵触媒１のＮＯｘ吸蔵特性を反映させた多項式を用いて推定するとともに、多項式の係数ｋ_ｉを、実測されたＮＯｘ浄化率に基づいて逐次補正するので、常に最新の触媒状態をモデル化した多項式とすることができ、高い精度でＮＯｘ吸蔵量を推定することができるという利点がある。また、ＮＯｘ吸蔵量を正確に推定することでＮＯｘ吸蔵量に基づいた最適なリッチ運転制御を実行することができ、燃費を向上させることができる。
【００３９】
また、ＮＯｘ吸蔵触媒１を変更した場合でも、変更後の触媒１の特性に応じた各係数ｋ_ｉに更新されるので、モデル式を変更する必要がないという利点がある。これに対して、従来のように数学的なモデル式を適用した場合には、各触媒の特性に応じて係数を特定する作業（合わせ込み）が必要となり多大な労力を要するという課題があるが、本願発明ではこのような合わせ込みを行なうことなく、正確にＮＯｘ吸蔵量を推定することができる。
【００４０】
ＮＯｘ放出量を還元剤利用率に基づいて算出するので、比較的高い精度でＮＯｘ放出量を推定することができる。
また、リッチ運転中は多項式の係数ｋ_ｉを保持するとともに、リーン運転開始時に保持された係数を用いて推定されるＮＯｘ浄化率と、実測されたＮＯｘ浄化率との差がしきい値以上であるとＮＯｘ吸蔵量を補正するので、やはりＮＯｘ吸蔵量の推定精度を高めることができる。
【００４１】
各係数ｋ_ｉの所定期間における平均値が所定の範囲を逸脱すると、触媒１が異常であると判定するので、高い精度で触媒１の異常を判定できる。また、このような異常判定時にはドライバに警報を発するので、速やかにドライバは触媒１の異常を認識でき、早期にドライバに触媒１の交換を促すことができる。また、このような触媒１の異常時にはリッチ運転が禁止されるのでＣＯがそのまま排出されるような事態を回避することができる。
【００４２】
なお、本発明のＮＯｘ吸蔵量の推定方法は、上述の実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々の変形が可能である。例えば、ＮＯｘの放出量については、ＮＯｘ吸蔵量の推定手法と同様の手法により求めてもよい。
つまり、上述ではＮＯｘ放出量を計算する際に用いる還元剤利用率［ｒ′］を図３に示す還元剤利用率設定マップから求めているが、下式（６）に示すように、触媒入口還元剤濃度［ｘ′］，排ガス温度［ｙ］，ＳＶ値［ｚ］の多項式から求めてもよい。

なお、ｍ_ｉ（ｉ＝１，２，・・・）は係数である。また、式（６）はＮＯｘ吸蔵触媒１のＮＯｘ放出特性を反映させた多項式であって、各係数ｍ_ｉには、式（１）と同様に実験等から得られた適当な値が初期値として入力されている。
【００４３】
また、この場合には少なくとも排気通路２上に還元剤（具体的にはＣＯ）濃度を検出するセンサ（ＣＯセンサ）を設け、このＣＯセンサから得られる還元剤濃度［ｘ′］の実測値と、式（６）から得られる還元剤濃度［ｘ′］の推定値とから各係数ｍ_ｉを逐次更新して還元剤利用率ｒ′を求め、この値を上述した式（４）に代入することでＮＯｘ放出量を求めることができる。
【００４４】
そして、このようにしてＮＯｘ放出量を求めることでＮＯｘ放出量の推定精度を高めることができる。
【００４５】
【発明の効果】
以上詳述したように、本発明のＮＯｘ吸蔵量の推定方法によれば、ＮＯｘ吸蔵量を、ＮＯｘ吸蔵触媒のＮＯｘ吸蔵特性を反映させた多項式を用いて推定するとともに、多項式の係数を、実測されたＮＯｘ浄化率に基づいて逐次補正するので、高い精度でＮＯｘ吸蔵量を推定することができるという利点がある。また、ＮＯｘ吸蔵量を正確に推定することでＮＯｘ吸蔵量に基づいた最適なリッチ運転制御を実行することができ、燃費を向上させることができる。また、ＮＯｘ吸蔵触媒を変更した場合でも、変更後の触媒の特性に応じた各係数に更新されるので、モデル式を変更する必要がないという利点がある（請求項１）。
また、ＮＯｘ放出量を還元剤利用率に基づいて算出するので、比較的高い精度でＮＯｘ放出量を推定することができる（請求項２）。
また、還元剤利用率を触媒のＮＯｘ放出特性を反映させた多項式を用いて推定するとともに、多項式の係数を還元剤の濃度に基づいて逐次補正するので、ＮＯｘ放出量の推定精度をさらに高めることができるという利点がある（請求項３）。
また、リッチ運転中は多項式の係数を保持するとともに、リーン運転開始時に推定ＮＯｘ浄化率と、実ＮＯｘ浄化率との差がしきい値以上であるとＮＯｘ吸蔵量を補正することにより、ＮＯｘ吸蔵量の推定精度をさらに高めることができる（請求項４）。
また、多項式における各係数の所定期間の平均値が所定の範囲を逸脱すると、触媒が異常である判定するので、触媒の異常判定を高い精度で実行することができる利点がある（請求項５）。
【図面の簡単な説明】
【図１】本発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法について説明するための模式的なブロック図である。
【図２】本発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法について説明するためのマップであって、ＮＯｘ触媒のＮＯｘ吸蔵時の特性を示すマップである。
【図３】本発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法について説明するためマップであって、ＮＯｘ触媒のＮＯｘ放出時の特性を示すマップである。
【図４】本発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法について説明するためのフローチャートである。
【図５】本発明の一実施形態にかかるＮＯｘ吸蔵量の推定方法について説明するためのフローチャートである。
【図６】従来の技術について説明する図である。
【符号の説明】
１ＮＯｘ吸蔵触媒（触媒）
２排気通路
３ＮＯｘセンサ（ＮＯｘ濃度検出手段）
４ＥＣＵ[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a method for estimating the amount of NOx stored in a NOx storage-type reduction catalyst provided in an exhaust passage.
[0002]
[Prior art]
In general, in a diesel engine, due to its structural characteristics, exhaust gas contains a large amount of oxygen and the atmosphere becomes an oxidizing atmosphere (that is, the exhaust air-fuel ratio becomes lean). Also, if a gasoline engine is capable of lean burn, the exhaust air-fuel ratio will also be lean during lean burn. Therefore, a lean NOx catalyst has been conventionally developed so that NOx can be purified even in such an oxidizing atmosphere.
[0003]
As such a lean NOx catalyst, there is a NOx storage type reduction catalyst (hereinafter, referred to as a NOx storage catalyst or simply a catalyst). In the NOx storage catalyst, NOx in the exhaust gas is stored in an oxidizing atmosphere, and the stored NOx is released when CO is present under a low oxygen concentration. Most of the NOx released from the NOx storage catalyst is released. The gas is reduced to harmless N ₂ by a reducing agent such as unburned HC and CO discharged from the engine and then discharged to the atmosphere.
[0004]
In a vehicle equipped with a NOx storage catalyst, generally, so-called rich spike control is performed. In other words, if the operation is continued at a lean air-fuel ratio, the NOx storage catalyst is in a saturated state in which NOx can no longer be stored, so that the air-fuel ratio is temporarily forcibly made rich at an appropriate timing (rich spike control is performed. This reduces the concentration of oxygen in the exhaust gas, and releases NOx from the NOx storage catalyst by supplying a reducing agent, thereby regenerating the NOx storage catalyst to a state where NOx can be stored.
[0005]
By the way, in order to perform the rich spike control as described above, it is necessary to accurately estimate (or detect) the NOx storage amount in the NOx storage catalyst. As a method for estimating the NOx storage amount, for example, the following method is used. Techniques are known.
(1) First prior art: NOx occlusion amount estimation method using NOx sensor In this first prior art, as shown in FIG. 6, NOx concentration is detected upstream and downstream of the catalyst 101, respectively. In addition to providing the sensors 102 and 103, an air flow sensor (not shown) for detecting the exhaust gas flow rate is provided in the exhaust passage.
[0006]
Then, the ECU 104 estimates the NOx storage amount by the following equation based on the information from the NOx sensors 102 and 103 and the information from the air flow sensor.
NOx storage amount = exhaust gas flow rate x (catalyst upstream NOx concentration-catalyst downstream NOx concentration)
{Circle around (2)} Second prior art: estimation method using a mathematical model (for example, Patent Document 1)
In the second prior art, a mathematical catalyst model based on the chemical phenomena and physical phenomena of a catalyst such as an occlusion reaction, an oxidation-reduction reaction, and a desorption reaction is provided, and the NOx occlusion amount is estimated using the catalyst model formula.
[0007]
[Patent Document 1]
JP-A-9-72235
[Problems to be solved by the invention]
However, in the first prior art, the NOx storage amount can be calculated during the lean operation, but the NOx release amount (decrease in the storage amount) during the rich operation cannot be detected by the NOx sensor. For this reason, when the rich operation time is insufficient, there is a problem that the NOx remaining in the catalyst is accumulated and a calculation error of the NOx storage amount occurs. In addition, there is a problem that NOx sensors are required upstream and downstream of the catalyst, respectively, which causes an increase in cost.
[0009]
Further, in the second conventional technique, the NOx storage amount can be estimated only by providing at least one NOx sensor downstream of the catalyst. However, since this technique is mainly applied to the air-fuel ratio control of the three-way catalyst. However, there is a problem that the model structure is different from the NOx storage catalyst, and it is difficult to apply the NOx storage catalyst to the NOx storage catalyst. Also, when estimating the NOx occlusion amount using other catalyst models, it is necessary to change the model formula because the characteristic values are different for each catalyst type. There is a problem that requires modeling.
[0010]
The present invention has been made in view of such a problem, and has been proposed in order to estimate the NOx storage amount of the NOx storage catalyst with high accuracy by always reflecting the latest state in the catalyst model. The aim is to provide a method.
[0011]
[Means for Solving the Problems]
For this reason, the method for estimating the NOx storage amount according to the first aspect of the present invention is a method for estimating the NOx storage amount of the NOx storage catalyst provided in the exhaust passage of the engine. The estimation is performed using a polynomial that reflects the NOx storage characteristics of the NOx storage catalyst, and the coefficient of each polynomial is sequentially corrected based on the actually measured NOx purification rate.
[0012]
Further, the NOx release amount of the NOx storage catalyst may be calculated by the following equation. NOx emission amount = ∫ (catalyst inlet reducing agent concentration × reducing agent utilization rate−0.5 × catalyst inlet oxygen concentration) × exhaust gas flow rate Preferably, the reducing agent utilization rate reflects the NOx emission characteristics of the NOx storage catalyst. And the coefficients of the polynomial are sequentially corrected based on the concentration of the reducing agent.
[0013]
Preferably, the engine is configured to be capable of switching between a lean operation in which the exhaust air-fuel ratio is lean and a rich operation in which the exhaust air-fuel ratio is rich, and holds the coefficient of the polynomial during the rich operation. At the same time, when the difference between the NOx purification rate obtained by using the held coefficient at the start of the lean operation and the actually measured NOx purification rate is equal to or larger than a threshold value, the NOx storage amount is corrected.
[0014]
Furthermore, when any one of the average values of the coefficients deviates from the set predetermined range, the catalyst may be determined to be abnormal.
[0015]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, a method for estimating the NOx occlusion amount according to an embodiment of the present invention will be described with reference to the drawings. FIG. 1 is a schematic diagram showing the entire configuration of an exhaust system of a vehicle to which the present invention is applied.
As shown in the drawing, a NOx storage catalyst 1 is interposed in the exhaust passage 2, and a NOx sensor 3 for detecting the NOx concentration at the catalyst outlet is provided downstream of the NOx storage catalyst 1.
[0016]
Although not shown, an O ₂ sensor, a temperature sensor, an air flow sensor (AFS), and the like are connected to the exhaust passage 2, and the oxygen concentration, the exhaust gas temperature, and the exhaust gas flow rate in the exhaust passage 2 are detected by these sensors. Each is detected.
Each sensor 3 is connected to an ECU 4 as control means. Here, the ECU 4 includes an input / output device, a storage device (ROM, RAM, non-volatile RAM, etc.), an arithmetic device (CPU), a timer counter, and the like. The ECU 4 allows the NOx storage catalyst 1 to store NOx. An estimate of the quantity is made.
[0017]
The ECU 4 is provided with various maps. For example, the concentrations of NOx, CO, and HC at the catalyst inlet (hereinafter referred to as catalyst inlet [NOx] and catalyst inlet [ CO] and the catalyst inlet [HC]) are read from the map. In the present embodiment, these values are obtained from the map, but various sensors may be provided in the exhaust passage 2 and these values may be directly detected from the sensors.
[0018]
By the way, in the prior art 2 described above, the characteristic itself of the NOx storage catalyst is described by a mathematical model to determine the NOx storage amount, whereas in the present invention, the NOx storage amount is calculated using a four-dimensional linear polynomial. It is supposed to.
That is, the ECU 4 stores a polynomial as shown in the following equation (1), and the NOx purification rate and the NOx storage ratio in the NOx storage catalyst 1 are calculated based on this polynomial.

In the equation (1), r is the NOx purification rate, x is the NOx storage ratio, y is the exhaust gas temperature, z is the SV value (or the exhaust gas flow rate), and k ₀ , k ₁ , k ₂ ,. is there.
[0019]
Here, the NOx storage characteristics of the NOx storage catalyst 1 can be known in advance by experiments or the like, and can be approximated by a polynomial. For example, if the number of parameters is three (three-dimensional), the characteristics of the catalyst can be represented by a curved surface as shown in FIG. 2, and this can be represented by a polynomial. Of course, since equation (1) is four-dimensional, it does not have a curved surface as shown in FIG. 2, but the concept is the same as in the case of three-dimensional.
[0020]
For each coefficient k _i (i = 1, 2,...), An appropriate value is input as an initial value from characteristics previously obtained through experiments or the like.
Then, during engine operation, always by sequentially correcting the coefficients k _i based on the value detected by the catalyst downstream of the NOx sensor 3 by the equation (1), it is possible to represent the latest catalytic state.
[0021]
Hereinafter, the method of NOx estimation will be described in more detail. First, the NOx purification rate r in the equation (1) can also be obtained by the following equation (2).
r = NOx concentration at catalyst outlet / NOx concentration at catalyst inlet ... (2)
Here, the catalyst outlet NOx concentration is a value detected by the NOx sensor 3, and the catalyst inlet NOx concentration is a value obtained from a map. Therefore, the NOx purification rate r obtained by the equation (2) can be called a sensor value (actually measured value), and the NOx purification rate r obtained by the equation (1) can be called an estimated value (calculated value).
[0022]
As the exhaust gas temperature [y] and the SV value [z] in the equation (1), measured data or map values are applied. Thus, the only unknown value in the equation (1) is the NOx storage ratio [x], and the other values are known. Therefore, the NOx occlusion ratio [x] can be obtained by the following equation (3) obtained by modifying the equation (1).
x = _{_{[r- (k 0 + k 2 y}} + k 3 z ··) ] _{_{/ (k 1 + k 4 y}} + ··) ·· (3)
Then, the NOx storage ratio [x] obtained by the equation (3) is returned to the equation (1) and used for the calculation in the next calculation cycle. That is, in equation (1), the purification rate r is calculated using the NOx storage ratio [x] obtained last time and [y] and [z] newly detected in the current calculation cycle. It is.
[0023]
Here, a purification rate r is the estimated value calculated by Equation (1), if the coefficient k _i is the exact value, the estimated value and the measured value obtained by the equation (2) of the purification rate r Should match.
Therefore, in the present invention, the estimated value and the measured value of the purification rate r are compared, and if there is a difference between the estimated value and the measured value, each coefficient k _i is minimized so that the estimated value becomes the measured value. The correction is performed sequentially using the square method. Then, by repeatedly performing such calculation and sequentially correcting each coefficient k _i , the polynomial of Expression (1) can be changed to an expression that accurately reflects the state of the NOx storage catalyst 1 as needed.
[0024]
That, NOx occlusion catalyst 1 is one in which its properties by usage and aging and the like are changed by correcting repeatedly even coefficients k _i times represents formula (1) is exactly of the catalyst 1 state Thus, the NOx storage amount can be estimated with high accuracy. The NOx occlusion ratio [x] obtained by the equation (3) is an instantaneous NOx occlusion amount in the calculation cycle at the beginning, but this calculation result is repeatedly returned to the equation (1) to calculate the purification rate r. By doing so, the integrated value of NOx stored in the catalyst 1 can be calculated.
[0025]
As described above, during the lean operation for storing NOx, the catalyst outlet [NOx] is taken in from the NOx sensor 3 provided downstream of the catalyst, and is compared with the estimated value obtained by the equation (2). Correcting the coefficient k _i by the multiplicative method makes it possible to obtain an accurate NOx storage amount.
By the way, since NOx is released during the rich operation, if the accurate NOx release amount cannot be detected or estimated, the initial value of the NOx occlusion amount in the next lean operation cannot be accurately calculated.
[0026]
Therefore, along with being rich operation keeps the _{k i,} the following equation (4), (5) to calculate the NOx emission amount, equation (4) from the NOx occlusion amount at the time immediately before the start rich operation, calculated by (5) The calculated NOx release amount is subtracted to calculate the NOx storage amount during the rich operation.
NOx emission amount = ∫ (catalyst inlet reducing agent concentration × reducing agent utilization rate [r ′] − 0.5 × catalyst inlet [O ₂ ]) × exhaust gas flow rate (4)
Reducing agent utilization rate [r '] = f (y, z) (5)
Here, the ECU 1 is provided with a reducing agent utilization rate setting map as shown in FIG. This map stores characteristics of the reducing agent utilization based on the exhaust gas temperature [y] and the SV value [z], and sets the reducing agent utilization based on these parameters [y] and [z]. It has become so.
[0027]
Further, the catalyst inlet reducing agent concentration (catalyst inlet [CO], catalyst inlet [HC]) can be obtained from the map using the accelerator opening and the engine speed as parameters as described above. Value or detected value.
Accordingly, when switched to rich operation from the lean operation, stores the NOx occlusion amount holds each coefficient k _i during the lean operation end in ECU 4, when switched to the lean operation the rich operation again, the NOx storage The remaining amount of NOx is calculated by subtracting the NOx release amount at the end of the rich operation from the amount.
[0028]
The ECU 4 corrects the NOx storage amount when the difference between the measured value and the estimated value of the NOx purification rate r at the start of the lean operation is equal to or larger than a threshold value. That is, immediately after switching from the rich operation to the lean operation, if the difference between the NOx purification rate estimated by the equation (1) and the NOx purification rate obtained by the equation (2) is equal to or more than a predetermined value, the lean operation is started. It is determined that the calculated residual NOx amount (= NOx storage amount at the end of the previous lean operation-NOx release amount at the end of the rich operation) is not accurate, and this NOx residual amount [x] (that is, the lean operation start) NOx storage amount at the time) is corrected based on the equation (3).
[0029]
Further, ECU 4 also has unexpected when coefficients by abnormality of the catalyst 1 k _i is the value is not normal, functional informing the driver detects this. Specifically, for each coefficient k _i , a moving average in a predetermined continuous calculation cycle is constantly calculated, and when the value of the moving average deviates from a predetermined range set for each coefficient, the catalyst 1 is deteriorated or deteriorated. It is determined that an abnormality such as breakage has occurred. The predetermined range is, for example, ± alpha to the initial value of each coefficient k _i (alpha is a constant value) is adapted to be set by adding a.
[0030]
Further, when the ECU 4 determines that the catalyst 1 is abnormal, for example, the warning light on the instrument panel is turned on and the rich operation is prohibited. In this way, by prohibiting the rich operation when the catalyst 1 is abnormal, it is possible to avoid a situation in which CO is discharged.
Since the method for estimating the NOx storage amount according to one embodiment of the present invention is configured as described above, the NOx storage amount is estimated as follows.
[0031]
First, during the lean operation, the actual NOx purification rate r is calculated from the equation (2), and the exhaust gas temperature [y] and the SV value [z] are obtained from the measured data or the map value. In the polynomial of the equation (1), an appropriate value obtained by an experiment or the like is input as an initial value for each coefficient k _i . Then, the respective values are substituted into Expression (3) obtained by modifying Expression (1), and the NOx storage ratio [x] is calculated.
[0032]
Next, the NOx storage ratio [x] obtained by the equation (3) and the newly obtained exhaust gas temperature [y] and SV value [z] are substituted into the equation (1) to obtain the NOx purification rate r (estimated). Value). The minimum and the estimated NOx purification ratio r, is compared with the actual NOx purification rate r newly calculated by the equation (2), the estimated NOx purification rate, each coefficient k _i to match the actual NOx purification rate Correct by the square method.
[0033]
Then, such calculations the repeatedly executed, so that each coefficient k _i Equation (1) by updated at any time becomes equation accurately represent the state of the catalyst 1, it is possible to estimate the NOx adsorption amount with high accuracy become. Further, by repeatedly calculating the NOx storage ratio [x] using the equation (3), the integrated value of the NOx stored in the catalyst 1 can be calculated.
[0034]
Furthermore, since during rich operation NOx is released, along with storing the respective coefficients k _i to interrupt the calculation of the NOx occlusion amount by the formula (1) in this case, the above-mentioned formula (4), (5) To calculate the NOx release amount.
When the operation is switched from the rich operation to the lean operation again, the NOx remaining amount at the start of the lean operation is calculated by subtracting the NOx release amount at the end of the rich operation from the NOx storage amount at the end of the lean operation.
[0035]
If the difference between the measured value and the estimated value of the NOx purification rate r at the start of the lean operation is equal to or larger than the threshold value, the NOx storage amount is corrected. Here, it will be described with reference to the flowchart of FIG method for correcting the NOx occlusion amount, firstly, sequentially corrected and identifying each coefficient k _i of each coefficient k _i of the polynomial (1) during lean operation (step S1) .
When it is determined that the rich operation is started (step S2), and holds the respective coefficient k _i (step S3), and thereafter, when the rich operation is completed (step S4), and the measured value of the NOx purification rate The estimated values are compared with each other, and it is determined whether these differences are equal to or larger than a threshold value (step S5).
[0036]
If the difference is less than the threshold value between the actual measurement value and the estimated value of the NOx purification rate, and continue to update the coefficients k _i during lean operation returns to Step 1 again, NOx occlusion if equal to or greater than the threshold value The amount is corrected to the value calculated by equation (3) (step S6), and the routine returns.
By correcting the NOx storage amount in this way, the estimation accuracy of the NOx storage amount can be further improved.
[0037]
On the other hand, when it is determined that the catalyst 1 has deteriorated or is abnormal, the warning light on the instrument panel is turned on to inform the driver and rich operation is prohibited. This will be described with reference to the flowchart of FIG. 5, first, the polynomial (1) sequentially corrects each coefficient k _i of identifying the coefficients k _i (step S11), and then a period of time of each coefficient k _i A moving average is calculated (step S12).
[0038]
Next, it is determined whether or not the calculated moving average of each coefficient k _i is within a predetermined range (step S13). If the moving average of each coefficient k _i is within the predetermined range, step S13 is performed again. returning to S2, if there in a predetermined range in one of the moving average of the coefficients k _i, it is determined that the abnormality in the catalyst 1 resulted perform an alarm (step S14), and prohibits the rich operation (Step S15).
As described in detail above, according to the method for estimating the NOx storage amount according to one embodiment of the present invention, the NOx storage amount is estimated using a polynomial that reflects the NOx storage characteristics of the NOx storage catalyst 1, and Since the coefficient k _i of the polynomial is successively corrected based on the actually measured NOx purification rate, it is possible to always use a polynomial that models the latest catalyst state, and to estimate the NOx storage amount with high accuracy. There are advantages. Further, by accurately estimating the NOx storage amount, it is possible to execute optimal rich operation control based on the NOx storage amount, and to improve fuel efficiency.
[0039]
Also, when you change the NOx storage catalyst 1, therefore it is updated each coefficient k _i corresponding to the characteristics of the catalyst 1 after the change, there is an advantage that there is no need to change the model equation. On the other hand, when a mathematical model formula is applied as in the related art, there is a problem that an operation (adjustment) of specifying a coefficient according to the characteristics of each catalyst is required, and a great deal of labor is required. According to the present invention, the NOx storage amount can be accurately estimated without performing such adjustment.
[0040]
Since the NOx emission amount is calculated based on the reducing agent utilization rate, the NOx emission amount can be estimated with relatively high accuracy.
In addition, holds the coefficients k _i during rich operation polynomials, the NOx purification rate is estimated using a coefficient which is held at the lean operation start, the difference between the actually measured NOx purification rate than the threshold value If there is, the NOx storage amount is corrected, so that the estimation accuracy of the NOx storage amount can also be improved.
[0041]
If the average value in a predetermined period of each coefficient k _i deviates from the predetermined range, determines that the catalyst 1 is abnormal, it can determine an abnormality of the catalyst 1 with high accuracy. Further, when such an abnormality is determined, an alarm is issued to the driver, so that the driver can quickly recognize the abnormality of the catalyst 1 and prompt the driver to replace the catalyst 1 at an early stage. In addition, when the catalyst 1 is abnormal, the rich operation is prohibited, so that a situation in which CO is directly discharged can be avoided.
[0042]
Note that the method of estimating the NOx storage amount of the present invention is not limited to the above embodiment, and various modifications can be made without departing from the spirit of the present invention. For example, the amount of released NOx may be obtained by a method similar to the method of estimating the amount of stored NOx.
That is, in the above description, the reductant utilization rate [r '] used in calculating the NOx release amount is obtained from the reductant utilization rate setting map shown in FIG. 3, but as shown in the following equation (6), It may be obtained from a polynomial of the reducing agent concentration [x '], the exhaust gas temperature [y], and the SV value [z].

Note that _mi (i = 1, 2,...) Is a coefficient. Further, equation (6) is a polynomial which reflects the NOx emission characteristics of NOx storage catalyst 1, each coefficient m _i, equation (1) and the appropriate value is an initial value obtained from the same experiment, etc. Is entered as
[0043]
In this case, a sensor (CO sensor) for detecting the concentration of the reducing agent (specifically, CO) is provided at least on the exhaust passage 2, and the measured value of the reducing agent concentration [x '] obtained from this CO sensor is From the estimated value of the reducing agent concentration [x '] obtained from the equation (6), each coefficient _mi is successively updated to obtain the reducing agent utilization rate r', and this value is substituted into the above-described equation (4). Thus, the NOx release amount can be obtained.
[0044]
By obtaining the NOx release amount in this way, the estimation accuracy of the NOx release amount can be improved.
[0045]
【The invention's effect】
As described in detail above, according to the NOx storage amount estimation method of the present invention, the NOx storage amount is estimated using a polynomial that reflects the NOx storage characteristics of the NOx storage catalyst, and the coefficient of the polynomial is measured. Since the correction is sequentially performed based on the determined NOx purification rate, there is an advantage that the NOx storage amount can be estimated with high accuracy. Further, by accurately estimating the NOx storage amount, it is possible to execute optimal rich operation control based on the NOx storage amount, and to improve fuel efficiency. Further, even when the NOx storage catalyst is changed, the coefficient is updated to each coefficient according to the changed characteristics of the catalyst, so that there is an advantage that it is not necessary to change the model formula (claim 1).
Further, since the NOx emission amount is calculated based on the reducing agent utilization rate, the NOx emission amount can be estimated with relatively high accuracy (claim 2).
Further, since the reducing agent utilization rate is estimated using a polynomial that reflects the NOx emission characteristics of the catalyst, and the coefficient of the polynomial is sequentially corrected based on the concentration of the reducing agent, the estimation accuracy of the NOx emission amount is further improved. (Claim 3).
In addition, during the rich operation, the polynomial coefficient is held, and when the difference between the estimated NOx purification rate and the actual NOx purification rate at the start of the lean operation is equal to or larger than a threshold value, the NOx storage amount is corrected to thereby store NOx. The accuracy of estimating the amount can be further improved (claim 4).
Further, when the average value of each coefficient in the polynomial in a predetermined period deviates from a predetermined range, it is determined that the catalyst is abnormal. Therefore, there is an advantage that the abnormality determination of the catalyst can be performed with high accuracy. .
[Brief description of the drawings]
FIG. 1 is a schematic block diagram for explaining a method for estimating an NOx storage amount according to an embodiment of the present invention.
FIG. 2 is a map for explaining a method of estimating a NOx storage amount according to one embodiment of the present invention, and is a map showing characteristics of a NOx catalyst during NOx storage.
FIG. 3 is a map for explaining a method of estimating the NOx storage amount according to one embodiment of the present invention, and is a map showing characteristics of the NOx catalyst at the time of releasing NOx.
FIG. 4 is a flowchart illustrating a method for estimating the amount of stored NOx according to an embodiment of the present invention.
FIG. 5 is a flowchart illustrating a method for estimating an NOx storage amount according to an embodiment of the present invention.
FIG. 6 is a diagram illustrating a conventional technique.
[Explanation of symbols]
1 NOx storage catalyst (catalyst)
2 Exhaust passage 3 NOx sensor (NOx concentration detecting means)
4 ECU

Claims

A method for estimating the NOx storage amount of a NOx storage catalyst provided in an exhaust passage of an engine,
The NOx storage amount is estimated using a polynomial that reflects the NOx storage characteristics of the NOx storage catalyst,
A method for estimating a NOx storage amount, wherein each coefficient of the polynomial is sequentially corrected based on an actually measured NOx purification rate.

2. The method according to claim 1, wherein the NOx release amount in the NOx storage catalyst is calculated by the following equation.
NOx emission amount = ∫ (catalyst inlet reducing agent concentration x reducing agent utilization rate-0.5 x catalyst inlet oxygen concentration) x exhaust gas flow rate

Estimating the reducing agent utilization rate using a polynomial that reflects the NOx release characteristics of the NOx storage catalyst,
3. The method according to claim 2, wherein the coefficient of the polynomial is corrected sequentially based on the concentration of the reducing agent.

The engine is configured to be capable of switching between a lean operation in which the exhaust air-fuel ratio is lean and a rich operation in which the exhaust air-fuel ratio is rich,
During the rich operation, the coefficient of the polynomial is retained, and the difference between the NOx purification rate obtained using the retained coefficient at the start of the lean operation and the actually measured NOx purification rate is equal to or greater than a threshold value. 2. The method for estimating a NOx storage amount according to claim 1, wherein the NOx storage amount is corrected.

6. The NOx storage amount estimating method according to claim 1, wherein the catalyst is determined to be abnormal when the average value of each coefficient in a predetermined period deviates from a predetermined range.