(57)【特許請求の範囲】
【請求項１】
多気筒内燃機関の排気通路に設けられた触媒と、
前記触媒の上流側に設けられた空燃比センサとしての触媒前センサと、
前記触媒の下流側に設けられた空燃比センサとしての触媒後センサと、
前記触媒前センサの出力に基づく主空燃比制御と、前記触媒後センサの出力に基づく補助空燃比制御とを実行する空燃比制御手段と、
前記触媒後センサの出力に基づき前記補助空燃比制御のための制御量を算出する制御量算出手段と、
前記内燃機関の回転変動を検出する回転変動検出手段と、
前記回転変動の検出値に基づき、気筒間空燃比ばらつき異常を検出する異常検出手段と、
前記ばらつき異常の検出中に前記制御量のガード範囲を縮小するガード範囲縮小手段と、
を備えることを特徴とする気筒間空燃比ばらつき異常検出装置。
【請求項２】
前記異常検出手段は、前記回転変動の検出値に基づき、一部気筒の空燃比が他の気筒の空燃比に対しリーン側にずれているリーンずれ異常を検出し、
前記ガード範囲縮小手段は、前記リーンずれ異常の検出中に前記制御量のガード範囲を縮小する
ことを特徴とする請求項１に記載の気筒間空燃比ばらつき異常検出装置。
【請求項３】
前記ガード範囲縮小手段は、ばらつき異常検出開始と同時に前記制御量のガード範囲を縮小し、
前記空燃比制御手段は、ばらつき異常検出開始直前の制御量が縮小されたガード範囲から外れていたとき、その検出開始直前の制御量に最も近く且つ縮小されたガード範囲を定めるガード値を前記制御量として、前記ばらつき異常検出中の前記補助空燃比制御を実行する
ことを特徴とする請求項１または２に記載の気筒間空燃比ばらつき異常検出装置。
【請求項４】
前記ガード範囲縮小手段は、ばらつき異常検出終了と同時に前記制御量のガード範囲を元の規定のガード範囲に復帰し、
前記空燃比制御手段は、ばらつき異常検出開始直前の制御量を前記制御量として、ばらつき異常検出終了後の前記補助空燃比制御を開始する
ことを特徴とする請求項３に記載の気筒間空燃比ばらつき異常検出装置。
【請求項５】
前記異常検出手段は、前記回転変動の検出値が所定の異常判定値以上であるとき、ばらつき異常ありと判定する
ことを特徴とする請求項１〜４のいずれか一項に記載の気筒間空燃比ばらつき異常検出装置。
【請求項６】
前記回転変動検出手段は、前記内燃機関の気筒毎の回転変動を検出する
ことを特徴とする請求項１〜５のいずれか一項に記載の気筒間空燃比ばらつき異常検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、気筒間空燃比のばらつき異常を検出するための装置に係り、特に、多気筒内燃機関において気筒間の空燃比が比較的大きくばらついていることを検出する装置に関する。
【背景技術】
【０００２】
一般に、触媒を利用した排気浄化システムを備える内燃機関では、排気中有害成分の触媒による浄化を高効率で行うため、内燃機関で燃焼される混合気の空気と燃料との混合割合、すなわち空燃比のコントロールが欠かせない。こうした空燃比の制御を行うため、内燃機関の排気通路に空燃比センサを設け、これによって検出された空燃比を所定の目標空燃比に一致させるようフィードバック制御を実施している。
【０００３】
一方、多気筒内燃機関においては、通常全気筒に対し同一の制御量を用いて空燃比制御を行うため、空燃比制御を実行したとしても実際の空燃比が気筒間でばらつくことがある。このときばらつきの程度が小さければ、空燃比フィードバック制御で吸収可能であり、また触媒でも排気中有害成分を浄化処理可能なので、排気エミッションに影響を与えず、特に問題とならない。
【０００４】
しかし、例えば一部の気筒の燃料噴射系が故障するなどして、気筒間の空燃比が大きくばらつくと、排気エミッションを悪化させてしまい、問題となる。このような排気エミッションを悪化させる程の大きな空燃比ばらつきは異常として検出するのが望ましい。特に自動車用内燃機関の場合、排気エミッションの悪化した車両の走行を未然に防止するため、気筒間空燃比ばらつき異常を車載状態（オンボード）で検出することが要請されており、最近ではこれを法規制化する動きもある。
【０００５】
例えば特許文献１に記載の装置では、内燃機関の回転変動に基づいて各気筒の気筒間ばらつき値を算出するようにしている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００４−５２６２０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、空燃比フィードバック制御に関し、触媒の上流側と下流側とにそれぞれ空燃比センサとしての触媒前センサおよび触媒後センサを設け、触媒前センサの出力に基づく主空燃比制御と、触媒後センサの出力に基づく補助空燃比制御とを実行する場合がある。
【０００８】
そして補助空燃比制御のための制御量については、その値が過剰になるのを防止するため、規定のガード範囲が設けられているのが一般的である。
【０００９】
他方、内燃機関の回転変動を検出し、その検出値に基づいて気筒間空燃比ばらつき異常を検出することが考えられる。
【００１０】
こうした場合、本発明者の鋭意研究の結果によれば、補助空燃比制御量と回転変動との間に相関性があり、補助空燃比制御量の絶対値が、規定のガード範囲のガード値付近まで大きくなっていると、ばらつき異常検出時に誤検出が発生する虞があることが判明した。
【００１１】
そこで本発明は、上記の事情に鑑みて創案され、その目的は、ばらつき異常検出時における誤検出を防止できる気筒間空燃比ばらつき異常検出装置を提供することにある。
【課題を解決するための手段】
【００１２】
本発明の一の態様によれば、
多気筒内燃機関の排気通路に設けられた触媒と、
前記触媒の上流側に設けられた空燃比センサとしての触媒前センサと、
前記触媒の下流側に設けられた空燃比センサとしての触媒後センサと、
前記触媒前センサの出力に基づく主空燃比制御と、前記触媒後センサの出力に基づく補助空燃比制御とを実行する空燃比制御手段と、
前記触媒後センサの出力に基づき前記補助空燃比制御のための制御量を算出する制御量算出手段と、
前記内燃機関の回転変動を検出する回転変動検出手段と、
前記回転変動の検出値に基づき、気筒間空燃比ばらつき異常を検出する異常検出手段と、
前記ばらつき異常の検出中に前記制御量のガード範囲を縮小するガード範囲縮小手段と、
を備えることを特徴とする気筒間空燃比ばらつき異常検出装置が提供される。
【００１３】
好ましくは、前記異常検出手段は、前記回転変動の検出値に基づき、一部気筒の空燃比が他の気筒の空燃比に対しリーン側にずれているリーンずれ異常を検出し、
前記ガード範囲縮小手段は、前記リーンずれ異常の検出中に前記制御量のガード範囲を縮小する。
【００１４】
好ましくは、前記ガード範囲縮小手段は、ばらつき異常検出開始と同時に前記制御量のガード範囲を縮小し、
前記空燃比制御手段は、ばらつき異常検出開始直前の制御量が縮小されたガード範囲から外れていたとき、その検出開始直前の制御量に最も近いガード値を前記制御量として、前記ばらつき異常検出中の前記補助空燃比制御を実行する。
【００１５】
好ましくは、前記ガード範囲縮小手段は、ばらつき異常検出終了と同時に前記制御量のガード範囲を元の規定のガード範囲に復帰し、
前記空燃比制御手段は、ばらつき異常検出開始直前の制御量を前記制御量として、ばらつき異常検出終了後の前記補助空燃比制御を開始する。
【００１６】
好ましくは、前記異常検出手段は、前記回転変動の検出値が所定の異常判定値以上であるとき、ばらつき異常ありと判定する。
【００１７】
好ましくは、前記回転変動検出手段は、前記内燃機関の気筒毎の回転変動を検出する。
【発明の効果】
【００１８】
本発明によれば、ばらつき異常検出時における誤検出を防止できるという、優れた効果が発揮される。
【図面の簡単な説明】
【００１９】
【図１】本発明の実施形態に係る内燃機関の概略図である。
【図２】触媒前センサおよび触媒後センサの出力特性を示すグラフである。
【図３】空燃比制御のメインルーチンを示すフローチャートである。
【図４】主空燃比補正量の算出マップである。
【図５】補助空燃比補正量の設定のためのサブルーチンを示すフローチャートである。
【図６】触媒後センサ出力差とその積算の様子を示すタイムチャートである。
【図７】補助空燃比補正量の算出マップである。
【図８】リーンずれ異常が発生した場合の各気筒の燃料噴射量の推移を示すグラフである。
【図９】回転変動パラメータを説明するためのタイムチャートである。
【図１０】インバランス割合と回転変動パラメータの関係を示すグラフである。
【図１１】触媒後センサ学習値と回転変動の関係を示すグラフであり、規定のガード範囲を示す。
【図１２】触媒後センサ学習値と回転変動の関係を示すグラフであり、縮小されたガード範囲を示す。
【図１３】ばらつき異常検出を伴う補助空燃比制御のルーチンに関するフローチャートである。
【図１４】図１３のルーチンを実行した際の各値の推移を示すタイムチャートである。
【発明を実施するための形態】
【００２０】
以下、本発明の実施形態を添付図面に基づき説明する。
【００２１】
図１は、本実施形態に係る内燃機関の概略図である。図示されるように、内燃機関（エンジン）１は、シリンダブロック２に形成された燃焼室３の内部で燃料および空気の混合気を燃焼させ、燃焼室３内でピストンを往復移動させることにより動力を発生する。本実施形態の内燃機関１は自動車に搭載された多気筒内燃機関であり、より具体的には直列４気筒火花点火式内燃機関である。内燃機関１は＃１〜＃４気筒を備える。但し気筒数、形式等は特に限定されない。
【００２２】
図示しないが、内燃機関１のシリンダヘッドには吸気ポートを開閉する吸気弁と、排気ポートを開閉する排気弁とが気筒ごとに配設されており、各吸気弁および各排気弁は、カムシャフトを含む動弁機構によって開閉駆動される。シリンダヘッドの頂部には、燃焼室３内の混合気に点火するための点火プラグ７が気筒ごとに取り付けられている。
【００２３】
各気筒の吸気ポートは気筒毎の枝管４を介して吸気集合室であるサージタンク８に接続されている。サージタンク８の上流側には吸気管１３が接続されており、吸気管１３の上流端にはエアクリーナ９が設けられている。そして吸気管１３には、上流側から順に、吸入空気量を検出するためのエアフローメータ５と、電子制御式のスロットルバルブ１０とが組み込まれている。吸気ポート、枝管、サージタンク８及び吸気管１３により吸気通路が形成される。
【００２４】
吸気通路、特に吸気ポート内に燃料を噴射するインジェクタ（燃料噴射弁）１２が気筒ごとに配設されている。インジェクタ１２から噴射された燃料は吸入空気と混合されて混合気をなし、この混合気が吸気弁の開弁時に燃焼室３に吸入され、ピストンで圧縮され、点火プラグ７で点火燃焼させられる。
【００２５】
一方、各気筒の排気ポートは排気マニフォールド１４に接続される。排気マニフォールド１４は、その上流部をなす気筒毎の枝管１４ａと、その下流部をなす排気集合部１４ｂとからなる。排気集合部１４ｂの下流側には排気管６が接続されている。排気ポート、排気マニフォールド１４及び排気管６により排気通路が形成される。
【００２６】
排気管６の上流側と下流側にはそれぞれ三元触媒からなる触媒、すなわち上流触媒１１と下流触媒１９が直列に取り付けられている。上流触媒１１の上流側及び下流側にそれぞれ排気ガスの空燃比を検出するための第１及び第２の空燃比センサ、即ち触媒前センサ１７及び触媒後センサ１８が設置されている。これら触媒前センサ１７及び触媒後センサ１８は、上流触媒１１の直前及び直後の位置に設置され、排気中の酸素濃度に基づいて空燃比を検出する。このように上流触媒１１の上流側の排気合流部に単一の触媒前センサ１７が設置されている。
【００２７】
上述の点火プラグ７、スロットルバルブ１０及びインジェクタ１２等は、制御手段としての電子制御ユニット（以下ＥＣＵと称す）２０に電気的に接続されている。ＥＣＵ２０は、何れも図示されないＣＰＵ、ＲＯＭ、ＲＡＭ、入出力ポート、および記憶装置等を含むものである。またＥＣＵ２０には、図示されるように、前述のエアフローメータ５、触媒前センサ１７、触媒後センサ１８のほか、内燃機関１のクランク角を検出するクランク角センサ１６、アクセル開度を検出するアクセル開度センサ１５、その他の各種センサが図示されないＡ／Ｄ変換器等を介して電気的に接続されている。ＥＣＵ２０は、各種センサの検出値等に基づいて、所望の出力が得られるように、点火プラグ７、スロットルバルブ１０、インジェクタ１２等を制御し、点火時期、燃料噴射量、燃料噴射時期、スロットル開度等を制御する。なおスロットル開度は通常アクセル開度に応じた開度に制御される。
【００２８】
触媒前センサ１７は所謂広域空燃比センサからなり、比較的広範囲に亘る空燃比を連続的に検出可能である。図２に触媒前センサ１７の出力特性を示す。図示するように、触媒前センサ１７は、排気空燃比に比例した大きさの電圧信号Ｖｆを出力する。排気空燃比がストイキ（理論空燃比、例えばＡ／Ｆ＝１４．６）であるときの出力電圧はＶｒｅｆｆ（例えば約３．３Ｖ）である。
【００２９】
他方、触媒後センサ１８は所謂Ｏ2センサからなり、ストイキを境に出力値が急変する特性を持つ。図２に触媒後センサ１８の出力特性を示す。図示するように、排気空燃比がストイキであるときの出力電圧、すなわちストイキ相当値はＶｒｅｆｒ（例えば０．４５Ｖ）である。触媒後センサ１８の出力電圧は所定の範囲（例えば０〜１Ｖ）内で変化する。排気空燃比がストイキよりリーンのとき、触媒後センサの出力電圧はストイキ相当値Ｖｒｅｆｒより低くなり、排気空燃比がストイキよりリッチのとき、触媒後センサの出力電圧はストイキ相当値Ｖｒｅｆｒより高くなる。
【００３０】
上流触媒１１及び下流触媒１９は、それぞれに流入する排気ガスの空燃比Ａ／Ｆがストイキ近傍のときに排気中の有害成分であるＮＯｘ，ＨＣおよびＣＯを同時に浄化する。この三者を同時に高効率で浄化できる空燃比の幅（ウィンドウ）は比較的狭い。
【００３１】
そこで上流触媒１１に流入する排気ガスの空燃比がストイキ近傍に制御されるように、空燃比フィードバック制御（ストイキ制御）がＥＣＵ２０により実行される。この空燃比制御は、触媒前センサ１７によって検出された排気空燃比を所定の目標空燃比であるストイキに一致させるような主空燃比制御（主空燃比フィードバック制御）と、触媒後センサ１８によって検出された排気空燃比をストイキに一致させるような補助空燃比制御（補助空燃比フィードバック制御）とからなる。
【００３２】
ここでストイキ制御について説明する。図３にはストイキ制御のメインルーチンを示す。このメインルーチンはＥＣＵ２０により１エンジンサイクル（＝７２０°ＣＡ）毎に繰り返し実行される。
【００３３】
まずステップＳ１０１では、筒内混合気の空燃比をストイキとするような基本の燃料噴射量即ち基本噴射量Ｑｂが算出される。基本噴射量Ｑｂは例えば、エアフローメータにより検出された吸入空気量Ｇａに基づき、式：Ｑｂ＝Ｇａ／１４．６により算出される。
【００３４】
ステップＳ１０２では触媒前センサ２０の出力Ｖｆが取得される。ステップＳ１０３では、このセンサ出力Ｖｆとストイキ相当センサ出力Ｖｒｅｆｆ（図２参照）との差、即ち触媒前センサ出力差ΔＶｆ＝Ｖｆ−Ｖｒｅｆｆが算出される。
【００３５】
ステップＳ１０４では、この触媒前センサ出力差ΔＶｆに基づき、図４に示したようなマップ（関数でもよい、以下同様）から主空燃比補正量（補正係数）Ｋｆが算出される。触媒前センサ出力差ΔＶｆ及び主空燃比補正量Ｋｆは、主空燃比制御のための制御量をなす。例えばゲインをＰｆとするとＫｆ＝Ｐｆ×ΔＶｆで表される。
【００３６】
そしてステップＳ１０５では、図５に示すサブルーチンで設定された補助空燃比補正量Ｋｒの値が取得される。最後に、ステップＳ１０６にて、インジェクタ２から噴射すべき最終的な燃料噴射量即ち最終噴射量Ｑｆｎｌが式：Ｑｆｎｌ＝Ｋｆ×Ｑｂ＋Ｋｒにより算出される。
【００３７】
図４のマップによれば、触媒前センサ出力Ｖｆがストイキ相当センサ出力Ｖｒｅｆｆより大きい（ΔＶｆ＞０）ほど、即ち実際の触媒前空燃比がストイキからリーン側に離れるほど、１に対しより大きな補正量Ｋｆが得られ、基本噴射量Ｑｂは増量補正される。反対に、触媒前センサ出力Ｖｆがストイキ相当センサ出力Ｖｒｅｆｆより小さい（ΔＶｆ＜０）ほど、即ち実際の触媒前空燃比がストイキからリッチ側に離れるほど、１に対しより小さな補正量Ｋｆが得られ、基本噴射量Ｑｂは減量補正される。こうして、触媒前センサ２０によって検出された空燃比をストイキに一致させるような主空燃比フィードバック制御が実行される。
【００３８】
ステップＳ１０６で得られた最終噴射量Ｑｆｎｌの値は、制御対象となる全気筒に対し一律に用いられる。即ち、１エンジンサイクルの間、最終噴射量Ｑｆｎｌに等しい量の燃料が各気筒のインジェクタ２から順次噴射され、次のエンジンサイクルでは新たに計算された最終噴射量Ｑｆｎｌの燃料が各気筒のインジェクタ２から順次噴射される。
【００３９】
なお、周知のように、最終噴射量Ｑｆｎｌの算出に当たっては他の補正（水温補正、バッテリ電圧補正等）を追加することも可能である。
【００４０】
図５には、補助空燃比補正量の設定のためのサブルーチンを示す。このサブルーチンはＥＣＵ２０により所定の演算周期で繰り返し実行される。
【００４１】
まずステップＳ２０１では、ＥＣＵ２０に装備されたタイマのカウントが実行され、ステップＳ２０２では、触媒後センサ１７の出力Ｖｒが取得される。ステップＳ２０３では、このセンサ出力Ｖｒとストイキ相当センサ出力Ｖｒｅｆｒ（図２参照）との差、即ち触媒後センサ出力差ΔＶｒ＝Ｖｒｅｆｒ−Ｖｒが算出され、この触媒後センサ出力差ΔＶｒが前回積算値に積算される。図６には触媒後センサ出力差ΔＶｒとその積算の様子を示す。
【００４２】
ステップＳ２０４では、タイマ値が所定値ｔｓを超えたか否かが判断される。所定値ｔｓを超えていなければルーチンが終了される。
【００４３】
タイマ値が所定値ｔｓを超えている場合、ステップＳ２０５で、この時点での触媒後センサ出力差積算値ΣΔＶｒが、触媒後センサ学習値ΔＶｒｇとして更新記憶すなわち学習される。そしてステップＳ２０６で、この触媒後センサ学習値ΔＶｒｇに基づき、図７に示したようなマップから、補助空燃比補正量Ｋｒが算出され、この補助空燃比補正量Ｋｒが更新記憶すなわち学習される。触媒後センサ学習値ΔＶｒｇ及び補助空燃比補正量Ｋｒは、補助空燃比制御のための制御量をなす。例えばゲインをＰｒとするとＫｒ＝Ｐｒ×ΔＶｒｇで表される。最後に、ステップＳ２０７にて、触媒後センサ出力差積算値ΣΔＶｒ及びタイマがリセットされる。
【００４４】
触媒後センサ出力差ΔＶｒを所定時間ｔｓの間積算する理由は、触媒後センサ出力Ｖｒのストイキ相当センサ出力Ｖｒｅｆｒに対する時間平均的なズレ量を検知するためである。積算時間を規定する所定値ｔｓは１エンジンサイクルより遙かに長い時間であり、よって触媒後センサ学習値ΔＶｒｇ及び補助空燃比補正量Ｋｒの更新は１エンジンサイクルより長い周期で行われる。
【００４５】
図７のマップによれば、触媒後センサ出力Ｖｒが時間平均的にストイキ相当センサ出力Ｖｒｅｆｒより小さい（ΔＶｒｇ＞０）ほど、即ち実際の触媒後空燃比がストイキからリーン側に離れるほど、０に対しより大きな補正量Ｋｒが得られ、最終噴射量算出の際に基本噴射量Ｑｂは増量補正される。反対に、触媒後センサ出力Ｖｒが時間平均的にストイキ相当センサ出力Ｖｒｅｆｒより大きい（ΔＶｒｇ＜０）ほど、即ち実際の触媒後空燃比がストイキからリッチ側に離れるほど、０に対しより小さな補正量Ｋｒが得られ、基本噴射量Ｑｂは減量補正される。こうして、触媒後センサ２１によって検出された触媒後空燃比をストイキに一致させるような補助空燃比フィードバック制御が実行される。
【００４６】
エンジン及び車両の個体バラツキや触媒前センサ２０の劣化等の理由で、主空燃比フィードバック制御を実行してもその結果がストイキからズレることがあるので、このズレを補正する目的で、補助空燃比フィードバック制御が実行される。
【００４７】
なお、この例では新たな学習値ΔＶｒｇ及び補正量Ｋｒが算出される度にこれらの値自身で更新を行うようにしたが、なまし等の平均化処理を行って更新速度を遅らせるようにしてもよい。
【００４８】
さて、例えば全気筒のうちの一部の気筒（特に１気筒）のインジェクタ１２が故障し、気筒間に空燃比のばらつき（インバランス：imbalance）が発生したとする。例えば＃１気筒が他の＃２、＃３及び＃４気筒よりも燃料噴射量が少なくなり、その空燃比が大きくリーン側にずれる場合等である。このときでも前述の主空燃比フィードバック制御により比較的大きな補正量を与えれば、触媒前センサ１７に供給されるトータルガスの空燃比をストイキに制御できる場合がある。しかし、気筒別に見ると、＃１気筒がストイキより大きくリーン、＃２、＃３及び＃４気筒がストイキよりリッチであり、全体のバランスとしてストイキとなっているに過ぎず、エミッション上好ましくないことは明らかである。そこで本実施形態では、かかる気筒間空燃比ばらつき異常を検出する装置が装備されている。
【００４９】
ここで便宜上、ばらつき度合いの大きさを表すパラメータであるインバランス割合（％）なる用語について説明する。インバランス割合とは、全気筒のうちある１気筒のみが燃料噴射量ズレを起こしている場合に、その燃料噴射量ズレを起こしている気筒（インバランス気筒）の燃料噴射量がどれくらいの割合で、燃料噴射量ズレを起こしていない気筒（バランス気筒）の燃料噴射量即ち基準噴射量からズレているかを示す値である。インバランス割合をＩＢ、インバランス気筒の燃料噴射量をＱｉｂ、バランス気筒の燃料噴射量即ち基準噴射量をＱｓとすると、ＩＢ＝（Ｑｉｂ−Ｑｓ）／Ｑｓで表される。インバランス割合ＩＢが大きいほど、インバランス気筒のバランス気筒に対する燃料噴射量ズレが大きく、空燃比ばらつき度合いは大きい。
【００５０】
図８には、空燃比ばらつき異常として、１気筒のみの燃料噴射量が過少となるリーンずれ異常が発生した場合の各気筒の燃料噴射量の推移を示す。
【００５１】
図８（Ａ）は、異常発生前で且つ各気筒の燃料噴射量が正常な場合の各気筒の燃料噴射量を示す。この場合、前述のストイキ制御の結果、各気筒の燃料噴射量はストイキ相当量Ｑｓｔとなっている。このストイキ相当量Ｑｓｔが前述の基準噴射量Ｑｓであり、＃１，＃２，＃３，＃４気筒のインバランス割合は図示されるように±０％である。
【００５２】
図８（Ｂ）は、異常発生時点における各気筒の燃料噴射量を示す。この場合、＃１気筒の燃料噴射量がストイキ相当量Ｑｓｔより少なく、他の＃２，＃３，＃４気筒の燃料噴射量はストイキ相当量Ｑｓｔとなっている。例えば図示されるように、＃１気筒のインバランス割合は−４０％であり、＃２，＃３，＃４気筒のインバランス割合は±０％である。＃１気筒が、ばらつき異常の原因となっている異常気筒であり、＃２，＃３，＃４気筒は正常気筒である。
【００５３】
この状態で前述のストイキ制御が実行されると、触媒前センサ１７に供給されるトータルガス（全気筒分の排気ガス）の空燃比がストイキになるように、全気筒の燃料噴射量が一律に１０％分ずつ増量補正される。すなわち、＃１気筒のみ４０％のリーンずれが生じても、触媒前センサ１７に対する影響は４０％を気筒数４で除した１０％であり、触媒前センサ１７からは１０％のリーンずれが生じたという情報が得られ、ＥＣＵ２０は全気筒の燃料噴射量を一律に１０％ずつ増量補正する。
【００５４】
この増量補正が開始されてある程度時間が経過すると、各気筒の燃料噴射量は図８（Ｃ）に示す如く変化し、＃１気筒の燃料噴射量はストイキ相当量Ｑｓｔに対し３０％相当少ないリーン状態、他の＃２，＃３，＃４気筒の燃料噴射量はストイキ相当量Ｑｓｔに対し１０％相当多いリッチ状態となる。言い換えれば図示されるように、＃１気筒のインバランス割合は−３０％であり、＃２，＃３，＃４気筒のインバランス割合は＋１０％である。
【００５５】
このように、異常発生時点からある程度の時間を経過しても、気筒間の空燃比ばらつきあるいは不均衡は解消されず、ストイキに対しリッチな気筒とリーンな気筒が混在する。
【００５６】
ばらつき異常の種類として、このようなリーンずれ異常のほか、１気筒のみの燃料噴射量が過大となるリッチずれ異常もある。本実施形態はこれらのうちリーンずれ異常を対象とする。しかしながら本発明はリッチずれ異常に対しても適用可能である。リーンずれ異常は次の方法によりエンジン１の回転変動に基づいて検出される。
【００５７】
図９において、（ａ）はクランク角（°ＣＡ）、（ｂ）は３０°ＣＡ時間Ｔ３０（ｓ）、（ｃ）は回転変動パラメータＹを示す。回転変動パラメータＹとは、エンジン１の回転変動の大きさを表すパラメータのことをいい、以下に述べる方法で検出される。点火順序は＃１、＃３、＃４、＃２の各気筒順である。図中、「正常」とは、いずれの気筒にも空燃比ズレが生じていない正常な場合を示し、「リーンずれ異常」とは、＃１気筒のみにインバランス割合ＩＢ＝−３０（％）のリーンずれ異常が生じている場合を示す。
【００５８】
３０°ＣＡ時間Ｔ３０とは、クランクシャフトが３０°ＣＡ回転するのに要した時間のことをいい、クランク角センサ１６からの出力信号に基づき検出される。３０°ＣＡ時間Ｔ３０が長い（大きい）ほど、クランクシャフトないしエンジンの回転速度は遅い。回転変動パラメータＹとは、前回点火気筒ＴＤＣ（圧縮上死点）から今回点火気筒ＴＤＣまでの間の３０°ＣＡ時間Ｔ３０の差のことをいう。
【００５９】
まずリーンずれ異常の場合に着目する。図示例では＃１気筒にリーンずれ異常が生じているので、（ｂ）に示すように、＃１気筒を点火しても十分な燃焼およびトルクを得られず、次回点火気筒である＃３気筒のＴＤＣの時点において３０°ＣＡ時間Ｔ３０が長くなっている。これに対応して、＃１気筒ＴＤＣから＃３気筒ＴＤＣまでの３０°ＣＡ時間Ｔ３０の差は大きく、（ｃ）に示す回転変動パラメータＹの値は大きくなっている。
【００６０】
これに対し、＃１以外の気筒、例えば＃３気筒の場合は正常であるので、（ｂ）に示すように、＃３気筒ＴＤＣから＃４気筒ＴＤＣまでの３０°ＣＡ時間Ｔ３０の差は小さく、（ｃ）に示す回転変動パラメータＹの値も小さい。
【００６１】
次に正常の場合に着目すると、いずれの気筒も正常であるので、気筒間の燃焼バラツキは少なく、（ｂ）に示すように、前回点火気筒ＴＤＣから今回点火気筒ＴＤＣまでの３０°ＣＡ時間Ｔ３０の差はいずれも小さくほぼ一定で、（ｃ）に示す回転変動パラメータＹの値も常に小さくほぼゼロ付近の値となっている。
【００６２】
つまり、正常の場合には回転変動が小さく、回転変動パラメータＹの値が常に小さいが、異常の場合には、回転変動が大きく、異常気筒に対応した回転変動パラメータＹの値が大きくなる。従って、気筒毎の回転変動パラメータＹを検出し、各気筒の回転変動パラメータＹの検出値を所定の異常判定値と比較することで、リーンずれ異常の有無を検出することが可能である。
【００６３】
本実施形態においては、ＥＣＵ２０が、クランク角センサ１６の出力に基づき気筒毎の回転変動パラメータＹを検出する。そして各気筒の回転変動パラメータＹの検出値が異常判定値以上とならない場合、リーンずれ異常なし、即ち正常と判定する。他方、異常判定値以上となる回転変動パラメータＹの検出値がある場合、リーンずれ異常ありと判定し、当該検出値に対応した気筒を異常気筒と特定する。回転変動パラメータＹの検出値として、複数の検出値の平均値を用いてもよい。
【００６４】
なお、ある１気筒でリーンずれ異常が生じると、当該気筒で爆発後のトルクが十分に出なくなり、回転の立ち上がりが悪化する。こうしたことから気筒毎の爆発後の回転の立ち上がりにムラが生じ、全気筒ないしエンジンサイクル単位で見た場合の回転変動が大きくなる。
【００６５】
図１０には、インバランス割合ＩＢと回転変動パラメータＹの関係を示す。図示されるように、インバランス割合ＩＢと回転変動パラメータＹの間には強い相関性があり、インバランス割合ＩＢが減少（マイナス側に増加）するほど回転変動パラメータＹは増加する。ここで図中のＩＢ２は、異常判定値に相当するインバランス割合ＩＢの値であり、例えば−３０（％）である。
【００６６】
ところで、主空燃比制御および補助空燃比制御のための制御量（以下それぞれ主空燃比制御量および補助空燃比制御量ともいう）については、それぞれ、各値が過剰になるのを防止するため、規定のガード範囲が設けられている。
【００６７】
一方、本発明者の鋭意研究の結果によれば、補助空燃比制御量と回転変動との間に相関性があり、補助空燃比制御量の絶対値が、規定のガード範囲のガード値付近まで大きくなっていると、ばらつき異常検出時に誤検出が発生する虞があることが判明した。
【００６８】
この点を図１１に基づき説明する。図１１においては、横軸に補助空燃比制御量としての触媒後センサ学習値ΔＶｒｇがとってあり、縦軸に回転変動がとってある。
【００６９】
補助空燃比制御において、触媒後センサ学習値ΔＶｒｇがゼロのとき、これによる燃料噴射量補正はなしであり、触媒後センサ学習値ΔＶｒｇがゼロからプラス側に大きくなるほど燃料噴射量はより増量側（リッチ側）に補正され、触媒後センサ学習値ΔＶｒｇがゼロからマイナス側に大きくなるほど燃料噴射量はより減量側（リーン側）に補正される。
【００７０】
他方、縦軸の回転変動は、回転変動パラメータＹの大きさと等価であり、図の上側に行くほど大きくなる。
【００７１】
触媒後センサ学習値ΔＶｒｇに対する規定のガード範囲は、図中の矢印ａで示す範囲であり、マイナス側の下限ガード値ΔＶｒｇ１以上で且つプラス側の上限ガード値ΔＶｒｇ２以下の範囲である。但しΔＶｒｇ１＜０＜ΔＶｒｇ２である。
【００７２】
図中の白丸プロットは、リーンずれ異常が生じていない正常状態のときのデータである。他方、図中の黒四角プロットは、１気筒のみにリーンずれ異常が生じている異常状態のときのデータである。
【００７３】
図示されるように、正常状態のときも異常状態のときも、触媒後センサ学習値ΔＶｒｇが変化すると回転変動も変化し、触媒後センサ学習値ΔＶｒｇがマイナス側からプラス側へと大きくなるにつれ、回転変動は小さくなる傾向にある。
【００７４】
ストイキに対しリッチ側となっている気筒では、燃焼状態が安定し十分なトルクが出るため、気筒毎の回転変動は小さい。他方、ストイキに対しリーン側となっている気筒では、逆に燃焼状態が不安定となり、十分なトルクが出ないため、気筒毎の回転変動が大きい。触媒後センサ学習値ΔＶｒｇが変化すると、系全体がリッチ側またはリーン側に寄るので、燃焼状態が変化し、同一条件でも回転変動が変化する。
【００７５】
正常状態のときに最も大きな回転変動を示すデータｂは、触媒後センサ学習値ΔＶｒｇが下限ガード値ΔＶｒｇ１に等しいときのデータであり、他方、異常状態のときに最も小さな回転変動を示すデータｃは、触媒後センサ学習値ΔＶｒｇが上限ガード値ΔＶｒｇ２に等しいときのデータである。
【００７６】
触媒後センサ学習値ΔＶｒｇの値は、エンジン及び車両の個体バラツキ、触媒前センサ２０の劣化度合い、エンジンの運転状態（定常、過渡等）の違い等により、図示の如く規定のガード範囲ａ内でばらつき得るものである。するとこのばらつきに対応して回転変動もばらつく。
【００７７】
一方、ばらつき異常検出時には、このような触媒後センサ学習値ΔＶｒｇに応じた回転変動のばらつきを考慮することなく、一律に、回転変動パラメータＹの大小のみで異常の有無を判定する。
【００７８】
そこで、触媒後センサ学習値ΔＶｒｇに応じた回転変動のばらつきの影響をできるだけ受けずに異常検出するためには、異常判定値を、データｂ、ｃ間の範囲ｄ内の回転変動に対応した値、好ましくは当該範囲ｄの中間値ｅに設定する必要がある。
【００７９】
しかし、この範囲ｄが比較的狭いため、当該範囲ｄ内の中間値ｅに異常判定値を設定したとしても、誤検出ないし誤判定する虞があった。例えば、データｂが図示の位置から何等かの理由で回転変動大側にずれ、中間値ｅを超えると、本来正常と判定すべきなのに、誤って異常と誤判定してしまう（これを誤異常判定という）。また、データｃが図示の位置から何等かの理由で回転変動小側にずれ、中間値ｅを超えると、本来異常と判定すべきなのに、誤って正常と誤判定してしまう（これを誤正常判定という）。
【００８０】
そこで本実施形態では、図１２に示すように、異常検出中に触媒後センサ学習値ΔＶｒｇのガード範囲を縮小する。すなわち、異常検出中のガード範囲は図中の矢印ａ’で示す範囲であり、この範囲ａ’は規定のガード範囲ａよりも狭い。そしてリーン側の下限ガード値は規定のΔＶｒｇ１からΔＶｒｇ１’に増大変更され、リッチ側の上限ガード値も規定のΔＶｒｇ２からΔＶｒｇ２’に減少変更される。ΔＶｒｇ１＜ΔＶｒｇ１’＜０＜ΔＶｒｇ２’＜ΔＶｒｇ１である。以下、変更後の下限ガード値および上限ガード値をそれぞれ制限下限ガード値および制限上限ガード値という。
【００８１】
すると、正常状態のときに最も大きな回転変動を示すデータは、触媒後センサ学習値ΔＶｒｇが制限下限ガード値ΔＶｒｇ１’に等しいときのデータｆとなり、他方、異常状態のときに最も小さな回転変動を示すデータは、触媒後センサ学習値ΔＶｒｇが制限上限ガード値ΔＶｒｇ２’に等しいときのデータｇとなる。
【００８２】
これらデータｆ，ｇ間の回転変動の範囲ｄ’は、図１１に示した規定の範囲ｄよりも広い。従ってかかる範囲ｄ’内の値、好ましくは当該範囲ｄ’の中間値ｅ’に異常判定値を設定することにより、誤検出すなわち誤異常判定を防止することができる。
【００８３】
すなわち、正常の場合に異常判定値ｅ’に最も近いデータｆは、データｂよりも、異常判定値ｅ’からより離れた小さい値となる。また、異常の場合に異常判定値ｅ’に最も近いデータｇは、データｃよりも、異常判定値ｅ’からより離れた大きい値となる。
【００８４】
よって、正常の場合にも異常の場合にも、触媒後センサ学習値ΔＶｒｇに応じた回転変動のばらつきの影響を受け難くなり、結果として誤検出ないし誤判定を防止することができる。例えば、データｆが図示の位置から何等かの理由で回転変動大側にずれても、異常判定値ｅ’に対するマージンが大きいので、これを超えづらくなり、誤異常判定を防止できる。またデータｇが図示の位置から何等かの理由で回転変動小側にずれても、異常判定値ｅ’に対するマージンが大きいので、これを超えづらくなり、誤正常判定を防止できる。
【００８５】
このように本実施形態によれば、ばらつき異常検出時における誤検出を防止し、検出精度を向上することが可能となる。
【００８６】
上記の考察に従い、本実施形態の異常判定値は、範囲ｄ’の中間値ｅ’に設定される。
【００８７】
次に、図１３および図１４を参照して、ばらつき異常検出を伴う補助空燃比制御の方法をより具体的に説明する。図１３は補助空燃比制御のルーチンに関するフローチャートを示し、図１４は当該ルーチンを実行した際の各値の推移を示す。
【００８８】
図１３のルーチンは、ＥＣＵ２０により、所定の演算周期毎に繰り返し実行されることができる。図１４において、（Ａ）は検出フラグのオンオフ状態を示し、（Ｂ）は学習フラグのオンオフ状態を示し、（Ｃ）は触媒後センサ学習値ΔＶｒｇの値を示し、（Ｄ）はバッファ内に格納された触媒後センサ学習値ΔＶｒｇの値を示す。これらについては後に説明する。
【００８９】
まず図１３に着目して、ステップＳ３０１では、今回の演算タイミングにおける触媒後センサ学習値ΔＶｒｇの値が取得される。
【００９０】
次にステップＳ３０２では、検出フラグがオンか否かが判断される。検出フラグは、ばらつき異常検出が実行されるときにオン、それ以外のときにオフとなるフラグである。検出フラグは、所定の実行条件が成立するとオフからオンに切り替わり、検出が終了するとオンからオフに切り替わる。
【００９１】
検出フラグがオンの場合、ばらつき異常検出が開始または実行され、ステップＳ３０３に進んで学習フラグが強制的にオフされる。学習フラグは、補助空燃比制御に関する学習値である触媒後センサ学習値ΔＶｒｇ（および補助空燃比補正量Ｋｒ）の学習を許可するためのフラグであり、オンのとき学習が許可され、オフのとき学習が禁止される。それ故、学習フラグがオンのときのみ触媒後センサ学習値ΔＶｒｇが更新され、学習フラグがオフのときには更新が停止されると共にオフになる直前の学習値が保持される（図１４（Ｂ）、（Ｃ）参照）。
【００９２】
なお、学習フラグは基本的に図示しない別ルーチンによりオンオフされ、所定条件が成立したときにオン、当該所定条件が非成立になったときにオフされる。ステップＳ３０３では、当該所定条件の成立または非成立に拘わらず、学習フラグが強制的にオフされる。
【００９３】
ステップＳ３０３において学習フラグがオフされるので、ばらつき異常検出中には触媒後センサ学習値ΔＶｒｇの学習および更新が停止される。
【００９４】
次いで、ステップＳ３０４において、今回のタイミングが、検出フラグがオフからオンに切り替わったタイミングであるか否か、言い換えれば、ばらつき異常検出の開始時に該当するか否かが判断される。
判定結果がイエスであるとき、ステップＳ３０５において、前回のタイミングで取得したばらつき異常検出開始直前の触媒後センサ学習値ΔＶｒｇの値が専用のバッファに格納、記憶される。この様子は図１４（Ｃ）、（Ｄ）に示されている。図１４において、時刻ｔ１が今回のタイミングすなわちばらつき異常検出開始時点であるが、この時、その直前の触媒後センサ学習値の値ΔＶｒｇｘが、（Ｄ）に示すように専用のバッファ内に格納される。このように格納されたバッファ内の触媒後センサ学習値の値を「バッファ内ΔＶｒｇ」という。バッファ内ΔＶｒｇの初期値はゼロである。
【００９５】
次いで、ステップＳ３０６において、触媒後センサ学習値ΔＶｒｇのガード範囲が、図１２に示した如く規定のガード範囲から縮小される。すなわち触媒後センサ学習値ΔＶｒｇのガード範囲は、ばらつき異常検出開始と同時に縮小される。
【００９６】
そして、ステップＳ３０７において、縮小後すなわちばらつき異常検出中の補助空燃比制御に用いられる触媒後センサ学習値の値ΔＶｒｇ’が決定される。この決定方法は次の通りである。
【００９７】
まず、ばらつき異常検出開始直前（前回タイミング）の触媒後センサ学習値ΔＶｒｇが縮小後のガード範囲内の値であったとき、その値がそのままばらつき異常検出中の触媒後センサ学習値の値ΔＶｒｇ’として、決定される。
【００９８】
他方、ばらつき異常検出開始直前の触媒後センサ学習値ΔＶｒｇが縮小後のガード範囲から外れた値であったとき、その検出開始直前の触媒後センサ学習値ΔＶｒｇに最も近い制限下限ガード値ΔＶｒｇ１’または制限上限ガード値ΔＶｒｇ２’が、ばらつき異常検出中の触媒後センサ学習値の値ΔＶｒｇ’として、決定される。
【００９９】
図１２を参照して、例えばばらつき異常検出開始直前の触媒後センサ学習値ΔＶｒｇがデータｂで示されるようなΔＶｒｇ１であったとき、これに最も近い制限下限ガード値ΔＶｒｇ１’が、ばらつき異常検出中の触媒後センサ学習値の値ΔＶｒｇ’として用いられる。また、例えばばらつき異常検出開始直前の触媒後センサ学習値ΔＶｒｇがΔＶｒｇ２であったとき、これに最も近い制限上限ガード値ΔＶｒｇ２’が、ばらつき異常検出中の触媒後センサ学習値の値ΔＶｒｇ’として用いられる。
【０１００】
図１４（Ｃ）には前者の例を示す。すなわち、ばらつき異常検出開始直前の触媒後センサ学習値ΔＶｒｇｘが制限下限ガード値ΔＶｒｇ１’より小さい値であったため、ばらつき異常検出開始と同時に、触媒後センサ学習値はより大きい制限下限ガード値ΔＶｒｇ１’に変更される。
【０１０１】
このように、検出開始直前の触媒後センサ学習値が縮小後のガード範囲から外れていたとき、その学習値に最も近いガード値を用いて検出中の補助空燃比制御を行うので、検出中の補助空燃比制御の精度を担保できる。
【０１０２】
図１３に戻って、この後ステップＳ３０８では、ステップＳ３０７で決定された触媒後センサ学習値の値ΔＶｒｇ’を使用して補助空燃比制御が実行される。この決定された値ΔＶｒｇ’は、ばらつき異常検出中に常時使用される。以上でルーチンが終了される。
【０１０３】
他方、ステップＳ３０４の判定結果がノーであるときには、ステップＳ３０５，Ｓ３０６，Ｓ３０７がスキップされ、ステップＳ３０８に進む。
【０１０４】
以上のステップＳ３０２〜Ｓ３０８によれば、検出開始と同時に（ステップＳ３０２：イエス）、学習フラグがオフされ（ステップＳ３０２）、触媒後センサ学習値ΔＶｒｇの更新が停止される。そして開始時点（ステップＳ３０４：イエス）において、その直前の触媒後センサ学習値の値ΔＶｒｇがバッファに記憶され（ステップＳ３０５）、ガード範囲が縮小され（ステップＳ３０６）、新たな触媒後センサ学習値の値ΔＶｒｇ’が決定され（ステップＳ３０７）、このΔＶｒｇ’に基づき補助空燃比制御が実行される（ステップＳ３０８）。検出開始後は、開始時点で決定された触媒後センサ学習値の値ΔＶｒｇ’が保持され、このΔＶｒｇ’に基づき補助空燃比制御が実行される（ステップＳ３０８）。
【０１０５】
他方、ステップＳ３０２において検出フラグがオフの場合、ばらつき異常検出が終了され、ステップＳ３０９に進む。
【０１０６】
ステップＳ３０９においては、今回のタイミングが、検出フラグがオンからオフに切り替わったタイミングであるか否か、言い換えればばらつき異常検出の終了時に該当するか否かが判断される。
判定結果がイエスであるとき、ステップＳ３１０において、学習フラグの強制オフが解除される。これにより触媒後センサ学習値は更新学習可能な元の状態に復帰する。
【０１０７】
そして、ステップＳ３１１において、専用バッファに格納、記憶されていた検出開始直前の触媒後センサ学習値の値が専用バッファから取得される。この様子が図１４（Ｃ）、（Ｄ）に示されている。図１４において、時刻ｔ２が今回のタイミングすなわちばらつき異常検出終了時点であるが、この時、バッファに格納されていた触媒後センサ学習値の値ΔＶｒｇｘがバッファから取り出され、触媒後センサ学習値ΔＶｒｇとして用いられる。
【０１０８】
次いでステップＳ３１２において、ガード範囲が元の規定のガード範囲に復帰される。これによりガード範囲は図１１に示した如く拡大される。
【０１０９】
そしてステップＳ３１３において、通常の触媒後センサ学習値の値ΔＶｒｇを使用して補助空燃比制御が実行される。以上でルーチンが終了される。
【０１１０】
他方、ステップＳ３０９の判定結果がノーであるときには、ステップＳ３１０，Ｓ３１１，Ｓ３１２がスキップされ、ステップＳ３１３に進む。
【０１１１】
以上のステップＳ３０９〜Ｓ３１３によれば、検出終了と同時に（ステップＳ３０９：イエス）、学習フラグの強制オフが解除され（ステップＳ３１０）、学習フラグのオンオフ状態に応じた触媒後センサ学習値ΔＶｒｇの更新学習が再開される。そしてバッファから検出開始直前の触媒後センサ学習値の値が取り出され（ステップＳ３１１）、ガード範囲が復帰、拡大され（ステップＳ３１２）、取り出された触媒後センサ学習値の値から補助空燃比制御が実行される（ステップＳ３１３）。図１４（Ｃ）に示すように、ガード範囲復帰の結果、制限下限ガード値ΔＶｒｇ１’より小さい触媒後センサ学習値の値ΔＶｒｇｘをとることが可能である。
【０１１２】
その後は通常通り、ステップＳ３０１で取得された触媒後センサ学習値の値ΔＶｒｇに基づき補助空燃比制御が実行される（ステップＳ３１３）。
【０１１３】
このように、検出開始時にその直前の触媒後センサ学習値を一旦記憶し、検出終了時にそれを取り出して検出終了後の補助空燃比制御を開始するので、検出終了後の補助空燃比制御の精度を担保できる。
【０１１４】
このルーチンによれば、ばらつき異常検出中に触媒後センサ学習値の学習が強制的に禁止される（ステップＳ３０３）。ばらつき異常検出中には、それまで学習してきた学習値とは異なる別の学習値（ΔＶｒｇ’）に基づき補助空燃比制御が行われる場合がある。よってばらつき異常検出中に学習を禁止することで、別の学習値に基づく不適切な学習を行うことを防止できる。
【０１１５】
以上、本発明の好適な実施形態を詳細に述べたが、本発明の実施形態は他にも様々なものが考えられる。例えば上記実施形態では補助空燃比制御量として触媒後センサ学習値ΔＶｒｇを用いたが、これと等価な補助空燃比補正量Ｋｒを用いてもよい。この場合、触媒後センサ学習値ΔＶｒｇのガード範囲に対応した補助空燃比補正量Ｋｒのガード範囲も存在するので、補助空燃比制御量のガード範囲とは補助空燃比補正量Ｋｒのガード範囲を指す。本発明はリッチずれ異常に対しても適用可能である。リッチずれ異常が生じたときもエンジンの回転変動が大きくなり得るからである。
【０１１６】
本発明の実施形態は前述の実施形態のみに限らず、特許請求の範囲によって規定される本発明の思想に包含されるあらゆる変形例や応用例、均等物が本発明に含まれる。従って本発明は、限定的に解釈されるべきではなく、本発明の思想の範囲内に帰属する他の任意の技術にも適用することが可能である。
【符号の説明】
【０１１７】
１  内燃機関
１１  上流触媒
１７  触媒前センサ
１８  触媒後センサ
１６  クランク角センサ
２０  電子制御ユニット（ＥＣＵ）
【図１】
0005278466.tif 000002
【図２】
0005278466.tif 000003
【図３】
0005278466.tif 000004
【図４】
0005278466.tif 000005
【図５】
0005278466.tif 000006
【図６】
0005278466.tif 000007
【図７】
0005278466.tif 000008
【図８】
0005278466.tif 000009
【図９】
0005278466.tif 000010
【図１０】
0005278466.tif 000011
【図１１】
0005278466.tif 000012
【図１２】
0005278466.tif 000013
【図１３】
0005278466.tif 000014
【図１４】
0005278466.tif 000015