(57)【特許請求の範囲】
【請求項１】
流入する排気ガスの空燃比がリーンであるときにはＮＯXを吸収し、流入する排気ガスの空燃比がリッチになると吸収したＮＯXを排気ガス中に含まれる還元剤によって放出し還元するＮＯX吸収剤を機関排気通路内に配置し、該ＮＯX吸収剤下流の機関排気通路内に排気ガス中のＮＯX濃度を検出しうるセンサを配置した内燃機関の排気浄化装置であって、前記センサの出力を利用してＮＯX吸収剤に吸収しうると推定される最大ＮＯX吸収量の推定値と、ＮＯX吸収剤に流入する排気ガスの空燃比がリーンであるときにＮＯX吸収剤にて達成されると推定されるＮＯX吸収剤へのＮＯX吸収速度の推定値との少なくとも一方を算出し、斯くして算出された最大ＮＯX吸収量の推定値またはＮＯX吸収速度の推定値が判定値よりも小さいときにはＮＯX吸収剤が劣化したと判定するようにした内燃機関の排気浄化装置において、ＮＯX吸収剤下流に流出するＮＯXの量と、ＮＯX吸収剤の最大ＮＯX吸収量と、ＮＯX吸収剤におけるＮＯX吸収速度との間に成り立つ関係式を予め求めておき、ＮＯX吸収剤下流に流出するＮＯXの量を上記センサの出力に基づいて算出し、斯くして算出されたＮＯXの量を上記関係式に代入することにより最大ＮＯX吸収量の推定値と、ＮＯX吸収速度の推定値とを算出するようにした内燃機関の排気浄化装置。
【請求項２】
ＮＯX吸収剤に流入する排気ガスの空燃比がリーンであるときに機関運転状態に基づいて算出されるＮＯX吸収剤に吸収されるＮＯXの量を積算することによりＮＯX吸収剤に吸収されているＮＯXの量を推定するＮＯX吸収量推定手段を具備し、該ＮＯX吸収量推定手段により推定されたＮＯX吸収量が許容値を超えたときにＮＯX吸収剤に流入する排気ガスの空燃比をリーンからリッチに切換えるようにした請求項１に記載の内燃機関の排気浄化装置。
【請求項３】
上記ＮＯX吸収量推定手段がＮＯX吸収剤に吸収されているＮＯXの量を推定するときに機関運転状態として最大ＮＯX吸収量とＮＯX吸収速度とを用いる請求項２に記載の内燃機関の排気浄化装置。
【請求項４】
流入する排気ガスの空燃比がリーンであるときにはＮＯXを吸収し、流入する排気ガスの空燃比がリッチになると吸収したＮＯXを排気ガス中に含まれる還元剤によって放出し還元するＮＯX吸収剤を機関排気通路に配置し、該ＮＯX吸収剤下流の機関排気通路内に排気ガス中のアンモニア濃度を検出しうるセンサを配置し、ＮＯX吸収剤に流入する排気ガスの空燃比がリーンからリッチにされたときにＮＯX吸収剤からＮＯXを放出し還元するために使用されなかった余剰の還元剤がアンモニアの形でＮＯX吸収剤から流出し、上記センサによりＮＯX吸収剤から流出する排気ガス中のアンモニア濃度が判定値を超えたときにＮＯX吸収剤からのＮＯXの放出・還元が完了したと判定するようにした内燃機関の排気浄化装置において、ＮＯX吸収剤に流入する排気ガスの空燃比がリーンであるときには機関運転状態に基づいて算出されるＮＯX吸収剤に吸収されるＮＯXの量を積算し、ＮＯX吸収剤に流入する排気ガスの空燃比がリッチであるときにはＮＯX吸収剤におけるＮＯXの還元速度に基づいて算出されるＮＯX吸収剤において還元されるＮＯXの量を減算することによりＮＯX吸収剤に吸収されているＮＯXの量を推定するＮＯX吸収量推定手段を具備し、ＮＯX吸収剤からのＮＯXの放出・還元が完了したと判定されたときに上記ＮＯX吸収量推定手段により推定されるＮＯX吸収量が基準値からずれているときには上記ＮＯX吸収量推定手段にて使用されるＮＯX還元速度を補正し、斯くして補正されたＮＯX還元速度が判定値よりも遅いときにはＮＯX吸収剤が劣化したと判定するようにした内燃機関の排気浄化装置。
【請求項５】
ＮＯX吸収剤からのＮＯXの放出・還元が完了したと判定されたときにＮＯX吸収量推定手段により推定されるＮＯX吸収量が基準値よりも多いときにはＮＯX吸収量推定手段にて使用されるＮＯX還元速度を速くするように補正し、基準値よりも少ないときには該ＮＯX還元速度を遅くするように補正する請求項４に記載の内燃機関の排気浄化装置。
【請求項６】
上記ＮＯX吸収量推定手段がＮＯX吸収剤に吸収されているＮＯXの量を推定するときに機関運転状態としてＮＯX吸収剤の最大ＮＯX吸収量とＮＯX吸収剤へのＮＯXの吸収速度とを用いる請求項４に記載の内燃機関の排気浄化装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は内燃機関の排気浄化装置に関する。
【０００２】
【従来の技術】
内燃機関の燃焼室から排出される排気ガス中の窒素酸化物（ＮＯX ）を浄化するためのＮＯX 浄化触媒を機関排気通路内に配置した排気浄化装置が公知である。ＮＯX 浄化触媒としては流入する排気ガスの空燃比がリーンであるときにはＮＯX を吸収し、流入する排気ガスの空燃比がリッチになると吸収したＮＯX を排気ガス中に含まれる還元剤によって放出し還元するＮＯX 吸収剤が知られている。このようなＮＯX 吸収剤においては吸収することができるＮＯX の量、すなわち最大ＮＯX 吸収量に限界がある。そしてＮＯX 吸収剤には吸収されているＮＯX の量がこの最大ＮＯX 吸収量を越えるとＮＯX 吸収剤はもはやＮＯX を吸収することはできず、このためＮＯX 吸収剤下流へとＮＯX が流出してしまう。そこでＮＯX 吸収剤に吸収されているＮＯX の量が最大ＮＯX 吸収量を越える前にＮＯX 吸収剤に流入する排気ガスの空燃比をリーンからリッチに切換えてＮＯX 吸収剤からＮＯX を放出し還元する必要がある。しかしながら最大ＮＯX 吸収量が非常に少なくなったときには還元剤として炭化水素（燃料）を頻繁にＮＯX 吸収剤に供給しなければならず、燃費が悪化する。この場合において燃費の悪化を防止するためには何らかの処置をＮＯX 吸収剤に施こさなければならないが、何らかの処置を施すべきか否かを判断するためにはＮＯX 吸収剤が劣化しているか否かを判定する必要がある。
【０００３】
ＮＯX 吸収剤の劣化を判定するシステムとしては特開２０００−１０４５３６号公報に開示されたシステムが公知である。当該公報によればＮＯX 吸蔵触媒下流に排気ガス中のＮＯX 濃度を検出することができるＮＯX センサを配置し、該ＮＯX センサにより検出されるＮＯX 濃度が所定値に達したときにＮＯX 吸蔵触媒に流入する排気ガスの空燃比をリーンからリッチに切換えるようにしたシステムにおいて、排気ガスの空燃比をリーンからリッチに切換える処理を実行する間隔が短かくなったことをもってＮＯX 吸蔵触媒が劣化したと判定するようにしている。
【０００４】
【発明が解決しようとする課題】
ところで上述した特開２０００−１０４５３６号公報によれば排気ガスの空燃比をリーンからリッチに切換える処理を実行する間隔に基づいてＮＯX 吸蔵触媒の劣化の有無を判定しているがこの排気ガスの空燃比をリーンからリッチに切換える処理の実行間隔は例えばＮＯX 吸蔵触媒に流入するＮＯX の量が多くなればたとえＮＯX 吸蔵触媒が劣化していなくても短くなる。したがって当該公報のシステムではＮＯX 吸蔵触媒の劣化を正確に判定することはできない。こうした事情に鑑み本発明の目的はＮＯX 吸収剤の劣化を正確に判定することにある。
【０００５】
【課題を解決するための手段】
上記課題を解決するために１番目の発明では流入する排気ガスの空燃比がリーンであるときにはＮＯXを吸収し、流入する排気ガスの空燃比がリッチになると吸収したＮＯXを排気ガス中に含まれる還元剤によって放出し還元するＮＯX吸収剤を機関排気通路内に配置し、該ＮＯX吸収剤下流の機関排気通路内に排気ガス中のＮＯX濃度を検出しうるセンサを配置した内燃機関の排気浄化装置であって、前記センサの出力を利用してＮＯX吸収剤に吸収しうると推定される最大ＮＯX吸収量の推定値と、ＮＯX吸収剤に流入する排気ガスの空燃比がリーンであるときにＮＯX吸収剤にて達成されると推定されるＮＯX吸収剤へのＮＯX吸収速度の推定値との少なくとも一方を算出し、斯くして算出された最大ＮＯX吸収量の推定値またはＮＯX吸収速度の推定値が判定値よりも小さいときにはＮＯX吸収剤が劣化したと判定するようにした内燃機関の排気浄化装置において、ＮＯX吸収剤下流に流出するＮＯXの量と、ＮＯX吸収剤の最大ＮＯX吸収量と、ＮＯX吸収剤におけるＮＯX吸収速度との間に成り立つ関係式を予め求めておき、ＮＯX吸収剤下流に流出するＮＯXの量を上記センサの出力に基づいて算出し、斯くして算出されたＮＯXの量を上記関係式に代入することにより最大ＮＯX吸収量の推定値と、ＮＯX吸収速度の推定値とを算出する。すなわち排気ガス中のＮＯX濃度を検出することができるセンサの出力を利用して算出されるＮＯX吸収剤の最大ＮＯX吸収量またはＮＯX吸収剤へのＮＯX吸収速度に基づいてＮＯX吸収剤の劣化が判定される。
【０００７】
２番目の発明では１番目の発明においてＮＯX吸収剤に流入する排気ガスの空燃比がリーンであるときに機関運転状態に基づいて算出されるＮＯX吸収剤に吸収されるＮＯXの量を積算することによりＮＯX吸収剤に吸収されているＮＯXの量を推定するＮＯX吸収量推定手段を具備し、該ＮＯX吸収量推定手段により推定されたＮＯX吸収量が許容値を超えたときにＮＯX吸収剤に流入する排気ガスの空燃比をリーンからリッチに切換える。
【０００８】
３番目の発明では２番目の発明において上記ＮＯX吸収量推定手段がＮＯX吸収剤に吸収されているＮＯXの量を推定するときに機関運転状態として最大ＮＯX吸収量とＮＯX吸収速度とを用いる。
上記課題を解決するために４番目の発明では流入する排気ガスの空燃比がリーンであるときにはＮＯXを吸収し、流入する排気ガスの空燃比がリッチになると吸収したＮＯXを排気ガス中に含まれる還元剤によって放出し還元するＮＯX吸収剤を機関排気通路に配置し、該ＮＯX吸収剤下流の機関排気通路内に排気ガス中のアンモニア濃度を検出しうるセンサを配置し、ＮＯX吸収剤に流入する排気ガスの空燃比がリーンからリッチにされたときにＮＯX吸収剤からＮＯXを放出し還元するために使用されなかった余剰の還元剤がアンモニアの形でＮＯX吸収剤から流出し、上記センサによりＮＯX吸収剤から流出する排気ガス中のアンモニア濃度が判定値を超えたときにＮＯX吸収剤からのＮＯXの放出・還元が完了したと判定するようにした内燃機関の排気浄化装置において、ＮＯX吸収剤に流入する排気ガスの空燃比がリーンであるときには機関運転状態に基づいて算出されるＮＯX吸収剤に吸収されるＮＯXの量を積算し、ＮＯX吸収剤に流入する排気ガスの空燃比がリッチであるときにはＮＯX 吸収剤におけるＮＯX の還元速度に基づいて算出されるＮＯX吸収剤において還元されるＮＯXの量を減算することによりＮＯX吸収剤に吸収されているＮＯXの量を推定するＮＯX吸収量推定手段を具備し、ＮＯX吸収剤からのＮＯXの放出・還元が完了したと判定されたときに上記ＮＯX吸収量推定手段により推定されるＮＯX吸収量が基準値からずれているときには上記ＮＯX吸収量推定手段にて使用されるＮＯX還元速度を補正し、斯くして補正されたＮＯX還元速度が判定値よりも遅いときにはＮＯX吸収剤が劣化したと判定する。すなわちアンモニア濃度を検出しうるセンサの出力に基づいてＮＯX吸収剤におけるＮＯXの放出・還元の完了を判定し、このときにＮＯX吸収量推定手段により推定されるＮＯX吸収量に基づいて補正されたＮＯX還元速度に基づいてＮＯX吸収剤の劣化が判定される。
【０００９】
５番目の発明では４番目の発明においてＮＯX吸収剤からのＮＯXの放出・還元が完了したと判定されたときにＮＯX吸収量推定手段により推定されるＮＯX吸収量が基準値よりも多いときにはＮＯX吸収量推定手段にて使用されるＮＯX還元速度を速くするように補正し、基準値よりも少ないときには該ＮＯX還元速度を遅くするように補正する。
【００１０】
６番目の発明では４番目の発明において上記ＮＯX吸収量推定手段がＮＯX吸収剤に吸収されているＮＯXの量を推定するときに機関運転状態としてＮＯX吸収剤の最大ＮＯX吸収量とＮＯX吸収剤へのＮＯXの吸収速度とを用いる。
【００１１】
【発明の実施の形態】
図１は本発明を筒内噴射式火花点火機関に適用した場合を示している。しかしながら本発明は圧縮着火式内燃機関に適用することもできる。図１を参照すると１は機関本体、２はシリンダブロック、３はシリンダブロック２内で往復動するピストン、４はシリンダブロック２上に固定されたシリンダヘッド、５はピストン３とシリンダヘッド４との間に形成された燃焼室、６は吸気弁、７は吸気ポート、８は排気弁、９は排気ポートを夫々示す。図１に示したようにシリンダヘッド４の内壁面の中央部には点火栓１０が配置され、シリンダヘッド４内壁面周辺部には燃料噴射弁１１が配置される。またピストン３の頂面上には燃料噴射弁１１の下方から点火栓１０の下方まで延びるキャビティ１２が形成されている。
【００１２】
各気筒の吸気ポート７は夫々対応する吸気枝管１３を介してサージタンク１４に連結され、サージタンク１４は吸気ダクト１５およびエアフロメータ１６を介してエアクリーナ（図示せず）に連結される。吸気ダクト１５内にはステップモータ１７によって駆動されるスロットル弁１８が配置される。一方、各気筒の排気ポート９は排気マニホルド１９に連結され、この排気マニホルド１９は酸化触媒または三元触媒２０を内臓した触媒コンバータ２１および排気管２２を介してＮＯX 吸収剤２３を内臓したケーシング２４に連結される。排気マニホルド１９とサージタンク１４とは再循環排気ガス（以下、ＥＧＲガスという）導管２６を介して互いに連結され、このＥＧＲガス導管２６内にはＥＧＲガス制御弁２７が配置される。
【００１３】
電子制御ユニット３１はディジタルコンピュータからなり、双方向性バス３２を介して相互に接続されたＲＡＭ（ランダムアクセスメモリ）３３、ＲＯＭ（リードオンリメモリ）３４、ＣＰＵ（マイクロプロセッサ）３５、入力ポート３６および出力ポート３７を具備する。エアフロメータ１６は吸入空気量に比例した出力電圧を発生し、この出力電圧が対応するＡＤ変換器３８を介して入力ポート３６に入力される。排気マニホルド１９には空燃比を検出するための空燃比センサ２８が取付けられ、この空燃比センサ２８の出力信号が対応するＡＤ変換器３８を介して入力ポート３６に入力される。またＮＯX 吸収剤２３を内臓したケーシング２４の出口に接続された排気管２５内には排気ガス中のＮＯX 濃度およびアンモニア濃度を共に検出可能なＮＯX アンモニアセンサ２９と、空燃比センサ３０とが配置され、これらＮＯX アンモニアセンサ２９および空燃比センサ３０の出力信号が対応するＡＤ変換器３８を介して入力ポート３６に入力される。
【００１４】
またアクセルペダル４０にはアクセルペダル４０の踏込み量に比例した出力電圧を発生する負荷センサ４１が接続され、負荷センサ４１の出力電圧は対応するＡＤ変換器３８を介して入力ポート３６に入力される。クランク角センサ４２は例えばクランクシャフトが３０度回転する毎に出力パルスを発生し、この出力パルスが入力ポート３６に入力される。ＣＰＵ３５ではこのクランク角センサ４２の出力パルスから機関回転数が計算される。一方、出力ポート３７は対応する駆動回路３９を介して点火栓１０、燃料噴射弁１１、ステップモータ１７およびＥＧＲ制御弁２７に接続される。
【００１５】
次に図２を参照しつつ図１に示したＮＯX アンモニアセンサ２９のセンサ部の構造について簡単に説明する。図２を参照するとＮＯX アンモニアセンサ２９のセンサ部は互いに積層された６つの酸化ジルコニア等の酸素イオン伝導性固体電解質層からなり、これら６つの固体電解質層を以下、上から順に第１層Ｌ1 、第２層Ｌ2 、第３層Ｌ3 、第４層Ｌ4 、第５層Ｌ5 、第６層Ｌ6 と称する。
【００１６】
図２を参照すると第１層Ｌ1 と第３層Ｌ3 との間に例えば多孔質のまたは細孔が形成されている第１の拡散律速部材５０と第２の拡散律速部材５１とが配置されており、これら拡散律速部材５０，５１間には第１室５２が形成され、第２の拡散律速部材５１と第２層Ｌ2 との間には第２室５３が形成されている。また第３層Ｌ3 と第５層Ｌ5 との間には外気に連通している大気室５４が形成されている。一方、第１の拡散律速部材５０の外端面は排気ガスと接触している。したがって排気ガスは第１の拡散律速部材５０を介して第１室５２内に流入し、斯くして第１室５２内は排気ガスで満たされている。
【００１７】
一方、第１室５２に面する第１層Ｌ1 の内周面上には陰極側第１ポンプ電極５５が形成され、第１層Ｌ1 の外周面上には陽極側第１ポンプ電極５６が形成されており、これら第１ポンプ電極５５，５６間には第１ポンプ電圧源５７により電圧が印加される。第１ポンプ電極５５，５６間に電圧が印加されると第１室５２内の排気ガス中に含まれる酸素が陰極側第１ポンプ電極５５と接触して酸素イオンとなり、この酸素イオンは第１層Ｌ1 内を陽極側第１ポンプ電極５６に向けて流れる。したがって第１室５２内の排気ガス中に含まれる酸素は第１層Ｌ1 内を移動して外部に汲み出されることになり、このとき外部に汲み出される酸素量は第１ポンプ電圧源５７の電圧が高くなるほど多くなる。
【００１８】
一方、大気室５４に面する第３層Ｌ3 の内周面上には基準電極５８が形成されている。ところで酸素イオン伝導性固体電解質では固体電解質層の両側において酸素濃度に差があると酸素濃度の高い側から酸素濃度の低い側に向けて固体電解質層内を酸素イオンが移動する。図２に示した例では大気室５４内の酸素濃度の方が第１室５２内の酸素濃度よりも高いので大気室５４内の酸素は基準電極５８と接触することにより電荷を受け取って酸素イオンとなり、この酸素イオンは第３層Ｌ3 、第２層Ｌ2 および第１層Ｌ1 内を移動し、陰極側第１ポンプ電極５５において電荷を放出する。その結果、基準電極５８と陰極側第１ポンプ電極５５との間に符号５９で示した電圧Ｖ0 が発生する。この電圧Ｖ0 は大気圧室５４内と第１室５２内との酸素濃度差に比例する。
【００１９】
図２に示した例ではこの電圧Ｖ0 が第１室５２内の酸素濃度が１p.p.m.のときに生ずる電圧に一致するように第１ポンプ電圧源５７の電圧がフィードバック制御される。すなわち第１室５２内の酸素は第１室５２内の酸素濃度が１p.p.m.となるように第１層Ｌ1 を通って汲み出され、それによって第１室５２内の酸素濃度が１p.p.m.に維持される。
【００２０】
なお陰極側第１ポンプ電極５５はＮＯX に対しては還元性の低い材料、例えば金Ａｕと白金Ｐｔとの合金から形成されており、したがって排気ガス中に含まれるＮＯX は第１室５２内ではほとんど還元されない。したがってこのＮＯX は第２の拡散律速部材５１を通って第２室５３内に流入する。
一方、第２室５３に面する第１層Ｌ1 の内周面上には陰極側第２ポンプ電極６０が形成されており、この陰極側第２ポンプ電極６０と陽極側第１ポンプ電極５６との間には第２ポンプ電圧源６１により電圧が印加される。これらポンプ電極６０，５６間に電圧が印加されると第２室５３内の排気ガス中に含まれる酸素が陰極側第２ポンプ電極６０と接触して酸素イオンとなり、この酸素イオンは第１層Ｌ1 内を陽極側第１ポンプ電極５６に向けて流れる。したがって第２室５３内の排気ガス中に含まれる酸素は第１層Ｌ1 内を移動して外部に汲み出されることになり、このとき外部に汲み出される酸素量は第２ポンプ電圧源６１の電圧が高くなるほど多くなる。
【００２１】
一方、前述したように酸素イオン伝導性固体電解質では固体電解質層の両側において酸素濃度に差があると酸素濃度の高い側から酸素濃度の低い側に向けて固体電解質層内を酸素イオンが移動する。図２に示した例では大気室５４内の酸素濃度の方が第２室５３内の酸素濃度よりも高いので大気室５４内の酸素は基準電極５８と接触することにより電荷を受け取って酸素イオンとなり、この酸素イオンは第３層Ｌ3 、第２層Ｌ2 および第１層Ｌ1 内を移動し、陰極側第２ポンプ電極６０において電荷を放出する。その結果、基準電極５８と陰極側第２ポンプ電極６０との間に符号６２で示した電圧Ｖ1 が発生する。この電圧Ｖ1 は大気圧室５４内と第２室５３内との酸素濃度差に比例する。
【００２２】
図２に示した例ではこの電圧Ｖ1 が第２室５３内の酸素濃度が０．０１p.p.m.のときに生ずる電圧に一致するように第２ポンプ電圧源６１の電圧がフィードバック制御される。すなわち第２室５３内の酸素は第２室５３内の酸素濃度が０．０１p.p.m.となるように第１層Ｌ1 を通って汲み出され、それによって第２室５３内の酸素濃度が０．０１p.p.m.に維持される。
【００２３】
なお陰極側第２ポンプ電極６０もＮＯX に対しては還元性の低い材料、例えば金Ａｕと白金Ｐｔとの合金から形成されており、したがって排気ガス中に含まれるＮＯX は陰極側第２ポンプ電極６０と接触してもほとんど還元されない。
一方、第２室５３に面する第３層Ｌ3 の内周面上にはＮＯX 検出用の陰極側ポンプ電極６３が形成されている。この陰極側ポンプ電極６３はＮＯX に対して強い還元性を有する材料、例えばロジウムＲｈや白金Ｐｔから形成されている。したがって第２室５３内のＮＯX 、実際には大部分を占めるＮＯが陰極側ポンプ電極６３上においてＮ2 とＯ2 とに分解される。図２に示したようにこの陰極側ポンプ電極６３と基準電極５８との間には一定電圧６４が印加されており、したがって陰極側ポンプ電極６３上において分解生成されたＯ2 は酸素イオンとなって第３層Ｌ3 内を基準電極５８に向けて移動する。このとき陰極側ポンプ電極６３と基準電極５８との間にはこの酸素イオン量に比例した符号６５で示した電流Ｉ1 が流れる。
【００２４】
前述したように第１室５２内ではＮＯX はほとんど還元されず、また第２室５３内には酸素はほとんど存在しない。したがって電流Ｉ1 は排気ガス中に含まれるＮＯX 濃度に比例することになり、斯くして電流Ｉ1 から排気ガス中のＮＯX 濃度を検出できることになる。
一方、排気ガス中に含まれるアンモニアＮＨ3 は第１室５２内においてＮＯとＨ2 Ｏとに分解され（４ＮＨ3 ＋５Ｏ2 →４ＮＯ＋６Ｈ2 Ｏ）、この分解されたＮＯは第２の拡散律速部材５１を通って第２室５３内に流入する。このＮＯは陰極側ポンプ電極６３上においてＮ2 とＯ2 とに分解され、分解生成されたＯ2 は酸素イオンとなって第３層Ｌ3 内を基準電極５８に向けて移動する。このときにも電流Ｉ1 は排気ガス中に含まれるＮＨ3 濃度に比例し、斯くして電流Ｉ1 から排気ガス中のＮＨ3 濃度を検出できることになる。
【００２５】
図３は電流Ｉ1 と排気ガス中のＮＯX 濃度およびＮＨ3 濃度との関係を示している。図３から電流Ｉ1 は排気ガス中のＮＯX 濃度およびＮＨ3 濃度に比例していることがわかる。
一方、排気ガス中の酸素濃度が高いほど、すなわち空燃比がリーンであるほど第１室５２から外部に汲み出される酸素量が多くなり、符号６６で示した電流Ｉ2 が増大する。したがってこの電流Ｉ2 から排気ガスの空燃比を検出することができる。
【００２６】
なお第５層Ｌ5 と第６層Ｌ6 との間にはＮＯX アンモニアセンサ２９のセンサ部を加熱するための電気ヒータ６７が配置されており、この電気ヒータ６７によってＮＯX アンモニアセンサ２９のセンサ部は７００℃から８００℃に加熱される。
図４はＮＯX 吸収剤２３下流の排気管２５内に配置された空燃比センサ３０の出力電圧Ｅ（Ｖ）、すなわち一般的な表現を用いると空燃比検出手段の出力信号レベルを示している。図４からわかるように空燃比センサ３０は排気ガスの空燃比がリッチのときには０．９（Ｖ）程度の出力電圧を発生し、排気ガスの空燃比がリーンのときには０．１（Ｖ）程度の出力電圧を発生する。すなわち図４に示した例ではリッチであることを示す出力信号レベルは０．９（Ｖ）であり、リーンであることを示す出力信号レベルは０．１（Ｖ）である。
【００２７】
一方、前述したようにＮＯX アンモニアセンサ２９の電流Ｉ2 から排気ガスの空燃比を検出することができ、したがって空燃比検出手段としてＮＯX アンモニアセンサ２９を用いることもできる。この場合には空燃比センサ３０を設ける必要がない。
次に図５（Ａ）を参照しつつ図１に示した内燃機関の燃料噴射制御について説明する。なお図５（Ａ）において縦軸は機関負荷Ｑ／Ｎ（吸入空気量Ｑ／機関回転数Ｎ）を表しており、横軸は機関回転数Ｎを表している。
【００２８】
図５（Ａ）において実線Ｘ1 よりも低負荷側の運転領域では成層燃焼が行われる。すなわちこのときには図１に示したように圧縮行程末期に燃料噴射弁１１からキャビティ１２内に向けて燃料Ｆが噴射される。この燃料はキャビティ１２の内周面により案内されて点火栓１０周りに混合気を形成し、この混合気が点火栓１０によって着火燃焼せしめられる。このとき燃焼室５内における平均空燃比はリーンとなっている。
【００２９】
一方、図５（Ａ）において実線Ｘ1 よりも高負荷側の領域では吸気行程中に燃料噴射弁１１から燃料が噴射され、このときには均一混合気燃焼が行われる。なお実線Ｘ1 と鎖線Ｘ2 との間ではリーン空燃比のもとで均一混合気燃焼が行われ、鎖線Ｘ2 と鎖線Ｘ3 との間では理論空燃比のもとで均一混合気燃焼が行われ、鎖線Ｘ3 よりも高負荷側ではリッチ空燃比のもとで均一混合気燃焼が行われる。
【００３０】
本発明では空燃比を理論空燃比とするのに必要な基本燃料噴射量ＴＡＵが図５（Ｂ）に示したように機関負荷Ｑ／Ｎおよび機関回転数Ｎの関数としてマップの形で予めＲＯＭ３４内に記憶されており、基本的にはこの基本燃料噴射量ＴＡＵに補正係数Ｋを乗算することによって（場合によっては後述するように修正係数ＫS をさらに乗算することによって）最終的な燃料噴射量ＴＡＵＯ（＝Ｋ・ＴＡＵ）が算出される。この補正係数Ｋは図５（Ｃ）に示したように機関負荷Ｑ／Ｎおよび機関回転数Ｎの関数としてマップの形で予めＲＯＭ３４内に記憶されている。
【００３１】
この補正係数Ｋの値はリーン空燃比のもとで燃焼が行われる図５（Ａ）の鎖線Ｘ2 よりも低負荷側の運転領域では１．０よりも小さく、リッチ空燃比のもとで燃焼が行われる図５（Ａ）の鎖線Ｘ3 よりも高負荷側の運転領域では１．０よりも大きくなる。またこの補正係数Ｋは鎖線Ｘ2 と鎖線Ｘ3 との間の運転領域では１．０とされ、このとき空燃比は理論空燃比となるように空燃比センサ２８の出力信号に基づいてフィードバック制御される。
【００３２】
機関排気通路内に配置されたＮＯX 吸収剤２３は例えばアルミナを担体とし、この担体上に例えばカリウムＫ、ナトリウムＮａ、リチウムＬｉ、セシウムＣｓのようなアルカリ金属、バリウムＢａ、カルシウムＣａのようなアルカリ土類、ランタンＬａ、イットリウムＹのような希土類から選ばれた少くとも一つと、白金Ｐｔのような貴金属とが担持されている。この場合、ケーシング２４内に例えばコージライトからなるパティキュレートフィルタを配置し、このパティキュレートフィルタ上にアルミナを担体とするＮＯX 吸収剤２３を担持させることもできる。
【００３３】
いずれの場合であっても機関吸気通路、燃焼室５およびＮＯX 吸収剤２３上流の排気通路内に供給された燃料（炭化水素）の量に対する空気の量の比をＮＯX 吸収剤２３への流入排気ガスの空燃比と称するとこのＮＯX 吸収剤２３は流入排気ガスの空燃比がリーンのときにはＮＯX を吸収し、流入排気ガスの空燃比が理論空燃比またはリッチになると吸収したＮＯX を放出するＮＯX の吸放出作用を行う。
【００３４】
このＮＯX 吸収剤２３を機関排気通路内に配置すればＮＯX 吸収剤２３は実際にＮＯX の吸放出作用を行うがこの吸放出作用の詳細なメカニズムについては明らかでない部分もある。しかしながらこの吸放出作用は図６に示したようなメカニズムで行われているものと考えられる。次にこのメカニズムについて担体上に白金ＰｔおよびバリウムＢａを担持させた場合を例にとって説明するが他の貴金属、アルカリ金属、アルカリ土類、希土類を用いても同様なメカニズムとなる。
【００３５】
図１に示した内燃機関では使用頻度の高い大部分の運転状態において空燃比がリーンの状態で燃焼が行われる。このように空燃比がリーンの状態で燃焼が行われている場合には排気ガス中の酸素濃度は高く、このときには図６（Ａ）に示したようにこれら酸素Ｏ2 がＯ2-またはＯ2-の形で白金Ｐｔの表面に付着する。一方、流入排気ガス中のＮＯは白金Ｐｔの表面上でＯ2-またはＯ2-と反応し、ＮＯ2 となる（２ＮＯ＋Ｏ2 →２ＮＯ2 ）。次いで生成されたＮＯ2 の一部は白金Ｐｔ上で酸化されつつ吸収剤内に吸収されて酸化バリウムＢａＯと結合しながら図６（Ａ）に示したように硝酸イオンＮＯ3-の形で吸収剤内に拡散する。このようにしてＮＯX がＮＯX 吸収剤２３内に吸収される。流入排気ガス中の酸素濃度が高い限り白金Ｐｔの表面でＮＯ2 が生成され、吸収剤のＮＯX 吸収能力が飽和しない限りＮＯ2 が吸収剤内に吸収されて硝酸イオンＮＯ3-が生成される。
【００３６】
一方、流入排気ガスの空燃比がリッチにされると流入排気ガス中の酸素濃度が低下し、その結果、白金Ｐｔの表面でのＮＯ2 の生成量が低下する。ＮＯ2 の生成量が低下すると反応が逆方向（ＮＯ3-→ＮＯ2 ）に進み、斯くして吸収剤内の硝酸イオンＮＯ3-がＮＯ2 の形で吸収剤から放出される。このときＮＯX 吸収剤２３から放出されたＮＯX は図６（Ｂ）に示したように流入排気ガス中に含まれる多量の未燃ＨＣ，ＣＯと反応して還元せしめられる。このようにして白金Ｐｔの表面上にＮＯ2 が存在しなくなると吸収剤から次から次へとＮＯ2 が放出される。したがって流入排気ガスの空燃比がリッチにされると短時間のうちにＮＯX 吸収剤２３からＮＯX が放出され、しかもこの放出されたＮＯX が還元されるために大気中にＮＯX が排出されることはない。
【００３７】
なおこの場合、流入排気ガスの空燃比を理論空燃比にしてもＮＯX 吸収剤２３からＮＯX が放出される。しかしながら流入排気ガスの空燃比を理論空燃比にした場合にはＮＯX 吸収剤２３からＮＯX が徐々にしか放出されないためにＮＯX 吸収剤２３に吸収されている全ＮＯX を放出させるには若干長い時間を要する。ところでＮＯX 吸収剤２３のＮＯX 吸収能力には限界があり、したがってＮＯX 吸収剤２３のＮＯX 吸収能力が飽和する前にＮＯX 吸収剤２３からＮＯX を放出させる必要がある。ところがＮＯX 吸収剤２３はＮＯX 吸収能力が十分なうちは排気ガス中に含まれるほとんど全てのＮＯX を吸収するがＮＯX 吸収能力の限界に近づくと一部のＮＯX を吸収しえなくなり、斯くしてＮＯX 吸収剤２３がＮＯX 吸収能力の限界に近づくとＮＯX 吸収剤２３から下流へ流出するＮＯX 量が増大しはじめる。
【００３８】
そこで本発明の第１実施例ではＮＯX 吸収剤２３に吸収されているトータルのＮＯX 吸収量を推定し、このＮＯX 吸収量が最大ＮＯX 吸収量に近づいたときにＮＯX 吸収剤２３に流入する排気ガスの空燃比を一時的にリッチにしてＮＯX 吸収剤２３からＮＯX を放出させるようにしている。この場合、ＮＯX 吸収剤２３に流入する排気ガスの空燃比をリッチにする方法は種々の方法がある。例えば燃焼室５内における混合気の平均空燃比をリッチにすることにより排気ガスの空燃比をリッチにすることもできるし、膨張行程末期または排気行程中に追加の燃料を噴射することによって排気ガスの空燃比をリッチにすることもできるし、またはＮＯX 吸収剤２３上流の排気通路内に追加の燃料を噴射することによって排気ガスの空燃比をリッチにすることもできる。本発明の実施例では１番目の方法、すなわちリッチ空燃比のもとで均一混合気燃焼を行わせることによって排気ガスの空燃比をリッチにするようにしている。
【００３９】
ところで排気ガス中にはＳＯX が含まれており、ＮＯX 吸収剤２３にはＮＯX ばかりでなくＳＯX も吸収される。このＮＯX 吸収剤２３へのＳＯX の吸収メカニズムはＮＯX の吸収メカニズムと同じであると考えられる。すなわちＮＯX の吸収メカニズムを説明したときと同様に担体上に白金ＰｔおよびバリウムＢａを担持させた場合を例にとって説明すると前述したように流入排気ガスの空燃比がリーンのときには酸素Ｏ2 がＯ2-またはＯ2-の形で白金Ｐｔの表面に付着しており、流入排気ガス中のＳＯ2 は白金Ｐｔの表面上でＯ2-またはＯ2-と反応してＳＯ3 となる。次いで生成されたＳＯ3 の一部は白金Ｐｔ上でさらに酸化されつつ吸収剤内に吸収されて酸化バリウムＢａＯと結合しながら、硫酸イオンＳＯ42- の形で吸収剤内に拡散し、安定した硫酸塩ＢａＳＯ4 を生成する。
【００４０】
しかしながらこの硫酸塩ＢａＳＯ4 は安定していて分解しづらく、流入排気ガスの空燃比を単にリッチにしただけでは硫酸塩ＢａＳＯ4 は分解されずにそのまま残る。したがってＮＯX 吸収剤２３内には時間が経過するにつれて硫酸塩ＢａＳＯ4 が増大することになり、斯くして時間が経過するにつれてＮＯX 吸収剤２３が吸収しうるＮＯX 量が低下することになる。すなわち時間が経過するにつれてＮＯX 吸収剤２３が劣化することになる。
【００４１】
ところがこの場合、ＮＯX 吸収剤２３の温度が一定温度、例えば６００℃以上になるとＮＯX 吸収剤２３内において硫酸塩ＢａＳＯ4 が分解し、このときＮＯX 吸収剤２３に流入する排気ガスの空燃比をリッチにするとＮＯX 吸収剤２３からＳＯX を放出させることができる。そこで本発明の実施例ではＮＯX 吸収剤２３からＳＯX を放出すべきときにＮＯX 吸収剤２３の温度が高い場合にはＮＯX 吸収剤２３に流入する排気ガスの空燃比をリッチにしてＮＯX 吸収剤２３からＳＯX を放出させ、ＳＯX を放出すべきときにＮＯX 吸収剤２３の温度が低い場合にはＮＯX 吸収剤２３の温度を上昇させると共にＮＯX 吸収剤２３に流入する排気ガスの空燃比をリッチにするようにしている。
【００４２】
次にＮＯX 吸収剤２３からＮＯX を放出すべくＮＯX 吸収剤２３に流入する排気ガスの空燃比をリッチにしたときの還元剤の量とＮＯX 吸収剤２３から下流へ流出する排気ガス中のアンモニアＮＨ3 の濃度との関係について説明する。まず初めに還元剤の量について説明する。ＮＯX 吸収剤２３に流入する排気ガスの空燃比を理論空燃比にするのに必要な燃料量に対して過剰な燃料はＮＯX の放出および還元のために使用されるのでこの過剰な燃料の量がＮＯX の放出および還元に使用される還元剤の量に一致する。このことはＮＯX 吸収剤２３からＮＯX を放出すべきときに燃焼室５内における混合気の空燃比をリッチにした場合でも、膨張行程末期または排気行程中に追加の燃料を噴射した場合でも、ＮＯX 吸収剤２３上流の排気通路内に追加の燃料を噴射した場合でも当てはまる。
【００４３】
次にアンモニアの濃度について説明する。空燃比がリーンのとき、すなわち酸化雰囲気のときにはアンモニアＮＨ3 はほとんど発生しない。ところが空燃比がリッチになると、すなわち還元雰囲気になると吸入空気中または排気ガス中の窒素Ｎ2 が酸化触媒または三元触媒２０において炭化水素ＨＣにより還元され、アンモニアＮＨ3 が生成される。しかしながら空燃比がリッチになるとＮＯX 吸収剤２３からＮＯX が放出され、生成されたアンモニアＮＨ3 はこのＮＯX を還元するために使用されるのでＮＯX 吸収剤２３からＮＯX が放出されている間、正確には供給された還元剤がＮＯX の放出および還元のために使用されている間はＮＯX 吸収剤２３から下流へアンモニアＮＨ3 は流出しない。これに対してＮＯX 吸収剤２３からのＮＯX の放出が完了した後も空燃比がリッチにされていると、より正確に言うとＮＯX 吸収剤２３からＮＯX を放出し還元するために使用されない余剰の還元剤が供給されるとアンモニアＮＨ3 はもはやＮＯX の還元のために消費されることがなくなり、斯くしてこのときにはＮＯX 吸収剤２３から下流へアンモニアＮＨ3 が流出することになる。
【００４４】
このことはＮＯX 吸収剤２３の上流に酸化触媒または三元触媒２０が設けられていない場合でも生ずる。すなわちＮＯX 吸収剤２３も還元機能を有する白金Ｐｔ等の触媒を具えているので空燃比がリッチになるとＮＯX 吸収剤２３においてアンモニアＮＨ3 が生成される可能性がある。しかしながらたとえアンモニアＮＨ3 が生成されたとしてもこのアンモニアＮＨ3 はＮＯX 吸収剤２３から放出されたＮＯX を還元するために使用されるためにＮＯX 吸収剤２３から下流へはアンモニアＮＨ3 が流出しない。ところがＮＯX 吸収剤２３からＮＯX を放出し還元するために使用されない余剰の還元剤が供給されると前述したようにＮＯX 吸収剤２３から下流へアンモニアＮＨ3 が流出することになる。
【００４５】
このようにＮＯX 吸収剤２３に流入する排気ガスの空燃比がリッチにされたときにＮＯX 吸収剤２３からＮＯX を放出し還元するために使用されない余剰の還元剤が供給されるとこの余剰の還元剤はアンモニアＮＨ3 の形でＮＯX 吸収剤２３から下流へ流出し、このとき流出するアンモニア量は余剰の還元剤の量に比例する。したがってこのとき流出するアンモニア量から余剰の還元剤量がわかることになる。このアンモニア量はアンモニア濃度を検出可能なＮＯX アンモニアセンサ２９によって検出される。この場合、このアンモニア濃度の積算値は余剰の還元剤量を表していると考えられ、したがってアンモニア濃度の積算値は余剰の還元剤量を表わす代表値であると言える。またこのアンモニア濃度の最大値が余剰の還元剤量を表していると考えることもでき、したがってアンモニア濃度の最大値は余剰の還元剤量を表わす代表値であると言える。
【００４６】
次に図７を参照しつつ還元剤の供給制御の第１実施例について説明する。図７を参照するとΣＮＯＸはＮＯX 吸収剤２３に吸収されているトータルのＮＯX 量（総ＮＯX 吸収量）を示しており、Ｉ1 はＮＯX アンモニアセンサ２９の検出電流を示している。なお図７においてＮＯX およびＮＨ3 は排気ガス中のＮＯX 濃度およびＮＨ3 濃度の変化によるＮＯX アンモニアセンサ２９の検出電流の変化を夫々示しており、これら検出電流の双方がＮＯX アンモニアセンサ２９の検出電流Ｉ1 に表われる。またＥは空燃比センサ３０の出力電圧を示しており、Ａ／Ｆは燃焼室５内における混合気の平均空燃比を示している。
【００４７】
図７に示したように総ＮＯX 吸収量ΣＮＯＸが増大してＮＯX 吸収剤２３の吸収能力限界に近づくとＮＯX 吸収剤２３から下流へとＮＯX が流出しはじめるのでＮＯX アンモニアセンサ２９の検出電流Ｉ1 が上昇しはじめる。図７に示した実施例ではＮＯX 吸収剤２３の総ＮＯX 吸収量を推定し、ＮＯX 吸収剤２３に流入する排気ガスの空燃比をリッチにした後に再びリッチにするまでのリッチ時間間隔をこの総ＮＯX 吸収量の推定値に基づき制御する。すなわちこの第１実施例ではＮＯX 吸収剤２３に吸収されているトータルのＮＯX 量を推定するための総ＮＯX 吸収量推定手段を具備しており、図７に示したように総ＮＯX 吸収量推定手段により推定された総ＮＯX 吸収量ΣＮＯＸが許容値ＮＯＸmax −αを越えたときに空燃比をリーンからリッチに一時的に切換えるようにしている。ここでＮＯＸmax はＮＯX 吸収剤２３に吸収しうる最大ＮＯX 吸収量である。
【００４８】
空燃比Ａ／Ｆがリーンからリッチに切換えられてもリッチ空燃比の排気ガスがＮＯX 吸収剤２３に到達するには時間を要するので空燃比Ａ／Ｆがリッチに切換えられた直後はＮＯX 吸収剤２３から下流へ流出するＮＯX 量が増大し続ける。次いでリッチ空燃比の排気ガス中に含まれる還元剤によるＮＯX の還元作用が開始されるためにＮＯX 吸収剤２３から下流へはＮＯX が流出しなくなる。したがって空燃比がリーンからリッチに切換えられるとＮＯX アンモニアセンサ２９の検出電流Ｉ1 は短時間上昇した後、零まで低下する。一方、空燃比がリーンからリッチに切換えられるとＮＯX 吸収剤２３からのＮＯX の放出作用が開始され、斯くしてＮＯX 吸収剤２３に吸収されているＮＯX 量ΣＮＯＸは徐々に減少する。
【００４９】
次に第１実施例における総ＮＯX 吸収量の算出方法について説明する。ＮＯX 吸収剤２３に流入する排気ガス（以下、流入排気ガスと称す。）の空燃比がリーンであるときに単位時間当たりにＮＯX 吸収剤２３に吸収されるＮＯX の量（以下、単位ＮＯX 吸収量と称す。）はＮＯX 吸収剤２３に最大限に吸収しうるＮＯX の量（以下、最大ＮＯX 吸収量と称す。）と、ＮＯX 吸収剤２３へＮＯX が吸収される速度に関する係数（以下、ＮＯX 吸収速度）と、流入排気ガス中のＮＯX 濃度（以下、流入ＮＯX 濃度と称す。）と、ＮＯX 吸収剤２３に現在吸収されているトータルのＮＯX の量（総ＮＯX 吸収量）とに基づいて算出可能である。
【００５０】
そしてこれらパラメータ、すなわち単位ＮＯX 吸収量と、最大ＮＯX 吸収量と、ＮＯX 吸収速度と、流入ＮＯX 濃度と、総ＮＯX 吸収量との間には或る一定の関係が成り立つ。そこで本実施例ではこれパラメータ間に成り立つ関係式を実験により予め求めておき、所定の時間間隔ごとに当該関係式に最大ＮＯX 吸収量と、ＮＯX 吸収速度と、流入ＮＯX 濃度と、総ＮＯX 吸収量とを代入して単位ＮＯX 吸収量を算出し、斯くして算出された単位ＮＯX 吸収量を積算することにより総ＮＯX 吸収量を算出する。なお本実施例において実験により求められる関係式は図８（Ａ）に示した式である。図示した関係式においてＡａｂｃは単位ＮＯX 吸収量、ＮＯＸmax は最大ＮＯX 吸収量、ＫａｂはＮＯX 吸収速度、Ｃnox は流入ＮＯX 濃度、ΣＮＯＸは総ＮＯX 吸収量である。
【００５１】
ここで最大ＮＯX 吸収量とＮＯX 吸収速度との初期値としては実験により予め求めた値を用いる。これら最大ＮＯX 吸収量およびＮＯX 吸収速度はＮＯX 吸収剤２３が劣化しない限りにおいては機関運転状態の変化とは関係なくほぼ一定である。しかしながらＮＯX 吸収剤２３の使用が開始されるとＮＯX 吸収剤２３はＳＯX の吸収や排気ガスの熱などにより徐々に劣化し、最大ＮＯX 吸収量およびＮＯX 吸収速度は変化する。したがって本実施例では後述する方法によりこれら最大ＮＯX 吸収量およびＮＯX 吸収速度を適宜修正する。一方、流入ＮＯX 濃度は機関運転状態の変化に伴って変化し、総ＮＯX 吸収量も時間の経過と共に変化するのでこれら流入ＮＯX 濃度および総ＮＯX 吸収量としては単位ＮＯX 吸収量算出時にその都度算出した値を用いる。
【００５２】
また流入ＮＯX 濃度Ｃnox は次のようにして算出される。すなわちＮＯX 吸収剤２３に単位時間当たりに流入するＮＯX の量（以下、単位ＮＯX 流入量と称す。）は機関回転数と機関負荷との関数となるので単位ＮＯX 流入量ＮＡを図９に示したようにマップの形で機関回転数Ｎと機関負荷Ｑ／Ｎとの関数で予めＲＯＭ３４に記憶しておき、当該マップに基づいて算出される単位ＮＯX 流入量ＮＡを単位時間当たりの吸気量にて除算することにより流入ＮＯX 濃度を算出することができる。
【００５３】
ところで上述したように最大ＮＯX 吸収量はＮＯX 吸収剤２３の劣化に伴って徐々に少なくなり、ＮＯX 吸収速度は徐々に遅くなる。この場合には最大ＮＯX 吸収量およびＮＯX 吸収速度をそのまま用いて単位ＮＯX 吸収量を算出しても正確な単位ＮＯX 吸収量は算出されない。そこで本実施例では以下の方法により最大ＮＯX 吸収量とＮＯX 吸収速度とを適宜修正し、これにより正確な最大ＮＯX 吸収量とＮＯX 吸収速度とを獲得する。
【００５４】
すなわちＮＯX 吸収剤２３下流にはＮＯX 吸収剤２３から流出する排気ガス中のＮＯX 濃度（以下、流出ＮＯX 濃度と称す。）を検出しうるＮＯX アンモニアセンサ２９が配置されているのでこのＮＯX アンモニアセンサ２９の出力を用いれば単位時間当たりにＮＯX 吸収剤２３から流出するＮＯX の量（以下、単位ＮＯX 流出量と称す。）Ａｏｕｔｓを算出することができる。ここで上述したように単位ＮＯX 流入量ＮＡは図９のマップから算出可能であるので図８（Ｂ）に示した関係式に従って斯くして算出した単位ＮＯX 流入量ＮＡから単位ＮＯX 流出量Ａｏｕｔｓを差し引けば単位ＮＯX 吸収量Ａａｂｓが算出される。なお単位ＮＯX 流出量ＡｏｕｔはＮＯX アンモニアセンサ２９の出力電流Ｉ1 に所定の係数Ｋ1 を掛けて算出される。
【００５５】
斯くして算出した単位ＮＯX 吸収量（以下、実測値と称す。）Ａａｂｓと上記関係式に基づいて算出した単位ＮＯX 吸収量（以下、理論値と称す。）Ａａｂｃとを比較したときにＮＯX 吸収剤２３が劣化していなければこれら実測値Ａａｂｓと理論値Ａａｂｃはほぼ等しくなる。しかしながらＮＯX 吸収剤２３が劣化しているときにはこれら実測値Ａａｂｓと理論値Ａａｂｃとは相当にずれている。したがってこの場合には今回、上記関係式を用いて理論値Ａａｂｃを算出する以前に理論値Ａａｂｃを算出したときに用いた単位ＮＯX 吸収量、流入ＮＯX 濃度、および総ＮＯX 吸収量を記憶しておき、これらパラメータと、今回、上記関係式を用いて理論値Ａａｂｃを算出したときに用いた流入ＮＯX 濃度、総ＮＯX 吸収量、および単位ＮＯX 吸収量の実測値をそれぞれ上記関係式に代入して最大ＮＯX 吸収量ＮＯＸmax とＮＯX 吸収速度Ｋａｂとを改めて算出する。すなわち本実施例では最大ＮＯX 吸収量ＮＯＸmax とＮＯX 吸収速度ＫａｂとはＮＯX アンモニアセンサ２９の出力により修正されることになる。これによればＮＯX 吸収剤２３の最大ＮＯX 吸収量とＮＯX 吸収速度とを正確に把握することができる。
【００５６】
ところで空燃比がリーンからリッチに切換えられると過剰な燃料、すなわち還元剤はＮＯX を還元するために消費されるのでＮＯX 吸収剤２３から下流へ流出する排気ガスの空燃比はほぼ理論空燃比となる。この場合、なぜそうなるのかという理由は明確ではないがＮＯX 吸収剤２３が劣化していないときにはＮＯX 吸収剤２３から下流へ流出する排気ガスの空燃比はわずかばかりリーンとなり、ＮＯX 吸収剤２３が劣化するとＮＯX 吸収剤２３から下流へ流出する排気ガスの空燃比はわずかばかりリッチとなる傾向がある。しかしながらいずれの場合でもＮＯX 吸収剤２３からのＮＯX の放出作用が完了する頃になるとＮＯX 吸収剤２３から下流へ流出する排気ガスの空燃比はリッチとなる。
【００５７】
図７は空燃比がリーンからリッチに切換えられたときにＮＯX 吸収剤２３から下流へ流出する排気ガスの空燃比がわずかばかりリーンとなっている場合を示しており、ＮＯX 吸収剤２３からのＮＯX の放出作用が完了する頃になると、すなわち総ＮＯX 吸収量ΣＮＯＸが零に近づくと空燃比センサ３０の出力電圧Ｅがリッチであることを示す出力信号レベルに向けて変化し、すなわち上昇することがわかる。この出力信号レベルＥの変化は応答性がよく、したがってこの出力信号レベルＥの変化に基づいて空燃比をリッチからリーンに切換えればＮＯX 吸収剤２３からのＮＯX の放出作用が完了したときに空燃比をリッチからリーンに切換えることができる。
【００５８】
そこで図７に示した実施例では空燃比センサ３０の出力電圧Ｅに対して基準電圧ＥS を設定しておき、すなわち一般的な表現を用いると空燃比検出手段の出力信号レベルＥに対して基準レベルＥS を設定しておき、出力信号レベルＥが基準レベルＥS を越えたときに空燃比をリッチからリーンに切換えるようにしている。
【００５９】
ところで空燃比センサ３０の出力電圧ＥはＮＯX の放出作用の完了に対して応答性よく変化するが空燃比センサ３０やＮＯX 吸収剤２３の性能のばらつき、或いは経時変化によって出力電圧Ｅは変化のしかたがいろいろと変わる。したがって基準レベルＥS を一定値に固定しておくとＮＯX の放出完了時に空燃比をリッチからリーンに切換えられなく場合が生じてくる。
【００６０】
一方、空燃比がリーンからリッチに切換えられたときにＮＯX 吸収剤２３からＮＯX を放出し還元するために使用されない余剰の還元剤が供給されたとするとこのときにはＮＯX 吸収剤２３から下流へアンモニアＮＨ3 が流出するので図７に示したようにＮＯX アンモニアセンサ２９の検出電流Ｉ1 が上昇する。この場合、図７においてハッチングで示される検出電流Ｉ1 の積算値ΣＩおよび検出電流Ｉ1 の最大値Ｉmax は余剰の還元剤量を表している。
【００６１】
このＮＯX アンモニアセンサ２９の検出電流Ｉ1 はＮＯX の放出完了に対して応答遅れがあるが余剰の還元剤量は検出電流Ｉ1 から正確に求まる。そこで本発明ではＮＯX アンモニアセンサ２９の検出電流Ｉ1 の変化に基づいて、すなわちアンモニア濃度の変化に基づいてＮＯX 吸収剤２３からのＮＯX の放出が完了したときに排気ガスの空燃比がリッチからリーンに切換えられるように基準レベルＥS を変化させるようにしている。
【００６２】
具体的に言うと検出電流Ｉ1 の積算値ΣＩ、或いは検出電流Ｉ1 の最大値Ｉmax に対して小さな値の目標値を予め設定しておき、ΣＩまたはＩmax が目標値よりも大きくなったとき、すなわち余剰の還元剤量が多いときにはリッチからリーンへの空燃比の切換時期を早めて余剰の還元剤量が少くなくなるように基準レベルＥS を低下させ、すなわち基準レベルＥS をリーンであることを示す出力信号レベル側へ変化させ、ΣＩまたはＩmax が目標値よりも小さくなったとき、すなわち余剰の還元剤量が零または零に近いときにはリッチからリーンへの空燃比の切換時期を遅らせて余剰の還元剤量が多くなくなるように基準レベルＥS を上昇させ、すなわち基準レベルＥS をリッチであることを示す出力信号レベル側へ変化させるようにしている。
【００６３】
ところでＮＯX 吸収剤２３に流入する排気ガスの空燃比がリーンからリッチに切換えられたときにはＮＯX 吸収剤２３に吸収されているＮＯX が排気ガス中の還元剤（すなわち炭化水素）によって放出・還元される。したがって総ＮＯX 吸収量ΣＮＯＸを徐々に減少すべきである。そこで次に本実施例においてＮＯX 吸収剤２３に流入する排気ガスの空燃比がリッチであるときの総ＮＯX 吸収量ΣＮＯＸの算出方法について説明する。
【００６４】
ＮＯX 吸収剤２３に流入する排気ガスの空燃比がリーンであるときに単位時間当たりにＮＯX 吸収剤２３から放出・還元されるＮＯX の量（以下、単位ＮＯX 還元量と称す。）はＮＯX 吸収剤２３からＮＯX が放出・還元される速度に関する係数（以下、ＮＯX 還元速度と称す。）と、流入排気ガス中の還元剤濃度（以下、流入還元剤濃度と称す。）と、ＮＯX 吸収剤２３に現在吸収されているトータルのＮＯX の量（総ＮＯX 吸収量）とに基づいて算出可能である。
【００６５】
そしてこれらパラメータ、すなわち単位ＮＯX 還元量と、ＮＯX 還元速度と、流入還元剤濃度と、総ＮＯX 吸収量との間には或る一定の関係が成り立つ。そこで本実施例ではこれらパラメータ間に成り立つ関係式を実験式により予め求めておき、所定の時間間隔ごとに当該関係式にＮＯX 還元速度と、流入還元剤濃度と、総ＮＯX 吸収量とを代入して単位ＮＯX 還元量を算出し、斯くして算出された単位ＮＯX 還元量をその時の総ＮＯX 吸収量から減算することにより総ＮＯX 吸収量を算出する。なお本実施例において実験により求められる関係式は図８（Ｄ）に示した式である。図示した関係式においてＡｒｅは単位ＮＯX 還元量、ＫｒｅはＮＯX 還元速度、Ｃｈｃは流入還元剤濃度、ΣＮＯＸは総ＮＯX 吸収量である。
【００６６】
ここでＮＯX 還元速度の初期値としては実験により予め求めた値を用いる。このＮＯX 還元速度はＮＯX 吸収剤２３が劣化しない限りにおいては機関運転状態の変化とは関係なくほぼ一定である。しかしながらＮＯX 吸収剤２３の使用が開始されるとＮＯX 吸収剤２３はＳＯX の吸収や排気ガスの熱などにより徐々に劣化し、ＮＯX 還元速度は変化する。したがって本実施例では後述する方法によりこのＮＯX 還元速度を適宜修正する。一方、流入還元剤濃度は機関運転状態、具体的には燃焼室５内に導入せしめられる空気の量と燃料噴射弁１０から噴射される燃料の量との変化に伴って変化し、総ＮＯX 吸収量も時間の経過と共に変化するのでこれら流入還元剤濃度および総ＮＯX 吸収量としては単位ＮＯX 還元量算出時にその都度算出した値を用いる。
【００６７】
ところで上述したようにＮＯX 還元速度ＫｒｅはＮＯX 吸収剤２３の劣化に伴って徐々に遅くなる。この場合にはＮＯX 還元速度Ｋｒｅをそのまま用いて単位ＮＯX 還元量Ａｒｅを算出しても正確な単位ＮＯX 還元量は算出されない。そこで本実施例では以下の方法によりＮＯX 還元速度Ｋｒｅを適宜修正し、これにより正確なＮＯX 還元速度を獲得する。
【００６８】
すなわちＮＯX 吸収剤２３下流にはＮＯX 吸収剤２３から流出する排気ガス中のアンモニア濃度（以下、流出アンモニア濃度と称す。）を検出しうるＮＯX アンモニアセンサ２９が配置されているのでこのＮＯX アンモニアセンサ２９の出力を用いればＮＯX 吸収剤２３におけるＮＯX の放出・還元が完了したことを知ることができる。ここで流入排気ガスの空燃比がリッチである間に総ＮＯX 吸収量を算出するために用いられるＮＯX 還元速度が真の値であればＮＯX アンモニアセンサ２９がＮＯX 吸収剤２３におけるＮＯX の放出・還元が完了したことを検出する前においては総ＮＯX 吸収量は零とはならず、少なくとも或る一定の値よりも多く、一方、ＮＯX アンモニアセンサ２９がＮＯX 吸収剤２３におけるＮＯX の放出・還元が完了したことを検出したときにおいては総ＮＯX 吸収量は零となるか、或いは少なくとも或る一定の値よりも少なくなるはずである。
【００６９】
言い換えればＮＯX アンモニアセンサ２９がＮＯX 吸収剤におけるＮＯX 放出・還元が完了したことを検出する前において総ＮＯX 吸収量が零となるか、或いは或る一定の値よりも少なくなったときには総ＮＯX 吸収量を算出するために用いたＮＯX 還元速度が速すぎたと判断することができるので本実施例ではこの場合には現在使用しているＮＯX 還元速度を所定値だけ遅くする。ここでの所定値は一定値でもよいし、或いは総ＮＯX 吸収量が零となったとき、或いは或る一定の値よりも少なくなったときからＮＯX アンモニアセンサ２９によりＮＯX 吸収剤２３におけるＮＯX の放出・還元が完了したことが検出されるまでの時間に基づいて当該時間が長いほど大きくなる値でもよい。
【００７０】
一方、ＮＯX アンモニアセンサ２９がＮＯX 吸収剤２３におけるＮＯX の放出・還元が完了したことを検出したときに総ＮＯX 吸収量が零となっていないか、或いは依然として或る一定の値よりも多いときには総ＮＯX 吸収量を算出するために用いたＮＯX 還元速度が遅すぎたと判断することができるので本実施例ではこの場合には現在使用しているＮＯX 還元速度を所定値だけ速くする。ここでの所定値は一定値でもよいし、或いはＮＯX アンモニアセンサ２９によりＮＯX 吸収剤２３におけるＮＯX の放出・還元が完了したことが検出されたときにおける総ＮＯX 吸収量に基づいて当該総ＮＯX 吸収量が多いほど大きくなる値でもよい。このように本実施例ではＮＯX 還元速度はＮＯX アンモニアセンサ２９の出力により修正されることになる。これによればＮＯX 還元速度を正確に把握することができる。
【００７１】
ところで上述したように本実施例では上記関係式にて使用される最大ＮＯX 吸収量、ＮＯX 吸収速度、およびＮＯX 還元速度はＮＯX 吸収剤２３の劣化度合に応じて修正せしめられる。言い換えればこれら最大ＮＯX 吸収量、ＮＯX 吸収速度、およびＮＯX 還元速度からＮＯX 吸収剤２３の劣化度合を知ることができる。すなわち修正せしめられた最大ＮＯX 吸収量がＮＯX 吸収剤２３の劣化が許容範囲内にあるならば達成されたであろう量よりも少ないとき、或いは修正せしめられたＮＯX 吸収速度またはＮＯX 還元速度が同様にＮＯX 吸収剤２３の劣化度合が許容範囲内にあるならば達成されたであろう速度よりも遅いときには本実施例ではＮＯX 吸収剤２３が劣化したと診断する。これによればＮＯX 吸収剤２３の劣化を正確に診断することができる。
【００７２】
図１０は第１実施例を実行するためのルーチンを示している。図１０を参照するとまず初めにステップ１００において図５（Ｂ）に示したマップから基本燃料噴射量ＴＡＵが算出される。次いでステップ１０１ではＮＯX 吸収剤２３からＮＯX を放出すべきことを示すＮＯX 放出フラグがセットされているか否かが判別される。ＮＯX 放出フラグがセットされていないときにはステップ１０２に進んで後述する図１１および図１２のルーチンにより算出される総ＮＯX 吸収量ΣＮＯＸが最大ＮＯX 吸収量ＮＯＸmax よりも値αだけ少ない値を越えたか否かが判別される。ΣＮＯＸ≦ＮＯＸmax −αのとき、すなわちＮＯX 吸収剤２３のＮＯX 吸収能力に未だ余裕があるときにはステップ１０４にジャンプする。ステップ１０４では図５（Ｃ）に示したマップから補正係数Ｋが算出される。次いでステップ１０５では基本燃料噴射量ＴＡＵに補正係数Ｋを乗算することによって最終的な燃料噴射量ＴＡＵＯ（＝Ｋ・ＴＡＵ）が算出され、この噴射量ＴＡＵＯでもって燃料噴射が行われる。次いでステップ１０６ではＮＯX 吸収剤２３からのＳＯX を放出するためのＳＯX 放出処理を行うべきか否かが判断される。ＳＯX 放出処理を行う必要のないときには処理サイクルを完了する。
【００７３】
一方、ステップ１０２においてΣＮＯＸ＞ＮＯＸmax −αになったと判断されたときにはステップ１０３に進んでＮＯX 放出フラグがセットされ、次いでステップ１０３ａに進んでＮＨ3 検出フラグがセットされる。次いでステップ１０４に進む。
ＮＯX 放出フラグがセットされると次の処理サイクルではステップ１０１からステップ１０８に進んでリッチ補正係数ＫR が算出される。次いでステップ１０９では基本燃料噴射量ＴＡＵにリッチ補正係数ＫR とリッチ修正係数ＫS とを乗算することによって最終的な燃料噴射量ＴＡＵＯ（＝ＫR ・ＴＡＵ・ＫS ）が算出され、この噴射量ＴＡＵＯでもって燃料噴射が行われる。なおリッチ修正係数ＫS は図１１および図１２のルーチンにおいてＮＯX 還元速度を修正するのに伴って算出される係数である。ステップ１０９によればリーン空燃比のもとでの成層燃焼またはリーン空燃比のもとでの均一混合気燃焼からリッチ空燃比のもとでの均一混合気燃焼に切換えられ、それによってＮＯX 吸収剤２３からのＮＯX の放出作用が開始される。
【００７４】
次いでステップ１１０では空燃比センサ３０の出力電圧Ｅが基準電圧ＥS を越えたか否かが判別される。Ｅ≦ＥS のときにはステップ１０６に進む。これに対してＥ＞ＥS になるとステップ１１１に進んでＮＯX 放出フラグがリセットされる。ＮＯX 放出フラグがリセットされると空燃比がリッチからリーンに切換えられる。一方、ステップ１０６においてＳＯX 放出処理を行うべきであると判断されるとステップ１０７に進んでＮＯX 吸収剤２３からＳＯX を放出させる処理が行われる。すなわちＮＯX 吸収剤２３の温度をほぼ６００℃以上に維持しつつ空燃比がリッチとされる。
【００７５】
図１１および図１２は本実施例の総ＮＯX 吸収量ΣＮＯＸを算出するためのルーチンを示している。図１１において初めにステップ２００においてＮＯX 吸収剤２３に流入する排気ガスの空燃比がリーンであるか否かが判別される。ステップ２００において流入排気ガスの空燃比がリーンであると判別されるとステップ２０１において機関回転数Ｎと機関負荷Ｌと吸入空気量Ｑとが算出され、次いでステップ２０２において図９に示したマップから機関回転数Ｎと機関負荷Ｌとに基づいてＮＯX 吸収剤２３に流入するＮＯX の量ＮＡが算出される。次いでステップ２０３においてＮＯX 吸収剤２３に流入する排気ガス中のＮＯX 濃度Ｃnox が算出され、次いでステップ２０４において図８（Ａ）に示した関係式から単位時間当たりにＮＯX 吸収剤２３に吸収されるＮＯX の量（単位ＮＯX 吸収量）の理論値Ａａｂｃが算出される。次いでステップ２０５において図８（Ｂ）および図８（Ｃ）の関係式から単位ＮＯX 吸収量の実測値Ａａｂｓが算出される。
【００７６】
次いでステップ２０７において理論値Ａａｂｃと実測値Ａａｂｓとの間のずれが所定値βの範囲にあるか否かが判別される。ステップ２０７においてＡａｂｓ−β＜Ａａｂｃ＜Ａａｂｓ＋βであると判別されたときには図８（Ａ）に示した関係式にて用いられた最大ＮＯX 吸収量ＮＯＸmax およびＮＯX 吸収速度Ｋａｂの値が真の値であると判断し、ステップ２０８において単位ＮＯX 吸収量Ａａｂとしてその理論値Ａａｂｃを採用し、ステップ２０９においてこの単位ＮＯX 吸収量Ａａｂが現在の総ＮＯX 吸収量ΣＮＯＸに加算され、新たな総ＮＯX 吸収量ΣＮＯＸが算出され、次いでステップ２０９ａにおいて今回のルーチン実行時における最大ＮＯX 吸収量ＮＯＸmax と、流入排気ガス中のＮＯX 濃度Ｃnox と、単位ＮＯX 吸収量Ａａｂとが保存される。
【００７７】
一方、ステップ２０７においてＡａｂｓ−β≧ＡａｂｃまたはＡａｂｃ≧Ａａｂｓ＋βであると判別されたときには図８（Ａ）に示した関係式にて用いられた最大ＮＯX 吸収量ＮＯＸmax またはＮＯX 吸収速度Ｋａｂの値が真の値ではないと判断し、ステップ２１１において単位ＮＯX 吸収量Ａａｂとしてその実測値Ａａｂｓを採用し、次いでステップ２１２において今回のルーチン以前にステップ２０９ａにおいて保存された単位ＮＯX 吸収量Ａａｂと、流入排気ガス中のＮＯX 濃度Ｃnox と、総ＮＯX 吸収量ΣＮＯＸとを図８（Ａ）に示した関係式に代入して１つの方程式を獲得し、さらに今回のルーチン実行時における単位ＮＯX 吸収量Ａａｂと、流入排気ガス中のＮＯX 濃度Ｃnox と、総ＮＯX 吸収量ΣＮＯＸとを同様に図８（Ａ）に示した関係式に代入して別の１つの方程式を獲得し、これら２つの方程式から求められる２つの値ＮＯＸmax ，Ｋａｂをそれぞれ新たな最大ＮＯX 吸収量およびＮＯX 吸収速度とすることでこれらパラメータを修正し、ステップ２０９に進む。
【００７８】
ところでステップ２００において流入排気ガスの空燃比がリッチであると判別されると図１２のステップ２１３において機関回転数Ｎと機関負荷Ｌと吸入空気量Ｑとが算出され、次いでステップ２１４において機関回転数Ｎと機関負荷Ｌと吸入空気量Ｑとに基づいてＮＯX 吸収剤２３に流入する排気ガス中の還元剤、すなわち炭化水素の濃度Ｃｈｃが算出される。すなわち燃料噴射弁１０から噴射される燃料の量は機関回転数Ｎと機関負荷Ｌとに基づいて決定され、斯くして決定される燃料噴射量のうち燃焼室５内にて燃焼しない燃料量を吸入空気量Ｑで割れば流入排気ガス中の還元剤濃度が算出されるので結局のところ機関回転数Ｎと、機関負荷Ｌと、吸入空気量Ｑとに基づいて流入排気ガス中の還元剤濃度Ｃｈｃを算出することができる。
【００７９】
次いでステップ２１５において図８（Ｄ）に示した関係式から単位ＮＯX 還元量Ａｒｅが算出され、次いでステップ２１６において現在の総ＮＯX 吸収量ΣＮＯＸから単位ＮＯX 還元量Ａｒｅが差し引かれて新たに総ＮＯX 吸収量ΣＮＯＸが算出され、ステップ２１７に進む。
ステップ２１７ではＮＯX アンモニアセンサ２９の出力電流Ｉ1 が基準値Ｉｔを越えたか否かが判別される。ステップ２１７においてＩ1 ＞Ｉｔであるときにはステップ２１８に進んでＮＯX 放出フラグがリセットされる。ここで図１０のステップ１１１においてＮＯX 放出フラグがリセットされていなければステップ２１８においてＮＯX 放出フラグがリセットされることによりＮＯX 吸収剤２３に流入する排気ガスの空燃比がリッチからリーンとされ、したがってＮＯX 吸収剤２３からのＮＯX の放出・還元が図１０のフローチャートによるのではなく強制的に終了せしめられる。
【００８０】
次いでステップ２１９において総ＮＯX 吸収量ΣＮＯＸが判定値Ａよりも多いか否かが判別される。ステップ２１９においてΣＮＯＸ≦Ａであれば図８（Ｄ）の関係式にて用いられたＮＯX 還元速度Ｋｒｅの値が真の値であると判断し、現在使用しているＮＯX 還元速度を修正することなくそのままルーチンを終了する。一方、ステップ２１９においてΣＮＯＸ＞Ａであると判別されたときにはステップ２２０に進んで図８（Ｄ）の関係式にて用いられるＮＯX 還元速度Ｋｒｅが所定の値だけ速くせしめられ、次いでステップ２２１において流入排気ガスの空燃比をリッチとするときに図１０のステップ１０９において使用されるリッチ修正係数ＫS がリッチ度合が大きくなるように修正される。
【００８１】
一方、ステップ２１７においてＩ1 ≦Ｉｔであると判別されたときにはステップ２２２において総ＮＯX 吸収量ΣＮＯＸが判定値Ａよりも少ないか否かが判別される。ステップ２２２においてΣＮＯＸ≧Ａであると判別されたときには図８（Ｄ）の関係式にて用いられたＮＯX 還元速度Ｋｒｅの値が真の値であると判断し、現在使用しているＮＯX 還元速度を修正することなくそのままルーチンを終了する。一方、ステップ２２２においてΣＮＯＸ＜Ａであると判別されたときにはステップ２２３に進んで図８（Ｄ）の関係式にて用いられるＮＯX 還元速度Ｋｒｅが所定の値だけ遅くせしめられ、次いでステップ２２４において流入排気ガスの空燃比をリッチとするときに図１０のステップ１０９において使用されるリッチ修正係数ＫS がリッチ度合が小さくなるように修正される。
【００８２】
図１３は目標レベルＥS を算出するためのルーチンを示している。図１３を参照するとまず初めにステップ３００においてＮＨ3 検出フラグがセットされているか否かが判別される。このＮＨ3 検出フラグは図１０のステップ１０２においてΣＮＯＸ＞ＮＯＸmax −αとなったときにセットされる。ＮＨ3 検出フラグがセットされているときにはステップ３０１に進んでＮＨ3 検出フラグがセットされてからの経過時間ｔが一定時間ｔ1 を越えたか否かが判別される。この一定時間ｔ1 は空燃比がリーンからリッチにされた後にＮＯX アンモニアセンサ２９の検出電流Ｉ1 が零まで低下し終わるまでの時間である。ｔ＞ｔ1 になるとステップ３０２に進んでＮＨ3 検出フラグがセットされてからの経過時間ｔが一定時間ｔ2 を越えたか否かが判別される。この一定時間ｔ2 はＮＯX 吸収剤２３から下流へアンモニアが流出したときにどのようなアンモニア量であってもＮＯX アンモニアセンサ２９がアンモニア濃度を検出しうるのに十分な時間である。ｔ≦ｔ2 のときにはステップ３０３に進む。
【００８３】
ステップ３０３ではＮＯX アンモニアセンサ２９の検出電流Ｉ1 が算出される。次いでステップ３０４ではこの検出電流Ｉ1 をΣＩに加算することによって検出電流の積算値ΣＩが算出される。次いでステップ３０２においてｔ＞ｔ2 になったと判別されたときにはステップ３０５に進んで検出電流の積算値ΣＩが目標値Ｓr よりも大きいか否かが判別される。ΣＩ＞Ｓr のときにはステップ３０６に進んで目標レベルＥS が予め定められた設定値ａだけ減少せしめられ、次いでステップ３０８に進む。これに対してΣＩ≦Ｓr のときにはステップ３０７に進んで目標レベルＥS が予め定められた設定値ａだけ増大せしめられ、次いでステップ３０８に進む。ステップ３０８ではΣＩがクリアされ、ＮＨ3 検出フラグがリセットされる。
【００８４】
図１４は目標レベルＥS を算出するためのルーチンの別の例を示している。図１４を参照するとまず初めにステップ４００においてＮＨ3 検出フラグがセットされているか否かが判別される。このＮＨ3 検出フラグは図１０のステップ１０２においてΣＮＯＸ＞ＮＯＸmax −αとなったときにセットされる。ＮＨ3 検出フラグがセットされているときにはステップ４０１に進んでＮＨ3 検出フラグがセットされてからの経過時間ｔが一定時間ｔ1 を越えたか否かが判別される。この一定時間ｔ1 は前述したように空燃比がリーンからリッチにされた後にＮＯX アンモニアセンサ２９の検出電流Ｉ1 が零まで低下し終わるまでの時間である。ｔ＞ｔ1 になるとステップ４０２に進んでＮＨ3 検出フラグがセットされてからの経過時間ｔが一定時間ｔ2 を越えたか否かが判別される。この一定時間ｔ2 は前述したようにＮＯX 吸収剤２３から下流へアンモニアが流出したときにどのようなアンモニア量であってもＮＯX アンモニアセンサ２９がアンモニア濃度を検出しうるのに十分な時間である。ｔ≦ｔ2 のときにはステップ４０３に進む。
【００８５】
ステップ４０３では検出電流Ｉ1 がＩmax よりも大きいか否かが判別される。Ｉ1 ＞Ｉmax のときにはステップ４０４に進んでＩ1 が検出電流の最大値Ｉmax とされる。次いでステップ４０２においてｔ＞ｔ2 になったと判別されたときにはステップ４０５に進んで検出電流の最大値Ｉmax が目標値Ｉmaxrよりも大きいか否かが判別される。Ｉmax ＞Ｉmaxrのときにはステップ４０６に進んで目標レベルＥS が予め定められた設定値ａだけ減少せしめられ、次いでステップ４０８に進む。これに対してＩmax ≦Ｉmaxrのときにはステップ４０７に進んで目標レベルＥS が予め定められた設定値ａだけ増大せしめられ、次いでステップ４０８に進む。ステップ４０８ではＩmax がクリアされ、ＮＨ3 検出フラグがリセットされる。
【００８６】
図１５はＮＯX 吸収剤２３の劣化を判定するためのルーチンを示している。図１５では初めにステップ５００において最大ＮＯX 吸収量ＮＯＸmax が判定値Ｃよりも大きいか否かが判別され、ステップ５０１ではＮＯX 吸収速度Ｋａｂが判定値Ｄよりも大きいか否かが判別され、ステップ５０２ではＮＯX 還元速度Ｋｒｅが所定値Ｅよりも大きいか否かが判別される。なおここで用いられるＮＯＸmax ，Ｋａｂ，Ｋｒｅは図８（Ａ）または図８（Ｄ）に示した関係式にて用いられるパラメータである。
【００８７】
ステップ５００においてＮＯＸmax ＞Ｃと判別され、且つステップ５０１においてＫａｂ＞Ｄと判別され、且つステップ５０２においてＫｒｅ＞Ｅと判別されたときにはステップ５０３に進んでＮＯX 吸収剤２３が劣化していることを表す劣化フラグがオフとされる。一方、ステップ５００においてＮＯＸmax ≦Ｃと判別されるか、或いはステップ５０１においてＫａｂ≦Ｄと判別されるか、或いはステップ５０２においてＫｒｅ≦Ｅと判別されるかしたときにはステップ５０４に進んで劣化フラグがオンされる。
【００８８】
なお上述した実施例では最大ＮＯX 吸収量ＮＯＸmax ，ＮＯX 吸収速度Ｋａｂ，ＮＯX 還元速度Ｋｒｅといったパラメータの初期値として未使用状態のＮＯX 吸収剤における値を採用し、使用が開始されてから修正を重ねるようにしているがＮＯＸmax ，Ｋａｂ，ＫｒｅはＮＯX 吸収剤２３の温度によっても変化するので上述した実施例においてＮＯX 吸収剤の温度を加味するようにしてもよい。具体的には例えばＮＯX 吸収剤２３の温度に応じて定まる各パラメータの初期値を予め求めてＲＯＭ３４に記憶しておき、上述した各パラメータに対する修正に関しては各パラメータに対する修正係数を算出し、単位ＮＯX 吸収量Ａａｂまたは単位ＮＯX 還元量Ａｒｅを算出するときにＮＯX 吸収剤２３の温度により算出される各パラメータの初期値をこの修正係数により修正し、斯くして修正したパラメータを図８（Ａ）または図（Ｄ）に示した関係式において使用するようにすればよい。
【００８９】
また上述した実施例ではＮＯX 吸収剤２３に流入する排気ガス中のＮＯX 濃度Ｃnox を機関回転数Ｎと、機関負荷Ｌと、吸入空気量Ｑとから算出するようにしているがＮＯX 吸収剤２３上流側にＮＯX アンモニアセンサを配置し、このＮＯX アンモニアセンサによりＮＯX 濃度Ｃnox を直接検出するようにしてもよい。またＮＯX 吸収剤２３のＮＯX 還元能力が不十分である場合にはＮＯX アンモニアセンサ２９下流にＮＯX 還元能力が高い触媒を配置するようにしてもよい。
【００９０】
また修正されたパラメータ（すなわち最大ＮＯX 吸収量、ＮＯX 吸収速度）を判定値と比較するのではなく、初期値をそのまま利用して図８（Ａ）の関係式を用いて算出した単位ＮＯX 吸収量に基づいてＮＯX 吸収剤下流に流出するＮＯX の量を算出し、斯くして算出したＮＯX 量とＮＯX アンモニアセンサ２９の出力から算出したＮＯX 吸収剤下流に流出するＮＯX の量とを比較し、両者の間に大きく隔りがあるときにＮＯX 吸収剤２３が劣化したと判定するようにしてもよい。
【００９１】
【発明の効果】
１番目の発明によれば排気ガス中のＮＯX 濃度を検出することができるセンサの出力を利用してＮＯX 吸収剤の最大ＮＯX 吸収量、またはＮＯX 吸収剤へのＮＯX の吸収速度が算出され、これら最大ＮＯX 吸収量またはＮＯX 吸収速度に基づいてＮＯX 吸収剤の劣化が判定される。このように排気ガス中のＮＯX 濃度の実測値、すなわち推定値ではない値を用いて算出される最大ＮＯX 吸収量またはＮＯX 吸収速度といったパラメータを利用するのでＮＯX 吸収剤の劣化を正確に判定することができる。
【００９２】
４番目の発明によれば排気ガス中のアンモニア濃度を検出しうるセンサの出力に基づいて決定されるＮＯX吸収剤におけるＮＯXの放出・還元の完了タイミングにおいてＮＯX吸収量推定手段により推定されるＮＯX吸収量に基づいて補正されるＮＯX吸収剤におけるＮＯX 還元速度に基づいてＮＯX吸収剤の劣化が判定される。このように排気ガス中のアンモニア濃度の実測値、すなわち推定値ではない値を用いて算出されるＮＯX還元速度といったパラメータを利用するのでＮＯX吸収剤の劣化を正確に判定することができる。
【図面の簡単な説明】
【図１】内燃機関の全体図である。
【図２】ＮＯX アンモニアセンサのセンサ部の構造を示す図である。
【図３】ＮＯX アンモニアセンサによる検出電流を示す図である。
【図４】空燃比センサの出力電圧を示す図である。
【図５】基本燃料噴射量、補正係数等を示す図である。
【図６】ＮＯX 吸収剤のＮＯX 吸放出作用を説明するための図である。
【図７】空燃比センサの出力電圧、ＮＯX アンモニアセンサの検出電流等を示すタイムチャートである。
【図８】単位ＮＯX 吸収量等を算出するための関係式を示す図である。
【図９】ＮＯX 流入量のマップを示す図である。
【図１０】機関運転を制御するためのフローチャートである。
【図１１】総ＮＯX 吸収量を算出するためのフローチャートである。
【図１２】総ＮＯX 吸収量を算出するためのフローチャートである。
【図１３】目標レベルＥS を算出するためのフローチャートである。
【図１４】目標レベルＥS を算出するためのフローチャートである。
【図１５】ＮＯX 吸収剤の劣化を判定するためのフローチャートである。
【符号の説明】
１１…燃料噴射弁
２３…ＮＯX 吸収剤
２９…ＮＯX アンモニアセンサ
３０…空燃比センサ
【図１】
0003775229.tif 000002
【図２】
0003775229.tif 000003
【図３】
0003775229.tif 000004
【図４】
0003775229.tif 000005
【図５】
0003775229.tif 000006
【図６】
0003775229.tif 000007
【図７】
0003775229.tif 000008
【図８】
0003775229.tif 000009
【図９】
0003775229.tif 000010
【図１０】
0003775229.tif 000011
【図１１】
0003775229.tif 000012
【図１２】
0003775229.tif 000013
【図１３】
0003775229.tif 000014
【図１４】
0003775229.tif 000015
【図１５】
0003775229.tif 000016