(57)【特許請求の範囲】
【請求項１】
内燃機関の運転状態に基づいて、触媒へ流入するガスの特性を推定する触媒入りガス推定手段と、
予め構築された触媒反応の総括式を用いて、前記触媒へ流入するガスの特性に基づいて、触媒から排出されるガスの特性を推定する触媒出ガス推定手段と、
触媒に与えられた熱負荷に基づいて、前記総括式を補正する補正手段と、を備え、
前記総括式は、触媒内で行われる反応の反応速度定数を表す数式であり、前記数式中には、触媒中の活性点の数を表す頻度因子の項と、反応の際のエネルギーを表す活性化エネルギーの項と、を含み、
前記補正手段は、
前記触媒入りガス推定手段で推定された前記触媒へ流入するガスの温度及び酸素濃度に基づいて、触媒に与えられた熱負荷の総量に応じて変化する前記総括式の補正項を算出する手段と、
触媒へ流入するガスの空燃比に基づいて、前記総括式の補正項を可変する手段と、
前記補正項を用いて前記総括式を補正する手段と、
を含むことを特徴とする排気ガスの推定装置。
【請求項２】
内燃機関の運転状態に基づいて、触媒へ流入するガスの特性を推定する触媒入りガス推定手段と、
予め構築された触媒反応の総括式を用いて、前記触媒へ流入するガスの特性に基づいて、触媒から排出されるガスの特性を推定する触媒出ガス推定手段と、
触媒に与えられた熱負荷に基づいて、前記総括式を補正する補正手段と、を備え、
前記総括式は、触媒内で行われる反応の反応速度定数を表す数式であり、前記数式中には、触媒中の活性点の数を表す頻度因子の項と、反応の際のエネルギーを表す活性化エネルギーの項と、を含み、
前記補正手段は、
前記触媒入りガス推定手段で推定された前記触媒へ流入するガスの温度及び酸素濃度に基づいて、触媒に与えられた熱負荷の総量に応じて変化する前記総括式の補正項を算出する手段と、
触媒へ流入するガスの空燃比及び温度に基づいて、前記総括式の補正項を可変する手段と、
前記補正項を用いて前記総括式を補正する手段と、
を含むことを特徴とする排気ガスの推定装置。
【請求項３】
前記触媒出ガス推定手段は、前記総括式によって触媒におけるガスの浄化率を規定し、当該浄化率に基づいて前記触媒から排出されるガスの特性を推定することを特徴とする請求項１または２記載の排出ガスの推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、排気ガスの推定装置に関する。
【背景技術】
【０００２】
従来、例えば特開２００５−１０５９８１号公報には、内燃機関の触媒内の浄化反応を表したモデルに、実際に検出、取得される特定ガスの触媒に流入する際の濃度、特定ガスが触媒から流出する際の濃度、触媒に流入し、特定ガスと反応する際の反応ガスの濃度、およびその触媒の温度を利用することにより、頻度因子を算出することが記載されている。
【０００３】
【特許文献１】特開２００５−１０５９８１号公報
【特許文献２】特開２００４−８９０８号公報
【特許文献３】特開平９−２０９７４０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、実際の触媒は、高温の排気ガスが通過するため、触媒の温度、および排気ガス中の酸素濃度によって劣化する。上記従来の技術では、このような熱負荷による触媒劣化を考慮していないため、触媒から排出されたガスのエミッション値を精度良く推定することは困難である。
【０００５】
この発明は、上述のような問題を解決するためになされたものであり、モデル計算による排出ガス成分の推定精度を向上することを目的とする。
【課題を解決するための手段】
【０００６】
第１の発明は、上記の目的を達成するため、内燃機関の運転状態に基づいて、触媒へ流入するガスの特性を推定する触媒入りガス推定手段と、予め構築された触媒反応の総括式を用いて、前記触媒へ流入するガスの特性に基づいて、触媒から排出されるガスの特性を推定する触媒出ガス推定手段と、触媒に与えられた熱負荷に基づいて、前記総括式を補正する補正手段と、を備え、前記総括式は、触媒内で行われる反応の反応速度定数を表す数式であり、前記数式中には、触媒中の活性点の数を表す頻度因子の項と、反応の際のエネルギーを表す活性化エネルギーの項と、を含み、前記補正手段は、前記触媒入りガス推定手段で推定された前記触媒へ流入するガスの温度及び酸素濃度に基づいて、触媒に与えられた熱負荷の総量に応じて変化する前記総括式の補正項を算出する手段と、触媒へ流入するガスの空燃比に基づいて、前記総括式の補正項を可変する手段と、前記補正項を用いて前記総括式を補正する手段と、を含むことを特徴とする。
【０００７】
第２の発明は、上記の目的を達成するため、内燃機関の運転状態に基づいて、触媒へ流入するガスの特性を推定する触媒入りガス推定手段と、予め構築された触媒反応の総括式を用いて、前記触媒へ流入するガスの特性に基づいて、触媒から排出されるガスの特性を推定する触媒出ガス推定手段と、触媒に与えられた熱負荷に基づいて、前記総括式を補正する補正手段と、を備え、前記総括式は、触媒内で行われる反応の反応速度定数を表す数式であり、前記数式中には、触媒中の活性点の数を表す頻度因子の項と、反応の際のエネルギーを表す活性化エネルギーの項と、を含み、前記補正手段は、前記触媒入りガス推定手段で推定された前記触媒へ流入するガスの温度及び酸素濃度に基づいて、触媒に与えられた熱負荷の総量に応じて変化する前記総括式の補正項を算出する手段と、触媒へ流入するガスの空燃比及び温度に基づいて、前記総括式の補正項を可変する手段と、前記補正項を用いて前記総括式を補正する手段と、を含むことを特徴とする。
【０００８】
第３の発明は、第１又は第２の発明において、前記触媒出ガス推定手段は、前記総括式によって触媒におけるガスの浄化率を規定し、当該浄化率に基づいて前記触媒から排出されるガスの特性を推定することを特徴とする。
【発明の効果】
【００１２】
第１の発明によれば、予め構築された触媒反応の総括式を用いて、触媒から流出するガスの特性を推定するモデルにおいて、触媒に与えられた熱負荷に基づいて、総括式を補正するため、触媒から排出されるガスの特性を高精度に推定することが可能となる。従って、触媒から排出されるガスのエミッション値を精度良く求めることが可能となる。
また、本発明によれば、熱負荷により触媒劣化が進行すると、触媒中の活性点の数を表す頻度因子の値と、反応の際のエネルギーを表す活性化エネルギーの値の少なくとも一方が変化するため、頻度因子の項と、活性化エネルギーの項とを含む総括式を用いて排出ガスの特性を推定するモデルにおいて、熱負荷に応じて総括式を補正することで、排出ガスの特性を精度良く推定することが可能となる。
また、本発明によれば、触媒へ流入するガスの温度及び酸素濃度から熱負荷が求まるため、触媒へ流入するガスの温度及び酸素濃度に基づいて補正項を算出し、総括式を補正することで、熱負荷による総括式の補正を行うことができる。
また、本発明によれば、触媒における浄化率は、触媒へ流入するガスの空燃比に応じて変化するため、触媒へ流入するガスの空燃比に基づいて総括式の補正項を可変することで、触媒から排出されるガスの特性を精度良く推定することができる。
【００１３】
第２の発明によれば、予め構築された触媒反応の総括式を用いて、触媒から流出するガスの特性を推定するモデルにおいて、触媒に与えられた熱負荷に基づいて、総括式を補正するため、触媒から排出されるガスの特性を高精度に推定することが可能となる。従って、触媒から排出されるガスのエミッション値を精度良く求めることが可能となる。
また、本発明によれば、熱負荷により触媒劣化が進行すると、触媒中の活性点の数を表す頻度因子の値と、反応の際のエネルギーを表す活性化エネルギーの値の少なくとも一方が変化するため、頻度因子の項と、活性化エネルギーの項とを含む総括式を用いて排出ガスの特性を推定するモデルにおいて、熱負荷に応じて総括式を補正することで、排出ガスの特性を精度良く推定することが可能となる。
また、本発明によれば、触媒へ流入するガスの温度及び酸素濃度から熱負荷が求まるため、触媒へ流入するガスの温度及び酸素濃度に基づいて補正項を算出し、総括式を補正することで、熱負荷による総括式の補正を行うことができる。
また、本発明によれば、触媒における浄化率は、触媒へ流入するガスの空燃比及び温度に応じて変化するため、触媒へ流入するガスの空燃比に基づいて総括式の補正項を可変することで、触媒から排出されるガスの特性を精度良く推定することができる。
【００１４】
第３の発明によれば、総括式によって触媒におけるガスの浄化率を規定するため、浄化率に基づいて触媒から排出されるガスの特性を推定することができる。
【発明を実施するための最良の形態】
【００１８】
以下、図面に基づいてこの発明の一実施形態について説明する。尚、各図において共通する要素には、同一の符号を付して重複する説明を省略する。なお、以下の実施の形態によりこの発明が限定されるものではない。
【００１９】
図１は、本発明の一実施形態に係る触媒モデル１０及びその周辺の構成を示す模式図である。図１の構成は、実機の内燃機関システムが排気通路に備える触媒と、その周辺の構成をモデル化したものである。実機の内燃機関システムにおいては、吸入空気量、燃料噴射量などが車載のＥＣＵ（Electronic Control Unit）によって制御される。図１に示すモデルは、例えばＥＣＵ内に構築されるものである。
【００２０】
図１に示すように、触媒モデル１０には、触媒入りガスモデル１２、触媒劣化モデル１４が接続されている。触媒入りガスモデル１２は、内燃機関の機関回転数Ne、負荷率KL、空燃比A/F、点火時期SA、バルブタイミングVTなどの運転条件の入力値に基づいて、触媒に入るガス（触媒入りガス）の成分、温度、流量などを算出するモデルである。
【００２１】
触媒モデル１０は、触媒入りガスモデル１２で算出された触媒入りガスの成分、温度、流量等を入力として、触媒から排出されるガス（触媒出ガス）の成分、温度を出力するモデルである。触媒モデル１０は、実機の触媒で行われる反応を、三元反応と酸素吸脱の総括反応でモデル化したものである。触媒モデル１０によれば、触媒出ガスの成分、温度に基づいて、排気ガスのエミッション値が明らかになる。従って、排気ガスの成分を実測することなく、エミッション値をモデル計算により求めることが可能である。
【００２２】
図１に示すように、触媒モデル１０には、触媒形状の物性値が入力される。触媒形状の物性値としては、例えば、触媒のセルの開口面積、セル数、熱伝達率などがある。触媒モデル１０は、これらの特性値を考慮した上で、触媒出ガスの成分、温度を出力する。
【００２３】
実機システムの触媒内では、一酸化炭素ＣＯと一酸化窒素ＮＯとの反応、一酸化炭素ＣＯと酸素Ｏ２との反応、炭化水素ＨＣと酸素Ｏ２との反応など、複数の反応が行われる。触媒モデル１０では、これらの全ての反応がモデル化される。そして、触媒モデル１０は、これらの反応のそれぞれについて、触媒が有するガスの浄化率に基づいて、触媒入りガスの成分、温度、流量から触媒出ガスの成分、温度を算出する。
【００２４】
以下、一酸化炭素ＣＯと酸素Ｏ２との反応を例に挙げて説明する。触媒内において、一酸化炭素ＣＯは酸素Ｏ２と反応して二酸化炭素ＣＯ２になる。図２は、この反応による一酸化炭素ＣＯの浄化率（触媒の反応効率）を示す特性図である。図２に示すように、一酸化炭素ＣＯの浄化率は、触媒へ流入する排気ガスの空燃比A/Fに応じて変化する。浄化率は、以下の（１）式で表すことができる。
浄化率＝（［ＣＯ］入−［ＣＯ］出）／［ＣＯ］入        ・・・（１）
上式において、［ＣＯ］入は、触媒入りガス中における一酸化炭素ＣＯの濃度である。また、［ＣＯ］出は、触媒出ガス中における一酸化炭素ＣＯの濃度である。
【００２５】
図２に示すように、空燃比A/Fがリッチの場合は、排気ガスに含まれる酸素量が少ないため、一酸化炭素ＣＯの浄化率は低下する。空燃比A/Fがリーンになると、排気ガス中に含まれる酸素量が増加するため、一酸化炭素ＣＯの浄化率は増加する。
【００２６】
触媒モデル１０では、一酸化炭素ＣＯと酸素Ｏ２との反応において、反応速度定数Ｋを同定することによって、予め図２の特性がモデル化されている。従って、触媒モデル１０によれば、図２の特性を用いることで、触媒入りガスモデル１２から入力された触媒入りガス中の一酸化炭素濃度に基づいて、触媒出ガス中の一酸化炭素濃度を算出することができる。
【００２７】
以下、触媒モデル１０において、図２の特性をモデル化する手法を説明する。一酸化炭素ＣＯと酸素Ｏ２が反応する際の速度は、反応速度Ｒによって表される。そして、反応速度Ｒは、反応速度定数Ｋによって決定される。ここで、反応速度Ｒと反応速度定数Ｋとの間には、以下の（２）式の関係が成立する。
Ｒ＝Ｋ＊［ＣＯ］入＊［Ｏ２］入／Ｇ        ・・・（２）
【００２８】
また、反応速度定数Ｋは、通常、以下の（３）式より算出される。
Ｋ＝ａ＊ｅｘｐ（−Ｅ／ＲＴ）        ・・・（３）
これらの式において、ａは頻度因子、Ｇは抑制項、Ｅは活性化エネルギーを示している。また、［Ｏ２］入は、触媒入りガス中における酸素Ｏ２の濃度である。また、Ｒは気体定数、Ｔはガスの温度である。ここで、頻度因子ａは、触媒における活性ポイントの数を表すパラメータである。頻度因子ａの値が増加するほど活性ポイントの数が増加するため、触媒における反応が促進される。従って、頻度因子ａの値が大きくなるほど、反応速度定数Ｋの値は大きくなる。また、活性化エネルギーＥは、反応が生じる際に必要なエネルギーを示している。活性化エネルギーＥの値が小さくなるほど、反応が生じ易くなるため、反応速度定数Ｋの値は大きくなる。
【００２９】
（３）式によれば、頻度因子ａ、抑制項Ｇ、活性化エネルギーＥの値を適正値に設定することで、図２の特性を同定することができる。触媒モデル１０は、（３）式を総括式として触媒反応をモデル化することで、触媒入りガス中の成分等に基づいて、触媒出ガス中のエミッション値を出力する。
【００３０】
図２の特性をモデル化する際には、実機の触媒から得られる浄化率の実測値に基づいて、反応速度定数Ｋを同定する処理が行われる。以下、図３及び図４に基づいて、反応速度定数Ｋを同定する手法を説明する。図３は、空燃比と浄化率との関係を実測する方法を示す模式図である。図３（Ａ）は、触媒入りガス中および触媒出ガス中の一酸化炭素の濃度の実測値を示している。図３（Ａ）において、破線で示す特性は触媒入りガス中の一酸化炭素濃度を示しており、実線で示す特性は触媒出ガス中の一酸化炭素濃度を示している。また、図３（Ｂ）は、実機の触媒へ送られる触媒入りガスの空燃比を示している。
【００３１】
図３（Ｂ）に示すように、反応速度定数Ｋを同定する際には、触媒入りガスの空燃比がリッチ−リーンの間で変化するように、内燃機関へ供給される混合気の空燃比が制御される。図３（Ｂ）の例では、時間Ｔを周期として、空燃比が１４と１５の間で変化するように制御が行われる。なお、図３（Ｂ）に示す時刻ｔ０以前には、触媒入りガスの空燃比はリーンに設定されている。
【００３２】
図３（Ａ）に示す一酸化炭素濃度は、触媒入りガスの空燃比を図３（Ｂ）に示す特性で可変させた場合に実測されたものである。図３（Ｂ）中に破線で示すように、時刻ｔ０で空燃比が１４（リッチ）に設定されると、触媒入りガス中の一酸化炭素濃度は比較的高い値となり、時刻ｔ１で空燃比が１５（リーン）に設定されると、触媒入りガス中の一酸化炭素濃度は低下する。
【００３３】
図３（Ａ）に実線で示すように、時刻ｔ０で空燃比が１４（リッチ）に設定された直後は、時刻ｔ０以前の空燃比がリーンであることから、触媒出ガス中の一酸化炭素濃度は低下している。その後、触媒出ガス中の一酸化炭素濃度は、時間の経過とともに上昇し、触媒入りガス中の一酸化炭素濃度と同程度の値となる。
【００３４】
時刻ｔ１で空燃比が１５に設定されると、触媒入りガス中の一酸化炭素濃度が低下し、触媒入りガス中に含まれる一酸化炭素は触媒を通過する際に殆ど浄化される。従って、空燃比が１５に設定されている時刻ｔ１から時刻ｔ２の間では、触媒出ガス中の一酸化炭素濃度はほぼ０となる。
【００３５】
図３のように、空燃比がリッチの場合とリーンの場合のそれぞれについて、触媒入りガス中の一酸化炭素濃度と、触媒出ガス中の一酸化炭素濃度が求まると、（１）式に基づいて、空燃比がリッチの場合とリーンの場合のそれぞれについて、一酸化炭素の浄化率の実測値を求めることができる。
【００３６】
反応速度定数ｋの同定は、空燃比がリッチの場合における浄化率と、空燃比がリーンの場合の浄化率とに基づいて行われる。この際、一酸化炭素の濃度は、時刻ｔ１の直前の時刻ｔ３、および時刻ｔ２の直前の時刻ｔ４で安定状態となっているため、時刻ｔ３および時刻ｔ４のタイミングで実測された浄化率に基づいて、反応速度定数ｋの同定を行うことで、酸素の吸脱反応による影響を抑えることができ、精度良く同定を行うことができる。
【００３７】
上述のように、時刻ｔ３及び時刻ｔ４における浄化率は、（１）式から算出することができる。また、時刻ｔ３における触媒入りガスの空燃比は１４であり、時刻ｔ４における触媒出ガスの空燃比は１５である。従って、空燃比が１４と１５の場合における浄化率を算出することができる。
【００３８】
図４は、図３の手法により算出された空燃比と浄化率の関係に基づいて、反応速度定数Ｋを同定する手法を示す模式図である。図４に示すように、空燃比が１４の場合の浄化率と、空燃比が１５の場合の浄化率は図４中に黒丸でプロットされる。そして、図３の手法による浄化率の算出は、触媒入りガスの空燃比を例えば１４．２と１４．８の間で変化させた場合にも行われる。これにより、図４に示すように、空燃比が１４．２の場合の浄化率と、空燃比が１４．８の場合の浄化率は図４中に白丸でプロットされる。従って、触媒入りガスの空燃比を可変して、このような処理を繰り返すことにより、図２で示すにように、浄化率と空燃比との関係を取得することができる。
【００３９】
そして、このようにして求められた浄化率の特性に適合するように、反応速度定数Ｋの同定が行われる。具体的には、図４の浄化率の特性に適合するように、反応速度乗数Ｋを決定するパラメータのうち、頻度因子ａ、抑制項Ｇ、活性化エネルギーＥの値が決定される。同定された反応速度乗数Ｋによれば、図２の特性を再現することができる。これにより、触媒モデル１０によれば、触媒入りガスの一酸化炭素濃度に基づいて、触媒出ガスの一酸化炭素濃度を算出することが可能となる。
【００４０】
次に、触媒劣化モデル１４について説明する。実機の触媒は、排気ガスによって高温状態となり、触媒温度、酸素濃度（空燃比）などの熱負荷の要因によって劣化していく。触媒劣化モデル１４は、触媒モデル１０から触媒温度、酸素濃度などの熱負荷の入力を受け、熱負荷に応じて触媒反応速度定数Ｋ、最大酸素吸蔵量などの特性値が変化する様子をモデル化する。
【００４１】
そして、触媒劣化モデル１４から触媒モデル１０には、これらの特性値の変化量を表す入力値が贈られる。従って、触媒モデル１０は、反応速度定数Ｋを同定することで得られた上述した浄化率のモデルと、触媒劣化モデル１４から入力された反応速度定数Ｋの変化量を表すパラメータとに基づいて、触媒劣化を考慮した上で触媒出ガスのエミッションを正確に出力することができる。
【００４２】
以下、触媒劣化モデル１４で行われる処理について説明する。熱負荷によって触媒が劣化すると、（３）式で表される触媒反応速度定数の右辺の項において、主として、頻度因子ａ、活性化エネルギーＥの項が熱負荷に応じて変化する。図６は、頻度因子ａ、活性化エネルギーＥが熱負荷に応じて変化する様子の一例を示す特性図である。図６において、横軸は触媒に与えられた熱負荷の積算値を示している。図６に示すように、頻度因子ａ、活性化エネルギーＥは、触媒に与えられた熱負荷の総量が多くなるほど変化する。
【００４３】
図７は、図３及び図４に示す方法で取得された浄化率の特性が、熱負荷に応じて変化した様子を示す特性図である。ここで、図７中に示す破線の特性は、図３及び図４の手法で反応速度定数Ｋを同定することによって得られた浄化率の特性を示しており、触媒に熱負荷による劣化が生じていない状態で取得した特性を示している。また、図７中に実線で示す特性は、熱負荷の総量が図６中に示す所定値Ｓに到達した時点で、実機の触媒から得られた実測値の浄化率を示している。
【００４４】
図７に示すように、触媒は熱負荷の総量に応じて劣化するため、熱負荷の総量が図６中に示す所定値Ｓに到達した時点では、浄化率の実測値と、モデル化により得られた浄化率との間に乖離が生じる。従って、熱負荷を考慮しない場合は、触媒モデル１０により算出された触媒出ガスのエミッション値と実際のエミッション値との間に誤差が生じる。
【００４５】
図７に示すように、浄化率の実測値と、モデル化により得られた浄化率との間の乖離量は、触媒入りガスの空燃比によって変化する。このため、本実施形態では、（３）式の右辺に、空燃比A/Fに応じて変化する補正項f（Ａ／Ｆ）を追加することとしている。以下の（３）’式は、（３）式の右辺に補正項ｆ（Ａ／Ｆ）を乗算して得られた式である。
Ｋ＝ｆ（Ａ／Ｆ）＊ａ＊ｅｘｐ（−Ｅ／ＲＴ）        ・・・（３）’
【００４６】
補正項ｆ（Ａ／Ｆ）は、空燃比Ａ／Ｆの関数であり、また、触媒に与えられた熱負荷に応じて変化する係数である。（３）’式によれば、（３）式の右辺に空燃比Ａ／Ｆに応じた補正係数ｆ（Ａ／Ｆ）を乗算しているため、図７中に矢印で示すように、熱負荷の総量に応じて、破線で示されるモデル算出値を実測値に補正することが可能となる。
【００４７】
具体的には、図３及び図４の手法により、リッチ−リーンの間で反応速度定数Ｋを同定する際に、ストイキでの補正係数ｆ（Ａ／Ｆ）は１として、その他のＡ／Ｆで補正係数ｆ（Ａ／Ｆ）の値を変化させて、（３）式の右辺に補正係数ｆ（Ａ／Ｆ）を乗算するようにする。または、リッチとリーンで反応速度定数Ｋを同定する際に、実測値との乖離が少ない空燃比側での補正係数ｆ（Ａ／Ｆ）は１として、その他のＡ／Ｆで補正係数ｆ（Ａ／Ｆ）の値を変化させて、（３）式の右辺に補正係数ｆ（Ａ／Ｆ）を乗算するようにする。
【００４８】
なお、補正方法としては、補正係数ｆ（Ａ／Ｆ）のように関数の値を乗算しても良いし、マップから算出した補正係数を用いて補正を行うなど、マップにより反応速度定数Ｋを補正しても良い。マップの場合は、補間などの手法を用いることで、算出値の精度を高めることが好適である。
【００４９】
図１に示すように、触媒劣化モデル１４は、触媒モデル１０から触媒温度、酸素濃度などの値を取得し、触媒に与えられた熱負荷の総量を算出する。そして、触媒劣化モデル１４は、熱負荷の総量に応じて、補正係数ｆ（Ａ／Ｆ）を算出し、算出した補正係数ｆ（Ａ／Ｆ）を触媒モデル１０へ入力する。触媒モデル１０では、（３）式の右辺に補正係数ｆ（Ａ／Ｆ）を乗算し、（３）’式を構築する。このような手法によれば、触媒モデル１０において、熱負荷による触媒劣化を考慮した上で、触媒出ガスの成分を高精度に算出することが可能となる。
【００５０】
なお、（３）’式では、（３）式の右辺に補正係数ｆ（Ａ／Ｆ）を乗算しているが、空燃比Ａ／Ｆをストイキの空燃比１４．７で除算して得られる空気過剰率λ（＝Ａ／Ｆ／１４．７）を変数とし、補正係数ｆ’（λ）を（３）式の右辺に乗算しても良い。
【００５１】
また、図７で説明した、浄化率の実測値と、モデル化による浄化率との間の乖離量は、触媒入りガスの温度に応じて変化する場合がある。このため、補正係数ｆ（Ａ／Ｆ）を触媒入りガスの温度Ｔに応じて変化させるようにしても良い。すなわち、この場合は、補正係数ｆ（Ａ／Ｆ，Ｔ）を（３）式の右辺に乗算して得られた以下の（３）”式により、触媒反応をモデル化する。
Ｋ＝ｆ（Ａ／Ｆ，Ｔ）＊ａ＊ｅｘｐ（−Ｅ／ＲＴ）        ・・・（３）”
これにより、モデル算出値をより正確に補正することが可能となり、エミッション値の算出精度を向上することが可能となる。
【００５２】
上述したように、触媒内では、一酸化炭素ＣＯと酸素Ｏ２との反応以外にも複数の反応が行われる。触媒モデル１０では、これらの全ての反応についても、一酸化炭素ＣＯと酸素Ｏ２との反応と同様に、反応速度定数Ｋが同定される。そして、触媒劣化モデル１４では、これらの全ての反応のそれぞれについて、熱負荷に応じて変化する補正係数ｆ（Ａ／Ｆ）を算出する。従って、モデル化された複数の反応の全てにおいて、熱負荷による触媒劣化を考慮した上で、触媒出ガスの成分を高精度に算出することが可能である。
【００５３】
以上説明したように本実施形態によれば、熱負荷による触媒劣化を考慮した上で、触媒モデル１０から触媒出ガスの成分、温度などを算出することが可能となり、触媒出ガスのエミッションを高精度に算出することが可能となる。従って、算出した触媒出ガスのエミッションに基づいて、内燃機関システムを高精度に制御することが可能となる。
【図面の簡単な説明】
【００５４】
【図１】本発明の一実施形態に係る触媒モデル１０及びその周辺の構成を示す模式図である。
【図２】一酸化炭素ＣＯと酸素Ｏ２との反応において、触媒による一酸化炭素ＣＯの浄化率（触媒の反応効率）を示す特性図である。
【図３】図３は、空燃比と浄化率との関係を実測する方法を示す模式図である。
【図４】図３の手法により算出された空燃比と浄化率の関係に基づいて、反応速度定数Ｋを同定する手法を示す模式図である。
【図５】触媒入りガス及び触媒出ガスの濃度（浄化特性）が、温度及びＡ／Ｆに応じて変化する様子を示す特性図である。
【図６】頻度因子ａ、活性化エネルギーＥが熱負荷に応じて変化する様子の一例を示す特性図である。
【図７】図３及び図４に示す方法で取得された浄化率の特性が、熱負荷に応じて変化した様子を示す特性図である。
【符号の説明】
【００５５】
１０  触媒モデル
１２  触媒入りガスモデル
１４  触媒劣化モデル
【図１】
0004816164.tif 000002
【図２】
0004816164.tif 000003
【図３】
0004816164.tif 000004
【図４】
0004816164.tif 000005
【図５】
0004816164.tif 000006
【図６】
0004816164.tif 000007
【図７】
0004816164.tif 000008