JP2019070376A - エンジン制御装置 - Google Patents

Abstract

【課題】触媒の昇温制御において、触媒での還元反応を考慮することにより、脱離したＳＯｘに起因した白煙の発生の抑制しつつ昇温制御に要する時間の長期化を抑制することのできるエンジン制御装置を提供する。【解決手段】ＥＣＵ３０は、要素Ｍ１〜Ｍ１０によって制御周期ごとにＰＭ再生制御中のＤＯＣ２２ａの床温の目標温度を算出する。これらの要素のうち、推定部Ｍ７は、流入ＳＯｘ量と代表温度とを用いて、すり抜けＳＯ３量を制御周期ごとに推定する。推定部Ｍ８は、ＤＯＣ２２ａにおいてＳＯ３からＳＯ２へと還元される割合であるＳＯ２還元率を推定する。そして、算出部Ｍ９は、サルフェート白煙に関する制約に相当する制約ＳＯ３量と、すり抜けＳＯ３量と、ＳＯ２還元率とを用いて、ＤＯＣ２２ａから脱離することが許されるＳＯ３量である許容脱離ＳＯ３量を制御周期ごとに算出する。【選択図】図３

Description

本発明は、エンジン制御装置に関し、詳しくは、排気管に触媒を備えるディーゼルエンジンの制御装置に関する。

従来、ディーゼルエンジンの排気管に設けられる触媒に吸着した硫黄酸化物（ＳＯ_２またはＳＯ_３をいい、以下においてこれらを区別しない場合には「ＳＯｘ」と総称する。）を、定期的に脱離させる昇温制御を行うことが知られている。昇温制御に関連する文献として、例えば特開２０１７−１０６３８１号公報が挙げられる。この公報には、触媒におけるＳＯｘの反応をモデル化することにより、触媒からＳＯ３の状態で排出されるＳＯｘ量を推定する技術が開示されている。ＳＯ_３はＨ_２０に溶解することでミスト化し、サルフェート白煙として視認される。上記の技術では、モデル演算によって推定された触媒から排出されるＳＯ_３量がサルフェート白煙に関する制約に抵触しないように触媒の目標床温が決定される。

特開２０１７−１０６３８１号公報 特開２０１５−１６９１０５号公報

上記公報の技術では、具体的には、触媒に流入したＳＯｘが触媒に吸着・脱離する反応と、触媒をすり抜けるＳＯｘのうちＳＯ_２が酸化してＳＯ_３に転化する反応をモデル化している。しかしながら、触媒におけるＳＯｘの反応は、触媒への吸着・脱離やＳＯ_２の酸化反応だけではない。すなわち、触媒の昇温制御において排気に燃料が添加されると、触媒ではＳＯ_３からＳＯ_２への還元反応が起きる。上記公報の技術では、触媒でのＳＯ_３の還元反応をモデル化していない。このため、上記公報の技術では、触媒からＳＯ_３の状態で排出されるＳＯｘ量が実際よりも多く推定されるおそれがある。この場合、触媒の目標床温が実際に許容される温度よりも低く設定されてしまうため、床温制御の期間が長くなり燃費の悪化を招いてしまう。

本発明は、上述のような課題に鑑みてなされたもので、触媒の昇温制御において、触媒での還元反応を考慮することにより、脱離したＳＯｘに起因した白煙の発生の抑制しつつ昇温制御に要する時間の長期化を抑制することのできるエンジン制御装置を提供することを目的とする。

上記の課題を解決するため、本発明は、ディーゼルエンジンの排気管に設けられる浄化装置と、浄化装置に未燃燃料を供給する燃料供給手段と、燃料供給手段から未燃燃料を供給することにより浄化装置の温度を微粒子が燃焼する温度域の目標温度まで上昇させる昇温制御を実行する制御手段と、を有するエンジン制御装置に適用される。エンジン制御装置の制御手段は、所定の制御周期ごとに浄化装置の温度の代表値である代表温度を取得する温度取得手段と、浄化装置に流入するＳＯｘ量を流入ＳＯｘ量として制御周期ごとに推定する流入ＳＯｘ量推定手段と、流入ＳＯｘ量と代表温度とを用いて、浄化装置の温度上昇中の各温度において浄化装置に最終的に吸着するＳＯｘ量を代表温度に関連付けたグラフとして表される最終吸着ＳＯｘ分布を制御周期ごとに推定する最終吸着ＳＯｘ分布推定手段と、流入ＳＯｘ量と代表温度とを用いて、浄化装置にＳＯｘの状態で流入して浄化装置に吸着することなくすり抜けるＳＯｘのうちＳＯ_３に転化するＳＯｘ量をすり抜けＳＯ_３量として制御周期ごとに推定するすり抜けＳＯ_３量推定手段と、浄化装置においてＳＯ_３からＳＯ_２へと還元される割合であるＳＯ_２還元率を推定するＳＯ_２還元率推定手段と、サルフェート白煙に関する制約に相当する浄化装置の下流における制約ＳＯ_３量と、すり抜けＳＯ_３量と、ＳＯ_２還元率とを用いて、浄化装置から脱離することが許されるＳＯ_３量を許容脱離ＳＯ_３量として制御周期ごとに算出する許容脱離ＳＯ_３量算出手段と、最終吸着ＳＯｘ分布と許容脱離ＳＯ_３量とを用いて、目標温度を制御周期ごとに算出する目標温度算出手段と、を備えている。

また、本発明のエンジン制御装置において、ＳＯ_２還元率推定手段は、ＳＯ_２還元率と、燃料供給手段による未燃燃料の供給量と、浄化装置に流入するガス量との関係に基づいて、ＳＯ_２還元率を推定するように構成されていてもよい。

また、本発明のエンジン制御装置において、すり抜けＳＯ_３量推定手段は、以下のように構成されていてもよい。すなわち、すり抜けＳＯ_３量推定手段は、浄化装置の温度上昇中の各温度において浄化装置に吸着するＳＯｘ量を代表温度に関連付けたグラフとして表される吸着ＳＯｘ分布と、浄化装置の温度上昇中の各温度において浄化装置に吸着するＳＯｘ最大量を代表温度に関連付けたグラフとして表される飽和ＳＯｘ分布と、を用いて、浄化装置におけるＳＯｘ飽和率を制御周期ごとに推定し、流入ＳＯｘ量とＳＯｘ飽和率とを用いて、浄化装置に流入して浄化装置に新たに吸着するＳＯｘ量を新規吸着ＳＯｘ量として制御周期ごとに推定し、新規吸着ＳＯｘ量を用いて、浄化装置に流入して浄化装置に吸着することなくすり抜けるＳＯｘ量をすり抜けＳＯｘ量として制御周期ごとに推定し、そして、浄化装置においてＳＯ_３に転化するＳＯ_２の転化率と代表温度との関係を表した転化率マップと、代表温度と、すり抜けＳＯｘ量と、を用いて、すり抜けＳＯ_３量を制御周期ごとに推定するように構成されている。

また、本発明のエンジン制御装置において、最終吸着ＳＯｘ分布推定手段は、以下のように構成されていてもよい。すなわち、最終吸着ＳＯｘ分布推定手段は、流入ＳＯｘ量と代表温度とを用いて、浄化装置から新たに脱離するＳＯｘ量を新規脱離ＳＯｘ量として制御周期ごとに推定し、そして、最終吸着ＳＯｘ分布推定手段は、新規脱離ＳＯｘ量を用いて、最終吸着ＳＯｘ分布を制御周期ごとに推定するように構成されている。

また、本発明のエンジン制御装置において、温度取得手段は、排気管における浄化装置の下流側に流れたガスの温度を代表温度として取得するように構成されていてもよい。

