(57)【要約】【課題】排出されるＮＯｘの状態を高精度に推定する。【解決手段】エンジンシステム（１０）において、ＥＣＵ（１００）の内部には、エンジン（２００）の三元触媒（２２３）から排出される触媒排出ガスの状態を推定するエミッション推定モデル（３００）が構築されている。エミッション推定モデル（３００）における触媒モデル（３１０）は、触媒流入モデル（３２０）及び触媒劣化モデル（３３０）からの指標値の入力を受け、予め実際の浄化対象物質の浄化率に近似するように係数が同定されたモデル式に従って浄化対象物質の各々について浄化特性を推定し、最終的にエミッション値を出力する。この際、ＮＯｘについては、触媒内の酸素濃度の影響が考慮され、当該モデル式に代えて、触媒内の酸素濃度がある程度低下した後にＮＯｘの還元反応が開始されるように設定されたモデル式が採用される。【選択図】図２
2008069708.tif 000002
【特許請求の範囲】
【請求項１】
触媒を備えた内燃機関において該触媒から排出される排出ガスの状態を推定する排出ガスの状態推定装置であって、
前記触媒へ流入する流入ガスの状態を特定する特定手段と、
予め前記触媒における浄化対象物質の各々の浄化特性を規定するものとして該各々について構築された、該各々の濃度に対応する項と、前記触媒における触媒反応に係る反応速度定数との積を含んでなる基準モデル式に従い、前記特定された流入ガスの状態に基づいて前記排出ガスの状態を推定する推定手段と、
前記浄化対象物質が窒素酸化物である場合に、酸素濃度が減少するのに応じて前記窒素酸化物の浄化率が上昇するように前記基準モデル式を補正する第1の補正手段と
を具備することを特徴とする排出ガスの状態推定装置。
【請求項２】
前記触媒に付与される熱負荷に基づいて前記基準モデル式を補正する第２の補正手段を更に具備する
ことを特徴とする請求項１に記載の排出ガスの状態推定装置。
【請求項３】
前記反応速度定数は、前記触媒中の活性点の数を表す頻度因子に対応する項と、前記触媒反応が生じる際のエネルギを表す活性化エネルギに対応する項とを含む
ことを特徴とする請求項１又は２に記載の排出ガスの状態推定装置。
【請求項４】
前記基準モデル式は、前記積と前記浄化対象物質に対応する還元剤の濃度との積を、該還元剤の種類毎に加算してなる加算項を含む
ことを特徴とする請求項１から３のいずれか一項に記載の排出ガスの状態推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、触媒を備えた内燃機関において当該触媒から排出される排出ガスの状態を推定する排出ガスの状態推定装置の技術分野に関する。
【背景技術】
【０００２】
この種の技術分野において、触媒反応に係る反応速度の解析手法が提案されている（例えば、特許文献１参照）。特許文献１に開示された、流体中の所定成分の濃度変化予測方法及び反応速度解析法（以下、「従来の技術」と称する）によれば、排ガスの成分濃度が触媒との接触によって変化する反応モデルを設定し、当該反応モデルに基づいて、所定成分の濃度変化量を表す速度式及び所定成分の濃度変化を表す位置・濃度特性を求めることにより、所定成分の濃度変化が予測される。このため、触媒担体の長さや所定成分の初期濃度など条件を変化させたときの所定成分の濃度変化に与える影響を見積もることが可能であるとされている。
【０００３】
【特許文献１】特開２００４−８９０８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
触媒において、ＮＯｘ（Nitrogen Oxide：窒素酸化物）には還元反応が生じ、酸素が脱離するが、触媒内に酸素が共存する場合、酸素濃度が相対的に高い領域では、係る還元反応の進行が遅れ易い。ところが、従来の技術では、このような酸素濃度が還元反応に与える影響が考慮されていないため、触媒内の酸素濃度が高い場合には特に、ＮＯｘの濃度変化の予測精度が低下し易い。即ち、従来の技術には、内燃機関の運転条件によっては、特にＮＯｘの濃度変化が十分に推定され難いという技術的な問題点がある。
【０００５】
本発明は上述した問題点に鑑みてなされたものであり、排出されるＮＯｘの状態を高精度に推定し得る排出ガスの状態推定装置を提供することを課題とする。
【課題を解決するための手段】
【０００６】
上述した課題を解決するため、本発明に係る排出ガスの状態推定装置は、触媒を備えた内燃機関において該触媒から排出される排出ガスの状態を推定する排出ガスの状態推定装置であって、前記触媒へ流入する流入ガスの状態を特定する特定手段と、予め前記触媒における浄化対象物質の各々の浄化特性を規定するものとして該各々について構築された、該各々の濃度に対応する項と、前記触媒における触媒反応に係る反応速度定数との積を含んでなる基準モデル式に従い、前記特定された流入ガスの状態に基づいて前記排出ガスの状態を推定する推定手段と、前記浄化対象物質が窒素酸化物である場合に、酸素濃度が減少するのに応じて前記窒素酸化物の浄化率が上昇するように前記基準モデル式を補正する第1の補正手段とを具備することを特徴とする。
【０００７】
本発明に係る「内燃機関」とは、例えば複数の気筒を有し、当該複数の気筒の各々における燃焼室において燃料が燃焼した際に発生する爆発力を、例えばピストン及びコネクティングロッド等の機械的な伝達経路を経て、例えばクランク軸等の入出力軸を介して動力として取り出すことが可能な機関を包括する概念であり、例えば２サイクル或いは４サイクルレシプロエンジン等を指す。
【０００８】
本発明に係る排出ガスの状態推定装置によれば、その動作時には、例えばＥＣＵ（Electronic Control Unit：電子制御ユニット）等の各種処理ユニット、各種コントローラ或いはマイコン装置等各種コンピュータシステム等として構成される特定手段により、触媒へ流入する流入ガスの状態が特定される。
【０００９】
ここで、「流入ガスの状態」とは、後述する排出ガスの状態の推定に供し得る指標値により定量的に表された状態を包括する概念であり、例えば、流入ガスにおける浄化対象物質の成分、成分比率、成分濃度及び成分量等の少なくとも一部、並びに流入ガスの流量及び温度等を含む趣旨である。
【００１０】
尚、本発明における「特定」とは、例えば、何らかの検出手段を介して直接的に又は間接的に物理的数値又は物理的数値に対応する電気信号等として検出すること、何らかの検出手段を介して直接的に又は間接的に例えば電気信号等の形で検出された、特定対象と対応関係を有する物理的数値に基づいて予め然るべき記憶手段等に記憶されたマップ等から該当する数値を選択すること、このような物理的数値又は選択された数値等から、予め設定されたアルゴリズムや計算式に従って導出すること、或いはこのように検出、選択又は導出された数値等を、例えば電気信号等の形で単に取得すること等を包括する広い概念である。
【００１１】
例えば特定手段は、例えば内燃機関の機関回転数、負荷率、空燃比、点火時期及びバルブタイミング等、上述した流入ガスの状態との直接的又は間接的な対応関係が例えば予め実験的に、経験的に、理論的に或いはシミュレーション等に基づいて実践上不利益が顕在化しない程度の信頼性を伴って確定されている内燃機関の運転状態を規定する各種指標値に基づいて、例えば予め設定されたアルゴリズムや計算式等に従って流入ガスの状態を然るべき指標値として算出すること、或いは当該各種指標値に基づいて然るべき値をマップ等から選択すること等により流入ガスの状態を特定する。
【００１２】
このように特定手段により流入ガスの状態が特定されると、例えばＥＣＵ等の各種処理ユニット、各種コントローラ或いはマイコン装置等各種コンピュータシステム等として構成される推定手段により、予め構築された基準モデル式に従い、特定手段によって特定された流入ガスの状態に基づいて、排出ガスの状態が推定される。
【００１３】
ここで、「基準モデル式」とは、予め触媒における浄化対象物質の各々の浄化特性を規定すべく、当該各々について構築された式であり、浄化対象物質の濃度に対応する項と、触媒における触媒反応に係る反応速度定数との積を含む式である。この基準モデル式は、推定された流入ガスの状態の関数であり、推定手段は、この基準モデル式に例えば特定された流入ガスの状態（即ち、形態としては各種指標値）を代入して得られる、例えば、浄化対象物質各々についての浄化率、浄化量或いは浄化速度（触媒反応の反応速度）等に基づいて、排出ガスの状態を推定する。
【００１４】
ここで、「排出ガスの状態」とは、内燃機関におけるエミッションの種類や量を含む概念としてのエミッション特性の評価に供し得る指標値によって表された状態を包括する概念であり、例えば上述した流入ガスに係る状態と、意味合いとして少なくとも一部が同等であってもよく、例えば、排出ガスにおける浄化対象物質の成分、成分比率、成分量及び成分濃度、並びに排出ガスの温度等を含む趣旨である。
【００１５】
尚、「反応速度定数」とは、例えば予め実験的に、経験的に、理論的に或いはシミュレーション等に基づいて、実際の触媒における浄化特性を近似的に或いは代替的に表し得るように例えば各種係数又は特性項の値等が決定（同定）されてなる定数であり、浄化対象物質の濃度に乗じられることに鑑みれば、実質的には浄化対象物質の浄化特性（例えば浄化量、浄化率若しくは浄化速度又は触媒における反応速度等）を決定付ける定数である。
【００１６】
尚、「浄化対象物質の濃度に対応する項」とは、当該濃度を含み、且つ基準モデル式を実際の触媒における浄化特性に近似させ得る或いは代替させ得る項を包括する概念であり、当該濃度のみに限定されず、例えば、当該濃度のべき乗項等であってもよい趣旨である。
【００１７】
ここで特に、浄化対象物質がＮＯ（Nitric Oxide：一酸化窒素）等を含む窒素酸化物（以下、ＮＯｘ）である場合、触媒内の酸素濃度が相対的に高い領域では、酸素の脱離が阻害され、還元反応の進行が遅れ易い。従って、上述した基準モデル式によってＮＯｘの浄化特性を表そうとした場合、触媒内の酸素濃度が相対的に高くなり易い、例えば内燃機関の高負荷領域（高流量域）等において、排出ガスの状態に係る推定精度が低下し易い。
【００１８】
そこで、本発明に係る排出ガスの状態推定装置では、例えばＥＣＵ等の各種処理ユニット、各種コントローラ或いはマイコン装置等各種コンピュータシステム等として構成される第１の補正手段の作用により、浄化対象物がＮＯｘである場合に基準モデル式が補正される。より具体的には、第１の補正手段は、酸素濃度が減少するのに応じてＮＯｘの浄化率が上昇するように基準モデル式を補正する。
【００１９】
尚、「浄化率が上昇するように」とは、浄化率が上昇することが判明している、或いは浄化率が上昇すると推定される限りにおいて、基準モデル式自体がＮＯｘの浄化率を直接規定する式であるか否かとは無関係であってよい趣旨である。また、第１の補正手段に係る補正の態様は、ＮＯｘの浄化率を酸素濃度の減少に応じて上昇せしめ、基準モデル式により規定されるＮＯｘの浄化特性を、何らこの種の補正がなされない場合と比較して幾らかなりとも実際のＮＯｘの浄化特性に近似させ得る得る限りにおいて限定されない趣旨である。例えば基準モデル式に対し、何らかの補正項を加算、減算、乗算又は除算すること等により係る補正がなされてもよい。また、「酸素濃度が減少するのに応じて」とは、必ずしも酸素濃度の変化が浄化率の変化と一対一に対応しておらずともよい趣旨であり、浄化率の変化と酸素濃度の変化とは、相互に一対多の関係であってもよい。
【００２０】
本発明に係る排出ガスの状態推定装置の一の態様では、前記触媒に付与される熱負荷に基づいて前記基準モデル式を補正する第２の補正手段を更に具備する。
【００２１】
この態様によれば、例えばＥＣＵ等の各種処理ユニット、各種コントローラ或いはマイコン装置等各種コンピュータシステム等として構成される第２の補正手段の作用により、基準モデル式が、触媒に付与される熱負荷に基づいて補正される。
【００２２】
ここで、触媒に付与される熱負荷とは、好適には触媒内に流入する流入ガスの温度、触媒の温度、流入ガスの流量、流入ガスにおける酸素濃度、或いは空燃比等の関数であり、第２の補正手段は、例えば、これらを表す指標値（好適には、特定手段により特定される、流入ガスの状態の一部である）に応じて可変な補正係数、補正式或いは補正関数を、予め然るべき記憶手段に記憶されたマップから選択する、或いは予め設定されたアルゴリズムや算出式に従って算出又は導出すること等を介して取得し、基準モデル式を補正する。
【００２３】
触媒における浄化対象物質の浄化特性は、触媒に付与される熱負荷に応じて相応に劣化するが、この態様によれば、熱負荷に応じて、基準モデル式、即ち触媒における浄化特性が補正されるため、最終的に推定手段によって推定される排出ガスの状態が、係る熱負荷を考慮したものとなり、一層正確に排出ガスの状態を推定することが可能となる。
【００２４】
本発明に係る排出ガスの状態推定装置の他の態様では、前記反応速度定数は、前記触媒中の活性点の数を表す頻度因子に対応する項と、前記触媒反応が生じる際のエネルギを表す活性化エネルギに対応する項とを含む。
【００２５】
この態様によれば、基準モデル式を構成する反応速度定数に、活性点の数を表す頻度因子及び触媒反応が生じる際に必要となるエネルギを表す活性化エネルギの項が含まれるため、触媒における浄化対象物質の浄化特性を、実際の浄化特性により近付けることが可能となり、高精度に排出ガスの状態を推定することが可能となる。
【００２６】
尚、この態様によれば、前述した第２の補正手段に係る補正に際し、熱負荷に応じて頻度因子及び活性化エネルギ各々の項を補正すること等により、触媒の劣化が好適に近似され得る。
【００２７】
本発明に係る排出ガスの状態推定装置の他の態様では、前記基準モデル式は、前記積と前記浄化対象物質に対応する還元剤の濃度との積を、該還元剤の種類毎に加算してなる加算項を含む。
【００２８】
この態様によれば、基準モデル式が、前述した反応速度定数と浄化対象物質の濃度との積に更に還元剤（例えば、浄化対象物質がＮＯｘであればＣＯ（一酸化炭素）やＨＣ（炭化水素）等）の濃度を乗じた項を、還元剤の種類毎に加算してなる加算項を含む。即ち、この態様においては、基準モデル式により、浄化対象物質の濃度のみならず還元剤の影響をも考慮した上での浄化対象物質の浄化特性が規定される。従って、より精細に排出ガスの状態を推定することが可能となる。
【００２９】
尚、還元剤の影響が考慮される場合、上述した加算項を構成する、浄化対象物質の濃度に対応する項と還元剤の濃度に対応する項と反応速度定数との積の各々について、例えば反応速度定数に係る係数や特性項（即ち、例えば、上述した頻度因子や活性化エネルギ或いは更に抑制項等）の同定が必要となり、同定に係る時間的負荷が高くなり易い。従って、還元剤の影響を考慮するか否かは、例えば、車両の仕様、仕向け、内燃機関の要求性能等に応じて、その都度個別具体的に決定されてよい。
【００３０】
本発明のこのような作用及び他の利得は次に説明する実施形態から明らかにされる。
【発明を実施するための最良の形態】
【００３１】
＜発明の実施形態＞
以下、図面を参照して、本発明の好適な各種実施形態について説明する。
【００３２】
＜１：第１実施形態＞
＜１−１：実施形態の構成＞
＜１−１−１：エンジンシステムの構成＞
始めに、図１を参照して、本発明の第１実施形態に係るエンジンシステム１０の構成について説明する。ここに、図１は、エンジンシステム１０の模式図である。
【００３３】
図１において、エンジンシステム１０は、ＥＣＵ１００及びエンジン２００を備える。
【００３４】
ＥＣＵ１００は、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）及びＲＡＭ（Random Access Memory）等を備え、エンジンシステム１０の動作全体を制御する電子制御ユニットである。また、ＥＣＵ１００内部には、エンジン２００のエミッション状態（即ち、本発明に係る「排出ガスの状態」の一例）を推定するための、後述するエミッション推定モデル３００が構築されており、本実施形態において、ＥＣＵ１００は、本発明に係る「排出ガスの状態推定装置」の一例としても機能するように構成されている。
【００３５】
エンジン２００は、本発明に係る「内燃機関」の一例である。エンジン２００は、気筒２０１内にその一部たる点火プラグの一部が露出してなる点火装置２０２の点火動作により混合気を爆発させると共に、爆発力に応じて生じるピストン２０３の往復運動を、コネクションロッド２０４を介してクランクシャフト２０５の回転運動に変換することが可能に構成されている。また、クランクシャフト２０５近傍には、クランクシャフト２０５の回転位置（即ち、クランク角）を検出するクランクポジションセンサ２０６が設置されている。クランクポジションセンサ２０６は、ＥＣＵ１００と電気的に接続されており、ＥＣＵ１００は、クランクポジションセンサ２０６によって検出されたクランク角に基づいて、点火装置２０２の点火時期等を制御することが可能に構成されている。また、ＥＣＵ１００は、クランクシャフト２０５の回転位置に基づいてエンジン２００の機関回転数Ｎｅを算出することが可能に構成されている。以下に、エンジン２００の要部構成を、その動作の一部と共に説明する。
【００３６】
気筒２０１内における燃料の燃焼に際し、外部から吸入された空気は吸気管２０７を通過し、吸気ポート２１３において、インジェクタ２１４から噴射された燃料と混合されて前述の混合気となる。燃料は、燃料タンク２１５に貯留されており、低圧ポンプ２１７の作用により、デリバリパイプ２１６を介してインジェクタ２１４に圧送供給されている。インジェクタ２１４は、ＥＣＵ１００と電気的に接続されており、この供給される燃料を、ＥＣＵ１００の制御に従って吸気ポート２１３に噴射することが可能に構成されている。
【００３７】
気筒２０１内部と吸気管２０７とは、吸気バルブ２１８の開閉によって連通状態が制御されている。気筒２０１内部で燃焼した混合気は排気となり吸気バルブ２１８の開閉に連動して開閉する排気バルブ２１９の開弁時に排気ポート２２０を介して排気管２２１に導かれる。
【００３８】
吸気管２０７上には、クリーナ２０８が配設されており、外部から吸入される空気が浄化される構成となっている。また、クリーナ２０８の下流側（気筒側）には更に、エアフローメータ２０９が配設されている。エアフローメータ２０９は、ホットワイヤー式と称される形態を有しており、吸入された空気の質量流量を直接検出することが可能に構成されている。尚、エアフローメータ２０９は、ＥＣＵ１００と電気的に接続されており、検出された吸入空気の質量流量は、ＥＣＵ１００によって絶えず把握される構成となっている。
【００３９】
吸気管２０７におけるエアフローメータ２０９の下流側には、気筒２０１内部へ吸入される空気に係る吸入空気量を調節するスロットルバルブ２１０が配設されている。このスロットルバルブ２１０には、スロットルポジションセンサ２１２が電気的に接続されており、その開度であるスロットル開度を検出することが可能に構成されている。
【００４０】
スロットバルブモータ２１１は、ＥＣＵ１００と電気的に接続され、スロットルバルブ２１０を駆動することが可能に構成されたモータである。ＥＣＵ１００は、不図示のアクセルポジションセンサによって検出されるアクセル開度に基づいてスロットルバルブモータ２１１の駆動状態を制御することが可能に構成されており、これによりスロットルバルブ２１０の開閉状態（即ち、スロットル開度）が制御される構成となっている。
【００４１】
尚、スロットルバルブ２１０は、上述したように一種の電子制御式スロットルバルブであり、スロットル開度は、ＥＣＵ１００により運転者の意思（即ち、アクセル開度）とは無関係に制御され得る。
【００４２】
排気管２２１には、三元触媒２２３が設置されている。三元触媒２２３は、エンジン２００から排出される排出ガス中のＣＯ及びＨＣを酸化し、且つＮＯｘを還元することにより夫々浄化することが可能に構成された触媒であり、本発明に係る「触媒」の一例である。
【００４３】
排気管２２１における三元触媒２２３の上流側には、空燃比センサ２２２が配設されている。空燃比センサ２２２は、排気ポート２２０を介して排出される排出ガスからエンジン２００の空燃比を検出することが可能に構成されている。空燃比センサ２２２は、ＥＣＵ１００と電気的に接続されており、検出された空燃比は、絶えずＥＣＵ１００によって把握される構成となっている。
【００４４】
また、気筒２０１を収容するシリンダブロックに設置されたウォータージャケットには、エンジン２００を冷却するための冷却水の温度を検出するための温度センサ２２４が配設されている。温度センサ２２４は、ＥＣＵ１００と電気的に接続されており、検出された冷却水の温度は、ＥＣＵ１００によって絶えず把握される構成となっている。
【００４５】
＜１−１−２：エミッション推定モデルの構成＞
次に、図２を参照し、ＥＣＵ１００内部に構築されたエミッション推定モデル３００の構成について説明する。ここに、図２は、エミッション推定モデル３００の構成を概念的に表してなる模式図である。
【００４６】
図２において、エミッション推定モデル３００は、触媒モデル３１０、触媒流入ガスモデル３２０及び触媒劣化モデル３３０を含んで構成される。
【００４７】
触媒モデル３１０は、三元触媒２２３から排出される触媒排出ガスのエミッション状態（例えば、触媒排出ガスの成分、成分比率、成分濃度又は成分量等）を表すエミッション値を推定結果として出力することが可能に構成されたモデルである。触媒モデル３１０は、実機の触媒（例えば三元触媒２２３或いはそれと同等な構成を有する触媒等）で生じる触媒反応をモデル化したものであり、本発明に係る「推定手段」、「第１の補正手段」及び「第２の補正手段」の一例として機能するように構成されている。
【００４８】
触媒流入ガスモデル３２０は、エンジン２００における機関回転数Ｎｅ、負荷率ＫＬ、空燃比Ａ／Ｆ、点火時期ＳＡ及びバルブタイミングＶＴ等の指標値（図示「エンジン状態指標値」に相当）を入力値として受け取り、予め設定されたアルゴリズム及び算出式に従って三元触媒２２３に流入する触媒流入ガスの状態（例えば、触媒流入ガスの成分、成分比率、成分濃度又は成分量、触媒流入ガスの温度及び流量等）を表す指標値を算出することが可能に構成されている。触媒流入ガスモデル３２０は、触媒モデル３１０と電気的に接続されており、算出された触媒流入ガスの状態を表す指標値は、触媒モデル３１０に入力される。触媒流入ガスモデル３２０は、本発明に係る「特定手段」の一例である。
【００４９】
触媒劣化モデル３３０は、触媒モデル３１０と電気的に接続され、熱負荷によって生じる、三元触媒２２３における浄化特性の劣化をモデル化したものである。触媒劣化モデル３３０は、触媒モデル３１０から触媒温度及び酸素濃度等を含む熱負荷指標値を入力として受け取り、受け取った熱負荷指標値に基づいて、熱負荷に応じた補正値を算出して触媒モデル３１０に出力することが可能に構成されている。尚、触媒劣化モデル３３０は、触媒モデル３１０と共に本発明に係る「第２の補正手段」の一例をなしている。
【００５０】
＜１−２：実施形態の動作＞
以下、本実施形態の動作として、上記エミッション推定モデル３００の動作について説明する。尚、三元触媒２２３では、主としてＮＯｘの還元反応（ＮＯｘからの酸素の脱離反応）、ＣＯの酸化反応、及びＨＣの酸化反応等の各種触媒反応が発生するが、以下の説明では、ＮＯｘの還元反応を例に挙げて説明することとする。
【００５１】
＜１−２−１：エミッション状態推定の基本動作＞
触媒モデル３１０では、基本的に下記（１）式（本発明に係る「基準モデル式」の一例）により、浄化対象物質の各々について三元触媒２２３における浄化率が規定されている。尚、（１）式における［ＮＯｘ］とは、浄化対象物質の一つたるＮＯｘの濃度を表し、初期値となる値が触媒流入ガスモデル３２０より、前述した触媒流入ガスの状態を規定する指標値の一部として取得されている。また、ｋは基本的に「１」に設定され、浄化率の同定の過程で適宜調整される、べき乗指数である。
【００５２】
ｄ［ＮＯｘ］／ｄｔ＝Ｋ・［ＮＯｘ］＾ｋ・・・（１）
ここで、上記（１）式におけるＫは反応速度定数であり、下記（２）式により規定される。
【００５３】
Ｋ＝Ａ／Ｇ・ｅｘｐ（―Ｅ／ＲＴ）・・・（２）
ここで、Ａは頻度因子、Ｇは抑制項、Ｅは活性化エネルギを表す。また、Ｒは公知の気体定数であり、Ｔは触媒流入ガスの絶対温度である。ここで、頻度因子Ａは、三元触媒２２３における活性点の数を表すパラメータであり、値が大きい程、活性点の数が増えることを表す。従って、頻度因子Ａの値が大きい程、三元触媒２２３における触媒反応が促進されている状態が表される。
【００５４】
また、活性化エネルギＥは、触媒反応が生じる際に必要となるエネルギを表しており、小さい程、触媒反応が生じ易い状態を表す。即ち、活性化エネルギＥが小さい程、反応速度定数Ｋの値は大きくなる。また、抑制項Ｇは、触媒反応を抑制する要因を表し、大きい程触媒反応が抑制された状態を表す。即ち反応速度定数Ｋの値が小さくなる。
【００５５】
上記（１）式に従って、三元触媒２２３における浄化率が、浄化対象物質の各々について算出される際、触媒モデル３１０には先ず、ＲＯＭ等に予め記憶されている、三元触媒２２３の形状に関する物性値、例えばセル開口面積、セル数及び熱伝達率等の指標値が入力される。そして、触媒流入ガスモデル３２０から、前述した、触媒流入ガスの状態を規定する指標値が入力される。
【００５６】
これら物性値及び指標値が入力されると、触媒モデル３１０は、（１）式に係る数値演算処理を、設定されたセル数に対応する頻度で繰り返す。より具体的には、触媒モデル３１０は、例えば第１のセルについて、初期値となるＮＯｘ濃度（即ち、触媒流入ガスのＮＯｘ濃度）及び（１）式に基づいた数値演算により第１のセルについてのＮＯｘ浄化率を取得し、係る取得された浄化率に基づいて、第２のセルにおける浄化率推定に係る数値演算処理に供すべきＮＯｘ濃度を取得する。
【００５７】
触媒モデル３１０は、更に同様の処理を設定されたセル数分繰り返し、最終的に三元触媒２２３の出口における、即ち触媒排出ガスにおけるＮＯｘ濃度を算出する。触媒排出ガスにおけるＮＯｘ濃度を算出すると、触媒モデル３１０は、係るＮＯｘ濃度に基づいてエミッション値を算出し、エミッション状態推定結果として出力する。触媒モデル３１０では、浄化対象物質の各々について、このような処理が並行して実行されている。
【００５８】
一方、三元触媒２２３は、気筒２０１から排出される排出ガス（実質的には、触媒流入ガス）によって高温となり、触媒温度や酸素濃度（或いは、酸素濃度を規定するものとしての空燃比又は吸入空気量等）等の熱負荷要因により劣化する。触媒劣化モデル３３０は、触媒モデル３１０から入力される、触媒温度及び酸素濃度等の熱負荷指標値に基づいて、上記（１）式を補正するための熱負荷補正係数を算出する。算出された熱負荷補正係数は、触媒モデル３１０に入力される。
【００５９】
尚、（１）式において、熱負荷の影響は、主として頻度因子A及び活性化エネルギＥに現れる。より具体的には、熱負荷が大きい程、頻度因子Ａは小さく、且つ活性化エネルギＥは大きくなる。即ち、熱負荷が大きい程、（１）式により規定されるＮＯｘ浄化率は低下する。触媒劣化モデル３３０は、熱負荷に応じた頻度因子Ａ及び活性化エネルギＥの変化特性を総括的に表すものとして熱負荷補正係数を算出し、触媒モデル３１０に出力する。
【００６０】
触媒モデル３１０は、上述した触媒排出ガスのエミッション特性の推定に際し、触媒劣化モデル３３０から入力される熱負荷補正係数を（１）式に乗じることによって（１）式を補正する。従って触媒モデル３１０は、三元触媒２２３に付与される熱負荷を考慮した、正確なエミッション値を出力することができる。
【００６１】
上述したように、上記（１）式を構成する頻度因子Ａ、抑制項Ｇ、活性化エネルギＥ、及びべき乗指数ｋの値を実情に適した値に設定することが出来れば、（１）式に基づいて触媒モデル３１０が実行する数値演算処理によって、実際の三元触媒２２３における浄化対象物質の浄化率を正確に近似することが可能となり、最終的に触媒モデル３１０から出力される触媒排出ガスのエミッション値を、実際に触媒排出ガスのエミッションを実測することなく正確に推定することが可能となる。
【００６２】
尚、係る事情に鑑みれば、触媒モデル３１０におけるエミッション状態の推定精度は、反応速度定数Ｋの同定精度によって大きく影響される。従って、反応速度定数Ｋの同定は、実機の触媒（例えば、エンジン２００に搭載される前の三元触媒２２３或いはそれと同等の触媒等）から得られる、例えば空燃比等の運転条件に応じた浄化率の実測値に適合するように、予めエミッション推定モデル３００がＥＣＵ１００内に構築される以前に行われている。
【００６３】
＜１−２−２：酸素濃度を考慮した補正＞
一方、三元触媒２２３の内部に酸素が十分に存在する場合、係る潤沢に存在する酸素によってＮＯｘからの酸素の脱離が阻害されるため、ＮＯｘの触媒反応は遅滞し易い。即ち、浄化率が低下し易い。他方、触媒流入ガスの流速が相対的に高くなるエンジン２００の高負荷領域（高流量域）では、触媒流入ガス中の絶対的な酸素量が多いことに加え、触媒流入ガスの流速が高いことに起因する他の酸化反応の反応速度の低下等により、三元触媒２２３内部の酸素濃度が高くなり易い。
【００６４】
従って、例えば、エンジン２００の高負荷領域では、顕著には触媒入り口に近い程ＮＯｘの浄化が進行し難い状況が発生する。このため、高負荷領域では、酸素濃度によってはＮＯｘが三元触媒２２３内を浄化されないまま通過する、所謂「ＮＯｘの吹き抜け」と称されるような事態が発生する。
【００６５】
ところが、上述した（１）式によって規定される浄化率をＮＯｘに適用する場合、（１）式にこのような酸素濃度の影響が記述されていないため、係るＮＯｘの吹き抜け現象を表すことが難しい。即ち、（１）式に従ってＮＯｘの浄化率を規定する場合、高負荷領域において顕著に実際の浄化率から乖離し易く、触媒排出ガスのエミッション特性の推定精度が極端に低下し易い。
【００６６】
ここで、図３を参照して、三元触媒２２３におけるＮＯｘの浄化特性について説明する。ここに、図３は、三元触媒２２３におけるＮＯｘ及びＯ２の濃度変化を概念的に表してなる模式図である。
【００６７】
図３において、縦軸には三元触媒２２３内のガス濃度が、また横軸には三元触媒２２３の入り口（即ち、三元触媒２２３のフロント側端部であり、以下、適宜「触媒入り口」と称する）からの長さＬが表されている。
【００６８】
図３において、酸素濃度及びＮＯｘ濃度の夫々について、エンジン２００の相互に異なる二種類の負荷領域、即ち、低負荷（低流量）領域及び高負荷（高流量）領域に対応する濃度変化の特性（即ち、浄化特性の一例）が示される。
【００６９】
低負荷領域における酸素濃度の変化特性は、図示ＰＲＦＯ２１（細実線）として示される。即ち、触媒入り口（図示長さＬが０の位置）において、酸素濃度は初期値Ｂを採り、三元触媒２２３の出口（即ち、リア側端部であり、以下、適宜「触媒出口」と称する）、即ち、図示長さＬｍａｘの位置において濃度は略ゼロとなる。
【００７０】
低負荷領域におけるＮＯｘ濃度の変化特性は、図示ＰＲＦＮＯｘ１（太実線）として表される。ここで、上述したように、触媒内部では、酸素濃度がある程度減少した後にＮＯｘの浄化が開始される。このＮＯｘの浄化が開始される酸素濃度をＤ（Ｄ＜Ｂ）とすると、ＮＯｘの浄化が開始される触媒内の位置は、図示長さＬ１に相当する位置となる。即ち、触媒入り口から長さＬ１に至るまでの領域では、ＮＯｘは浄化されずに（実践上有効であるとみなされる程度には浄化されずに）触媒内を通過する。また、触媒入り口からの長さがＬ１となる位置からＮＯｘの浄化が開始されると、ＮＯｘ濃度は徐々に減少し、触媒出口で概ねゼロとなる。
【００７１】
一方、高負荷領域における酸素濃度の変化特性は、図示ＰＲＦＯ２２（図示細破線）として示される。図示の通り、高負荷領域では触媒流入ガスの流速が相対的に高いため、酸素の浄化速度が相対的に緩やかとなり、上述した濃度Ｄまで酸素濃度が低下するのは、触媒入り口から長さＬ２の位置である。
【００７２】
高負荷領域におけるＮＯｘ濃度の特性は、図示ＰＲＦＮＯｘ２（図示太破線）となる。即ち、低負荷領域の場合と同様の考え方により、ＮＯｘの濃度が低下し始めるのは触媒入り口から長さＬ２の位置であり、触媒出口においてＮＯｘ濃度はＣ（Ｂ＜Ｃ＜Ｄ）となる。このように、高負荷領域では、酸素濃度の減少速度が緩慢になるため、触媒内部でＮＯｘが十分に浄化されず、相当程度のＮＯｘが浄化されぬまま排出され、触媒排出ガスのエミッションが悪化するといった事態が発生する。
【００７３】
このようなＮＯｘの吹き抜け現象は、先に述べた（１）式では十分に表し得ないため、顕著には高負荷領域において、触媒モデル３１０に係るエミッションの推定結果と、実際のエミッション特性とが大きく乖離することになる。
【００７４】
そこで、本実施形態では、ＮＯｘの浄化率を規定する式が、下記（３）式によって表される。即ち、触媒モデル３１０は、ＮＯｘの浄化率を推定する場合に、上記（１）式を、酸素濃度の影響を考慮すべく（３）式に補正する。尚、（３）式において、ｎは、ｋと同様に実際の浄化率に適合するように同定されるべき乗指数であり、基本的には１に設定される。
【００７５】
ｄ［ＮＯｘ］／ｄｔ＝Ｋ・［ＮＯｘ］＾ｋ／（１＋［Ｏ２］＾ｎ）・・・（３）
係る（３）式によれば、酸素濃度が高い領域ではＮＯｘの浄化率が低くなり、また酸素濃度が十分に低くなれば、実質的にＮＯｘの浄化が開始される。即ち、図３に概念的に示される触媒内のＮＯｘの濃度変化特性を高精度に近似することが可能となる。
【００７６】
尚、酸素濃度の影響を考慮する際の補正の態様は無論これに限定されず、例えば（３）式において、酸素濃度の項を、酸素濃度に例えば「Ａ１ｅｘｐ（−Ｅ１／ＲＴ）」なる定数を乗じてなる項に置換してもよい。ここで、Ａ１及びＥ１は、上述した反応速度定数と同様に頻度因子と活性化エネルギを表し、事前に実際のＮＯｘの浄化率と適合するように同定される。
【００７７】
ここで、図４を参照し、（３）式に基づいて得られるＮＯｘ濃度の推定結果について説明する。ここに、図４は、吸入空気量Ｇａに対する触媒排出ガス中のＮＯｘ濃度の特性を概念的に表す模式図である。
【００７８】
図４において、図３に示した実際のＮＯｘ濃度の変化特性に相当する特性が、図示ＰＲＦＮＯｘ３（実線）として示される。即ち、既に説明したように、吸入空気量Ｇａが高い領域（高負荷領域）において、触媒排出ガス中のＮＯｘ濃度は急激に上昇する。
【００７９】
一方、式（３）に従って触媒モデル３１０により推定されたＮＯｘ濃度は、図示ＰＲＦＮＯｘ４（図示二点鎖線）として表される。即ち、式（１）に更に酸素濃度の補正項が適用されるため、ＰＲＦＮＯｘ４は、ＰＲＦＮＯｘ３に漸近し、吸入空気量Ｇａが高い領域において急激に上昇する。
【００８０】
ここで、図４には更に、本実施形態に係る比較例として、図示ＰＲＦＮＯｘ５（図示破線）及びＰＲＦＮＯｘ６（図示一点鎖線）が示される。両特性は、いずれも式（１）に従って触媒モデル３１０により推定されたＮＯｘ濃度の変化特性であり、いずれも酸素濃度の影響が考慮されていないため、吸入空気量Ｇａの変化に対しＮＯｘ濃度の変化がリニアとなる。
【００８１】
従って、ＰＲＦＮＯｘ５を用いて高負荷領域におけるＮＯｘ濃度の推定精度を高めようとすれば、低負荷領域（吸入空気量Ｇａが低い領域）におけるＮＯｘ濃度の推定精度が低下し、ＰＲＦＮＯｘ６を用いて低負荷領域におけるＮＯｘ濃度の推定精度を高めようとすれば、高負荷領域（吸入空気量Ｇａが高い領域）におけるＮＯｘ濃度の推定精度が低下する。即ち、式（１）を用いた場合、酸素濃度に対するＮＯｘ濃度の変化を十分に推定することが難しいことが明らかである。
【００８２】
その点、ＰＲＦＮＯｘ４として示される通り、酸素濃度の補正項により酸素濃度の影響が考慮された場合には、酸素濃度（近似的には吸入空気量Ｇａ）に応じたＮＯｘの濃度変化を高精度に近似することが可能となり、触媒モデル３１０によって、触媒排出ガス中のエミッション特性を高精度に推定することが可能となるのである。
【００８３】
＜２：第２実施形態＞
第１実施形態では、三元触媒２２３におけるＮＯｘの浄化率が、（３）式により、ＮＯｘ濃度と反応速度定数Ｋとの積を含む形で定義された。しかしながら、ＮＯｘの浄化率は、還元剤として機能するＣＯやＨＣ等の濃度とも関係し得るため、ＮＯｘの濃度変化を推定する場合、これら還元剤の影響を考慮することにより一層推定精度が向上する可能性がある。
【００８４】
そこで、第２実施形態として、触媒モデル３１０は、第１実施形態に係る（３）式に代えて、下記（４）式により表されるＮＯｘ浄化率に基づいて触媒排出ガス中のＮＯｘ濃度を推定する。
【００８５】
浄化率＝ｄ［ＮＯｘ］／ｄｔ＝Ｋ１・［ＮＯｘ］＾ｋ１・［第１還元剤］＾ｍ１／（１＋［Ｏ２］＾ｎ１）＋Ｋ２・［ＮＯｘ］＾ｋ２・［第２還元剤］＾ｍ２／（１＋［Ｏ２］＾ｎ２）・・・（４）
ここで、［第１還元剤］及び［第２還元剤］とは、夫々第１還元剤（例えば、ＣＯ又はＨＣ）の濃度及び第２還元剤（例えば、ＨＣ又はＣＯ）の濃度を表す。また、Ｋ１及びＫ２は夫々反応速度定数であり、ｋ１、ｋ２、ｍ１、ｍ２はｎ１及びｎ２は、夫々べき乗指数である。
【００８６】
係る（４）式に従えば、第１実施形態と同様に酸素濃度の影響を考慮しつつ、更に還元剤毎に還元剤の濃度が考慮されるため、実際のＮＯｘの浄化率が一層正確に近似される。但し、この際、ＥＣＵ１００内に構築される前段階である係数同定に係る負荷が増大するため、エミッション特性の推定に際し、触媒モデル３１０がいずれの式を用いるべきかは、エンジン２００の仕様、仕向け及び要求性能等に応じてエンジン２００毎に個別具体的に決定されてよい。いずれの場合であっても、（１）式と比較すればＮＯｘの浄化率を正確に表すことが可能となるから、三元触媒２２３におけるＮＯｘの浄化率を高精度に推定することが可能となり、最終的に触媒排出ガスにおけるエミッション特性、顕著にはＮＯｘの排出量の推定に係る精度が飛躍的に向上する。
【００８７】
本発明は、上述した実施形態に限られるものではなく、請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う排出ガスの状態推定装置もまた本発明の技術的範囲に含まれるものである。
【図面の簡単な説明】
【００８８】
【図１】本発明の第１実施形態に係るエンジンシステムを表す模式図である。
【図２】図１のエンジンシステムにおいてＥＣＵ内に構築されたエミッション推定モデルの構成を概念的に表してなる模式図である。
【図３】三元触媒におけるＮＯｘ及びＯ２の濃度特性を概念的に表してなる模式図である。
【図４】吸入空気量に対する触媒排出ガス中のＮＯｘ濃度の特性を概念的に表す模式図である。
【符号の説明】
【００８９】
１０…エンジンシステム、１００…ＥＣＵ、２００…エンジン、２０１…気筒、２０３…ピストン、２０５…クランクシャフト、２２２…空燃比センサ、２２３…三元触媒、３００…エミッション推定モデル、３１０…触媒モデル、３２０…触媒流入ガスモデル、３３０…触媒劣化モデル。
【図１】
2008069708.tif 000003
【図２】
2008069708.tif 000004
【図３】
2008069708.tif 000005
【図４】
2008069708.tif 000006