(57)【要約】【課題】常に最新の触媒状態をモデル化した多項式を用いつつ、ＮＯｘ吸蔵量を高精度で推定するＮＯｘ吸蔵量の推定方法を提供する。【解決手段】エンジンの排気通路に介装されたＮＯｘ吸蔵触媒のＮＯｘ吸蔵量を推定する方法であって、ＮＯｘ吸蔵触媒のＮＯｘ吸蔵特性を反映させた多項式を用いてＮＯｘ吸蔵量が推定される工程(12)を備え、この工程は、排気通路にて実測されたＮＯｘ浄化率と、ＮＯｘ吸蔵触媒の触媒特性が変化してもその値が変化しないポイントのＮＯｘ浄化率とに基づいて多項式の各係数が逐次補正される工程(18)を有する。【選択図】図６
2007263048.tif 000002
【特許請求の範囲】
【請求項１】
エンジンの排気通路に介装されたＮＯｘ吸蔵触媒のＮＯｘ吸蔵量を推定する方法であって、
前記ＮＯｘ吸蔵触媒のＮＯｘ吸蔵特性を反映させた多項式を用いて前記ＮＯｘ吸蔵量が推定される工程を備え、
該工程は、前記排気通路にて実測されたＮＯｘ浄化率と、前記ＮＯｘ吸蔵触媒の触媒特性が変化してもその値が変化しないポイントのＮＯｘ浄化率とに基づいて前記多項式の各係数が逐次補正される工程を有すること
を特徴とするＮＯｘ吸蔵量の推定方法。
【請求項２】
前記ポイントのＮＯｘ浄化率は、略零の値であることを特徴とする請求項１に記載のＮＯｘ吸蔵量の推定方法。
【請求項３】
前記各係数は、以下の適応推定ロジックによって逐次補正されることを特徴とする請求項１又は２に記載のＮＯｘ吸蔵量の推定方法。
ｋｎｅｗ＝ｋｏｌｄ＋δ｛α０（ｒ０−ｋｏｌｄｐ０Ｔ）ｐ０＋α１（ｒ１−ｋｏｌｄｐ１Ｔ）ｐ１＋・・・＋αＮ（ｒＮ−ｋｏｌｄｐＮＴ）ｐＮ｝
ここで、ｋｎｅｗは多項式の各項の係数の更新値のベクトル、ｋｏｌｄは多項式の各項の係数の前回計算値のベクトル、ｒｉはＮＯｘ浄化率、ｐｉはＮＯｘ吸蔵割合ｘや排ガス温度ｙやＳＶ値ｚ等の各入力値のベクトル、Ｔは転置行列、δはゲイン、αｉは実測によって収集されたデータやＮ個の特定された前記ポイントのデータの各重み付け係数を示す。
【請求項４】
エンジンの排気通路に介装され、リーン運転時に排ガス中のＮＯｘを吸蔵するとともに該吸蔵したＮＯｘをリッチ運転時に放出還元するＮＯｘ吸蔵触媒と、
該ＮＯｘ吸蔵触媒のＮＯｘ吸蔵特性を反映させた多項式を用いてＮＯｘ吸蔵量を推定する吸蔵量推定部とを具備し、
該吸蔵量推定部は、前記排気通路にて実測されたＮＯｘ浄化率と、前記ＮＯｘ吸蔵触媒の触媒特性が変化してもその値が変化しないポイントのＮＯｘ浄化率とに基づいて前記多項式の各係数を逐次補正する係数同定部を含む
ことを特徴とするＮＯｘ吸蔵触媒の吸蔵量推定装置。
【請求項５】
前記ポイントのＮＯｘ浄化率は、略零の値であることを特徴とする請求項４に記載のＮＯｘ吸蔵触媒の吸蔵量推定装置。
【請求項６】
前記各係数は、以下の適応推定ロジックによって逐次補正されることを特徴とする請求項４又は５に記載のＮＯｘ吸蔵触媒の吸蔵量推定装置。
ｋｎｅｗ＝ｋｏｌｄ＋δ｛α０（ｒ０−ｋｏｌｄｐ０Ｔ）ｐ０＋α１（ｒ１−ｋｏｌｄｐ１Ｔ）ｐ１＋・・・＋αＮ（ｒＮ−ｋｏｌｄｐＮＴ）ｐＮ｝
ここで、ｋｎｅｗは多項式の各項の係数の更新値のベクトル、ｋｏｌｄは多項式の各項の係数の前回計算値のベクトル、ｒｉはＮＯｘ浄化率、ｐｉはＮＯｘ吸蔵割合ｘや排ガス温度ｙやＳＶ値ｚ等の各入力値のベクトル、Ｔは転置行列、δはゲイン、αｉは実測によって収集されたデータやＮ個の特定された前記ポイントのデータの各重み付け係数を示す。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、排気通路に介装されたＮＯｘ吸蔵型還元触媒（以下、ＮＯｘ吸蔵触媒と称する）に対するＮＯｘ吸蔵量の推定方法に関する。
【背景技術】
【０００２】
一般に、ＮＯｘ吸蔵触媒は、排気空燃比がリーンのときに排ガス中のＮＯｘ（窒素酸化物）を吸蔵し、排気空燃比がリッチのときに吸蔵したＮＯｘを放出還元する。具体的には、この触媒は、酸素過剰状態（酸化雰囲気）において排ガス中のＮＯｘを硝酸塩として吸蔵し、この吸蔵したＮＯｘを一酸化炭素過剰状態（還元雰囲気）において窒素に還元させる特性を有している。
【０００３】
そして、この触媒を装備したエンジンでは、ＮＯｘ吸蔵量が飽和に至る前にリッチ運転へ間欠的に切り換えるリッチスパイクを行い、ＮＯｘ吸蔵量の増加に伴う触媒の性能低下を抑制させる。これにより、触媒が再生され、排ガスが良好に浄化される。
ここで、このリッチスパイクを行うには、ＮＯｘ吸蔵量を正確に推定又は検出する必要がある。そのため、触媒の化学・物理反応に基づく数学的な触媒モデルを用い、このモデル式を用いてＮＯｘ吸蔵量を推定する技術が開示されている（例えば、特許文献１）。
【特許文献１】特開平９−７２２３５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところで、当該技術は三元触媒の空燃比制御に適用させることを主眼としている。つまり、モデル構造がＮＯｘ吸蔵触媒とは異なり、ＮＯｘ吸蔵触媒への応用が困難である。また、その他の触媒モデルを用いてＮＯｘ吸蔵量を推定する場合には、モデル式を逐次変更しなければならない。触媒種毎にその特性値が異なるからである。しかも、触媒の劣化等に対応させるためには、より詳細なモデル化が必要となる。このように、上述した技術では、ＮＯｘ吸蔵量の推定は可能であったとしても、ＮＯｘ吸蔵量を高精度で推定する点に関しては依然として課題が残されている。
【０００５】
ここで、ＮＯｘ吸蔵量を高精度で推定するためには、ＮＯｘ吸蔵特性を反映させた多項式を用い、そして、実測によって収集されたデータに基づいて常に最新の状態をモデル化し、ＮＯｘ吸蔵量を推定する方法が考えられる。しかしながら、この方法を用いる場合には、上記収集されたデータのみに拘泥されてはならない点に留意する必要がある。エンジン運転範囲のうち狭い領域のデータが収集され、当該データのみに基づいて触媒モデルを最新化すると、触媒モデルがこの収集された特定の運転範囲へのフィッティングに支配され、他の運転範囲における触媒モデルの推定精度が低下するとの懸念があるからである。
【０００６】
本発明は、このような課題に鑑みてなされたもので、常に最新の触媒状態をモデル化した多項式を用いつつ、ＮＯｘ吸蔵量を高精度で推定するＮＯｘ吸蔵量の推定方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記の目的を達成すべく、請求項１記載のＮＯｘ吸蔵量の推定方法は、エンジンの排気通路に介装されたＮＯｘ吸蔵触媒のＮＯｘ吸蔵量を推定する方法であって、ＮＯｘ吸蔵触媒のＮＯｘ吸蔵特性を反映させた多項式を用いてＮＯｘ吸蔵量が推定される工程を備え、この工程は、排気通路にて実測されたＮＯｘ浄化率と、ＮＯｘ吸蔵触媒の触媒特性が変化してもその値が変化しないポイントのＮＯｘ浄化率とに基づいて多項式の各係数が逐次補正される工程を有することを特徴としている。
【０００８】
また、請求項４記載のＮＯｘ吸蔵触媒の吸蔵量推定装置は、エンジンの排気通路に介装され、リーン運転時に排ガス中のＮＯｘを吸蔵するとともに吸蔵したＮＯｘをリッチ運転時に放出還元するＮＯｘ吸蔵触媒と、ＮＯｘ吸蔵触媒のＮＯｘ吸蔵特性を反映させた多項式を用いてＮＯｘ吸蔵量を推定する吸蔵量推定部とを具備し、吸蔵量推定部は、排気通路にて実測されたＮＯｘ浄化率と、ＮＯｘ吸蔵触媒の触媒特性が変化してもその値が変化しないポイントのＮＯｘ浄化率とに基づいて多項式の各係数を逐次補正する係数同定部を含むことを特徴としている。
【０００９】
また、請求項２，５記載の発明では、ポイントのＮＯｘ浄化率は、略零の値であることを特徴としている。
更に、請求項３，６記載の発明では、各係数は、以下の適応推定ロジックによって逐次補正されることを特徴としている。
ｋｎｅｗ＝ｋｏｌｄ＋δ｛α０（ｒ０−ｋｏｌｄｐ０Ｔ）ｐ０＋α１（ｒ１−ｋｏｌｄｐ１Ｔ）ｐ１＋・・・＋αＮ（ｒＮ−ｋｏｌｄｐＮＴ）ｐＮ｝
ここで、ｋｎｅｗは多項式の各項の係数の更新値のベクトル、ｋｏｌｄは多項式の各項の係数の前回計算値のベクトル、ｒｉはＮＯｘ浄化率、ｐｉはＮＯｘ吸蔵割合ｘや排ガス温度ｙやＳＶ値ｚ等の各入力値のベクトル、Ｔは転置行列、δはゲイン、αｉは実測によって収集されたデータやＮ個の特定された前記ポイントのデータの各重み付け係数を示す。
【発明の効果】
【００１０】
従って、請求項１記載の本発明のＮＯｘ吸蔵量の推定方法や、請求項４記載のＮＯｘ吸蔵触媒の吸蔵量推定装置によれば、ＮＯｘ吸蔵特性を反映させた多項式を用いてＮＯｘ吸蔵量を推定するにあたり、この多項式の各係数が、実測によって収集されたＮＯｘ浄化率の他、触媒特性が変化してもその値が変化しないポイントのＮＯｘ浄化率に基づいて逐次補正されている。つまり、仮に、実測されたＮＯｘ浄化率が狭い領域に集中したデータであってとしても、この狭い領域から離れた他の領域のデータをも考慮した多項式の各係数が同定可能となり、この多項式の各係数は大域的なデータに基づいた値になる。
【００１１】
そして、上述した他の領域のデータには不変のデータが用いられており、当該他の領域における誤差は解消されるので、各係数の同定計算が安定するし、その計算精度も向上する。
これらの結果、多項式による触媒モデルが、収集された特定の運転範囲へのフィッティングに支配され難くなり、ＮＯｘ吸蔵量の推定精度の向上に寄与する。
【００１２】
また、請求項２や請求項５記載の発明によれば、上述した他の領域では、例えば排ガス温度が低温時或いは高温時の如く、ＮＯｘ吸蔵触媒が劣化してもＮＯｘ浄化率が略零のままであるポイントのデータが用いられている。従って、当該他の領域における誤差は、触媒特性が変化しても確実に解消されるので、ＮＯｘ吸蔵量の推定精度の向上により一層寄与する。
【００１３】
更に、請求項３や請求項６記載の発明によれば、右辺第１項は前回の係数であり、同第２項中の添字０で示される部分は実測によって収集されたデータに相当する。そして、同項中の添字１〜Ｎで示される部分は、上述した他の領域にＮ個の特定されたポイントのデータに相当する。従って、当該他の領域における誤差は、Ｎ個分のデータによって触媒特性が変化しても確実に解消されることから、ＮＯｘ吸蔵量の推定精度の向上により一層寄与する。
【発明を実施するための最良の形態】
【００１４】
以下、図面により本発明の実施形態について説明する。
図１は、本発明が適用される車両排気系の概略構成図を示している。同図に示されるように、排気通路２にはＮＯｘ吸蔵触媒４が介装されている。この触媒４は、排気空燃比がストイキオよりもリーン状態のときに排ガス中のＮＯｘを吸蔵し、排気空燃比がリッチ状態にて排ガス中に還元剤としての未燃燃料（ＨＣ）や一酸化炭素（ＣＯ）が存在するときには吸蔵したＮＯｘの放出還元を行う。なお、この触媒４は公知の構造である。
【００１５】
この触媒４の下流側には触媒出口側のＮＯｘ濃度を検出するＮＯｘセンサ６が設けられている。なお、図示はしないが、この排気通路２には、Ｏ２ センサ、温度センサ及びエアフローセンサ（ＡＦＳ）等がそれぞれ設けられており、これら各センサは、排気通路２内の酸素濃度［Ｏ２］、排ガス温度［ｙ］及びＳＶ値（又は排ガス流速）［ｚ］がそれぞれ検出される。
【００１６】
また、上述した各センサはＥＣＵ（電子コントロールユニット）１０に電気的に接続されている。このＥＣＵ１０は入出力装置、記憶装置（ＲＯＭ，ＲＡＭ，不揮発性ＲＡＭ等）、演算装置（ＣＰＵ）及びタイマカウンタ等を備えており、触媒４に対するＮＯｘ吸蔵量の推定の他、エンジンの総合的な制御が実施される。更に、ＥＣＵ１０には種々のマップが設けられている。例えばアクセル開度やエンジン回転速度をパラメータとして、それぞれ触媒入口側のＮＯｘ濃度［ＮＯｘ］、ＣＯ濃度［ＣＯ］及びＨＣ濃度［ＨＣ］が上記マップから読み出される。なお、これらの値は上記マップに代えて、排気通路２に設けた種々のセンサから直接に検出されても良い。
【００１７】
ここで、本実施形態のＥＣＵ１０は吸蔵量推定部１２を備えており、ＮＯｘ吸蔵量が触媒４のＮＯｘ吸蔵特性を反映させた４次元の線形多項式によって算出されている。
より詳しくは、図２に示される如く、吸蔵量推定部１２は推定ＮＯｘ浄化率演算部１４を有し、この演算部１４には式（１）に示される多項式が記憶されている。そして、係数同定部１８にて同定される当該多項式の各係数をｋｉ（ｉ＝０，１，２・・・）とし、触媒４におけるＮＯｘ吸蔵割合［ｘ］が特定されると、既知の排ガス温度［ｙ］及びＳＶ値［ｚ］からＮＯｘ浄化率［ｒＥ］が算出される。
【００１８】
ｒＥ＝ｆ（ｘ，ｙ，ｚ）
＝ｋ０＋ｋ１ｘ＋ｋ２ｙ＋ｋ３ｚ＋ｋ４ｘｙ＋ｋ５ｙｚ
＋ｋ６ｚｘ＋ｋ７ｘ２ｙ＋ｋ８ｘｙ２＋・・    ・・・（１）
このように、触媒４のＮＯｘ吸蔵特性は、実験等により予め知ることができ、多項式で近似可能である。例えば、パラメータが３つ（３次元）であるとすれば、触媒４の特性は図３に示された曲面で表され、この曲面を多項式で表現することができる。なお、式（１）は４次元であるので図３に示された曲面とは異なるものの、この次数は任意に設定可能であって、考え方は３次元の場合と同様である。
【００１９】
次いで、吸蔵量推定部１２は吸蔵量演算部２０を有し、この演算部２０には式（１）を変形した式（２）が記憶されており、係数同定部１８にて同定される各係数ｋｉと排ガス温度［ｙ］及びＳＶ値［ｚ］が入力されると、触媒４におけるＮＯｘ吸蔵割合［ｘ］つまり、現在の吸蔵量／最大の吸蔵量が算出される。
ｘ＝〔ｒＥ−（ｋ０＋ｋ２ｙ＋ｋ３ｚ＋・・）〕／（ｋ１＋ｋ４ｙ＋・・）  ・・・（２）
この式（２）で特定されたＮＯｘ吸蔵割合［ｘ］は、初期段階ではその計算周期における瞬間的なＮＯｘ吸蔵量であるが、図２に示されるように、このＮＯｘ吸蔵割合［ｘ］が演算部１４の式（１）に戻されており、このＮＯｘ吸蔵割合［ｘ］は次回の計算周期におけるＮＯｘ浄化率［ｒＥ］の計算に用いられている。換言すれば、式（１）では前回求めたＮＯｘ吸蔵割合［ｘ］と、今回の計算周期で新たに検出された排ガス温度［ｙ］及びＳＶ値［ｚ］とを用いて、今回の計算周期におけるＮＯｘ浄化率［ｒＥ］を算出している。このように、吸蔵量推定部１２ではＮＯｘ吸蔵割合［ｘ］を式（１）に繰り返し戻してＮＯｘ浄化率［ｒＥ］を算出する。
【００２０】
一方、吸蔵量推定部１２は実ＮＯｘ浄化率演算部１６を有しており、この演算部１６では、ＮＯｘ浄化率［ｒＲ］が式（３）を用いて求められている。
ｒＲ＝（［ＮＯｘ］− 触媒出口側のＮＯｘ濃度 ）／［ＮＯｘ］  ・・・（３）
この触媒出口側のＮＯｘ濃度はＮＯｘセンサ６による検出値であり、［ＮＯｘ］は上記マップから得られた値である。そして、式（１）にて得られたＮＯｘ浄化率［ｒＥ］は推定値（計算値）と定義されるのに対し、式（３）にて得られたＮＯｘ浄化率［ｒＲ］はセンサ値（実測値）と定義される。
【００２１】
ところで、式（１）の多項式は、各係数ｋｉがＮＯｘセンサ６の検出値等に基づいてリーン運転中に逐次補正され、常に最新の触媒状態をモデル化した式として構成されている。
より具体的には、各係数ｋｉ が正確な値であれば式（１）のＮＯｘ浄化率［ｒＥ］と式（３）のＮＯｘ浄化率［ｒＲ］とが一致することを鑑み、上記係数同定部１８では、これらＮＯｘ浄化率の推定値［ｒＥ］と実測値［ｒＲ］とを比較して、推定値［ｒＥ］と実測値［ｒＲ］との間に差が存在する場合には、推定値［ｒＥ］が実測値［ｒＲ］となるように各係数ｋｉ を補正している。この補正は、上述したＮＯｘ吸蔵割合［ｘ］を式（１）に戻してＮＯｘ浄化率［ｒＥ］を算出する度に実施される。つまり、触媒４は使用状況や経年劣化等によってその特性が変化するものであるが、各係数ｋｉ が何度も繰り返して補正されることにより、式（１）の多項式が触媒４の状態を正確に表す式に修正されている。
【００２２】
しかも、本実施形態の係数同定部１８では、ｐがＮＯｘ吸蔵割合［ｘ］や排ガス温度［ｙ］やＳＶ値［ｚ］等の入力値ベクトル[１，ｘ，ｙ，ｚ，ｘｙ，ｙｚ，・・・]、Ｔが転置行列、δがゲイン、及びαがデータの重み付け係数とすると、式（４）の適応推定ロジックを用いて各係数ｋｉを補正している。
【００２３】
【数１】
2007263048.tif 000003
・・・（４）
ここで、左辺のｋｎｅｗは各係数ｋｉの今回の値のベクトル、右辺第１項のｋｏｌｄは各係数ｋｉの前回の値のベクトルであり、右辺第２項中のｉ＝０で示される部分は実測によって収集された収集データに相当する。そして、右辺第２項中のｉ＝１〜Ｎで示される部分は、実測によって収集された領域以外の他の領域に特定されたＮ個の特定されたポイントの基準データに相当する。
【００２４】
そして、本実施形態の当該他の領域とは、触媒４の触媒特性が変化（例えば、劣化）してもその値が変化しないポイント（例えば、排ガス温度［ｙ］が低温時或いは高温時）であり、このポイントのＮＯｘ浄化率は排ガス温度［ｙ］が低温時或いは高温時には略零になる。そこで、係数同定部１８では、上記変化しないポイントを特定ポイントとして用いて不変のＮＯｘ浄化率を得ているのである。
【００２５】
このように、係数同定部１８にて同定された各係数ｋｉは推定ＮＯｘ浄化率演算部１４及び吸蔵量演算部２０に出力され、演算部２０にてＮＯｘ吸蔵割合［ｘ］が算出される。次いで、このＮＯｘ吸蔵割合［ｘ］を式（１）に繰り返し戻してＮＯｘ浄化率［ｒＥ］を算出することにより、触媒４に吸蔵されたＮＯｘの積算値（ＮＯｘ吸蔵量）が高い精度で算出可能となる。
【００２６】
一方、吸蔵量推定部１２は放出量演算部２４を有しており、この演算部２４では、リッチ運転中に式（５）、式（６）を用いてＮＯｘ放出量ｑｄを算出している。リッチ運転中はＮＯｘが触媒４から放出されるため、ＮＯｘ放出量を正確に推定又は検出し、次回のリーン運転時におけるＮＯｘ吸蔵量の初期値を正確に算出するためである。なお、リーン運転からリッチ運転に切り替えられると、係数同定部１８ではリーン運転終了時点の各係数ｋｉ が保持され、吸蔵量演算部２０ではリーン運転終了時点のＮＯｘ吸蔵割合［ｘ］、つまり、ＮＯｘ吸蔵量が記憶されている。
【００２７】
ｑｄ＝∫(［ＣＯ］,［ＨＣ］の濃度×［ｒ′］−０．５×［Ｏ２］)×［ｚ］  ・・・（５）
［ｒ′］＝ｆ（ｙ，ｚ）  ・・・（６）
ここで、［ｒ′］は還元剤利用率であり、図４に示されるマップが設けられている。このマップには、排ガス温度［ｙ］及びＳＶ値［ｚ］に基づく還元剤利用率［ｒ′］の特性が記憶されており、排ガス温度［ｙ］及びＳＶ値［ｚ］から還元剤利用率［ｒ′］が設定される。
【００２８】
次いで、吸蔵量推定部１２は最終吸蔵量演算部２２を有している。この演算部２２では、リッチ運転からリーン運転に再び切り替えられると、リーン運転終了時点のＮＯｘ吸蔵量からＮＯｘ放出量ｑｄを差し引き、リッチ運転終了時点にて触媒４内に残存しているＮＯｘ吸蔵量Ｑａ、すなわち、次回のリーン運転時におけるＮＯｘ吸蔵量の初期値を算出する。
【００２９】
更に、この演算部２２では、リッチ運転終了時点におけるＮＯｘ浄化率の推定値［ｒＥ］と実測値［ｒＲ］とを比較し、推定値［ｒＥ］と実測値［ｒＲ］との差が閾値Ａを超えている場合には、リッチ運転終了時点にて触媒４内に残存しているＮＯｘ吸蔵量が正確ではない旨を判定し、推定値［ｒＥ］が実測値［ｒＲ］となるように、式（２）に基づいてリーン運転時におけるＮＯｘ吸蔵量の初期値を補正している。
【００３０】
図５を参照すると、吸蔵量推定部１２によるＮＯｘ吸蔵量推定のフローチャートが示されている。同図のステップＳ５０１では、リーン運転が終了するまでは、式（１）によるＮＯｘ浄化率の推定値［ｒＥ］が式（３）によるＮＯｘ浄化率の実測値［ｒＲ］に一致するように、式（４）による基準データを同時使用した適応推定ロジックを用いて各係数ｋｉ を補正し、式（２）によるＮＯｘ吸蔵割合 [ｘ]を繰り返し算出して触媒４に吸蔵されたＮＯｘの積算値を求め、ステップＳ５０２に進む。
【００３１】
ステップＳ５０２では、リッチ運転が開始されるか否かが判別され、リッチ運転が開始、すなわち、ＹＥＳと判定された場合にはステップＳ５０３に進み、このリッチ運転が終了するまでは、リーン運転終了時点の各係数ｋｉ やＮＯｘ吸蔵量がそれぞれ保持されるとともに、式（５）や式（６）によるＮＯｘ放出量ｑｄを算出する。
続いて、ステップＳ５０４では、リッチ運転が終了するか否かが判別され、リッチ運転が終了、すなわち、ＹＥＳと判定された場合にはステップＳ５０５に進む。このステップＳ５０５では、リッチ運転終了時点にて触媒４内に残存しているＮＯｘ吸蔵量が正確であるか否かが判別され、リッチ運転終了時点におけるＮＯｘ浄化率の推定値［ｒＥ］と実測値［ｒＲ］との差が閾値Ａを超えている、すなわち、ＹＥＳと判定された場合にはステップＳ５０６に進み、リーン運転時におけるＮＯｘ吸蔵量の初期値を補正して一連のルーチンを抜ける。
【００３２】
一方、ステップＳ５０５にて、リッチ運転終了時点におけるＮＯｘ浄化率の推定値［ｒＥ］と実測値［ｒＲ］との差が閾値Ａを超えていない場合には、ステップＳ５０１に戻り、上述した各係数ｋｉ の補正等が順次実施される。
以上のように、本発明によれば、ＮＯｘ吸蔵特性を反映させた式（１）の多項式を用いてＮＯｘ吸蔵量を推定するにあたり、この多項式の各係数ｋｉが、実測によって収集されたＮＯｘ浄化率［ｒＲ］の他、触媒特性が劣化してもその値が変化しないポイント、例えば、排ガス温度［ｙ］が低温時或いは高温時のポイントのＮＯｘ浄化率（＝略零）に基づいて逐次補正されている。
【００３３】
この点につき詳述すると、まず、図６にはエンジンの運転状態に対応する排ガス温度の範囲に対し、触媒４の劣化前のＮＯｘ浄化率の実測値（図中、一点鎖線で示す：真値）と、劣化後のＮＯｘ浄化率の実測値（図中、点線で示す：真値）とがそれぞれ示されており、触媒４の劣化に伴って山型の形状が内側に向けて小さくなることが分かる。
ここで、同図（ａ）に示された劣化後のＮＯｘ浄化率の推定値（図中、実線で示す）は、実測による収集データ（図中●で示す）のみを用いて多項式の各係数ｋｉが補正されており、実線で示された当該推定値の曲線は、収集された運転範囲が排ガス温度［ｙ］が中温程度の領域に集中していた場合、その領域へのフィッティングに支配され、その他の運転範囲（排ガス温度［ｙ］が低温及び高温）では上記点線で示された実測値との差が非常に大きくなってしまう。
【００３４】
これに対し、本発明では同図（ｂ）に示されるように、実測による収集データが中温程度の排ガス温度の領域に集中したとしても、ＮＯｘ浄化率が略零になる低温や高温の排ガス温度の上記各実測値に基準データとして用いるためのポイントが設定されている（図中、更に●で示す）。よって、実線で示された当該推定値の曲線は、収集された運転範囲（排ガス温度［ｙ］が中温程度）では、点線で示された実測値にほぼ一致するし、その他の運転範囲（排ガス温度［ｙ］が低温及び高温）でもまた、収集された運転範囲へのフィッティングに偏ることなく、上記点線で示された実測値にほぼ一致する。
【００３５】
このように、本発明による基準データを同時使用した適応推定ロジックによれば、多項式の各係数ｋｉが大域的なデータに基づいた値になり、仮にエンジンの運転状態が収集された運転範囲（排ガス温度［ｙ］が中温程度）から他の運転範囲（排ガス温度［ｙ］が低温や高温）に移行した場合にも、移行前に実施された他の運転範囲の演算結果に高い信頼性が得られるので、各係数の同定計算が安定するし、その計算精度も向上する。
【００３６】
これらの結果、多項式によるモデルの推定誤差が最小限度に抑えられ、触媒モデルが収集された特定の運転範囲へのフィッティングに支配され難くなり、ＮＯｘ吸蔵量の推定精度の向上に寄与する。ひいては、リッチスパイクの実施時期はＮＯｘ吸蔵量に基づいて判断されることから、このリッチスパイクの間隔の最適化が図られ、ＮＯｘ浄化率の向上や燃料消費の抑制が達成可能となる。
【００３７】
しかも、触媒４の状態をモニタする精度が向上するので、Ｓパージの時期を決定する機能や、触媒の劣化を運転者に警告する機能等を車両に備えることが可能になる。
以上で本発明の一実施形態についての説明を終えるが、本発明は上記実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々の変更ができるものである。
例えば、上記実施形態では、２つのポイントを設定しているが、このポイントの数は任意に設定可能である。また、ＮＯｘ放出量ｑｄは図４に示されたマップから還元剤利用率［ｒ′］を求めているが、触媒入口側の還元剤濃度［ｘ′］、排ガス温度［ｙ］及びＳＶ値［ｚ］を用い、触媒４のＮＯｘ放出特性を反映させた多項式から求めても良い（式（７））。
【００３８】
ｒ′＝ｆ（ｘ′，ｙ，ｚ）
＝ｍ０＋ｍ１ｘ′＋ｍ２ｙ＋ｍ３ｚ＋ｍ４ｘ′ｙ＋ｍ５ｙｚ
＋ｍ６ｚｘ′＋ｍ７ｘ′２ｙ＋ｍ８ｘ′ｙ２＋・・    ・・・（７）
ここで、ｍｉ （ｉ＝１，２，・・・）は係数である。
この場合には少なくとも排気通路２に還元剤の濃度を検出するセンサを設け、このセンサから得られた還元剤濃度の実測値［ｘＲ′］と、式（７）から得られた還元剤濃度の推定値［ｘＥ′］とから各係数ｍｉを逐次補正して還元剤利用率ｒ′を求め、ここで求めた還元剤利用率ｒ′を上述した式（５）に代入することでＮＯｘ放出量ｑｄを求めることができる。これにより、ＮＯｘ放出量の推定精度が高められるとの効果を奏する。
【図面の簡単な説明】
【００３９】
【図１】本発明の一実施形態に係るＮＯｘ吸蔵量の推定方法を説明する排気系の概略構成図である。
【図２】図１のＥＣＵに構成された吸蔵量推定部のブロック図である。
【図３】図１のＮＯｘ吸蔵触媒によるＮＯｘ吸蔵特性を示すマップである。
【図４】図１のＮＯｘ吸蔵触媒によるＮＯｘ放出特性を示すマップである。
【図５】図１の吸蔵量推定部によるフローチャートである。
【図６】（ａ）は従来技術による実測ＮＯｘ浄化率と排ガス温度との関係図であり、（ｂ）は図１の推定方法による実測ＮＯｘ浄化率および推定ＮＯｘ浄化率と排ガス温度との関係図であり、本発明の作用効果を説明する図である。
【符号の説明】
【００４０】
２  排気通路
４  ＮＯｘ吸蔵触媒
１０  ＥＣＵ（電子コントロールユニット）
１２  吸蔵量推定部
１８  係数同定部
【図１】
2007263048.tif 000004
【図２】
2007263048.tif 000005
【図３】
2007263048.tif 000006
【図４】
2007263048.tif 000007
【図５】
2007263048.tif 000008
【図６】
2007263048.tif 000009