また、本発明のエンジン制御装置において、浄化装置は、排気管を流れる微粒子を捕集するフィルタを含んで構成されていてもよい。そして、このような構成では、制御手段は、昇温制御を、フィルタに捕集された微粒子量の推定値が除去要求量に到達したときに開始するように構成されていてもよい。

本発明によれば、許容脱離ＳＯ_３量の算出において、サルフェート白煙に関する制約に相当する浄化装置の下流における制約ＳＯ_３量と、すり抜けＳＯ_３量と、ＳＯ_３からＳＯ_２へと還元される割合であるＳＯ_２還元率とが考慮される。ＳＯ_２還元率が考慮されると、ＳＯ_２還元率が考慮されない場合よりも許容脱離ＳＯ_３量の推定精度を高めることができる。これにより、床温制御の目標温度を実際に許容される床温に近づけることができるので、脱離したＳＯｘに起因した白煙の発生の抑制しつつ昇温制御に要する時間の長期化を抑制することが可能となる。

実施の形態のシステム構成を示す図である。 ＤＯＣ２２ａにおけるＳＯｘの吸着と脱離を説明するための図である。 目標床温Ｔｔｒｇを算出するためのロジックを示す機能ブロック図である。 吸着ＳＯｘ分布と飽和ＳＯｘ分布を説明するための図である。 基準飽和ＳＯｘ分布と補正後の飽和ＳＯｘ分布の関係を説明するための図である。 総吸着余裕ＳＯ_２量を説明するための図である。 新規吸着ＳＯｘ量とすり抜けＳＯｘ量の関係を説明するための図である。 吸着率ｍａｐを説明するための図である。 吸着後ＳＯｘ分布を説明するための図である。 脱離可能総ＳＯｘ量を説明するための図である。 最終吸着ＳＯｘ分布と吸着後ＳＯｘ分布の関係を説明するための図である。 ＳＯ_３転化率ｍａｐを説明するための図である。 ＳＯ_３転化率補正ｍａｐを説明するための図である。 ＳＯ_２還元率ｍａｐを説明するための図である。 許容脱離ＳＯ_３量を説明するための図である。 目標床温Ｔｔｒｇを説明するための図である。 実施の形態のシステムが備えるＥＣＵのハードウェア構成の例を示す図である。 実施の形態のシステムが備えるＥＣＵのハードウェア構成の他の例を示す図である。

以下、図面を参照して本発明の実施の形態について説明する。ただし、以下に示す実施の形態において各要素の個数、数量、量、範囲等の数に言及した場合、特に明示した場合や原理的に明らかにその数に特定される場合を除いて、その言及した数に、この発明が限定されるものではない。また、以下に示す実施の形態において説明する構造やステップ等は、特に明示した場合や明らかに原理的にそれに特定される場合を除いて、この発明に必ずしも必須のものではない。

実施の形態１．
本発明の実施の形態１について図を参照して説明する。

［システム構成の説明］
図１は、実施の形態のシステム構成を示す図である。図１に示すシステムは、車両に搭載されるディーゼルエンジン１０（以下単に「エンジン１０」ともいう。）を備えたエンジン制御装置として構成されている。エンジン１０の各気筒には、燃料としての軽油を噴射するインジェクタ１２が設けられている。なお、図１に描かれるエンジン１０は直列４気筒エンジンであるが、エンジン１０の気筒数および気筒配列は特に限定されない。また、図１には、４つのインジェクタ１２のうちの１つが描かれている。

エンジン１０の排気マニホールド１４には、ターボチャージャ１６の排気タービン１６ａの入口が接続されている。排気タービン１６ａは吸気管１８に設けられたコンプレッサ１６ｂに連結されている。コンプレッサ１６ｂは、排気タービン１６ａの回転により駆動して吸気を過給する。排気タービン１６ａの出口には排気管２０が接続されている。排気管２０には排気浄化装置２２が設けられている。排気浄化装置２２は、ＤＯＣ（Diesel Oxidation Catalyst）２２ａと、ＤＰＦ（Diesel Particulate Filter）２２ｂと、を備えている。ＤＯＣ２２ａは、排気中の炭化水素（ＨＣ）や一酸化炭素（ＣＯ）を酸化して、水（Ｈ_２Ｏ）や二酸化炭素（ＣＯ_２）に転化させる機能を有する触媒である。ＤＰＦ２２ｂは、排気中に含まれる微粒子（ＰＭ）を捕集するフィルタである。排気浄化装置２２の上流には、インジェクタ１２と共通の未燃燃料を排気管２０に添加する燃料添加弁２４が設けられている。燃料添加弁２４は、燃料供給手段の例示である。

図１に示すシステムは、制御手段としてのＥＣＵ（Electronic Control Unit）３０を備えている。ＥＣＵ３０は、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、ＣＰＵ（Central Processing Unit）及び入出力インターフェースを備えるマイクロコンピュータである。ＥＣＵ３０は、入出力インターフェースを通して車両に搭載された各種センサの信号を取り込み処理する。各種センサには、吸気管１８の入口付近に設けられたエアフローメータ３２、ＤＯＣ２２ａの出口温度を検出する温度センサ３４、ＤＰＦ２２ｂの上下流における圧力差を検出する差圧センサ３６が含まれている。ＥＣＵ３０は、取り込んだ各センサの信号を処理して所定の制御プログラムに従って各種アクチュエータを操作する。ＥＣＵ３０によって操作されるアクチュエータには、上述したインジェクタ１２、燃料添加弁２４が含まれている。

［ＤＰＦ２２ｂの再生制御］
本実施の形態では、ＥＣＵ３０によるエンジン制御として、ＤＰＦ２２ｂの昇温制御（以下「ＰＭ再生制御」ともいう。）が行われる。ＰＭ再生制御は、ＤＰＦ２２ｂで捕集したＰＭの推定値が除去要求量に到達したときに、燃料添加弁２４から燃料を添加する制御である。例えば差圧センサ３６で検出した圧力差が所定値に到達したときに、ＰＭの推定値が除去要求量に到達したと判断することができる。燃料添加弁２４から燃料を添加することでＤＯＣ２２ａにおいて添加燃料を酸化し、この酸化反応熱によりＤＰＦ２２ｂの床温を６００℃以上まで上昇させる。これにより、ＤＰＦ２２ｂで捕集したＰＭを燃焼除去できるので、ＤＰＦ２２ｂの捕集機能を回復させることができる。なお、ＤＰＦ２２ｂの床温を６００℃以上に上昇させるための燃料添加弁２４からの添加燃料量（以下「ＤＰＦ用燃料量）は、ＤＰＦの床温と関連付けたマップに基づき決定されるものとする。このようなマップは例えばＥＣＵ３０のＲＯＭに記憶させておき、ＤＰＦ２２ｂの実際の床温に応じて適宜読み出すことができる。

［ＰＭ再生制御における問題点］
ところで、ディーゼルエンジンの燃料や潤滑油には一般に硫黄が含まれており、燃料の燃焼に伴ってこうした硫黄からＳＯｘが生成する。本実施の形態においても同様で、エンジン１０での燃料の燃焼に伴ってＳＯｘが生成する。生成したＳＯｘはエンジン１０から排出されて排気浄化装置２２に流入し、主にＤＯＣ２２ａに吸着する。但し、ＤＯＣ２２ａの床温が高くなると、ここに吸着していたＳＯｘが脱離し始める。ＤＯＣ２２ａの組成等により多少変動するものの、ＰＭ再生制御を行う温度域ではＤＯＣ２２ａからＳＯｘが脱離して下流側に放出される。

