(57)【要約】【課題】触媒の昇温制御において、触媒での還元反応を考慮することにより、脱離したＳＯｘに起因した白煙の発生の抑制しつつ昇温制御に要する時間の長期化を抑制することのできるエンジン制御装置を提供する。【解決手段】ＥＣＵ３０は、要素Ｍ１〜Ｍ１０によって制御周期ごとにＰＭ再生制御中のＤＯＣ２２ａの床温の目標温度を算出する。これらの要素のうち、推定部Ｍ７は、流入ＳＯｘ量と代表温度とを用いて、すり抜けＳＯ３量を制御周期ごとに推定する。推定部Ｍ８は、ＤＯＣ２２ａにおいてＳＯ３からＳＯ２へと還元される割合であるＳＯ２還元率を推定する。そして、算出部Ｍ９は、サルフェート白煙に関する制約に相当する制約ＳＯ３量と、すり抜けＳＯ３量と、ＳＯ２還元率とを用いて、ＤＯＣ２２ａから脱離することが許されるＳＯ３量である許容脱離ＳＯ３量を制御周期ごとに算出する。【選択図】図３
2019070376.tif 000002
【特許請求の範囲】
【請求項１】
ディーゼルエンジンの排気管に設けられる浄化装置と、前記浄化装置に未燃燃料を供給する燃料供給手段と、前記燃料供給手段から未燃燃料を供給することにより前記浄化装置の温度を微粒子が燃焼する温度域の目標温度まで上昇させる昇温制御を実行する制御手段と、を有するエンジン制御装置であって、
前記制御手段は、
所定の制御周期ごとに前記浄化装置の温度の代表値である代表温度を取得する温度取得手段と、
前記浄化装置に流入するＳＯｘ量を流入ＳＯｘ量として前記制御周期ごとに推定する流入ＳＯｘ量推定手段と、
前記流入ＳＯｘ量と前記代表温度とを用いて、前記浄化装置の温度上昇中の各温度において前記浄化装置に最終的に吸着するＳＯｘ量を前記代表温度に関連付けたグラフとして表される最終吸着ＳＯｘ分布を前記制御周期ごとに推定する最終吸着ＳＯｘ分布推定手段と、
前記流入ＳＯｘ量と前記代表温度とを用いて、前記浄化装置にＳＯｘの状態で流入して前記浄化装置に吸着することなくすり抜けるＳＯｘのうちＳＯ３に転化するＳＯｘ量をすり抜けＳＯ３量として前記制御周期ごとに推定するすり抜けＳＯ３量推定手段と、
前記浄化装置においてＳＯ３からＳＯ２へと還元される割合であるＳＯ２還元率を推定するＳＯ２還元率推定手段と、
サルフェート白煙に関する制約に相当する前記浄化装置の下流における制約ＳＯ３量と、前記すり抜けＳＯ３量と、前記ＳＯ２還元率とを用いて、前記浄化装置から脱離することが許されるＳＯ３量を許容脱離ＳＯ３量として前記制御周期ごとに算出する許容脱離ＳＯ３量算出手段と、
前記最終吸着ＳＯｘ分布と前記許容脱離ＳＯ３量とを用いて、前記目標温度を前記制御周期ごとに算出する目標温度算出手段と、
を備えることを特徴とするエンジン制御装置。
【請求項２】
前記ＳＯ２還元率推定手段は、
前記ＳＯ２還元率と、前記燃料供給手段による未燃燃料の供給量と、前記浄化装置に流入するガス量との関係に基づいて、前記ＳＯ２還元率を推定する
ように構成されていることを特徴とする請求項１に記載のエンジン制御装置。
【請求項３】
前記すり抜けＳＯ３量推定手段は、
前記浄化装置の温度上昇中の各温度において前記浄化装置に吸着するＳＯｘ量を前記代表温度に関連付けたグラフとして表される吸着ＳＯｘ分布と、前記浄化装置の温度上昇中の各温度において前記浄化装置に吸着するＳＯｘ最大量を前記代表温度に関連付けたグラフとして表される飽和ＳＯｘ分布と、を用いて、前記浄化装置におけるＳＯｘ飽和率を前記制御周期ごとに推定し、
前記流入ＳＯｘ量と前記ＳＯｘ飽和率とを用いて、前記浄化装置に流入して前記浄化装置に新たに吸着するＳＯｘ量を新規吸着ＳＯｘ量として前記制御周期ごとに推定し、
前記新規吸着ＳＯｘ量を用いて、前記浄化装置に流入して前記浄化装置に吸着することなくすり抜けるＳＯｘ量をすり抜けＳＯｘ量として前記制御周期ごとに推定し、
前記浄化装置においてＳＯ３に転化するＳＯ２の転化率と前記代表温度との関係を表した転化率マップと、前記代表温度と、前記すり抜けＳＯｘ量と、を用いて、前記すり抜けＳＯ３量を前記制御周期ごとに推定する
ように構成されていることを特徴とする請求項１又は２に記載のエンジン制御装置。
【請求項４】
前記最終吸着ＳＯｘ分布推定手段は、
前記流入ＳＯｘ量と前記代表温度とを用いて、前記浄化装置から新たに脱離するＳＯｘ量を新規脱離ＳＯｘ量として前記制御周期ごとに推定し、
前記新規脱離ＳＯｘ量を用いて、前記最終吸着ＳＯｘ分布を前記制御周期ごとに推定する
ように構成されていることを特徴とする請求項１乃至３の何れか１項に記載のエンジン制御装置。
【請求項５】
前記温度取得手段は、
前記排気管における前記浄化装置の下流側に流れたガスの温度を前記代表温度として取得する
ように構成されていることを特徴とする請求項１乃至４の何れか１項に記載のエンジン制御装置。
【請求項６】
前記浄化装置は、前記排気管を流れる微粒子を捕集するフィルタを含み、
前記制御手段は、
前記昇温制御を、前記フィルタに捕集された微粒子量の推定値が除去要求量に到達したときに開始する
ように構成されていることを特徴とする請求項１乃至５の何れか１項に記載のエンジン制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジン制御装置に関し、詳しくは、排気管に触媒を備えるディーゼルエンジンの制御装置に関する。
【背景技術】
【０００２】
従来、ディーゼルエンジンの排気管に設けられる触媒に吸着した硫黄酸化物（ＳＯ２またはＳＯ３をいい、以下においてこれらを区別しない場合には「ＳＯｘ」と総称する。）を、定期的に脱離させる昇温制御を行うことが知られている。昇温制御に関連する文献として、例えば特開２０１７−１０６３８１号公報が挙げられる。この公報には、触媒におけるＳＯｘの反応をモデル化することにより、触媒からＳＯ３の状態で排出されるＳＯｘ量を推定する技術が開示されている。ＳＯ３はＨ２０に溶解することでミスト化し、サルフェート白煙として視認される。上記の技術では、モデル演算によって推定された触媒から排出されるＳＯ３量がサルフェート白煙に関する制約に抵触しないように触媒の目標床温が決定される。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１７−１０６３８１号公報
【特許文献２】特開２０１５−１６９１０５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記公報の技術では、具体的には、触媒に流入したＳＯｘが触媒に吸着・脱離する反応と、触媒をすり抜けるＳＯｘのうちＳＯ２が酸化してＳＯ３に転化する反応をモデル化している。しかしながら、触媒におけるＳＯｘの反応は、触媒への吸着・脱離やＳＯ２の酸化反応だけではない。すなわち、触媒の昇温制御において排気に燃料が添加されると、触媒ではＳＯ３からＳＯ２への還元反応が起きる。上記公報の技術では、触媒でのＳＯ３の還元反応をモデル化していない。このため、上記公報の技術では、触媒からＳＯ３の状態で排出されるＳＯｘ量が実際よりも多く推定されるおそれがある。この場合、触媒の目標床温が実際に許容される温度よりも低く設定されてしまうため、床温制御の期間が長くなり燃費の悪化を招いてしまう。
【０００５】
本発明は、上述のような課題に鑑みてなされたもので、触媒の昇温制御において、触媒での還元反応を考慮することにより、脱離したＳＯｘに起因した白煙の発生の抑制しつつ昇温制御に要する時間の長期化を抑制することのできるエンジン制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記の課題を解決するため、本発明は、ディーゼルエンジンの排気管に設けられる浄化装置と、浄化装置に未燃燃料を供給する燃料供給手段と、燃料供給手段から未燃燃料を供給することにより浄化装置の温度を微粒子が燃焼する温度域の目標温度まで上昇させる昇温制御を実行する制御手段と、を有するエンジン制御装置に適用される。エンジン制御装置の制御手段は、所定の制御周期ごとに浄化装置の温度の代表値である代表温度を取得する温度取得手段と、浄化装置に流入するＳＯｘ量を流入ＳＯｘ量として制御周期ごとに推定する流入ＳＯｘ量推定手段と、流入ＳＯｘ量と代表温度とを用いて、浄化装置の温度上昇中の各温度において浄化装置に最終的に吸着するＳＯｘ量を代表温度に関連付けたグラフとして表される最終吸着ＳＯｘ分布を制御周期ごとに推定する最終吸着ＳＯｘ分布推定手段と、流入ＳＯｘ量と代表温度とを用いて、浄化装置にＳＯｘの状態で流入して浄化装置に吸着することなくすり抜けるＳＯｘのうちＳＯ３に転化するＳＯｘ量をすり抜けＳＯ３量として制御周期ごとに推定するすり抜けＳＯ３量推定手段と、浄化装置においてＳＯ３からＳＯ２へと還元される割合であるＳＯ２還元率を推定するＳＯ２還元率推定手段と、サルフェート白煙に関する制約に相当する浄化装置の下流における制約ＳＯ３量と、すり抜けＳＯ３量と、ＳＯ２還元率とを用いて、浄化装置から脱離することが許されるＳＯ３量を許容脱離ＳＯ３量として制御周期ごとに算出する許容脱離ＳＯ３量算出手段と、最終吸着ＳＯｘ分布と許容脱離ＳＯ３量とを用いて、目標温度を制御周期ごとに算出する目標温度算出手段と、を備えている。
【０００７】
また、本発明のエンジン制御装置において、ＳＯ２還元率推定手段は、ＳＯ２還元率と、燃料供給手段による未燃燃料の供給量と、浄化装置に流入するガス量との関係に基づいて、ＳＯ２還元率を推定するように構成されていてもよい。
【０００８】
また、本発明のエンジン制御装置において、すり抜けＳＯ３量推定手段は、以下のように構成されていてもよい。すなわち、すり抜けＳＯ３量推定手段は、浄化装置の温度上昇中の各温度において浄化装置に吸着するＳＯｘ量を代表温度に関連付けたグラフとして表される吸着ＳＯｘ分布と、浄化装置の温度上昇中の各温度において浄化装置に吸着するＳＯｘ最大量を代表温度に関連付けたグラフとして表される飽和ＳＯｘ分布と、を用いて、浄化装置におけるＳＯｘ飽和率を制御周期ごとに推定し、流入ＳＯｘ量とＳＯｘ飽和率とを用いて、浄化装置に流入して浄化装置に新たに吸着するＳＯｘ量を新規吸着ＳＯｘ量として制御周期ごとに推定し、新規吸着ＳＯｘ量を用いて、浄化装置に流入して浄化装置に吸着することなくすり抜けるＳＯｘ量をすり抜けＳＯｘ量として制御周期ごとに推定し、そして、浄化装置においてＳＯ３に転化するＳＯ２の転化率と代表温度との関係を表した転化率マップと、代表温度と、すり抜けＳＯｘ量と、を用いて、すり抜けＳＯ３量を制御周期ごとに推定するように構成されている。
【０００９】
また、本発明のエンジン制御装置において、最終吸着ＳＯｘ分布推定手段は、以下のように構成されていてもよい。すなわち、最終吸着ＳＯｘ分布推定手段は、流入ＳＯｘ量と代表温度とを用いて、浄化装置から新たに脱離するＳＯｘ量を新規脱離ＳＯｘ量として制御周期ごとに推定し、そして、最終吸着ＳＯｘ分布推定手段は、新規脱離ＳＯｘ量を用いて、最終吸着ＳＯｘ分布を制御周期ごとに推定するように構成されている。
【００１０】
また、本発明のエンジン制御装置において、温度取得手段は、排気管における浄化装置の下流側に流れたガスの温度を代表温度として取得するように構成されていてもよい。
【００１１】
また、本発明のエンジン制御装置において、浄化装置は、排気管を流れる微粒子を捕集するフィルタを含んで構成されていてもよい。そして、このような構成では、制御手段は、昇温制御を、フィルタに捕集された微粒子量の推定値が除去要求量に到達したときに開始するように構成されていてもよい。
【発明の効果】
【００１２】
本発明によれば、許容脱離ＳＯ３量の算出において、サルフェート白煙に関する制約に相当する浄化装置の下流における制約ＳＯ３量と、すり抜けＳＯ３量と、ＳＯ３からＳＯ２へと還元される割合であるＳＯ２還元率とが考慮される。ＳＯ２還元率が考慮されると、ＳＯ２還元率が考慮されない場合よりも許容脱離ＳＯ３量の推定精度を高めることができる。これにより、床温制御の目標温度を実際に許容される床温に近づけることができるので、脱離したＳＯｘに起因した白煙の発生の抑制しつつ昇温制御に要する時間の長期化を抑制することが可能となる。
【図面の簡単な説明】
【００１３】
【図１】実施の形態のシステム構成を示す図である。
【図２】ＤＯＣ２２ａにおけるＳＯｘの吸着と脱離を説明するための図である。
【図３】目標床温Ｔｔｒｇを算出するためのロジックを示す機能ブロック図である。
【図４】吸着ＳＯｘ分布と飽和ＳＯｘ分布を説明するための図である。
【図５】基準飽和ＳＯｘ分布と補正後の飽和ＳＯｘ分布の関係を説明するための図である。
【図６】総吸着余裕ＳＯ２量を説明するための図である。
【図７】新規吸着ＳＯｘ量とすり抜けＳＯｘ量の関係を説明するための図である。
【図８】吸着率ｍａｐを説明するための図である。
【図９】吸着後ＳＯｘ分布を説明するための図である。
【図１０】脱離可能総ＳＯｘ量を説明するための図である。
【図１１】最終吸着ＳＯｘ分布と吸着後ＳＯｘ分布の関係を説明するための図である。
【図１２】ＳＯ３転化率ｍａｐを説明するための図である。
【図１３】ＳＯ３転化率補正ｍａｐを説明するための図である。
【図１４】ＳＯ２還元率ｍａｐを説明するための図である。
【図１５】許容脱離ＳＯ３量を説明するための図である。
【図１６】目標床温Ｔｔｒｇを説明するための図である。
【図１７】実施の形態のシステムが備えるＥＣＵのハードウェア構成の例を示す図である。
【図１８】実施の形態のシステムが備えるＥＣＵのハードウェア構成の他の例を示す図である。
【発明を実施するための形態】
【００１４】
以下、図面を参照して本発明の実施の形態について説明する。ただし、以下に示す実施の形態において各要素の個数、数量、量、範囲等の数に言及した場合、特に明示した場合や原理的に明らかにその数に特定される場合を除いて、その言及した数に、この発明が限定されるものではない。また、以下に示す実施の形態において説明する構造やステップ等は、特に明示した場合や明らかに原理的にそれに特定される場合を除いて、この発明に必ずしも必須のものではない。
【００１５】
実施の形態１．
本発明の実施の形態１について図を参照して説明する。
【００１６】
［システム構成の説明］
図１は、実施の形態のシステム構成を示す図である。図１に示すシステムは、車両に搭載されるディーゼルエンジン１０（以下単に「エンジン１０」ともいう。）を備えたエンジン制御装置として構成されている。エンジン１０の各気筒には、燃料としての軽油を噴射するインジェクタ１２が設けられている。なお、図１に描かれるエンジン１０は直列４気筒エンジンであるが、エンジン１０の気筒数および気筒配列は特に限定されない。また、図１には、４つのインジェクタ１２のうちの１つが描かれている。
【００１７】
エンジン１０の排気マニホールド１４には、ターボチャージャ１６の排気タービン１６ａの入口が接続されている。排気タービン１６ａは吸気管１８に設けられたコンプレッサ１６ｂに連結されている。コンプレッサ１６ｂは、排気タービン１６ａの回転により駆動して吸気を過給する。排気タービン１６ａの出口には排気管２０が接続されている。排気管２０には排気浄化装置２２が設けられている。排気浄化装置２２は、ＤＯＣ（Diesel Oxidation Catalyst）２２ａと、ＤＰＦ（Diesel Particulate Filter）２２ｂと、を備えている。ＤＯＣ２２ａは、排気中の炭化水素（ＨＣ）や一酸化炭素（ＣＯ）を酸化して、水（Ｈ２Ｏ）や二酸化炭素（ＣＯ２）に転化させる機能を有する触媒である。ＤＰＦ２２ｂは、排気中に含まれる微粒子（ＰＭ）を捕集するフィルタである。排気浄化装置２２の上流には、インジェクタ１２と共通の未燃燃料を排気管２０に添加する燃料添加弁２４が設けられている。燃料添加弁２４は、燃料供給手段の例示である。
【００１８】
図１に示すシステムは、制御手段としてのＥＣＵ（Electronic Control Unit）３０を備えている。ＥＣＵ３０は、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、ＣＰＵ（Central Processing Unit）及び入出力インターフェースを備えるマイクロコンピュータである。ＥＣＵ３０は、入出力インターフェースを通して車両に搭載された各種センサの信号を取り込み処理する。各種センサには、吸気管１８の入口付近に設けられたエアフローメータ３２、ＤＯＣ２２ａの出口温度を検出する温度センサ３４、ＤＰＦ２２ｂの上下流における圧力差を検出する差圧センサ３６が含まれている。ＥＣＵ３０は、取り込んだ各センサの信号を処理して所定の制御プログラムに従って各種アクチュエータを操作する。ＥＣＵ３０によって操作されるアクチュエータには、上述したインジェクタ１２、燃料添加弁２４が含まれている。
【００１９】
［ＤＰＦ２２ｂの再生制御］
本実施の形態では、ＥＣＵ３０によるエンジン制御として、ＤＰＦ２２ｂの昇温制御（以下「ＰＭ再生制御」ともいう。）が行われる。ＰＭ再生制御は、ＤＰＦ２２ｂで捕集したＰＭの推定値が除去要求量に到達したときに、燃料添加弁２４から燃料を添加する制御である。例えば差圧センサ３６で検出した圧力差が所定値に到達したときに、ＰＭの推定値が除去要求量に到達したと判断することができる。燃料添加弁２４から燃料を添加することでＤＯＣ２２ａにおいて添加燃料を酸化し、この酸化反応熱によりＤＰＦ２２ｂの床温を６００℃以上まで上昇させる。これにより、ＤＰＦ２２ｂで捕集したＰＭを燃焼除去できるので、ＤＰＦ２２ｂの捕集機能を回復させることができる。なお、ＤＰＦ２２ｂの床温を６００℃以上に上昇させるための燃料添加弁２４からの添加燃料量（以下「ＤＰＦ用燃料量）は、ＤＰＦの床温と関連付けたマップに基づき決定されるものとする。このようなマップは例えばＥＣＵ３０のＲＯＭに記憶させておき、ＤＰＦ２２ｂの実際の床温に応じて適宜読み出すことができる。
【００２０】
［ＰＭ再生制御における問題点］
ところで、ディーゼルエンジンの燃料や潤滑油には一般に硫黄が含まれており、燃料の燃焼に伴ってこうした硫黄からＳＯｘが生成する。本実施の形態においても同様で、エンジン１０での燃料の燃焼に伴ってＳＯｘが生成する。生成したＳＯｘはエンジン１０から排出されて排気浄化装置２２に流入し、主にＤＯＣ２２ａに吸着する。但し、ＤＯＣ２２ａの床温が高くなると、ここに吸着していたＳＯｘが脱離し始める。ＤＯＣ２２ａの組成等により多少変動するものの、ＰＭ再生制御を行う温度域ではＤＯＣ２２ａからＳＯｘが脱離して下流側に放出される。
【００２１】
ＤＯＣ２２ａにおけるＳＯｘの吸着と脱離について、図２を参照しながら説明する。図２は、ＤＯＣ２２ａにおけるＳＯｘの吸着と脱離を説明するための図である。この図に示すように、ＤＯＣ２２ａは、基材（不図示）の表面を覆うコート材２２ｃと、貴金属２２ｄ（Ｐｔ，Ｐｄ等）とを備えている。貴金属２２ｄはコート材２２ｃに分散担持されており、ＨＣやＣＯを酸化する際の活性点となる。但し、排気中のＳＯ２が貴金属２２ｄに吸着し、または、排気中のＳＯ３がコート材２２ｃに吸着する。貴金属２２ｄに吸着したＳＯ２の一部は貴金属２２ｄから脱離して排気中に戻り、または、貴金属２２ｄ上で酸化されてＳＯ３となり、ＳＯ３の状態でコート材２２ｃに吸着する。即ち、貴金属２２ｄにはＳＯ２が吸着し、コート材２２ｃには排気由来のＳＯ３とＳＯ２由来のＳＯ３が吸着する。何れにせよ、ＳＯｘが吸着することで、ＤＯＣ２２ａにおけるＨＣ等の酸化機能が阻害されることになる。
【００２２】
上述した２経路によりコート材２２ｃに吸着したＳＯ３は、コート材２２ｃの床温が高くなることで脱離する。また、コート材２２ｃの床温が高くなることで貴金属２２ｄ上でのＳＯ２からＳＯ３への転化が促進されるので、このようなＳＯ３もコート材２２ｃから脱離する。従って、ＰＭ再生制御を行うことで、上述したＤＰＦ２２ｂの捕集機能だけでなく、ＤＯＣ２２ａにおけるＨＣ等の酸化機能を回復させることもできる。ところが、図２に示すように、コート材２２ｃから脱離したＳＯ３が排気管２０に存在するＨ２Ｏと反応することでＨ２ＳＯ４が発生する。そして、このＨ２ＳＯ４の濃度が一定濃度を超えると視認可能な白煙（サルフェート白煙）となるので、エンジン１０を搭載した車両の商品価値を損ねてしまうおそれがある。
【００２３】
［ＤＯＣ２２ａの温度制御］
ＤＯＣ２２ａの下流における排気中のＨ２ＳＯ４の濃度が高くなり過ぎないように燃料添加弁２４から燃料を添加すれば、ＰＭ再生制御中のサルフェート白煙の発生を抑制できる。そこで、本実施の形態のＥＣＵ３０は、制御モデルを用いてＤＯＣ２２ａの下流におけるＳＯ３の濃度を算出する。そして、ＥＣＵ３０は、算出したＳＯ３の濃度がサルフェート白煙に関する制約を満たすように、ＰＭ再生制御中のＤＯＣ２２ａの床温の目標温度（以下「目標床温Ｔｔｒｇ」ともいう。）を算出する。なお、このような制約ＳＯ３濃度（ＤＯＣ２２ａの下流におけるＳＯ３濃度の上限値）は、例えばＥＣＵ３０のＲＯＭに記憶させておくことができる。
【００２４】
ＥＣＵ３０は、目標床温Ｔｔｒｇに基づいて、燃料添加弁２４から添加する燃料量を算出する。以下の説明では、目標床温Ｔｔｒｇを実現するための燃料量を、制約充足用燃料量と称する。制約充足用燃料量よりもＤＰＦ用燃料量の方が多い場合に、ＤＰＦ用燃料量ではなく制約充足用燃料量を採用する。これにより、サルフェート白煙に関する制約を満たしつつ、ＤＯＣ２２ａにおけるＨＣ等の酸化機能を回復させることができる。
【００２５】
［本実施の形態の特徴］
本願の発明者は、ＰＭ再生制御中のＤＯＣ２２ａでのＳＯｘ反応について鋭意研究を重ねた。そして、ＰＭ再生制御中のＤＯＣ２２ａでは、ＳＯ２からＳＯ３への酸化反応だけでなく、ＳＯ３からＳＯ２への還元反応も活発に起きていることに着眼した。従来の技術のように、還元反応を考慮しない制御モデルでは、ＤＯＣ２２ａから排出されるＳＯ３量が実際よりも多く算出されてしまう。その結果、目標床温Ｔｔｒｇは、実際に許容される床温よりも低い温度に設定される。この場合、ＰＭ再生制御の実行期間が長期化し、燃費の悪化を招いてしまう。
【００２６】
本実施の形態のシステムは、ＰＭ再生制御中のＤＯＣ２２ａでのＳＯ３の還元反応を制御モデルに盛り込むこととした。これにより、ＤＯＣ２２ａから排出されるＳＯ３量の推定精度を高めることができるので、目標床温Ｔｔｒｇを実際に許容される床温に近づけることができる。その結果、ＰＭ再生制御の実行期間が長期化を抑制することができるので、燃費の悪化が抑制される。以下、本実施の形態のシステムにおいて実行される目標床温Ｔｔｒｇの算出ロジックについて詳細に説明する。
【００２７】
［目標床温Ｔｔｒｇの算出ロジック］
図３は、目標床温Ｔｔｒｇを算出するためのロジックを示す機能ブロック図であり、これはＥＣＵ３０により実現される。この図に示すように、ＥＣＵ３０は、流入ＳＯｘ量推定部Ｍ１と、ＳＯｘ飽和率推定部Ｍ２と、新規吸着ＳＯｘ量およびすり抜けＳＯｘ量推定部Ｍ３と、吸着後ＳＯｘ分布推定部Ｍ４と、新規脱離ＳＯｘ量推定部Ｍ５と、最終吸着ＳＯｘ分布推定部Ｍ６と、すり抜けＳＯ３量推定部Ｍ７と、ＳＯ２還元率推定部Ｍ８と、許容脱離ＳＯ３量算出部Ｍ９と、白煙抑制目標床温算出部Ｍ１０と、を備えており、これらの要素Ｍ１〜Ｍ９によって制御周期ごとに（具体的にはエンジン１０の燃焼サイクルごとに）目標床温Ｔｔｒｇを算出するものとする。なお、以下の説明においては、要素Ｍ１〜Ｍ９を簡略化するものとし、例えば流入ＳＯｘ量推定部Ｍ１を「推定部Ｍ１」ともいう。
【００２８】
推定部Ｍ１は、ＤＯＣ２２ａに流入するＳＯｘの量（以下「流入ＳＯｘ量」ともいう。）を推定する。なお、本明細書でいう「ＤＯＣ２２ａに流入するＳＯｘ」には、エンジン１０で生成し、ここから排出されてＤＯＣ２２ａに流入するＳＯｘだけでなく、燃料添加弁２４から添加された燃料のＤＯＣ２２ａでの酸化反応に伴って生成し、ＤＯＣ２２ａ上を流れるＳＯｘも含まれるものとする。
【００２９】
推定部Ｍ１は、具体的に、インジェクタ１２からの噴射燃料量（筒内噴射量）および燃料添加弁２４からの添加燃料量（排気添加量）を変数とする次式（１）により、第ｔ番目のサイクルにおける流入ＳＯｘ量を推定する。なお、式（１）の燃料Ｓ濃度は、燃料中の硫黄濃度であり、燃料供給系に別途設けた硫黄濃度センサの検出値を用いてもよく、設定値を用いてもよい。
【００３０】
【数１】
2019070376.tif 000003
【００３１】
式（１）の流入燃料量（排気添加量（ｔ），筒内噴射量（ｔ））は、「ＤＯＣ２２ａに流入するＳＯｘ」の由来となった燃料の第ｔ番目のサイクルにおける量であり、燃料の比重（軽油比重）を用いて次式（２）により算出される。
【００３２】
【数２】
2019070376.tif 000004
【００３３】
なお、以下の説明においては、流入ＳＯｘ量（排気添加量（ｔ），筒内噴射量（ｔ））を流入ＳＯｘ量（ｔ）ともいう。また、流入燃料量（排気添加量（ｔ），筒内噴射量（ｔ））を流入燃料量（ｔ）ともいう。
【００３４】
推定部Ｍ２は、ＤＯＣ２２ａにおけるＳＯｘの飽和率（以下「ＳＯｘ飽和率」ともいう。）を推定する。ＳＯｘ飽和率の推定には、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯｘの量（以下「吸着ＳＯｘ量」ともいう。）を、ＤＯＣ２２ａの床温に関連付けたグラフとして表される分布（以下「吸着ＳＯｘ分布」ともいう。）と、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯｘの最大量（以下「飽和ＳＯｘ量」ともいう。）を、ＤＯＣ２２ａの床温に関連付けたグラフとして表される分布（以下「飽和ＳＯｘ分布」ともいう。）と、が用いられる。先ず、吸着ＳＯｘ分布と飽和ＳＯｘ分布について、ＳＯ３を例として図４を参照しながら説明する。
【００３５】
図４は、吸着ＳＯｘ分布と飽和ＳＯｘ分布を説明するための図である。図４に「吸着ＳＯ３量」として示したデータは、次の手法により収集したものである。具体的には先ず、図４に「現在温度」として示した床温において、十分な量のＳＯｘをＤＯＣ２２ａに吸着させる。続いて、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａから脱離したＳＯ３の量を、上昇速度を一定とした条件のもとで測定する。そして、この脱離ＳＯ３量をＤＯＣ２２ａの床温に関連付けてグラフを作成する。これにより、脱離ＳＯ３量を表した分布（以下「脱離ＳＯ３分布」ともいう。）を得ることができる。これと同様の手法により、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａから脱離したＳＯ２の量を、ＤＯＣ２２ａの床温に関連付けたグラフ（以下「脱離ＳＯ２分布」ともいう。）を得ることもできる。なお、ＤＯＣ２２ａから脱離するＳＯ３については、これをセンサで直接的に測定してもよいし、ＳＯｘまたはＳＯ２を検出するセンサを用いて両者を測定し、これらの差から算出してもよい（ＳＯ３＝ＳＯｘ−ＳＯ２）。
【００３６】
ここで、ＤＯＣ２２ａの床温上昇中にＤＯＣ２２ａから脱離するＳＯ３は、実際には図４に「現在温度」として示した床温でＤＯＣ２２ａに吸着させたＳＯ３である。しかし、ある床温でＤＯＣ２２ａから脱離するＳＯ３は、その床温に至るまでＤＯＣ２２ａに吸着し続けることのできたＳＯ３であり、更に言えば、その床温でＤＯＣ２２ａに吸着することができたＳＯ３であると考えることもできる。このような仮定に基づいて、上述した脱離ＳＯ３分布の縦軸を、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯ３の量に置き換えると、図４に示した「吸着ＳＯ３量」のデータのグラフ、即ち、吸着ＳＯ３分布を得ることができる。そして、これと同様の手法によって、吸着ＳＯ２分布を得ることもできる。
【００３７】
また、図４に「飽和ＳＯ３量」として示したデータは、「吸着ＳＯ３量」のデータと同様の手法により収集したものである。この「飽和ＳＯ３量」のデータは、具体的に、上昇速度を極低速とした条件のもと、ＤＯＣ２２ａの床温上昇中の各床温（例えば５℃間隔）においてＤＯＣ２２ａから脱離したＳＯ３の量に相当している。ＤＯＣ２２ａの床温の上昇速度が極低速であることから、この「飽和ＳＯ３量」のデータは、ＤＯＣ２２ａから脱離したＳＯ３の量の最大値であると考えることができる。また、この最大値に対しては、上述した仮定を適用することができる。即ち、ある床温でＤＯＣ２２ａから脱離するＳＯ３の最大量は、その床温でＤＯＣ２２ａに吸着することができたＳＯ３の最大量に等しいと考えることができる。このような仮定に基づいて上述した脱離ＳＯ３分布の縦軸を上記ＳＯ３の最大量に置き換えると、図４に示した「飽和ＳＯ３量」のデータのグラフ、即ち、飽和ＳＯ３分布を得ることができる。そして、これと同様の手法によって、飽和ＳＯ２分布を得ることもできる。
【００３８】
推定部Ｍ２は、第ｔ番目のサイクルにおけるＤＯＣ２２ａの現在床温Ｔ２を変数とする次式（３）により、第ｔ番目のサイクルにおけるＳＯｘ飽和率（Ｔ２（ｔ），ｔ）を推定する。なお、現在床温Ｔ２は、ＤＯＣ２２ａの現在の床温の代表温度として、例えば温度センサ３４の検出値を使用することができる。
【００３９】
【数３】
2019070376.tif 000005
【００４０】
式（３）のＳＯｘ飽和率（Ｔ２（ｔ），ｔ）の算出過程は次のとおりである。先ず、ＤＯＣ２２ａの床温上昇中における床温Ｔ１と、現在床温Ｔ２とを変数とする次式（４）および（５）により、第ｔ番目のサイクルにおける飽和ＳＯ２分布（Ｔ１，Ｔ２（ｔ），ｔ）および飽和ＳＯ３分布（Ｔ１，Ｔ２（ｔ），ｔ）をそれぞれ算出する。
【００４１】
【数４】
2019070376.tif 000006
【００４２】
式（４）の基準飽和ＳＯ２分布は、十分な量のＳＯｘをＤＯＣ２２ａに吸着させるときの床温（図４の「現在温度」）を基準床温（例えば上述した吸着限界量が最大となる３００℃付近の床温）として作成した飽和ＳＯ２分布である。式（５）の基準飽和ＳＯ３分布もこれと同様である。式（４）の床温補正ＳＯ２ｍａｐ（Ｔ２（ｔ））は、基準飽和ＳＯ３分布を現在床温Ｔ２の飽和ＳＯ２分布に変換するための補正値を定めたマップである。式（５）の床温補正ＳＯ３ｍａｐ（Ｔ２（ｔ））もこれと同様である。このような基準飽和ＳＯｘ分布と補正マップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、現在床温Ｔ２に応じて適宜読み出すことができる。
【００４３】
基準飽和ＳＯｘ分布と補正後の飽和ＳＯｘ分布の関係を、ＳＯ２を例として図５を参照しながら説明する。図５は、基準飽和ＳＯｘ分布と補正後の飽和ＳＯｘ分布の関係を説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、ＤＯＣ２２ａの床温上昇中にＤＯＣ２２ａからＳＯ２が脱離し始める温度（下限温度）と、ＤＯＣ２２ａからＳＯ２が脱離し終わる温度（上限温度）に、それぞれ相当している。この図に示す３種類の分布の違いは、現在床温Ｔ２にある。即ち、現在床温Ｔ２が基準温度と等しい場合は、補正後の飽和ＳＯ２分布の形状が基準飽和ＳＯ２分布の形状と一致する（中央）。一方、現在床温Ｔ２が基準温度よりも低い場合（左方）や、現在床温Ｔ２が基準温度よりも高い場合（右方）は、補正後の飽和ＳＯ２分布の形状が基準飽和ＳＯ２分布の形状と一致しなくなる。なお、現在床温Ｔ２が基準温度よりも高い場合（右方）には、補正後の飽和ＳＯ２分布の形状が、現在床温Ｔ２よりも低温側のデータが欠落しているような形状となる。この理由は、現在床温Ｔ２よりも低温側では、本来であればこの床温域においてＤＯＣ２２ａに吸着し続けることのできたはずのＳＯｘが、ＤＯＣ２２ａから既に脱離していると考えられるためである。
【００４４】
続いて、式（４）により算出した飽和ＳＯ２分布（Ｔ１，Ｔ２（ｔ），ｔ）を次式（６）に代入して、第ｔ番目のサイクルにおける総飽和ＳＯ２量（Ｔ２（ｔ），ｔ）を算出する。また、式（５）により算出した飽和ＳＯ３分布（Ｔ１，Ｔ２（ｔ），ｔ）を次式（７）に代入して、第ｔ番目のサイクルにおける総飽和ＳＯ３量を算出する。
【００４５】
【数５】
2019070376.tif 000007
【００４６】
総飽和ＳＯ２量（Ｔ２（ｔ），ｔ）および総飽和ＳＯ３量（Ｔ２（ｔ），ｔ）を算出したら、次式（８）にこれらを代入して、第ｔ番目のサイクルにおける総飽和量（Ｔ２（ｔ），ｔ）を算出する。
【００４７】
【数６】
2019070376.tif 000008
【００４８】
なお、以下の説明においては、総飽和ＳＯ２量（Ｔ２（ｔ），ｔ）を単に総飽和ＳＯ２量（ｔ）ともいう。また、総飽和ＳＯ３量（Ｔ２（ｔ），ｔ）を単に総飽和ＳＯ３量（ｔ）ともいう。また、総飽和量（Ｔ２（ｔ），ｔ）を単に総飽和量（ｔ）ともいう。
【００４９】
式（８）により総飽和量（ｔ）を算出したら、飽和ＳＯ２分布（Ｔ１，Ｔ２（ｔ），ｔ）と、推定部Ｍ６で推定した第ｔ番目のサイクルにおける最終吸着ＳＯ２分布（Ｔ１，ｔ）とを次式（９）に代入して、第ｔ番目のサイクルにおける吸着余裕ＳＯ２分布（Ｔ１，Ｔ２（ｔ），ｔ）を算出する。また、飽和ＳＯ３分布（Ｔ１，Ｔ２（ｔ），ｔ）と、推定部Ｍ６で推定した第ｔ番目のサイクルにおける最終吸着ＳＯ３分布（Ｔ１，ｔ）とを次式（１０）に代入して、第ｔ番目のサイクルにおける吸着余裕ＳＯ３分布（Ｔ１，Ｔ２（ｔ），ｔ）を算出する。
【００５０】
【数７】
2019070376.tif 000009
【００５１】
続いて、式（９）により算出した吸着余裕ＳＯ２分布（Ｔ１，Ｔ２（ｔ），ｔ）を次式（１１）に代入して、第ｔ番目のサイクルにおける総吸着余裕ＳＯ２量（Ｔ２（ｔ），ｔ）を算出する。また、式（１０）により算出した吸着余裕ＳＯ３分布（Ｔ１，Ｔ２（ｔ），ｔ）を次式（１２）に代入して、第ｔ番目のサイクルにおける総吸着余裕ＳＯ３量（Ｔ２（ｔ），ｔ）を算出する。
【００５２】
【数８】
2019070376.tif 000010
【００５３】
なお、以下の説明においては、吸着余裕ＳＯ２分布（Ｔ１，Ｔ２（ｔ），ｔ）を単に吸着余裕ＳＯ２分布（ｔ）ともいう。また、吸着余裕ＳＯ３分布（Ｔ１，Ｔ２（ｔ），ｔ）を単に吸着余裕ＳＯ３分布（ｔ）ともいう。また、総吸着余裕ＳＯ２量（Ｔ２（ｔ），ｔ）を単に総吸着余裕ＳＯ２量（ｔ）ともいう。また、総吸着余裕ＳＯ３量（Ｔ２（ｔ），ｔ）を単に総吸着余裕ＳＯ３量（ｔ）ともいう。
【００５４】
図６を参照して、総吸着余裕ＳＯ２量（ｔ）を説明する。図６は、総吸着余裕ＳＯ２量を説明するための図である。なお、総吸着余裕ＳＯ３量（ｔ）についてはこれと同様である。この図に示すように、総吸着余裕ＳＯ２量（ｔ）は、飽和ＳＯ２分布から、飽和ＳＯ２分布と吸着ＳＯ２分布の重複部分を除いた面積として表すことができる。なお、この図の右方の分布に領域Ａとして示すように、ＤＯＣ２２ａの床温上昇中の各床温においてＤＯＣ２２ａに吸着するＳＯ２の量、即ち、吸着ＳＯ２量がその最大量、即ち、飽和ＳＯ２量を上回る場合は、ＤＯＣ２２ａが飽和していると考えられることから、総吸着余裕ＳＯ２量（ｔ）の算出から除外される。また、この右方の分布において、現在床温Ｔ２よりも低温側のデータが欠落している理由については、図５で説明した通りである。
【００５５】
そして、総吸着余裕ＳＯ２量（ｔ）および総吸着余裕ＳＯ３量（ｔ）を算出したら、これらを次式（１３）に代入して、第ｔ番目のサイクルにおける総吸着余裕量（Ｔ２（ｔ），ｔ）を算出する。
【００５６】
【数９】
2019070376.tif 000011
【００５７】
更に、式（８）により算出した総飽和量（ｔ）と、式（１３）により算出した総吸着余裕量（ｔ）とを式（３）に代入すれば、飽和率（Ｔ２（ｔ），ｔ）を算出できる。なお、以下の説明においては、飽和率（Ｔ２（ｔ），ｔ）を単に飽和率（ｔ）ともいう。
【００５８】
図３に戻り目標床温Ｔｔｒｇの算出ロジックの説明を続ける。推定部Ｍ３は、「ＤＯＣ２２ａに流入するＳＯｘ」であってＤＯＣ２２ａに新たに吸着するＳＯｘの量（以下「新規吸着ＳＯｘ量」ともいう。）、および、「ＤＯＣ２２ａに流入するＳＯｘ」であってＤＯＣ２２ａに吸着することなくすり抜けるＳＯｘの量（以下「すり抜けＳＯｘ量」ともいう。）を推定する。先ず、新規吸着ＳＯｘ量とすり抜けＳＯｘ量の関係について、図７を参照して説明する。図７は、新規吸着ＳＯｘ量とすり抜けＳＯｘ量の関係を説明するための図である。この図に矢印で示すように、新規吸着ＳＯｘ量とすり抜けＳＯｘ量の和が、流入ＳＯｘ量に等しくなる。この理由は、「ＤＯＣ２２ａに流入するＳＯｘ」のうちの一部がＤＯＣ２２ａに吸着し、残りがＤＯＣ２２ａに吸着することなくすり抜けるためである。
【００５９】
推定部Ｍ３は、具体的に、推定部Ｍ１で推定した流入ＳＯｘ量（ｔ）と、推定部Ｍ２で推定した飽和率（ｔ）とを変数とする次式（１４）により新規吸着ＳＯｘ量を推定し、次式（１５）によりすり抜けＳＯｘ量を推定する。
【００６０】
【数１０】
2019070376.tif 000012
【００６１】
なお、以下の説明においては、新規吸着ＳＯｘ量（流入ＳＯｘ量（ｔ），飽和率（ｔ））を単に新規吸着ＳＯｘ量（ｔ）ともいう。また、すり抜けＳＯｘ量（流入ＳＯｘ量（ｔ），飽和率（ｔ））を単にすり抜けＳＯｘ量（ｔ）ともいう。
【００６２】
式（１４）および（１５）の吸着率ｍａｐは、第ｔ番目のサイクルにおいて「ＤＯＣ２２ａに流入するＳＯｘ」のうち、ＤＯＣ２２ａに吸着するＳＯｘの割合（即ち、吸着率）が、飽和率（ｔ）によって変わるという特性に基づいて作成されたマップである。図８は、吸着率ｍａｐを説明するための図である。この図に示すように、吸着率ｍａｐの特性は、飽和率（ｔ）が低い領域では吸着率が高く、飽和率（ｔ）が高くなるにつれて吸着率が徐々に低下するというものである。なお、このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、現在床温Ｔ２に応じて適宜読み出すことができる。
【００６３】
図３に戻り、推定部Ｍ４は、推定部Ｍ３で推定した新規吸着ＳＯｘ量を吸着ＳＯｘ分布に反映させて吸着後ＳＯｘ分布を推定する。吸着後ＳＯｘ分布について、ＳＯ２を例とした図９を参照しながら説明する。図９は、吸着後ＳＯｘ分布を説明するための図である。この図に示すように、吸着後ＳＯ２分布は、前回サイクル（例えば、第ｔ−１番目のサイクル）における最終吸着ＳＯ２分布に、今回サイクル（例えば、第ｔ番目のサイクル）においてＤＯＣ２２ａに新たに吸着するＳＯ２の量を表した分布（以下「新規吸着ＳＯ２分布」ともいう。）を加算することで推定される。
【００６４】
推定部Ｍ４は、具体的には先ず、新規吸着ＳＯｘ量（ｔ）、総吸着余裕量（ｔ）および吸着余裕ＳＯ２分布（ｔ）を変数とする次式（１６）により、第ｔ番目のサイクルにおける新規吸着ＳＯ２分布を算出する。新規吸着ＳＯ２分布と同様に、推定部Ｍ４は、ＤＯＣ２２ａに新たに吸着するＳＯ３の量を表した分布（以下「新規吸着ＳＯ３分布」ともいう。）を、次式（１７）により算出する。なお、吸着余裕ＳＯ２分布（ｔ）と総吸着余裕量（ｔ）には、推定部Ｍ２において算出されたものが使用される。
【００６５】
【数１１】
2019070376.tif 000013
【００６６】
なお、以下の説明においては、新規吸着ＳＯ２分布（新規吸着ＳＯｘ量（ｔ），吸着余裕ＳＯ２分布（ｔ），総吸着余裕量（ｔ））を、単に新規吸着ＳＯ２分布（ｔ）ともいう。また、新規吸着ＳＯ３分布（新規吸着ＳＯｘ量（ｔ），吸着余裕ＳＯ３分布（ｔ），総吸着余裕量（ｔ））を、単に新規吸着ＳＯ３分布（ｔ）ともいう。
【００６７】
推定部Ｍ４は、続いて、算出した新規吸着ＳＯ２分布と、第ｔ−１番目のサイクルにおける最終吸着ＳＯ２分布（ｔ−１）とを次式（１８）に代入して、吸着後ＳＯ２分布を算出する。また、算出した新規吸着ＳＯ３分布と、第ｔ−１番目のサイクルにおいて推定部Ｍ６で推定した吸着ＳＯ３分布（ｔ−１）とを次式（１９）に代入して、吸着後ＳＯ３分布を算出する。
【００６８】
【数１２】
2019070376.tif 000014
【００６９】
図３に戻り、推定部Ｍ５は、推定部Ｍ４で推定した吸着後ＳＯｘ分布に基づいて、ＤＯＣ２２ａから新たに脱離するＳＯｘの量（以下「新規脱離ＳＯｘ量」ともいう。）を推定する。
【００７０】
推定部Ｍ５は、具体的には先ず、ＤＯＣ２２ａから脱離することのできるＳＯｘの総量（以下「脱離可能総ＳＯｘ量」ともいう。）を推定する。脱離可能総ＳＯｘ量について、ＳＯ２を例とした図１０を参照しながら説明する。図１０は、脱離可能総ＳＯｘ量を説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、脱離可能総ＳＯｘ量は、現在床温Ｔ２よりも低温側で、尚且つ、下限温度ＴＬよりも高温側の吸着後ＳＯｘ分布の面積に相当する。
【００７１】
ＤＯＣ２２ａから脱離することのできるＳＯ２の総量、即ち、脱離可能総ＳＯ２量は、現在床温Ｔ２を変数とする次式（２０）により算出される。ＤＯＣ２２ａから脱離することのできるＳＯ３の総量、即ち、脱離可能総ＳＯ３量は、現在床温Ｔ２を変数とする次式（２１）により算出される。
【００７２】
【数１３】
2019070376.tif 000015
【００７３】
推定部Ｍ５は、算出した脱離可能総ＳＯ２量を次式（２２）に代入して、第ｔ番目のサイクルにおいてＤＯＣ２２ａから新たに脱離するＳＯ２の量、即ち、新規脱離ＳＯ２量を算出する。また、算出した脱離可能総ＳＯ３量を次式（２３）に代入して、第ｔ番目のサイクルにおいてＤＯＣ２２ａから新たに脱離するＳＯ３の量、即ち、新規脱離ＳＯ３量を算出する。なお、式（２２）および（２３）の脱離率には設定値が使用され、例えばＥＣＵ３０のＲＯＭに記憶させておくことができる。
【００７４】
【数１４】
2019070376.tif 000016
【００７５】
図３に戻り、推定部Ｍ６は、推定部Ｍ５で推定した新規脱離ＳＯｘ量を吸着後ＳＯｘ分布に反映させて、最終吸着ＳＯｘ分布を推定する。
【００７６】
推定部Ｍ６は具体的に、推定部Ｍ５で推定した新規脱離ＳＯｘ量の分だけＳＯｘが脱離し、吸着後ＳＯｘ分布の形状が変形すると仮定して、最終吸着ＳＯｘ分布（脱離後ＳＯｘ分布）を推定する。最終吸着ＳＯｘ分布と吸着後ＳＯｘ分布の関係について、ＳＯ２を例とした図１１を参照しながら説明する。図１１は、最終吸着ＳＯｘ分布と吸着後ＳＯｘ分布の関係を説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、吸着後ＳＯ２分布の下限温度ＴＬからの積分値が新規脱離ＳＯ２量に一致するときの吸着後ＳＯ２分布の面積、即ち、下限温度ＴＬから温度ＴｄＳＯ２までの面積を、吸着後ＳＯｘ分布から削った後に残る分布が、最終吸着ＳＯ２分布となる。
【００７７】
図１１の温度ＴｄＳＯ２が床温Ｔ１を上回る場合には、ＤＯＣ２２ａからＳＯ２が全て脱離していることを意味する。これを考慮すると、第ｔ番目のサイクルにおける最終吸着ＳＯ２分布は床温Ｔ１を変数とする次式（２４）で表され、第ｔ番目のサイクルにおける最終吸着ＳＯ３分布は次式（２５）で表されることになる。なお、式（２５）の温度ＴｄＳＯ３は、吸着後ＳＯ３分布の下限温度ＴＬからの積分値が新規脱離ＳＯ３量に一致するときの床温Ｔ１に相当している。
【００７８】
【数１５】
2019070376.tif 000017
【００７９】
なお、新規脱離ＳＯ２量と温度ＴｄＳＯ２の関係は次式（２６）で表すことができ、新規脱離ＳＯ３量と温度ＴｄＳＯ３の関係は次式（２７）で表すことができる。
【００８０】
【数１６】
2019070376.tif 000018
【００８１】
図３に戻り、推定部Ｍ７は、上述したすり抜けＳＯｘのうちＳＯ３に転化するＳＯｘの量（以下「すり抜けＳＯ３量」ともいう。）を推定する。
【００８２】
図２で説明したように、ＤＯＣ２２ａでは貴金属２２ｄに吸着したＳＯ２の一部がＳＯ３に転化する。この転化がすり抜けＳＯｘ中のＳＯ２にも起こると仮定した場合、すり抜けＳＯｘがＳＯ３に転化する割合は現在床温Ｔ２とエンジン１０の排気流量（ガス流量）ＧＡに依存する。推定部Ｍ７では、すり抜け量と現在床温Ｔ２とガス流量ＧＡを変数とする次式（２８）により、第ｔ番目のサイクルにおけるすり抜けＳＯ３量を推定する。また、すり抜けＳＯｘのうちＳＯ３に転化しないＳＯ２の量（以下「すり抜けＳＯ２量」ともいう。）は、次式（２９）で表すことができる。
【００８３】
【数１７】
2019070376.tif 000019
【００８４】
式（２８）および（２９）のＳＯ３転化率ｍａｐ（Ｔ２（ｔ））は、第ｔ番目のサイクルにおいて「ＤＯＣ２２ａに流入するＳＯｘ」のうち、ＤＯＣ２２ａからＳＯ３の状態で排出されるＳＯｘの割合（即ち、ＳＯ３転化率）が、ＤＯＣ２２ａの現在床温Ｔ２によって変わるという特性に基づいて作成されたマップである。図１２は、ＳＯ３転化率ｍａｐを説明するための図である。ＳＯ３転化率ｍａｐの特性は、例えば図１２に示すように、現在床温Ｔ２がある温度域αにある場合はＳＯ３転化率が高くなり、この温度域αよりも低温側では、ＳＯ２からＳＯ３への転化が起こり難くなるというものである。このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、現在床温Ｔ２に応じて適宜読み出すことができる。
【００８５】
また、式（２８）および（２９）のＳＯ３転化率補正ｍａｐ（ＧＡ（ｔ））は、第ｔ番目のサイクルにおいて、ＳＯ３転化率がガス流量ＧＡによって変わるという特性に基づいて作成されたマップである。図１３は、ＳＯ３転化率補正ｍａｐを説明するための図である。ＳＯ３転化率補正ｍａｐの特性は、例えば図１３に示すように、ガス流量ＧＡがあるガス流量域βにある場合にはＳＯ３転化率が高くなり、このガス流量域βよりも交流量側では、ガス流量ＧＡが高流量になるほどＳＯ２からＳＯ３への転化が起こり難くなるというものである。このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、ガス流量ＧＡに応じて適宜読み出すことができる。なお、ガス流量ＧＡは、例えばエアフローメータ３２の検出値を用いることができる。
【００８６】
図３に戻り、推定部Ｍ８は、上述したすり抜けＳＯｘのうちＤＯＣ２２ａにおいてＳＯ３からＳＯ２へと還元される割合（以下、「ＳＯ２還元率」ともいう）を推定する。
【００８７】
ＤＯＣ２２ａではＳＯ３の一部が還元反応によってＳＯ２に転化する。ＳＯ２への還元反応は、ＤＯＣ２２ａの還元雰囲気の影響を受ける。推定部Ｍ８では、エンジン１０の排気流量（ガス流量）ＧＡに対する排気添加量の割合ΔＦ／Ａを変数とする次式（３０）により、第ｔ番目のサイクルにおけるＳＯ２還元率を推定する。
【００８８】
【数１８】
2019070376.tif 000020
【００８９】
式（３０）のＳＯ２還元率ｍａｐ（ΔＦ／Ａ（ｔ））は、第ｔ番目のサイクルにおいてＳＯ２還元率が、ＤＯＣ２２ａの還元雰囲気の強さを示すΔＦ／Ａによって変わるという特性に基づいて作成されたマップである。図１４は、ＳＯ２還元率ｍａｐを説明するための図である。ＳＯ２還元率ｍａｐの特性は、例えば図１４に示すように、ΔＦ／Ａが０から高くなるにつれて、つまり還元雰囲気が強くなるにつれてＳＯ２還元率が０から１まで次第に高くなり、それ移行は１近傍に維持される傾向を示している。このようなマップは、例えばＥＣＵ３０のＲＯＭに記憶させておくことができ、ΔＦ／Ａに応じて適宜読み出すことができる。なお、ΔＦ／Ａは、なお、すり抜け量とＤＯＣ２２ａに流入するガス流量ＧＡを変数とする次式（３１）により、第ｔ番目のサイクルにおけるΔＦ／Ａを算出することができる。
【００９０】
【数１９】
2019070376.tif 000021
【００９１】
図３に戻り、算出部Ｍ９は、ＤＯＣ２２ａの床温上昇中にＤＯＣ２２ａから脱離してもよいＳＯ３の量（以下「許容脱離ＳＯ３量」ともいう。）を算出する。許容脱離ＳＯ３量について、図１５を参照して説明する。図１５は、許容脱離ＳＯ３量を説明するための図である。この図に示す制約ＳＯ３量は、サルフェート白煙に関する制約に相当しており、この図においてはすり抜けＳＯ３量と許容脱離ＳＯ３量の和が制約ＳＯ３量と等しくなっている。すり抜けＳＯ３量と許容脱離ＳＯ３量の和は、ＤＯＣ２２ａの下流におけるＳＯ３の量であることから、この和が制約ＳＯ３量となれば制約が満たされることになる。
【００９２】
但し、すり抜けＳＯ３量と許容脱離ＳＯ３量の和のうちの一部は、ＤＯＣ２２ａを流通する過程でＳＯ２に還元される。従って、制約ＳＯ３量とすり抜けＳＯ３量とＳＯ２還元率を変数とする次式（３２）を許容脱離ＳＯ３量が満たせば、制約が満たされることになる。
【００９３】
【数２０】
2019070376.tif 000022
【００９４】
なお、以下の説明においては、許容脱離ＳＯ３量（制約ＳＯ３量（ガス流量（ｔ）），すり抜けＳＯ３量（ｔ），ＳＯ２還元率（ΔＦ／Ａ（ｔ）））を単に許容脱離ＳＯ３量（ｔ）ともいう。
【００９５】
図３に戻り、算出部Ｍ１０は、ＰＭ再生制御中のサルフェート白煙の発生を抑制するための第ｔ番目のサイクルにおける目標床温Ｔｔｒｇを算出する。目標床温Ｔｔｒｇについて、図１６を参照しながら説明する。図１６は、目標床温Ｔｔｒｇを説明するための図である。なお、この図の横軸のＴＬおよびＴＨは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、最終吸着ＳＯ３分布の低温側からの積分値に脱離率を乗じた値が、算出部Ｍ９で算出した許容脱離ＳＯ３量に一致するときの床温Ｔ１が目標床温Ｔｔｒｇに相当する。
【００９６】
なお、第ｔ番目のサイクルにおける許容脱離ＳＯ３量と目標床温Ｔｔｒｇの関係は次式（３３）で表すことができる。式（３３）の脱離率には設定値が使用され、例えばＥＣＵ３０のＲＯＭに記憶させておくことができる。
【００９７】
【数２１】
2019070376.tif 000023
【００９８】
このように、以上説明した実施の形態１のシステムによれば、ＤＯＣ２２ａでのＳＯ３の還元反応をモデル化することにより、ＤＯＣ２２ａから排出されるＳＯ３量の推定精度を高めることができる。これにより、制約ＳＯ３量を守りつつ目標床温Ｔｔｒｇを上げることができるので、白煙を抑えつつＰＭの生成を早期に完了させて燃費悪化を抑えることが可能となる。
【００９９】
なお、上述した実施の形態においては、推定部Ｍ１が本発明の「流入ＳＯｘ量推定手段」に相当し、推定部Ｍ２、推定部Ｍ３、及び推定部Ｍ７が本発明の「すり抜けＳＯ３量推定手段」に相当し、推定部Ｍ５及び推定部Ｍ６が本発明の「最終吸着ＳＯｘ分布推定手段」に相当し、算出部Ｍ８が本発明の「ＳＯ２還元率推定手段」に相当し、算出部Ｍ９が本発明の「許容脱離ＳＯ３量算出手段」に相当し、算出部Ｍ１０が本発明の「目標温度算出手段」に相当している。
【０１００】
また、上述した実施の形態においては、現在床温Ｔ２が本発明の「代表温度」に相当している。
【０１０１】
ところで、上述した実施の形態では、燃料添加弁２４からの燃料の添加によりＰＭ再生制御を行った。しかしこのＰＭ再生制御を、インジェクタ１２からの燃料の噴射（具体的には、メイン噴射よりも後のサブ噴射（例えばポスト噴射））により行ってもよい。この場合は、式（１）の排気添加量を、インジェクタ１２からのサブ噴射量に置き換えればよい。
【０１０２】
また、上述した実施の形態では、ＰＭ再生制御中を例としてＤＯＣ２２ａの床温の目標温度を算出した。しかし、ＤＯＣ２２ａからＳＯｘを脱離させる制御をＰＭ再生制御と併せて行うような場合に、この脱離制御中に上述した手法によりＤＯＣ２２ａの床温の目標温度を算出してもよい。このように、上述した目標温度の算出手法は、ＤＯＣ２２ａからＳＯｘが脱離する温度域までＤＯＣ２２ａの床温を上昇させる制御一般に適用することができる。
【０１０３】
また、上述した実施の形態では、ＤＯＣ２２ａとＤＰＦ２２ｂを備える排気浄化装置２２を例として説明した。しかし、ＤＯＣ２２ａにおけるＨＣ等の酸化機能をＤＰＦ２２ｂに付与して、排気浄化装置２２からＤＯＣ２２ａを省略してもよい。この場合は、上述した目標温度の算出手法を、酸化機能が付与されたＤＰＦ２２ｂに適用することで上述した実施の形態と同様の効果を得ることができる。
【０１０４】
また、上述した実施の形態ではエンジン１０がターボチャージャ１６を備えるとしたが、エンジン１０がターボチャージャ１６を備えていなくてもよい。即ち、上述した目標温度の算出手法は、非過給ディーゼルエンジンのシステムにも適用できる。
【０１０５】
さらに、実施の形態のシステムが備えるＥＣＵ３０は、以下のように構成されてもよい。図１７は、実施の形態の運転支援制御システムが備える制御装置のハードウェア構成の例を示す図である。ＥＣＵ３０の各機能は、処理回路により実現される。図１７に示す例では、ＥＣＵ３０の処理回路は、少なくとも１つのプロセッサ３０１と少なくとも１つのメモリ３０２とを備える。
【０１０６】
処理回路が少なくとも１つのプロセッサ３０１と少なくとも１つのメモリ３０２とを備える場合、ＥＣＵ３０の各機能は、ソフトウェア、ファームウェア、またはソフトウェアとファームウェアとの組み合わせにより実現される。ソフトウェアおよびファームウェアの少なくとも一方は、プログラムとして記述される。ソフトウェアおよびファームウェアの少なくとも一方は、少なくとも１つのメモリ３０２に格納される。少なくとも１つのプロセッサ３０１は、少なくとも１つのメモリ３０２に記憶されたプログラムを読み出して実行することにより、ＥＣＵ３０の各機能を実現する。少なくとも１つのプロセッサ３０１は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、ＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｏｒ）ともいう。例えば、少なくとも１つのメモリ３０２は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、フラッシュメモリ、ＥＰＲＯＭ（Ｅｒａｓａｂｌｅ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ  Ｅｒａｓａｂｌｅ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｒｅａｄ−Ｏｎｌｙ  Ｍｅｍｏｒｙ）等の、不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク等である。
【０１０７】
図１８は、実施の形態のシステムが備えるＥＣＵのハードウェア構成の他の例を示す図である。図１８に示す例では、ＥＣＵ３０の処理回路は、少なくとも１つの専用のハードウェア３０３を備える。
【０１０８】
処理回路が少なくとも１つの専用のハードウェア３０３を備える場合、処理回路は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）、またはこれらを組み合わせたものである。ＥＣＵ３０の各部の機能がそれぞれ処理回路で実現されても良い。また、ＥＣＵ３０の各部の機能がまとめて処理回路で実現されても良い。
【０１０９】
また、ＥＣＵ３０の各機能について、一部を専用のハードウェア３０３で実現し、他の一部をソフトウェアまたはファームウェアで実現してもよい。このように、処理回路は、ハードウェア３０３、ソフトウェア、ファームウェア、またはこれらの組み合わせによって、ＥＣＵ３０の各機能を実現する。
【符号の説明】
【０１１０】
１０  ディーゼルエンジン
１２  インジェクタ
２０  排気管
２２  排気浄化装置
２２ａ  ＤＯＣ
２２ｂ  ＤＰＦ
２２ｃ  コート材
２２ｄ  貴金属
２４  燃料添加弁
３０  ＥＣＵ
【図１】
2019070376.tif 000024
【図２】
2019070376.tif 000025
【図３】
2019070376.tif 000026
【図４】
2019070376.tif 000027
【図５】
2019070376.tif 000028
【図６】
2019070376.tif 000029
【図７】
2019070376.tif 000030
【図８】
2019070376.tif 000031
【図９】
2019070376.tif 000032
【図１０】
2019070376.tif 000033
【図１１】
2019070376.tif 000034
【図１２】
2019070376.tif 000035
【図１３】
2019070376.tif 000036
【図１４】
2019070376.tif 000037
【図１５】
2019070376.tif 000038
【図１６】
2019070376.tif 000039
【図１７】
2019070376.tif 000040
【図１８】
2019070376.tif 000041