(57)【特許請求の範囲】
【請求項１】
多気筒内燃機関の排気通路に設けられた空燃比センサと、
前記空燃比センサの出力を所定のサンプル周期毎に取得する取得手段と、
前記取得手段により取得された複数のセンサ出力から当該センサ出力の変動度合いに相関するパラメータを算出するパラメータ算出手段と、
前記パラメータ算出手段により算出されたパラメータを所定の判定値と比較して気筒間空燃比ばらつき異常の有無を判定する判定手段と、
筒内に流入している空気量である筒内空気量を前記サンプル周期毎に推定する推定手段と、
前記推定手段により推定された前記サンプル周期毎の筒内空気量の時系列データを記憶する記憶手段と、
前記取得手段が前記センサ出力を取得する毎に、この取得時点より輸送遅れ時間だけ前の時点での筒内空気量の値を前記時系列データに基づいて求め、この求めた筒内空気量の値に対応した対応判定値を算出する対応判定値算出手段と、
前記対応判定値算出手段によって算出された複数の対応判定値に基づき前記判定値を決定する判定値決定手段と、
を備え、
前記輸送遅れ時間が、前記サンプル周期の整数倍に等しい予め定められた時間であり、
前記サンプル周期が、クランク角の単位を有する
ことを特徴とする気筒間空燃比ばらつき異常検出装置。
【請求項２】
前記輸送遅れ時間が、前記サンプル周期の複数の整数倍に等しい予め定められた時間である
ことを特徴とする請求項１に記載の気筒間空燃比ばらつき異常検出装置。
【請求項３】
前記輸送遅れ時間が、前記筒内から前記空燃比センサに排気ガスが到達するまでの時間である
ことを特徴とする請求項１または２に記載の気筒間空燃比ばらつき異常検出装置。
【請求項４】
前記対応判定値算出手段は、前記筒内空気量の値に対応した対応判定値を所定のマップに従って算出する
ことを特徴とする請求項１〜３のいずれか一項に記載の気筒間空燃比ばらつき異常検出装置。
【請求項５】
前記判定値決定手段は、前記複数の対応判定値の平均値を前記判定値として決定する
ことを特徴とする請求項１〜４のいずれか一項に記載の気筒間空燃比ばらつき異常検出装置。
【請求項６】
前記内燃機関の吸気通路に設けられた吸入空気量検出手段をさらに備え、
前記推定手段は、前記吸入空気量検出手段により検出された吸入空気量に基づき前記筒内空気量を推定する
ことを特徴とする請求項１〜５のいずれか一項に記載の気筒間空燃比ばらつき異常検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、多気筒内燃機関の気筒間空燃比のばらつき異常を検出するための装置に係り、特に、多気筒内燃機関において気筒間の空燃比が比較的大きくばらついていることを検出する装置に関する。
【背景技術】
【０００２】
一般に、触媒を利用した排気浄化システムを備える内燃機関では、排気中有害成分の触媒による浄化を高効率で行うため、内燃機関で燃焼される混合気の空気と燃料との混合割合、すなわち空燃比のコントロールが欠かせない。こうした空燃比の制御を行うため、内燃機関の排気通路に空燃比センサを設け、これによって検出された空燃比を所定の目標空燃比に一致させるようフィードバック制御を実施している。
【０００３】
一方、多気筒内燃機関においては、通常全気筒に対し同一の制御量を用いて空燃比制御を行うため、空燃比制御を実行したとしても実際の空燃比が気筒間でばらつくことがある。このときばらつきの程度が小さければ、空燃比フィードバック制御で吸収可能であり、また触媒でも排気中有害成分を浄化処理可能なので、排気エミッションに影響を与えず、特に問題とならない。
【０００４】
しかし、例えば一部の気筒の燃料噴射系が故障するなどして、気筒間の空燃比が大きくばらつくと、排気エミッションを悪化させてしまい、問題となる。このような排気エミッションを悪化させる程の大きな空燃比ばらつきは異常として検出するのが望ましい。特に自動車用内燃機関の場合、排気エミッションが悪化した車両の走行を未然に防止するため、気筒間空燃比ばらつき異常を車載状態（オンボード）で検出することが要請されている。
【０００５】
例えば特許文献１には、空燃比センサ出力の単位時間当たりの変動量に基づいてばらつき異常を検出する装置が開示されている。この装置では、当該変動量が吸入空気量に応じて変化することに鑑み、吸入空気量に応じて変動量を補正している。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】国際公開第２０１１／０３０４５１号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、気筒間空燃比ばらつき異常を検出する場合、空燃比センサ出力の変動度合いに相関するパラメータを所定の判定値と比較してばらつき異常の有無を判定することが考えられる。
【０００８】
そして当該パラメータが吸入空気量に応じて変化することに鑑み、パラメータの比較対象である判定値を吸入空気量に応じて変更することが考えられる。
【０００９】
しかし、空燃比センサ出力の変動度合いに影響を及ぼすのは、吸気通路で検出された吸入空気量の大きさではなく、空燃比センサに実際に到達した排気ガスの流量である。また、空燃比センサは、これに実際に到達した排気ガスの空燃比を検出している。さらに、筒内で生成された排気ガスが空燃比センサに実際に到達するまでの間に、輸送遅れという時間差が存在する。
【００１０】
従来、これらの事項を十分考慮していなかったため、判定値を適切に設定することが困難であった。このため誤判定の虞があり、十分な検出精度を確保するのが困難であった。
【００１１】
そこで本発明は、上記事情に鑑みて創案され、その目的は、判定値を適切に設定して検出精度を向上し得る気筒間空燃比ばらつき異常検出装置を提供することにある。
【課題を解決するための手段】
【００１２】
本発明の一の態様によれば、
多気筒内燃機関の排気通路に設けられた空燃比センサと、
前記空燃比センサの出力を所定のサンプル周期毎に取得する取得手段と、
前記取得手段により取得された複数のセンサ出力から当該センサ出力の変動度合いに相関するパラメータを算出するパラメータ算出手段と、
前記パラメータ算出手段により算出されたパラメータを所定の判定値と比較して気筒間空燃比ばらつき異常の有無を判定する判定手段と、
筒内に流入している空気量である筒内空気量を前記サンプル周期毎に推定する推定手段と、
前記推定手段により推定された前記サンプル周期毎の筒内空気量の時系列データを記憶する記憶手段と、
前記取得手段が前記センサ出力を取得する毎に、この取得時点より輸送遅れ時間だけ前の時点での筒内空気量の値を前記時系列データに基づいて求め、この求めた筒内空気量の値に対応した対応判定値を算出する対応判定値算出手段と、
前記対応判定値算出手段によって算出された複数の対応判定値に基づき前記判定値を決定する判定値決定手段と、
を備えることを特徴とする気筒間空燃比ばらつき異常検出装置が提供される。
【００１３】
好ましくは、前記輸送遅れ時間が、前記サンプル周期の整数倍に等しい予め定められた時間である。
【００１４】
好ましくは、前記輸送遅れ時間が、前記サンプル周期の複数の整数倍に等しい予め定められた時間である。
【００１５】
好ましくは、前記サンプル周期が、クランク角の単位を有する。
【００１６】
好ましくは、前記輸送遅れ時間が、前記筒内から前記空燃比センサに排気ガスが到達するまでの時間である。
【００１７】
好ましくは、前記対応判定値算出手段は、前記筒内空気量の値に対応した対応判定値を所定のマップに従って算出する。
【００１８】
好ましくは、前記判定値決定手段は、前記複数の対応判定値の平均値を前記判定値として決定する。
【００１９】
好ましくは、前記気筒間空燃比ばらつき異常検出装置は、前記内燃機関の吸気通路に設けられた吸入空気量検出手段をさらに備え、前記推定手段は、前記吸入空気量検出手段により検出された吸入空気量に基づき前記筒内空気量を推定する。
【発明の効果】
【００２０】
本発明によれば、判定値を適切に設定して検出精度を向上することができるという、優れた効果が発揮される。
【図面の簡単な説明】
【００２１】
【図１】本発明の実施形態に係る内燃機関の概略図である。
【図２】触媒前センサおよび触媒後センサの出力特性を示すグラフである。
【図３】気筒間空燃比ばらつき度合いに応じた排気空燃比の変動を示すグラフである。
【図４】図３のＵ部に相当する拡大図である。
【図５】インバランス割合と出力変動パラメータの関係を示すグラフである。
【図６】出力変動パラメータと吸入空気量の関係を示すグラフである。
【図７】筒内空気量と、触媒前センサの位置における排気ガス流量との時間的推移を示すタイムチャートである。
【図８】メモリの記憶内容を示す図である。
【図９】ばらつき異常検出ルーチンのフローチャートである。
【発明を実施するための形態】
【００２２】
以下、本発明の実施形態を添付図面に基づき説明する。
【００２３】
図１は、本実施形態に係る内燃機関の概略図である。図示されるように、内燃機関（エンジン）１は、シリンダブロック２に形成された燃焼室３の内部で燃料および空気の混合気を燃焼させ、燃焼室３内でピストンを往復移動させることにより動力を発生する。本実施形態の内燃機関１は自動車に搭載された多気筒内燃機関であり、より具体的には直列４気筒火花点火式内燃機関である。内燃機関１は＃１〜＃４気筒を備える。但し気筒数、形式等は特に限定されない。
【００２４】
図示しないが、内燃機関１のシリンダヘッドには吸気ポートを開閉する吸気弁と、排気ポートを開閉する排気弁とが気筒ごとに配設されており、各吸気弁および各排気弁はカムシャフトによって開閉させられる。シリンダヘッドの頂部には、燃焼室３内の混合気に点火するための点火プラグ７が気筒ごとに取り付けられている。
【００２５】
各気筒の吸気ポートは気筒毎の枝管４を介して吸気集合室であるサージタンク８に接続されている。サージタンク８の上流側には吸気管１３が接続されており、吸気管１３の上流端にはエアクリーナ９が設けられている。そして吸気管１３には、上流側から順に、吸入空気量を検出するためのエアフローメータ５（吸入空気量検出手段）と、電子制御式のスロットルバルブ１０とが組み込まれている。吸気ポート、枝管４、サージタンク８及び吸気管１３により吸気通路が形成される。
【００２６】
吸気通路、特に吸気ポート内に燃料を噴射するインジェクタ（燃料噴射弁）１２が気筒ごとに配設されている。インジェクタ１２から噴射された燃料は吸入空気と混合されて混合気をなし、この混合気が吸気弁の開弁時に燃焼室３に吸入され、ピストンで圧縮され、点火プラグ７で点火燃焼させられる。なおインジェクタは燃焼室３内に燃料を直接噴射するものであってもよい。
【００２７】
一方、各気筒の排気ポートは排気マニフォールド１４に接続される。排気マニフォールド１４は、その上流部をなす気筒毎の枝管１４ａと、その下流部をなす排気集合部１４ｂとからなる。排気集合部１４ｂの下流側には排気管６が接続されている。排気ポート、排気マニフォールド１４及び排気管６により排気通路が形成される。
【００２８】
排気管６の上流側と下流側にはそれぞれ三元触媒からなる触媒、すなわち上流触媒１１と下流触媒１９が直列に取り付けられている。これら触媒１１，１９は酸素吸蔵能（Ｏ2ストレージ能）を有する。すなわち、触媒１１，１９は、排気ガスの空燃比がストイキ（理論空燃比、例えばＡ／Ｆ＝１４．６）より大きい（リーンな）ときに排気ガス中の過剰酸素を吸蔵し、ＮＯｘを還元する。また触媒１１，１９は、排気ガスの空燃比がストイキより小さい（リッチな）ときに吸蔵酸素を放出し、排気ガス中のＨＣ，ＣＯを酸化する。
【００２９】
上流触媒１１の上流側及び下流側にそれぞれ排気ガスの空燃比を検出するための第１及び第２の空燃比センサ、即ち触媒前センサ１７及び触媒後センサ１８が設置されている。これら触媒前センサ１７及び触媒後センサ１８は、上流触媒１１の直前及び直後の位置に設置され、排気中の酸素濃度に基づいて空燃比を検出する。このように上流触媒１１の上流側の排気合流部に単一の触媒前センサ１７が設置されている。本実施形態の場合、触媒前センサ１７が本発明にいう「空燃比センサ」に該当する。
【００３０】
上述の点火プラグ７、スロットルバルブ１０及びインジェクタ１２等は、制御手段としての電子制御ユニット（以下ＥＣＵと称す）２０に電気的に接続されている。ＥＣＵ２０は、何れも図示されないＣＰＵ、ＲＯＭ、ＲＡＭ、入出力ポート、および記憶装置等を含むものである。またＥＣＵ２０には、図示されるように、前述のエアフローメータ５、触媒前センサ１７、触媒後センサ１８のほか、内燃機関１のクランク角を検出するクランク角センサ１６、アクセル開度を検出するアクセル開度センサ１５、その他の各種センサが図示されないＡ／Ｄ変換器等を介して電気的に接続されている。ＥＣＵ２０は、各種センサの検出値等に基づいて、所望の出力が得られるように、点火プラグ７、スロットルバルブ１０、インジェクタ１２等を制御し、点火時期、燃料噴射量、燃料噴射時期、スロットル開度等を制御する。
【００３１】
スロットルバルブ１０にはスロットル開度センサ（図示せず）が設けられ、スロットル開度センサからの信号がＥＣＵ２０に送られる。ＥＣＵ２０は、通常、アクセル開度に応じて定まる目標スロットル開度に、スロットルバルブ１０の開度（スロットル開度）をフィードバック制御する。
【００３２】
ＥＣＵ２０は、エアフローメータ５からの信号に基づき、単位時間当たりの吸入空気の量である吸入空気量すなわち吸気流量を検出する。そしてＥＣＵ２０は、検出したアクセル開度、スロットル開度および吸入空気量の少なくとも一つに基づき、エンジン１の負荷を検出する。
【００３３】
ＥＣＵ２０は、クランク角センサ１６からのクランクパルス信号に基づき、クランク角自体を検出すると共にエンジン１の回転数を検出する。ここで「回転数」とは単位時間当たりの回転数のことをいい、回転速度と同義である。本実施形態では１分間当たりの回転数ｒｐｍのことをいう。
【００３４】
触媒前センサ１７は所謂広域空燃比センサからなり、比較的広範囲に亘る空燃比を連続的に検出可能である。図２に触媒前センサ１７の出力特性を示す。図示するように、触媒前センサ１７は、排気空燃比に比例した大きさの電圧信号Ｖｆを出力する。排気空燃比がストイキであるときの出力電圧はＶｒｅｆｆ（例えば約３．３Ｖ）である。
【００３５】
他方、触媒後センサ１８は所謂Ｏ2センサからなり、ストイキを境に出力値が急変する特性を持つ。図２に触媒後センサ１８の出力特性を示す。図示するように、排気空燃比がストイキであるときの出力電圧、すなわちストイキ相当値はＶｒｅｆｒ（例えば０．４５Ｖ）である。触媒後センサ１８の出力電圧は所定の範囲（例えば０〜１Ｖ）内で変化する。排気空燃比がストイキよりリーンのとき、触媒後センサの出力電圧はストイキ相当値Ｖｒｅｆｒより低くなり、排気空燃比がストイキよりリッチのとき、触媒後センサの出力電圧はストイキ相当値Ｖｒｅｆｒより高くなる。
【００３６】
上流触媒１１及び下流触媒１９は、それぞれに流入する排気ガスの空燃比Ａ／Ｆがストイキ近傍のときに排気中の有害成分であるＮＯｘ，ＨＣおよびＣＯを同時に浄化する。この三者を同時に高効率で浄化できる空燃比の幅（ウィンドウ）は比較的狭い。
【００３７】
そこで通常運転時、上流触媒１１に流入する排気ガスの空燃比がストイキ近傍に制御されるように、空燃比フィードバック制御がＥＣＵ２０により実行される。この空燃比フィードバック制御は、触媒前センサ１７によって検出された排気空燃比を所定の目標空燃比であるストイキに一致させるような主空燃比制御（主空燃比フィードバック制御）と、触媒後センサ１８によって検出された排気空燃比をストイキに一致させるような補助空燃比制御（補助空燃比フィードバック制御）とからなる。
【００３８】
なお、このように目標空燃比をストイキとする空燃比フィードバック制御をストイキ制御という。ストイキは基準空燃比をなす。
【００３９】
さて、例えば全気筒のうちの一部の気筒（特に１気筒）のインジェクタ１２が故障し、気筒間に空燃比のばらつき（インバランス：imbalance）が発生する場合がある。例えば＃１気筒のインジェクタ１２が故障し、＃１気筒が他の＃２、＃３及び＃４気筒よりも燃料噴射量が多くなり、その空燃比が大きくリッチ側にずれる場合等である。このときでも前述のストイキ制御により比較的大きな補正量を与えれば、触媒前センサ１７に供給されるトータルガスの空燃比をストイキに制御できる場合がある。しかし、気筒別に見ると、＃１気筒がストイキより大きくリッチ、＃２、＃３及び＃４気筒がストイキより若干リーンであり、全体のバランスとしてストイキとなっているに過ぎず、エミッション上好ましくないことは明らかである。そこで本実施形態では、かかる気筒間空燃比ばらつき異常を検出する装置が装備されている。
【００４０】
［ばらつき異常検出の態様］
以下、本実施形態におけるばらつき異常検出の態様を説明する。
【００４１】
図３に示すように、気筒間空燃比ばらつきが発生すると、１エンジンサイクル（＝７２０°ＣＡ）内での排気空燃比の変動が大きくなる。（Ｂ）の空燃比線図ａ，ｂ，ｃはそれぞればらつき無し、１気筒のみ２０％のインバランス割合でリッチずれ、及び１気筒のみ５０％のインバランス割合でリッチずれの場合の、触媒前センサ１７による検出空燃比Ａ／Ｆを示す。見られるように、ばらつき度合いが大きくなるほど空燃比変動の振幅が大きくなる。
【００４２】
ここでインバランス割合（％）とは、気筒間空燃比のばらつき度合いに相関する一つのパラメータである。即ち、インバランス割合とは、全気筒のうちある１気筒のみが燃料噴射量ズレを起こしている場合に、その燃料噴射量ズレを起こしている気筒（インバランス気筒）の燃料噴射量がどれくらいの割合で、燃料噴射量ズレを起こしていない気筒（バランス気筒）の燃料噴射量即ち基準噴射量からズレているかを示す値である。インバランス割合をＩＢ、インバランス気筒の燃料噴射量をＱｉｂ、バランス気筒の燃料噴射量即ち基準噴射量をＱｓとすると、ＩＢ＝（Ｑｉｂ−Ｑｓ）／Ｑｓで表される。インバランス割合ＩＢが大きいほど、インバランス気筒のバランス気筒に対する燃料噴射量ズレが大きく、空燃比ばらつき度合いは大きい。
【００４３】
図３から理解されるように、インバランス割合が大きいほど、すなわち気筒間空燃比のばらつき度合いが大きいほど、触媒前センサ１７の出力変動が大きくなる。
【００４４】
よってこの特性を利用し、本実施形態では、触媒前センサ１７の出力変動度合いに相関する出力変動パラメータＸを、気筒間空燃比ばらつき度合いに相関するパラメータとして用い、且つ出力変動パラメータＸを算出（あるいは検出）する。そしてこの算出された出力変動パラメータＸに基づき、ばらつき異常を検出する。なお、前述のインバランス割合は単に説明目的のためだけに用いる。
【００４５】
以下に出力変動パラメータＸの算出方法を説明する。図４は図３のＵ部に相当する拡大図であり、特に１エンジンサイクル内の触媒前センサ出力の変動を簡略的に示す。触媒前センサ出力としては、触媒前センサ１７の出力電圧Ｖｆを空燃比Ａ／Ｆに換算した値を用いる。但し触媒前センサ１７の出力電圧Ｖｆを直接用いることも可能である。
【００４６】
図４（Ｂ）に示すように、ＥＣＵ２０は、１エンジンサイクル内において、所定のサンプル周期τ毎に、触媒前センサ出力Ａ／Ｆの値を取得する。そして今回（ｎ）のタイミングで取得した値Ａ／Ｆnと、前回（ｎ−１）のタイミングで取得した値Ａ／Ｆn-1との差の絶対値（出力差という）ΔＡ／Ｆnを次式（１）により求める。この出力差ΔＡ／Ｆnは今回のタイミングにおける微分値あるいは傾きと言い換えることができる。
【００４７】
【数１】
0005397454.tif 000002
【００４８】
最も単純には、この出力差ΔＡ／Ｆnが触媒前センサ出力の変動を表す。変動度合いが大きくなるほど空燃比線図の傾きが大きくなり、出力差ΔＡ／Ｆnが大きくなるからである。そこで所定の１タイミングにおける出力差ΔＡ／Ｆnの値を出力変動パラメータとすることができる。
【００４９】
但し、本実施形態では精度向上のため、複数の出力差ΔＡ／Ｆnの平均値を出力変動パラメータとする。本実施形態では、１エンジンサイクルの間、各タイミング毎に出力差ΔＡ／Ｆnを積算し、最終積算値をサンプル数Ｎで除し、１エンジンサイクル内の差ΔＡ／Ｆnの平均値を求める。そしてさらに、Ｍエンジンサイクル分（例えばＭ＝１００）だけ出力差ΔＡ／Ｆnの平均値を積算し、最終積算値をサイクル数Ｍで除し、Ｍエンジンサイクル内の出力差ΔＡ／Ｆnの平均値を求める。こうして求められた最終的な平均値を出力変動パラメータＸとする。触媒前センサ出力の変動度合いが大きくなるほど出力変動パラメータＸは大きくなる。
【００５０】
ここで、サンプル周期τはクランク角の単位°ＣＡを有する。サンプル周期τは例えば１０°ＣＡとすることができ、この場合、１エンジンサイクル内のサンプル数Ｎは７２０／１０＝７２個となる。なおサンプル周期τは他のクランク角とすることもできる。
【００５１】
触媒前センサ出力Ａ／Ｆは増加する場合と減少する場合とがあるので、これら各場合の一方についてだけ上記出力差ΔＡ／Ｆnあるいはその平均値を求め、これを出力変動パラメータとしても良い。特に１気筒のみリッチずれの場合、当該１気筒に対応した排気ガスを触媒前センサが受けた時にその出力が急速に減少（リッチ側に変化）する傾向があるので、減少側のみの値をリッチずれ検出のために用いることも可能である。もっとも、これに限定されず、増加側の値のみを用いることも可能である。
【００５２】
また、触媒前センサ出力の変動度合いに相関する如何なる値をも出力変動パラメータとすることができる。例えば、１エンジンサイクル内における触媒前センサ出力の最大ピークと最小ピークの差（所謂ピークトゥピーク; peak to peak）、または２階微分値の最大ピークまたは最小ピークの絶対値に基づいて、出力変動パラメータを算出することもできる。触媒前センサ出力の変動度合いが大きいほど、触媒前センサ出力の最大ピークと最小ピークの差は大きくなり、また２階微分値の最大ピークまたは最小ピークの絶対値も大きくなるからである。
【００５３】
図５には、インバランス割合ＩＢ（％）と出力変動パラメータＸの関係を示す。図示されるように、インバランス割合ＩＢと出力変動パラメータＸの間には強い相関関係があり、インバランス割合ＩＢの絶対値が増加するほど空燃比変動パラメータＸも増加する。
【００５４】
算出された出力変動パラメータＸを、所定の判定値と比較して、ばらつき異常の有無が判定される。例えば、算出された出力変動パラメータＸが判定値以上であればばらつき異常あり（異常）、算出された出力変動パラメータＸが判定値より小さければばらつき異常なし（正常）と判定される。
【００５５】
ところで、出力変動パラメータＸの値は、気筒間空燃比のばらつき度合いのみならず、触媒前センサ１７に到達あるいは供給される排気ガス流量に応じても変化する。
【００５６】
図６には、出力変動パラメータＸと、エアフローメータ５により検出された吸入空気量Ｇａの関係を示す。図中、白三角および線ａはばらつき異常の無い正常な場合を示し、黒丸および線ｂはばらつき異常の有る異常な場合を示す。図示する関係は、吸入空気量一定という条件下で得られた関係であり、この場合、吸入空気量Ｇａは、触媒前センサ１７に実際に到達した排気ガス流量に等しいとみなせる。
【００５７】
図示されるように、正常な場合にも異常な場合にも、出力変動パラメータＸは吸入空気量Ｇａの増加につれ増加する傾向にある。但し、同一吸入空気量の場合、当然ながら、異常時の出力変動パラメータＸは正常時の出力変動パラメータＸより大きい。
【００５８】
そこで、正常時のパラメータ曲線ａと異常時のパラメータ曲線ｂとの間に、正常・異常の境界を定める判定値αを規定する判定値曲線ｃを設定する。この判定値曲線ｃも、出力変動パラメータＸと同様、吸入空気量Ｇａの増加につれ増加する傾向にある。
【００５９】
これにより、ある吸入空気量の下で算出された出力変動パラメータＸに対しては、同一の吸入空気量に対応する判定値曲線ｃ上の判定値αを用い、これらを比較することにより、吸入空気量の影響を加味した正確な判定を実行することができる。
【００６０】
本実施形態において、図６に示すような判定値曲線ｃは、予め実験的に求められ、マップ（または関数、以下同様）の形でＥＣＵ２０に記憶される。但し、詳しくは後述するが、本実施形態においては吸入空気量Ｇａの代わりに、筒内に流入している空気量である筒内空気量Ｇｃを用いるため、筒内空気量Ｇｃと判定値αの関係を規定する判定値曲線ｃがＥＣＵ２０に記憶されている。
【００６１】
ところで、出力変動パラメータＸに実質的に影響を及ぼすのは、吸気通路に設置されたエアフローメータ５により検出された吸入空気量Ｇａの大きさではなく、触媒前センサ１７に実際に到達した排気ガスの流量である。また触媒前センサ１７は、これに実際に到達した排気ガスの空燃比を検出している。さらに、筒内の燃焼室３で生成された排気ガスが触媒前センサ１７に実際に到達するまでの間に、輸送遅れという時間差が存在する。
【００６２】
従来、これらの事項を十分考慮していなかったため、判定値αを適切に設定することが困難であった。このため誤判定の虞があり、十分な検出精度を確保するのが困難であった。
【００６３】
以下、この点について図７を用いて説明する。図中（Ａ）は、筒内に流入している（あるいは存在する）空気量である筒内空気量Ｇｃを示す。また（Ｂ）は、触媒前センサ１７の位置に到達した排気ガスの流量（センサ位置ガス流量という）Ｇｅを示す。筒内に流入した空気が混合気となって燃焼し、排気ガスを生成し、この排気ガスが筒内から輸送遅れ時間Ｌだけ遅れて触媒前センサ１７に到達する。このため、筒内空気量Ｇｃとセンサ位置ガス流量Ｇｅとの間には、輸送遅れによる時間差があることを除いて、相関関係がある。
【００６４】
図示例において、あるサンプル時期ｎにおける筒内空気量はＧｃnであり、これより１サンプル周期（τ）前の筒内空気量はＧｃn-1であり、さらにこれより１サンプル周期前の筒内空気量はＧｃn-2である。
【００６５】
ここで本実施形態を説明する前に、便宜上比較例を説明する。サンプル時期ｎにおいて、触媒前センサ出力Ａ／Ｆnの値が取得されると共に、前記式（１）により出力差ΔＡ／Ｆnが算出される。
【００６６】
一方、サンプル時期ｎにおいて、触媒前センサ出力Ａ／Ｆnに対応する対応判定値αnも算出される。対応判定値αnは、詳しくは後述するが、判定値αを決定する際の基礎となる値である。具体的には、複数の対応判定値αnの平均値が判定値αをなす。
【００６７】
この比較例において、対応判定値αnは、サンプル時期ｎにおける筒内空気量Ｇｃnに基づき、図６に示したようなマップから算出される。すなわち、当該マップから筒内空気量Ｇｃnに対応した対応判定値αnが算出される。
【００６８】
しかし、サンプル時期ｎにおける筒内空気量Ｇｃnおよび対応判定値αnは、当該サンプル時期ｎにおけるセンサ位置ガス流量Ｇｅnに対応していない。当該センサ位置ガス流量Ｇｅnに対応するのは、サンプル時期ｎより輸送遅れ時間Ｌだけ前の筒内空気量であり、図示例の場合、サンプル時期ｎより２サンプル周期（２τ）前の筒内空気量Ｇｃn-2である。
【００６９】
従って結果的に、サンプル時期ｎにおいて算出された出力差ΔＡ／Ｆnと対応判定値αnは互いに対応しておらず、これが、最終的に算出される判定値αを不適切なものとする原因となる。
【００７０】
そこで本実施形態では、サンプル時期ｎより輸送遅れ時間Ｌだけ前の時点での筒内空気量、図示例ではサンプル時期ｎより２サンプル周期前の時点での筒内空気量Ｇｃn-2に基づき、図６に示したようなマップから、サンプル時期ｎにおける対応判定値αnを算出する。
【００７１】
これにより、サンプル時期ｎにおける出力差ΔＡ／Ｆnと対応判定値αnとを互いに対応したものとすることができる。そして最終的に算出される判定値αを適切なものとし、検出精度を向上することができる。
【００７２】
また本実施形態によれば、筒内空気量Ｇｃが変化し、これに起因して同一サンプル時期における筒内空気量Ｇｃとセンサ位置ガス流量Ｇｅとが対応しなくなる期間ＩＩにおいても、出力差ΔＡ／Ｆnに対応した適切な対応判定値αnを算出することができる。よって検出精度向上のみならず、ロバスト性向上、耐誤検出余裕確保、検出機会あるいは検出頻度確保に大変有利である。
【００７３】
すなわち、両者が対応する期間Ｉでは、比較例のように、サンプル時期ｎにおける筒内空気量Ｇｃnに基づき対応判定値αnを算出しても問題はない。しかし、両者が対応しなくなる期間ＩＩでは、比較例だと前述したようにサンプル時期ｎの出力差ΔＡ／Ｆnに対応した適切な対応判定値αnを得られない。
【００７４】
その結果、例えば筒内空気量Ｇｃが増加する加速時には、あるサンプル時期ｎにおけるセンサ位置ガス流量Ｇｅnに対して筒内空気量Ｇｃnが相対的に多くなり、これに伴って対応判定値αnも大きくなり、本来異常なのに誤って正常と判定してしまう虞がある。
【００７５】
逆に、筒内空気量Ｇｃが減少する減速時には、サンプル時期ｎにおけるセンサ位置ガス流量Ｇｅnに対して筒内空気量Ｇｃnが相対的に少なくなり、これに伴って対応判定値αnも小さくなり、本来正常なのに誤って異常と判定してしまう虞がある。
【００７６】
しかし、本実施形態だと、両者が対応しなくなる期間ＩＩにおいても両者の対応関係を保つことができ、適切な対応判定値αnを算出することが可能である。その結果、加減速に伴って筒内空気量Ｇｃが変動しても、適切な対応判定値αnを算出することができる。よってロバスト性を向上できると共に、耐誤検出余裕を確保できる。また、ばらつき異常検出は、所定期間内における吸入空気量Ｇａ（これは筒内空気量Ｇｃに相関する）の変動幅が所定値以内という定常運転条件下に限って行われる。本実施形態の場合だと、その吸入空気量Ｇａの変動幅を比較例よりも大きくすることができ、すなわち筒内空気量Ｇｃの変動を比較例よりも許容することができ、結果的に多くの検出機会あるいは検出頻度を確保できるようになる。
【００７７】
ところで、ＥＣＵ２０は、筒内空気量Ｇｃnをサンプル周期τ毎に推定し、この推定したサンプル周期τ毎の筒内空気量Ｇｃnの時系列データを自身のメモリに記憶するようになっている。
【００７８】
筒内空気量Ｇｃの推定に関して、ＥＣＵ２０は、エアフローメータ５により検出された吸入空気量Ｇａに基づいて筒内空気量Ｇｃを推定する。例えばＥＣＵ２０は、あるサンプル時期ｎにおいて、エアフローメータ５により検出された吸入空気量Ｇａの値を取得すると共に、この取得した吸入空気量Ｇａの値に基づいて、当該サンプル時期ｎまたはその後の所定サンプル時期における筒内空気量Ｇｃを推定または予測する。この推定方法については、マップを使用した方法、モデルを使用して計算する方法など、様々な方法が利用可能であり、当該方法は既に公知あるいは周知であるので詳細な説明を省略する。なお推定に際しては吸入空気量Ｇａ以外のパラメータ（例えばエンジン回転数等）を併せて使用してもよい。
【００７９】
図８は、メモリの記憶内容を示す。図示するように、メモリのアドレスＡＤ０には、今回のサンプル時期ｎに推定された筒内空気量Ｇｃnの値が格納される。メモリのアドレスＡＤ１には、１サンプル周期前の前回のサンプル時期ｎ−１に推定された筒内空気量Ｇｃn-1の値が格納される。同様に、メモリのアドレスＡＤ２には２回前のサンプル時期ｎ−２に推定された筒内空気量Ｇｃn-2の値が格納され、メモリのアドレスＡＤｉにはｉ回前のサンプル時期ｎ−ｉに推定された筒内空気量Ｇｃn-iの値が格納される。こうして格納された筒内空気量Ｇｃn〜Ｇｃn-iの各値が時系列データをなす。
【００８０】
このメモリの記憶内容はサンプル周期τ毎に更新され、常に最新の筒内空気量データがアドレスＡＤ０に格納される。そして最も古いアドレスＡＤｉの筒内空気量データはサンプル周期τ毎に消去される。
【００８１】
本実施形態では、今回のサンプル時期ｎにおいて、触媒前センサ出力Ａ／Ｆnが取得され且つ出力差ΔＡ／Ｆnが算出されると共に、Ｋ回前のサンプル時期ｎ−Ｋにおける筒内空気量Ｇｃn-Kの値が時系列データから読み出される。この読み出された筒内空気量Ｇｃn-Kが、今回のサンプル時期ｎにおける対応空気量ＧｃＬnとされる。そしてこの対応空気量ＧｃＬnの値に基づいて、今回のサンプル時期ｎにおける対応判定値αnが算出される。すなわち図６に示したマップから、対応空気量ＧｃＬnに対応した対応判定値αnが算出される。
【００８２】
Ｋは例えば２であり、この場合、図８に示すように、２回前のサンプル時期ｎ−２における筒内空気量Ｇｃn-2が対応空気量ＧｃＬnとされ、この対応空気量ＧｃＬnに基づいて今回のサンプル時期ｎにおける対応判定値αnが算出される。
【００８３】
Ｋは前記輸送遅れ時間Ｌを規定する定数である。すなわち、本実施形態では、サンプル周期τと定数Ｋの積が輸送遅れ時間Ｌに等しいとみなされる。定数Ｋは１以上の整数である。こうして輸送遅れ時間Ｌは、サンプル周期τの整数倍に等しい予め定められた時間とされる。また定数Ｋを２以上とした場合、輸送遅れ時間Ｌは、サンプル周期τの複数の整数倍に等しい予め定められた時間とされる。
【００８４】
定数Ｋは、排気ポート入口から触媒前センサ１７までの排気通路の長さ、容積、形状等や、排気量等のハード因子に起因してエンジン毎に定まる固有値ないし適合値であり、予め実験的に求められＥＣＵ２０に記憶されている。
【００８５】
適合値である定数Ｋを用い、サンプル周期τのＫ倍として輸送遅れ時間Ｌを定めることにより、非常にシンプルな方法で、輸送遅れ時間Ｌを規定し、輸送遅れ時間Ｌだけ前の対応空気量ＧｃＬnを求め、且つこれを用いて今回サンプル時期ｎの対応判定値αnを算出することができる。
【００８６】
なお、比較例においては、今回のサンプル時期ｎに当該時期の筒内空気量Ｇｃnに基づき対応判定値αnが算出される。これによって、今回のサンプル時期ｎにおける触媒前センサ出力Ａ／Ｆnに対して筒内空気量Ｇｃが対応しなくなり、不適切な対応判定値αnが算出されてしまう。本実施形態においては、今回のサンプル時期ｎより輸送遅れ時間Ｌ＝Ｋτだけ前の時点の筒内空気量Ｇｃn-Kに基づき対応判定値αnが算出されるので、適切な対応判定値αnが算出可能である。
【００８７】
ここで、エンジン回転数が高くなるほど、ガス流速が増大し、時間単位での輸送遅れ時間Ｌが短くなる。しかし、本実施形態ではサンプル周期τをクランク角単位で設定し、輸送遅れ時間Ｌをサンプル周期τのＫ倍として規定している。よってこのような回転上昇によるガス流速増大が生じても、これに対応して時間単位での輸送遅れ時間Ｌを短くすることができる。結局、エンジン回転数やガス流速の変化があっても常に正確な対応空気量ＧｃＬnを得、正確な対応判定値αnを算出することができる。
【００８８】
なお、定数Ｋは２に限らず、試験結果に応じて任意な値に設定できる。メモリ中には少なくともＫサンプル周期前までのデータを記憶しておけばよく、必ずしも（Ｋ＋１）サンプル周期以上前のデータを記憶する必要はないが、本実施形態では便宜上（Ｋ＋１）サンプル周期以上前のデータも記憶している。
【００８９】
メモリ中の時系列データに基づき、補完計算により、例えば２．５サンプル周期前の筒内空気量を対応空気量ＧｃＬnとして求めることも可能である。この場合、２サンプル周期前の筒内空気量Ｇｃn-2と３サンプル周期前の筒内空気量Ｇｃn-3との平均値を対応空気量ＧｃＬnとして求め、これを今回サンプル時期ｎの対応判定値αnの算出の基礎とする。結局、輸送遅れ時間Ｌは、必ずしもサンプル周期τの整数倍に等しくなくてもよい。
【００９０】
以下、図９を参照して、本実施形態に係るばらつき異常検出ルーチンを説明する。このルーチンはＥＣＵ２０によりサンプル周期τ（°ＣＡ）毎に繰り返し実行される。
【００９１】
まずステップＳ１０１では、異常検出を行うのに適した所定の前提条件が成立しているか否かが判断される。この前提条件は、例えば次の各条件が成立したときに成立する。
（１）エンジンの暖機が終了している。ＥＣＵ２０は、図示しない水温センサで検出された水温が所定値（例えば７５℃）以上であるとき暖機終了と判断する。
（２）触媒前センサ１７および触媒後センサ１８が活性化している。ＥＣＵ２０は、両センサのインピーダンスがそれぞれ所定の活性温度相当の値になっているとき、両センサが活性化していると判断する。
（３）上流触媒１１および下流触媒１９が活性化している。ＥＣＵ２０は、別途推定した上流触媒１１および下流触媒１９の温度がそれぞれ所定の活性温度相当の値になっているとき、両触媒が活性化していると判断する。
（４）エンジンが定常運転中である。
（５）エアフローメータ５により検出された吸入空気量Ｇａが所定値（例えば２０（ｇ／ｓ））以上である。
（６）ストイキ制御中である。
【００９２】
前提条件が成立していない場合にはルーチンが終了される。他方、前提条件が成立している場合にはステップＳ１０２に進んで異常検出が実質的に開始または実行される。
【００９３】
条件（４）に関し、ＥＣＵ２０は、前述したように、所定期間内における吸入空気量Ｇａの変動幅が所定値以内であるときエンジンが定常運転中であると判断する。この条件（４）が非成立だと異常検出が開始されず、また異常検出の途中で条件（４）が非成立になると異常検出が途中で中止されてしまう。
【００９４】
本実施形態によれば、前述したように、吸入空気量Ｇａの変動に伴って筒内空気量Ｇｃが変動したとしても、その変動中、対応判定値αnを精度良く算出できる。よって吸入空気量Ｇａの変動幅を比較的大きく設定することができ、異常検出の開始条件および中止条件を緩和し、検出機会あるいは検出頻度をより多く確保することができる。
【００９５】
ステップＳ１０２において、今回のサンプル時期ｎないしタイミングにおける触媒前センサ出力Ａ／Ｆnが取得される。なお触媒前センサ出力Ａ／Ｆnは触媒前センサ１７の出力電圧Ｖｆを空燃比に換算した値である。
【００９６】
ステップＳ１０３において、今回のタイミングにおけるセンサ出力差ΔＡ／Ｆnが前式（１）より算出される。
【００９７】
ステップＳ１０４において、センサ出力差ΔＡ／Ｆnが積算され、すなわち今回のタイミングにおける積算センサ出力差ΣΔＡ／Ｆnが次式（２）より算出される。
【００９８】
【数２】
0005397454.tif 000003
【００９９】
ステップＳ１０５において、今回のタイミングにおける対応空気量ＧｃＬnが、図８を参照して説明したように、メモリから読み出される。本実施形態の場合、メモリ中のアドレスＡＤ２に格納された２サンプル周期前の筒内空気量Ｇｃn-2が、今回タイミングの対応空気量ＧｃＬnとして読み出される。
【０１００】
ステップＳ１０６において、今回タイミングの対応空気量ＧｃＬnに対応した対応判定値αnが図６に示したマップから算出される。
【０１０１】
ステップＳ１０７において、対応判定値αnが積算され、すなわち今回のタイミングにおける積算対応判定値Σαnが次式（３）より算出される。
【０１０２】
【数３】
0005397454.tif 000004
【０１０３】
ステップＳ１０８において、１エンジンサイクルが終了したか否かが判断される。終了してなければルーチンが終了され、終了した場合にはステップＳ１０９に進む。
【０１０４】
ステップＳ１０９では、今回の１エンジンサイクル終了時点における最終的な積算センサ出力差ΣΔＡ／Ｆnがサンプル数Ｎで除して平均化され、平均センサ出力差Ｒmが算出される。
【０１０５】
そしてステップＳ１１０において、平均センサ出力差Ｒmが積算され、すなわち今回のエンジンサイクル終了時点における積算平均センサ出力差ΣＲmが次式（４）より算出される。
【０１０６】
【数４】
0005397454.tif 000005
【０１０７】
同様に、ステップＳ１１１において、今回の１エンジンサイクル終了時点における最終的な積算対応判定値Σαnがサンプル数Ｎで除して平均化され、平均対応判定値Ｓmが算出される。
【０１０８】
そしてステップＳ１１２において、平均対応判定値Ｓmが積算され、すなわち今回のエンジンサイクル終了時点における積算平均対応判定値ΣＳmが次式（５）より算出される。
【０１０９】
【数５】
0005397454.tif 000006
【０１１０】
次に、ステップＳ１１３において、Ｍエンジンサイクル（但しＭは２以上の整数で、例えば１００）が終了したか否かが判断される。終了してなければルーチンが終了され、終了した場合にはステップＳ１１４に進む。
【０１１１】
ステップＳ１１４では、Ｍエンジンサイクル終了時点における最終的な積算平均センサ出力差ΣＲmがサイクル数Ｍで除して平均化され、出力変動パラメータＸが算出される。すなわち、複数（Ｎ×Ｍ個）のセンサ出力差ΔＡ／Ｆnの平均値が最終的な出力変動パラメータＸとして算出される。
【０１１２】
同様に、ステップＳ１１５では、Ｍエンジンサイクル終了時点における最終的な積算平均対応判定値ΣＳmがサイクル数Ｍで除して平均化され、判定値αが算出される。すなわち、複数（Ｎ×Ｍ個）の対応判定値αnの平均値が最終的な判定値αとして決定される。
【０１１３】
このように、各サンプル時期のセンサ出力Ａ／Ｆnおよびセンサ出力差ΔＡ／Ｆnに対応した対応判定値αnを平均化して判定値αを決定するので、適切な判定値αを決定することができる。
【０１１４】
次いで、ステップＳ１１６において、出力変動パラメータＸが判定値αと比較される。
【０１１５】
出力変動パラメータＸが判定値α以上である場合、ステップＳ１１７に進んでばらつき異常有り、すなわち異常と判定され、ルーチンが終了される。なお異常判定と同時に、異常の事実をユーザに知らせるべくチェックランプ等の警告装置を起動するのが好ましい。
【０１１６】
他方、出力変動パラメータＸが判定値α未満である場合、ステップＳ１１８に進んでばらつき異常無し、すなわち正常と判定され、ルーチンが終了される。
【０１１７】
以上、本発明の好適な実施形態を詳細に述べたが、本発明の実施形態は他にも様々なものが考えられる。上記の数値はあくまで例示であり、適宜変更が可能である。
【０１１８】
本発明の実施形態は前述の実施形態のみに限らず、特許請求の範囲によって規定される本発明の思想に包含されるあらゆる変形例や応用例、均等物が本発明に含まれる。従って本発明は、限定的に解釈されるべきではなく、本発明の思想の範囲内に帰属する他の任意の技術にも適用することが可能である。
【符号の説明】
【０１１９】
１  内燃機関（エンジン）
２  インジェクタ
３  燃焼室
５  エアフローメータ
６  排気管
１１  上流触媒
１４  排気マニフォールド
１７  触媒前センサ
１８  触媒後センサ
２０  電子制御ユニット（ＥＣＵ）
τ  サンプル周期
Ｘ  出力変動パラメータ
α  判定値
Ｇｃ  筒内空気量
Ｌ  輸送遅れ時間
αn  対応判定値
【図１】
0005397454.tif 000007
【図２】
0005397454.tif 000008
【図３】
0005397454.tif 000009
【図４】
0005397454.tif 000010
【図５】
0005397454.tif 000011
【図６】
0005397454.tif 000012
【図７】
0005397454.tif 000013
【図８】
0005397454.tif 000014
【図９】
0005397454.tif 000015