ＤＯＣ２２ａにおけるＳＯｘの吸着と脱離について、図２を参照しながら説明する。図２は、ＤＯＣ２２ａにおけるＳＯｘの吸着と脱離を説明するための図である。この図に示すように、ＤＯＣ２２ａは、基材（不図示）の表面を覆うコート材２２ｃと、貴金属２２ｄ（Ｐｔ，Ｐｄ等）とを備えている。貴金属２２ｄはコート材２２ｃに分散担持されており、ＨＣやＣＯを酸化する際の活性点となる。但し、排気中のＳＯ_２が貴金属２２ｄに吸着し、または、排気中のＳＯ_３がコート材２２ｃに吸着する。貴金属２２ｄに吸着したＳＯ_２の一部は貴金属２２ｄから脱離して排気中に戻り、または、貴金属２２ｄ上で酸化されてＳＯ_３となり、ＳＯ_３の状態でコート材２２ｃに吸着する。即ち、貴金属２２ｄにはＳＯ_２が吸着し、コート材２２ｃには排気由来のＳＯ_３とＳＯ_２由来のＳＯ_３が吸着する。何れにせよ、ＳＯｘが吸着することで、ＤＯＣ２２ａにおけるＨＣ等の酸化機能が阻害されることになる。

上述した２経路によりコート材２２ｃに吸着したＳＯ_３は、コート材２２ｃの床温が高くなることで脱離する。また、コート材２２ｃの床温が高くなることで貴金属２２ｄ上でのＳＯ_２からＳＯ_３への転化が促進されるので、このようなＳＯ_３もコート材２２ｃから脱離する。従って、ＰＭ再生制御を行うことで、上述したＤＰＦ２２ｂの捕集機能だけでなく、ＤＯＣ２２ａにおけるＨＣ等の酸化機能を回復させることもできる。ところが、図２に示すように、コート材２２ｃから脱離したＳＯ_３が排気管２０に存在するＨ_２Ｏと反応することでＨ_２ＳＯ_４が発生する。そして、このＨ_２ＳＯ_４の濃度が一定濃度を超えると視認可能な白煙（サルフェート白煙）となるので、エンジン１０を搭載した車両の商品価値を損ねてしまうおそれがある。

［ＤＯＣ２２ａの温度制御］
ＤＯＣ２２ａの下流における排気中のＨ_２ＳＯ_４の濃度が高くなり過ぎないように燃料添加弁２４から燃料を添加すれば、ＰＭ再生制御中のサルフェート白煙の発生を抑制できる。そこで、本実施の形態のＥＣＵ３０は、制御モデルを用いてＤＯＣ２２ａの下流におけるＳＯ_３の濃度を算出する。そして、ＥＣＵ３０は、算出したＳＯ_３の濃度がサルフェート白煙に関する制約を満たすように、ＰＭ再生制御中のＤＯＣ２２ａの床温の目標温度（以下「目標床温Ｔｔｒｇ」ともいう。）を算出する。なお、このような制約ＳＯ_３濃度（ＤＯＣ２２ａの下流におけるＳＯ_３濃度の上限値）は、例えばＥＣＵ３０のＲＯＭに記憶させておくことができる。

ＥＣＵ３０は、目標床温Ｔｔｒｇに基づいて、燃料添加弁２４から添加する燃料量を算出する。以下の説明では、目標床温Ｔｔｒｇを実現するための燃料量を、制約充足用燃料量と称する。制約充足用燃料量よりもＤＰＦ用燃料量の方が多い場合に、ＤＰＦ用燃料量ではなく制約充足用燃料量を採用する。これにより、サルフェート白煙に関する制約を満たしつつ、ＤＯＣ２２ａにおけるＨＣ等の酸化機能を回復させることができる。

［本実施の形態の特徴］
本願の発明者は、ＰＭ再生制御中のＤＯＣ２２ａでのＳＯｘ反応について鋭意研究を重ねた。そして、ＰＭ再生制御中のＤＯＣ２２ａでは、ＳＯ_２からＳＯ_３への酸化反応だけでなく、ＳＯ_３からＳＯ_２への還元反応も活発に起きていることに着眼した。従来の技術のように、還元反応を考慮しない制御モデルでは、ＤＯＣ２２ａから排出されるＳＯ_３量が実際よりも多く算出されてしまう。その結果、目標床温Ｔｔｒｇは、実際に許容される床温よりも低い温度に設定される。この場合、ＰＭ再生制御の実行期間が長期化し、燃費の悪化を招いてしまう。

本実施の形態のシステムは、ＰＭ再生制御中のＤＯＣ２２ａでのＳＯ_３の還元反応を制御モデルに盛り込むこととした。これにより、ＤＯＣ２２ａから排出されるＳＯ_３量の推定精度を高めることができるので、目標床温Ｔｔｒｇを実際に許容される床温に近づけることができる。その結果、ＰＭ再生制御の実行期間が長期化を抑制することができるので、燃費の悪化が抑制される。以下、本実施の形態のシステムにおいて実行される目標床温Ｔｔｒｇの算出ロジックについて詳細に説明する。

［目標床温Ｔｔｒｇの算出ロジック］
図３は、目標床温Ｔｔｒｇを算出するためのロジックを示す機能ブロック図であり、これはＥＣＵ３０により実現される。この図に示すように、ＥＣＵ３０は、流入ＳＯｘ量推定部Ｍ１と、ＳＯｘ飽和率推定部Ｍ２と、新規吸着ＳＯｘ量およびすり抜けＳＯｘ量推定部Ｍ３と、吸着後ＳＯｘ分布推定部Ｍ４と、新規脱離ＳＯｘ量推定部Ｍ５と、最終吸着ＳＯｘ分布推定部Ｍ６と、すり抜けＳＯ_３量推定部Ｍ７と、ＳＯ_２還元率推定部Ｍ８と、許容脱離ＳＯ_３量算出部Ｍ９と、白煙抑制目標床温算出部Ｍ１０と、を備えており、これらの要素Ｍ１〜Ｍ９によって制御周期ごとに（具体的にはエンジン１０の燃焼サイクルごとに）目標床温Ｔｔｒｇを算出するものとする。なお、以下の説明においては、要素Ｍ１〜Ｍ９を簡略化するものとし、例えば流入ＳＯｘ量推定部Ｍ１を「推定部Ｍ１」ともいう。

推定部Ｍ１は、ＤＯＣ２２ａに流入するＳＯｘの量（以下「流入ＳＯｘ量」ともいう。）を推定する。なお、本明細書でいう「ＤＯＣ２２ａに流入するＳＯｘ」には、エンジン１０で生成し、ここから排出されてＤＯＣ２２ａに流入するＳＯｘだけでなく、燃料添加弁２４から添加された燃料のＤＯＣ２２ａでの酸化反応に伴って生成し、ＤＯＣ２２ａ上を流れるＳＯｘも含まれるものとする。

推定部Ｍ１は、具体的に、インジェクタ１２からの噴射燃料量（筒内噴射量）および燃料添加弁２４からの添加燃料量（排気添加量）を変数とする次式（１）により、第ｔ番目のサイクルにおける流入ＳＯｘ量を推定する。なお、式（１）の燃料Ｓ濃度は、燃料中の硫黄濃度であり、燃料供給系に別途設けた硫黄濃度センサの検出値を用いてもよく、設定値を用いてもよい。

式（１）の流入燃料量（排気添加量（ｔ），筒内噴射量（ｔ））は、「ＤＯＣ２２ａに流入するＳＯｘ」の由来となった燃料の第ｔ番目のサイクルにおける量であり、燃料の比重（軽油比重）を用いて次式（２）により算出される。

なお、以下の説明においては、流入ＳＯｘ量（排気添加量（ｔ），筒内噴射量（ｔ））を流入ＳＯｘ量（ｔ）ともいう。また、流入燃料量（排気添加量（ｔ），筒内噴射量（ｔ））を流入燃料量（ｔ）ともいう。

