現在の車両の運転状態で、窒素酸化物低減触媒を被毒する硫黄酸化物の質量を計算するステップ、現在の車両の運転状態で、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を計算するステップ、及び前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量から、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を引いた値を積分して、窒素酸化物低減触媒を被毒する硫黄酸化物の量を計算するステップ、を含むことを特徴とする窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量は、排気ガス内の硫黄酸化物の質量、窒素酸化物低減触媒を通過する排気ガスの体積速度、触媒温度、及び現在貯蔵された硫黄酸化物の量に基づいて、計算されることを特徴とする請求項1に記載の窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量を計算するステップは、触媒温度に応じた体積当たりの硫黄酸化物の貯蔵能力を計算するステップ、前記体積当たりの硫黄酸化物の貯蔵能力と、窒素酸化物低減触媒の有効体積を用いて、現在の硫黄酸化物貯蔵能力を計算するステップ、現在の硫黄酸化物貯蔵能力と、現在貯蔵された硫黄酸化物の量を用いて、相対的な硫黄酸化物貯蔵レベルを計算するステップ、相対的な硫黄酸化物貯蔵レベルに応じた基本の硫黄酸化物貯蔵効率を計算するステップ、及び前記基本の硫黄酸化物貯蔵効率と、排気ガス内の硫黄酸化物の質量を用いて、窒素酸化物低減触媒に貯蔵される硫黄酸化物の質量を計算するステップ、を含むことを特徴とする請求項1に記載の窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。前記基本の硫黄酸化物貯蔵効率は、窒素酸化物低減触媒を通過する排気ガスの体積速度と、窒素酸化物低減触媒の温度に応じて、補正されることを特徴とする請求項3に記載の窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量は、窒素酸化物低減触媒の前端の空燃比lambda、実際の還元剤の質量、エンジンの回転数、現在の燃料噴射量、窒素酸化物低減触媒の後端の空燃比、現在貯蔵された硫黄酸化物の量、及び窒素酸化物低減触媒の温度に基づいて、計算されることを特徴とする請求項1に記載の窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量を計算するステップは、エンジンの回転数と現在の燃料噴射量に応じた脱着効率を計算するステップ、前記脱着効率と還元剤の質量を用いて、基本の硫黄酸化物の脱着質量を計算するステップ、及び前記基本の硫黄酸化物の脱着質量を補正して、脱着される硫黄酸化物の質量を計算するステップ、を含むことを特徴とする請求項1に記載の窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。前記基本の硫黄酸化物の脱着質量は、窒素酸化物低減触媒の前端の空燃比に応じた第1補正係数、窒素酸化物低減触媒の後端の空燃比に応じた第2補正係数、及び窒素酸化物低減触媒の温度と現在貯蔵された硫黄酸化物の量に応じた第3補正係数の少なくとも一つを用いて、補正されることを特徴とする請求項6に記載の窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量を積分することによって、窒素酸化物低減触媒から除去される硫黄酸化物の量を計算するステップをさらに含むことを特徴とする請求項1に記載の窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法。燃焼室内に燃料を噴射する第1インジェクタを有するエンジンで発生した排気ガスが流れる排気パイプ、前記排気パイプに取り付けられ、還元剤を噴射する第2インジェクタ、前記第2インジェクタの後端の排気パイプに装着され、前記第2インジェクタで噴射された還元剤を用いて、排気ガスに含まれている窒素酸化物を低減させる窒素酸化物低減触媒、及び、エンジンの運転条件に応じて、前記窒素酸化物低減触媒に貯蔵される硫黄酸化物の量と、窒素酸化物低減触媒から除去される硫黄酸化物の量を予測する制御部、を含み、前記制御部は、排気ガス内の硫黄酸化物の質量、窒素酸化物低減触媒を通過する排気ガスの体積速度、触媒温度、及び現在貯蔵された硫黄酸化物の量に基づいて、窒素酸化物低減触媒を被毒する硫黄酸化物の質量を計算し、窒素酸化物低減触媒の前端の空燃比lambda、実際の還元剤の質量、エンジンの回転数、現在の燃料噴射量、窒素酸化物低減触媒の後端の空燃比、現在貯蔵された硫黄酸化物の量、及び窒素酸化物低減触媒の温度に基づいて、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を計算することを特徴とする排気装置。前記制御部は、前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量から、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を引いた値を積分することによって、窒素酸化物低減触媒を被毒する硫黄酸化物の量を計算することを特徴とする請求項9に記載の排気装置。前記制御部は、前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量を積分することによって、窒素酸化物低減触媒から除去される硫黄酸化物の量を計算することを特徴とする請求項9に記載の排気装置。前記制御部は、触媒温度に応じた体積当たりの硫黄酸化物の貯蔵能力を計算し、前記体積当たりの硫黄酸化物の貯蔵能力と窒素酸化物低減触媒の有効体積を用いて、現在の硫黄酸化物貯蔵能力を計算し、現在の硫黄酸化物貯蔵能力と現在貯蔵された硫黄酸化物の量を用いて、相対的な硫黄酸化物貯蔵レベルを計算し、相対的な硫黄酸化物貯蔵レベルに応じた基本の硫黄酸化物貯蔵効率を計算し、前記基本の硫黄酸化物貯蔵効率と排気ガス内の硫黄酸化物の質量を用いて、窒素酸化物低減触媒に貯蔵される硫黄酸化物の質量を計算することを特徴とする請求項9に記載の排気装置。前記基本の硫黄酸化物貯蔵効率は、窒素酸化物低減触媒を通過する排気ガスの体積速度と、窒素酸化物低減触媒の温度に応じて、補正されることを特徴とする請求項12に記載の排気装置。前記制御部は、エンジンの回転数と現在の燃料噴射量に応じた脱着効率を計算し、前記脱着効率と還元剤の質量を用いて、基本の硫黄酸化物の脱着質量を計算し、前記基本の硫黄酸化物の脱着質量を補正して、脱着される硫黄酸化物の質量を計算することを特徴とする請求項9に記載の排気装置。前記基本の硫黄酸化物の脱着質量は、窒素酸化物低減触媒の前端の空燃比に応じた第1補正係数、窒素酸化物低減触媒の後端の空燃比に応じた第2補正係数、及び窒素酸化物低減触媒の温度と現在貯蔵された硫黄酸化物の量に応じた第3補正係数の少なくとも一つを用いて、補正されることを特徴とする請求項14に記載の排気装置。前記還元剤は、燃料であり、前記排気装置は、第2インジェクタと窒素酸化物低減触媒との間の前記排気パイプ上に装着され、前記燃料を分解する燃料分解触媒をさらに含むことを特徴とする請求項9に記載の排気装置。本発明は、上記問題点を解決するためになされたものであって、本発明の目的は、最適の窒素酸化物低減触媒の温度、最適の窒素酸化物低減触媒を通過する体積速度、最適の窒素酸化物低減触媒の前端の空燃比、及び最適の実際の還元剤質量を定義することによって、窒素酸化物低減触媒に実際に貯蔵される硫黄酸化物の量を正確に予測する方法を提供することにある。本発明の他の目的は、正確に予測された硫黄酸化物の貯蔵量から、窒素酸化物低減触媒の再生時期及び還元剤投入量を正確に予測する排気装置を提供することにある。本発明は、現在の車両の運転状態で、窒素酸化物低減触媒を被毒する硫黄酸化物の質量を計算するステップ、現在の車両の運転状態で、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を計算するステップ、及び前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量から、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を引いた値を積分して、窒素酸化物低減触媒を被毒する硫黄酸化物の量を計算するステップ、を含むことを特徴とする。前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量は、排気ガス内の硫黄酸化物の質量、窒素酸化物低減触媒を通過する排気ガスの体積速度、触媒温度、及び現在貯蔵された硫黄酸化物の量に基づいて、計算されることを特徴とする。前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量を計算するステップは、触媒温度に応じた体積当たりの硫黄酸化物の貯蔵能力を計算するステップ、前記体積当たりの硫黄酸化物の貯蔵能力と、窒素酸化物低減触媒の有効体積を用いて、現在の硫黄酸化物貯蔵能力を計算するステップ、現在の硫黄酸化物貯蔵能力と、現在貯蔵された硫黄酸化物の量を用いて、相対的な硫黄酸化物貯蔵レベルを計算するステップ、相対的な硫黄酸化物貯蔵レベルに応じた基本の硫黄酸化物貯蔵効率を計算するステップ、及び前記基本の硫黄酸化物貯蔵効率と、排気ガス内の硫黄酸化物の質量を用いて、窒素酸化物低減触媒に貯蔵される硫黄酸化物の質量を計算するステップ、を含むことを特徴とする。前記基本の硫黄酸化物貯蔵効率は、窒素酸化物低減触媒を通過する排気ガスの体積速度と、窒素酸化物低減触媒の温度に応じて、補正されることを特徴とする。前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量は、窒素酸化物低減触媒の前端の空燃比lambda、実際の還元剤の質量、エンジンの回転数、現在の燃料噴射量、窒素酸化物低減触媒の後端の空燃比、現在貯蔵された硫黄酸化物の量、及び窒素酸化物低減触媒の温度に基づいて、計算されることを特徴とする。前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量を計算するステップは、エンジンの回転数と現在の燃料噴射量に応じた脱着効率を計算するステップ、前記脱着効率と還元剤の質量を用いて、基本の硫黄酸化物の脱着質量を計算するステップ、及び前記基本の硫黄酸化物の脱着質量を補正して、脱着される硫黄酸化物の質量を計算するステップ、を含むことを特徴とする。前記基本の硫黄酸化物の脱着質量は、窒素酸化物低減触媒の前端の空燃比に応じた第1補正係数、窒素酸化物低減触媒の後端の空燃比に応じた第2補正係数、及び窒素酸化物低減触媒の温度と現在貯蔵された硫黄酸化物の量に応じた第3補正係数の少なくとも一つを用いて、補正されることを特徴とする。前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量を積分することによって、窒素酸化物低減触媒から除去される硫黄酸化物の量を計算するステップをさらに含むことを特徴とする。また、本発明は、燃焼室内に燃料を噴射する第1インジェクタを有するエンジンで発生した排気ガスが流れる排気パイプ、前記排気パイプに取り付けられ、還元剤を噴射する第2インジェクタ、前記第2インジェクタの後端の排気パイプに装着され、前記第2インジェクタで噴射された還元剤を用いて、排気ガスに含まれている窒素酸化物を低減させる窒素酸化物低減触媒、及び、エンジンの運転条件に応じて、前記窒素酸化物低減触媒に貯蔵される硫黄酸化物の量と、窒素酸化物低減触媒から除去される硫黄酸化物の量を予測する制御部、を含み、前記制御部は、排気ガス内の硫黄酸化物の質量、窒素酸化物低減触媒を通過する排気ガスの体積速度、触媒温度、及び現在貯蔵された硫黄酸化物の量に基づいて、窒素酸化物低減触媒を被毒する硫黄酸化物の質量を計算し、窒素酸化物低減触媒の前端の空燃比lambda、実際の還元剤の質量、エンジンの回転数、現在の燃料噴射量、窒素酸化物低減触媒の後端の空燃比、現在貯蔵された硫黄酸化物の量、及び窒素酸化物低減触媒の温度に基づいて、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を計算することを特徴とする。前記制御部は、前記窒素酸化物低減触媒を被毒する硫黄酸化物の質量から、窒素酸化物低減触媒から脱着される硫黄酸化物の質量を引いた値を積分することによって、窒素酸化物低減触媒を被毒する硫黄酸化物の量を計算することを特徴とする。前記制御部は、前記窒素酸化物低減触媒から脱着される硫黄酸化物の質量を積分することによって、窒素酸化物低減触媒から除去される硫黄酸化物の量を計算することを特徴とする。前記制御部は、触媒温度に応じた体積当たりの硫黄酸化物の貯蔵能力を計算し、前記体積当たりの硫黄酸化物の貯蔵能力と窒素酸化物低減触媒の有効体積を用いて、現在の硫黄酸化物貯蔵能力を計算し、現在の硫黄酸化物貯蔵能力と現在貯蔵された硫黄酸化物の量を用いて、相対的な硫黄酸化物貯蔵レベルを計算し、相対的な硫黄酸化物貯蔵レベルに応じた基本の硫黄酸化物貯蔵効率を計算し、前記基本の硫黄酸化物貯蔵効率と排気ガス内の硫黄酸化物の質量を用いて、窒素酸化物低減触媒に貯蔵される硫黄酸化物の質量を計算することを特徴とする。前記基本の硫黄酸化物貯蔵効率は、窒素酸化物低減触媒を通過する排気ガスの体積速度と、窒素酸化物低減触媒の温度に応じて、補正されることを特徴とする。前記制御部は、エンジンの回転数と現在の燃料噴射量に応じた脱着効率を計算し、前記脱着効率と還元剤の質量を用いて、基本の硫黄酸化物の脱着質量を計算し、前記基本の硫黄酸化物の脱着質量を補正して、脱着される硫黄酸化物の質量を計算することを特徴とする。前記基本の硫黄酸化物の脱着質量は、窒素酸化物低減触媒の前端の空燃比に応じた第1補正係数、窒素酸化物低減触媒の後端の空燃比に応じた第2補正係数、及び窒素酸化物低減触媒の温度と現在貯蔵された硫黄酸化物の量に応じた第3補正係数の少なくとも一つを用いて、補正されることを特徴とする。前記還元剤は、燃料であり、前記排気装置は、第2インジェクタと窒素酸化物低減触媒との間の前記排気パイプ上に装着され、前記燃料を分解する燃料分解触媒をさらに含むことを特徴とする。本発明によれば、窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を正確に予測することによって、排気装置の脱硫効率を向上することができる。また、窒素酸化物低減触媒に貯蔵された硫黄酸化物の正確な量に応じて再生時期及び還元剤の噴射量を調節することによって、燃費を向上することができる。以下に、添付図面に基づき、本発明の好ましい実施形態について詳細に説明する。図1は、本発明の実施形態による窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法が適用できる排気装置の一例を示す概略図である。図1に示すように、内燃機関の排気装置は、エンジン10、排気パイプ20、排気ガス再循環ExhaustGasRecirculation;EGR装置80、燃料分解触媒32、媒煙フィルタParticulateFilter30、窒素酸化物低減触媒40、及び制御部50を含む。エンジン10は、燃料及び空気が混合された混合気を燃焼させて、化学的エネルギーを機械的エネルギーに変換する。エンジン10は、吸気マニホールド18に連結されて、燃焼室12の内部に空気が流入され、燃焼過程で発生した排気ガスは、排気マニホールド16に集められた後、エンジン外部に排出される。燃焼室12には第1インジェクタ14が取り付けられていて、燃料を燃焼室12の内部に噴射する。ここでは、ディゼルエンジンを例示したが、希薄燃焼leanburnガソリンエンジンを用いてもよい。ガソリンエンジンを用いる場合、吸気マニホールド18を通じて混合気が燃焼室12の内部に流入され、燃焼室12の上部には、点火するための点火プラグ図示せずが装着される。また、多様な圧縮比、好ましくは16.5以下の圧縮比を有するエンジンを用いてもよい。排気パイプ20は、排気マニホールド16に連結されて、排気ガスを車両の外部に排出させる。排気パイプ20上には、媒煙フィルタ30及び窒素酸化物低減触媒40が装着されていて、排気ガス内に含まれている炭化水素、一酸化炭素、粒子状物質、及び窒素酸化物などが除去される。排気ガス再循環装置80は、排気パイプ20上に装着されて、エンジン10から排出される排気ガス一部を排気ガス再循環装置80を介してエンジンに再供給する。また、排気ガス再循環装置80は、吸気マニホールド18に連結されて、排気ガスの一部を空気に混合して、燃焼温度を制御する。このような燃焼温度の制御は、制御部50の制御によって、吸気マニホールド18に供給される排気ガスの量を調節することによって行われる。排気ガス再循環装置80の後方の排気パイプ20には、第1酸素センサ25が取り付けられていて、排気ガス再循環装置80を通過した排気ガス内の酸素量を検出する。本明細書では、第1酸素センサの測定値をエンジン出口の空燃比lambdaと称する。第2インジェクタ90は、排気ガス再循環装置80の後方の排気パイプ20に取り付けられ、制御部50に電気的に連結されて、制御部50の制御によって排気パイプ20内に燃料の追加噴射を行う。媒煙フィルタ30は、第2インジェクタ90の後方の排気パイプ20に装着されている。媒煙フィルタ30の前端部には、燃料分解触媒FuelCrackingCatalystが備えられている。この場合、燃料分解触媒32は、第2インジェクタ90と窒素酸化物低減触媒40との間に配置される。ここでは、媒煙フィルタ30とは別に燃料分解触媒32が備えられていることを示したが、燃料分解触媒32を媒煙フィルタ30の前端部にコーティングすることもできる。燃料分解触媒32は、触媒反応によって、燃料内に含まれている炭素化合物のチェ-ン環を切って分解させる。つまり、燃料分解触媒32は、熱分解ThermalCracking機能により、炭化水素を構成する連結環を切って分解する。これにより、追加噴射された燃料の有効反応表面積が増加し、高反応性の酸素が含まれている炭化水素OxygenatedHC、CO、H2等を生成する。サーマルクラッキングは、下記のような手順で進められる。C16H34→2n-C8H17→n-C6H13→n-C4H9→C2H5→C2H4C16H34→8C2H4H2ここで、は、ラジカルを意味する。ここで、炭化水素は、排気ガス及び燃料に含まれている炭素と水素とから構成された化合物全てを称する。燃料分解触媒32の後方には、媒煙フィルタ30の一種である媒煙濾過装置ParticulateFilter30が取り付けられていて、排気パイプ20を通じて排出される排気ガスに含まれている粒子状物質ParticulateMatters:PMを捕集する。ここで、媒煙フィルタ30は、媒煙濾過装置30と同様の意味である。しかし、媒煙濾過装置30の代わりに、他の種類の媒煙フィルタ30例えば、触媒媒煙フィルタcatalyzedparticulatefilter:CPFを用いてもよい。また、媒煙濾過装置30には、酸化触媒OxidationCatalystがコーティングされてもよい。このような酸化触媒は、排気ガスに含まれている炭化水素と一酸化炭素を二酸化炭素に酸化し、排気ガスに含まれている一酸化窒素を二酸化窒素に酸化させる。酸化触媒は、媒煙フィルタ30の所定部門に多くコーティングされていてもよく、媒煙フィルタ30の全領域に均等にコーティングされていてもよい。燃料分解触媒32の前方排気パイプには第1温度センサ35が取り付けられ、燃料分解触媒32の入口温度を測定し、燃料分解触媒32の後方には第2温度センサ36が取り付けられ、燃料分解触媒32の出口温度または媒煙フィルタ30の入口温度を測定する。一方、排気パイプ20には差圧センサ55が取り付けられている。差圧センサ55は、媒煙フィルタ30の前端部と後端部との圧力差を測定し、これに対する信号を制御部50に伝達する。制御部50は、差圧センサ55で測定された圧力差が第1設定圧力以上である場合、媒煙フィルタ30を再生するように制御することができる。この場合、第1インジェクタ14で燃料を後噴射することによって、媒煙フィルタ30の内部に捕集された粒子状物質を燃焼させることができる。これとは異なり、第2インジェクタ90で燃料を追加噴射することによって、媒煙フィルタ30を再生させることもできる。窒素酸化物低減触媒40は、媒煙フィルタ30の後方の排気パイプ20上に装着されて、排気ガスに含まれている窒素酸化物を吸着し、燃料の追加噴射によって吸着した窒素酸化物を脱着して還元反応を進めることにより、排気ガスに含まれている窒素酸化物を浄化する。窒素酸化物低減触媒40の前方及び後方には、第3温度センサ60及び第4温度センサ65がそれぞれ取り付けられていて、窒素酸化物低減触媒40の入口部の温度と出口部の温度とを測定する。ここでは、窒素酸化物低減触媒40が2つの部分に分かれていることを示した。これは、担体にコーティングされた金属比率を異にすることによって特定機能を遂行するためのものである。例えば、エンジン10に近い第1部分には、パラジウムPdの比率を高めることによって耐熱機能を強化し、テールパイプに近い第2部分には、白金Ptの比率を高めることによって炭化水素のスリップを防止することができる。これとは異なり、担体にコーティングされた金属比率が全領域において同様である窒素酸化物低減触媒40を用いても良い。また、窒素酸化物低減触媒40の前方排気パイプ20には第2酸素センサ62が取り付けられており、窒素酸化物低減触媒40の後方の排気パイプ20には第3酸素センサ70が取り付けられる。第2酸素センサ62は、窒素酸化物低減触媒40に流入する排気ガスに含まれている酸素量を測定して、これに対する信号を制御部50に伝達する。これにより、制御部50が排気ガスのリーン/リッチ制御lean/richcontrolを行うことを助ける。また、第3酸素センサ70は、本発明の実施形態による内燃機関の排気装置が排気ガスに含まれている有害物質を正常に除去しているか否かをモニターするためのものである。ここでは、排気パイプ20に第2酸素センサ62を追加装着したことを示した。しかし、排気パイプ20に第2酸素センサ62を取り付けずに、第1酸素センサ25及び第3酸素センサ70の測定値、燃料消耗量、及びエンジンの可動ヒストリーhistoryの少なくとも一つに基づいて、窒素酸化物低減触媒40に流入される排気ガスに含まれている酸素量を推定することもできる。本明細書では、第2酸素センサ62の測定値を触媒の前端の空燃比lambdaと称する。制御部50は、各センサから検出された信号に基づいてエンジンの運転条件を判断し、エンジンの運転条件に基づいて燃料の追加噴射量及び追加噴射時期を制御することによって、窒素酸化物低減触媒40に吸着した窒素酸化物を脱着する。一例として、制御部50は、窒素酸化物低減触媒40に吸着した窒素酸化物の量が設定値以上である場合、燃料を追加噴射するように制御する。また、制御部50は、窒素酸化物低減触媒40で窒素酸化物の還元反応が活性化するよう、排気ガス内で窒素酸化物に対する炭化水素の比率が設定比率以上になるように制御する。前記設定された比率は、5である。一方、制御部50は、エンジンの運転条件に基づいて窒素酸化物低減触媒40に吸着した窒素酸化物の量、窒素酸化物低減触媒の後端の窒素酸化物のスリップ量、及び窒素酸化物に対する炭化水素の比率を計算する。このような計算は数多い実験によって定められたマップに基づいて行われる。また、制御部50は、エンジンの運転条件、エンジンの状態または窒素酸化物低減触媒の状態に応じて第2インジェクタ90の燃料噴射パターンを変化させる。ここで、エンジンの状態は、エンジンの作動期間を考慮して推定し、窒素酸化物低減触媒の状態は、窒素酸化物低減触媒の劣化を考慮して推定される。さらに、制御部50は、窒素酸化物低減触媒40に被毒された硫黄酸化物の量が設定量以上になると、第1インジェクタ14と第2インジェクタ90を制御して、燃料を追加噴射することができる。この場合、追加噴射された燃料によって、窒素酸化物低減触媒40に被毒された硫黄酸化物が除去される。燃料の追加噴射によって、窒素酸化物低減触媒40に被毒された硫黄酸化物が除去されることを、窒素酸化物低減触媒40の脱硫または再生と言う。したがって、本明細書において、窒素酸化物低減触媒40の再生とは、窒素酸化物低減触媒40に吸着した窒素酸化物または硫黄酸化物の脱着を意味する。一方、制御部50は、媒煙フィルタ30の再生を進行する。また、制御部50は、第2インジェクタ90での追加噴射を制御する代わりに、第1インジェクタ14での後噴射を制御することによって、窒素酸化物低減触媒40で窒素酸化物の還元反応を活性化することができる。この場合、後噴射された燃料は、燃料分解触媒32で高反応性の還元剤に変化し、窒素酸化物低減触媒40で窒素酸化物の還元反応を促進させる。したがって、本明細書及び特許請求の範囲における追加噴射は、後噴射を含むと解釈すべきである。ここでは、窒素酸化物低減触媒40としてLNT触媒を用いた例を示した。以下、窒素酸化物低減触媒40の一例について、詳しく説明する。窒素酸化物低減触媒40は、担体にコーティングされた第1、2触媒層を含む。第1触媒層は排気ガスに近接して配置され、第2触媒層は担体に近接して配置される。第1触媒層は。排気ガスに含まれている窒素酸化物を酸化させ、酸化された窒素酸化物の一部を、燃えない燃料または排気ガスに含まれている炭化水素との酸化-還元反応によって、還元させる。また、酸化された窒素酸化物の他の一部は第2触媒層に拡散される。第2触媒層は、第1触媒層から拡散した窒素酸化物を吸着し、追加噴射される燃料によって、吸着した窒素酸化物を脱着して、第1触媒層で還元されるようにする。第2触媒層に拡散された窒素酸化物は、硝酸塩の形態で第2触媒層に吸着される。また、第2触媒層で脱着された窒素酸化物は、二酸化窒素の形態で第1触媒層へ移動し、二酸化窒素の一部は還元され、他の一部はスリップし、残りは第1触媒層に吸着される。第2触媒層は吸着物質を含む。このような吸着物質には弱塩基性酸化物が用いられる。弱塩基性酸化物には、アルカリまたはアルカリ土類金属を含む酸化物が用いられ、より具体的に、バリウムを含む酸化物が用いられる。以下、窒素酸化物低減触媒40の作動原理について、詳しく説明する。第2インジェクタ90で燃料が追加噴射されない場合、排気ガスに含まれている窒素酸化物は第1触媒層で酸化される。酸化された窒素酸化物の一部は、排気ガスに含まれている炭化水素と酸化-還元反応を行い、窒素気体に還元される。この過程で、排気ガスに含まれている炭化水素は二酸化炭素に酸化される。また、酸化された窒素酸化物の他の一部と、排気ガスに含まれている窒素酸化物は、第2触媒層に拡散して吸着される。第2インジェクタ90で燃料が追加噴射される場合、追加噴射された燃料が燃料分解触媒を通過し、この過程で、燃料が低分子の炭化水素に変換される。また、低分子の炭化水素の一部は酸素と結合された炭化水素に変換されて、窒素酸化物低減触媒40を通過する。このとき、第2触媒層では、窒素酸化物が前記炭化水素との置換反応によって脱着される。また、第1触媒層では、前記脱着された窒素酸化物と、炭化水素/酸素と結合した炭化水素との間の酸化-還元反応によって、窒素酸化物は窒素気体に還元され、炭化水素/酸素と結合した炭化水素は二酸化炭素に酸化される。よって、排気ガスに含まれている窒素酸化物と炭化水素が浄化される。図2は、本発明の実施形態による窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法に用いられる制御部における入出力関係を示すブロック図である。図2に示すように、SOx質量検出部100、排気ガス質量検出部110、第3温度センサ60、第4温度センサ65、第1酸素センサ25、吸気質量検出部120、還元剤噴射量検出部130、燃料噴射量検出部140、エンジン回転数センサ150、及び第3酸素センサ70は、制御部50に電気的に連結されており、検出した値を制御部50に伝達する。SOx質量検出部100は、窒素酸化物低減触媒40の前端を流れる排気ガスに含まれているSOxの質量を検出する。通常、制御部50は、混合気の燃焼状態、排気ガスの温度、エンジン出口の空燃比、窒素酸化物低減触媒の前端の空燃比、窒素酸化物低減触媒の入口部と出口部との温度差などを用いて、排気ガスに含まれているSOxの質量を予測する。排気ガス質量検出部110は、排気パイプ20を流れる排気ガスの質量を検出する。第3温度センサ60は、窒素酸化物低減触媒40の入口部の温度を検出する。第4温度センサ65は、窒素酸化物低減触媒40の出口部の温度を検出する。第3温度センサ60及び第4温度センサ65で検出された温度は、設定された計算により、窒素酸化物低減触媒40の温度を決定するのに用いられる。このとき、温度に応じた硫黄酸化物の貯蔵特性を用いてもよい。これとは異なって、窒素酸化物低減触媒40の入口部の温度または窒素酸化物低減触媒40の出口部の温度を窒素酸化物低減触媒40の温度で決定することもできる。第1酸素センサ25は、エンジン10の後端の空燃比を検出する。吸気質量検出部120は吸気パイプに装着され、1サイクルにおいて吸入される吸入空気の質量を検出する。還元剤噴射量検出部130は、現在噴射している還元剤の噴射量を検出する。還元剤の噴射量は、制御部50によってデューティ制御されるため、現在のデューティ値を読み取ることによって検出することができる。また、制御部50は、還元剤の噴射量から還元剤の質量を推定することができる。燃料噴射量検出部140は、現在噴射している燃料噴射量を検出する。最近は、燃料は主噴射とパイロット噴射によって噴射される。したがって、燃料噴射量検出部140は一周期において噴射される主噴射量とパイロット噴射量を検出する。また、燃料噴射量は、制御部50によってデューティ制御されるため、現在のデューティ値を読み取ることによって検出することができる。エンジン回転数センサ150は、クランクシャフト図示せずの位相変化からエンジンの回転数を検出する。第3酸素センサ70は、窒素酸化物低減触媒40の後端の空燃比を検出する。制御部50では、前記伝えられた値に基づき、エンジンの運転条件、燃料の噴射量、燃料の噴射時期、燃料の噴射パターン、燃料の追加噴射量すなわち、還元剤の噴射量、追加噴射時期すなわち、再生時期、及び追加噴射パターンを決定し、第1、2インジェクタ14、90を制御するための信号を第1、2インジェクタ14、90に出力する。また、制御部50は、差圧センサ55で測定された値に基づいて媒煙フィルタ30の再生を制御する。上述したように、媒煙フィルタ30の再生は、第1インジェクタ14による後噴射または第2インジェクタ90による追加噴射によって行われる。さらに、制御部50は、窒素酸化物低減触媒40に貯蔵される硫黄酸化物の質量と、窒素酸化物低減触媒40から脱着される硫黄酸化物の質量を計算し、これらに基づいて、窒素酸化物低減触媒40に貯蔵される硫黄酸化物の量と、脱硫によって除去される硫黄酸化物の量を計算する。一方、本発明の実施形態による内燃機関の排気装置には、図2に記載したセンサ以外に複数のセンサが取り付けられても良いが、説明の便宜のため省略する。図3は、本発明の実施形態による窒素酸化物低減触媒に貯蔵される硫黄酸化物の量を予測する方法を行うフローチャートであり、図4は、本発明の実施形態において、窒素酸化物低減触媒に被毒される硫黄酸化物の質量の計算を行うフローチャートであり、図5は、本発明の実施形態において、窒素酸化物低減触媒から脱着される硫黄酸化物の質量の計算を行うフローチャートである。図3に示すように、SOx質量検出部100は、窒素酸化物低減触媒40の前端で排気ガスに含まれているSOxの質量を検出しS200、排気ガス質量検出部110は排気ガスの質量を検出しS210、第3温度センサ60は窒素酸化物低減触媒40の入口部の温度を検出しS220、第4温度センサ65は窒素酸化物低減触媒40の出口部の温度を検出しS230、第1酸素センサ25はエンジン出口の空燃比を検出しS240、吸気質量検出部120は吸気質量を検出するS250。また、還元剤噴射量検出部130は排気ガスに噴射される還元剤の噴射量を検出し、制御部50は前記還元剤の噴射量に基づいて還元剤の質量を推定するS260。燃料噴射量検出部140は検出された燃料噴射量を積算することによって総燃料噴射量を検出しS270、エンジン回転数センサ150はエンジンの回転数を検出しS280、燃料噴射量検出部140は現在噴射される燃料量を検出しS290、第3酸素センサ70は窒素酸化物低減触媒40の後端の空燃比を検出するS300。制御部50は、排気ガスの質量に基づいて窒素酸化物低減触媒40を通過する排気ガスの体積速度を計算するS310。排気ガスの体積速度は、本発明の実施形態による方法を行うのに大変重要な因子であるが、どのような状態を基準として計算するかが問題となる。数多い実験を通じて、正常状態を基準として排気ガスの体積速度を計算すると、適切な結果が得られることが分かった。正常状態における排気ガスの体積速度SVは、次の[式1]によって計算される。制御部50は、窒素酸化物低減触媒40の入口部の温度及び出口部の温度に基づいて窒素酸化物低減触媒40の温度を計算するS320。また、制御部50は、窒素酸化物低減触媒40の前端の空燃比を計算しS330、実際の還元剤質量を計算するS340。窒素酸化物低減触媒40で除去される硫黄酸化物の質量を計算するために、窒素酸化物低減触媒40の前端の空燃比は大変重要な因子である。窒素酸化物低減触媒40の前端の空燃比は、第2酸素センサ62によって検出することができる。しかし、第2酸素センサ62の不正確性及び時間遅延の理由により、窒素酸化物低減触媒40の前端の空燃比は、次の数2により計算できる。また、燃料質量は数3で表される。ここで、qFuel、internalは、燃料噴射量であり、iは、シリンダの回転当りの燃焼数であり、zは、シリンダの個数であり、nは、エンジン速度である。したがって、数2は、次の数4で表される。また、窒素酸化物低減触媒40の前端の空燃比は、次の数5により計算される。制御部50は、窒素酸化物低減触媒40に貯蔵された硫黄酸化物の量を計算しS350、これに基づいて、窒素酸化物低減触媒40を被毒する硫黄酸化物の質量を計算するS360。窒素酸化物低減触媒40に貯蔵された硫黄酸化物の量は、以前再生後に窒素酸化物低減触媒40に残る硫黄酸化物の量と、以前再生後に窒素酸化物低減触媒40に新しく貯蔵される硫黄酸化物の量から計算される。窒素酸化物低減触媒40に貯蔵される硫黄酸化物の質量を計算する過程は、図4を参照しながら説明する。図4に示すように、制御部50は、窒素酸化物低減触媒40の温度に応じた体積当たりの硫黄酸化物の貯蔵能力を計算するS500。窒素酸化物低減触媒40の温度に応じた体積当たりの硫黄酸化物の貯蔵能力は、設定されたマップに保存されている。制御部50は、体積当たりの硫黄酸化物の貯蔵能力と有効体積を用いて、現在の硫黄酸化物貯蔵能力を計算するS510。その後、制御部50は、現在の硫黄酸化物貯蔵能力と現在貯蔵された硫黄酸化物の量を用いて、相対的な硫黄酸化物貯蔵レベルを計算するS520。相対的な硫黄酸化物の貯蔵レベルは、現在貯蔵された硫黄酸化物の量を現在の硫黄酸化物貯蔵能力で割ることによって、計算される。制御部50は、相対的な硫黄酸化物貯蔵レベルに応じた基本の硫黄酸化物貯蔵効率を計算するS530。相対的な硫黄酸化物貯蔵レベルに応じた基本の硫黄酸化物貯蔵効率は、設定されたマップに保存されている。その後、制御部50は、排気ガスの体積速度と窒素酸化物低減触媒40の温度に応じた補正の硫黄酸化物の貯蔵効率を計算するS540。体積速度と触媒温度に応じた補正係数が設定されたマップに保存されており、制御部50は、補正係数を基本の硫黄酸化物貯蔵効率にかけることによって、補正の硫黄酸化物の貯蔵効率を計算する。次に、制御部50は、補正の硫黄酸化物の貯蔵効率と排気ガス内の硫黄酸化物の質量を用いて、窒素酸化物低減触媒40に貯蔵される硫黄酸化物の質量を計算するS550。つまり、補正の硫黄酸化物の貯蔵効率を排気ガス内の硫黄酸化物の質量にかけることによって、窒素酸化物低減触媒40に貯蔵される硫黄酸化物の質量を計算する。それから、制御部50は、窒素酸化物低減触媒40から脱着される硫黄酸化物の質量を計算するS370。窒素酸化物低減触媒40から脱着される硫黄酸化物の質量を計算するステップは、図5を参照しながら説明する。図5に示すように、制御部50は、エンジンの回転数と現在の燃料噴射量を用いて、脱着効率を計算するS600。エンジンの回転数と現在の燃料噴射量に応じた脱着効率は、設定されたマップに保存されている。制御部50は脱着効率と還元剤の質量を用いて、基本の硫黄酸化物の脱着質量を計算するS610。例えば、基本の硫黄酸化物の脱着質量は、脱着効率と還元剤の質量をかけることによって計算できる。その後、制御部50は、窒素酸化物低減触媒の前端の空燃比に応じた第1補正係数を計算しS620、窒素酸化物低減触媒の後端の空燃比に応じた第2補正係数を計算しS630、窒素酸化物低減触媒の温度と現在の硫黄酸化物被毒量に応じた第3補正係数を計算するS640。その後、制御部50は、基本の硫黄酸化物の脱着質量と第1、2、3補正係数を用いて、脱着される硫黄酸化物の質量を計算するS370。例えば、脱着される硫黄酸化物の質量は、基本の硫黄酸化物の脱着質量に第1、2、3補正係数をかけることによって計算できる。制御部50は、窒素酸化物低減触媒40に被毒される硫黄酸化物の質量から、窒素酸化物低減触媒40から脱着される硫黄酸化物の質量を引いた値を積分しS380、窒素酸化物低減触媒40を被毒する硫黄酸化物の量を計算するS390。また、制御部50は、窒素酸化物低減触媒40から脱着される硫黄酸化物の質量を積分しS400、脱硫によって窒素酸化物低減触媒40で除去される硫黄酸化物の質量を計算するS410。以上、本発明に関する好ましい実施形態を説明したが、本発明は前記実施形態に限定されず、本発明の属する技術範囲を逸脱しない範囲での全ての変更が含まれる。