(57)【要約】  推定値に基づいてプラントの制御量を制御するプラントの制御装置であって、プラントの固体ばらつきや経年劣化による推定値の誤差を抑制できるプラントの制御装置を提供すること。  排気浄化システムの制御装置は、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴを、入力ベクトルＵに基づきニューラルネットワークにより算出するＩｎｅｒｔ−ＥＧＲ推定値算出部７１１と、Ｉｎｅｒｔ−ＥＧＲ量と相関のある排気空燃比の推定値ΦＨＡＴを、入力ベクトルＵに基づきニューラルネットワークにより算出するＬＡＦセンサ出力推定値算出部７１２と、排気空燃比を検出するＬＡＦセンサ３４と、ＬＡＦセンサ３４の検出値ΦＡＣＴとＬＡＦセンサ３４の出力の推定値ΦＨＡＴとの推定誤差ＥＨＡＴが最小になるように、適応入力ＵＶＮＳを算出する非線形適応修正器７１３と、を備える。
2010109667.tif 000002
【特許請求の範囲】
【請求項１】
プラントの制御装置であって、
前記プラントの状態を示す複数の物理量のうちの少なくとも１つである第１物理量の推定値を、複数の入力に基づき所定のアルゴリズムにより算出する第１推定値算出手段と、
前記第１物理量と相関のある第２物理量の推定値を、複数の入力に基づき所定のアルゴリズムにより算出する第２推定値算出手段と、
前記第２物理量を検出する検出手段と、
前記検出手段により検出された第２物理量の検出値と前記第２推定値算出手段により算出された第２物理量の推定値との偏差が最小になるように、前記第１推定値算出手段及び前記第２推定値算出手段に入力される適応入力を算出する適応入力算出手段と、を備え、
前記第１物理量の推定値に基づいて、前記プラントの所定の制御量を制御することを特徴とするプラントの制御装置。
【請求項２】
前記第１推定値算出手段のアルゴリズム及び前記第２推定値算出手段のアルゴリズムは、それぞれ、所定の関数に従って出力する複数のニューロンを結合して構成されたニューラルネットワークであることを特徴とする請求項１に記載のプラントの制御装置。
【請求項３】
前記第１推定値算出手段への複数の入力、及び、前記第２推定値算出手段への複数の入力には、それぞれ、複数の異なる時刻の物理量に関するデータが含まれることを特徴とする請求項２に記載のプラントの制御装置。
【請求項４】
前記適応入力算出手段は、
前記第１推定値算出手段への複数の入力、及び、前記第２推定値算出手段への複数の入力のうちの少なくとも１つを参照パラメータとし、当該参照パラメータを基底とする空間に、互いに重複する複数の領域を定義するとともに、各領域にそれぞれ「０」でない値を持つ正規化された複数の重み関数を設定する重み関数設定手段と、
前記重み関数の値と前記偏差との積が最小になるように、前記領域ごとに修正値（Ｕｉｊ）を算出する修正値算出手段と、
前記重み関数の値と前記修正値との積の全領域に亘る総和に基づいて適応入力を決定する決定手段と、を備えることを特徴とする請求項１から３の何れかに記載のプラントの制御装置。
【請求項５】
前記プラントは、
内燃機関の排気系を流通する排気の一部を、前記内燃機関の吸気系に還流する排気還流装置を備えた内燃機関の排気浄化システムであり、
前記プラントの第１物理量は、前記排気還流装置により前記内燃機関に還流される排気に関するパラメータを含むことを特徴とする請求項１から４の何れかに記載のプラントの制御装置。
【請求項６】
前記排気還流装置は、排気系を流通する排気の一部を吸気系に還流する排気還流通路と、当該排気還流通路に設けられた排気還流制御弁とを備え、
前記制御装置は、前記還流される排気に関するパラメータの推定値が所定の目標値に一致するように、前記排気還流制御弁の操作量を決定するコントローラをさらに備えることを特徴とする請求項５に記載のプラントの制御装置。
【請求項７】
前記プラントは、
内燃機関の排気系に設けられ、還元剤の存在下で前記排気系を流通するＮＯｘを還元する選択還元触媒と、
前記排気系のうち前記選択還元触媒の上流側に、還元剤又は還元剤の元となる添加剤を供給する還元剤供給手段と、を備える内燃機関の排気浄化システムであり、
前記プラントの第１物理量は、前記選択還元触媒に流入する排気中のＮＯｘに関するパラメータを含むことを特徴とする請求項１から４の何れかに記載のプラントの制御装置。
【請求項８】
前記排気中のＮＯｘに関するパラメータの推定値に基づいて、前記還元剤供給手段による還元剤又は添加剤の供給量を決定するコントローラをさらに備えることを特徴とする請求項７に記載のプラントの制御装置。
【請求項９】
前記プラントは、
内燃機関の排気系に設けられ、前記内燃機関で燃焼する混合気を理論空燃比よりもリーンにしたときに排気中のＮＯｘを吸着又は吸蔵し、還元雰囲気下で前記吸着又は吸蔵したＮＯｘを還元するＮＯｘ浄化触媒と、
前記ＮＯｘ浄化触媒に流入する排気を還元雰囲気にする還元制御処理を実行する還元化手段と、を備える内燃機関の排気浄化システムであり、
前記プラントの第１物理量は、前記ＮＯｘ浄化触媒に流入する排気中のＮＯｘに関するパラメータ、及び、前記ＮＯｘ浄化触媒に流入する排気中の還元成分に関するパラメータを含むことを特徴とする請求項１から４の何れかに記載の制御装置。
【請求項１０】
前記ＮＯｘに関するパラメータ及び前記還元成分に関するパラメータの推定値に基づいて、前記還元制御処理の実行を指令するコントローラをさらに備えることを特徴とする請求項９に記載の制御装置。
【請求項１１】
前記第２物理量は、前記排気系を流通する排気の空燃比であることを特徴とする請求項５から１０の何れかに記載のプラントの制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プラントの制御装置に関する。特に、プラントの状態を示す所定の物理量の推定値を算出し、この推定値に基づいてプラントの制御量を制御するプラントの制御装置に関する。
【背景技術】
【０００２】
内燃機関やその排気浄化システムなどのプラントにおいて、このプラントの状態を示す複数の物理量は、センサにより検出されている物理量（検出物理量）と、センサにより検出されていない物理量（非検出物理量）とに分けられる。ここで、非検出物理量には、原理的にセンサで直に検出できない物理量や、耐久性やコストなどの様々な理由からセンサを用いて直に検出しない物理量などが含まれる。このような非検出物理量に基づいてプラントを制御する必要が生じた場合、制御装置では、基本的には他の検出物理量に基づいて、その推定値を算出する。【０００３】
非検出物理量の具体例の１つとして、例えば内燃機関の排気浄化システムでは、ＥＧＲ量やＥＧＲ率などの排気還流装置に関わる物理量が挙げられる。特許文献１や特許文献２には、ＥＧＲ量やＥＧＲ率の推定値の算出や、この推定値に基づいた制御に関する技術が示されている。【０００４】
特許文献１には、各種バルブの開度と、吸入新気量と、機関回転数と、ＥＧＲ率との関係が定められたマップにより、ＥＧＲ率の推定値を算出する制御装置が示されている。この制御装置では、ＥＧＲ率の推定値が目標ＥＧＲ率に維持するようにディーゼルパティキュレートフィルタ（以下、「ＤＰＦ（Ｄｉｅｓｅｌ  Ｐａｒｔｉｃｕｌａｔｅ  Ｆｉｌｔｅｒ）」という）を再生することにより、内燃機関から排出されるＮＯｘの量の変動を抑制している。【０００５】
特許文献２には、ＥＧＲ弁の開度やエアフローメータで検出された吸入空気量などに基づいて、ＥＧＲ量の推定値を算出する制御装置が示されている。この制御装置では、推定したＥＧＲ量に基づいて、燃料噴射時期やパイロット噴射量を補正することにより、内燃機関の過渡運転状態における燃焼騒音を低減したり、ＮＯｘの排出量を低減したりしている。
【特許文献１】特開２００８−１０６７１７号公報
【特許文献２】特開２００８−１９７８２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
ところで、特許文献１及び特許文献２に示された技術や、その他従来から知られている技術などにおいて、ＥＧＲ量やＥＧＲ率を推定する際には、予め設定されたマップや、物理モデルに基づいて予め設定された演算式が用いられる。しかしながら、このような予め設定されたマップや演算式を用いた場合、例えば、システムの経年劣化や固体ばらつきが生じると、推定値と実際の値との間に誤差が生じてしまい、結果として適切な制御が行えなくなってしまう。【０００７】
本発明は上述した点を考慮してなされたものであり、推定値に基づいてプラントの制御量を制御するプラントの制御装置であって、プラントの固体ばらつきや経年劣化による推定値の誤差を抑制できるプラントの制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するため本発明は、プラント（２，２Ａ，２Ｂ）の制御装置を提供する。前記プラントの制御装置は、前記プラントの状態を示す複数の物理量のうちの少なくとも１つである第１物理量の推定値（ＩＥＧＲＨＡＴ，ＮＯＸＨＡＴ，ＲＥＤＨＡＴ）を、複数の入力に基づき所定のアルゴリズムにより算出する第１推定値算出手段（７１１，８１１，９１１，９１５）と、前記第１物理量と相関のある第２物理量の推定値（ΦＨＡＴ）を、複数の入力に基づき所定のアルゴリズムにより算出する第２推定値算出手段（７１２，８１２，９１２）と、前記第２物理量を検出する検出手段（３４）と、前記検出手段により検出された第２物理量の検出値（ΦＡＣＴ）と前記第２推定値算出手段により算出された第２物理量の推定値（ΦＨＡＴ）との偏差（ＥＨＡＴ）が最小になるように、前記第１推定値算出手段及び前記第２推定値算出手段に入力される適応入力（ＵＶＮＳ）を算出する適応入力算出手段（７１３，８１３，９１３）と、を備える。前記プラントの制御装置は、前記第１物理量の推定値に基づいて、前記プラントの所定の制御量を制御する。【０００９】
この構成によれば、第１推定値算出手段で所定のアルゴリズムに基づいて第１物理量の推定値を算出し、第２推定値算出手段で所定のアルゴリズムに基づいて第２物理量の推定値を算出する。ここで、第１推定値算出手段及び第２推定値算出手段には、第２物理量の推定値と検出手段の出力値との偏差が最小になるように算出された適応入力が入力される。さらにこの第１物理量の推定値に基づいて、プラントの所定の制御量を制御する。【００１０】
ここで、例えば、プラントに固体ばらつきや経年劣化が生じることにより、第１物理量の推定値に誤差が発生したとする。この場合、第１物理量と相関のある第２物理量の推定値にも、誤差が発生すると考えられる。上記構成によれば、この第１物理量の推定値に発生する誤差は、第２物理量の推定値と検出手段の検出値との偏差として検出される。さらにこの偏差が最小になるように適応入力が算出され、第１推定値算出手段及び第２推定手段に入力される。これにより、第２物理量の推定値とともに第１物理量の推定値の誤差が抑制される。また、このような第１物理量の推定値に基づいてプラントの所定の制御量を制御することにより、プラントの制御量を適切な状態に制御することができる。【００１１】
好ましくは、前記第１推定値算出手段のアルゴリズム及び前記第２推定値算出手段のアルゴリズムは、それぞれ、所定の関数（ｆ（ｘ）、ｇ（ｘ））に従って出力する複数のニューロンを結合して構成されたニューラルネットワークである。【００１２】
この構成によれば、非線形な動特性の再現性に優れたニューラルネットワークにより第１物理量及び第２物理量の推定値を算出する。これにより、例えば、実際の第１物理量が非線形な挙動を示したとしても、これを精度良く推定することができる。【００１３】
好ましくは、前記第１推定値算出手段への複数の入力（Ｕ）、及び、前記第２推定値算出手段への複数の入力（Ｕ）には、それぞれ、複数の異なる時刻の物理量に関するデータが含まれる。【００１４】
この構成によれば、第１推定値算出手段及び第２推定値算出手段への入力に、複数の異なる時刻の物理量に関するデータを含めることで、推定値の動的挙動の再現性をより向上することができる。【００１５】
好ましくは、前記適応入力算出手段は、前記第１推定値算出手段への複数の入力、及び、前記第２推定値算出手段への複数の入力のうちの少なくとも１つを参照パラメータとし、当該参照パラメータを基底とする空間に、互いに重複する複数の領域を定義するとともに、各領域にそれぞれ「０」でない値を持つ正規化された複数の重み関数（Ｗｉｊ）を設定する重み関数設定手段（７１５）と、前記重み関数の値と前記偏差との積が最小になるように、前記領域ごとに修正値（Ｕｉｊ）を算出する修正値算出手段（７１６）と、前記重み関数の値と前記修正値との積の全領域に亘る総和（ΣΣＷｉｊＵｉｊ）に基づいて適応入力を決定する決定手段（７１７）と、を備える。【００１６】
この構成によれば、参照パラメータを基底とする空間に複数の領域を定義するとともに、各領域に重み関数を設定する。そして、重み関数と上記偏差との積が最小になるように、領域ごとに修正値を算出する。さらに、重み関数と上記修正値との積の全領域に亘る総和に基づいて適応入力を決定する。
ところで、プラントにおける経年劣化や固体ばらつきが、第１物理量の推定値の誤差に及ぼす影響は、プラントの状態に応じて異なったものになると考えられる。この構成によれば、プラントの状態を示す参照パラメータを基底とした空間内の領域ごとに修正値を算出することにより、プラントの状態ごとに異なる誤差への影響を考慮して適応入力を算出することができる。【００１７】
好ましくは、前記プラントは、内燃機関（１）の排気系を流通する排気の一部を、前記内燃機関の吸気系に還流する排気還流装置（４０，４５）を備えた内燃機関の排気浄化システム（２）である。前記プラントの第１物理量は、前記排気還流装置により前記内燃機関に還流される排気に関するパラメータを含む。【００１８】
この構成によれば、排気還流装置により内燃機関に還流される排気に関するパラメータを第１物理量とし、この第１物理量の推定値を上述の手順で算出する。排気還流装置により還流される排気に関するパラメータは、現存するセンサでは精度良く検出することができない。そこで、このパラメータの推定値を第１推定値算出手段で精度良く算出することにより、排気浄化システムを、その固体ばらつきや経年劣化に応じて適切な状態に制御することができる。
特に、排気還流装置により還流される排気に関するパラメータは、過渡時において非線形的な挙動を示す。このため、上述のようなニューラルネットワークに基づいて、上述のパラメータの推定値を算出することで、このような非線形的な挙動も再現することができる。【００１９】
好ましくは、前記排気還流装置は、排気系を流通する排気の一部を吸気系に還流する排気還流通路（４１，４６）と、当該排気還流通路に設けられた排気還流制御弁（４２，４７）とを備える。前記制御装置は、前記還流される排気に関するパラメータの推定値（ＩＥＧＲＨＡＴ）が所定の目標値（ＩＥＧＲＩＤＥＡＬ＿ＣＭＤ）に一致するように、前記排気還流制御弁の操作量（ＬＨＰ，ＬＬＰ）を決定するコントローラ（７，７２，７３，７４）をさらに備える。【００２０】
この構成によれば、コントローラは、排気還流装置により還流される排気に関するパラメータの推定値が所定の目標値に一致するように、排気還流制御弁の操作量を決定する。これにより、排気浄化システムの固体ばらつきや経年劣化に合わせて排気還流制御弁の操作量を適切に決定することができる。したがって、排気浄化システムの固体ばらつきや経年劣化によるＮＯｘの排出量の増加を抑制することができる。【００２１】
好ましくは、前記プラントは、内燃機関の排気系に設けられ、還元剤の存在下で前記排気系を流通するＮＯｘを還元する選択還元触媒（６１）と、前記排気系のうち前記選択還元触媒の上流側に、還元剤又は還元剤の元となる添加剤を供給する還元剤供給手段（６２）と、を備える内燃機関の排気浄化システム（２Ａ）である。前記プラントの第１物理量は、前記選択還元触媒に流入する排気中のＮＯｘに関するパラメータを含む。【００２２】
この構成によれば、選択還元触媒に流入する排気中のＮＯｘに関するパラメータを第１物理量として、この第１物理量の推定値を上述の手順で算出する。排気中のＮＯｘを検出するセンサは、現存するものでは、検出分解能や応答性が低く、また出力ばらつきが大きい。このため、排気中のＮＯｘに関するパラメータを、精度良く検出することができない。そこで、このような排気中のＮＯｘに関するパラメータの推定値を第１推定値算出手段で精度良く算出することにより、排気浄化システムを、その固体ばらつきや経年劣化に応じて適切な状態に制御することができる。【００２３】
好ましくは、前記プラントの制御装置は、前記排気中のＮＯｘに関するパラメータの推定値（ＮＯＸＨＡＴ）に基づいて、前記還元剤供給手段による還元剤又は添加剤の供給量（ＧＵＲＥＡ）を決定するコントローラ（８、８２，８３，８４，８５）をさらに備える。【００２４】
この構成によれば、コントローラは、排気中のＮＯｘに関するパラメータの推定値に基づいて、還元剤又は添加剤の供給量を決定する。これにより、排気浄化システムの固体ばらつきや経年劣化に合わせて、選択還元触媒におけるＮＯｘ浄化率を高く維持することができる。【００２５】
好ましくは、前記プラントは、内燃機関の排気系に設けられ、前記内燃機関で燃焼する混合気を理論空燃比よりもリーンにしたときに排気中のＮＯｘを吸着又は吸蔵し、還元雰囲気下で前記吸着又は吸蔵したＮＯｘを還元するＮＯｘ浄化触媒（６５）と、前記ＮＯｘ浄化触媒に流入する排気を還元雰囲気にする還元制御処理を実行する還元化手段（９）と、を備える内燃機関の排気浄化システム（２Ｂ）である。前記プラントの第１物理量は、前記ＮＯｘ浄化触媒に流入する排気中のＮＯｘに関するパラメータ、及び、前記ＮＯｘ浄化触媒に流入する排気中の還元成分に関するパラメータを含む。【００２６】
この構成によれば、ＮＯｘ浄化触媒に流入する排気中のＮＯｘ及び還元成分に関するパラメータを第１物理量として、この第１物理量の推定値を上述の手順で算出する。このような排気中のＮＯｘや還元成分は、上述のように現存するセンサでは精度良く検出することができない。そこで、このような排気中のＮＯｘや還元成分に関するパラメータの推定値を第１推定値算出手段で精度良く算出することにより、排気浄化システムを、その固体ばらつきや経年劣化に応じて適切な状態に制御することができる。【００２７】
好ましくは、前記プラントの制御装置は、前記ＮＯｘに関するパラメータ及び前記還元成分に関するパラメータの推定値（ＮＯＸＨＡＴ，ＲＥＤＨＡＴ）に基づいて、前記還元制御処理の実行を指令するコントローラ（９，９２，９３）をさらに備える。【００２８】
この構成によれば、コントローラは、排気中のＮＯｘ及び還元成分に関するパラメータの推定値に基づいて、還元制御処理の実行を指令する。これにより、排気浄化システムの固体ばらつきや経年劣化に合わせて、ＮＯｘ浄化触媒におけるＮＯｘ浄化率を高く維持することができる。【００２９】
好ましくは、前記第２物理量は、前記排気系を流通する排気の空燃比である。【００３０】
この構成によれば、排気の空燃比を第２物理量として、この第２物理量の推定値を算出する。上述の排気還流装置により還流される排気に関するパラメータ、排気中のＮＯｘに関するパラメータ、及び排気中の還元成分に関するパラメータは、何れも排気の空燃比と強い相関がある。このような排気の空燃比を第２物理量の推定値するとともに、この推定値と検出値との偏差が最小になるように適応入力を算出することにより、第１物理量の推定値の精度を向上することができる。
【図面の簡単な説明】
【００３１】
【図１】本発明の第１実施形態に係るエンジン及びその排気浄化システムと、その制御装置との構成を示す模式図である。
【図２】シリンダに吸入された吸気の成分を示す図である。
【図３】Ｉｎｅｒｔ−ＥＧＲ量と、エンジンから排出されるＮＯｘ、ＨＣ、及びスートの量との関係を示す図である。
【図４】上記実施形態に係る排気浄化システムの制御装置の構成を示すブロック図である。
【図５】上記実施形態に係る適応バーチャルセンサシステムの構成を示すブロック図である。
【図６】上記実施形態に係るＩｎｅｒｔ−ＥＧＲ推定値算出部のニューラルネットワーク構造を示す図である。
【図７】上記実施形態に係るシグモイド関数を示す図である。
【図８】上記実施形態に係る非線形適応修正器の構成を示すブロック図である。
【図９】上記実施形態に係るエンジン回転数を定義域とした４つの第１重み関数を示す図である。
【図１０】上記実施形態に係る燃料噴射量を定義域とした４つの第２重み関数を示す図である。
【図１１】上記実施形態に係る２つの参照パラメータを定義域とした１６個の重み関数を示す図である。
【図１２】上記実施形態に係るＩｎｅｒｔ−ＥＧＲ量制御の手順を示すフローチャートである。
【図１３】上記実施形態に係るシミュレーション結果を示す図である。
【図１４】上記実施形態に係るシミュレーション結果を示す図である。
【図１５】上記実施形態に係るシミュレーション結果を示す図である。
【図１６】上記実施形態に係るシミュレーション結果を示す図である。
【図１７】本発明の第２実施形態に係るエンジン及びその排気浄化システムの構成を示す模式図である。
【図１８】上記実施形態に係る選択還元触媒におけるＮＯｘ浄化率と、ＮＯ２とＮＯｘの比率との関係を示す図である。
【図１９】上記実施形態に係る排気浄化システムの制御装置の構成を示すブロック図である。
【図２０】上記実施形態に係る適応バーチャルセンサシステムの構成を示すブロック図である。
【図２１】上記実施形態に係る最大ストレージ容量と選択還元触媒温度との関係を示す図である。
【図２２】上記実施形態に係る選択還元触媒のストレージモデルの概念を示す模式図である。
【図２３】上記実施形態に係るシミュレーション結果を示す図である。
【図２４】上記実施形態に係るシミュレーション結果を示す図である。
【図２５】上記実施形態に係るシミュレーション結果を示す図である。
【図２６】上記実施形態に係るシミュレーション結果を示す図である。
【図２７】本発明の第３実施形態に係るエンジン及びその排気浄化システムの構成を示す模式図である。
【図２８】上記実施形態に係るＮＯｘ浄化触媒におけるＮＯｘの吸着／吸蔵効率と、ＮＯｘの吸着／吸蔵量との関係を示す図である。
【図２９】上記実施形態に係る排気浄化システムの制御装置の構成を示すブロック図である。
【図３０】上記実施形態に係る適応バーチャルセンサシステムの構成を示すブロック図である。
【図３１】上記実施形態に係るＮＯｘ浄化触媒における最大ＮＯｘ吸着／吸蔵量とＮＯｘ浄化触媒温度との関係を示す図である。
【図３２】上記実施形態に係るシミュレーション結果を示す図である。
【図３３】上記実施形態に係るシミュレーション結果を示す図である。
【図３４】上記実施形態に係るシミュレーション結果を示す図である。
【図３５】上記実施形態に係るシミュレーション結果を示す図である。
【符号の説明】
【００３２】
１…エンジン（内燃機関）
２，２Ａ，２Ｂ…排気浄化システム（プラント、排気浄化システム）
２０…吸気管（吸気系）
２１…吸気マニホールド（吸気系）
３０…排気管（排気系）
３１…排気マニホールド（排気系）
３４…ＬＡＦセンサ（検出手段）
４０…高圧ＥＧＲ装置（排気還流装置）
４１…高圧ＥＧＲ管（排気還流通路）
４２…高圧ＥＧＲバルブ（排気還流制御弁）
４５…低圧ＥＧＲ装置（排気還流装置）
４６…低圧ＥＧＲ管（排気還流通路）
４７…低圧ＥＧＲバルブ（排気還流制御弁）
６１…選択還元触媒
６２…ユリア噴射装置（還元剤供給手段）
６５…ＮＯｘ浄化触媒
７…ＥＣＵ
７１…適応バーチャルセンサシステム
７１１…Ｉｎｅｒｔ−ＥＧＲ推定値算出部（第１推定値算出手段）
７１２…ＬＡＦセンサ出力推定値算出部（第２推定値算出手段）
７１３…非線形適応修正器（適応入力算出手段）
７１５…重み関数設定部（重み関数設定手段）
７１６…局所適応入力算出部（修正値算出手段）
７１７…適応係数算出部（決定手段）
７２…Ｉｎｅｒｔ−ＥＧＲ量目標値算出部（コントローラ）
７３…Ｉｎｅｒｔ−ＥＧＲコントローラ（コントローラ）
７４…リフト量算出部（コントローラ）
８…ＥＣＵ
８１…適応バーチャルセンサシステム
８１１…ＮＯｘ量推定値算出部（第１推定値算出手段）
８１２…ＬＡＦセンサ出力推定値算出部（第２推定値算出手段）
８１３…非線形適応修正器（適応入力算出手段）
８２…フィードフォワード噴射量決定部（コントローラ）
８３…ストレージ量目標値設定部（コントローラ）
８４…フィードバック噴射量決定部（コントローラ）
８５…加算器（コントローラ）
９…ＥＣＵ
９１…適応バーチャルセンサシステム
９１１…ＮＯｘ量推定値算出部（第１推定値算出手段）
９１２…ＬＡＦセンサ出力推定値算出部（第２推定値算出手段）
９１３…非線形適応修正器（適応入力算出手段）
９１５…還元剤量推定値算出部（第１推定値算出手段）
９２…ＮＯｘ吸着／吸蔵量推定部（コントローラ）
９３…リッチモードコントローラ（コントローラ）
【発明を実施するための形態】
【００３３】
［第１実施形態］
以下、本発明の第１実施形態を、図面を参照して説明する。
図１は、本実施形態に係る内燃機関（以下「エンジン」という）１及びその排気浄化システム２と、その制御装置との構成を示す模式図である。
【００３４】
エンジン１は、リーンバーン運転方式のガソリンエンジン又はディーゼルエンジンであり、図示しない車両に搭載されている。エンジン１には、各シリンダの燃焼室内に燃料を噴射する燃料噴射弁が設けられている。これら燃料噴射弁は、電子制御ユニット（以下「ＥＣＵ」という）７により電気的に接続されており、燃料噴射弁の開弁時間及び閉弁時間は、ＥＣＵ７により制御される。
【００３５】
排気浄化システム２は、エンジン１に接続され吸気が流通する吸気管２０と、エンジン１の排気が流通する排気管３０と、排気の一部を吸気に還流する高圧排気還流装置（以下、「高圧ＥＧＲ装置」という）４０及び低圧排気還流装置（以下、「低圧ＥＧＲ装置」という）４５と、排気を浄化する酸化触媒３５及びＤＰＦ３６と、エンジン１に吸気を圧送するプライマリターボチャージャ５０及びセカンダリターボチャージャ５５と、を含んで構成される。
【００３６】
吸気管２０は、吸気マニホールド２１の複数の分岐部を介してエンジン１の各気筒の吸気ポートに接続されている。排気管３０は、排気マニホールド３１の複数の分岐部を介してエンジン１の各気筒の排気ポートに接続されている。
【００３７】
吸気管２０には、プライマリターボチャージャ５０と、セカンダリターボチャージャ５５と、インタークーラ５９とが上流側からこの順で設けられている。
【００３８】
プライマリターボチャージャ５０は、排気管３０に設けられたタービン５１と、吸気管２０に設けられたコンプレッサ５２と、を備える。タービン５１は、排気管３０を流通する排気の運動エネルギにより駆動される。コンプレッサ５２は、タービン５１の回転により駆動され、吸気を加圧する。さらに、プライマリターボチャージャ５０は、開閉動作によりタービン５１の回転速度を変更する図示しない可変ベーンを備える。
【００３９】
セカンダリターボチャージャ５５は、排気管３０のうちプライマリターボチャージャ５０のタービン５１よりも上流側に設けられたタービン５６と、吸気管２０のうちプライマリターボチャージャ５０のコンプレッサ５２よりも下流側に設けられたコンプレッサ５７と、を備える。タービン５６は、排気管３０を流通する排気の運動エネルギにより駆動される。コンプレッサ５７は、タービン５６の回転により駆動され、吸気を加圧する。さらに、セカンダリターボチャージャ５５は、開閉動作によりタービン５６の回転速度を変更する図示しない可変ベーンを備える。また、吸気管２０には、セカンダリターボチャージャ５５のコンプレッサ５７を迂回するバイパス管２２が設けられている。バイパス管２２には、このバイパス管２２を開閉するバイパス弁２３が設けられている。
【００４０】
インタークーラ５９は、２つのターボチャージャ５０，５５により加圧された吸気を冷却する。
【００４１】
排気管３０のうち、プライマリターボチャージャ５０のタービン５１の下流には、酸化触媒３５とＤＰＦ３６とが、上流側からこの順で設けられている。
【００４２】
酸化触媒３５は、排気との反応により発生する熱で排気を昇温する。この酸化触媒３５には、例えば、触媒として作用する白金（Ｐｔ）を、アルミナ（Ａｌ２Ｏ３）担体に担持させたものに、ＨＣの吸着作用に優れたゼオライトと、ＨＣの水蒸気改質作用に優れたロジウム（Ｒｈ）を加えて構成されたものが用いられる。
【００４３】
ＤＰＦ３６は、排気がフィルタ壁の微細な孔を通過する際、排気中の炭素を主成分とする粒子状物質（以下、「ＰＭ（Ｐａｒｔｉｃｕｌａｔｅ  Ｍａｔｔｅｒ）」という）を、フィルタ壁の表面及びフィルタ壁中の孔に堆積させることによって捕集する。フィルタ壁の構成材料としては、例えば、炭化珪素（ＳｉＣ）などのセラミックスの多孔体が使用される。
【００４４】
高圧ＥＧＲ装置４０は、高圧ＥＧＲ管４１と、高圧ＥＧＲバルブ４２と、高圧ＥＧＲクーラ４３と、高圧ＥＧＲインテークシャッタ４４Ｉと、高圧ＥＧＲエキゾーストシャッタ４４Ｅと、を含んで構成される。
【００４５】
高圧ＥＧＲ管４１は、排気マニホールド３１と吸気マニホールド２１とを接続する。高圧ＥＧＲバルブ４２は、高圧ＥＧＲ管４１に設けられ、この高圧ＥＧＲ管４１を介して還流される排気の流量を制御する。高圧ＥＧＲクーラ４３は、高圧ＥＧＲ管４１を介して還流される排気を冷却する。高圧ＥＧＲインテークシャッタ４４Ｉは、吸気管２０のうちインタークーラ５９の下流に設けられ、高圧ＥＧＲエキゾーストシャッタ４４Ｅは、排気管３０のうちタービン５６の上流に設けられている。
【００４６】
高圧ＥＧＲバルブ４２、高圧ＥＧＲインテークシャッタ４４Ｉ、及び高圧ＥＧＲエキゾーストシャッタ４４Ｅは、図示しないアクチュエータを介してＥＣＵ７に接続されており、その開度（リフト量）はＥＣＵ７により電磁的に制御される。
【００４７】
低圧ＥＧＲ装置４５は、低圧ＥＧＲ管４６と、低圧ＥＧＲバルブ４７と、低圧ＥＧＲクーラ４８と、低圧ＥＧＲインテークシャッタ４９Ｉと、低圧ＥＧＲエキゾーストシャッタ４９Ｅと、を含んで構成される。
【００４８】
低圧ＥＧＲ管４６は、排気管３０のうちＤＰＦ３６の下流側と吸気管２０のうちコンプレッサ５２の上流側とを接続する。低圧ＥＧＲバルブ４７は、低圧ＥＧＲ管４６に設けられ、この低圧ＥＧＲ管４６を介して還流される排気の流量を制御する。低圧ＥＧＲクーラ４８は、低圧ＥＧＲ管４６を介して還流される排気を冷却する。低圧ＥＧＲインテークシャッタ４９Ｉは、吸気管２０のうち低圧ＥＧＲ管４６の接続部よりも上流側に設けられ、低圧ＥＧＲエキゾーストシャッタ４９Ｅは、排気管３０のうち低圧ＥＧＲ管４６の接続部よりも下流側に設けられている。
【００４９】
低圧ＥＧＲバルブ４７、低圧ＥＧＲインテークシャッタ４９Ｉ、及び低圧ＥＧＲエキゾーストシャッタ４９Ｅは、図示しないアクチュエータを介してＥＣＵ７に接続されており、その開度（リフト量）はＥＣＵ７により電磁的に制御される。
【００５０】
ＥＣＵ７には、エンジン１のクランク軸の回転角度を検出するクランク角度位置センサ１１、及びエンジン１により駆動される車両のアクセルペダルの踏み込み量を検出するアクセルセンサ１２が接続されており、これらセンサの検出信号は、ＥＣＵ７に供給される。ここで、エンジン１の回転数ＮＥは、クランク角度位置センサ１１の出力に基づいてＥＣＵ７により算出される。エンジン１の負荷を示す燃料噴射量ＧＦＵＥＬは、アクセルセンサ１２の出力に基づいてＥＣＵ７により算出される。
【００５１】
これらセンサ１１，１２に加えて、ＥＣＵ７には、排気浄化システム２の各部分における物理量を検出する吸気圧力センサ２４、第１排気圧力センサ３２、第２排気圧力センサ３３、ＬＡＦセンサ３４、第１リフトセンサ１３、及び第２リフトセンサ１４が接続されている。
吸気圧力センサ２４は、吸気管２０のうちインタークーラ５９と高圧ＥＧＲインテークシャッタ４９Ｉとの間の吸気圧力Ｐ２を検出し、検出値に略比例した信号をＥＣＵ７に送信する。第１排気圧力センサ３２は、高圧ＥＧＲ管４１のうち高圧ＥＧＲクーラ４３の上流側の排気圧力Ｐ３を検出し、検出値に略比例した信号をＥＣＵ７に送信する。第２排気圧力センサ３３は、排気管３０のうちタービン５１と酸化触媒３５との間の排気圧力Ｐ４Ｌを検出し、検出値に略比例した信号をＥＣＵ７に送信する。ＬＡＦセンサ３４は、排気管３０のうちタービン５１と酸化触媒３５との間の排気の空燃比ΦＡＣＴを検出し、検出値に略比例した信号をＥＣＵ７に送信する。第１リフトセンサ１３は、高圧ＥＧＲバルブ４２のリフト量ＬＨＰ＿ＡＣＴを検出し、検出値に略比例した信号をＥＣＵ７に送信する。第２リフトセンサ１４は、低圧ＥＧＲバルブ４７のリフト量ＬＬＰ＿ＡＣＴを検出し、検出値に略比例した信号をＥＣＵ７に送信する。
【００５２】
ＥＣＵ７は、各種センサからの入力信号波形を整形し、電圧レベルを所定のレベルに修正し、アナログ信号値をデジタル信号値に変換するなどの機能を有する入力回路と、中央演算処理ユニット（以下「ＣＰＵ」という）とを備える。この他、ＥＣＵ７は、ＣＰＵで実行される各種演算プログラム及び演算結果などを記憶する記憶回路と、高圧ＥＧＲバルブ４２、低圧ＥＧＲバルブ４７、ターボチャージャ５０，５５、及びエンジン１の燃料噴射弁などに制御信号を出力する出力回路と、を備える。
【００５３】
次に、以上のような排気浄化システム２を制御するＥＣＵを構成するにあたり、本願発明者が着目した課題について、図２及び図３を参照して説明する。
【００５４】
図２は、シリンダに吸入された吸気の成分を示す図である。
図２に示すように、吸気は、新たに吸入された分（新気成分）と、高圧ＥＧＲ装置及び低圧ＥＧＲ装置により排気から還流された分（ＥＧＲ成分）とで構成される。ところで、ディーゼルエンジンやリーンバーンガソリンエンジンは、リーン燃焼を行うことにより、排気中には不活性ガスの他残留酸素が多く存在する。そこで、図２に示すように、ＥＧＲ成分のうち、この不活性ガスの成分のみを指してＩｎｅｒｔ−ＥＧＲという。また、この不活性ガスの量をＩｎｅｒｔ−ＥＧＲ量という。
【００５５】
図３は、以上のように定義されたＩｎｅｒｔ−ＥＧＲ量と、エンジンから排出される排気のＮＯｘ、ＨＣ、及びスートの濃度との関係を示す図である。図３に示すように、Ｉｎｅｒｔ−ＥＧＲ量と排気のＮＯｘ濃度、ＨＣ濃度、及びスート濃度との間には強い相関がある。
【００５６】
より具体的には、Ｉｎｅｒｔ−ＥＧＲ量を多くすると、混合気の比熱量が増加することにより燃焼温度が低下し、ＮＯｘ濃度は減少する。また、このような燃焼温度の低下に伴い、ＨＣやスートの濃度が増加する。また、さらにＩｎｅｒｔ−ＥＧＲ量を多くすると、燃焼温度の低下により未燃のＨＣが増加するため、スート濃度は減少する。そこで、これらＮＯｘ濃度、ＨＣ濃度、及びスート濃度を、できるだけ全て低減するためには、図３中、ハッチングで示す領域内に、Ｉｎｅｒｔ−ＥＧＲ量を制御することが好ましい。
【００５７】
図１に示すような排気浄化システム２では、還流される排気の温度をＥＧＲクーラ４３，４８で低下させることにより、不活性ガスの濃度を増加し、Ｉｎｅｒｔ−ＥＧＲ量を増加することができる。また、例えば、タービン５１，５６を駆動することで熱エネルギーを失い温度が低下した排気を、低圧ＥＧＲ管４６を介して還流することによっても、Ｉｎｅｒｔ−ＥＧＲ量を増加することができる。
【００５８】
ここで、排気のＮＯｘ濃度、ＨＣ濃度、及びスート濃度が最適になるように、図１の排気浄化システム２に対し、Ｉｎｅｒｔ−ＥＧＲ量のフィードバック制御を行うことを検討する。この場合、各燃焼サイクルにおいてシリンダに吸入されるＩｎｅｒｔ−ＥＧＲ量を検出する必要があるが、シリンダ内のＩｎｅｒｔ−ＥＧＲ量を必要な精度で検出でき、かつ、実際のエンジンの運転状況下での使用に耐えうる耐久性を有するセンサは存在しない。このため、上述の吸気圧力Ｐ２、排気圧力Ｐ３，Ｐ３Ｌや、その他の排気浄化システム２の状態を示す物理量に基づいて推定する必要がある。
【００５９】
従来では、物理量を推定する場合、この物理量の振る舞いを再現する物理モデルを構築し、このモデルに基づく演算式が用いられる。しかしながら、Ｉｎｅｒｔ−ＥＧＲ量を物理モデルに基づいて推定する場合、以下のような課題がある。
【００６０】
（１）先ず、実際のＩｎｅｒｔ−ＥＧＲ量は非線形で予測のしにくい挙動を示す。特に、過渡時における挙動は複雑であり、このような挙動を十分な精度で再現できる簡易な物理モデルを構築するのは容易ではない。
【００６１】
（２）また、このような物理モデルが構築できたとしても、排気浄化システム２の経年変化や固体ばらつきに対応することは困難である。したがって、例えば、ＥＧＲバルブやポート部、各インテークシャッタなどへのデポジットが生じた場合には、Ｉｎｅｒｔ−ＥＧＲ量の推定値が実際の値からずれてしまい、結果として排気の浄化性能が低下するおそれがある。
【００６２】
以下では、以上のような２つの課題（１）及び（２）に鑑みてなされた、排気浄化システム２の制御装置の構成について説明する。
【００６３】
図４は、排気浄化システム２の制御装置の構成を示すブロック図である。なお、図４には、排気浄化システム２におけるＩｎｅｒｔ−ＥＧＲ量のフィードバック制御に係る構成のみを図示する。より具体的には、高圧ＥＧＲバルブのリフト量ＬＨＰ及び低圧ＥＧＲバルブのリフト量ＬＬＰの決定に関する、ＥＣＵにより構成されるモジュールのみを図示する。
【００６４】
このモジュールは、適応バーチャルセンサシステム７１と、Ｉｎｅｒｔ−ＥＧＲ量目標値算出部７２と、Ｉｎｅｒｔ−ＥＧＲコントローラ７３と、リフト量算出部７４とを含んで構成される。
【００６５】
適応バーチャルセンサシステム７１は、複数のセンサ２４，３２，３３，３４，１３，１４の検出値Ｐ２，Ｐ３，Ｐ３Ｌ，ΦＡＣＴ，ＬＨＰ＿ＡＣＴ，ＬＬＰ＿ＡＣＴ、並びに燃料噴射量ＧＦＵＥＬ及びエンジン回転数ＮＥに基づいてＩｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴを算出する。なお、この適応バーチャルセンサシステム７１の詳細な構成については、後に図５〜図１１を参照して説明する。
【００６６】
Ｉｎｅｒｔ−ＥＧＲ量目標値算出部７２は、複数のセンサ２４，３２，３３，１３，１４の検出値Ｐ２，Ｐ３，Ｐ３Ｌ，ＬＨＰ＿ＡＣＴ，ＬＬＰ＿ＡＣＴ、並びに、エンジン回転数ＮＥ及び燃料噴射量ＧＦＵＥＬに基づいて、Ｉｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤを算出する。ここで、Ｉｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤは、上述のように排気のＮＯｘ、ＨＣ、及びスートの濃度をともに低減するように予め設定されたマップに基づいて決定される。
【００６７】
リフト量算出部７４は、Ｉｎｅｒｔ−ＥＧＲ量の目標値、並びに、エンジン回転数ＮＥ及び燃料噴射量ＧＦＵＥＬに対する高圧ＥＧＲバルブリフト量ＬＨＰ及び低圧ＥＧＲバルブリフト量ＬＬＰが設定されたマップを備えており、このマップに基づいてバルブリフト量ＬＨＰ，ＬＬＰを決定する。
【００６８】
Ｉｎｅｒｔ−ＥＧＲコントローラ７３は、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴとＩｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤとの偏差ＥＩＥ（下記式（１）参照）が「０」となるように、リフト量算出部７４の上述のマップに対する入力を修正する。より具体的には、以下に示すような応答指定型制御アルゴリズムに基づいて、Ｉｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤに対する修正目標値ＩＥＧＲＣＭＤを算出する。本実施形態では、以下に示すような応答指定型制御アルゴリズムに基づいて、修正目標値ＩＥＧＲＣＭＤを算出する。
【数１】
2010109667.tif 000003
【００６９】
ここで、記号（ｋ）は、離散化した時間を示す記号であり、所定の制御周期ごとに検出又は算出されたデータであることを示す。すなわち、記号（ｋ）が今回の制御タイミングにおいて検出又は算出されたデータであるとした場合、記号（ｋ−１）は前回の制御タイミングにおいて検出又は算出されたデータであることを示す。なお、以下の説明においては、記号（ｋ）を適宜、省略する。
【００７０】
先ず、下記式（２）に示すように、切換関数設定パラメータＰＯＬＥＩＥと前回制御時の偏差ＥＩＥ（ｋ−１）との積と、今回制御時の偏差ＥＩＥ（ｋ）との和を算出し、これを切換関数σＩＥ（ｋ）として定義する。なお、切換関数設定パラメータＰＯＬＥＩＥは、所定の設定テーブルに基づいて、−１から０の間で設定されたものが用いられる。
【数２】
2010109667.tif 000004
【００７１】
次に、切換関数σＩＥ（ｋ）に基づいて、到達則入力ＵＲＣＨ＿ＩＥ（ｋ）、及び適応則入力ＵＡＤＰ＿ＩＥ（ｋ）を算出する。より具体的には、到達則入力ＵＲＣＨ＿ＩＥ（ｋ）は、偏差状態量を切換直線上に載せるための入力であり、下記式（３）に示すように、切換関数σＩＥ（ｋ）に所定の到達則制御ゲインＫＲＣＨ＿ＩＥを乗算することで算出される。
【数３】
2010109667.tif 000005
【００７２】
適応則入力ＵＡＤＰ＿ＩＥ（ｋ）は、モデル化誤差や外乱の影響を抑制し、偏差状態量を切換直線に載せるための入力であり、下記式（４）に示すように、切換関数σＩＥ（ｋ）と所定の適応則ゲインＫＡＤＰ＿ＩＥを乗算したものと、前回制御時の適応則入力ＵＡＤＰ＿ＩＥ（ｋ−１）との和により算出される。
【数４】
2010109667.tif 000006
【００７３】
そして、下記式（５）に示すように、これらＵＲＣＨ（ｋ）、及びＵＡＤＰ（ｋ）の和を算出し、これを、補正係数ＫＥＧＲ（ｋ）として定義する。
【数５】
2010109667.tif 000007
【００７４】
さらに、下記式（６）に示すように、Ｉｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤ（ｋ）に、算出した補正係数ＫＥＧＲ（ｋ）を乗算することにより、Ｉｎｅｒｔ−ＥＧＲ量の修正目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤ（ｋ）を算出する。
【数６】
2010109667.tif 000008
【００７５】
以上のような構成により、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴが目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤに一致するように、高圧ＥＧＲバルブリフト量ＬＨＰ及び低圧ＥＧＲバルブリフト量ＬＬＰが決定される。
【００７６】
以下、図５〜図１１を参照して、適応バーチャルセンサシステム７１の構成について詳細に説明する。
図５は、適応バーチャルセンサシステム７１の構成を示すブロック図である。
【００７７】
適応バーチャルセンサシステム７１は、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴを算出するＩｎｅｒｔ−ＥＧＲ推定値算出部７１１と、ＬＡＦセンサ３４の出力（排気空燃比）の推定値ΦＨＡＴを算出するＬＡＦセンサ出力推定値算出部７１２と、適応入力ＵＶＮＳを算出する非線形適応修正器７１３とを含んで構成される。
【００７８】
Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１は、上述の課題（１）を解決するため、図５及び後述の図６に示すように、非線形な動特性の再現性に優れたニューラルネットワークを用いてＩｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴを算出する。
【００７９】
また、さらに、この適応バーチャルセンサシステム７１では、上述の課題（２）を解決するため、すなわち排気浄化システムの劣化や固体ばらつきに対応するため、Ｉｎｅｒｔ−ＥＧＲ量とは別の物理量であり、かつ、Ｉｎｅｒｔ−ＥＧＲ量と相関のある物理量である排気空燃比の推定値ΦＨＡＴを、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１と同様のニューラルネットワーク構造を有するＬＡＦセンサ出力推定値算出部７１２で算出する。
【００８０】
さらに、下記式（７）に示すように、算出した排気空燃比の推定値ΦＨＡＴ（ｋ）とＬＡＦセンサ３４の検出値ΦＡＣＴ（ｋ）との間の推定誤差ＥＨＡＴ（ｋ）を加算器７１４により算出する。
【数７】
2010109667.tif 000009
【００８１】
さらに、非線形適応修正器７１３では、後に図８〜図１１を参照して詳述するように、算出した推定誤差ＥＨＡＴが最小になるように、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１及びＬＡＦセンサ出力推定値算出部７１２に対し共通に入力される適応入力ＵＶＮＳを算出する。
【００８２】
この適応入力ＵＶＮＳは、例えば、「０」〜「１」の間に設定されるものであり、後に詳述するように、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１及びＬＡＦセンサ出力推定値算出部７１２において、排気浄化システムの劣化や固体ばらつきによる基準品からのずれを示す値として、それぞれのニューラルネットワークの学習時に設定される入力である。
【００８３】
すなわち、この適応バーチャルセンサシステム７１では、直感的には、排気浄化システムの劣化や固体ばらつきにより生じるＩｎｅｒｔ−ＥＧＲ量の推定値の誤差を、Ｉｎｅｒｔ−ＥＧＲ量と相関のある物理量の推定誤差ＥＨＡＴにより間接的に検出する。そして、この誤差が最小になるように、基準品からのずれを示す入力として予め用意しておいた適応入力ＵＶＮＳを算出し、この適応入力ＵＶＮＳをＩｎｅｒｔ−ＥＧＲ推定値算出部７１１及びＬＡＦセンサ出力推定値算出部７１２に入力する。これにより、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１及びＬＡＦセンサ出力推定値算出部７１２のニューラルネットワーク構造において、排気浄化システムの劣化や固体ばらつきに対する適応特性を実現することができる。
【００８４】
ここで、本実施形態の適応バーチャルセンサシステム７１におけるＩｎｅｒｔ−ＥＧＲ量と、排気空燃比との関係について説明する。この適応バーチャルセンサシステム７１では、出力として必要となるのはＩｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴである。これに対して、排気空燃比の推定値ΦＨＡＴは、上述の課題（２）を解決しＩｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴの真値との誤差を小さくするために、補助的に算出されるものである。
【００８５】
すなわち、排気浄化システムが劣化したり固体ばらつきが発生したりした場合には、Ｉｎｅｒｔ−ＥＧＲ量に及ぼされる影響と、補助的に検出される物理量に及ぼされる影響とがほぼ等しいことが好ましい。
【００８６】
したがって、このような補助的に検出する物理量としては、上述のように推定する必要のある物理量と上述のような相関があること、並びに、センサにより常時検出できる物理量であることが好ましい。このような条件を満たす物理量であれば、補助的に検出する物理量は、排気空燃比に限られない。
【００８７】
以下では、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１、ＬＡＦセンサ出力推定値算出部７１２、及び非線形適応修正器７１３の構成について順に説明する。
【００８８】
［Ｉｎｅｒｔ−ＥＧＲ推定値算出部］
図６は、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１のニューラルネットワーク構造を示す図である。
このニューラルネットワークは、所定の関数に従って出力する複数のニューロンを結合して構成され、ｍ成分の入力ベクトルＵ（ｋ）に応じて、値Ｙ（ｋ）を出力する。図６に示すように、このニューラルネットワークは、ｍ個のニューロンＷ１ｊ（ｊ＝１〜ｍ）で構成された入力層と、ｍ×（ｎ−１）個のニューロンＷｉｊ（ｉ＝２〜ｎ，ｊ＝１〜ｍ）で構成された中間層と、１個のニューロンＹで構成された出力層との３つの層を含んで構成された階層型である。
入力層：Ｗ１ｊ      （ｊ＝１，２，…，ｍ）
中間層：Ｗｉｊ      （ｉ＝２，３，…，ｎ，ｊ＝１，２，…，ｍ）
出力層：Ｙ
【００８９】
入力層のｍ個のニューロンＷ１ｊ（ｊ＝１〜ｍ）の動作について説明する。
入力層のニューロンＷ１ｊには、信号Ｔ１ｊ（ｋ）が入力される。この入力信号Ｔ１ｊ（ｋ）には、それぞれ、下記式（８）に示すように入力ベクトルＵ（ｋ）のｊ番目の成分Ｕｊ（ｋ）が用いられる。
【数８】
2010109667.tif 000010
【００９０】
入力層のニューロンＷ１ｊは、中間層のｍ個のニューロンＷ２ｊ（ｊ＝１〜ｍ）に所定の重みで結合しており、これら結合したｍ個のニューロンＷ２ｊへ信号Ｖ１ｊ（ｋ）を出力する。すなわち、このニューロンＷ１ｊは、下記式（９），（１０）に示すように、シグモイド関数ｆ（ｘ）に従って、入力信号Ｔ１ｊ（ｋ）に応じた信号Ｖ１ｊ（ｋ）をｍ個のニューロンＷ２ｊに出力する。
【数９】
2010109667.tif 000011
【数１０】
2010109667.tif 000012
【００９１】
図７は、シグモイド関数ｆ（ｘ）を示す図である。この図７には、上記式（１０）において、ε＝０とし、β＝０．５，１．０，２．０，３．０とした場合を示す。
シグモイド関数ｆ（ｘ）の値域は、［ε，ε＋１］となっている。また、図７に示すように、シグモイド関数ｆ（ｘ）は、βを大きくするに従い、ｘ＝０を中心としたステップ関数に近づく。
【００９２】
上記式（１０）において、係数βはシグモイド関数ｆ（ｘ）の傾きゲインを示し、係数εはシグモイド関数ｆ（ｘ）のオフセット値を示す。傾きゲインβは、後述のニューラルネットワークの学習により設定する。オフセット値εは、後述のニューラルネットワークの学習により設定するか、又は所定の値に設定しておく。
【００９３】
次に中間層の（ｎ−１）×ｍ個のニューロンＷｉｊ（ｉ＝２〜ｎ，ｊ＝１〜ｍ）の動作について説明する。
中間層のニューロンＷｉｊ（ｉ＝２〜ｎ，ｊ＝１〜ｍ）には、結合するニューロンから出力されたｍ個の信号Ｖｉ−１，ｊ（ｊ＝１〜ｍ）のそれぞれに所定の重みωｉ−１，ｊ（ｊ＝１〜ｍ）を乗じた信号の和が入力される。したがって、中間層のニューロンＷｉｊには、下記式（１１）に示すような信号Ｔｉｊ（ｋ）が入力される。
【数１１】
2010109667.tif 000013
【００９４】
中間層のニューロンのうち出力層に結合するｍ個を除いたニューロン、すなわち、（ｎ−２）×ｍ個のニューロンＷｉｊ（ｉ＝２〜ｎ−１，ｊ＝１〜ｍ）は、中間層のｍ個のニューロンＷｉ＋１，ｊ（ｊ＝１〜ｍ）に重みωｉｊで結合しており、これら結合したニューロンＷｉ＋１，ｊへ信号Ｖｉｊ（ｋ）を出力する。すなわち、このニューロンＷｉｊ（ｉ＝２〜ｎ−１，ｊ＝１〜ｍ）は、下記式（１２）に示すように、シグモイド関数ｆ（ｘ）に従って、入力信号Ｔｉｊ（ｋ）に応じた信号Ｖｉｊ（ｋ）をｍ個のニューロンＷｉ＋１，ｊに出力する。
【数１２】
2010109667.tif 000014
【００９５】
また、中間層のｍ個のニューロンＷｎｊ（ｊ＝１〜ｍ）は、出力層のニューロンＹに重みωｎｊで結合しており、この出力層のニューロンＹへ信号Ｖｎｊ（ｋ）を出力する。すなわち、これらニューロンＷｎｊ（ｊ＝１〜ｍ）は、下記式（１３）に示すように、シグモイド関数ｆ（ｘ）に従って、入力信号Ｔｎｊ（ｋ）に応じた信号Ｖｎｊ（ｋ）をニューロンＹに出力する。
【数１３】
2010109667.tif 000015
【００９６】
次に出力層のニューロンＹの動作について説明する。
出力層のニューロンＹには、結合する中間層のニューロンから出力されたｍ個の信号Ｖｎ，ｊ（ｊ＝１〜ｍ）に所定の重みωｎ，ｊ（ｊ＝１〜ｍ）を乗じた信号の和が入力される。したがって、出力層のニューロンＹには、下記式（１４）に示すような信号Ｔ（ｋ）が入力される。
【数１４】
2010109667.tif 000016
【００９７】
出力層のニューロンＹは、下記式（１５），（１６）に示すように、シグモイド関数ｇ（ｘ）に従って、入力信号Ｔ（ｋ）に応じた信号Ｙ（ｋ）を出力する。
【数１５】
2010109667.tif 000017
【数１６】
2010109667.tif 000018
【００９８】
シグモイド関数ｇ（ｘ）は、上述の図７に示す関数ｆ（ｘ）と、定性的には同じ振る舞いを示すが、値域が［δ，δ＋α］である点でシグモイド関数ｆ（ｘ）と異なる。上記式（１６）において、係数γはシグモイド関数ｇ（ｘ）の傾きゲインを示し、係数δはシグモイド関数ｇ（ｘ）のオフセット値を示す。また、係数αはニューラルネットワークの出力の取り得る自由度を設定するための出力ゲインを示す。傾きゲインγ及び出力ゲインαは、後述のニューラルネットワークの学習により設定する。オフセット値δは、後述のニューラルネットワークの学習により設定するか、又は所定の値に設定しておく。
【００９９】
次に、Ｉｎｅｒｔ−ＥＧＲ量を推定するための、ニューラルネットワークの学習について説明する。
【０１００】
先ず、ニューラルネットワークに対する入力ベクトルＵ（ｋ）の成分を、下記式（１７）に示すように定義する。このように、入力ベクトルＵ（ｋ）の成分には、Ｉｎｅｒｔ−ＥＧＲ量を推定するために必要となる複数の物理量（燃料噴射量ＧＦＵＥＬ、吸気圧力Ｐ２、排気圧力Ｐ３、排気圧力Ｐ３Ｌ、高圧ＥＧＲバルブリフト量の検出値ＬＨＰ＿ＡＣＴ、低圧ＥＧＲバルブリフト量の検出値ＬＬＰ＿ＡＣＴ、エンジン回転数ＮＥ）と、上述の適応入力ＵＶＮＳとが含まれる。また、入力ベクトルの成分には、このように異なる種類の物理量に関するデータが含まれているとともに、異なる時刻の物理量に関するデータも含まれている。
【数１７】
2010109667.tif 000019
【０１０１】
また、このような入力ベクトルＵ（ｋ）に対するニューラルネットワークの出力Ｙ（ｋ）を、下記式（１８）に示すように、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴ（ｋ）として定義する。
【数１８】
2010109667.tif 000020
【０１０２】
次に、Ｉｎｅｒｔ−ＥＧＲ量を実際に検出できるセンサと、状態の異なる排気浄化システムを少なくとも２組準備する。１つは、例えば基準となる新品の排気浄化システム（以下、単に「基準品」という）であり、もう１つは、固体ばらつきや劣化などの理由により、基準品に対し大きく異なった特性を有する排気浄化システム（以下、単に「劣化品」という）である。そして、実際に準備した排気浄化システムを運転することにより、上記式（１７）の入力ベクトルＵの成分（ＧＦＵＥＬ，Ｐ２，Ｐ３，Ｐ３Ｌ，ＬＨＰ＿ＡＣＴ，ＬＬＰ＿ＡＣＴ，ＮＥ）と、上述のセンサにより検出されたＩｎｅｒｔ−ＥＧＲ量の検出値との関係を記録する。なお、この入力ベクトルの成分とＩｎｅｒｔ−ＥＧＲ量の検出値との関係を示すデータは、準備した状態の異なる排気浄化システムごとに取得する。
【０１０３】
次に、取得したデータに基づいてニューラルネットワークの学習を行う。すなわち、入力ベクトルＵの成分（ＧＦＵＥＬ，Ｐ２，Ｐ３，Ｐ３Ｌ，ＬＨＰ＿ＡＣＴ，ＬＬＰ＿ＡＣＴ，ＮＥ）とＩｎｅｒｔ−ＥＧＲ量の検出値との関係が、ニューラルネットワークにより再現されるようにニューロンの関数ｆ（ｘ），ｇ（ｘ）の各種ゲイン（α，β，γ，δ，ε）、並びに、各ニューロンの結合の強さを示す重みωｉｊ（ｉ＝１〜ｎ，ｊ＝１〜ｍ）を設定する。なお、ニューラルネットワークの学習のアルゴリズムには、既知の方法が用いられる。具体的には、例えば、逆誤差伝播法などの学習アルゴリズムの他、遺伝的アルゴリズムなどの最適化アルゴリズムが挙げられる。
【０１０４】
ここで、上述のような学習を行う際における適応入力ＵＶＮＳの設定について説明する。例えば、データを取得するための排気浄化システムとして、基準品と劣化品との２つを準備した場合について説明する。先ず、基準品のデータに基づいて学習を行う際には、適応入力ＵＶＮＳを「１」と設定する。また、劣化品のデータに基づいて学習を行う際には、適応入力ＵＶＮＳを「０」と設定する。
基準品：適応入力ＵＶＮＳ←１
劣化品：適応入力ＵＶＮＳ←０
【０１０５】
以上のように設定された適応入力ＵＶＮＳを含む入力ベクトルＵを用いてニューラルネットワークの学習を行うことにより、適応入力ＵＶＮＳを「０」〜「１」の間で連続的に変化させた場合に、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴ及びＬＡＦセンサの出力の推定値ΦＨＡＴを基準品から劣化品へ連続的に変化させることができる。
【０１０６】
以上のようにして適応入力ＵＶＮＳを設定した上でニューラルネットワークの学習を行うことにより、排気浄化システムの劣化や固体ばらつきに対する適応特性が実現された適応バーチャルセンサシステムを構築することができる。
【０１０７】
また、例えば、基準品と劣化品との間の中間的な特性を有する排気浄化システムを準備した場合には、適応入力ＵＶＮＳを「１」と「０」との間の値、例えば「０．３」や「０．６」等の値に設定した上で、ニューラルネットワークの学習を行う。これにより、排気浄化システムの劣化や固体ばらつきに対する適応特性を、より現実的に実現することができる。
【０１０８】
［ＬＡＦセンサ出力推定値算出部］
図５に戻って、ＬＡＦセンサ出力推定値算出部７１２の構成について説明する。
ＬＡＦセンサ出力推定値算出部７１２は、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１と同様に、ニューラルネットワークによりＬＡＦセンサ３４の出力の推定値ΦＨＡＴを算出する。なお、ＬＡＦセンサ出力推定値算出部７１２のニューラルネットワーク構造は、図６及び図７を参照して詳述したＩｎｅｒｔ−ＥＧＲ推定値算出部７１１のニューラルネットワーク構造とほぼ同じであり、その詳細な説明を省略する。
【０１０９】
また、入力ベクトルＵの成分も、上記式（１７）と同様に定義する。すなわち、入力ベクトルＵの成分には、排気空燃比を推定するために必要となる複数の物理量（燃料噴射量ＧＦＵＥＬ、吸気圧力Ｐ２、排気圧力Ｐ３、排気圧力Ｐ３Ｌ、高圧ＥＧＲバルブリフト量の検出値ＬＨＰ＿ＡＣＴ、低圧ＥＧＲバルブリフト量の検出値ＬＬＰ＿ＡＣＴ、エンジン回転数ＮＥ）と、適応入力ＵＶＮＳとが含まれる。なお、本実施形態では、ＬＡＦセンサ出力推定値算出部のニューラルネットワークに対する入力ベクトルＵとしてＩｎｅｒｔ−ＥＧＲ推定値算出部と同じ入力ベクトルＵを用いたが、これに限るものではない。これら入力ベクトルの成分には互いに異なる物理量を含むようにしてもよいが、適応入力ＵＶＮＳは共通のものを用いる。
【０１１０】
さらに、下記式（１９）に示すように、ニューラルネットワークの出力Ｙ（ｋ）をＬＡＦセンサの出力の推定値ΦＨＡＴ（ｋ）として定義する。
【数１９】
2010109667.tif 000021
【０１１１】
また、ニューラルネットワークに対する学習も、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１のニューラルネットワークと同様の手順により行う。したがって、その詳細な説明を省略する。
【０１１２】
これにより、上述のＩｎｅｒｔ−ＥＧＲ推定値算出部７１１と同様に、排気浄化システムの劣化や固体ばらつきに対する適応特性が実現された適応バーチャルセンサシステムを構成することができる。
【０１１３】
［非線形適応修正器］
次に、図８〜図１１を参照して、非線形適応修正器７１３の構成について説明する。
図８は、非線形適応修正器７１３の構成を示すブロック図である。
【０１１４】
非線形適応修正器７１３は、重み関数設定部７１５と、局所適応入力算出部７１６と、適応係数算出部７１７と、を備える。非線形適応修正器７１３は、このような構成により推定誤差ＥＨＡＴが最小となるように適応入力ＵＶＮＳを算出する。
【０１１５】
より具体的には、この非線形適応修正器７１３では、２つの参照パラメータ（エンジン回転数ＮＥ及び燃料噴射量ＧＦＵＥＬ）を基底とする空間を定義するとともに、この空間を複数の領域に分ける。さらに、領域ごとに後述の局所適応入力Ｕｉｊ（ｉ＝１〜４，ｊ＝１〜４）を算出し、これら局所適応入力Ｕｉｊを後述の重み関数Ｗｉｊ（ｉ＝１〜４，ｊ＝１〜４）で重み結合することにより、適応入力ＵＶＮＳを算出する。
【０１１６】
ところで、排気浄化システムにおける経年劣化や固体ばらつきが、Ｉｎｅｒｔ−ＥＧＲ量や排気空燃比の推定値の誤差に及ぼす影響は、エンジンの運転条件、すなわち参照パラメータの値ごとに異なったものになると考えられる。この非線形適応修正器７１３では、参照パラメータを基底とする空間内の領域ごとに局所適応入力Ｕｉｊを算出することにより、参照パラメータの値ごとに異なる上述の誤差への影響を考慮して適応入力ＵＶＮＳを算出することができる。
【０１１７】
図９は、エンジン回転数ＮＥを定義域とした４つの第１重み関数ＷＮｉ（ｉ＝１〜４）を示す図である。図９に示すように、４つの第１重み関数ＷＮｉは、それぞれ、定義域に互いに重複した４つの領域を定義し、これら領域において「０」でない値を持つように設定される。
【０１１８】
より具体的には、定義域は、［Ｎ０，Ｎ２］と、第２領域［Ｎ１，Ｎ３］と、第３領域［Ｎ２，Ｎ４］と、第４領域［Ｎ３，Ｎ５］とに分けられる。ここで、図９に示すように、Ｎ０＜Ｎ１＜Ｎ２＜Ｎ３＜Ｎ４＜Ｎ５とする。したがって、第１領域と第２領域は区間［Ｎ１，Ｎ２］で重複し、第２領域と第３領域は区間［Ｎ２，Ｎ３］で重複し、第３領域と第４領域は区間［Ｎ３，Ｎ４］で重複する。
【０１１９】
関数ＷＮ１は、第１領域［Ｎ０，Ｎ２］において「０」でない値を持つように設定される。より具体的には、関数ＷＮ１は、区間［Ｎ０，Ｎ１］において「１」に設定され、区間［Ｎ１，Ｎ２］において「１」から「０」に減少するように設定される。
関数ＷＮ２は、第２領域［Ｎ１，Ｎ３］において「０」でない値を持つように設定される。より具体的には、関数ＷＮ２は、区間［Ｎ１，Ｎ２］において「１」から「０」に上昇するように設定され、区間［Ｎ２，Ｎ３］において「１」から「０」に減少するように設定される。したがって、関数ＷＮ１と関数ＷＮ２は、区間［Ｎ１，Ｎ２］の中心で交差する。
関数ＷＮ３は、第３領域［Ｎ２，Ｎ４］において「０」でない値を持つように設定される。より具体的には、関数ＷＮ３は、区間［Ｎ２，Ｎ３］において「１」から「０」に上昇するように設定され、区間［Ｎ３，Ｎ４］において「１」から「０」に減少するように設定される。したがって、関数ＷＮ２と関数ＷＮ３は、区間［Ｎ２，Ｎ３］の中心で交差する。
関数ＷＮ４は、第４領域［Ｎ３，Ｎ５］において「０」でない値を持つように設定される。より具体的には、関数ＷＮ４は、区間［Ｎ３，Ｎ４］において「１」から「０」に上昇するように設定され、区間［Ｎ４，Ｎ５］において「１」に設定される。したがって、関数ＷＮ３と関数ＷＮ４は、区間［Ｎ３，Ｎ４］の中心で交差する。
【０１２０】
また、図９には、Ｉｎｅｒｔ−ＥＧＲ量のエンジン回転数ＮＥに対する定性的な振る舞いを破線で示す。この図に示すように、低い回転数におけるＩｎｅｒｔ−ＥＧＲ量は、高い回転数におけるＩｎｅｒｔ−ＥＧＲ量よりも多くなる傾向がある。そこで、このようなＩｎｅｒｔ−ＥＧＲ量の非線形性に合わせて、本実施形態では、区間［Ｎ１，Ｎ２］は、区間［Ｎ３，Ｎ４］よりも狭くなるように設定する。
【０１２１】
また、以上のように構成された第１重み関数ＷＮｉは、下記式（２０）に示すように、その総和関数がエンジン回転数ＮＥによらず「１」となるように正規化される。
【数２０】
2010109667.tif 000022
【０１２２】
図１０は、燃料噴射量ＧＦＵＥＬを定義域とした４つの第２重み関数ＷＧｊ（ｊ＝１〜４）を示す図である。
図１０に示すように、４つの第２重み関数ＷＧｊは、それぞれ、定義域に互いに重複した４つの領域を定義し、これら領域において「０」でない値を持つように設定される。
【０１２３】
より具体的には、定義域は、［Ｇ０，Ｇ２］と、第２領域［Ｇ１，Ｇ３］と、第３領域［Ｇ２，Ｇ４］と、第４領域［Ｇ３，Ｇ５］とに分けられる。ここで、図１０に示すように、Ｇ０＜Ｇ１＜Ｇ２＜Ｇ３＜Ｇ４＜Ｇ５とする。したがって、第１領域と第２領域は区間［Ｇ１，Ｇ２］で重複し、第２領域と第３領域は区間［Ｇ２，Ｇ３］で重複し、第３領域と第４領域は区間［Ｇ３，Ｇ４］で重複する。
【０１２４】
関数ＷＧ１は、第１領域［Ｇ０，Ｇ２］において「０」でない値を持つように設定される。より具体的には、関数ＷＧ１は、区間［Ｇ０，Ｇ１］において「１」に設定され、区間［Ｇ１，Ｇ２］において「１」から「０」に減少するように設定される。
関数ＷＧ２は、第２領域［Ｇ１，Ｇ３］において「０」でない値を持つように設定される。より具体的には、関数ＷＧ２は、区間［Ｇ１，Ｇ２］において「１」から「０」に上昇するように設定され、区間［Ｇ２，Ｇ３］において「１」から「０」に減少するように設定される。したがって、関数ＷＧ１と関数ＷＧ２は、区間［Ｇ１，Ｇ２］の中心で交差する。
関数ＷＧ３は、第３領域［Ｇ２，Ｇ４］において「０」でない値を持つように設定される。より具体的には、関数ＷＧ３は、区間［Ｇ２，Ｇ３］において「１」から「０」に上昇するように設定され、区間［Ｇ３，Ｇ４］において「１」から「０」に減少するように設定される。したがって、関数ＷＧ２と関数ＷＧ３は、区間［Ｇ２，Ｇ３］の中心で交差する。
関数ＷＧ４は、第４領域［Ｇ３，Ｇ５］において「０」でない値を持つように設定される。より具体的には、関数ＷＧ４は、区間［Ｇ３，Ｇ４］において「１」から「０」に上昇するように設定され、区間［Ｇ４，Ｇ５］において「１」に設定される。したがって、関数ＷＧ３と関数ＷＧ４は、区間［Ｇ３，Ｇ４］の中心で交差する。
【０１２５】
また、図１０には、エンジンのＮＯｘ排出量の燃料噴射量ＧＦＵＥＬに対する定性的な振る舞いを破線で示す。この図に示すように、ＮＯｘ排出量は、燃料噴射量が所定量を超えると急激に増加する傾向がある。そこで、このようなＮＯｘ排出量の非線形性に合わせて、本実施形態では、区間［Ｇ３，Ｇ４］は、区間［Ｇ１，Ｇ２］よりも狭くなるように設定する。
【０１２６】
また、以上のように構成された第２重み関数ＷＮｉは、下記式（２１）に示すように、その総和関数が燃料噴射量ＧＦＵＥＬによらず「１」となるように正規化される。
【数２１】
2010109667.tif 000023
【０１２７】
図１１は、２つの参照パラメータ（ＮＥ，ＧＦＵＥＬ）を定義域とした１６個の重み関数Ｗｉｊ（ｉ＝１〜４，ｊ＝１〜４）を示す図である。図１１において、横軸はエンジン回転数ＮＥを示し、縦軸は燃料噴射量ＧＦＵＥＬを示す。図１１に示すように、２つの参照パラメータ（ＮＥ，ＧＦＵＥＬ）の定義域には、互いに重複する１６個の領域が定義される。
【０１２８】
１６個の重み関数Ｗｉｊは、下記式（２２）に示すように、第１重み関数ＷＮｉの各成分と第２重み関数ＷＧｊの各成分との積により定義される。これにより、１６個の領域においてそれぞれ「０」でない値を持つ重み関数Ｗｉｊが定義される。なお、図１１には、４つの重み関数Ｗ１１，Ｗ２２，Ｗ３３，Ｗ４４のみを図示する。
【数２２】
2010109667.tif 000024
【０１２９】
また、上記式（２０）及び（２１）と同様に、重み関数Ｗｉｊの総和関数は、下記式（２３）に示すように、２つの参照パラメータ（ＮＥ，ＧＦＵＥＬ）によらず「１」となるように正規化される。
【数２３】
2010109667.tif 000025
【０１３０】
図８に戻って、重み関数設定部７１５は、複数の第１重み関数ＷＮｉが設定された第１重み関数算出部７１５１と、複数の第２重み関数ＷＧｊが設定された第２重み関数算出部７１５２と、これら第１重み関数ＷＮｉ及び第２重み関数ＷＧｊに基づいて重み関数Ｗｉｊを算出する乗算器７１５３と、推定誤差ＥＨＡＴに対し領域ごとに重み付けする乗算器７１５４と、を含んで構成される。
【０１３１】
第１重み関数算出部７１５１は、図９に示すような制御マップを検索することにより、エンジン回転数ＮＥ（ｋ）に応じた第１重み関数の値ＷＮｉ（ｋ）を算出する。
【０１３２】
第２重み関数算出部７１５２は、図１０に示すような制御マップを検索することにより、燃料噴射量ＧＦＵＥＬ（ｋ）に応じた第２重み関数の値ＷＧｊ（ｋ）を算出する。
【０１３３】
乗算器７１５３は、下記式（２４）に示すように、第１重み関数算出部７１５１により算出された第１重み関数の値ＷＮｉ（ｋ）と第２重み関数の値ＷＧｊ（ｋ）との各成分を乗算することにより、重み関数の値Ｗｉｊ（ｋ）を算出する。
【数２４】
2010109667.tif 000026
【０１３４】
乗算器７１５４は、下記式（２５）に示すように、算出された重み関数の値Ｗｉｊ（ｋ）の各成分を推定誤差ＥＨＡＴ（ｋ）に乗算することにより、領域ごとに重み付けされた誤差信号ＷＥＶＮＳｉｊ（ｋ）を算出する。
【数２５】
2010109667.tif 000027
【０１３５】
局所適応入力算出部７１６は、領域ごとに重み付けされた誤差信号ＷＥＶＮＳｉｊが「０」になるように、局所適応入力Ｕｉｊ（ｉ＝１〜４，ｊ＝１〜４）を領域ごとに算出する。
【０１３６】
本実施形態では、誤差信号ＷＥＶＮＳｉｊの収束速度を設定できる応答指定型制御アルゴリズムにより、局所適応入力Ｕｉｊを算出する。この応答指定型制御アルゴリズムとは、偏差の収束挙動を規定した関数に基づいて、偏差の収束速度と収束挙動の両方を指定できる制御アルゴリズムのことをいう。
【０１３７】
局所適応入力算出部７１６は、この応答指定型制御アルゴリズムが実行可能に構成された複数のスライディングモードコントローラを備える。以下では、これらスライディングモードコントローラの動作について説明する。
【０１３８】
先ず、下記式（２６）に示すように、切換関数設定パラメータＰＯＬＥＶと前回制御時の誤差信号ＷＥＶＮＳｉｊ（ｋ−１）との積と、ＷＥＶＮＳｉｊ（ｋ）との和を算出し、これを切換関数σＶ＿ｉｊ（ｋ）として定義する。なお、切換関数設定パラメータＰＯＬＥＶは、所定の設定テーブルに基づいて、−１から０の間で設定されたものが用いられる。
【数２６】
2010109667.tif 000028
【０１３９】
次に、切換関数σＶ＿ｉｊ（ｋ）に基づいて、到達則入力ＵＲＣＨ＿Ｖ＿ｉｊ（ｋ）、及び適応則入力ＵＡＤＰ＿Ｖ＿ｉｊ（ｋ）を算出し、さらに下記式（２７）に示すように、これらＵＲＣＨ（ｋ）、及びＵＡＤＰ（ｋ）の和を算出し、これを局所適応入力Ｕｉｊ（ｋ）として定義する。
【数２７】
2010109667.tif 000029
【０１４０】
到達則入力ＵＲＣＨ＿Ｖ＿ｉｊ（ｋ）は、偏差状態量を切換直線上に載せるための入力であり、下記式（２８）に示すように、切換関数σＶ＿ｉｊ（ｋ）に所定の到達則制御ゲインＫＲＣＨ＿Ｖを乗算することで算出される。
【数２８】
2010109667.tif 000030
【０１４１】
適応則入力ＵＡＤＰ＿Ｖ＿ｉｊ（ｋ）は、モデル化誤差や外乱の影響を抑制し、偏差状態量を切換直線に載せるための入力であり、下記式（２９）に示すように、切換関数σＶ＿ｉｊ（ｋ）と所定の適応則ゲインＫＡＤＰ＿Ｖを乗算したものと、前回制御時の適応則入力ＵＡＤＰ＿Ｖ＿ｉｊ（ｋ−１）との和により算出される。
【数２９】
2010109667.tif 000031
【０１４２】
適応係数算出部７１７は、領域ごとに算出された局所適応入力Ｕｉｊを重み関数Ｗｉｊによる重み結合したものに、「１」を加算することにより、適応入力ＵＶＮＳを算出する。すなわち、適応係数算出部７１７は、下記式（３０）に示すように、局所適応入力Ｕｉｊ（ｋ）と重み関数の値Ｗｉｊとの積の全領域（ｉ＝１〜４，ｊ＝１〜４）に亘る総和に、「１」を加算することにより、適応入力ＵＶＮＳ（ｋ）を算出する。
【数３０】
2010109667.tif 000032
【０１４３】
ここで、上記式（３０）において、「１」を加算した理由は、上述のように適応入力ＵＶＮＳ＝１を基準品の状態とした上で、適応則入力ＵＡＤＰ＿Ｖ＿ｉｊの初期値を「０」とし、かつ、適応入力ＵＶＮＳの初期値を「１」とするためである。なお、適応則入力ＵＡＤＰ＿Ｖ＿ｉｊの初期値を「１」とした場合には、上記式（３０）において「１」を加算する必要は無い。また、適応入力ＵＶＮＳの初期値を、劣化品を示す「０」にする場合にも、上記式（３０）において「１」を加算する必要は無い。
【０１４４】
また、上述のように重み関数Ｗｉｊは、各領域においてのみ「０」でない値を持つ関数ＷＮｉ及びＷＧｊの積に基づいて算出されるものであるため、重み関数の値が「０」になる領域も存在する。したがって、上記式（３０）において、全領域（ｉ＝１〜４，ｊ＝１〜４）に亘る総和を演算する際には、このような重み関数の値Ｗｉｊ（ｋ）が「０」になる領域に関する演算を除外してもよい。これにより、演算負荷を軽減することができる。
【０１４５】
次に、Ｉｎｅｒｔ−ＥＧＲ量制御の具体的な手順について、図１２を参照して説明する。
図１２は、ＥＣＵにより実行されるＩｎｅｒｔ−ＥＧＲ量制御の手順を示すフローチャートである。この処理は、所定の制御周期（例えば、５０ｍｓｅｃ）ごとに実行される。
【０１４６】
ステップＳ１では、ＥＧＲ弁故障フラグＦＥＧＲＮＧが「１」であるか否かを判別する。このＥＧＲ弁故障フラグＦＥＧＲＮＧは、図示しない判定処理において高圧ＥＧＲ弁又は低圧ＥＧＲ弁が故障したと判定されたときに「１」に設定され、それ以外のときには「０」に設定される。この判別がＹＥＳの場合には、ステップＳ１１に移り、Ｉｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＣＭＤを強制的に「０」にした後に、ステップＳ１２に移る。この判別がＮＯの場合には、ステップＳ２に移る。
【０１４７】
ステップＳ２では、ＬＡＦセンサ故障フラグＦＬＡＦＮＧが「１」であるか否かを判別する。このＬＡＦセンサ故障フラグＦＬＡＦＮＧは、図示しない判定処理においてＬＡＦセンサが故障したと判定されたときに「１」に設定され、それ以外のときには「０」に設定される。この判別がＹＥＳの場合には、ステップＳ１１に移る。この判別がＮＯの場合には、ステップＳ３に移る。
【０１４８】
ステップＳ３では、圧力センサ故障フラグＦＣＮＳＮＧが「１」であるか否かを判別する。この圧力センサ故障フラグＦＣＮＳＮＧは、図示しない判定処理において圧力センサの何れかが故障したと判定されたときに「１」に設定され、それ以外のときには「０」に設定される。この判別がＹＥＳの場合には、ステップＳ１１に移る。この判別がＮＯの場合には、ステップＳ３に移る。
【０１４９】
ステップＳ４では、適応バーチャルセンサシステムにより、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴを算出する（上記式（８）〜（１８）参照）。
【０１５０】
ステップＳ５では、ＬＡＦセンサ活性フラグＦＬＡＦＡＣＴが「１」であるか否かを判別する。このＬＡＦセンサ活性フラグＦＬＡＦＡＣＴは、図示しない判定処理においてＬＡＦセンサが活性状態に達したと判定されたときに「１」に設定され、それ以外のときには「０」に設定される。この判別がＹＥＳの場合にはステップＳ６に移り、推定誤差ＥＨＡＴを算出する（上記式（７）参照）。また、この判別がＮＯの場合にはステップＳ７に移り、推定誤差ＥＨＡＴを強制的に「０」にする。
【０１５１】
ステップＳ８では、適応入力ＵＶＮＳを算出する（上記式（２０）〜（３０）参照）。
ステップＳ９では、Ｉｎｅｒｔ−ＥＧＲ量目標値算出部７２により、Ｉｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤを算出する。
【０１５２】
ステップＳ１０では、目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤに対する修正目標値ＩＥＧＲＣＭＤを算出し（上記式（１）〜（６）参照）、ステップＳ１２では、修正目標値ＩＥＧＲＣＭＤに基づいて高圧ＥＧＲバルブリフト量ＬＨＰ及び低圧ＥＧＲバルブリフト量ＬＬＰを算出し、この処理を終了する。
【０１５３】
次に、以上のように構成された本実施形態のＩｎｅｒｔ−ＥＧＲ量制御のシミュレーション結果について、図１３〜図１６を参照して詳述する。
【０１５４】
図１３は、排気浄化システムの高圧ＥＧＲ弁及び低圧ＥＧＲ弁を基準品とし、かつ、適応入力ＵＶＮＳを「１」に固定した場合におけるシミュレーション結果を示す図である。
この場合、各ＥＧＲ弁を基準品としたため、適応入力ＵＶＮＳを強制的に「１」に固定しても、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴは実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴに一致する。このため、実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴを目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤに精度良く制御することができる。したがって、エンジンから排出されるＮＯｘ量は、可能な限り抑制される。
【０１５５】
図１４は、排気浄化システムの高圧ＥＧＲ弁及び低圧ＥＧＲ弁を基準品とし、かつ、適応入力ＵＶＮＳを、非線形適応修正器により算出させた場合におけるシミュレーション結果を示す図である。
この場合、各ＥＧＲ弁を基準品としたため、適応入力ＵＶＮＳを非線形適応器により算出させても、初期値の「１」のまま変動しない。したがって、上述の図１３に示す結果とほぼ同様となり、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴは実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴに一致する。また、実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴを目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤに精度良く制御することができる。したがって、エンジンから排出されるＮＯｘ量は、可能な限り抑制される。
【０１５６】
図１５は、排気浄化システムの高圧ＥＧＲ弁及び低圧ＥＧＲ弁を時刻Ａにおいて基準品から劣化品へ仮想的に変化させ、かつ、適応入力ＵＶＮＳを「１」に固定した場合におけるシミュレーション結果を示す図である。
この場合、時刻Ａ以降、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴと実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴとの間で誤差が発生する。また、ＬＡＦセンサの出力ΦＡＣＴと推定値ΦＨＡＴとの間にも誤差が発生する。しかしながら、適応入力ＵＶＮＳを「１」に固定したため、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴには誤差が生じたままである。このため、Ｉｎｅｒｔ−ＥＧＲ量の目標値ＩＥＧＲＩＤＥＡＬ＿ＣＭＤは、Ｉｎｅｒｔ−ＥＧＲコントローラにより修正されることはない。この結果、エンジンから排出されるＮＯｘ量は増加する。なお、このＮＯｘ量の排出量の増加は、エンジンの高負荷側において特に顕著である。
【０１５７】
図１６は、排気浄化システムの高圧ＥＧＲ弁及び低圧ＥＧＲ弁を時刻Ａにおいて基準品から劣化品へ仮想的に変化させ、かつ、適応入力ＵＶＮＳを非線形適応修正器により算出させた場合におけるシミュレーション結果を示す図である。
この場合、時刻Ａ以降、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴと実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴとの間で誤差が発生する。また、ＬＡＦセンサの出力ΦＡＣＴと推定値ΦＨＡＴとの間にも誤差が発生する。これに伴い、非線形適応修正器は、発生した誤差を最小にするように適応入力ＵＶＮＳを「１」から小さな値へ修正する。また、この適応入力ＵＶＮＳの修正により、Ｉｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴとＬＡＦセンサの出力の推定値ΦＨＡＴに発生した誤差が次第に小さくなる。
Ｉｎｅｒｔ−ＥＧＲコントローラは、実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴの変化を、推定値量ＩＥＧＲＨＡＴを介して検出し、実Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＡＣＴが目標Ｉｎｅｒｔ−ＥＧＲ量ＩＥＧＲＩＤＥＡＬ＿ＣＭＤに収束するように補正係数ＫＥＧＲを修正する。これにより、各ＥＧＲ弁が劣化し、Ｉｎｅｒｔ−ＥＧＲ導入量が低下した場合であっても、これを修正し、適正なＩｎｅｒｔ−ＥＧＲ量に維持できるため、排出されるＮＯｘ量の増加を抑制することができる。
【０１５８】
以上により、本実施形態の適応バーチャルセンサシステムが、排気浄化システムの劣化に対して優れたロバスト性を発揮し、精度良くＩｎｅｒｔ−ＥＧＲ量の推定値ＩＥＧＲＨＡＴを算出できることが確認された。また、この適応バーチャルセンサシステムが、排気浄化システムの固体ばらつきに対しても優れたロバスト性を発揮できることは明らかである。
【０１５９】
本実施形態では、例えば、ＬＡＦセンサ３４により検出手段が構成され、ＥＣＵ７により、第１推定値算出手段、第２推定値算出手段、適応入力算出手段、重み関数設定手段、修正値算出手段、決定手段、及びコントローラが構成される。
【０１６０】
具体的には、例えば、Ｉｎｅｒｔ−ＥＧＲ推定値算出部７１１により第１推定値算出手段が構成され、ＬＡＦセンサ出力推定値算出部７１２により第２推定値算出手段が構成され、非線形適応修正器７１３により適応入力算出手段が構成され、重み関数設定部７１５により重み関数設定手段が構成され、局所適応入力算出部７１６により修正値算出手段が構成され、適応係数算出部７１７により決定手段が構成される。また、例えば、Ｉｎｅｒｔ−ＥＧＲ量目標値算出部７２、Ｉｎｅｒｔ−ＥＧＲコントローラ７３、及びリフト量算出部７４により、コントローラが構成される。
【０１６１】
［第２実施形態］
次に、本発明の第２実施形態を、図面を参照して説明する。
以下の第２実施形態の説明にあたって、第１実施形態と同一構成要件については同一符号を付し、その説明を省略化又は簡略化する。
【０１６２】
図１７は、本実施形態に係るエンジン１及びその排気浄化システム２Ａの構成を示す模式図である。
図１７に示すように、本実施形態は、選択還元触媒６１及びユリア噴射装置６２を備える点と、ＥＣＵ８の構成とが、第１実施形態と異なる。
【０１６３】
排気浄化システム２Ａは、排気管３０のうち酸化触媒３５の下流側に設けられ、この排気管３０を流通する排気中の窒素酸化物（以下、「ＮＯｘ」という）を、還元剤としてのアンモニアの存在下で浄化する選択還元触媒６１と、排気管３０のうち選択還元触媒６１の上流側に、還元剤の元となる尿素水を供給するユリア噴射装置６２と、を備える。
【０１６４】
ユリア噴射装置６２は、ユリアタンク６２１と、ユリア噴射弁６２３とを備える。
ユリアタンク６２１は、尿素水を貯蔵する。ユリア噴射弁６２３は、ＥＣＵ８に接続されており、ＥＣＵ８からの制御信号により動作し、この制御信号に応じた量の尿素水を排気管３０内の酸化触媒３５と選択還元触媒６１との間に噴射する。すなわち、ユリア噴射制御が実行される。
【０１６５】
酸化触媒３５は、排気管３０のうち選択還元触媒６１及びユリア噴射弁６２３よりも上流側に設けられ、排気中のＮＯをＮＯ２に変換し、これにより、選択還元触媒６１におけるＮＯｘの還元を促進する。
【０１６６】
選択還元触媒６１は、アンモニア等の還元剤が存在する雰囲気下で、排気中のＮＯｘを選択的に還元する。具体的には、ユリア噴射装置６２により尿素水を噴射すると、この尿素水は、排気の熱により熱分解又は加水分解されて還元剤としてのアンモニアが生成される。生成されたアンモニアは、選択還元触媒６１に供給され、これらアンモニアにより、排気中のＮＯｘは選択的に還元される。
【０１６７】
ところで、この選択還元触媒６１は、尿素水から生成したアンモニアで排気中のＮＯｘを還元する機能を有するとともに、生成したアンモニアを所定の量だけ貯蔵する機能も有する。以下では、選択還元触媒６１において貯蔵されたアンモニア量をストレージ量とし、選択還元触媒６１において貯蔵できるアンモニア量を最大ストレージ容量とする。
【０１６８】
このようにして貯蔵されたアンモニアは、排気中のＮＯｘの還元にも適宜消費される。このため、ストレージ量が大きくなるに従い、選択還元触媒６１におけるＮＯｘ還元率は高くなる。また、エンジンから排出されたＮＯｘの量に対し尿素水の供給量が少ない場合等には、貯蔵されたアンモニアが、この尿素水の不足分を補うようにしてＮＯｘの還元に消費される。
【０１６９】
ここで、選択還元触媒６１において、最大ストレージ容量を超えてアンモニアが生成された場合、生成されたアンモニアは、選択還元触媒６１の下流側へ排出される。このようにしてアンモニアが選択還元触媒６１に貯蔵されず、その下流側へ排出されることを、以下では「アンモニアスリップ」という。
【０１７０】
後に詳述するように、選択還元触媒６１のストレージ量が所定の目標値に維持されるようにユリア噴射制御を行うことにより、選択還元触媒６１におけるＮＯｘ浄化率を高く維持しながら、アンモニアスリップの発生も極力抑制することができる。
【０１７１】
ＥＣＵ８には、クランク角度位置センサ１１、アクセルセンサ１２、吸気圧力センサ２４、第１排気圧力センサ３２、第２排気圧力センサ３３、ＬＡＦセンサ３４、第１リフトセンサ１３、及び第２リフトセンサ１４の他、酸化触媒温度センサ３７、及び選択還元触媒温度センサ３８が接続されている。
【０１７２】
酸化触媒温度センサ３７は、酸化触媒３５の温度ＴＤＯＣを検出し、検出値に略比例した信号をＥＣＵ８に送信する。選択還元触媒温度センサ３８は、選択還元触媒６１の温度ＴＳＣＲを検出し、検出値に略比例した信号をＥＣＵ８に送信する。
【０１７３】
次に、以上のような排気浄化システム２Ａを制御するＥＣＵを構成するにあたり、本願発明者が着目した課題について説明する。
【０１７４】
このような排気浄化システムを制御する場合、従来では、ユリア噴射装置と選択還元触媒との間に排気中のＮＯｘ量を検出するＮＯｘセンサを設け、このＮＯｘセンサの検出値に基づいてユリア噴射量を決定することが知られている。しかしながら、このようにＮＯｘセンサを用いた場合には、以下のような課題がある。
【０１７５】
（３）先ず、現存するＮＯｘセンサは、高精度でユリア噴射制御を行うにはＮＯｘの観測分解能が十分ではなく、また固体ばらつきが大きい。このため、ユリア噴射量が不足してＮＯｘ浄化率が低下したり、逆にユリア噴射量が過多になってしまい過剰なアンモニアスリップが発生したりする虞がある。
【０１７６】
（４）また、現存するＮＯｘセンサは、高精度でユリア噴射制御を行うには応答性能が十分では無い。このため、特に過渡時において大きなセンシング遅れが生じてしまい、結果としてユリア噴射量が不足しＮＯｘ浄化率が低下する虞がある。
【０１７７】
（５）また、現存するＮＯｘセンサは、センサ素子の割れを防止するために、急速昇温することができない。このため、エンジンの始動後、ＮＯｘセンサが活性に達するまで、数百秒程度かかってしまう場合がある。したがって、この間は、ＮＯｘセンサの出力を利用できないため、ＮＯｘ浄化率が低下したりアンモニアスリップが発生したりする虞がある。
【０１７８】
（６）図１８は、選択還元触媒におけるＮＯｘ浄化率と、ＮＯ２とＮＯｘの比率（ＮＯ２／ＮＯｘ＝ＲＮＯ２）との関係を示す図である。
図１８に示すように、選択還元触媒におけるＮＯｘ浄化率は、ＮＯ２とＮＯｘの比率によって大きく異なる。特に、選択還元触媒の温度が低くなるに従い、ＮＯ２とＮＯｘの比率に対するＮＯｘ浄化率の変動は大きくなる。
選択還元触媒にはこのような特性があるにもかかわらず、現存するＮＯｘセンサではＮＯ２とＮＯｘの比率を検出することは出来ない。したがって、このようなＮＯｘセンサの出力に基づいてユリア噴射量を決定した場合、噴射した尿素水が選択還元触媒におけるＮＯｘ浄化率の向上に寄与しない余剰分がＮＯ２とＮＯｘの比率に応じて発生してしまい、アンモニアスリップが発生するおそれがある。
【０１７９】
以上のように、ユリア噴射制御を行うためにＮＯｘセンサを設けることの利点は多くない。また、ＮＯｘセンサの出力を用いずに、予め設定したマップを用いてユリア噴射量を決定することも考えられるが、この場合も、高い精度でユリア噴射量を決定することは困難であり、また、エンジンや排気浄化システムの劣化や固体ばらつきに対応することも困難である。
【０１８０】
以下では、以上のような課題に鑑みてなされた、排気浄化システム２Ａの制御装置の構成について説明する。以下、詳細に説明するように、本実施形態では、第１実施形態の排気浄化システムにおけるＩｎｅｒｔ−ＥＧＲ量の推定値と同様にして、排気中のＮＯｘ量の推定値を算出する適応バーチャルセンサシステムを構築する。すなわち、本実施形態では、ニューラルネットワークを用いて排気中のＮＯｘ量を推定し、さらにこのＮＯｘ量の推定値に基づいてユリア噴射制御を行う。
【０１８１】
図１９は、排気浄化システム２Ａの制御装置の構成を示すブロック図である。なお、図１９には、排気浄化システム２Ａにおけるユリア噴射制御に係る構成のみを図示する。より具体的には、ユリア噴射装置のユリア噴射量ＧＵＲＥＡの決定に関する、ＥＣＵにより構成されるモジュールのみを図示する。
【０１８２】
このモジュールは、適応バーチャルセンサシステム８１と、フィードフォワード噴射量決定部８２と、ストレージ量目標値設定部８３と、フィードバック噴射量決定部８４と、を含んで構成される。
【０１８３】
このモジュールにおいて、ユリア噴射量ＧＵＲＥＡ（ｋ）は、下記式（３１）に示すように、加算器８５によりフィードフォワード噴射量ＧＵＲＥＡ＿ＦＦ（ｋ）とフィードバック噴射量ＧＵＲＥＡ＿ＳＴ（ｋ）との和を算出することにより決定される。
【数３１】
2010109667.tif 000033
【０１８４】
以下、詳細に説明するように、フィードフォワード噴射量ＧＵＲＥＡ＿ＦＦ（ｋ）は、適応バーチャルセンサシステム８１により算出された排気中のＮＯｘ量の推定値ＮＯＸＨＡＴに基づいて、フィードフォワード噴射量決定部８２により決定される。
また、フィードバック噴射量ＧＵＲＥＡ＿ＳＴ（ｋ）は、ストレージ量目標値設定部８３により設定された目標値ＳＴＵＲＥＡ＿ＣＭＤ（ｋ）にストレージ量を維持するように、フィードバック噴射量決定部８４により決定される。
【０１８５】
適応バーチャルセンサシステム８１は、複数のセンサ２４，３２，３３，３４，１３，１４，３７，３８の検出値Ｐ２，Ｐ３，Ｐ３Ｌ，ΦＡＣＴ，ＬＨＰ＿ＡＣＴ，ＬＬＰ＿ＡＣＴ，ＴＤＯＣ，ＴＳＣＲに基づいて、酸化触媒と選択還元触媒との間の排気のＮＯｘ量の推定値ＮＯＸＨＡＴを算出する。
【０１８６】
図２０は、適応バーチャルセンサシステム８１の構成を示すブロック図である。
適応バーチャルセンサシステム８１は、ＮＯｘ量の推定値ＮＯＸＨＡＴを算出するＮＯｘ量推定値算出部８１１と、ＬＡＦセンサ３４の出力（排気空燃比）の推定値ΦＨＡＴを算出するＬＡＦセンサ出力推定値算出部８１２と、適応入力ＵＶＮＳを算出する非線形適応修正器８１３とを含んで構成される。
【０１８７】
この適応バーチャルセンサシステム８１では、上述の課題（３）〜（５）を解決するため、上述の第１実施形態におけるＩｎｅｒｔ−ＥＧＲ量の推定と同様にしてＮＯｘ量の推定値ＮＯＸＨＡＴを算出する。
【０１８８】
すなわち、本実施形態の適応バーチャルセンサシステム８１では、ニューラルネットワークが構築されたＮＯｘ量推定値算出部８１１により、ＮＯｘ量の推定値ＮＯＸＨＡＴを算出する。さらに、ＮＯｘ量とは別の物理量であり、かつ、ＮＯｘ量と相関のある物理量である排気空燃比の推定値ΦＨＡＴを、ＮＯｘ量推定値算出部８１１と同様のニューラルネットワーク構造を有するＬＡＦセンサ出力推定値算出部８１２で算出する。なお、これらＮＯｘ量推定値算出部８１１及びＬＡＦセンサ出力推定値算出部８１２のニューラルネットワーク構造は、それぞれ、第１実施形態のＩｎｅｒｔ−ＥＧＲ推定値算出部７１１及びＬＡＦセンサ出力推定値算出部７１２のニューラルネットワーク構造とほぼ同様の構成であるので、その詳細な説明を省略する。
【０１８９】
さらに、下記式（３２）に示すように、算出した排気空燃比の推定値ΦＨＡＴ（ｋ）とＬＡＦセンサ３４の検出値ΦＡＣＴ（ｋ）との間の推定誤差ＥＨＡＴ（ｋ）を加算器８１４により算出する。
【数３２】
2010109667.tif 000034
【０１９０】
さらに、非線形適応修正器８１３では、算出した推定誤差ＥＨＡＴが最小になるように、ＮＯｘ量推定値算出部８１１及びＬＡＦセンサ出力推定値算出部８１２に対し共通に入力される適応入力ＵＶＮＳを算出する。なお、この非線形適応修正器８１３は、第１実施形態の非線形適応修正器７１３とほぼ同様の構成であるので、その詳細な説明を省略する。
【０１９１】
次に、ＮＯｘ量を推定するための、ニューラルネットワークの学習について説明する。
先ず、ニューラルネットワークに対する入力ベクトルＵ（ｋ）の成分を、下記式（３３）に示すように定義する。このように、入力ベクトルＵ（ｋ）の成分には、ＮＯｘ量を推定するために必要となる複数の物理量（酸化触媒温度ＴＤＯＣ、選択還元触媒温度ＴＳＣＲ、燃料噴射量ＧＦＵＥＬ、吸気圧力Ｐ２、排気圧力Ｐ３、排気圧力Ｐ３Ｌ、高圧ＥＧＲバルブリフト量の検出値ＬＨＰ＿ＡＣＴ、低圧ＥＧＲバルブリフト量の検出値ＬＬＰ＿ＡＣＴ、エンジン回転数ＮＥ）と、適応入力ＵＶＮＳとが含まれる。また、入力ベクトルの成分には、このように異なる種類の物理量に関するデータが含まれているとともに、異なる時刻の物理量に関するデータも含まれている。また、酸化触媒の下流側のＮＯｘ量は、酸化触媒の酸化能力にも大きく依存するため、入力ベクトルＵの成分には、酸化触媒温度ＴＤＯＣが含まれる。
【数３３】
2010109667.tif 000035
【０１９２】
なお、入力ベクトルＵの成分には、例えば、エンジンの冷却水の温度を含めてもよい。これにより、エンジンの暖機過程におけるＮＯｘ量の推定精度を向上することができる。
【０１９３】
また、このような入力ベクトルＵ（ｋ）に対するニューラルネットワークの出力Ｙ（ｋ）を、下記式（３４）に示すように、ＮＯｘ量の推定値ＮＯＸＨＡＴ（ｋ）として定義する。
【数３４】
2010109667.tif 000036
【０１９４】
ニューラルネットワークの学習は、以上のようにして入力ベクトルＵ（ｋ）と出力Ｙ（ｋ）とを設定した上で、上述の第１実施形態と同様の手順に従って行う。すなわち、排気浄化システムの基準品と劣化品との少なくとも２つを準備し、入力ベクトルＵの成分と、酸化触媒と選択還元触媒との間の排気のＮＯｘ量との関係を示すデータを、準備した排気浄化システムごとに取得する。そして、この取得したデータに基づいて、ニューラルネットワークの学習を行う。
【０１９５】
ここで、ニューラルネットワークの学習に用いる上述のデータの教師データについて説明する。
教師データ、すなわち、ニューラルネットワークの出力ＮＯＸＨＡＴにより再現されるデータには、下記式（３５）に示すように、酸化触媒と選択還元触媒との間のＮＯｘの量を検出するＮＯｘセンサの出力ＮＯＸ（ｋ）に補正係数ＫＭＯＤ＿ＮＯ２（ｋ）を乗算した修正ＮＯｘ量ＮＯＸＭＯＤ（ｋ）のデータを用いる。
【数３５】
2010109667.tif 000037
【０１９６】
ここで、補正係数ＫＭＯＤ＿ＮＯ２（ｋ）は、下記式（３６）及び（３７）に示すように、上述の図１８に示すマップに基づいて決定された選択還元触媒におけるＮＯｘ浄化率ＩＴＡ＿ＮＯＸ（ｋ）を、最大浄化率ＩＴＡ＿ＮＯＸ＿ＭＡＸで割ったものである。なお、ＮＯｘに対するＮＯ２の比率ＲＮＯ２（ｋ）は、フーリエ変換赤外分光光度計（ＦＴＩＲ）などにより測定された値が用いられる。
【数３６】
2010109667.tif 000038
【数３７】
2010109667.tif 000039
【０１９７】
すなわち、上述の修正ＮＯｘ量ＮＯＸＭＯＤは、酸化触媒と選択還元触媒との間の排気のうち、選択還元触媒において浄化可能なＮＯｘ量に相当する。本実施形態では、このような修正ＮＯｘ量が記録された教師データを用いてニューラルネットワークの学習を行い、その出力ＮＯＸＨＡＴを設定する。このようなＮＯｘに対するＮＯ２の比率による影響が予め反映された出力ＮＯＸＨＡＴに基づいてユリア噴射量を決定することにより、余分な尿素水の噴射を抑制し、アンモニアスリップの発生を抑制しながら、選択還元触媒におけるＮＯｘ浄化率を高く維持することができる。したがって、上述の課題（６）を解決することができる。
【０１９８】
また、ＬＡＦセンサ出力推定値算出部８１２のニューラルネットワークに対する学習も、上述のＮＯｘ量推定値算出部８１１のニューラルネットワークと同様の手順により行う。
【０１９９】
フィードフォワード噴射量決定部８２は、下記式（３８）に示すように、適応バーチャルセンサシステム８１により算出されたＮＯｘ量の推定値ＮＯＸＨＡＴ（ｋ）に、変換係数ＫＣＯＮＶ＿ＮＯＸ＿ＵＲＥＡを乗算することにより、フィードフォワード噴射量ＧＵＲＥＡ＿ＦＦ（ｋ）を決定する。下記式（３８）において、変換係数ＫＣＯＮＶ＿ＮＯＸ＿ＵＲＥＡは、ＮＯｘ量からユリア噴射量に変換する変換係数である。より具体的には、変換係数ＫＣＯＮＶ＿ＮＯＸ＿ＵＲＥＡは、所定の量のＮＯｘを還元するために必要なユリア噴射量である。
【数３８】
2010109667.tif 000040
【０２００】
ここで、上述のように、適応バーチャルセンサシステム８１により算出されたＮＯｘ量の推定値ＮＯＸＨＡＴには、ＮＯｘに対するＮＯ２の比率や選択還元触媒温度ＴＳＣＲに応じて変化するＮＯｘ浄化率の影響が反映されている。したがって、上記式（３８）に示すように、推定値に係数を乗算するだけで、過分の無い適切なフィードフォワード噴射量ＧＵＲＥＡ＿ＦＦを決定することができる。
【０２０１】
ストレージ量目標値設定部８３は、選択還元触媒温度の検出値ＴＳＣＲ（ｋ）に基づいて、ストレージ量の目標値ＳＴＵＲＥＡ＿ＣＭＤ（ｋ）を設定する。
【０２０２】
図２１は、最大ストレージ容量ＳＴＵＲＥＡ＿ＭＡＸと選択還元触媒温度ＴＳＣＲとの関係を示す図であり、検出値ＴＳＣＲ（ｋ）に基づいて、ストレージ量の目標値ＳＴＵＲＥＡ＿ＣＭＤ（ｋ）を設定するためのマップを示す図である。
【０２０３】
この図に示すように、選択還元触媒温度ＴＳＣＲが上昇するに従い最大ストレージ容量ＳＴＵＲＥＡ＿ＭＡＸは減少する。そこで、ストレージ量の目標値ＳＴＵＲＥＡ＿ＣＭＤは、アンモニアスリップが発生しないように、最大ストレージ容量ＳＴＵＲＥＡ＿ＭＡＸよりもやや小さな値に設定される。
【０２０４】
図１９に戻って、ストレージ量目標値設定部８３では、図２１に示すようなマップに基づいて、検出値ＴＳＣＲ（ｋ）に応じた目標値ＳＴＵＲＥＡ＿ＣＭＤ（ｋ）を設定する。
【０２０５】
フィードバック噴射量決定部８４は、選択還元触媒の所定のストレージモデルに基づいてストレージ量ＳＴＵＲＥＡを推定しつつ、このストレージ量の推定値ＳＴＵＲＥＡが設定された目標値ＳＴＵＲＥＡ＿ＣＭＤに一致するように、フィードバック噴射量ＧＵＲＥＡ＿ＳＴを決定する。
【０２０６】
図２２は、選択還元触媒のストレージモデルの概念を示す模式図である。
このアンモニアストレージモデルは、選択還元触媒に流入する排気のＮＯｘ量に対するユリア噴射量に応じて、選択還元触媒におけるアンモニアのストレージ量の変化を推定するモデルである。具体的には、選択還元触媒におけるストレージ量の変化の状態を、所定のＮＯｘ量に対してユリア噴射量が適切な状態（図２２の（ａ）参照）と、ユリア噴射量が過剰な状態（図２２の（ｂ）参照）と、ユリア噴射量が不足した状態（図２２の（ｃ）参照）との、３つの状態に分類する。
【０２０７】
図２２の（ａ）に示すように、選択還元触媒に流入するＮＯｘに対して、ユリア噴射量が適切な状態である場合、すなわち、排気中のＮＯｘを最も効率良く還元できるアンモニアの量と、供給した尿素水から生成されるアンモニアの量とが略一致した場合には、ストレージ量の変化はない。
【０２０８】
図２２の（ｂ）に示すように、選択還元触媒に流入するＮＯｘに対して、ユリア噴射量が過剰な状態である場合、すなわち、供給した尿素水から生成されたアンモニアの量が、排気中のＮＯｘを最も効率良く還元できる量より多い場合には、この余剰分のアンモニアが選択還元触媒に貯蔵される。したがって、このような供給過剰（Ｏｖｅｒ−ｄｏｓｉｎｇ）状態では、ストレージ量は増加する。
【０２０９】
図２２の（ｃ）に示すように、選択還元触媒に流入するＮＯｘに対して、ユリア噴射量が不足した状態である場合、すなわち、供給した尿素水から生成されたアンモニアの量が、排気中のＮＯｘを最も効率良く還元できる量より少ない場合には、この不足分は貯蔵されたアンモニアから補われる。したがって、このような供給不足（Ｕｎｄｅｒ−ｄｏｓｉｎｇ）状態では、ストレージ量は減少する。
【０２１０】
フィードバック噴射量決定部８４では、以上のようなストレージモデルに基づいてストレージ量の推定値ＳＴＵＲＥＡを算出する。より具体的には、下記式（３９）〜（４２）に基づいて算出する。
【０２１１】
先ず、選択還元触媒に流入したＮＯｘを還元するために必要な量のユリア噴射量ＧＵＲＥＡ＿ＩＤＥＡＬ（ｋ）は、下記式（３９）に示すように、ＮＯｘ量の推定値ＮＯＸＨＡＴに基づいて算出される。
【数３９】
2010109667.tif 000041
【０２１２】
ストレージ量を増減する要因となるユリア噴射量の余剰分ＤＵＲＥＡ（ｋ）は、下記式（４０）に示すように、実際のユリア噴射量ＧＵＲＥＡ（ｋ）から還元に必要なユリア噴射量ＧＵＲＥＡ＿ＩＤＥＡＬ（ｋ）を減算することにより算出される。
【数４０】
2010109667.tif 000042
【０２１３】
したがって、ストレージ量の推定値ＳＴＵＲＥＡ（ｋ）は、最大ストレージ容量ＳＴＵＲＥＡ＿ＭＡＸ（ｋ）を上限値として、下記式（４１）及び（４２）に示すように、ユリア噴射量の余剰分ＤＵＲＥＡ（ｋ）に基づいて算出される。
【数４１】
2010109667.tif 000043
【数４２】
2010109667.tif 000044
【０２１４】
ここで、最大ストレージ容量ＳＴＵＲＥＡ＿ＭＡＸ（ｋ）は、選択還元触媒温度ＴＳＣＲ（ｋ）に応じて、上述の図２１に示すようなマップを検索することにより設定される。
【０２１５】
フィードバック噴射量決定部８４は、以上のように算出されたストレージ量の推定値ＳＴＵＲＥＡ（ｋ）が目標値ＳＴＵＲＥＡ＿ＣＭＤ（ｋ）に一致するように、下記式（４３）〜（４６）に以下に示すような拡大系Ｉ−Ｐ制御によりフィードバック噴射量ＧＵＲＥＡ＿ＳＴ（ｋ）を決定する。
【０２１６】
先ず、下記式（４３）に示すように、ストレージ量の推定値ＳＴＵＲＥＡ（ｋ）と目標値ＳＴＵＲＥＡ＿ＣＭＤ（ｋ）との偏差ＥＳＴ（ｋ）を算出する。
【数４３】
2010109667.tif 000045
【０２１７】
次に、偏差ＥＳＴ（ｋ）に積分ゲインＫＩＳＴを乗算したものを、下記式（４４）に示すように、積分項ＧＵＲＥＡ＿ＳＴ＿Ｉ（ｋ）として定義する。
【数４４】
2010109667.tif 000046
【０２１８】
一方、ストレージ量の推定値の微分値ＳＴＵＲＥＡ（ｋ）−ＳＴＵＲＥＡ（ｋ−１）を算出し、この微分値に比例ゲインＫＰＳＴを乗算したものを、下記式（４５）に示すように、比例項ＧＵＲＥＡ＿ＳＴ＿Ｐ（ｋ）として定義する。
【数４５】
2010109667.tif 000047
【０２１９】
次に、下記式（４６）に示すように、比例項ＧＵＲＥＡ＿ＳＴ＿Ｐ（ｋ）と積分項ＧＵＲＥＡ＿ＳＴ＿Ｉ（ｋ）の和を算出し、これをフィードバック噴射量ＧＵＲＥＡ＿ＳＴ（ｋ）として決定する。
【数４６】
2010109667.tif 000048
【０２２０】
次に、以上のように構成された本実施形態のユリア噴射制御のシミュレーション結果について、図２３〜図２６を参照して詳述する。
【０２２１】
図２３は、排気浄化システムを基準品とし、かつ、適応入力ＵＶＮＳを「１」に固定した場合におけるシミュレーション結果を示す図である。
この場合、排気浄化システムを基準品としたため、適応入力ＵＶＮＳを強制的に「１」に固定しても、ＮＯｘ量の推定値ＮＯＸＨＡＴは実ＮＯｘ量ＮＯＸＡＣＴに一致する。このため、ストレージ量の推定値ＳＴＵＲＥＡは、目標値ＳＴＵＲＥＡ＿ＣＭＤに維持される。この際、実際のストレージ量も目標値ＳＴＵＲＥＡ＿ＣＭＤに維持されているので、選択還元触媒におけるＮＯｘ浄化率を高く維持することができる。したがって、選択還元触媒の下流側のＮＯｘ量は、可能な限り抑制される。
【０２２２】
図２４は、排気浄化システムを基準品とし、かつ、適応入力ＵＶＮＳを非線形適応器により算出させた場合におけるシミュレーション結果を示す図である。
この場合、排気浄化システムを基準品としたため、適応入力ＵＶＮＳを非線形適応器により算出させても、初期値の「１」のまま変動しない。したがって、上述の図２３に示す結果とほぼ同様となり、ＮＯｘ量の推定値ＮＯＸＨＡＴは実ＮＯｘ量ＮＯＸＡＣＴに一致する。また、ストレージ量の推定値ＳＴＵＲＥＡは、目標値ＳＴＵＲＥＡ＿ＣＭＤに維持される。この際、実際のストレージ量も目標値ＳＴＵＲＥＡ＿ＣＭＤに維持されているので、選択還元触媒におけるＮＯｘ浄化率を高く維持することができる。したがって、選択還元触媒の下流側のＮＯｘ量は、可能な限り抑制される。
【０２２３】
図２５は、排気浄化システムを時刻Ａにおいて基準品から劣化品へ仮想的に変化させ、かつ、適応入力ＵＶＮＳを「１」に固定した場合におけるシミュレーション結果を示す図である。
この場合、時刻Ａ以降、ＮＯｘ量の推定値ＮＯＸＨＡＴと実ＮＯｘ量ＮＯＸＡＣＴとの間で誤差が発生する。また、ＬＡＦセンサの出力ΦＡＣＴと推定値ΦＨＡＴとの間にも誤差が発生する。しかしながら、適応入力ＵＶＮＳを「１」に固定したため、ＮＯｘ量の推定値ＮＯＸＨＡＴには誤差が生じたままである。このため、フィードフォワード噴射量ＧＵＲＥＡ＿ＦＦは、実ＮＯｘ量ＮＯＸＡＣＴに対して不足した状態となり、また、フィードバック噴射量ＧＵＲＥＡ＿ＳＴは、ストレージ量を目標値ＳＴＵＲＥＡ＿ＣＭＤに維持するように適切に決定されなくなる。このため、選択還元触媒におけるストレージ量が急激に減少してしまい、その後、回復することもない。結果として、選択還元触媒におけるＮＯｘ浄化率が低下してしまい、選択還元触媒の下流側のＮＯｘ量が増加してしまう。
【０２２４】
図２６は、排気浄化システムを時刻Ａにおいて基準品から劣化品へ仮想的に変化させ、かつ、適応入力ＵＶＮＳを非線形適応修正器により算出させた場合におけるシミュレーション結果を示す図である。
この場合、時刻Ａ以降、ＮＯｘ量の推定値ＮＯＸＨＡＴと実ＮＯｘ量ＮＯＸＡＣＴとの間で誤差が発生する。また、ＬＡＦセンサの出力ΦＡＣＴと推定値ΦＨＡＴとの間にも誤差が発生する。これに伴い、非線形適応修正器は、発生した誤差を最小にするように適応入力ＵＶＮＳを「１」から修正する。また、この適応入力ＵＶＮＳの修正により、ＮＯｘ量の推定値ＮＯＸＨＡＴとＬＡＦセンサの出力の推定値ΦＨＡＴに発生した誤差が次第に小さくなる。このため、フィードフォワード噴射量ＧＵＲＥＡ＿ＦＦは、実ＮＯｘ量ＮＯＸＡＣＴに対して適切に決定され、選択還元触媒の下流側のＮＯｘ量の増加が抑制される。また、フィードバック噴射量ＧＵＲＥＡ＿ＳＴもストレージ量を目標値ＳＴＵＲＥＡ＿ＣＭＤに維持するように決定される。これにより、ストレージ量は、時刻Ａにおいて、一旦は急激に減少するものの、ストレージ量の目標値ＳＴＵＲＥＡ＿ＣＭＤに向けて次第に回復し始める。
【０２２５】
以上より、本実施形態の適応バーチャルセンサシステムが、排気浄化システムの劣化に対して優れたロバスト性を発揮し、精度良くＮＯｘ量の推定値ＮＯＸＨＡＴを算出できることが確認された。また、この適応バーチャルセンサシステムが、排気浄化システムの固体ばらつきに対しても優れたロバスト性を発揮できることは明らかである。
【０２２６】
本実施形態では、例えば、ＬＡＦセンサ３４により検出手段が構成され、ＥＣＵ８により、第１推定値算出手段、第２推定値算出手段、適応入力算出手段、及びコントローラが構成される。
【０２２７】
具体的には、例えば、ＮＯｘ量推定値算出部８１１により第１推定値算出手段が構成され、ＬＡＦセンサ出力推定値算出部８１２により第２推定値算出手段が構成され、非線形適応修正器８１３により適応入力算出手段が構成される。また、例えば、フィードフォワード噴射量決定部８２、ストレージ量目標値設定部８３、フィードバック噴射量決定部８４、及び加算器８５により、コントローラが構成される。
【０２２８】
［第３実施形態］
次に、本発明の第３実施形態を、図面を参照して説明する。
以下の第３実施形態の説明にあたって、第１実施形態と同一構成要件については同一符号を付し、その説明を省略化又は簡略化する。
【０２２９】
図２７は、本実施形態に係るエンジン１及びその排気浄化システム２Ｂの構成を示す模式図である。
図２７に示すように、本実施形態は、ＮＯｘ浄化触媒６５を備える点と、ＥＣＵ９の構成とが、第１実施形態と異なる。
【０２３０】
排気浄化システム２Ｂは、排気管３０のうち酸化触媒３５の下流側に設けられ、この排気管３０を流通する排気中のＮＯｘを浄化するＮＯｘ浄化触媒６５を備える。
【０２３１】
ＮＯｘ浄化触媒６５は、エンジン１で燃焼する混合気を理論空燃比よりもリーンにしたときに排気中のＮＯｘを吸着又は吸蔵し、還元雰囲気下で吸着又は吸蔵したＮＯｘを還元する。
【０２３２】
ここで、還元雰囲気下とは、エンジン１で燃焼する混合気を理論空燃比よりもリッチにしたときにおける排気の雰囲気下（リッチ燃焼雰囲気下）、又は、還元剤が存在する排気の雰囲気下（還元剤雰囲気下）をいう。また、この還元剤雰囲気は、例えば、エンジン１で燃焼する混合気を理論空燃比よりもリーンにしたまま、排気工程や膨張行程において燃料を噴射（ポスト噴射）したり、排気管３０内に燃料を直接噴射したりすることで生成することができる。
以下では、上述のような方法により、ＮＯｘ浄化触媒６５に流入する排気を還元雰囲気にすることで排気中のＮＯｘを浄化する処理を、還元制御処理という。
【０２３３】
このＮＯｘ浄化触媒６５としては、例えば、アルミナ（Ａｌ２Ｏ３）、セリア（ＣｅＯ２）、及びセリウムと希土類の複合酸化物（以下、「セリア系複合酸化物」という）の担体に担持された、触媒として作用する白金（Ｐｔ）と、ＮＯｘ吸着能力を有するセリアもしくはセリア系複合酸化物と、触媒に生成されたアンモニアアンモニア（ＮＨ３）を、アンモニウムイオン（ＮＨ４＋）として保持する機能を有するゼオライトとを備えるものが用いられる。
【０２３４】
ＥＣＵ９には、クランク角度位置センサ１１、アクセルセンサ１２、吸気圧力センサ２４、第１排気圧力センサ３２、第２排気圧力センサ３３、ＬＡＦセンサ３４、第１リフトセンサ１３、及び第２リフトセンサ１４の他、冷却水温度センサ１５、及びＮＯｘ浄化触媒温度センサ３９が接続されている。
【０２３５】
冷却水温度センサ１５は、エンジン１の冷却水の温度ＴＷを検出し、検出値に略比例した信号をＥＣＵ９に送信する。ＮＯｘ浄化触媒温度センサ３９は、ＮＯｘ浄化触媒６５の温度ＴＬＮＣを検出し、検出値に略比例した信号をＥＣＵ９に送信する。
【０２３６】
このような排気浄化システム２Ｂを制御するＥＣＵを構成するにあたり、本願発明者が着目した課題について説明する。
【０２３７】
図２８は、ＮＯｘ浄化触媒におけるＮＯｘの吸着／吸蔵効率と、ＮＯｘの吸着／吸蔵量との関係を示す図である。
図２８に示すように、ＮＯｘの吸着／吸蔵効率は、ＮＯｘの吸着／吸蔵量が所定の量に達するまではほぼ一定であるが、所定の量を超えると急激に減少する。このため、高いＮＯｘ吸着／吸蔵効率を維持し続けるためには、ＮＯｘの吸着／吸蔵量が過大にならないように適切なタイミングで還元制御処理を行うことが好ましい。
【０２３８】
この場合、ＮＯｘ浄化触媒におけるＮＯｘの吸着／吸蔵量を推定する必要がある。従来では、ＮＯｘ浄化触媒の上流側に排気中のＮＯｘ量を検出するＮＯｘセンサを設け、このＮＯｘセンサの検出値に基づいてＮＯｘの吸着／吸蔵量を推定することが行われている。しかしながら、このようにＮＯｘセンサを用いた場合には、以下のような課題がある。
【０２３９】
（７）先ず、現存するＮＯｘセンサは、高精度でＮＯｘ吸着／吸蔵量を推定するにはＮＯｘの観測分解能が十分ではなく、また固体ばらつきが大きい。このため、ＮＯｘ吸着／吸蔵量の推定値が実際の値からずれてしまい、還元制御処理を実行するタイミングが適切でなくなる場合がある。したがって、ＮＯｘ浄化触媒からＮＯｘが飽和し、ＮＯｘ浄化触媒の下流側へ排出されるＮＯｘの量が増加する虞がある。
【０２４０】
（８）また、現存するＮＯｘセンサは、高精度でＮＯｘ吸着／吸蔵量を推定するには応答性能が十分では無い。このため、特に過渡時において大きなセンシング遅れが生じてしまい、結果として、ＮＯｘ吸着／吸蔵量の推定値が実際の値からずれてしまい、還元制御処理を実行するタイミングが適切でなくなる場合がある。したがって、ＮＯｘ浄化触媒が飽和状態となり、ＮＯｘ浄化触媒の下流側へ排出されるＮＯｘ量が増加する虞がある。
【０２４１】
（９）また、現存するＮＯｘセンサは、センサ素子の割れを防止するために、急速昇温することができない。このため、エンジンの始動後、ＮＯｘセンサが活性に達するまで、数百秒程度かかってしまう場合がある。したがって、この間は、ＮＯｘセンサの出力を利用できないため、ＮＯｘ吸着／吸蔵量を推定することができない。
【０２４２】
（１０）上述の第２実施形態では、還元剤、すなわち尿素水をユリア噴射装置により排気に直接噴射していたため、排気中の還元剤量を推定することは比較的容易である。これに対して本実施形態では、燃料を還元剤として用いるため、排気中の還元剤量は、燃焼室や排気管での燃焼、すなわちエンジンの運転条件により複雑に変化する。このため、従来の方法では還元剤量を推定することは困難である。
【０２４３】
以下では、以上のような課題に鑑みてなされた、排気浄化システム２Ｂの制御装置の構成について説明する。以下、詳細に説明するように、本実施形態では、第１実施形態の排気浄化システムにおけるＩｎｅｒｔ−ＥＧＲ量の推定値と同様にして、排気中のＮＯｘ量の推定値と、排気中の還元剤量の推定値とを算出する適応バーチャルセンサシステムを構築する。すなわち、本実施形態では、ニューラルネットワークを用いて排気中のＮＯｘ量及び還元剤量を推定し、さらにこれらＮＯｘ量及び還元剤量の推定値に基づいて還元制御処理を実行する。
【０２４４】
図２９は、排気浄化システム２Ｂの制御装置の構成を示すブロック図である。なお、図２９には、排気浄化システム２Ｂにおける還元制御処理の実行に係る構成のみを図示する。より具体的には、リッチモードフラグＦＲＩＣＨＭＯＤＥの決定に関する、ＥＣＵにより構成されるモジュールのみを図示する。
【０２４５】
リッチモードフラグＦＲＩＣＨＭＯＤＥは、上述の還元制御処理を実行する時期であることを示すフラグである。すなわち、リッチモードフラグＦＲＩＣＨＭＯＤＥを「１」にすることにより、還元制御処理の実行が指令される。この指令に基づいて、図示しないモジュールにより還元制御処理が実行される。また、このリッチモードフラグＦＲＩＣＨＭＯＤＥを「０」に戻すことにより、還元制御処理の停止が指令される。この指令に基づいて、上述の還元制御処理が停止される。
【０２４６】
このモジュールは、適応バーチャルセンサシステム９１と、ＮＯｘ吸着／吸蔵量推定部９２と、リッチモードコントローラ９３と、を含んで構成される。
【０２４７】
適応バーチャルセンサシステム９１は、複数のセンサ２４，３２，３３，３４，１３，１４，１５の検出値Ｐ２，Ｐ３，Ｐ３Ｌ，ΦＡＣＴ，ＬＨＰ＿ＡＣＴ，ＬＬＰ＿ＡＣＴ，ＴＷ並びにリッチモードフラグＦＲＩＣＨＭＯＤＥに基づいて、ＮＯｘ浄化触媒に流入する排気のＮＯｘ量の推定値ＮＯＸＨＡＴ及びＮＯｘ浄化触媒に流入する排気の還元剤量の推定値ＲＥＤＨＡＴを算出する。
【０２４８】
図３０は、適応バーチャルセンサシステム９１の構成を示すブロック図である。
適応バーチャルセンサシステム９１は、ＮＯｘ量の推定値ＮＯＸＨＡＴを算出するＮＯｘ量推定値算出部９１１と、還元剤量の推定値ＲＥＤＨＡＴを算出する還元剤量推定値算出部９１５と、ＬＡＦセンサ３４の出力（排気空燃比）の推定値ΦＨＡＴを算出するＬＡＦセンサ出力推定値算出部９１２と、適応入力ＵＶＮＳを算出する非線形適応修正器９１３とを含んで構成される。
【０２４９】
この適応バーチャルセンサシステム９１では、上述の課題（７）〜（１０）を解決するため、上述の第１実施形態におけるＩｎｅｒｔ−ＥＧＲ量の推定と同様にして、ＮＯｘ量の推定値ＮＯＸＨＡＴ及び還元剤量の推定値ＲＥＤＨＡＴを算出する。
【０２５０】
すなわち、本実施形態の適応バーチャルセンサシステム９１では、ニューラルネットワークが構築されたＮＯｘ量推定値算出部９１１及び還元剤量推定値算出部９１５により、ＮＯｘ量の推定値ＮＯＸＨＡＴ及び還元剤量の推定値ＲＥＤＨＡＴを算出する。さらに、これらＮＯｘ量及び還元剤量とは別の物理量であり、かつ、ＮＯｘ量及び還元剤量と相関のある物理量である排気空燃比の推定値ΦＨＡＴを、ＮＯｘ量推定値算出部９１１及び還元剤量推定値算出部９１５と同様のニューラルネットワーク構造を有するＬＡＦセンサ出力推定値算出部９１２で算出する。なお、これらＮＯｘ量推定値算出部８１１、還元剤量推定値算出部９１５、及びＬＡＦセンサ出力推定値算出部９１２のニューラルネットワーク構造は、それぞれ、第１実施形態のＩｎｅｒｔ−ＥＧＲ推定値算出部７１１及びＬＡＦセンサ出力推定値算出部７１２のニューラルネットワーク構造とほぼ同様の構成であるので、その詳細な説明を省略する。
【０２５１】
さらに、下記式（４７）に示すように、算出した排気空燃比の推定値ΦＨＡＴ（ｋ）とＬＡＦセンサ３４の検出値ΦＡＣＴ（ｋ）との間の推定誤差ＥＨＡＴ（ｋ）を加算器９１４により算出する。
【数４７】
2010109667.tif 000049
【０２５２】
さらに、非線形適応修正器９１３では、算出した推定誤差ＥＨＡＴが最小になるように、ＮＯｘ量推定値算出部９１１、還元剤量推定値算出部９１５、及びＬＡＦセンサ出力推定値算出部９１２に対し共通に入力される適応入力ＵＶＮＳを算出する。なお、この非線形適応修正器９１３は、第１実施形態の非線形適応修正器７１３とほぼ同様の構成であるので、その詳細な説明を省略する。
【０２５３】
次に、ＮＯｘ量及び還元剤量を推定するための、ニューラルネットワークの学習について説明する。
先ず、ニューラルネットワークに対する入力ベクトルＵ（ｋ）の成分を、下記式（４８）に示すように定義する。このように、入力ベクトルＵ（ｋ）の成分には、ＮＯｘ量及び還元剤量を推定するために必要となる複数の物理量（冷却水温度ＴＷ、リッチモードフラグＦＲＩＣＨＭＯＤＥ、燃料噴射量ＧＦＵＥＬ、吸気圧力Ｐ２、排気圧力Ｐ３、排気圧力Ｐ３Ｌ、高圧ＥＧＲバルブリフト量の検出値ＬＨＰ＿ＡＣＴ、低圧ＥＧＲバルブリフト量の検出値ＬＬＰ＿ＡＣＴ、エンジン回転数ＮＥ）と、適応入力ＵＶＮＳとが含まれる。また、入力ベクトルの成分には、このように異なる種類の物理量に関するデータが含まれているとともに、異なる時刻の物理量に関するデータも含まれている。
【数４８】
2010109667.tif 000050
【０２５４】
また、このような入力ベクトルＵ（ｋ）に対するニューラルネットワークの出力Ｙ（ｋ）を、下記式（４９）及び（５０）に示すように、ＮＯｘ量の推定値ＮＯＸＨＡＴ（ｋ）及び還元剤量の推定値ＲＥＤＨＡＴ（ｋ）として定義する。
【数４９】
2010109667.tif 000051
【数５０】
2010109667.tif 000052
【０２５５】
ニューラルネットワークの学習は、以上のようにして入力ベクトルＵ（ｋ）と出力Ｙ（ｋ）とを設定した上で、上述の第１実施形態と同様の手順に従って行う。すなわち、排気浄化システムの基準品と劣化品との少なくとも２つを準備し、入力ベクトルＵの成分と、排気中のＮＯｘ量及び排気中の還元剤量との関係を示すデータを、準備した排気浄化システムごとに取得する。さらに取得したデータに基づいて、ニューラルネットワークの学習を行う。
【０２５６】
ここで特に、エンジンの冷却水温度ＴＷを変化させた場合のデータも取得しておく。このようなデータに基づいてニューラルネットワークの学習を行うことにより、エンジン始動後の暖機過程での排気中のＮＯｘ量の推定精度を向上し、上述の課題（９）を解決することができる。
【０２５７】
また、ＬＡＦセンサ出力推定値算出部９１２のニューラルネットワークに対する学習も、上述のＮＯｘ量推定値算出部９１１及び還元剤量推定値算出部９１５のニューラルネットワークと同様の手順により行う。
【０２５８】
ＮＯｘ吸着／吸蔵量推定部９２は、適応バーチャルセンサシステム９１により算出されたＮＯｘ量の推定値ＮＯＸＨＡＴ及び還元剤量の推定値ＲＥＤＨＡＴに基づいて、ＮＯｘ浄化触媒におけるＮＯｘの吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴを算出する。
【０２５９】
このＮＯｘ吸着／吸蔵量推定部９２では、下記式（５１）及び（５２）に示すような演算を所定の演算周期ごとに行うことにより、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴ（ｋ）を算出する。なお、下記式（５１）において、ＳＴＮＯＸ＿ＭＡＸ（ｋ）は、ＮＯｘ浄化触媒における最大ＮＯｘ吸着／吸蔵量を示す。
【数５１】
2010109667.tif 000053
【数５２】
2010109667.tif 000054
【０２６０】
図３１は、ＮＯｘ浄化触媒における最大ＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＭＡＸとＮＯｘ浄化触媒温度ＴＬＮＣとの関係を示す図である。
この図に示すように、最大ＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＭＡＸは、ＮＯｘ浄化触媒温度ＴＬＮＣに応じて変化する。上記式（５１）において、最大ＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＭＡＸ（ｋ）は、ＮＯｘ浄化触媒温度ＴＬＮＣ（ｋ）に基づいて、この図に示すようなマップを検索することにより算出される。
【０２６１】
また、上記式（５２）において、右辺第２項は、排気中のＮＯｘを吸着／吸蔵したことによる増量分を示し、右辺第３項は、ＮＯｘを還元したことによる減量分を示す。
【０２６２】
すなわち、右辺第２項のＩＴＡ＿ＮＯＸ（ｋ）は、ＮＯｘ浄化触媒のＮＯｘ吸着／吸蔵効率の推定値を示す。このＮＯｘ吸着／吸蔵効率の推定値ＩＴＡ＿ＮＯＸ（ｋ）は、例えば、ＮＯｘ浄化触媒温度ＴＬＮＣ（ｋ）及びＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴ（ｋ）に基づいて、上述の図２８に示すようなマップを検索することにより算出される。
また、右辺第３項の係数ＫＲＥＤは、還元剤量からＮＯｘ量へ変換する変換係数である。より具体的には、所定の量の還元剤で還元化能なＮＯｘの量を示す。
【０２６３】
リッチモードコントローラ９３は、２つの閾値ＳＴＲＩＣＨ＿ＯＮ（ｋ）及びＳＴＲＩＣＨ＿ＯＦＦ（ｋ）と、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴ（ｋ）との比較に基づいて、リッチモードフラグＦＲＩＣＨＭＯＤＥ（ｋ）を設定する。
【０２６４】
より具体的には、リッチモードコントローラ９３は、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴ（ｋ）がＳＴＲＩＣＨ＿ＯＮ（ｋ）を超えた場合に、リッチモードフラグＦＲＩＣＨＭＯＤＥ（ｋ）を「０」から「１」にする。また、リッチモードコントローラ９３は、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴ（ｋ）がＳＴＲＩＣＨ＿ＯＦＦ（ｋ）を下回った場合に、リッチモードフラグＦＲＩＣＨＭＯＤＥ（ｋ）を「１」から「０」に戻す。
【０２６５】
また、これら閾値ＳＴＲＩＣＨ＿ＯＮ（ｋ）及びＳＴＲＩＣＨ＿ＯＦＦ（ｋ）は、それぞれ、下記式（５３）及び（５４）に示すように、図３１に示すマップを検索して設定された最大ＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＭＡＸ（ｋ）に係数ＫＲＯＮ及びＫＲＯＦＦを乗算することで算出される。ここで、係数ＫＲＯＮは例えば０．８０に設定され、係数ＫＲＯＦＦは例えば０．１５に設定される。
【数５３】
2010109667.tif 000055
【数５４】
2010109667.tif 000056
【０２６６】
次に、以上のように構成された本実施形態の還元制御処理のシミュレーション結果について、図３２〜図３５を参照して詳述する。
【０２６７】
図３２は、排気浄化システムを基準品とし、かつ、適応入力ＵＶＮＳを「１」に固定した場合におけるシミュレーション結果を示す図である。
この場合、排気浄化システムを基準品としたため、適応入力ＵＶＮＳを強制的に「１」に固定しても、ＮＯｘ量の推定値ＮＯＸＨＡＴは実ＮＯｘ量ＮＯＸＡＣＴに一致する。このため、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴは、実際のＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＡＣＴに一致する。このため、還元制御処理を実行するタイミングは適切に維持される。したがって、ＮＯｘ浄化触媒の下流側のＮＯｘ量は、可能な限り抑制される。
【０２６８】
図３３は、排気浄化システムを基準品とし、かつ、適応入力ＵＶＮＳを非線形適応器により算出させた場合におけるシミュレーション結果を示す図である。
この場合、排気浄化システムを基準品としたため、適応入力ＵＶＮＳを非線形適応器により算出させても、初期値の「１」のまま変動しない。したがって、上述の図３２に示す結果とほぼ同様となり、ＮＯｘ量の推定値ＮＯＸＨＡＴは実ＮＯｘ量ＮＯＸＡＣＴに一致する。また、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴは、実際のＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＡＣＴに一致する。このため、還元制御処理を実行するタイミングは適切に維持される。したがって、ＮＯｘ浄化触媒の下流側のＮＯｘ量は、可能な限り抑制される。
【０２６９】
図３４は、排気浄化システムを時刻Ａにおいて基準品から劣化品へ仮想的に変化させ、かつ、適応入力ＵＶＮＳを「１」に固定した場合におけるシミュレーション結果を示す図である。
この場合、時刻Ａ以降、ＮＯｘ量の推定値ＮＯＸＨＡＴと実ＮＯｘ量ＮＯＸＡＣＴとの間で誤差が発生する。また、ＬＡＦセンサの出力ΦＡＣＴと推定値ΦＨＡＴとの間にも誤差が発生する。しかしながら、適応入力ＵＶＮＳを「１」に固定したため、ＮＯｘ量の推定値ＮＯＸＨＡＴには誤差が生じたままである。このため、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴと、実際のＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＡＣＴとの間の誤差は次第に大きくなる。このため、還元制御処理を実行するタイミングが適切なタイミングからずれ始める。すなわち、実際のＮＯｘ吸着／吸蔵量ＳＴＮＯＸ＿ＡＣＴが閾値ＳＴＲＩＣＨ＿ＯＮを上回っているにもかかわらず、還元制御処理が実行されなくなってしまう。このため、ＮＯｘ浄化触媒が飽和状態になってしまい、結果として、ＮＯｘ浄化触媒の下流側のＮＯｘ量が増加してしまう。
【０２７０】
図３５は、排気浄化システムを時刻Ａにおいて基準品から劣化品へ仮想的に変化させ、かつ、適応入力ＵＶＮＳを非線形適応修正器により算出させた場合におけるシミュレーション結果を示す図である。
この場合、時刻Ａ以降、ＮＯｘ量の推定値ＮＯＸＨＡＴと実ＮＯｘ量ＮＯＸＡＣＴとの間で誤差が発生する。また、ＬＡＦセンサの出力ΦＡＣＴと推定値ΦＨＡＴとの間にも誤差が発生する。これに伴い、非線形適応修正器は、発生した誤差を最小にするとように適応入力ＵＶＮＳを「１」から修正する。また、この適応入力ＵＶＮＳの修正により、ＮＯｘ量の推定値ＮＯＸＨＡＴとＬＡＦセンサの出力の推定値ΦＨＡＴに発生した誤差が次第に小さくなる。このため、ＮＯｘ吸着／吸蔵量の推定値ＳＴＮＯＸ＿ＨＡＴが実際のＮＯｘ吸着／吸蔵量から大きくずれるのを抑制することができる。したがって、ＮＯｘ浄化触媒が飽和状態になるのを防止できる、また、これにより、ＮＯｘ浄化触媒の下流側のＮＯｘ量の増加を抑制することができる。
【０２７１】
以上より、本実施形態の適応バーチャルセンサシステムが、排気浄化システムの劣化に対して優れたロバスト性を発揮し、精度良くＮＯｘ量の推定値ＮＯＸＨＡＴを算出できることが確認された。また、この適応バーチャルセンサシステムが、排気浄化システムの固体ばらつきに対しても優れたロバスト性を発揮できることは明らかである。
【０２７２】
本実施形態では、例えば、ＬＡＦセンサ３４により検出手段が構成され、ＥＣＵ９により、第１推定値算出手段、第２推定値算出手段、適応入力算出手段、及びコントローラが構成される。
【０２７３】
具体的には、例えば、ＮＯｘ量推定値算出部９１１及び還元剤量推定値算出部９１５により第１推定値算出手段が構成され、ＬＡＦセンサ出力推定値算出部９１２により第２推定値算出手段が構成され、非線形適応修正器９１３により適応入力算出手段が構成される。また、例えば、ＮＯｘ吸着／吸蔵量推定部９２、及びリッチモードコントローラ９３により、コントローラが構成される。
【０２７４】
なお本発明は上述した実施形態に限るものではなく、種々の変形が可能である。
【０２７５】
例えば、上記実施形態では、ニューラルネットワークとして、図６に示すような階層型のパーセプトロンを用いたが、これに限らない。ニューラルネットワーク構造としては、この他、網目状のボルツマンマシン、出力の過去値を入力として用いるリカレント型ニューラルネットワーク、入力に各データの時系列データを用いるカオス型ニューラルネットワーク、多くの階層やニューロンを用いずに高精度の推定が可能な自己組織化マップ（ＳＯＭ）、並びに、ＲＢＦニューラルネットワークを用いてもよい。
【０２７６】
上記実施形態では、非線形適応修正器により、推定誤差及び参照パラメータに応じた適応入力を算出した（上述の図８〜図１１参照）が、これに限らない。例えば、スライディングモード制御、バックステッピング制御、ＰＩＤ制御、及び最適制御アルゴリズムなどの既知のフィードバックアルゴリズムを用いて適応入力を算出することも効果的である。
【０２７７】
上記実施形態の非線形適応修正器では、２つの参照パラメータ（エンジン回転数、燃料噴射量）を基底とした２次元の空間に重み関数を定義し、この重み関数を用いて適応入力ＵＶＮＳを算出したが、これに限らない。すなわち、重み関数を定義する空間の次元数、及び、参照パラメータとして用いる物理量の種類は、上記実施形態に示した例に限らない。空間の次元数は、例えば、１次元又は３次元以上であってもよい。また、参照パラメータとして用いる物理量には、例えば、吸気圧力Ｐ２、排気圧力Ｐ３，Ｐ３Ｌ、及びＥＧＲバルブリフト量ＬＬＰ，ＬＨＰなどを用いてもよい。
【０２７８】
上記実施形態の非線形適応修正器では、スライディングモード制御アルゴリズムを用いて領域ごとの局所適応入力Ｕｉｊを算出した（上記式（２６）〜（３０）参照）が、これに限らない。例えば、バックステッピング制御、ＰＩＤ制御、及び最適制御アルゴリズムなどの既知のフィードバックアルゴリズムを用いて領域ごとの局所適応入力Ｕｉｊを算出してもよい。
【０２７９】
上記第１実施形態では、Ｉｎｅｒｔ−ＥＧＲ量を推定したが、これに限らない。全体のＥＧＲ量のうち不活性分を示すＩｎｅｒｔ−ＥＧＲ量だけでなく、全体のＥＧＲ量もセンサを用いて検出することが難しい物理量である。このため、ＥＧＲ量を適応バーチャルセンサシステムにより推定し、この推定に基づいて各ＥＧＲ弁の操作量を決定することも効果的である。
【０２８０】
上記第２実施形態の適応バーチャルセンサシステムでは、排気のＮＯｘに対するＮＯ２の比率に応じて選択還元触媒におけるＮＯｘ浄化率が変化することを考慮して、選択還元触媒に流入する排気のうち選択還元触媒で浄化可能なＮＯｘ量の推定値をＮＯＸＨＡＴとして算出したが、これに限らない。すなわち、ＮＯｘに対するＮＯ２の比率による影響を、ニューラルネットワークの出力の段階で考慮したが、これに限らない。
例えば、排気中のＮＯｘ量の推定値を算出するニューラルネットワークと、同じ排気中のＮＯ２量の推定値を算出するニューラルネットワークとの２つを準備し、これらの出力に基づいて、ユリア噴射量ＧＵＲＥＡを決定する段階において、ＮＯｘに対するＮＯ２の比率による影響を考慮してもよい。
【０２８１】
上記第２実施形態では、アンモニアを還元剤とし、かつ、この還元剤の元となる添加剤として尿素水を供給する尿素添加式の排気浄化システムに、本発明を適用した例を示したが、これに限るものではない。
例えば、尿素水を供給しこの尿素水からアンモニアを生成せずに、直接アンモニアを供給してもよい。また、アンモニアの元となる添加剤としては、尿素水に限らず他の添加剤を用いてもよい。また、ＮＯｘを還元するための還元剤はアンモニアに限るものではない。本発明は、ＮＯｘを還元するための還元剤として、アンモニアの代わりに、例えば炭化水素を用いた排気浄化システムに適用することもできる。
【図１】
2010109667.tif 000057
【図２】
2010109667.tif 000058
【図３】
2010109667.tif 000059
【図４】
2010109667.tif 000060
【図５】
2010109667.tif 000061
【図６】
2010109667.tif 000062
【図７】
2010109667.tif 000063
【図８】
2010109667.tif 000064
【図９】
2010109667.tif 000065
【図１０】
2010109667.tif 000066
【図１１】
2010109667.tif 000067
【図１２】
2010109667.tif 000068
【図１３】
2010109667.tif 000069
【図１４】
2010109667.tif 000070
【図１５】
2010109667.tif 000071
【図１６】
2010109667.tif 000072
【図１７】
2010109667.tif 000073
【図１８】
2010109667.tif 000074
【図１９】
2010109667.tif 000075
【図２０】
2010109667.tif 000076
【図２１】
2010109667.tif 000077
【図２２】
2010109667.tif 000078
【図２３】
2010109667.tif 000079
【図２４】
2010109667.tif 000080
【図２５】
2010109667.tif 000081
【図２６】
2010109667.tif 000082
【図２７】
2010109667.tif 000083
【図２８】
2010109667.tif 000084
【図２９】
2010109667.tif 000085
【図３０】
2010109667.tif 000086
【図３１】
2010109667.tif 000087
【図３２】
2010109667.tif 000088
【図３３】
2010109667.tif 000089
【図３４】
2010109667.tif 000090
【図３５】
2010109667.tif 000091
【国際調査報告】
2010109667.tif 000092
2010109667.tif 000093