推定部Ｍ２は、ＤＯＣ２２ａにおけるＳＯｘの飽和率（以下「ＳＯｘ飽和率」ともいう。）を推定する。ＳＯｘ飽和率の推定には、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯｘの量（以下「吸着ＳＯｘ量」ともいう。）を、ＤＯＣ２２ａの床温に関連付けたグラフとして表される分布（以下「吸着ＳＯｘ分布」ともいう。）と、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯｘの最大量（以下「飽和ＳＯｘ量」ともいう。）を、ＤＯＣ２２ａの床温に関連付けたグラフとして表される分布（以下「飽和ＳＯｘ分布」ともいう。）と、が用いられる。先ず、吸着ＳＯｘ分布と飽和ＳＯｘ分布について、ＳＯ_３を例として図４を参照しながら説明する。

図４は、吸着ＳＯｘ分布と飽和ＳＯｘ分布を説明するための図である。図４に「吸着ＳＯ_３量」として示したデータは、次の手法により収集したものである。具体的には先ず、図４に「現在温度」として示した床温において、十分な量のＳＯｘをＤＯＣ２２ａに吸着させる。続いて、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａから脱離したＳＯ_３の量を、上昇速度を一定とした条件のもとで測定する。そして、この脱離ＳＯ_３量をＤＯＣ２２ａの床温に関連付けてグラフを作成する。これにより、脱離ＳＯ_３量を表した分布（以下「脱離ＳＯ_３分布」ともいう。）を得ることができる。これと同様の手法により、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａから脱離したＳＯ_２の量を、ＤＯＣ２２ａの床温に関連付けたグラフ（以下「脱離ＳＯ_２分布」ともいう。）を得ることもできる。なお、ＤＯＣ２２ａから脱離するＳＯ_３については、これをセンサで直接的に測定してもよいし、ＳＯｘまたはＳＯ_２を検出するセンサを用いて両者を測定し、これらの差から算出してもよい（ＳＯ_３＝ＳＯｘ−ＳＯ_２）。

ここで、ＤＯＣ２２ａの床温上昇中にＤＯＣ２２ａから脱離するＳＯ_３は、実際には図４に「現在温度」として示した床温でＤＯＣ２２ａに吸着させたＳＯ_３である。しかし、ある床温でＤＯＣ２２ａから脱離するＳＯ_３は、その床温に至るまでＤＯＣ２２ａに吸着し続けることのできたＳＯ_３であり、更に言えば、その床温でＤＯＣ２２ａに吸着することができたＳＯ_３であると考えることもできる。このような仮定に基づいて、上述した脱離ＳＯ_３分布の縦軸を、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯ_３の量に置き換えると、図４に示した「吸着ＳＯ_３量」のデータのグラフ、即ち、吸着ＳＯ_３分布を得ることができる。そして、これと同様の手法によって、吸着ＳＯ_２分布を得ることもできる。

また、図４に「飽和ＳＯ_３量」として示したデータは、「吸着ＳＯ_３量」のデータと同様の手法により収集したものである。この「飽和ＳＯ_３量」のデータは、具体的に、上昇速度を極低速とした条件のもと、ＤＯＣ２２ａの床温上昇中の各床温（例えば５℃間隔）においてＤＯＣ２２ａから脱離したＳＯ_３の量に相当している。ＤＯＣ２２ａの床温の上昇速度が極低速であることから、この「飽和ＳＯ_３量」のデータは、ＤＯＣ２２ａから脱離したＳＯ_３の量の最大値であると考えることができる。また、この最大値に対しては、上述した仮定を適用することができる。即ち、ある床温でＤＯＣ２２ａから脱離するＳＯ_３の最大量は、その床温でＤＯＣ２２ａに吸着することができたＳＯ_３の最大量に等しいと考えることができる。このような仮定に基づいて上述した脱離ＳＯ_３分布の縦軸を上記ＳＯ_３の最大量に置き換えると、図４に示した「飽和ＳＯ_３量」のデータのグラフ、即ち、飽和ＳＯ_３分布を得ることができる。そして、これと同様の手法によって、飽和ＳＯ_２分布を得ることもできる。

推定部Ｍ２は、第ｔ番目のサイクルにおけるＤＯＣ２２ａの現在床温Ｔ_２を変数とする次式（３）により、第ｔ番目のサイクルにおけるＳＯｘ飽和率（Ｔ_２（ｔ），ｔ）を推定する。なお、現在床温Ｔ_２は、ＤＯＣ２２ａの現在の床温の代表温度として、例えば温度センサ３４の検出値を使用することができる。

式（３）のＳＯｘ飽和率（Ｔ_２（ｔ），ｔ）の算出過程は次のとおりである。先ず、ＤＯＣ２２ａの床温上昇中における床温Ｔ_１と、現在床温Ｔ_２とを変数とする次式（４）および（５）により、第ｔ番目のサイクルにおける飽和ＳＯ_２分布（Ｔ_１，Ｔ_２（ｔ），ｔ）および飽和ＳＯ_３分布（Ｔ_１，Ｔ_２（ｔ），ｔ）をそれぞれ算出する。

式（４）の基準飽和ＳＯ_２分布は、十分な量のＳＯｘをＤＯＣ２２ａに吸着させるときの床温（図４の「現在温度」）を基準床温（例えば上述した吸着限界量が最大となる３００℃付近の床温）として作成した飽和ＳＯ_２分布である。式（５）の基準飽和ＳＯ_３分布もこれと同様である。式（４）の床温補正ＳＯ_２ｍａｐ（Ｔ_２（ｔ））は、基準飽和ＳＯ_３分布を現在床温Ｔ_２の飽和ＳＯ_２分布に変換するための補正値を定めたマップである。式（５）の床温補正ＳＯ_３ｍａｐ（Ｔ_２（ｔ））もこれと同様である。このような基準飽和ＳＯｘ分布と補正マップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、現在床温Ｔ_２に応じて適宜読み出すことができる。

基準飽和ＳＯｘ分布と補正後の飽和ＳＯｘ分布の関係を、ＳＯ_２を例として図５を参照しながら説明する。図５は、基準飽和ＳＯｘ分布と補正後の飽和ＳＯｘ分布の関係を説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、ＤＯＣ２２ａの床温上昇中にＤＯＣ２２ａからＳＯ_２が脱離し始める温度（下限温度）と、ＤＯＣ２２ａからＳＯ_２が脱離し終わる温度（上限温度）に、それぞれ相当している。この図に示す３種類の分布の違いは、現在床温Ｔ_２にある。即ち、現在床温Ｔ_２が基準温度と等しい場合は、補正後の飽和ＳＯ_２分布の形状が基準飽和ＳＯ_２分布の形状と一致する（中央）。一方、現在床温Ｔ_２が基準温度よりも低い場合（左方）や、現在床温Ｔ_２が基準温度よりも高い場合（右方）は、補正後の飽和ＳＯ_２分布の形状が基準飽和ＳＯ_２分布の形状と一致しなくなる。なお、現在床温Ｔ_２が基準温度よりも高い場合（右方）には、補正後の飽和ＳＯ_２分布の形状が、現在床温Ｔ_２よりも低温側のデータが欠落しているような形状となる。この理由は、現在床温Ｔ_２よりも低温側では、本来であればこの床温域においてＤＯＣ２２ａに吸着し続けることのできたはずのＳＯｘが、ＤＯＣ２２ａから既に脱離していると考えられるためである。

続いて、式（４）により算出した飽和ＳＯ_２分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を次式（６）に代入して、第ｔ番目のサイクルにおける総飽和ＳＯ_２量（Ｔ_２（ｔ），ｔ）を算出する。また、式（５）により算出した飽和ＳＯ_３分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を次式（７）に代入して、第ｔ番目のサイクルにおける総飽和ＳＯ_３量を算出する。

