(57)【要約】  【課題】  隣接する燃焼気筒の間隔による影響を気筒別空燃比の推定値に反映させる。  【解決手段】  気筒別空燃比推定モデルは、隣接する燃焼気筒の間隔（燃焼間隔）による影響を気筒別空燃比の推定値に反映させるために、隣接する燃焼気筒から排出されるガスの混合と、混合したガスが空燃比センサ１６の位置に到達するまでの移動とを考慮するように設計されている。隣接する燃焼気筒から排出されるガスの混合は、各気筒の排気マニホールド１２が合流する合流部２２で発生し、この合流部２２で混合したガスは空燃比センサ１６の位置まで流れて空燃比が検出される。隣接する燃焼気筒から排出されるガスの混合を評価する際に、各燃焼気筒の排気マニホールド１２の長さと形状を考慮する。混合ガスが空燃比センサ１６の位置に到達するまでの移動を評価する際に合流部２２から空燃比センサ１６の位置までの距離もしくは容積と各燃焼気筒の排出ガス量とを考慮する。  【選択図】  図１
2005207354.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の複数気筒の排気マニホールドが接続された集合排気管に、各気筒から排出されて合流部にて混合したガスの空燃比を検出する空燃比センサを設置し、この空燃比センサで検出した混合ガスの空燃比に基づいて当該ガスを排出した各気筒の空燃比（以下「気筒別空燃比」という）を推定する気筒別空燃比推定手段を備えた内燃機関の気筒別空燃比推定装置において、
前記気筒別空燃比推定手段は、隣接する燃焼気筒の間隔（以下「燃焼間隔」という）による影響を気筒別空燃比の推定値に反映させることを特徴とする内燃機関の気筒別空燃比推定装置。
【請求項２】
前記気筒別空燃比推定手段は、隣接する燃焼気筒から排出されるガスの混合と、混合したガスが前記空燃比センサの位置に到達するまでの移動とを前記燃焼間隔による影響として考慮することを特徴とする請求項１に記載の内燃機関の気筒別空燃比推定装置。
【請求項３】
前記気筒別空燃比推定手段は、前記隣接する燃焼気筒から排出されるガスの混合を評価する際に、各燃焼気筒の排気マニホールドの長さを考慮することを特徴とする請求項２に記載の内燃機関の気筒別空燃比推定装置。
【請求項４】
前記気筒別空燃比推定手段は、前記隣接する燃焼気筒から排出されるガスの混合を評価する際に、各気筒の排気マニホールドの形状を考慮することを特徴とする請求項２又は３に記載の内燃機関の気筒別空燃比推定装置。
【請求項５】
前記気筒別空燃比推定手段は、前記混合したガスが前記空燃比センサの位置に到達するまでの移動を評価する際に、各燃焼気筒のガスの合流部から前記空燃比センサの位置までの距離もしくは容積と各燃焼気筒の排出ガス量とを考慮することを特徴とする請求項２乃至４のいずれかに記載の内燃機関の気筒別空燃比推定装置。
【請求項６】
複数の気筒グループからなる内燃機関であって、各気筒グループ毎にそれぞれ複数気筒の排気マニホールドが接続された集合排気管が設けられていると共に、各集合排気管にそれぞれ前記空燃比センサが設置され、
前記気筒別空燃比推定手段は、各気筒グループ毎に燃焼間隔による影響を気筒別空燃比の推定値に反映させることを特徴とする請求項１乃至５のいずれかに記載の内燃機関の気筒別空燃比推定装置。
【請求項７】
請求項１乃至６のいずれかに記載の内燃機関の気筒別空燃比推定装置と、この気筒別空燃比推定装置により推定した気筒別空燃比の気筒間ばらつきを小さくする方向に各気筒の空燃比を制御する気筒別空燃比制御手段とを備えていることを特徴とする内燃機関の気筒別空燃比制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の複数気筒の排気マニホールドが接続された集合排気管に設置した空燃比センサで検出したガスの空燃比に基づいて当該ガスを排出した気筒の空燃比を推定する内燃機関の気筒別空燃比推定装置及び気筒別空燃比制御装置に関する発明である。
【背景技術】
【０００２】
近年、内燃機関の気筒間の空燃比ばらつきを少なくして空燃比制御精度を向上させるために、特許文献１（特許第３２１７６８０号公報）に記載されているように、内燃機関の排気系の挙動を記述するモデルを設定して、集合排気管に設置した１つの空燃比センサの検出値（集合排気管を流れるガスの空燃比）を該モデルに入力し、その内部状態を観測するオブザーバによって各気筒の空燃比（気筒別空燃比）を推定すると共に、その推定値に基づいて各気筒の空燃比を目標値にフィードバック制御するようにしたものがある。
【特許文献１】特許第２７１７７４４号公報（第１頁〜第２頁等）
【発明の開示】
【発明が解決しようとする課題】
【０００３】
例えば、Ｖ型エンジンのように、複数のバンク（気筒グループ）からなる内燃機関においては、各バンク毎にそれぞれ集合排気管を設け、各集合排気管にそれぞれ空燃比センサを設置したものがある。この構成では、各バンク毎に空燃比センサの検出値に基づいて気筒別空燃比を推定することになるが、１つのバンクに設けられた複数の気筒の燃焼間隔（排気行程の間隔）は等間隔とはならない。この理由を、Ｖ型８気筒エンジンを例にして説明する。Ｖ型８気筒エンジンは、２つのバンクからなり、各バンクにそれぞれ４気筒ずつ設けられている。エンジン全体（８気筒全体）で見れば、燃焼間隔は等間隔（９０℃Ａ間隔）であるが、図２に示すように、片方のバンクの４つの気筒＃１，＃３，＃５，＃７についてのみ見れば、燃焼間隔（排気行程の間隔）が９０℃Ａ、１８０℃Ａ、２７０℃Ａの３通りに変化するため、燃焼間隔が不等間隔になる。燃焼間隔が長い場合（２７０℃Ａの場合）は、空燃比センサの位置に到達するガスの中に、他の燃焼気筒から排出されるガスが混じっていないが、燃焼間隔が短い場合（９０℃Ａの場合）は、空燃比センサの位置に到達するガスの中に他の燃焼気筒から排出されるガスが混じり込んで空燃比が変化しているものと思われる。
【０００４】
しかし、従来の気筒別空燃比推定モデルは、排気系が１系統のエンジンのように、燃焼間隔が等間隔になるエンジンの排気系の挙動をモデル化したものであるため、このモデルを燃焼間隔が不等間隔になるＶ型８気筒エンジン等に適用しても、気筒別空燃比を精度良く推定できないという問題があった。
【０００５】
そこで、本発明の第１の目的は、燃焼間隔が等間隔、不等間隔のいずれの場合でも、気筒別空燃比を精度良く推定できる内燃機関の気筒別空燃比推定装置を提供することであり、また、本発明の第２の目的は、燃焼間隔が等間隔、不等間隔のいずれの場合でも、気筒別空燃比制御を精度良く実施することができる内燃機関の気筒別空燃比制御装置を提供することである。
【課題を解決するための手段】
【０００６】
上記第１の目的を達成するために、請求項１に係る発明は、内燃機関の複数気筒の排気マニホールドが接続された集合排気管に設置した空燃比センサで検出した混合ガスの空燃比に基づいて当該ガスを排出した各気筒の空燃比（以下「気筒別空燃比」という）を推定する気筒別空燃比推定手段を備えた内燃機関の気筒別空燃比推定装置において、前記気筒別空燃比推定手段は、隣接する燃焼気筒の間隔（以下「燃焼間隔」という）による影響を気筒別空燃比の推定値に反映させるようにしたことを特徴とするものである。このようにすれば、燃焼間隔が等間隔、不等間隔のいずれの場合でも、気筒別空燃比を精度良く推定することができる。
【０００７】
この場合、請求項２のように、隣接する燃焼気筒から排出されるガスの混合と、混合したガスが空燃比センサの位置に到達するまでの移動とを前記燃焼間隔による影響として考慮するようにすると良い。例えば、燃焼間隔が短くなるほど、隣接する燃焼気筒から排出されるガスの混合（重なり）の度合いが大きくなり、先の燃焼気筒のガスの空燃比が後の燃焼気筒のガスの空燃比の方向へ変化する度合いが大きくなる。更に、隣接する燃焼気筒から排出されるガスの混合（重なり）の度合いが大きくなるほど、集合排気管に流入するガス量が増加してガスの流速が速くなり、ガスが空燃比センサの位置に到達するまでの時間が短くなる。従って、隣接する燃焼気筒から排出されるガスの混合と、混合したガスが空燃比センサの位置に到達するまでの移動とを評価すれば、燃焼間隔による影響を精度良く評価することができる。
【０００８】
更に、請求項３のように、隣接する燃焼気筒から排出されるガスの混合を評価する際に、各燃焼気筒の排気マニホールドの長さを考慮するようにすると良い。このようにすれば、不等長の排気マニホールドを持つ内燃機関でも、隣接する燃焼気筒から排出されるガスの混合を精度良く評価することができる。
【０００９】
また、請求項４のように、隣接する燃焼気筒から排出されるガスの混合を評価する際に、各気筒の排気マニホールドの形状を考慮するようにすると良い。このようにすれば、排気マニホールドの形状がガス当り（空燃比センサ周辺のガス流動挙動）に及ぼす影響を考慮してガスの混合を精度良く評価することができる。
【００１０】
また、請求項５のように、混合したガスが空燃比センサの位置に到達するまでの移動を評価する際に、各燃焼気筒のガスの合流部から空燃比センサの位置までの距離もしくは容積と各燃焼気筒の排出ガス量とを考慮するようにすると良い。このようにすれば、各燃焼気筒のガスの合流部から空燃比センサの位置までガスが流動するのに要する時間を精度良く判定することができ、それによって、空燃比センサの検出値がいずれの時期に排出されたガスの空燃比であるかを精度良く特定することができる。
【００１１】
以上説明した請求項１乃至５に係る発明は、複数の気筒グループからなる内燃機関の各気筒グループ毎にそれぞれ集合排気管と空燃比センサが設けられた構成のものに適用すると良い。複数の気筒グループからなる内燃機関では、１つの気筒グループについてのみ見れば、燃焼間隔が不等間隔になるため、従来の気筒別空燃比推定方法では、気筒別空燃比を精度良く推定できないが、本発明を適用すれば、燃焼間隔が不等間隔になる場合でも、気筒別空燃比を精度良く推定できる。
【００１２】
また、前記第２の目的を達成するために、請求項７のように、請求項１乃至６のいずれかに記載の内燃機関の気筒別空燃比推定装置と、この気筒別空燃比推定装置により推定した気筒別空燃比の気筒間ばらつきを小さくする方向に各気筒の空燃比を制御する気筒別空燃比制御手段を備える構成としても良い。このようにすれば、燃焼間隔が等間隔、不等間隔のいずれの場合でも、気筒別空燃比制御を精度良く実施することができる。
【発明を実施するための最良の形態】
【００１３】
以下、本発明を例えばＶ型８気筒エンジンに適用した一実施例を説明する。
まず、図１に基づいてＶ型８気筒エンジンの排気系の構成を説明する。Ｖ型８気筒エンジン１１は、２つのバンク（ＡバンクとＢバンク）をＶ字型に配置し、Ａバンクには、４つの気筒＃１，＃３，＃５，＃７を直列に配置し、Ｂバンクには、残りの４つの気筒＃２，＃４，＃６，＃８を直列に配置した構成となっている。ＡバンクとＢバンクには、それぞれ別々の排気系が構成され、各バンクの４本の排気マニホールド１２は、それぞれ別々の集合排気管１４に接続されている。各バンクの集合排気管１４には、それぞれ排出ガスの空燃比を検出する空燃比センサ１６が設置され、各空燃比センサ１６の下流側に排出ガス浄化用の触媒１８が設置されている。
【００１４】
上記空燃比センサ１６等の各種センサの出力は、エンジン制御回路（ＥＣＵ）２０に入力される。このＥＣＵ２０は、マイクロコンピュータを主体として構成され、内蔵されたＲＯＭ（記憶媒体）に記憶された各種のエンジン制御プログラムを実行することで、エンジン運転状態に応じて各気筒の燃料噴射量や点火時期を制御する。
【００１５】
本実施例では、ＥＣＵ２０は、後述する気筒別空燃比制御用の各ルーチンを実行することで、後述する気筒別空燃比推定モデルを用いて各バンクの空燃比センサ１６の検出値（各バンクの集合排気管１４を流れる排出ガスの実空燃比）に基づいて各バンクの気筒別空燃比（各気筒毎の空燃比）を推定し、各バンク毎に気筒別空燃比推定値の平均値を算出して、その平均値を基準空燃比（各バンクの目標空燃比）に設定する。そして、各バンク毎に気筒別空燃比推定値と基準空燃比との偏差を各気筒毎に算出して、その偏差（気筒間の空燃比ばらつき）が小さくなるように気筒別補正量（各気筒の燃料補正量）を算出し、その算出結果に基づいて気筒別燃料噴射量を補正することで、各気筒に供給する混合気の空燃比を各気筒毎に補正して気筒間の空燃比ばらつきを少なくするように制御する（以下、この制御を気筒別空燃比制御という）。
【００１６】
ここで、各バンクの空燃比センサ１６の検出値（各バンクの集合排気管１４を流れる排出ガスの実空燃比）に基づいて各バンクの気筒別空燃比を推定するモデル（以下「気筒別空燃比推定モデル」という）の具体例を説明する。
【００１７】
Ｖ型８気筒エンジン１１は、エンジン全体（８気筒全体）で見れば、隣接する燃焼気筒の間隔（以下「燃焼間隔」という）は、等間隔（９０℃Ａ間隔）であるが、図２に示すように、片方のバンク（Ａバンク）の４つの気筒＃１，＃３，＃５，＃７についてのみ見れば、燃焼間隔（排気行程の間隔）が９０℃Ａ、１８０℃Ａ、２７０℃Ａの３通りに変化するため、燃焼間隔が不等間隔になる。燃焼間隔が長い場合（２７０℃Ａの場合）は、空燃比センサ１６の位置に到達するガスの中に、他の燃焼気筒から排出されるガスが混じっていないが、燃焼間隔が短い場合（９０℃Ａの場合）は、空燃比センサ１６の位置に到達するガスの中に他の燃焼気筒から排出されるガスが混じり込んで空燃比が変化しているものと思われる。
【００１８】
そこで、本実施例の気筒別空燃比推定モデルは、燃焼間隔による影響を気筒別空燃比の推定値に反映させるために、隣接する燃焼気筒から排出されるガスの混合（重なり）と、混合したガスが空燃比センサ１６の位置に到達するまでの移動とを考慮するように設計されている。隣接する燃焼気筒から排出されるガスの混合は、各バンクの４本の排気マニホールド１２が合流する合流部２２で発生し、この合流部２２で混合したガスは、空燃比センサ１６の位置まで流れて空燃比が検出されることになる。
【００１９】
例えば、燃焼間隔が短くなるほど、隣接する燃焼気筒から排出されるガスの混合（重なり）の度合いが大きくなり、先の燃焼気筒のガスの空燃比が後の燃焼気筒のガスの空燃比の方向へ変化する度合いが大きくなる。更に、隣接する燃焼気筒から排出されるガスの混合（重なり）の度合いが大きくなるほど、集合排気管１４を流れるガス量が増加してガスの流速が速くなり、混合ガスが合流部２２から空燃比センサ１６の位置に到達するまでの時間が短くなる。従って、隣接する燃焼気筒から排出されるガスの混合と、混合したガスが合流部２２から空燃比センサ１６の位置に到達するまでの移動とを評価すれば、燃焼間隔による影響を精度良く評価することができる。
【００２０】
隣接する燃焼気筒から排出されるガスの混合を評価する際に、各燃焼気筒の排気マニホールド１２の長さと形状を考慮するものとする。排気マニホールド１２の長さを考慮すれば、不等長の排気マニホールド１２を持つエンジンでも、隣接する燃焼気筒から排出されるガスの混合を精度良く評価することができる。また、排気マニホールド１２の形状を考慮すれば、排気マニホールド１２の形状がガス当り（空燃比センサ１６周辺のガス流動挙動）に及ぼす影響を考慮してガスの混合を精度良く評価することができる。
【００２１】
具体的には、合流部２２におけるガスの混合は、次式によりデル化している。
λa (j+1) ＝α(j) ・λ(j) ＋｛１−α(j) ｝・λa (j)     ……（１）
ここで、λa は合流部２２における混合ガスの空燃比、λは合流部２２のガスに混合する気筒の排出ガスの空燃比、αは合流部２２のガスに混合する気筒の排出ガスの混合割合である。(j) は今回の演算タイミングにおける値を意味し、(j+1) は次回の演算タイミングにおける値を意味している。
【００２２】
また、混合ガスが合流部２２から空燃比センサ１６の位置に到達するまでの移動を評価する際に、合流部２２から空燃比センサ１６の位置までの距離もしくは容積と各燃焼気筒の排出ガス量とを考慮するものとする。このようにすれば、合流部２２から空燃比センサ１６の位置までガスが流動するのに要する時間を精度良く判定することができ、それによって、空燃比センサ１６の検出値がいずれの時期に排出されたガスの空燃比であるかを精度良く特定することができる。
【００２３】
具体的には、混合ガスが合流部２２から空燃比センサ１６の位置に到達するまでの移動は、次式によりデル化している。
λs (i) ＝λa （i −Ｖex／Ｖcy−β）      ……（２）
ここで、λs は空燃比センサ１６の検出値、βは燃焼間隔によるガス量の重なり（混合）を考慮に入れるパラメータ、Ｖexは合流部２２から空燃比センサ１６の位置までの容積、Ｖcyは各気筒の排出ガス量（シリンダ容積）、i は今回の演算タイミングを表している。Ｖex／Ｖcyは、合流部２２から空燃比センサ１６の位置までの容積と各燃焼気筒の排出ガス量とを考慮に入れるパラメータとなる。
【００２４】
λa （i −Ｖex／Ｖcy−β）は、現在(i) よりも（Ｖex／Ｖcy＋β）だけ過去に溯った時点のλa であることを意味している。燃焼間隔によるガス量の重なりを考慮に入れるパラメータβは、予め隣接する燃焼気筒のガス量の重なり度合いに応じて設定されている。本実施例では、図２に示すように、燃焼間隔によるガス量の重なりを３種類（重なり“大”、重なり“中”、重なり“無”）に分類し、重なり“大”の場合は、β＝−１、重なり“中”の場合は、β＝０、重なり“無”の場合は、β＝１に設定している。これは、隣接する燃焼気筒のガス量の重なりが大きくなるほど、合流部２２から集合排気管１４に流入するガス量が増加してガスの流速が速くなり、ガスが合流部２２から空燃比センサ１６の位置に到達するまでの時間が短くなるという事情を考慮するためである。
【００２５】
上記合流部２２におけるガスの混合と空燃比センサ１６の位置への移動のモデルを次式にまとめて気筒別空燃比推定モデルを構築し、この気筒別空燃比推定モデルを用いて気筒別空燃比を推定する。尚、オブザーバにはカルマンフィルタを用いる。
【００２６】
【数１】
2005207354.tif 000003
【００２７】
ここで、λn は＃ｎ気筒の気筒別空燃比、αn は合流部２２のガスに混合する＃ｎ気筒の排出ガスの混合割合、Ｃn は＃ｎ気筒の排気マニホールド１２の形状による混合ガスへの影響を考慮に入れるパラメータ、ｒn は＃ｎ気筒の排気マニホールド１２の長さによる混合ガスへの影響を考慮に入れるパラメータである。
【００２８】
上記（３）式を状態空間モデルに変換すると、次の（４ａ）、（４ｂ）式が導き出される。
Ｘ(i+1) ＝Ａ・Ｘ(i) ＋Ｂ・ｕ(i) ＋Ｗ(i)       ……（４ａ）
Ｙ(i) ＝Ｃ・Ｘ(i) ＋Ｄ・ｕ(i)                 ……（４ｂ）
ここで、Ａ，Ｂ，Ｃ，Ｄはモデルのパラメータ、Ｙは空燃比センサ１６の検出値、Ｘは状態変数としての気筒別空燃比の影響の総和、Ｗはノイズである。
【００２９】
更に、上記（４ａ）、（４ｂ）式によりカルマンフィルタを設計すると、次の（５）式が得られる。
Ｘ＾(k+1｜k)＝Ａ・Ｘ＾(k｜k-1)＋Ｋ｛Ｙ(k) −Ｃ・Ａ・Ｘ＾(k｜k-1)｝  ……（５）  ここで、Ｘ＾（エックスハット）は気筒別空燃比の影響の総和の推定値、Ｋはカルマンゲインである。Ｘ＾(k+1｜k)の意味は、時間(k) の推定値により時間(k+1) の推定値を求めることを表す。
【００３０】
以上のようにして、気筒別空燃比推定モデルをカルマンフィルタ型オブザーバにて構成することにより、燃焼サイクルの進行に伴い気筒別空燃比の影響の総和を順次推定することができる。気筒別空燃比は、（３）式の逆変換により推定できる。
【００３１】
ＥＣＵ２０は、図３乃至図５に示す気筒別空燃比制御用の各ルーチンを実行することで、上記気筒別空燃比推定モデルにより各バンクの空燃比センサ１６の検出値に基づいて各バンクの気筒別空燃比を推定し、各バンク毎に気筒間の空燃比ばらつきを少なくするように各気筒の燃料噴射量を補正する気筒別空燃比制御を実行する。以下、各ルーチンの処理内容を説明する。
【００３２】
［気筒別空燃比制御メインルーチン］
図３の気筒別空燃比制御メインルーチンは、クランク角センサ（図示せず）の出力パルスに同期して所定クランク角毎（例えば３０℃Ａ毎）に起動される。本ルーチンが起動されると、まずステップ１０１で、後述する図４の実行条件判定ルーチンを実行して、気筒別空燃比制御の実行条件が成立しているか否かを判定する。この後、ステップ１０２に進み、図４の実行条件判定ルーチンでセットされた実行フラグがＯＮであるか否かで、気筒別空燃比制御の実行条件が成立しているか否かを判定する。その結果、実行フラグがＯＦＦ（実行条件が不成立）と判定された場合は、以降の処理を行うことなく、本ルーチンを終了する。
【００３３】
一方、実行フラグがＯＮ（実行条件成立）と判定された場合は、ステップ１０３に進み、現在のクランク角が各気筒の空燃比検出タイミング（空燃比センサ１６の出力のサンプルタイミング）であるか否かを判定し、空燃比検出タイミングでなければ、以降の処理を行うことなく、本ルーチンを終了する。
【００３４】
これに対して、現在のクランク角が空燃比検出タイミングであれば、ステップ１０４に進み、後述する図５の気筒別空燃比制御実行ルーチンを起動して、気筒別空燃比制御を実行する。
【００３５】
［実行条件判定ルーチン］
図４の実行条件判定ルーチンは、図３の気筒別空燃比制御メインルーチンのステップ１０１で実行されるサブルーチンである。本ルーチンが起動されると、まずステップ２０１で、空燃比センサ１６が使用可能な状態であるか否かを判定する。ここで、使用可能な状態とは、例えば、空燃比センサ１６が活性状態で、且つ、故障していない状態であることである。空燃比センサ１６が使用可能な状態でなければ、以降の処理をを行うことなく、本ルーチンを終了する。
【００３６】
一方、空燃比センサ１６が使用可能な状態であれば、ステップ２０２に進み、冷却水温が所定温度以上（エンジン１１が暖機状態）であるか否かを判定し、所定温度未満であれば、以降の処理をを行うことなく、本ルーチンを終了する。冷却水温が所定温度以上であれば、ステップ２０３に進み、エンジン回転速度と負荷（例えば吸気管圧力）とをパラメータとする運転領域マップを参照して、現在のエンジン運転領域が気筒別空燃比制御の実行領域であるか否かを判定する。高回転域や低負荷域では、気筒別空燃比の推定精度が悪化されるため、気筒別空燃比制御が禁止される。
【００３７】
現在のエンジン運転領域が気筒別空燃比制御の実行領域であれば、ステップ２０４に進み、実行フラグをＯＮにセットし、気筒別空燃比制御の実行領域でなければ、ステップ２０５に進み、実行フラグをＯＦＦにセットする。
【００３８】
［気筒別空燃比制御実行ルーチン］
図５の気筒別空燃比制御実行ルーチンは、図３の気筒別空燃比制御メインルーチンのステップ１０４で実行されるサブルーチンである。本ルーチンが起動されると、まずステップ３０１で、空燃比センサ１６の出力（空燃比検出値）を読み込み、次のステップ３０２で、前記気筒別空燃比推定モデルを用いて今回の空燃比推定対象となる気筒の空燃比を空燃比センサ１６の検出値に基づいて推定する。このステップ３０２の処理が特許請求の範囲でいう気筒別空燃比推定手段としての役割を果たす。この後、ステップ３０３に進み、全気筒の推定空燃比の平均値を算出して、その平均値を基準空燃比（全気筒の目標空燃比）に設定する。
【００３９】
この後、ステップ３０４に進み、各気筒の推定空燃比と基準空燃比との偏差を算出して、その偏差が小さくなるように気筒別補正量を算出した後、ステップ３０５に進み、気筒別補正量に基づいて気筒別燃料噴射量を補正することで、各気筒に供給する混合気の空燃比を各気筒毎に補正して気筒間の空燃比ばらつきを少なくするように制御する。
【００４０】
以上説明した本実施例では、気筒別空燃比を推定する際に、隣接する燃焼気筒から排出されるガスの混合と、混合したガスが空燃比センサ１６の位置に到達するまでの移動とを考慮するようにしたので、気筒別空燃比の推定値に燃焼間隔による影響を精度良く反映させることができて、気筒別空燃比を精度良く推定することができる。
【００４１】
尚、本発明は、Ｖ型エンジンに限定されず、直列エンジン、水平対向エンジン等にも適用することができ、また、排気系が２系統のエンジンの他に、排気系が１系統のエンジンにも適用することができる。
【図面の簡単な説明】
【００４２】
【図１】本発明の一実施例におけるエンジン排気系の概略構成図である。
【図２】隣接する燃焼気筒の排出ガスの重なりを説明する図である。
【図３】気筒別空燃比制御メインルーチンの処理の流れを示すフローチャートである。
【図４】実行条件判定ルーチンの処理の流れを示すフローチャートである。
【図５】気筒別空燃比制御実行ルーチンの処理の流れを示すフローチャートである。
【符号の説明】
【００４３】
１１…エンジン（内燃機関）、１２…排気マニホールド、１４…集合排気管、１６…空燃比センサ、１８…触媒、２０…ＥＣＵ（気筒別空燃比推定手段，気筒別空燃比制御手段）、２２…合流部
【図１】
2005207354.tif 000004
【図２】
2005207354.tif 000005
【図３】
2005207354.tif 000006
【図４】
2005207354.tif 000007
【図５】
2005207354.tif 000008