(57)【特許請求の範囲】
【請求項１】
エンジンの運転状態を検出する運転状態検出手段と、
上記エンジンの排気通路に設けられた排ガス浄化用触媒に吸蔵されているＯ２量と相関するＯ２ストレージ指標を算出するＯ２ストレージ指標算出手段と、
上記Ｏ２ストレージ指標算出手段により算出されたＯ２ストレージ指標に基づき上記触媒上の空燃比を推定する空燃比推定手段と、
上記エンジンの燃料カットからの復帰時において、上記運転状態検出手段により検出されたエンジン運転状態及び上記空燃比推定手段により推定された空燃比に基づき上記触媒の浄化反応により発生する反応熱量を算出する浄化反応熱量算出手段と、
上記浄化反応熱量算出手段により算出された反応熱量に基づき上記燃料復帰時の触媒温度を推定する触媒温度推定手段と
を備えた触媒温度推定装置において、
上記エンジンの燃料カット中に上記触媒の酸化反応により消費された排ガス中のＯ２消費量を算出するＯ２消費量算出手段と、
上記燃料カットからの復帰時に、上記エンジンの排気空燃比に基づき排ガス中のＯ２過不足量を算出し、該Ｏ２過不足量から上記Ｏ２消費量算出手段により算出されたＯ２消費量を減算して上記排ガスを理論空燃比にするために要するＯ２必要量を算出するＯ２必要量算出手段と
をさらに備え、
上記Ｏ２ストレージ指標算出手段は、上記Ｏ２必要量算出手段により算出されたＯ２必要量に基づき上記Ｏ２ストレージ指標を算出することを特徴とする触媒温度推定装置。
【請求項２】
上記燃料復帰時に上記エンジンから排出される排ガスと上記触媒との間の熱伝達量を算出する熱伝達量算出手段を備え、
上記触媒温度推定手段は、上記熱伝達量算出手段により算出された熱伝達量と上記浄化反応熱量算出手段により算出された反応熱量とから上記触媒の温度を推定することを特徴とする請求項１記載の触媒温度推定装置。
【請求項３】
上記空燃比推定手段は、上記Ｏ２ストレージ指標算出手段により算出されたＯ２ストレージ指標に基づき上記触媒がＯ２吸蔵限界の上限近傍または下限近傍にあると判定したときには、排気空燃比検出手段により検出された排気空燃比を上記触媒上の空燃比と見なす一方、上記触媒が吸蔵限界の上限近傍または下限近傍にないと判定したときには、理論空燃比を上記触媒上の空燃比と見なすことを特徴とする請求項１または２記載の触媒温度推定装置。
【請求項４】
上記触媒の活性状態に基づき該触媒上のＯ２放出及び吸蔵能力を算出する触媒能力算出手段を備え、
上記Ｏ２ストレージ指標算出手段は、上記触媒能力算出手段により算出された触媒能力が上記Ｏ２必要量算出手段により算出されたＯ２必要量より小のときには、該触媒能力に基づき上記Ｏ２ストレージ指標を算出することを特徴とする請求項１乃至３の何れかに記載の触媒温度推定装置。
【請求項５】
上記触媒温度推定手段は、推定した触媒温度を上記触媒中の熱伝導を模擬した３次フィルタにより補正することを特徴とする請求項１乃至４の何れかに記載の触媒温度推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、排気浄化用触媒の温度を推定する触媒温度推定装置に関するものである。
【背景技術】
【０００２】
エンジンの排気系には排ガス中に含まれるＨＣ（炭化水素）、ＣＯ（一酸化炭素）、ＮＯｘ（窒素酸化物）などの有害物質を浄化するための排気浄化用触媒が設けられている。この種の触媒の温度はエンジンの運転状態に応じて大きく変化するため、いかなる運転状態でも触媒の耐熱温度を越えないように、触媒温度の推定値に基づいて適切にエンジンを制御する必要がある。触媒温度の推定手法としては、定常運転時の排気温度を記憶しておき、過渡運転時には記憶している定常値をフィルタ処理して排気温度を推定している。
【０００３】
ところで、車両減速時に実行される燃料カット時には、排ガス中のＯ２（酸素）と触媒上に吸着されているＣＯやＨＣとが反応することで多量の酸化反応熱が発生し、燃料復帰時には、燃料カット中に触媒上に吸蔵されたＯ２により排ガス中のＣＯやＨＣが浄化されることで多量の浄化反応熱が発生する。従って、何れの場合も触媒の急激な温度上昇を引き起こすため、特に触媒温度を推定し難い条件となる。上記触媒温度の推定手法では、フィルタ処理により温度変化率を定常運転時より小さく推定しているだけのため、特に燃料カット及び燃料復帰が頻繁に繰り返される状況では、図７のタイムチャートに破線で示すように、実際の触媒温度からかけ離れた触媒温度が推定されてしまうという問題があった。
【０００４】
一方、上記とは別の触媒温度の推定手法も提案されている（例えば、特許文献１参照）。当該特許文献１の技術では、燃料カット以外の運転状態においては、エンジンから排出される排気温度を定常触媒温度と見なして、定常触媒温度から排気管壁温度を推定し、これらの定常触媒温度と排気管壁温度とから仮触媒流入排気温度を推定した上で、走行風による排気温度の低下量を見込んで仮触媒流入排気温度から推定した触媒流入排気温度に基づき触媒温度を推定する一方、燃料カット時には、以上の推定処理を中止して触媒温度の前回値をそのまま今回値としている。
【特許文献１】特開２００４−２６３６０６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上記特許文献１の触媒温度の推定手法のように燃料カット時に前回値をそのまま今回値として適用することは、燃料カット及び燃料復帰直時に触媒温度を変化しないものと推定していることを意味し、その推定値は例えば図７に示す従来手法と同様に、ＨＣやＣＯの酸化反応により急激に上昇する実際の触媒温度に対して大きく相違してしまう。よって、当該推定手法により推定された触媒温度に基づき燃料カット時や燃料復帰時の触媒温度の上昇を抑制できるようにエンジン制御を実行することは到底不可能であり、結果として耐熱温度を越えた温度上昇により触媒が破損する虞があった。
【０００６】
本発明はこのような問題点を解決するためになされたもので、その目的とするところは、エンジンの燃料復帰時の触媒温度の挙動を正確に推定でき、もって耐熱温度を越えた温度上昇による触媒の破損を未然に防止することができる触媒温度推定装置を提供することにある。
【課題を解決するための手段】
【０００７】
上記目的を達成するため、請求項１の発明は、エンジンの運転状態を検出する運転状態検出手段と、エンジンの排気通路に設けられた排ガス浄化用触媒に吸蔵されているＯ２量と相関するＯ２ストレージ指標を算出するＯ２ストレージ指標算出手段と、Ｏ２ストレージ指標算出手段により算出されたＯ２ストレージ指標に基づき触媒上の空燃比を推定する空燃比推定手段と、エンジンの燃料カットからの復帰時において、運転状態検出手段により検出されたエンジン運転状態及び空燃比推定手段により推定された空燃比に基づき触媒の浄化反応により発生する反応熱量を算出する浄化反応熱量算出手段と、浄化反応熱量算出手段により算出された反応熱量に基づき燃料復帰時の触媒温度を推定する触媒温度推定手段とを備えた触媒温度推定装置において、エンジンの燃料カット中に触媒の酸化反応により消費された排ガス中のＯ２消費量を算出するＯ２消費量算出手段と、燃料カットからの復帰時に、エンジンの排気空燃比に基づき排ガス中のＯ２過不足量を算出し、Ｏ２過不足量からＯ２消費量算出手段により算出されたＯ２消費量を減算して排ガスを理論空燃比にするために要するＯ２必要量を算出するＯ２必要量算出手段とをさらに備え、Ｏ２ストレージ指標算出手段が、Ｏ２必要量算出手段により算出されたＯ２必要量に基づきＯ２ストレージ指標を算出するものである。
【０００８】
従って、エンジンの燃料カット中に触媒の酸化反応により消費された排ガス中のＯ２消費量がＯ２消費量算出手段により算出され、燃料カットからの復帰時には、Ｏ２必要量算出手段によりエンジンの排気空燃比に基づき排ガス中のＯ２過不足量が算出されると共に、このＯ２過不足量からＯ２消費量算出手段により算出されたＯ２消費量が減算されて排ガスを理論空燃比にするために要するＯ２必要量が算出される。
そして、Ｏ２必要量に基づき、触媒に吸蔵されているＯ２量と相関するＯ２ストレージ指標がＯ２ストレージ指標算出手段により算出され、このＯ２ストレージ指標から推定された触媒上の空燃比及び運転状態検出手段により検出されたエンジンの運転状態に基づき触媒の浄化反応で発生する反応熱量が浄化反応熱量算出手段により算出され、算出された反応熱量に基づき触媒温度推定手段により触媒温度が推定される。
【０００９】
燃料復帰時には、燃料カット中に触媒上に吸蔵されたＯ２により排ガス中のＣＯやＨＣが浄化されることで多量の反応熱が発生し、この反応熱が触媒の温度を急激に上昇させる要因となるが、このような現象を反映した反応熱量に基づき触媒温度の推定処理が実行されることで、燃料復帰時の触媒温度の挙動を正確に推定可能となる。
請求項２の発明は、請求項１において、燃料復帰時にエンジンから排出される排ガスと触媒との間の熱伝達量を算出する熱伝達量算出手段を備え、触媒温度推定手段が、熱伝達量算出手段により算出された熱伝達量と浄化反応熱量算出手段により算出された反応熱量とから触媒の温度を推定するものである。
【００１０】
従って、浄化反応熱量に加えて排ガスと触媒との間の熱伝達量を考慮して触媒の温度が推定されることから、燃料復帰時の触媒温度の挙動を一層正確に推定可能となる。
請求項３の発明は、請求項１または２において、空燃比推定手段が、Ｏ２ストレージ指標算出手段により算出されたＯ２ストレージ指標に基づき触媒がＯ２吸蔵限界の上限近傍または下限近傍にあると判定したときには、排気空燃比検出手段により検出された排気空燃比を触媒上の空燃比と見なす一方、触媒が吸蔵限界の上限近傍または下限近傍にないと判定したときには、理論空燃比を触媒上の空燃比と見なすものである。
【００１１】
従って、Ｏ２ストレージ指標に基づき触媒がＯ２吸蔵限界の上限近傍または下限近傍にあると判定したときには、排気空燃比が触媒上の空燃比と見なされて反応熱量の算出に適用される一方、触媒が吸蔵限界の上限近傍または下限近傍にないと判定したときには、理論空燃比が触媒上の空燃比と見なされて反応熱量の算出に適用される。
触媒がＯ２ストレージ機能を発揮しているときには、排気空燃比Ａ／Ｆがストイキオ以外であっても触媒上の空燃比は理論空燃比に保持される一方、触媒がＯ２ストレージ機能を発揮していないときには、触媒上の空燃比は排気空燃比となることから、このような現象を反映した実状に則した正確な触媒上の空燃比に基づいて反応熱量、ひいては燃料復帰時の触媒温度の挙動を一層正確に推定可能となる。
【００１３】
請求項４の発明は、請求項１乃至３において、触媒の活性状態に基づき触媒上のＯ２放出及び吸蔵能力を算出する触媒能力算出手段を備え、Ｏ２ストレージ指標算出手段が、触媒能力算出手段により算出された触媒能力がＯ２必要量算出手段により算出されたＯ２必要量より小のときには、触媒能力に基づき上記Ｏ２ストレージ指標を算出するものである。
【００１４】
従って、排ガスを理論空燃比にするために要するＯ２必要量に対して触媒能力が小のときには触媒能力に基づきＯ２ストレージ指標が算出される。触媒能力以上のＯ２必要量が要求されても、実際に使用されるＯ２量は触媒能力相当分に制限されることから、この場合には触媒能力に基づき実情に即したＯ２ストレージ量の算出が可能となる。
請求項５の発明は、請求項１乃至４において、触媒温度推定手段が、推定した触媒温度を触媒中の熱伝導を模擬した３次フィルタにより補正するものである。
【００１５】
従って、触媒中の熱伝導を模擬した３次フィルタ補正により触媒温度の挙動が一層正確に推定される。
【発明の効果】
【００１６】
以上説明したように請求項１の発明の触媒温度推定装置によれば、エンジン運転状態、Ｏ２ストレージ指標、触媒上の空燃比に基づき触媒の浄化反応で発生する反応熱量を算出して触媒の温度を推定するため、エンジンの燃料復帰時の触媒温度の挙動を正確に推定でき、もって耐熱温度を越えた温度上昇による触媒の破損を未然に防止できると共に、排ガスを理論空燃比にするためのＯ２必要量に基づきＯ２ストレージ指標を算出することで、さらに正確に反応熱量を算出することができる。
請求項２の発明の触媒温度推定装置によれば、請求項１に加えて、排ガスと触媒との間の熱伝達量を考慮することで、触媒温度の挙動を一層正確に推定することができる。
【００１７】
請求項３の発明の触媒温度推定装置によれば、請求項１または２に加えて、触媒のＯ２ストレージ機能の有無に応じた正確な触媒上の空燃比に基づき反応熱量を算出するため、燃料復帰時の触媒温度の挙動を一層正確に推定することができる。
【００１８】
請求項４の発明の触媒温度推定装置によれば、請求項１乃至３に加えて、触媒能力を考慮して実情に即したＯ２ストレージ指標を算出することで、さらに正確に反応熱量を算出することができる。
請求項５の発明の触媒温度推定装置によれば、請求項１乃至４に加えて、３次フィルタ補正により触媒温度の挙動を一層正確に推定することができる。
【発明を実施するための最良の形態】
【００１９】
以下、本発明を具体化した触媒温度推定装置の一実施例を説明する。
図１は本実施形態の触媒温度推定装置を示す全体構成図であり、エンジン１は筒内噴射型火花点火式の直列４気筒ガソリンエンジンとして構成されている。エンジン１のシリンダヘッドには、各気筒毎に点火プラグ２と共に電磁式の燃料噴射弁３が取り付けられ、この燃料噴射弁３から燃焼室内に直接燃料が噴射される。シリンダヘッドには各気筒毎に略直立方向に吸気ポート４ａが形成されており、これらの吸気ポート４ａは吸気マニホールド４を介してスロットルバルブ５と接続され、スロットルバルブ５は図示しない吸気通路と接続されている。
【００２０】
又、シリンダヘッドには略水平方向に排気ポート６ａが形成されており、これらの排気ポート６ａには排気マニホールド６を介して排気通路７が接続されている。排気通路７の上流側には前段触媒８が配置され、下流側には床下触媒９が配置されている。前段触媒８及び床下触媒９は白金（Ｐｔ）、ロジウム（Ｒｈ）などの貴金属を担持した三元触媒として構成され、本実施形態では床下触媒９を対象として温度推定の処理を実行している。なお、温度推定の対象となる触媒は床下触媒９に限ることはなく、例えば前段触媒８や図示しないＮＯx浄化用のＮＯx吸蔵触媒を対象として温度を推定するようにしてもよい。床下触媒９の上流側には空燃比センサ１０（排気空燃比検出手段）が配設され、空燃比センサ１０は床下触媒９に流入する排ガスの空燃比を検出する。
【００２１】
車室内には入出力装置、記憶装置（ＲＯＭ、ＲＡＭ、不揮発性ＲＡＭ等）、中央処理装置（ＣＰＵ）、タイマカウンタ等を備えたＥＣＵ（電子コントロールユニット）２１が設置されている。ＥＣＵ２１の入力側には上記空燃比センサ１０、スロットルバルブ５の開度を検出するスロットルセンサ２２、エンジン１の回転速度を検出する回転速度センサ２３、冷却水温を検出する水温センサ２４などの各種センサ類が接続されており、これらセンサ類からの検出情報が入力される。ＥＣＵ２１の出力側には、上記点火プラグ２、燃料噴射弁３などの各種デバイス類が接続されている。
【００２２】
ＥＣＵ２１は各センサからの検出情報に基づいて点火時期や燃料噴射時期などを決定し、決定した制御量に基づいて点火プラグ２や燃料噴射弁３を駆動制御してエンジン１を運転する。また、アクセルオフによる車両減速及び所定車速以上などの所定の燃料カット条件が成立すると燃料カットにより燃料噴射を中止し、アクセルオンなどの所定の燃料復帰条件が成立すると燃料噴射を再開する。
【００２３】
そして、[背景技術]でも述べたように、燃料カット時及び燃料復帰時には床下触媒９の急激な温度上昇が発生するため触媒温度Ｔcatfを推定し難い条件となり、必然的に触媒温度Ｔcatfに基づいて温度抑制のための適切なエンジン制御を実行することも困難となる。そこで、本実施形態では燃料カット時及び燃料復帰時の反応熱を反映した床下触媒９の温度推定処理を実行しており、以下、当該温度推定処理について詳述する。
【００２４】
具体的な説明に先立って、燃料カット時及び燃料復帰時に触媒上で生起される反応と、当該反応に起因する触媒温度Ｔcatfへの影響について述べる。
燃料カット時及び燃料復帰時においても他の運転状態と同様に、基本的に触媒の昇温は排ガスと床下触媒９との間の熱伝達量ΔＨtに応じて行われるため、触媒温度Ｔcatfの推定には熱伝達量ΔＨtを考慮する必要がある。この熱伝達量ΔＨtによる触媒昇温に加えて、燃料カット時には、燃料カット以前に触媒９上に吸着された未燃燃料であるＣＯやＨＣ（共にリッチ成分であり、以下の説明では代表としてＣＯと称するがＨＣも含むものとする）が排ガス中のＯ２により酸化反応したときの反応熱で触媒温度Ｔcatfが上昇する現象が生じることから、このときの酸化反応熱量ΔＨrを触媒温度Ｔcatfの推定に反映させる必要がある。
【００２５】
また、燃料復帰時には、燃料カット中に触媒９に吸蔵されたＯ２により排ガス中の未燃燃料であるＣＯやＨＣ（上記と同じく代表としてＣＯと称する）が浄化されたときの反応熱で触媒温度Ｔcatfが上昇する現象が生じることから、このときの浄化反応熱量ΔＨcを触媒温度Ｔcatfの推定に反映させる必要がある。従って、以下に述べるように、燃料カット時には、排ガスと触媒９との間の熱伝達量ΔＨtに加えて酸化反応熱量ΔＨrを考慮して触媒温度Ｔcatfを推定し、燃料復帰時には、熱伝達量ΔＨtに加えて浄化反応熱量ΔＨcを考慮して触媒温度Ｔcatfを推定している。
【００２６】
まず、上記燃料カット時のＣＯ酸化反応熱量ΔＨrの算出処理について説明する。
図２はＥＣＵ２１のＦ／Ｃ時反応熱量算出部５０の処理を示す制御ブロック図である。燃料カット時には排ガス中のＯ２との反応により触媒上からＣＯが脱離する現象のみならず、触媒上にＣＯが吸着する現象も並行して発生し、これらの現象が相俟って触媒９のＣＯ吸着率θco、ひいては燃料カット時の酸化反応熱量ΔＨrに影響を及ぼす。そこで、当該ＣＯ酸化反応熱量ΔＨrの算出処理では、触媒上からのＣＯの脱離及び触媒上へのＣＯの吸着を考慮して触媒９の吸着ＣＯ量CATcoを算出し、この吸着ＣＯ量CATcoをＣＯ酸化反応熱量ΔＨrの算出処理に反映させている。
【００２７】
図２においてＦ／Ｃ時反応熱量算出部５０の脱離定数算出部５１では、前回処理で推定された触媒温度Ｔcatfから予め設定されたマップに従って脱離速度定数Ｋredが求められる。触媒上からのＣＯの脱離速度Ｒredは触媒９の活性状態に応じて変化することから、活性状態と相関する触媒温度Ｔcatfに応じて脱離速度定数Ｋredが設定される。
吸着定数算出部５２では、触媒温度Ｔcatf及び空燃比センサ１０により検出された排気空燃比A/Fから予め設定されたマップに従って吸着速度定数Ｋadが求められる。触媒上へのＣＯの吸着速度Ｒadは触媒９の活性状態のみならず、排気空燃比A/Fに応じた排ガス組成の相違に影響されることから、触媒温度Ｔcatf及び排気空燃比A/Fに応じて吸着速度定数Ｋadが設定される。なお、排気空燃比A/Fは周知の演算手法によりエンジン運転状態などから推定してもよい。
【００２８】
吸着サイト算出部５３では、触媒温度Ｔcatfから予め設定されたマップに従って触媒９のＣＯ吸着のために機能する全ての吸着サイトρco（mol）が求められる。吸着サイトρcoは触媒固有の値であるものの、触媒９の活性状態に応じて変化することから、活性状態と相関する触媒温度Ｔcatfに応じて吸着サイトρcoが設定される。
ＣＯ濃度算出部５４では、排気空燃比A/Fから予め設定されたマップに従って排ガス中のＣＯ濃度Ｐco（atm）が求められる。ＣＯ濃度Ｐcoは排気空燃比A/Fと相関することから、排気空燃比A/Fに応じてＣＯ濃度Ｐcoが設定される。
【００２９】
一方、前回処理時の酸化反応熱量ΔＨrを算出する過程で触媒９のＣＯ吸着率θcoが算出されており、当該ＣＯ吸着率θcoは上記吸着サイトρcoとの間に次式(1)の関係が成立している。
θco＝rρco/ρco ………(1)
ここに、rρcoは実際の触媒９のＣＯ吸着サイトである。
【００３０】
上記ＣＯ吸着率θcoと、上記脱離定数算出部５１で算出された脱離速度定数Ｋred、及び吸着サイト算出部５３で算出された吸着サイトρcoとが脱離速度算出部５５に入力される。脱離速度算出部５５では、これらの入力値に基づき次式(2)に従って触媒上からのＣＯ脱離速度Ｒred（mol/sec）が算出される。
Ｒred＝Ｋred・ρco・θco ………(2)
また、上記ＣＯ吸着率θcoは減算部５６に入力され、減算部５６では１からＣＯ吸着率θcoが減算される（１−θco）。減算部５６による算出値と、上記吸着定数算出部５２で算出された吸着速度定数Ｋad、吸着サイト算出部５３で算出された吸着サイトρco、及びＣＯ濃度算出部５４で算出されたＣＯ濃度Ｐcoとが吸着速度算出部５７に入力され、吸着速度算出部５７では、これらの入力値に基づき次式(3)に従って触媒上へのＣＯ吸着速度Ｒad（mol/sec）が算出される。
【００３１】
Ｒad＝Ｋad・ρco・（１−θco）・Ｐco ………(3)
また、ＣＯ濃度算出部５４で算出されたＣＯ濃度Ｐco、全圧１（atm）、単位時間当たりの排ガス中の全物質量ｎ all（mol/sec）が分圧算出部５８に入力され、分圧算出部５８では次式(4)に従って分圧として排ガス中のＣＯ量GASco（mol/sec）が算出される。なお、詳細は説明しないが排ガス中の全物質量ｎ allは、吸入空気量Ｑ、空気及び燃料の分子量、排気空燃比A/Fなどに応じて周知の演算手法により算出される。
【００３２】
GASco＝Ｐco/１・ｎ all ………(4)
上記触媒上へのＣＯ吸着速度Ｒadと排ガス中のＣＯ量GAScoとは最小選択部５９に入力され、最小選択部５９では入力値の小の側を選択して最終的な触媒上へのＣＯ吸着速度Ｒadとして決定する。即ち、ＣＯ吸着速度Ｒadは触媒９のＣＯ吸着能力を意味する値であり、当該ＣＯ吸着速度Ｒad以上のＣＯ量GAScoを排ガスが有している場合には、ＣＯ吸着速度Ｒadが実際の触媒上へのＣＯ吸着速度Ｒadとなるが、排ガス中のＣＯ量GAScoがＣＯ吸着速度Ｒad未満の場合には、触媒上へのＣＯ吸着速度ＲadはＣＯ量GAScoに制限される。このため、最小選択部５９で何れか小の側を選択しているのである。
【００３３】
上記脱離速度算出部５５で算出されたＣＯ脱離速度Ｒred及び最小選択部５９で選択された触媒上へのＣＯ吸着速度Ｒadは単位時間吸着ＣＯ量算出部６０に入力され、単位時間吸着ＣＯ量算出部６０では次式(5)に従って単位時間当たりに触媒９に吸着する単位時間吸着ＣＯ量Δadθcoが算出される。
Δadθco＝Ｒad−Ｒred ………(5)
算出された単位時間吸着ＣＯ量Δadθcoは、触媒９の吸着サイトρco及びＥＣＵ２１の演算周期ｆ（例えば、0.1msec）と共に吸着率変化量算出部６１に入力され、吸着率変化量算出部６１では次式(6)に従って演算周期ｆ間のＣＯ吸着率θcoの変化量Δθcoが算出される。
【００３４】
Δθco＝Δadθco/ρco・ｆ ………(6)
ＣＯ吸着率θcoの変化量ΔθcoはＣＯ吸着率θcoと共に吸着率算出部６２に入力され、吸着率算出部６２ではＣＯ吸着率θcoを前回処理時の値θco(n-1)と見なし、次式(7)に従って今回のＣＯ吸着率θco(n)が算出される。
θco(n)＝θco(n-1)＋Δθco ………(7)
その後、ＣＯ吸着率θco(n)は吸着サイトρco及びＥＣＵ２１の演算周期fと共に吸着ＣＯ量算出部６３に入力され、吸着ＣＯ量算出部６３では次式(9)に従って現在触媒９に吸着されている吸着ＣＯ量の単位時間相当値CATco（mol/sec）が算出される。
【００３５】
CATco＝θco(n)・ρco / f ………(9)
一方、Ｏ２濃度算出部６４では、排気空燃比A/Fから予め設定されたマップに従ってＯ２濃度Ｐo２が求められる。上記ＣＯ濃度Ｐcoと同じくＯ２濃度Ｐo２は排気空燃比A/Fに応じて決定されるが、後述する反応速度式(10)に適用するときの便宜のために、このときマップからはＯ２濃度Ｐo２を1/2乗した値Ｐo２1/2が算出される。
【００３６】
反応速度定数算出部６５では、触媒温度Ｔcatfから予め設定されたマップに従ってＣＯとＯ２とが反応するときの反応速度定数Ｋrが求められる。ＣＯとＯ２との反応速度は触媒９の活性状態に応じて変化することから触媒温度Ｔcatfに応じて反応速度定数Ｋrが設定される。
Ｏ２濃度算出部６４で算出されたＯ２濃度Ｐo２1/2、反応速度定数算出部６５で算出された反応速度定数ＫrはＣＯ吸着率θco及び吸着サイトρcoと共に反応速度算出部６６に入力され、反応速度算出部６６では次式(10)に従ってＣＯとＯ２とが反応可能な最大の反応速度ｒ（mol/sec）が算出される。
【００３７】
ｒ＝Ｋr・ρco・θco（Ｐo２）1/2 ………(10)
即ち、触媒上でのＣＯとＯ２との最大の反応速度ｒはＣＯ量CATco（＝θco・ρco）とＯ２濃度Ｐo２との積により定まり、且つ次式(11)で示すようにＣＯとＯ２との反応モル比が「２」であることから、上式(10)ではＯ２濃度をＰo２1/2としているのである。
ＣＯ＋1/2Ｏ２→ＣＯ２ ………(11)
上式(11)に基づいて設定された反応モル比の「２」、Ｏ２濃度算出部６４で算出されたＯ２濃度Ｐo２1/2、全圧１（atm）、単位時間当たりの排ガス中の全物質量ｎ allは排気Ｏ２量算出部６７に入力され、排気Ｏ２量算出部６７では次式(12)に従って排ガス中のＯ２量GASo２（mol/sec）が算出される。
【００３８】
GASo２＝２・Ｐo２/１・ｎ all ………(12)
吸着ＣＯ量算出部６３で算出された触媒上の吸着ＣＯ量CATco、反応速度算出部６６で算出された反応速度ｒ、排気Ｏ２量算出部６７で算出された排気Ｏ２量GASo２は最小選択部６８に入力され、最小選択部６８ではこれらの値の最小値がＣＯの反応量Ｒct（mol/sec）として選択される。基本的にＣＯとＯ２との反応は、吸着ＣＯ量CATcoと排気Ｏ２量GASo２との何れか小の側の供給量により定まるが、反応速度ｒを越える反応は不可能であるため、反応速度ｒを加えた３つの要件の最小値がＣＯ反応量Ｒctとして決定される。
【００３９】
最小選択部６８で算出されたＣＯ反応量Ｒctはエンジン１の燃料カット中にセットされる燃料カットフラグと共にＦ／Ｃ判定部６９に入力され、Ｆ／Ｃ判定部６９では、燃料カットフラグのリセットにより燃料カット中でないと判定したときには０が選択され、燃料カットフラグのセットにより燃料カット中と判定したときにはＣＯ反応量Ｒctが選択される。Ｆ/Ｃ判定部６９で選択された値は燃料の反応熱量（例えば、283kJ/mol）及びＥＣＵ２１の演算周期ｆと共にΔＨr算出部７０に入力され、ΔＨr算出部７０ではこれらの値に基づき触媒上でのＣＯの酸化反応により発生するＣＯ酸化反応熱量ΔＨrが算出される。従って、Ｆ／Ｃ時反応熱量算出部５０からは燃料カット中でないときには０が出力され、燃料カット中にはＣＯ酸化反応熱量ΔＨrが出力される。
【００４０】
Ｆ/Ｃ判定部６９で選択された値（Ｒctまたは０）は吸着サイトρcoと共に吸着率減少分算出部７１に入力され、吸着率減少分算出部７１ではこれらの値に基づいてＣＯの反応により触媒上の吸着率の減少分Δrθcoが算出される。算出された吸着率減少分Δrθcoは上記吸着率算出部６２で算出された今回のＣＯ吸着率θco(n)と共に吸着率更新部７２に入力され、吸着率更新部７２では次式(13)に従って最新のＣＯ吸着率θcoが算出される。
【００４１】
θco＝θco(n)−Δrθco ………(13)
吸着率更新部７２で算出されたＣＯ吸着率θcoは上限クリップ部７３により上限を1.0に制限され、下限クリップ部７４により下限を０に制限される。各クリップ部７３，７４での処理は、何らかの要因により不適切なＣＯ吸着率θcoが算出された場合への対策であり、これらの処理後のＣＯ吸着率θcoが次回のＥＣＵ２１の処理で適用される。
【００４２】
次に、上記排ガスと触媒９との間の熱伝達量ΔＨt（kJ/sec）の算出処理について説明する。
図３はＥＣＵ２１の排ガス・触媒間の熱伝達量算出部８０の処理を示す制御ブロック図である。排ガスと触媒９との間の熱伝達量ΔＨtは、排ガスと触媒９との温度差、排ガスと触媒９との接触面積、及び熱伝達係数に依存し、これらの間には次式(14）が成立する。
【００４３】
ΔＨt＝h・Ｓv・（Ｔex−Ｔcatf）・Ｖcat ………(14)
ここに、Ｓvは触媒９の比表面積、ｈは熱伝達係数、Ｔexは推定処理により推定された排ガス温度、Ｔcatfは触媒温度、Ｖcatは触媒体積である。なお、ここではTexは、推定処理により推定された排ガス温度としているが、図１に示すように排気温度センサ１１を用いて検出するようにしてもよい。
【００４４】
図３において熱伝達量算出部８０の接触面積算出部８１では、予め設定されたマップに従って吸入空気量Ｑから触媒９の有効接触面積Ｓcatが算出される。有効接触面積Ｓcatは排ガスに対する触媒９の接触面積であり、上式(14)の触媒９の比表面積Ｓvと触媒体積Ｖcatとの積に相当する。
また、熱伝達係数算出部８２では、予め設定されたマップに従って吸入空気量Ｑ及び排ガス温度Ｔexから熱伝達係数ｈが算出される。温度差算出部８３では、排ガス温度Ｔexと仮触媒温度Ｔcat（触媒温度Ｔcatfのフィルタ処理前の値）との差ΔＴが算出される。
【００４５】
上記有効接触面積Ｓcat、熱伝達係数ｈ、温度差ΔＴはＥＣＵ２１の演算周期ｆと共にΔＨt算出部８４に入力され、ΔＨt算出部８４では各値が乗算され、結果として上式(14)に従って演算周期ｆ間の熱伝達量ΔＨt（kJ/sec）が算出される（熱伝達量算出手段）。
次に、上記燃料復帰時の触媒９に吸蔵されたＯ２によるＣＯ浄化反応熱量ΔＨcの算出処理について説明する。
【００４６】
図４はＥＣＵ２１のＦ／Ｃ復帰時浄化熱量算出部１１０の処理を示す制御ブロック図、図５はＥＣＵ２１のＯ２吸蔵率算出部９０の処理を示す制御ブロック図である。触媒上でのＣＯ浄化反応熱量ΔＨcはエンジン１の運転状態、即ち、エンジン回転速度Ｎｅ、充填効率Ｅｃ、触媒上の空燃比A/Fに依存し、これらの値から決定することができる。但し、触媒上の空燃比Ａ／Ｆは基本的に排気空燃比A/Fを適用できるが、触媒９によるＯ２ストレージ機能の影響を受けることからこの点を考慮する必要がある。即ち、触媒９は排気空燃比A/Fがリーン側のときにＯ２を吸蔵する一方、排気空燃比A/FがリッチのときにＯ２を放出することで、排気空燃比A/Fをストイキオ（A/F=14.7）に保持する作用を奏するため、このようなＯ２ストレージ機能が発揮されているときには、排気空燃比A/Fがストイキオ以外であっても触媒上の空燃比はストイキオに保持される一方、Ｏ２ストレージ機能が発揮されていないときには、触媒上の空燃比は排気空燃比A/Fとなる。
【００４７】
そこで、当該ＣＯ浄化反応熱量ΔＨcの算出処理では、触媒９のＯ２ストレージ機能と相関するＯ２吸蔵率θo２を算出し、算出したＯ２吸蔵率θo２をＣＯ浄化反応熱量ΔＨc（kJ/sec）の算出処理に反映させており、まず、Ｏ２吸蔵率θo２の算出処理について説明する。
上記燃料カット時の触媒上でのＣＯの挙動と同様に、燃料復帰時には排ガス中のＣＯとの反応により触媒上からＯ２が放出される現象のみならず、排ガス中のＯ２が触媒上に吸蔵される現象も並行して発生しており、これらの現象が相俟って触媒９のＯ２吸蔵率θo２、ひいては燃料復帰時のＣＯ浄化反応熱量ΔＨcに影響を及ぼす。そこで、当該ＣＯ浄化反応熱量ΔＨcの算出処理では、触媒上からのＯ２の放出及び触媒上へのＯ２の吸蔵を考慮して触媒９の吸蔵Ｏ２量を算出し、この吸蔵Ｏ２量をＣＯ浄化反応熱量ΔＨcの算出処理に反映させている。
【００４８】
図５においてＯ２吸蔵率算出部９０の吸蔵サイト算出部９１では、触媒温度Ｔcatfから予め設定されたマップに従って触媒９のＯ２吸蔵のために機能する吸蔵サイトρo２（mol）が求められる。吸蔵サイトρo２は触媒固有の値であるものの、上記ＣＯ吸着サイトと同じく触媒９の活性状態に応じて変化することから、活性状態と相関する触媒温度Ｔcatfに応じて吸蔵サイトρo２が設定される。吸蔵サイトρo２はＯ２吸蔵率θo２との間に次式(15)の関係が成立している。
【００４９】
θo２＝rρo２/ρo２ ………(15)
ここに、rρo２は実際の触媒９のＯ２吸蔵サイトである。
吸蔵サイト算出部９１で算出された吸蔵サイトρo２及び前回処理時のＯ２吸蔵率θo２はＯ２放出能力量算出部９２に入力され、Ｏ２放出能力量算出部９２では、これらの入力値に基づき次式(16)に従って現在の触媒９のＯ２放出能力量Ｃred(mol)が算出される（触媒能力算出手段）。
【００５０】
Ｃred＝ρo２・θo２・(−１) ………(16)
なお、上式(15)で−１の乗算によりＯ２放出能力量Ｃredを負の値として設定しているのは、Ｏ２放出がＯ２吸蔵率θo２に対して減少側に作用するためである。
また、減算部９３では１からＯ２吸蔵率θo２が減算され（１−θo２）、減算後の値と上記吸蔵サイト算出部９１で算出された吸蔵サイトρo２とがＯ２吸蔵能力量算出部９４に入力され、Ｏ２吸蔵能力量算出部９４では、これらの入力値に基づき次式(17)に従って現在の触媒９のＯ２吸蔵能力量Ｃad(mol)が算出される（触媒能力算出手段）。
【００５１】
Ｃad＝ρo２・（１−θo２） ………(17)
一方、Ｏ２過不足量算出部９５では、排気空燃比A/Fから予め設定されたマップに従って排ガス中のＯ２過不足量ΔＯ２(vol%)が求められる。当該Ｏ２過不足量ΔＯ２はストイキオを基準としたＯ２の過剰量またはＯ２の不足量を意味し、排気空燃比A/Fに応じて設定される。設定されたＯ２過不足量ΔＯ２は単位時間当たりの排ガス中の全物質量ｎ all（mol/sec）と共に換算部９６に入力され、単位時間当たりの値（mol/sec）に換算される。
【００５２】
Ｏ２消費量算出部９７では上記Ｆ／Ｃ時反応熱量算出部５０のＦ/Ｃ判定部６９から出力されたＣＯ反応量Ｒct（mol/sec）、及びＣＯとＯ２との反応モル比の「２」が入力され、これらの入力値に基づき次式(18)に従って燃料カット中に触媒上のＣＯと反応して消費された排ガス中のＯ２消費量expΔＯ２が算出される（Ｏ２消費量算出手段）。
expΔＯ２＝Ｒct/２ ………(18)
減算部９８には、換算部９６で換算後のＯ２過不足量ΔＯ２、及びＯ２消費量算出部９７で算出されたＯ２消費量expΔＯ２が入力され、次式(19)に従って、排ガスをストイキオにするために必要な（即ち、排ガスを浄化するために必要な）単位時間当たりのＯ２量としてＯ２ストイキオ必要量Δado２が算出される（Ｏ２必要量算出手段）。具体的には、排気空燃比A/Fがリッチ側で排ガス中のＣＯのＯ２による酸化反応を要するときには触媒上からの放出Ｏ２量（負の値）として、また、排気空燃比A/Fがリーン側で余剰Ｏ２の吸蔵を要するときには触媒上への吸蔵Ｏ２量（正の値）としてＯ２ストイキオ必要量Δadθo２が算出される。算出されたＯ２ストイキオ必要量Δado２は換算部９９に入力され、次式(20)に従って演算周期ｆ間の値Δadθo２に換算される。
【００５３】
Δado２＝ΔＯ２−expΔＯ２ ………(19)
Δadθo２＝Δado２・ｆ ………(20)
換算後のＯ２ストイキオ必要量Δadθo２は、上記触媒９のＯ２放出能力量Ｃred及びＯ２吸蔵能力量Ｃadと排気空燃比A/Fと共にＯ２放出吸蔵量算出部１００に入力される。Ｏ２放出吸蔵量算出部１００では、排気空燃比A/Fがストイキオ以下のリッチ側で、触媒上においてＯ２の放出現象が生起されていると推測されるときには、Ｏ２放出能力量ＣredとＯ２ストイキオ必要量Δadθo２との大の側が選択される（共に負の値のため絶対値では小の側となる）。また、排気空燃比A/Fがストイキオを越えてリーン側で、触媒上においてＯ２の吸蔵現象が生起されていると推測されるときには、Ｏ２吸蔵能力量ＣadとＯ２ストイキオ必要量Δadθo２との小の側が選択される。
【００５４】
即ち、排ガスのストイキオのために触媒能力（Ｃred，Ｃad）以上のＯ２量（Δadθo２）が要求されても、実際に使用されるＯ２量は触媒能力相当分に制限され、逆にストイキオのために要求されるＯ２量（Δadθo２）を越える触媒能力（Ｃred，Ｃad）が確保されていても、全触媒能力の内のＯ２量相当分だけが使用されることから、このような点を鑑みてＯ２放出吸蔵量算出部１００での処理が設定されている。
【００５５】
Ｏ２放出吸蔵量算出部１００で選択された値（ＣredまたはＣredまたはΔadθo２）は吸蔵サイトρo２と共に吸蔵率減少分算出部１０１に入力され、吸蔵率減少分算出部１０１ではこれらの値に基づいてＯ２の反応により触媒上の吸蔵率の減少分Δrθo２が算出される。算出された吸蔵率の減少分Δrθo２は前回処理時のＯ２吸蔵率θo２共に吸蔵率更新部１０２に入力され、吸蔵率更新部１０２では次式(21)に従って最新のＯ２吸蔵率θo２が算出される（Ｏ２ストレージ指標算出手段）。
【００５６】
θo２＝θo２＋Δrθo２ ………(21)
吸蔵率更新部１０２で算出されたＯ２吸蔵率θo２は上限クリップ部１０３により上限を1.0に制限され、下限クリップ部１０４により下限を０に制限される。
以上のように算出された触媒９のＯ２吸蔵率θo２が図４に示すＦ／Ｃ復帰時浄化熱量算出部１１０のＡ／Ｆ判定部１１１に入力される。Ａ／Ｆ判定部１１１では、Ｏ２吸蔵率θoが０（Ｏ２の放出限界）または1.0（Ｏ２の吸蔵限界）であり触媒９がＯ２ストレージ機能を奏していないと判定したときには、触媒上の空燃比として排気空燃比A/Fが選択される（空燃比推定手段）。また、Ｏ２吸蔵率θoが０＞θ0＞1.0の条件を満たして触媒９がＯ2ストレージ機能を奏していると判定したときには、触媒上の空燃比としてストイキオが選択される（空燃比推定手段）。
【００５７】
Ａ／Ｆ判定部１１１で判定された触媒上の空燃比は回転速度センサ２３により検出されたエンジン回転速度Ｎｅ（運転状態検出手段）及びエンジン１の運転状態に基づき算出された充填効率Ｅｃ（運転状態検出手段）と共にΔＨc算出部１１２に入力され、ΔＨc算出部１１２ではエンジン回転速度Ｎｅ、充填効率Ｅｃ及び触媒上の空燃比Ａ／Ｆに基づき予め設定されたマップに従って単位時間当たりのＣＯ浄化反応熱量ΔＨc（kJ/sec）が算出される（浄化反応熱量算出手段）。
【００５８】
算出されたＣＯ浄化反応熱量ΔＨcは燃料カットフラグと共にＦ／Ｃ判定部１１３に入力され、Ｆ／Ｃ判定部１１３では、燃料カットフラグのリセットにより燃料カット中でないと判定したときにはＣＯ浄化反応熱量ΔＨcが選択され、燃料カットフラグのセットにより燃料カット中と判定したときには０が選択される。Ｆ/Ｃ判定部１１３で選択された値はＥＣＵ２１の演算周期ｆと共に周期換算部１１４に入力され、周期換算部１１４ではこれらの値に基づき演算周期ｆ間のＣＯ浄化反応熱量ΔＨc（kJ）が算出される。従って、Ｆ／Ｃ復帰時浄化熱量算出部１１０からは燃料カット中でないときにはＣＯ浄化反応熱量ΔＨcが出力され、燃料カット中には０が出力される。
【００５９】
次に、最終的な触媒温度Ｔcatfの推定処理について説明する。
図６はＥＣＵ２１が実行する触媒温度推定部１２０の処理状況を示す制御ブロック図である。Ｆ／Ｃ時反応熱量算出部５０で算出された燃料カット時のＣＯ酸化反応熱量ΔＨr、上記Ｆ／Ｃ復帰時反応熱量算出部１１０で算出された燃料復帰時のＣＯ浄化反応熱量ΔＨc、熱伝達量算出部８０で算出された排ガス・触媒間の熱伝達量ΔＨtが総熱量算出部１２１に入力される。上記のように燃料カット時にはＣＯ酸化反応熱量ΔＨrと排ガス・触媒間の熱伝達量ΔＨtとが有効値として算出され、燃料復帰時にはＣＯ浄化反応熱量ΔＨcと排ガス・触媒間の熱伝達量ΔＨtとが有効値として算出されるため、総熱量算出部１２１ではそれぞれの場合に対応する２つの値が加算されて触媒温度上昇の要因となる総熱量ΣＨが算出される。
【００６０】
算出された総熱量ΣＨは温度上昇量算出部１２２に入力され、総熱量ΣＨと予め設定された触媒９の熱容量とに基づいて触媒９の温度上昇量ΔＴが算出されて触媒温度算出部１２３に入力される。触媒温度算出部１２３では前回処理時の仮触媒温度Ｔcat（以下に述べるフィルタ処理前の触媒温度）に温度上昇量ΔＴが加算されて仮触媒温度Ｔcatが算出される。この仮触媒温度Ｔcatが前回値として記憶されると共に、この仮触媒温度Ｔcatをベースとして触媒中の熱伝導を模擬した３次フィルタにより最終的な触媒温度Ｔcatfが算出される（触媒温度推定手段）。
【００６１】
３次フィルタのゲイン算出部１２４では予め設定されたマップにより吸入空気量ＱからフィルタゲインＫcaが設定される。図に示すように吸入空気量Ｑの増加に伴って次第にフィルタゲインＫcaも増加し、設定されたフィルタゲインＫcaが上記仮触媒温度Ｔcatと共に第１フィルタ部１２５に入力される。第１フィルタ部１２５では上記仮触媒温度Ｔcat、フィルタゲインＫca、前回処理時の第１触媒温度セルＴf1（第1フィルタ処理後の触媒温度）に基づき、今回の第１触媒温度セルＴf1が算出され、この第１触媒温度セルＴf1が記憶されると共に第２フィルタ部１２６に入力される。
【００６２】
第２フィルタ部１２６では上記第１触媒温度セルＴf1、フィルタゲインＫca、前回処理時の第２触媒温度セルＴf2（第２フィルタ処理後の触媒温度）に基づき、今回の第２触媒温度セルＴf2が算出され、この第２触媒温度セルＴf2が記憶されると共に第３フィルタ部１２７に入力される。
第３フィルタ部１２７では上記第２触媒温度セルＴf2、フィルタゲインＫca、前回処理時の第３触媒温度セルＴf3（第３フィルタ処理後の触媒温度）に基づき、今回の第３触媒温度セルＴf3が算出され、この第３触媒温度セルＴf3が記憶されると共に最終的な触媒温度Ｔcatfとして出力される。吸入空気量Ｑと共にフィルタゲインＫcaが増加されるほど、各フィルタ部１２５〜１２７の処理では前回値が反映され難くなって今回値に近い値が算出され、結果として最終的な触媒温度Ｔcatfは吸入空気量Ｑと相関する触媒中での熱伝導を模擬した値として算出される。
【００６３】
以上のように本実施形態の触媒温度推定装置では、燃料復帰時において触媒上に吸蔵されたＯ２により排ガス中のＣＯが浄化されたときのＣＯ浄化反応熱量ΔＨcを算出し、このＣＯ浄化反応熱量ΔＨcと熱伝達量ΔＨtとに基づいて床下触媒９の温度を推定している。従って、燃料復帰時の触媒温度の挙動を正確に推定でき、この触媒温度Ｔcatfに基づいて適切にエンジン制御を実行することにより、耐熱温度を越えた温度上昇による触媒の破損を未然に防止することができる。
【００６４】
図７は燃料カットと燃料復帰が短い周期で繰り返されたときの触媒温度の推定状況を従来技術と比較した試験結果を示し、このような運転状態は運転者が頻繁にアクセルをオン・オフした場合が考えられる。燃料カット時に発生するＣＯ酸化反応熱量ΔＨr及び燃料復帰時に発生するＣＯ浄化反応熱量ΔＨcにより実触媒温度は排ガス温度Ｔexに比較して急激に上昇し、従来技術による推定値ではこの温度上昇がほとんど反映されない。これに対して本実施形態の推定手法によれば、実触媒温度に極めて近似する推定値Ｔcatfが得られており、この試験結果からも上記本発明の作用効果が裏付けられる。
【００６５】
しかも、ＣＯ浄化反応熱量ΔＨcの算出に適用する触媒上の空燃比として、Ｏ２吸蔵率θoが０または1.0のときには（Ｏ２ストレージ機能無し）排気空燃比A/Fを選択し、Ｏ２吸蔵率θoが０＞θ0＞1.0のときには(Ｏ２ストレージ機能有り)ストイキオを選択するため、実状に則した正確な触媒上の空燃比に基づいてＣＯ浄化反応熱量ΔＨc、ひいては触媒温度Ｔcatfの挙動を一層正確に推定することができる。
【００６６】
また、排ガスをストイキオにするために要するＯ２ストイキオ必要量Δadθo２、即ち、排ガス浄化に際して触媒上で放出及び吸蔵されるＯ２量に基づいてＯ２吸蔵率θo２を算出し、しかも、Ｏ２ストイキオ必要量Δadθo２と触媒のＯ２放出能力量Ｃred及びＯ２吸蔵能力量Ｃadとを大小関係に応じて選択し、選択した値に基づいてＯ２吸蔵率θo２を算出しているため、正確なＯ２吸蔵率θo２、ひいては正確なＣＯ浄化反応熱量ΔＨcを算出することができる。
【００６７】
以上で実施形態の説明を終えるが、本発明の態様はこの実施形態に限定されるものではない。例えば、上記実施形態では燃料カット時と燃料復帰時との双方で共に触媒温度Ｔcatfを推定したが、必ずしも燃料カット時に触媒温度Ｔcatfを推定する必要はなく、例えば燃料復帰時のみにＣＯ浄化反応熱量ΔＨcrに基づいて触媒温度Ｔcatfを推定するようにしてもよい。また、上記実施形態ではエンジン１は筒内噴射型火花点火式のガソリンエンジンとして説明したが、これに限定されることはなく、吸気ポートに燃料を供給するガソリンエンジンであってもよく、またディーゼルエンジンでも本発明を適用可能である。
【図面の簡単な説明】
【００６８】
【図１】実施形態の触媒温度推定装置を示す全体構成図である。
【図２】ＥＣＵのＦ／Ｃ時反応熱量算出部の処理を示す制御ブロック図である。
【図３】ＥＣＵの排ガス・触媒間の熱伝達量算出部の処理を示す制御ブロック図である。
【図４】ＥＣＵのＦ／Ｃ復帰時浄化熱量算出部の処理を示す制御ブロック図である。
【図５】ＥＣＵのＯ２吸蔵率算出部の処理を示す制御ブロック図である。
【図６】ＥＣＵの触媒温度推定部の処理状況を示す制御ブロック図である。
【図７】燃料カットと燃料復帰が短い周期で繰り返されたときの触媒温度の推定状況の試験結果を示す図である。
【符号の説明】
【００６９】
１      エンジン
９      床下触媒
１０    空燃比センサ（排気空燃比検出手段）
２１    ＥＣＵ
２３    回転速度センサ（運転状態検出手段）
８４    ΔＨt算出部（熱伝達量算出手段）
９４    Ｏ２吸蔵能力量算出部（触媒能力算出手段）
９２    Ｏ２放出能力量算出部（触媒能力算出手段）
９７ Ｏ２消費量算出部（Ｏ２消費量算出手段）
９８    減算部（Ｏ２必要量算出手段）
１０２  吸蔵率更新部（Ｏ２ストレージ指標算出手段）
１１１  Ａ／Ｆ判定部（空燃比推定手段）
１１２  ΔＨc算出部（浄化反応熱量算出手段）
１２３  触媒温度算出部（触媒温度推定手段）
【図１】
0004535287.tif 000002
【図２】
0004535287.tif 000003
【図３】
0004535287.tif 000004
【図４】
0004535287.tif 000005
【図５】
0004535287.tif 000006
【図６】
0004535287.tif 000007
【図７】
0004535287.tif 000008