総飽和ＳＯ_２量（Ｔ_２（ｔ），ｔ）および総飽和ＳＯ_３量（Ｔ_２（ｔ），ｔ）を算出したら、次式（８）にこれらを代入して、第ｔ番目のサイクルにおける総飽和量（Ｔ_２（ｔ），ｔ）を算出する。

なお、以下の説明においては、総飽和ＳＯ_２量（Ｔ_２（ｔ），ｔ）を単に総飽和ＳＯ_２量（ｔ）ともいう。また、総飽和ＳＯ_３量（Ｔ_２（ｔ），ｔ）を単に総飽和ＳＯ_３量（ｔ）ともいう。また、総飽和量（Ｔ_２（ｔ），ｔ）を単に総飽和量（ｔ）ともいう。

式（８）により総飽和量（ｔ）を算出したら、飽和ＳＯ_２分布（Ｔ_１，Ｔ_２（ｔ），ｔ）と、推定部Ｍ６で推定した第ｔ番目のサイクルにおける最終吸着ＳＯ_２分布（Ｔ_１，ｔ）とを次式（９）に代入して、第ｔ番目のサイクルにおける吸着余裕ＳＯ_２分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を算出する。また、飽和ＳＯ_３分布（Ｔ_１，Ｔ_２（ｔ），ｔ）と、推定部Ｍ６で推定した第ｔ番目のサイクルにおける最終吸着ＳＯ_３分布（Ｔ_１，ｔ）とを次式（１０）に代入して、第ｔ番目のサイクルにおける吸着余裕ＳＯ_３分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を算出する。

続いて、式（９）により算出した吸着余裕ＳＯ_２分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を次式（１１）に代入して、第ｔ番目のサイクルにおける総吸着余裕ＳＯ_２量（Ｔ_２（ｔ），ｔ）を算出する。また、式（１０）により算出した吸着余裕ＳＯ_３分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を次式（１２）に代入して、第ｔ番目のサイクルにおける総吸着余裕ＳＯ_３量（Ｔ_２（ｔ），ｔ）を算出する。

なお、以下の説明においては、吸着余裕ＳＯ_２分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を単に吸着余裕ＳＯ_２分布（ｔ）ともいう。また、吸着余裕ＳＯ_３分布（Ｔ_１，Ｔ_２（ｔ），ｔ）を単に吸着余裕ＳＯ_３分布（ｔ）ともいう。また、総吸着余裕ＳＯ_２量（Ｔ_２（ｔ），ｔ）を単に総吸着余裕ＳＯ_２量（ｔ）ともいう。また、総吸着余裕ＳＯ_３量（Ｔ_２（ｔ），ｔ）を単に総吸着余裕ＳＯ_３量（ｔ）ともいう。

図６を参照して、総吸着余裕ＳＯ_２量（ｔ）を説明する。図６は、総吸着余裕ＳＯ_２量を説明するための図である。なお、総吸着余裕ＳＯ_３量（ｔ）についてはこれと同様である。この図に示すように、総吸着余裕ＳＯ_２量（ｔ）は、飽和ＳＯ_２分布から、飽和ＳＯ_２分布と吸着ＳＯ_２分布の重複部分を除いた面積として表すことができる。なお、この図の右方の分布に領域Ａとして示すように、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯ_２の量、即ち、吸着ＳＯ_２量がその最大量、即ち、飽和ＳＯ_２量を上回る場合は、ＤＯＣ２２ａが飽和していると考えられることから、総吸着余裕ＳＯ_２量（ｔ）の算出から除外される。また、この右方の分布において、現在床温Ｔ_２よりも低温側のデータが欠落している理由については、図５で説明した通りである。

そして、総吸着余裕ＳＯ_２量（ｔ）および総吸着余裕ＳＯ_３量（ｔ）を算出したら、これらを次式（１３）に代入して、第ｔ番目のサイクルにおける総吸着余裕量（Ｔ_２（ｔ），ｔ）を算出する。

更に、式（８）により算出した総飽和量（ｔ）と、式（１３）により算出した総吸着余裕量（ｔ）とを式（３）に代入すれば、飽和率（Ｔ_２（ｔ），ｔ）を算出できる。なお、以下の説明においては、飽和率（Ｔ_２（ｔ），ｔ）を単に飽和率（ｔ）ともいう。

図３に戻り目標床温Ｔｔｒｇの算出ロジックの説明を続ける。推定部Ｍ３は、「ＤＯＣ２２ａに流入するＳＯｘ」であってＤＯＣ２２ａに新たに吸着するＳＯｘの量（以下「新規吸着ＳＯｘ量」ともいう。）、および、「ＤＯＣ２２ａに流入するＳＯｘ」であってＤＯＣ２２ａに吸着することなくすり抜けるＳＯｘの量（以下「すり抜けＳＯｘ量」ともいう。）を推定する。先ず、新規吸着ＳＯｘ量とすり抜けＳＯｘ量の関係について、図７を参照して説明する。図７は、新規吸着ＳＯｘ量とすり抜けＳＯｘ量の関係を説明するための図である。この図に矢印で示すように、新規吸着ＳＯｘ量とすり抜けＳＯｘ量の和が、流入ＳＯｘ量に等しくなる。この理由は、「ＤＯＣ２２ａに流入するＳＯｘ」のうちの一部がＤＯＣ２２ａに吸着し、残りがＤＯＣ２２ａに吸着することなくすり抜けるためである。

推定部Ｍ３は、具体的に、推定部Ｍ１で推定した流入ＳＯｘ量（ｔ）と、推定部Ｍ２で推定した飽和率（ｔ）とを変数とする次式（１４）により新規吸着ＳＯｘ量を推定し、次式（１５）によりすり抜けＳＯｘ量を推定する。

なお、以下の説明においては、新規吸着ＳＯｘ量（流入ＳＯｘ量（ｔ），飽和率（ｔ））を単に新規吸着ＳＯｘ量（ｔ）ともいう。また、すり抜けＳＯｘ量（流入ＳＯｘ量（ｔ），飽和率（ｔ））を単にすり抜けＳＯｘ量（ｔ）ともいう。

式（１４）および（１５）の吸着率ｍａｐは、第ｔ番目のサイクルにおいて「ＤＯＣ２２ａに流入するＳＯｘ」のうち、ＤＯＣ２２ａに吸着するＳＯｘの割合（即ち、吸着率）が、飽和率（ｔ）によって変わるという特性に基づいて作成されたマップである。図８は、吸着率ｍａｐを説明するための図である。この図に示すように、吸着率ｍａｐの特性は、飽和率（ｔ）が低い領域では吸着率が高く、飽和率（ｔ）が高くなるにつれて吸着率が徐々に低下するというものである。なお、このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、現在床温Ｔ_２に応じて適宜読み出すことができる。

図３に戻り、推定部Ｍ４は、推定部Ｍ３で推定した新規吸着ＳＯｘ量を吸着ＳＯｘ分布に反映させて吸着後ＳＯｘ分布を推定する。吸着後ＳＯｘ分布について、ＳＯ_２を例とした図９を参照しながら説明する。図９は、吸着後ＳＯｘ分布を説明するための図である。この図に示すように、吸着後ＳＯ_２分布は、前回サイクル（例えば、第ｔ−１番目のサイクル）における最終吸着ＳＯ_２分布に、今回サイクル（例えば、第ｔ番目のサイクル）においてＤＯＣ２２ａに新たに吸着するＳＯ_２の量を表した分布（以下「新規吸着ＳＯ_２分布」ともいう。）を加算することで推定される。

