(57)【特許請求の範囲】
【請求項１】
内燃機関の排気系に設けられ、排気リーン状態で排気中のＮＯｘを吸収するＮＯｘ浄化手段と、該ＮＯｘ浄化手段に吸収された硫黄酸化物量を推定する硫黄酸化物量推定手段と、該硫黄酸化物量推定手段により推定された硫黄酸化物量が設定値に達したときに、前記硫黄酸化物を除去する処理を行う硫黄酸化物除去手段とを備える内燃機関の排気浄化装置において、
前記硫黄酸化物量推定手段は、前記機関に供給する混合気の空燃比及び前記機関の運転状態に応じて、前記硫黄酸化物量の単位時間当たりの変化量を推定し、該推定した変化量を積算することにより前記硫黄酸化物量の推定を行うものであり、
前記空燃比が理論空燃比近傍に設定されているとき、前記機関運転状態に応じて前記変化量を推定する第１推定手段と、前記空燃比が理論空燃比よりリーン側に設定されているとき、前記機関運転状態に応じて前記変化量を推定する第２推定手段と、前記空燃比が理論空燃比よりリッチ側に設定されているとき、前記機関運転状態に応じて前記変化量を推定する第３推定手段とを有し、
前記第１推定手段は、前記ＮＯｘ浄化手段の温度が第１所定温度以上となる機関運転状態においては第１の負の変化量を出力し、前記ＮＯｘ浄化手段の温度が前記第１所定温度より低くなる機関運転状態においては第１の正の変化量を出力し、
前記第２推定手段は、前記機関運転状態に応じて第２の正の変化量を出力し、
前記第３推定手段は、前記ＮＯｘ浄化手段の温度が前記第１所定温度より低い第２所定温度以上となる機関運転状態においては第２の負の変化量を出力し、前記ＮＯｘ浄化手段の温度が前記第２所定温度より低くなる機関運転状態においては第３の正の変化量を出力することを特徴とする内燃機関の排気浄化装置。
【請求項２】
前記第１推定手段は、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第１の負の変化量の絶対値が増加するように、前記第１の負の変化量を算出し、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第１の正の変化量が減少するように、前記第１の正の変化量を算出し、前記第２推定手段は、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第２の正の変化量が増加するように、前記第２の正の変化量を算出し、前記第３推定手段は、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第２の負の変化量の絶対値が増加するように、前記第２の負の変化量を算出し、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第３の正の変化量が減少するように、前記第３の正の変化量を算出することを特徴とする請求項１に記載の内燃機関の排気浄化装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、排気系に窒素酸化物（ＮＯｘ）浄化装置を備えた内燃機関の排気浄化装置に関し、特にＮＯｘ浄化装置が硫黄被毒により劣化した場合に再生処理を行うものに関する。
【０００２】
【従来の技術】
内燃機関の排気系にＮＯｘ浄化装置を設けて、空燃比を理論空燃比よりリーン側に設定するリーン運転中においてＮＯｘを吸収させるとともに、適宜空燃比のリッチ化を実行して、ＮＯｘ浄化装置に吸収されたＮＯｘを還元して放出させるようにした排気浄化装置は、従来より知られている。また、このＮＯｘ浄化装置は、燃料に含まれる硫黄の酸化物を吸収すること、すなわち硫黄被毒によって、ＮＯｘの吸収能力が低下するため、ＮＯｘ浄化装置の上流側に硫黄酸化物を吸収する能力を有するＳＯｘ触媒を配置した排気浄化装置が知られている（例えば特開平１１−２４７６５０号公報）。
【０００３】
この公報に示された装置では、ＳＯｘ触媒に吸収された硫黄酸化物量を推定し、該推定した硫黄酸化物量が設定値に達すると、その硫黄酸化物を除去する再生処理を実行するために、吸収された硫黄酸化物量を推定するカウンタを設け、このカウンタの単位時間当たりの加算値を機関回転数及び吸気管内圧に応じて設定して、カウンタをインクリメントすることにより、ＳＯｘ触媒に吸収された硫黄酸化物量を推定している。
【０００４】
【発明が解決しようとする課題】
排気系にＳＯｘ触媒を設けない場合には、硫黄酸化物はＮＯｘ浄化装置に吸収されるので、ＮＯｘ浄化装置の温度が高い状態で空燃比をリッチ化することにより、硫黄酸化物を放出させる再生処理を行う必要がある。その場合には、先ずＮＯｘ浄化装置に吸収された硫黄酸化物量を推定する必要があるが、上記公報に示された手法をそのまま適用すると以下のような問題があった。
【０００５】
すなわち、機関に供給する混合気の空燃比を理論空燃比よりリーン側に設定するリーン運転を多用する場合でも、空燃比を理論空燃比に設定するストイキ運転あるいは理論空燃比よりリッチ側に設定するリッチ運転を実行する期間もある。そのようなストイキ運転またはリッチ運転中において、ＮＯｘ浄化装置の温度が高くなったときには、ＮＯｘ浄化装置に吸収された硫黄酸化物が放出されることになる。ところが、上記公報に示された従来の装置では、通常のストイキ運転またはリッチ運転中におけるそのような硫黄酸化物の放出が考慮されていないため、ＮＯｘ浄化装置に吸収された硫黄酸化物量の推定が不正確となり、再生処理の実行時期が最適な時期からずれるという問題があった。
【０００６】
本発明はこの点に着目してなされたものであり、ＮＯｘ浄化装置に吸収される硫黄酸化物量をより正確に推定し、硫黄被毒したＮＯｘ浄化装置の再生処理を最適な時期に実行できるようにした排気浄化装置を提供することを目的とする。
【０００７】
【課題を解決するための手段】
上記目的を達成するため請求項１に記載の発明は、燃機関の排気系に設けられ、排気リーン状態で排気中のＮＯｘを吸収するＮＯｘ浄化手段と、該ＮＯｘ浄化手段に吸収された硫黄酸化物量を推定する硫黄酸化物量推定手段と、該硫黄酸化物量推定手段により推定された硫黄酸化物量が設定値に達したときに、前記硫黄酸化物を除去する処理を行う硫黄酸化物除去手段とを備える内燃機関の排気浄化装置において、前記硫黄酸化物量推定手段は、前記機関に供給する混合気の空燃比及び前記機関の運転状態に応じて、前記硫黄酸化物量の単位時間当たりの変化量を推定し、該推定した変化量を積算することにより前記硫黄酸化物量の推定を行うものであり、前記空燃比が理論空燃比近傍に設定されているとき、前記機関運転状態に応じて前記変化量を推定する第１推定手段と、前記空燃比が理論空燃比よりリーン側に設定されているとき、前記機関運転状態に応じて前記変化量を推定する第２推定手段と、前記空燃比が理論空燃比よりリッチ側に設定されているとき、前記機関運転状態に応じて前記変化量を推定する第３推定手段とを有し、前記第１推定手段は、前記ＮＯｘ浄化手段の温度が第１所定温度（７００℃）以上となる機関運転状態においては第１の負の変化量（ＳＤＥＳＳＴ）を出力し、前記ＮＯｘ浄化手段の温度が前記第１所定温度より低くなる機関運転状態においては第１の正の変化量（ＳＡＢＳＳＴ）を出力し、前記第２推定手段は、前記機関運転状態に応じて第２の正の変化量（ＳＡＢＳＬ）を出力し、前記第３推定手段は、前記ＮＯｘ浄化手段の温度が前記第１所定温度より低い第２所定温度（６００℃）以上となる機関運転状態においては第２の負の変化量（ＳＤＥＳＲ）を出力し、前記ＮＯｘ浄化手段の温度が前記第２所定温度より低くなる機関運転状態においては第３の正の変化量（ＳＡＢＳＲ）を出力することを特徴とする。
【０００８】
この構成によれば、機関に供給する混合気の空燃比及び機関の運転状態に応じて、ＮＯｘ浄化手段に吸収された硫黄酸化物量の、単位時間当たりの変化量（すなわち吸収量または放出量）が推定され、該推定された変化量を積算することにより、ＮＯｘ浄化手段に吸収された硫黄酸化物量の推定が行われる。具体的には、空燃比が理論空燃比近傍に設定されているときは、第１推定手段により、ＮＯｘ浄化手段の温度が第１所定温度以上となる機関運転状態においては第１の負の変化量が算出されるとともに、ＮＯｘ浄化手段の温度が第１所定温度より低くなる機関運転状態においては第１の正の変化量が算出される。また空燃比が理論空燃比よりリーン側に設定されているときは、第２推定手段により、機関運転状態に応じて第２の正の変化量が算出される。また空燃比が理論空燃比よりリッチ側に設定されているときは、第３推定手段により、ＮＯｘ浄化手段の温度が第１所定温度より低い第２所定温度以上となる機関運転状態においては第２の負の変化量が算出されるとともに、ＮＯｘ浄化手段の温度が第２所定温度より低くなる機関運転状態においては第３の正の変化量が算出される。ＮＯｘ浄化手段に吸収される硫黄酸化物量は、空燃比が理論空燃比よりリーン側に設定されているとき増加し、空燃比が理論空燃比近傍あるいは理論空燃比よりリッチ側に設定されているときは、機関運転状態により増加または減少するので、第１〜第３推定手段により、硫黄酸化物量の変化量を推定することにより、硫黄酸化物量の単位時間当たりの変化量、すなわち吸収量または放出量を正確に推定することができる。その結果、硫黄酸化物量の推定を従来に比べてより正確に行うことができ、硫黄被毒したＮＯｘ浄化手段の再生処理を最適な時期に実行することができる。
【０００９】
請求項２に記載の発明は、請求項１に記載の内燃機関の排気浄化装置において、前記第１推定手段は、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第１の負の変化量の絶対値が増加するように、前記第１の負の変化量を算出し、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第１の正の変化量が減少するように、前記第１の正の変化量を算出し、前記第２推定手段は、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第２の正の変化量が増加するように、前記第２の正の変化量を算出し、前記第３推定手段は、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第２の負の変化量の絶対値が増加するように、前記第２の負の変化量を算出し、前記機関の回転速度及び／または前記機関の吸気圧が高くなるほど、前記第３の正の変化量が減少するように、前記第３の正の変化量を算出することを特徴とする。
【００１０】
この構成によれば、機関の回転速度及び／または機関の吸気圧が高くなるほど、第１の負の変化量の絶対値が増加するように、第１の負の変化量が算出され、機関の回転速度及び／または機関の吸気圧が高くなるほど、第１の正の変化量が減少するように、第１の正の変化量が算出され、機関の回転速度及び／または機関の吸気圧が高くなるほど、第２の正の変化量が増加するように、第２の正の変化量が算出され、機関の回転速度及び／または機関の吸気圧が高くなるほど、第２の負の変化量の絶対値が増加するように、第２の負の変化量が算出され、機関の回転速度及び／または機関の吸気圧が高くなるほど、第３の正の変化量が減少するように、第３の正の変化量が算出される。
【００１４】
また前記硫黄酸化物除去手段は、前記硫黄酸化物を除去する処理を実行するときは、所定時間に亘って前記機関に供給する混合気の空燃比を理論空燃比近傍に設定し、その後前記空燃比を理論空燃比よりリッチ側の空燃比に設定することが望ましい。
【００１５】
また前記硫黄酸化物除去手段は、前記硫黄酸化物を除去する処理を実行するときは、前記機関の点火時期を通常設定値より遅角させ、排気還流を停止するとともに、前記機関の出力トルクが変化しないように前記機関の吸入空気量を制御することが望ましい。
【００１６】
【発明の実施の形態】
以下本発明の実施の形態を図面を参照して説明する。
図１は、本発明の実施の一形態にかかる排気浄化装置を含む、内燃機関（以下「エンジン」という）及びその制御装置の全体構成図であり、例えば４気筒のエンジン１の吸気管２の途中にはスロットル弁３が配されている。スロットル弁３にはスロットル弁開度（θＴＨ）センサ４が連結されており、当該スロットル弁３の開度に応じた電気信号を出力してエンジン制御用電子コントロールユニット（以下「ＥＣＵ」という）５に供給する。
【００１７】
吸気管２にはスロットル弁３をバイパスする補助空気通路１７が接続されており、補助空気通路１７の途中には補助空気量を制御する補助空気制御弁１８が設けられている。補助空気制御弁１８は、ＥＣＵ５に接続されており、ＥＣＵ５によりその開弁量が制御される。
【００１８】
燃料噴射弁６はエンジン１とスロットル弁３との間かつ吸気管２の図示しない吸気弁の少し上流側に各気筒毎に設けられており、各噴射弁は図示しない燃料ポンプに接続されていると共にＥＣＵ５に電気的に接続されて当該ＥＣＵ５からの信号により燃料噴射弁６の開弁時間が制御される。
【００１９】
一方、スロットル弁３の直ぐ下流には吸気管内絶対圧（ＰＢＡ）センサ８が設けられており、この絶対圧センサ８により電気信号に変換された絶対圧信号は前記ＥＣＵ５に供給される。また、その下流には吸気温（ＴＡ)センサ９が取付けられており、吸気温ＴＡを検出して対応する電気信号を出力してＥＣＵ５に供給する。
【００２０】
エンジン１の本体に装着されたエンジン水温（ＴＷ）センサ１０はサーミスタ等から成り、エンジン水温（冷却水温）ＴＷを検出して対応する温度信号を出力してＥＣＵ５に供給する。
ＥＣＵ５には、エンジン１のクランク軸（図示せず）の回転角度を検出するクランク角度位置センサ１０が接続されており、クランク軸の回転角度に応じた信号がＥＣＵ５に供給される。クランク角度位置センサ１０は、エンジン１の特定の気筒の所定クランク角度位置で信号パルス（以下「ＣＹＬ信号パルス」という）を出力する気筒判別センサ、各気筒の吸入行程開始時の上死点（ＴＤＣ）より所定クランク角度前のクランク角度位置で（４気筒エンジンではクランク角１８０度毎に）ＴＤＣ信号パルスを出力するＴＤＣセンサ及びＴＤＣ信号パルスより短い一定クランク角周期（例えば３０度周期）で１パルス（以下「ＣＲＫ信号パルス」という）を発生するＣＲＫセンサから成り、ＣＹＬ信号パルス、ＴＤＣ信号パルス及びＣＲＫ信号パルスがＥＣＵ５に供給される。これらの信号パルスは、燃料噴射時期、点火時期等の各種タイミング制御及びエンジン回転数（エンジン回転速度）ＮＥの検出に使用される。
【００２１】
排気管１２に、ＮＯｘ浄化手段としてのＮＯｘ浄化装置１５が設けられている。ＮＯｘ浄化装置１５は、ＮＯｘを吸収するＮＯｘ吸収剤及び酸化、還元を促進するための触媒を内蔵する。ＮＯｘ浄化装置１５は、エンジン１に供給される混合気の空燃比が理論空燃比よりリーン側に設定され、排気中の酸素濃度が比較的高い状態（以下「排気リーン状態」という）では、ＮＯｘを吸収する一方、エンジン１に供給される混合気の空燃比が理論空燃比近傍または理論空燃比よりリッチ側に設定され、排気中の酸素濃度が比較的低く、ＨＣ、ＣＯ成分が多い状態（以下「排気リッチ状態」という）では、吸収されたＮＯｘがＨＣ、ＣＯにより還元されて、窒素ガスとして排出され、またＨＣ、ＣＯは酸化されて水蒸気及び二酸化炭素として排出されるように構成されている。
【００２２】
ＮＯｘ吸収剤のＮＯｘ吸収能力の限界、すなわち最大ＮＯｘ吸収量まで、ＮＯｘを吸収すると、それ以上ＮＯｘを吸収できなくなるので、適時ＮＯｘを放出させて還元するために空燃比のリッチ化、すなわち還元リッチ化を実行する。
またＮＯｘ吸収剤に硫黄酸化物（以下「ＳＯｘ」という）が吸収され、その吸収されたＳＯｘ量が設定値に達したときは、ＳＯｘを除去（パージ）する再生処理を実行する。なお、本実施形態では、ＮＯｘを吸着するタイプのＮＯｘ吸収剤を用いており、ＮＯｘ吸収剤によるＳＯｘの「吸収」を「吸着」ともいう。また、ＮＯｘ吸収剤からのＳＯｘの「放出」を「脱離」ともいう。
【００２３】
ＮＯｘ浄化装置１５の上流位置には、比例型空燃比センサ１４（以下「ＬＡＦセンサ１４」という）が装着されており、このＬＡＦセンサ１４は排気中の酸素濃度（空燃比）にほぼ比例した電気信号を出力し、ＥＣＵ５に供給する。
エンジン１の各気筒毎に設けられた点火プラグ１１は、ＥＣＵ５に接続されており、点火プラグ１１の駆動信号、すなわち点火信号がＥＣＵ５から供給される。
【００２４】
吸気管２のスロットル弁３の下流側と、排気管１２のＮＯｘ浄化装置１５の上流側との間には、排気還流通路２１が設けられており、排気還流通路２１の途中には排気還流量を制御する排気還流弁（以下「ＥＧＲ弁」という）２２が設けられている。ＥＧＲ弁２２は、ソレノイドを有する電磁弁であり、その弁開度はＥＣＵ５により制御される。ＥＧＲ弁２２には、その弁開度（弁リフト量）ＬＡＣＴを検出するリフトセンサ２３が設けられており、その検出信号はＥＣＵ５に供給される。排気還流通路２１及びＥＧＲ弁２２より、排気還流機構が構成される。
【００２５】
エンジン１は、吸気弁及び排気弁のバルブタイミングを、エンジンの高速回転領域に適した高速バルブタイミングと、低速回転領域に適した低速バルブタイミングとの２段階に切換可能なバルブタイミング切換機構３０を有する。このバルブタイミングの切換は、弁リフト量の切換も含み、さらに低速バルブタイミング選択時は２つに吸気弁のうちの一方を休止させて、空燃比を理論空燃比よりリーン化する場合においても安定した燃焼を確保するようにしている。
【００２６】
バルブタイミング切換機構３０は、バルブタイミングの切換を油圧を介して行うものであり、この油圧切換を行う電磁弁及び油圧センサがＥＣＵ５に接続されている。油圧センサの検出信号はＥＣＵ５に供給され、ＥＣＵ５は電磁弁を制御してエンジン１の運転状態に応じたバルブタイミングの切換制御を行う。
【００２７】
ＥＣＵ５には、エンジン１によって駆動される車両の走行速度（車速）ＶＰを検出する車速センサ３１が接続されており、その検出信号がＥＣＵ５に供給される。
ＥＣＵ５は、各種センサからの入力信号波形を整形し、電圧レベルを所定レベルに修正し、アナログ信号値をデジタル信号値に変換する等の機能を有する入力回路５ａ、中央演算処理ユニット（以下「ＣＰＵ」という）５ｂ、ＣＰＵ５ｂで実行される各種演算プログラム及び演算結果等を記憶する記憶手段５ｃ、前記燃料噴射弁６に駆動信号を供給する出力回路５ｄ等から構成される。
【００２８】
ＥＣＵ５は、各種エンジンパラメータ信号に基づいてエンジン運転状態を判別し、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて設定されるＥＧＲ弁２２の弁開度指令値ＬＣＭＤと、リフトセンサ２３によって検出される実弁開度ＬＡＣＴとの偏差を零にするようにＥＧＲ弁２２のソレノイドに制御信号を供給する。
【００２９】
ＣＰＵ５ｂは、上述の各種エンジンパラメータ信号に基づいて、種々のエンジン運転状態を判別するとともに、該判別されたエンジン運転状態に応じて、次式（１）に基づき、前記ＴＤＣ信号パルスに同期して開弁作動する燃料噴射弁６の燃料噴射時間ＴＯＵＴを演算する。
ＴＯＵＴ＝ＴＩＭ×ＫＣＭＤ×ＫＬＡＦ×Ｋ１＋Ｋ２…（１）
【００３０】
ここに、ＴＩＭは基本燃料量、具体的には燃料噴射弁６の基本燃料噴射時間であり、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて設定されたＴＩマップを検索して決定される。ＴＩマップは、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに対応する運転状態において、エンジンに供給する混合気の空燃比がほぼ理論空燃比になるように設定されている。すなわち、基本燃料量ＴＩＭは、エンジンの単位時間当たりの吸入空気量（質量流量）にほぼ比例する値を有する。
【００３１】
ＫＣＭＤは目標空燃比係数であり、エンジン回転数ＮＥ、スロットル弁開度θＴＨ、エンジン水温ＴＷ等のエンジン運転パラメータに応じて設定される。目標空燃比係数ＫＣＭＤは、空燃比Ａ／Ｆの逆数、すなわち燃空比Ｆ／Ａに比例し、理論空燃比のとき値１．０をとるので、目標当量比ともいう。また目標空燃比係数ＫＣＭＤは、後述するように還元リッチ化あるいはＳＯｘ除去のためのリッチ化（以下「ＳＯｘ除去リッチ化」という）を実行するときは、１．０より大きな値に設定される。
【００３２】
ＫＬＡＦは、フィードバック制御の実行条件が成立するときは、ＬＡＦセンサ１４の検出値から算出される検出当量比ＫＡＣＴが目標当量比ＫＣＭＤに一致するようにＰＩＤ（比例、積分、微分）制御により算出される空燃比補正係数である。
【００３３】
Ｋ１及びＫ２は夫々各種エンジンパラメータ信号に応じて演算される他の補正係数および補正変数であり、エンジン運転状態に応じた燃費特性、エンジン加速特性等の諸特性の最適化が図れるような所定値に決定される。
ＣＰＵ５ｂは上述のようにして求めた燃料噴射時間ＴＯＵＴに基づいて燃料噴射弁６を開弁させる駆動信号、点火プラグ１１を駆動する点火信号、補助空気量制御弁１８の駆動信号、及びＥＧＲ弁２２の駆動信号を出力回路５ｄを介して出力する。
【００３４】
図２は、ＮＯｘ浄化装置１５のＮＯｘ吸収剤に吸着されたＳＯｘ量（以下「ＳＯｘ吸着量」という）が設定値に達したときに、その再生制御を行う処理のフローチャートである。本処理は、ＴＤＣ信号パルスの発生に同期して、ＣＰＵ５ｂで実行される。
【００３５】
ステップＳ１０では、図３に示すＳＯｘ飽和判断処理を実行する。この処理では、設定空燃比（目標空燃比係数ＫＣＭＤ）及びエンジン１の運転状態に応じて、ＳＯｘ吸着量が算出され、そのＳＯｘ吸収量が設定値に達し、且つエンジン１が再生処理を実行するのに適した運転状態にあるとき、再生モードフラグＦＳＲＣＭＯＤＥが「１」に設定される。
【００３６】
続くステップＳ１１では、再生モードフラグＦＳＲＣＭＯＤＥが「１」であるか否かを判別し、ＦＳＲＣＭＯＤＥ＝０であってＳＯｘ吸着量が設定値に達していないときは、通常の燃料供給制御、点火時期制御、補助空気量制御及び排気還流制御を実行する（ステップＳ１２）。
【００３７】
ステップＳ１０の処理において再生モードフラグＦＳＲＣＭＯＤＥが「１」に設定されると、ステップＳ１１からＳ１３に進み、ＮＯｘ浄化装置１５の昇温を促進する昇温モード制御を実行する。この昇温モード制御では、空燃比が理論空燃比となるように燃料供給制御を行い、点火時期を最適点火時期よりリタードさせる点火時期制御を行い、さらにそのような燃料供給制御及び点火時期制御に対応してエンジン１の出力トルクが通常制御時とほぼ同一となるように補助空気量制御を行うとともに、ＥＧＲ弁２２を閉弁して排気還流を停止させる制御を行う（ステップＳ１３）。このような昇温モード制御により、エンジンの出力トルクの変動を抑制しつつ、ＮＯｘ浄化装置１５の昇温を促進することできる。
【００３８】
ステップＳ１３における昇温モード制御に対応する目標空燃比係数ＫＣＭＤ設定処理が図９に示されている。この処理では、後述するように、目標空燃比係数ＫＣＭＤを「１．０」に維持する制御を所定時間ＴＭＴＣＡＴ実行した時点で、昇温が完了したことを「１」で示す昇温完了フラグＦＴＣＡＴＯＫが「１」に設定される。
【００３９】
ステップＳ１４では、昇温完了フラグＦＴＣＡＴＯＫが「１」であるか否かを判別し、ＦＴＣＡＴＯＫ＝０である間は、直ちに本処理を終了する。ＦＴＣＡＴＯＫ＝１となると、ＳＯｘ除去（パージ）モード処理を実行する（ステップＳ１５）。この処理では、空燃比を理論空燃比よりリッチ側とする燃料供給制御を実行し、点火時期を最適点火時期よりリタードさせる点火時期制御を実行し、さらにそのような燃料供給制御及び点火時期制御に対応してエンジン１の出力トルクが通常制御時とほぼ同一となるように補助空気量制御を行うとともに、排気還流の停止を継続する。このようなＳＯｘ除去モード処理を実行することにより、エンジンの出力トルクの変動を抑制しつつ、ＮＯｘ浄化装置１５に吸収されたＳＯｘを還元して、排気とともに放出させ、ＮＯｘ浄化装置１５のＮＯｘ吸収能力を回復させることができる。
【００４０】
図３は、図２のステップＳ１０で実行されるＳＯｘ飽和判断処理のフローチャートである。
ステップＳ２０では、図４及び５に示すＳＯｘ吸着量推定処理を実行する。この処理では、設定空燃比（目標空燃比係数ＫＣＭＤ）及びエンジン１の運転状態に応じて、ＳＯｘ吸着量推定値が算出される。より具体的には、ＳＯｘ吸着量カウンタＣＳＡＢＳの値をインクリメントまたはデクリメントすることにより、ＳＯｘ吸着量推定値が算出される。カウンタＣＳＡＢＳの値は、イグニッションスイッチがオフされたときも記憶内容が保持されるバックアップメモリに記憶される。
【００４１】
ステップＳ２１では、ＳＯｘ吸着量カウンタＣＳＡＢＳの値が設定値ＣＳＲＭＡＣ１以上か否かを判別し、ＣＳＡＢＳ＜ＣＳＲＭＡＣ１であるときは直ちに本処理を終了する。
ＣＳＡＢＳ≧ＣＳＲＭＡＣ１となると、ステップＳ２１からステップＳ２２に進み、エンジン運転状態及び車速ＶＰが所定の条件を満たすか否かを判別する。すなわち、エンジン回転数ＮＥが所定上下限値ＮＥＳＲＭＨ，ＮＥＳＲＭＬ（例えば３０００ｒｐｍ，１５００ｒｐｍ）の範囲内にあるか否か、吸気管内絶対圧ＰＢＡが所定上下限値ＰＢＳＲＭＨ，ＰＢＳＲＭＬ（例えば８１．３ｋＰａ，６１．３ｋＰａ）の範囲内にあるか否か、吸気温ＴＡが所定上下限値ＴＡＳＲＭＨ，ＴＡＳＲＭＬ（例えば１００℃，０℃）の範囲内にあるか否か、エンジン水温ＴＷが所定上下限値ＴＷＳＲＭＨ，ＴＷＳＲＭＬ（例えば１００℃，８０℃）の範囲内にあるか否か、車速ＶＰが所定上下限値ＶＳＲＭＨ，ＶＳＲＭＬ（例えば１２０ｋｍ／ｈ，６０ｋｍ／ｈ）の範囲内にあるか否かを判別する。そして、これら判別に対する答がすべて肯定（ＹＥＳ）のときは、再生モードフラグＦＳＲＣＭＯＤＥを「１」に設定する（ステップＳ２３）一方、いずれかの答が否定（ＮＯ）のときは、再生モードフラグＦＳＲＣＭＯＤＥを「０」に設定する（ステップＳ２４）。
【００４２】
図４及び５は、図３のステップＳ２０で実行されるＳＯｘ吸着量推定処理のフローチャートである。
ステップＳ３０では、始動モード、すなわちクランキング中か否かを判別し、始動モードでないときは、燃料供給遮断運転中であることを「１」で示すフュエルカットフラグＦＦＣが「１」であるか否かを判別する（ステップＳ３１）。そして、始動モードであるときまたは燃料供給遮断運転中であるときは、直ちに本処理を終了する。
【００４３】
燃料供給遮断運転中でないときは、さらにアイドル運転中であることを「１」で示すアイドルフラグＦＩＤＬＥが「１」であるか否かを判別し（ステップＳ３２）、ＦＩＤＬＥ＝１であってアイドル運転中であるときは、加算値ＳＡＢＳをアイドル用加算値ＳＡＢＳＩＤＬに設定する（ステップＳ３３）。そして、ＳＯｘ吸着量カウンタＣＳＡＢＳを加算値ＳＡＢＳだけインクリメントし（ステップＳ６１）、本処理を終了する。
【００４４】
ＦＩＤＬＥ＝０であってアイドル運転中でないときは、目標空燃比係数ＫＣＭＤがリーン判定閾値ＫＢＳＤＥＳＬ（例えば０．９）より小さいか否かを判別し（ステップＳ３４）、ＫＣＭＤ＜ＫＢＳＤＥＳＬであってリーン運転中であるときは、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて図６（ａ）に示すリーン領域吸着マップを検索し、リーン運転用加算値ＳＡＢＳＬを算出する（ステップＳ３６）。リーン領域吸着マップは、同図にハッチングを付して示す、リーン運転が実行される領域について設定されており、エンジン回転数ＮＥが増加するほど、また吸気管内絶対圧ＰＢＡが増加するほど、加算値ＳＡＢＳＬが増加するように設定されている。
【００４５】
次いで加算値ＳＡＢＳを、このリーン運転用加算値ＳＡＢＳＬに設定して（ステップＳ３７）、前記ステップＳ６１に進む。
一方ステップＳ３４でＫＣＭＤ≧ＫＢＳＤＥＳＬであるときは、さらに目標空燃比係数ＫＣＭＤがリッチ判定閾値ＫＢＳＤＥＳＨ（例えば１．１）より大きいか否かを判別する（ステップＳ３８）。ＫＣＭＤ≦ＫＢＳＤＥＳＨであって空燃比を理論空燃比に設定するストイキ運転中であるときは、エンジン回転数ＮＥに応じて図７（ａ）に示すストイキ境界テーブルを検索し、ストイキ境界圧ＰＢＳＴＧを算出する（ステップＳ３９）。このストイキ境界テーブルは、ＮＯｘ浄化装置１５の温度がほぼ７００℃となるようなエンジン回転数ＮＥと吸気管内絶対圧ＰＢＡの組が得られるように設定されている。このテーブルで定義されるラインＬＳＴより吸気管内絶対圧ＰＢＡが高いときは、ＮＯｘ浄化装置１５の温度は、７００℃より高くなり、ＮＯｘ吸収剤に吸着されたＳＯｘが脱離されるが、逆にラインＬＳＴより吸気管内絶対圧ＰＢＡが低いときは、ＮＯｘ浄化装置１５の温度は、７００℃より低くなり、ＮＯｘ吸収剤にＳＯｘが吸着される。そこで、以下のステップＳ４０〜Ｓ４８の処理では、ＳＯｘが吸着される運転状態では、ＳＯｘ吸着量カウンタＣＳＡＢＳの加算値ＳＡＢＳを算出する一方、ＳＯｘが脱離される運転状態では、ＳＯｘ吸着量カウンタＣＳＡＢＳの減算値ＳＤＥＳを算出する。
【００４６】
先ずステップＳ４０では、吸気管内絶対圧ＰＢＡがストイキ境界圧ＰＢＳＴＧより低いか否かを判別し、ＰＢＡ＜ＰＢＳＴＧであるときは、後述するステップＳ４６で参照されるダウンカウントタイマｔＤＥＳＴＬＨＤを所定遅延時間ＴＤＥＳＴＬＨＤ（例えば３秒）にセットしてスタートさせる（ステップＳ４１）。次いで後述するステップＳ４５でスタートされたダウンカウントタイマｔＤＥＳＴＨＬＤの値が「０」であるか否かを判別し（ステップＳ４２）、ｔＤＥＳＴＨＬＤ＞０である間はステップＳ４７に進み、ｔＤＥＳＴＨＬＤ＝０となると、ステップＳ４３に進んで、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて、図６（ｂ）に示すストイキ領域吸着マップを検索し、ストイキ運転用加算値ＳＡＢＳＳＴを算出する。ストイキ領域吸着マップは、同図に右上がりのハッチングを付して示す、ＳＯｘがＮＯｘ吸収剤に吸着される領域について設定されており、エンジン回転数ＮＥが増加するほど、また吸気管内絶対圧ＰＢＡが増加するほど、加算値ＳＡＢＳＳＴが減少するように設定されている。なお、同図において破線で囲まれた領域は、図６（ａ）に示すリーン運転領域である。
【００４７】
次いで加算値ＳＡＢＳを、このストイキ運転用加算値ＳＡＢＳＳＴに設定して（ステップＳ４４）、前記ステップＳ６１に進む。
またステップＳ４０でＰＢＡ≧ＰＢＳＴＧであるときは、ダウンカウントタイマｔＤＥＳＴＨＬＤを所定遅延時間ＴＤＥＳＴＨＬＤ（例えば３秒）にセットしてスタートさせる（ステップＳ４５）。次いでステップＳ４１でスタートされたダウンカウントタイマｔＤＥＳＴＬＨＤの値が「０」であるか否かを判別し（ステップＳ４６）、ｔＤＥＳＴＬＨＤ＞０である間はステップＳ４３に進み、ｔＤＥＳＴＬＨＤ＝０となると、ステップＳ４７に進んで、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて、図６（ｂ）に示すストイキ領域脱離マップを検索し、ストイキ運転用減算値ＳＤＥＳＳＴを算出する。ストイキ領域脱離マップは、同図に右下がりのハッチングを付して示す、ＳＯｘがＮＯｘ吸収剤から脱離される領域について設定されており、エンジン回転数ＮＥが増加するほど、また吸気管内絶対圧ＰＢＡが増加するほど、減算値ＳＤＥＳＳＴが増加するように設定されている。
【００４８】
次いで減算値ＳＤＥＳを、このストイキ運転用減算値ＳＤＥＳＳＴに設定して（ステップＳ４８）、ステップＳ６２に進み、ＳＯｘ吸着量カウンタＣＳＡＢＳの値を減算値ＳＤＥＳだけデクリメントして、本処理を終了する。
ステップＳ４１，Ｓ４２，Ｓ４５及びＳ４６は、エンジン運転状態が吸着領域から脱離領域へ、またはその逆に変化したときに、制御の切換をエンジン運転状態が安定するまで遅延させるために設けられている。
【００４９】
ステップＳ３８で、ＫＣＭＤ＞ＫＢＳＤＥＳＨであって空燃比を理論空燃比よりリッチ側に設定するリッチ運転中であるときは、エンジン回転数ＮＥに応じて図７（ｂ）に示すリッチ境界テーブルを検索し、リッチ境界圧ＰＢＲＩＣＨＧを算出する（図５、ステップＳ５１）。このリッチ境界テーブルは、ＮＯｘ浄化装置１５の温度がほぼ６００℃となるようなエンジン回転数ＮＥと吸気管内絶対圧ＰＢＡの組が得られるように設定されている。このテーブルで定義されるラインＬＲより吸気管内絶対圧ＰＢＡが高いときは、ＮＯｘ浄化装置１５の温度は、６００℃より高くなり、ＮＯｘ吸収剤からＳＯｘが脱離されるが、逆にラインＬＲより吸気管内絶対圧ＰＢＡが低いときは、ＮＯｘ浄化装置１５の温度は、６００℃より低くなり、ＮＯｘ吸収剤にＳＯｘが吸着される。そこで、以下のステップＳ５２〜Ｓ６０の処理では、ＳＯｘが吸着される運転状態では、ＳＯｘ吸着量カウンタＣＳＡＢＳの加算値ＳＡＢＳを算出する一方、ＳＯｘが脱離される運転状態では、ＳＯｘ吸着量カウンタＣＳＡＢＳの減算値ＳＤＥＳを算出する。
【００５０】
先ずステップＳ５２では、吸気管内絶対圧ＰＢＡがリッチ境界圧ＰＢＲＩＣＨＧより低いか否かを判別し、ＰＢＡ＜ＰＢＲＩＣＨＧであるときは、後述するステップＳ５８で参照されるダウンカウントタイマｔＤＥＳＲＬＨＤを所定遅延時間ＴＤＥＳＲＬＨＤ（例えば３秒）にセットしてスタートさせる（ステップＳ５３）。次いで後述するステップＳ５７でスタートされたダウンカウントタイマｔＤＥＳＲＨＬＤの値が「０」であるか否かを判別し（ステップＳ５４）、ｔＤＥＳＲＨＬＤ＞０である間はステップＳ５９に進み、ｔＤＥＳＲＨＬＤ＝０となると、ステップＳ５５に進んで、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて、図６（ｃ）に示すリッチ領域吸着マップを検索し、リッチ運転用加算値ＳＡＢＳＲを算出する。リッチ領域吸着マップは、同図に右上がりのハッチングを付して示す、ＳＯｘがＮＯｘ吸収剤に吸着される領域について設定されており、エンジン回転数ＮＥが増加するほど、また吸気管内絶対圧ＰＢＡが増加するほど、加算値ＳＡＢＳＲが減少するように設定されている。なお、同図において破線で囲まれた領域は、図６（ａ）に示すリーン運転領域である。
【００５１】
次いで加算値ＳＡＢＳを、このリッチ運転用加算値ＳＡＢＳＲに設定して（ステップＳ５６）、前記ステップＳ６１に進む。
またステップＳ５２でＰＢＡ≧ＰＢＲＩＣＨＧであるときは、ダウンカウントタイマｔＤＥＳＲＨＬＤを所定遅延時間ＴＤＥＳＲＨＬＤ（例えば３秒）にセットしてスタートさせる（ステップＳ５７）。次いでステップＳ５３でスタートされたダウンカウントタイマｔＤＥＳＲＬＨＤの値が「０」であるか否かを判別し（ステップＳ５８）、ｔＤＥＳＲＬＨＤ＞０である間はステップＳ５５に進み、ｔＤＥＳＲＬＨＤ＝０となると、ステップＳ５９に進んで、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて、図６（ｃ）に示すリッチ領域脱離マップを検索し、リッチ運転用減算値ＳＤＥＳＲを算出する。リッチ領域脱離マップは、同図に右下がりのハッチングを付して示す、ＳＯｘがＮＯｘ吸収剤から脱離される領域について設定されており、エンジン回転数ＮＥが増加するほど、また吸気管内絶対圧ＰＢＡが増加するほど、減算値ＳＤＥＳＲが増加するように設定されている。
【００５２】
次いで減算値ＳＤＥＳを、このリッチ運転用減算値ＳＤＥＳＲに設定して（ステップＳ６０）、ステップＳ６２に進み、ＳＯｘ吸着量カウンタＣＳＡＢＳの値を減算値ＳＤＥＳだけデクリメントして、本処理を終了する。
以上のように図４，５の処理によれば、ＳＯｘ吸着量カウンタＣＳＡＢＳの値は、設定空燃比及びエンジン運転状態に応じた、単位時間当たりのＳＯｘ吸着量に対応する加算値ＳＡＢＳと、単位時間当たりのＳＯｘ脱離量に対応する減算値ＳＤＥＳにより、インクリメントまたはデクリメントされる。すなわち、ＮＯｘ吸収剤に吸収されたＳＯｘ量の推定値が、空燃比及びエンジン運転状態に応じて設定される加算値ＳＡＢＳ及び減算値ＳＤＥＳを積算することにより算出される。その結果、ＮＯｘ浄化装置１５のＮＯｘ吸収剤に吸着されるＳＯｘ量を正確に推定することができ、ＳＯｘ除去処理を最適な時期に実行することが可能となる。
【００５３】
図８は、目標空燃比係数ＫＣＭＤを算出する処理のフローチャートであり、本処理は、ＣＰＵ５ｂにおいてＴＤＣ信号パルスの発生に同期して実行される。
ステップＳ１０１では、再生モードフラグＦＳＲＣＭＯＤＥが「１」であるか否かを判別し、ＦＳＲＣＭＯＤＥ＝０であるときは、通常制御、すなわちエンジン運転状態に応じた目標空燃比係数ＫＣＭＤの設定を行う（ステップＳ１０３）。目標空燃比係数ＫＣＭＤは、基本的には、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて算出し、エンジン水温ＴＷの低温状態や所定の高負荷運転状態では、それらの運転状態に応じた値に変更される。次いで、後述する図９の処理で参照されるカウンタＮＦＩＳＲＭＤ１，ＮＦＩＳＲＭＤ２及びＮＦＩＳＲＭＤ３に、それぞれ所定値Ｎ１，Ｎ２及びＮ３をセットし（ステップＳ１０４）、本処理を終了する。
ＦＳＲＣＭＯＤＥ＝１であって、ＮＯｘ浄化装置１５の再生処理を実行するときは、図９に示すＮＯｘＣＡＴ再生ＫＣＭＤ算出処理を実行する（ステップＳ１０２）。
【００５４】
図９のステップＳ７０では、カウンタＮＦＩＳＲＭＤ１の値が「０」か否かを判別する。最初は、ＮＦＩＳＲＭＤ１＞０であるので、このカウンタＮＦＩＳＲＭＤ１の値を「１」だけデクリメントする（ステップＳ７１）とともに、ダウンカウントタイマｔｍＴＣＡＴを昇温モードの実行時間ＴＭＴＣＡＴ（例えば６０秒）に設定してスタートさせる（ステップＳ７２）。さらに、昇温モードが完了したことを「１」で示す昇温完了フラグＦＴＣＡＴＯＫを「０」に設定し（ステップＳ７３）、ダウンカウントタイマｔｍＳＰＲＧをＳＯｘ除去モードの実行時間ＴＭＳＰＲＧ（例えば３０秒）に設定してスタートさせ（ステップＳ７４）、本処理を終了する。
【００５５】
カウンタＮＦＩＳＲＭＤ１の値が「０」となると、ステップＳ７０からステップＳ７５に進み、昇温完了フラグＦＴＣＡＴＯＫが「１」であるか否かを判別する。最初は、ＦＴＣＡＴＯＫ＝０であるので、目標空燃比係数ＫＣＭＤを昇温モード用の値、すなわち理論空燃比に対応する「１．０」に設定し（ステップＳ７６）、タイマｔｍＴＣＡＴの値が「０」であるか否かを判別する（ステップＳ７７）。ｔｍＴＣＡＴ＞０である間は前記ステップＳ７４に進み、ｔｍＴＣＡＴ＝０となると、昇温完了フラグＦＴＣＡＴＯＫを「１」に設定する（ステップＳ７８）。昇温完了フラグＦＴＣＡＴＯＫが「１」に設定されると、以後はステップＳ７５から直ちにステップＳ８０に進む。
【００５６】
続くステップＳ８０では、カウンタＮＦＩＳＲＭＤ２の値が「０」であるか否かを判別する。最初はＮＦＩＳＲＭＤ２＞０であるので、このカウンタＮＦＩＳＲＭＤ２の値を「１」だけデクリメントして（ステップＳ７９）、前記ステップＳ７４に進む。ＮＦＩＳＲＭＤ２＝０となると、ステップＳ８０からステップＳ８１に進み、エンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じてＫＣＭＤＳＰＲＧマップを検索し、ＳＯｘ除去モード用のリッチ化設定値ＫＣＭＤＳＰＲＧを算出する。ＫＣＭＤＳＰＲＧマップは、エンジン回転数ＮＥが増加するほど、また吸気管内絶対圧ＰＢＡが増加するほど、リッチ化設定値ＫＣＭＤＳＰＲＧが増加するように設定されている。なお、ＳＯｘ除去用のリッチ化設定値ＫＣＭＤＳＰＲＧは、固定値（例えばＡ／Ｆ１２．５相当の値）としてもよい。
【００５７】
続くステップＳ８２では、目標空燃比係数ＫＣＭＤをリッチ化設定値ＫＣＭＤＳＰＲＧに設定し、タイマｔｍＳＰＲＧの値が「０」か否かを判別する（ステップＳ８３）。ｔｍＳＰＲＧ＞０である間は、直ちに本処理を終了し、ｔｍＳＰＲＧ＝０となると、カウンタＮＦＩＳＲＭＤ３の値が「０」か否かを判別する（ステップＳ８４）。最初はＮＦＩＳＲＭＤ３＞０であるので、このカウンタＮＦＩＳＲＭＤ３の値を「１」だけデクリメントして（ステップＳ８５）、本処理を終了する。ＮＦＩＳＲＭＤ３＝０となると、再生モードフラグＦＳＲＣＭＯＤＥを「０」に戻す（ステップＳ８６）とともに、ＳＯｘ吸着量カウンタＣＳＡＢＳの値を「０」に戻して（ステップＳ８７）、本処理を終了する。
【００５８】
以上のように図９の処理により、目標空燃比係数ＫＣＭＤが、昇温モードにおいては「１．０」に設定され、ＮＯｘ浄化装置１５の昇温完了後のＳＯｘ除去モードにおいてはリッチ化設定値ＫＣＭＤＳＰＲＧに設定される。その結果、ＮＯｘ浄化装置１５の昇温を促進して、吸収されたＳＯｘを還元して放出させることができ、硫黄被毒したＮＯｘ浄化装置１５を再生させることができる。
【００５９】
本実施形態では、ＥＣＵ５が、硫黄酸化物量推定手段及び硫黄酸化物除去手段を構成する。より具体的には、図４及び５の処理が硫黄酸化物量推定手段に相当し、図２のステップＳ１３及びＳ１５が硫黄酸化物除去手段に相当する。
なお本発明は上述した実施形態に限るものではなく、種々の変形が可能である。例えば、上述した実施形態では、ＮＯｘ浄化装置１５に吸収されたＳＯｘ量の推定を行い、その推定ＳＯｘ量が設定値に達したときに、ＳＯｘ除去処理を行うようにしたが、特開平１１−２４７６５０号公報に示されるような、ＮＯｘ浄化装置の上流側に配置されたＳＯｘ浄化装置に吸収されたＳＯｘ量の推定に、本発明を適用してもよい。
【００６０】
また上述した実施形態では、目標空燃比係数ＫＣＭＤ並びにエンジン回転数ＮＥ及び吸気管内絶対圧ＰＢＡに応じて、加算値ＳＡＢＳまたは減算値ＳＤＥＳを算出するようにしたが、目標空燃比係数ＫＣＭＤを、ＬＡＦセンサ１４により検出される検出当量比ＫＡＣＴに代えてもよい。
【００６１】
【発明の効果】
以上詳述したように請求項１に記載した発明によれば、機関に供給する混合気の空燃比及び機関の運転状態に応じて、ＮＯｘ浄化手段に吸収された硫黄酸化物量の、単位時間当たりの変化量（すなわち吸収量または放出量）が推定され、該推定された変化量を積算することにより、ＮＯｘ浄化手段に吸収された硫黄酸化物量の推定が行われる。具体的には、空燃比が理論空燃比近傍に設定されているときは、第１推定手段により、ＮＯｘ浄化手段の温度が第１所定温度以上となる機関運転状態においては第１の負の変化量が算出されるとともに、ＮＯｘ浄化手段の温度が第１所定温度より低くなる機関運転状態においては第１の正の変化量が算出される。また空燃比が理論空燃比よりリーン側に設定されているときは、第２推定手段により、機関運転状態に応じて第２の正の変化量が算出される。また空燃比が理論空燃比よりリッチ側に設定されているときは、第３推定手段により、ＮＯｘ浄化手段の温度が第１所定温度より低い第２所定温度以上となる機関運転状態においては第２の負の変化量が算出されるとともに、ＮＯｘ浄化手段の温度が第２所定温度より低くなる機関運転状態においては第３の正の変化量が算出される。ＮＯｘ浄化手段に吸収される硫黄酸化物量は、空燃比が理論空燃比よりリーン側に設定されているとき増加し、空燃比が理論空燃比近傍あるいは理論空燃比よりリッチ側に設定されているときは、機関運転状態により増加または減少するので、第１〜第３推定手段により、硫黄酸化物量の変化量を推定することにより、硫黄酸化物量の単位時間当たりの変化量、すなわち吸収量または放出量を正確に推定することができる。その結果、硫黄酸化物量の推定を従来に比べてより正確に行うことができ、硫黄被毒したＮＯｘ浄化手段の再生処理を最適な時期に実行することができる。
【図面の簡単な説明】
【図１】本発明の一実施形態にかかる内燃機関及びその排気浄化装置の構成を示す図である。
【図２】硫黄被毒したＮＯｘ浄化装置を再生させる制御を行う処理のフローチャートである
【図３】図２のＳＯｘ飽和判断処理を詳細に示すフローチャートである。
【図４】図３のＳＯｘ吸着量推定処理を詳細に示すフローチャートである。
【図５】図３のＳＯｘ吸着量推定処理を詳細に示すフローチャートである。
【図６】図４または５の処理で使用するマップを示す図である。
【図７】図４または５の処理で使用するテーブルを示す図ある。
【図８】目標空燃比係数（ＫＣＭＤ）を算出する処理のフローチャートである。
【図９】ＮＯｘ浄化装置の再生処理実行中における目標空燃比係数の設定を行う処理のフローチャートである。
【符号の説明】
１  内燃機関
５  電子コントロールユニット（硫黄酸化物量推定手段、硫黄酸化物除去手段）
７  吸気管内絶対圧センサ
１０  クランク角度位置センサ
１２  排気管
１５  ＮＯｘ浄化装置（ＮＯｘ浄化手段）
【図１】
0004509435.tif 000002
【図２】
0004509435.tif 000003
【図３】
0004509435.tif 000004
【図４】
0004509435.tif 000005
【図５】
0004509435.tif 000006
【図６】
0004509435.tif 000007
【図７】
0004509435.tif 000008
【図８】
0004509435.tif 000009
【図９】
0004509435.tif 000010