推定部Ｍ４は、具体的には先ず、新規吸着ＳＯｘ量（ｔ）、総吸着余裕量（ｔ）および吸着余裕ＳＯ_２分布（ｔ）を変数とする次式（１６）により、第ｔ番目のサイクルにおける新規吸着ＳＯ_２分布を算出する。新規吸着ＳＯ_２分布と同様に、推定部Ｍ４は、ＤＯＣ２２ａに新たに吸着するＳＯ_３の量を表した分布（以下「新規吸着ＳＯ_３分布」ともいう。）を、次式（１７）により算出する。なお、吸着余裕ＳＯ_２分布（ｔ）と総吸着余裕量（ｔ）には、推定部Ｍ２において算出されたものが使用される。

なお、以下の説明においては、新規吸着ＳＯ_２分布（新規吸着ＳＯｘ量（ｔ），吸着余裕ＳＯ_２分布（ｔ），総吸着余裕量（ｔ））を、単に新規吸着ＳＯ_２分布（ｔ）ともいう。また、新規吸着ＳＯ_３分布（新規吸着ＳＯｘ量（ｔ），吸着余裕ＳＯ_３分布（ｔ），総吸着余裕量（ｔ））を、単に新規吸着ＳＯ_３分布（ｔ）ともいう。

推定部Ｍ４は、続いて、算出した新規吸着ＳＯ_２分布と、第ｔ−１番目のサイクルにおける最終吸着ＳＯ_２分布（ｔ−１）とを次式（１８）に代入して、吸着後ＳＯ_２分布を算出する。また、算出した新規吸着ＳＯ_３分布と、第ｔ−１番目のサイクルにおいて推定部Ｍ６で推定した吸着ＳＯ_３分布（ｔ−１）とを次式（１９）に代入して、吸着後ＳＯ_３分布を算出する。

図３に戻り、推定部Ｍ５は、推定部Ｍ４で推定した吸着後ＳＯｘ分布に基づいて、ＤＯＣ２２ａから新たに脱離するＳＯｘの量（以下「新規脱離ＳＯｘ量」ともいう。）を推定する。

推定部Ｍ５は、具体的には先ず、ＤＯＣ２２ａから脱離することのできるＳＯｘの総量（以下「脱離可能総ＳＯｘ量」ともいう。）を推定する。脱離可能総ＳＯｘ量について、ＳＯ_２を例とした図１０を参照しながら説明する。図１０は、脱離可能総ＳＯｘ量を説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、脱離可能総ＳＯｘ量は、現在床温Ｔ_２よりも低温側で、尚且つ、下限温度ＴＬよりも高温側の吸着後ＳＯｘ分布の面積に相当する。

ＤＯＣ２２ａから脱離することのできるＳＯ_２の総量、即ち、脱離可能総ＳＯ_２量は、現在床温Ｔ_２を変数とする次式（２０）により算出される。ＤＯＣ２２ａから脱離することのできるＳＯ_３の総量、即ち、脱離可能総ＳＯ_３量は、現在床温Ｔ_２を変数とする次式（２１）により算出される。

推定部Ｍ５は、算出した脱離可能総ＳＯ_２量を次式（２２）に代入して、第ｔ番目のサイクルにおいてＤＯＣ２２ａから新たに脱離するＳＯ_２の量、即ち、新規脱離ＳＯ_２量を算出する。また、算出した脱離可能総ＳＯ_３量を次式（２３）に代入して、第ｔ番目のサイクルにおいてＤＯＣ２２ａから新たに脱離するＳＯ_３の量、即ち、新規脱離ＳＯ_３量を算出する。なお、式（２２）および（２３）の脱離率には設定値が使用され、例えばＥＣＵ３０のＲＯＭに記憶させておくことができる。

図３に戻り、推定部Ｍ６は、推定部Ｍ５で推定した新規脱離ＳＯｘ量を吸着後ＳＯｘ分布に反映させて、最終吸着ＳＯｘ分布を推定する。

推定部Ｍ６は具体的に、推定部Ｍ５で推定した新規脱離ＳＯｘ量の分だけＳＯｘが脱離し、吸着後ＳＯｘ分布の形状が変形すると仮定して、最終吸着ＳＯｘ分布（脱離後ＳＯｘ分布）を推定する。最終吸着ＳＯｘ分布と吸着後ＳＯｘ分布の関係について、ＳＯ_２を例とした図１１を参照しながら説明する。図１１は、最終吸着ＳＯｘ分布と吸着後ＳＯｘ分布の関係を説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、吸着後ＳＯ_２分布の下限温度ＴＬからの積分値が新規脱離ＳＯ_２量に一致するときの吸着後ＳＯ_２分布の面積、即ち、下限温度ＴＬから温度Ｔｄ_ＳＯ２までの面積を、吸着後ＳＯｘ分布から削った後に残る分布が、最終吸着ＳＯ_２分布となる。

図１１の温度Ｔｄ_ＳＯ２が床温Ｔ_１を上回る場合には、ＤＯＣ２２ａからＳＯ_２が全て脱離していることを意味する。これを考慮すると、第ｔ番目のサイクルにおける最終吸着ＳＯ_２分布は床温Ｔ_１を変数とする次式（２４）で表され、第ｔ番目のサイクルにおける最終吸着ＳＯ_３分布は次式（２５）で表されることになる。なお、式（２５）の温度Ｔｄ_ＳＯ３は、吸着後ＳＯ_３分布の下限温度ＴＬからの積分値が新規脱離ＳＯ_３量に一致するときの床温Ｔ_１に相当している。

なお、新規脱離ＳＯ_２量と温度Ｔｄ_ＳＯ２の関係は次式（２６）で表すことができ、新規脱離ＳＯ_３量と温度Ｔｄ_ＳＯ３の関係は次式（２７）で表すことができる。

図３に戻り、推定部Ｍ７は、上述したすり抜けＳＯｘのうちＳＯ_３に転化するＳＯｘの量（以下「すり抜けＳＯ_３量」ともいう。）を推定する。

図２で説明したように、ＤＯＣ２２ａでは貴金属２２ｄに吸着したＳＯ_２の一部がＳＯ_３に転化する。この転化がすり抜けＳＯｘ中のＳＯ_２にも起こると仮定した場合、すり抜けＳＯｘがＳＯ_３に転化する割合は現在床温Ｔ_２とエンジン１０の排気流量（ガス流量）ＧＡに依存する。推定部Ｍ７では、すり抜け量と現在床温Ｔ_２とガス流量ＧＡを変数とする次式（２８）により、第ｔ番目のサイクルにおけるすり抜けＳＯ_３量を推定する。また、すり抜けＳＯｘのうちＳＯ_３に転化しないＳＯ_２の量（以下「すり抜けＳＯ_２量」ともいう。）は、次式（２９）で表すことができる。

式（２８）および（２９）のＳＯ_３転化率ｍａｐ（Ｔ_２（ｔ））は、第ｔ番目のサイクルにおいて「ＤＯＣ２２ａに流入するＳＯｘ」のうち、ＤＯＣ２２ａからＳＯ_３の状態で排出されるＳＯｘの割合（即ち、ＳＯ_３転化率）が、ＤＯＣ２２ａの現在床温Ｔ_２によって変わるという特性に基づいて作成されたマップである。図１２は、ＳＯ_３転化率ｍａｐを説明するための図である。ＳＯ_３転化率ｍａｐの特性は、例えば図１２に示すように、現在床温Ｔ_２がある温度域αにある場合はＳＯ_３転化率が高くなり、この温度域αよりも低温側では、ＳＯ_２からＳＯ_３への転化が起こり難くなるというものである。このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、現在床温Ｔ_２に応じて適宜読み出すことができる。

また、式（２８）および（２９）のＳＯ_３転化率補正ｍａｐ（ＧＡ（ｔ））は、第ｔ番目のサイクルにおいて、ＳＯ_３転化率がガス流量ＧＡによって変わるという特性に基づいて作成されたマップである。図１３は、ＳＯ_３転化率補正ｍａｐを説明するための図である。ＳＯ_３転化率補正ｍａｐの特性は、例えば図１３に示すように、ガス流量ＧＡがあるガス流量域βにある場合にはＳＯ_３転化率が高くなり、このガス流量域βよりも交流量側では、ガス流量ＧＡが高流量になるほどＳＯ_２からＳＯ_３への転化が起こり難くなるというものである。このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、ガス流量ＧＡに応じて適宜読み出すことができる。なお、ガス流量ＧＡは、例えばエアフローメータ３２の検出値を用いることができる。

図３に戻り、推定部Ｍ８は、上述したすり抜けＳＯｘのうちＤＯＣ２２ａにおいてＳＯ_３からＳＯ_２へと還元される割合（以下、「ＳＯ_２還元率」ともいう）を推定する。

ＤＯＣ２２ａではＳＯ_３の一部が還元反応によってＳＯ_２に転化する。ＳＯ_２への還元反応は、ＤＯＣ２２ａの還元雰囲気の影響を受ける。推定部Ｍ８では、エンジン１０の排気流量（ガス流量）ＧＡに対する排気添加量の割合ΔＦ／Ａを変数とする次式（３０）により、第ｔ番目のサイクルにおけるＳＯ_２還元率を推定する。

式（３０）のＳＯ_２還元率ｍａｐ（ΔＦ／Ａ（ｔ））は、第ｔ番目のサイクルにおいてＳＯ_２還元率が、ＤＯＣ２２ａの還元雰囲気の強さを示すΔＦ／Ａによって変わるという特性に基づいて作成されたマップである。図１４は、ＳＯ_２還元率ｍａｐを説明するための図である。ＳＯ_２還元率ｍａｐの特性は、例えば図１４に示すように、ΔＦ／Ａが０から高くなるにつれて、つまり還元雰囲気が強くなるにつれてＳＯ_２還元率が０から１まで次第に高くなり、それ移行は１近傍に維持される傾向を示している。このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、ΔＦ／Ａに応じて適宜読み出すことができる。なお、ΔＦ／Ａは、なお、すり抜け量とＤＯＣ２２ａに流入するガス流量ＧＡを変数とする次式（３１）により、第ｔ番目のサイクルにおけるΔＦ／Ａを算出することができる。

図３に戻り、算出部Ｍ９は、ＤＯＣ２２ａの床温上昇中にＤＯＣ２２ａから脱離してもよいＳＯ_３の量（以下「許容脱離ＳＯ_３量」ともいう。）を算出する。許容脱離ＳＯ_３量について、図１５を参照して説明する。図１５は、許容脱離ＳＯ_３量を説明するための図である。この図に示す制約ＳＯ_３量は、サルフェート白煙に関する制約に相当しており、この図においてはすり抜けＳＯ_３量と許容脱離ＳＯ_３量の和が制約ＳＯ_３量と等しくなっている。すり抜けＳＯ_３量と許容脱離ＳＯ_３量の和は、ＤＯＣ２２ａの下流におけるＳＯ_３の量であることから、この和が制約ＳＯ_３量となれば制約が満たされることになる。

但し、すり抜けＳＯ_３量と許容脱離ＳＯ_３量の和のうちの一部は、ＤＯＣ２２ａを流通する過程でＳＯ_２に還元される。従って、制約ＳＯ_３量とすり抜けＳＯ_３量とＳＯ_２還元率を変数とする次式（３２）を許容脱離ＳＯ_３量が満たせば、制約が満たされることになる。

なお、以下の説明においては、許容脱離ＳＯ_３量（制約ＳＯ_３量（ガス流量（ｔ）），すり抜けＳＯ_３量（ｔ），ＳＯ_２還元率（ΔＦ／Ａ（ｔ）））を単に許容脱離ＳＯ_３量（ｔ）ともいう。

図３に戻り、算出部Ｍ１０は、ＰＭ再生制御中のサルフェート白煙の発生を抑制するための第ｔ番目のサイクルにおける目標床温Ｔｔｒｇを算出する。目標床温Ｔｔｒｇについて、図１６を参照しながら説明する。図１６は、目標床温Ｔｔｒｇを説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、最終吸着ＳＯ_３分布の低温側からの積分値に脱離率を乗じた値が、算出部Ｍ９で算出した許容脱離ＳＯ_３量に一致するときの床温Ｔ_１が目標床温Ｔｔｒｇに相当する。

なお、第ｔ番目のサイクルにおける許容脱離ＳＯ_３量と目標床温Ｔｔｒｇの関係は次式（３３）で表すことができる。式（３３）の脱離率には設定値が使用され、例えばＥＣＵ３０のＲＯＭに記憶させておくことができる。

このように、以上説明した実施の形態１のシステムによれば、ＤＯＣ２２ａでのＳＯ３の還元反応をモデル化することにより、ＤＯＣ２２ａから排出されるＳＯ_３量の推定精度を高めることができる。これにより、制約ＳＯ_３量を守りつつ目標床温Ｔｔｒｇを上げることができるので、白煙を抑えつつＰＭの生成を早期に完了させて燃費悪化を抑えることが可能となる。

なお、上述した実施の形態においては、推定部Ｍ１が本発明の「流入ＳＯｘ量推定手段」に相当し、推定部Ｍ２、推定部Ｍ３、及び推定部Ｍ７が本発明の「すり抜けＳＯ_３量推定手段」に相当し、推定部Ｍ５及び推定部Ｍ６が本発明の「最終吸着ＳＯｘ分布推定手段」に相当し、算出部Ｍ８が本発明の「ＳＯ_２還元率推定手段」に相当し、算出部Ｍ９が本発明の「許容脱離ＳＯ_３量算出手段」に相当し、算出部Ｍ１０が本発明の「目標温度算出手段」に相当している。

また、上述した実施の形態においては、現在床温Ｔ_２が本発明の「代表温度」に相当している。

ところで、上述した実施の形態では、燃料添加弁２４からの燃料の添加によりＰＭ再生制御を行った。しかしこのＰＭ再生制御を、インジェクタ１２からの燃料の噴射（具体的には、メイン噴射よりも後のサブ噴射（例えばポスト噴射））により行ってもよい。この場合は、式（１）の排気添加量を、インジェクタ１２からのサブ噴射量に置き換えればよい。

また、上述した実施の形態では、ＰＭ再生制御中を例としてＤＯＣ２２ａの床温の目標温度を算出した。しかし、ＤＯＣ２２ａからＳＯｘを脱離させる制御をＰＭ再生制御と併せて行うような場合に、この脱離制御中に上述した手法によりＤＯＣ２２ａの床温の目標温度を算出してもよい。このように、上述した目標温度の算出手法は、ＤＯＣ２２ａからＳＯｘが脱離する温度域までＤＯＣ２２ａの床温を上昇させる制御一般に適用することができる。

また、上述した実施の形態では、ＤＯＣ２２ａとＤＰＦ２２ｂを備える排気浄化装置２２を例として説明した。しかし、ＤＯＣ２２ａにおけるＨＣ等の酸化機能をＤＰＦ２２ｂに付与して、排気浄化装置２２からＤＯＣ２２ａを省略してもよい。この場合は、上述した目標温度の算出手法を、酸化機能が付与されたＤＰＦ２２ｂに適用することで上述した実施の形態と同様の効果を得ることができる。

また、上述した実施の形態ではエンジン１０がターボチャージャ１６を備えるとしたが、エンジン１０がターボチャージャ１６を備えていなくてもよい。即ち、上述した目標温度の算出手法は、非過給ディーゼルエンジンのシステムにも適用できる。

さらに、実施の形態のシステムが備えるＥＣＵ３０は、以下のように構成されてもよい。図１７は、実施の形態の運転支援制御システムが備える制御装置のハードウェア構成の例を示す図である。ＥＣＵ３０の各機能は、処理回路により実現される。図１７に示す例では、ＥＣＵ３０の処理回路は、少なくとも１つのプロセッサ３０１と少なくとも１つのメモリ３０２とを備える。

処理回路が少なくとも１つのプロセッサ３０１と少なくとも１つのメモリ３０２とを備える場合、ＥＣＵ３０の各機能は、ソフトウェア、ファームウェア、またはソフトウェアとファームウェアとの組み合わせにより実現される。ソフトウェアおよびファームウェアの少なくとも一方は、プログラムとして記述される。ソフトウェアおよびファームウェアの少なくとも一方は、少なくとも１つのメモリ３０２に格納される。少なくとも１つのプロセッサ３０１は、少なくとも１つのメモリ３０２に記憶されたプログラムを読み出して実行することにより、ＥＣＵ３０の各機能を実現する。少なくとも１つのプロセッサ３０１は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、ＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）ともいう。例えば、少なくとも１つのメモリ３０２は、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、フラッシュメモリ、ＥＰＲＯＭ（Ｅｒａｓａｂｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ Ｅｒａｓａｂｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｒｅａｄ−Ｏｎｌｙ Ｍｅｍｏｒｙ）等の、不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク等である。

図１８は、実施の形態のシステムが備えるＥＣＵのハードウェア構成の他の例を示す図である。図１８に示す例では、ＥＣＵ３０の処理回路は、少なくとも１つの専用のハードウェア３０３を備える。

処理回路が少なくとも１つの専用のハードウェア３０３を備える場合、処理回路は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）、またはこれらを組み合わせたものである。ＥＣＵ３０の各部の機能がそれぞれ処理回路で実現されても良い。また、ＥＣＵ３０の各部の機能がまとめて処理回路で実現されても良い。

また、ＥＣＵ３０の各機能について、一部を専用のハードウェア３０３で実現し、他の一部をソフトウェアまたはファームウェアで実現してもよい。このように、処理回路は、ハードウェア３０３、ソフトウェア、ファームウェア、またはこれらの組み合わせによって、ＥＣＵ３０の各機能を実現する。

１０ ディーゼルエンジン
１２ インジェクタ
２０ 排気管
２２ 排気浄化装置
２２ａ ＤＯＣ
２２ｂ ＤＰＦ
２２ｃ コート材
２２ｄ 貴金属
２４ 燃料添加弁
３０ ＥＣＵ

Claims (6)

1. ディーゼルエンジンの排気管に設けられる浄化装置と、前記浄化装置に未燃燃料を供給する燃料供給手段と、前記燃料供給手段から未燃燃料を供給することにより前記浄化装置の温度を微粒子が燃焼する温度域の目標温度まで上昇させる昇温制御を実行する制御手段と、を有するエンジン制御装置であって、
   前記制御手段は、
   所定の制御周期ごとに前記浄化装置の温度の代表値である代表温度を取得する温度取得手段と、
   前記浄化装置に流入するＳＯｘ量を流入ＳＯｘ量として前記制御周期ごとに推定する流入ＳＯｘ量推定手段と、
   前記流入ＳＯｘ量と前記代表温度とを用いて、前記浄化装置の温度上昇中の各温度において前記浄化装置に最終的に吸着するＳＯｘ量を前記代表温度に関連付けたグラフとして表される最終吸着ＳＯｘ分布を前記制御周期ごとに推定する最終吸着ＳＯｘ分布推定手段と、
   前記流入ＳＯｘ量と前記代表温度とを用いて、前記浄化装置にＳＯｘの状態で流入して前記浄化装置に吸着することなくすり抜けるＳＯｘのうちＳＯ_３に転化するＳＯｘ量をすり抜けＳＯ_３量として前記制御周期ごとに推定するすり抜けＳＯ_３量推定手段と、
   前記浄化装置においてＳＯ_３からＳＯ_２へと還元される割合であるＳＯ_２還元率を推定するＳＯ_２還元率推定手段と、
   サルフェート白煙に関する制約に相当する前記浄化装置の下流における制約ＳＯ_３量と、前記すり抜けＳＯ_３量と、前記ＳＯ_２還元率とを用いて、前記浄化装置から脱離することが許されるＳＯ_３量を許容脱離ＳＯ_３量として前記制御周期ごとに算出する許容脱離ＳＯ_３量算出手段と、
   前記最終吸着ＳＯｘ分布と前記許容脱離ＳＯ_３量とを用いて、前記目標温度を前記制御周期ごとに算出する目標温度算出手段と、
   を備えることを特徴とするエンジン制御装置。
2. 前記ＳＯ_２還元率推定手段は、
   前記ＳＯ_２還元率と、前記燃料供給手段による未燃燃料の供給量と、前記浄化装置に流入するガス量との関係に基づいて、前記ＳＯ_２還元率を推定する
   ように構成されていることを特徴とする請求項１に記載のエンジン制御装置。
3. 前記すり抜けＳＯ_３量推定手段は、
   前記浄化装置の温度上昇中の各温度において前記浄化装置に吸着するＳＯｘ量を前記代表温度に関連付けたグラフとして表される吸着ＳＯｘ分布と、前記浄化装置の温度上昇中の各温度において前記浄化装置に吸着するＳＯｘ最大量を前記代表温度に関連付けたグラフとして表される飽和ＳＯｘ分布と、を用いて、前記浄化装置におけるＳＯｘ飽和率を前記制御周期ごとに推定し、
   前記流入ＳＯｘ量と前記ＳＯｘ飽和率とを用いて、前記浄化装置に流入して前記浄化装置に新たに吸着するＳＯｘ量を新規吸着ＳＯｘ量として前記制御周期ごとに推定し、
   前記新規吸着ＳＯｘ量を用いて、前記浄化装置に流入して前記浄化装置に吸着することなくすり抜けるＳＯｘ量をすり抜けＳＯｘ量として前記制御周期ごとに推定し、
   前記浄化装置においてＳＯ_３に転化するＳＯ_２の転化率と前記代表温度との関係を表した転化率マップと、前記代表温度と、前記すり抜けＳＯｘ量と、を用いて、前記すり抜けＳＯ_３量を前記制御周期ごとに推定する
   ように構成されていることを特徴とする請求項１又は２に記載のエンジン制御装置。
4. 前記最終吸着ＳＯｘ分布推定手段は、
   
   前記流入ＳＯｘ量と前記代表温度とを用いて、前記浄化装置から新たに脱離するＳＯｘ量を新規脱離ＳＯｘ量として前記制御周期ごとに推定し、
   前記新規脱離ＳＯｘ量を用いて、前記最終吸着ＳＯｘ分布を前記制御周期ごとに推定する
   ように構成されていることを特徴とする請求項１乃至３の何れか１項に記載のエンジン制御装置。
5. 前記温度取得手段は、
   前記排気管における前記浄化装置の下流側に流れたガスの温度を前記代表温度として取得する
   ように構成されていることを特徴とする請求項１乃至４の何れか１項に記載のエンジン制御装置。
6. 前記浄化装置は、前記排気管を流れる微粒子を捕集するフィルタを含み、
   前記制御手段は、
   前記昇温制御を、前記フィルタに捕集された微粒子量の推定値が除去要求量に到達したときに開始する
   ように構成されていることを特徴とする請求項１乃至５の何れか１項に記載のエンジン制御装置。

Images