(57)【特許請求の範囲】
【請求項１】
燃料タンクにて発生する蒸発燃料をキャニスタに蓄え、該キャニスタに蓄えられた蒸発燃料を空気と共にキャニスタから内燃機関の吸気側に接続された放出通路を介して放出するようにした内燃機関の空燃比制御装置において、
前記放出通路の経路途中に配設され、前記蒸発燃料を含む空気のパージ量を変化させる流量制御弁と、
前記内燃機関の運転状態に基づいて前記流量制御弁による蒸発燃料及び空気のパージを実行するパージ実行手段と、
パージ実行に伴い前記キャニスタより前記内燃機関の吸気系に放出される蒸発燃料及び空気のパージ量がその時々に許容される限界値としてのパージ限界になるか否かを、パージ非実行時に予測するパージ限界予測手段と、
該パージ限界予測手段による予測結果からパージ限界となる旨が予測された場合、前記パージ実行を禁止するパージ禁止手段と
を具備することを特徴とする内燃機関の空燃比制御装置。
【請求項２】
前記パージ限界予測手段は、少なくともパージ非実行状態からパージ実行状態へ移行するタイミングで前記パージ限界の予測を実行するものである請求項１に記載の内燃機関の空燃比制御装置。
【請求項３】
前記パージ限界予測手段は、前記蒸発燃料及び空気のパージ量が前記内燃機関の運転に支障を来す領域をパージ限界として予測する請求項１又は２に記載の内燃機関の空燃比制御装置。
【請求項４】
前記パージ限界は、電磁式インジェクタが燃料噴射量を制御できる最小通電時間を基に設定される請求項３に記載の内燃機関の空燃比制御装置。
【請求項５】
前記パージ限界は、アイドル回転数制御弁が負担空気量を制御できる最小制御値を基に設定される請求項３に記載の内燃機関の空燃比制御装置。
【請求項６】
多気筒内燃機関に適用され、各気筒に接続された吸気管の集合部に前記蒸発燃料及び空気を放出する内燃機関の空燃比制御装置において、
前記パージ限界は、前記蒸発燃料及び空気のパージ量が気筒分配時に気筒間誤差を生じうるレベルを基に設定される請求項３に記載の内燃機関の空燃比制御装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、燃料タンクに発生する蒸発燃料を一時的に蓄えるキャニスタを備え、当該キャニスタに蓄えられた蒸発燃料を空気と共に内燃機関の吸気系に放出するようにした内燃機関の空燃比制御装置に関するものである。
【０００２】
【従来の技術】
従来より、自動車用内燃機関においては、燃料タンクにて発生する蒸発燃料を一時的に蓄えるキャニスタが設けられており、このキャニスタは放出通路を介して内燃機関の吸気側に接続されている。放出通路の途中には、例えば電磁弁からなる流量制御弁が配設され、この流量制御弁の開度に応じてキャニスタに蓄えられた蒸発燃料が空気と共に放出通路を介して内燃機関に放出されるようになっていた。
【０００３】
また、この種の先行技術文献としては、特公平７−３２１１号公報にて開示された「燃料蒸発ガス排出抑止装置」が知られている。この公報の装置では、排気ガス中の残存酸素濃度を検出し、その検出結果から空燃比を補正するためのフィードバック補正量を算出する。そして、内燃機関に放出する蒸発燃料により生じるフィードバック補正量のズレ量が所定値以上となった時に、キャニスタからの蒸発燃料及び空気（蒸発ガス）のパージ流量を減量し、濃い蒸発ガスが導入されても所定の空燃比に制御することを可能とするものである。
【０００４】
【発明が解決しようとする課題】
ところが、上記従来技術では以下に示す問題が生ずる。つまり、高温環境下或いは揮発性の高い燃料を燃料タンク内に保有している場合等、濃い濃度の蒸発ガスがキャニスタに導入される場合には、パージ非実行状態からパージ実行状態への移行後においてそのパージ実行に伴い蒸発ガス中の多量の燃料が内燃機関に導入されるため（パージ燃料が過多になるため）、その移行後、直ちにパージ流量が減量される可能性がある。かかる場合、パージ非実行からパージ実行への移行時にはパージ流量が所定量に制御されても、その直後にはパージ中止によりパージ流量が「０」になり、さらにその後、再びパージ流量が所定量に制御される。即ち、「パージ流量＞０」→「パージ流量＝０」→「パージ流量＞０」→・・・となるようにパージの実行／非実行が繰り返されることとなる。その結果、空燃比が大きく乱れるようにして変動し、排気エミッションが悪化するなどの問題を招く。
【０００５】
本発明は、上記問題に着目してなされたものであって、その目的とするところは、空燃比を変動なく制御し、ひいては排気エミッションの悪化を抑制することができる内燃機関の空燃比制御装置を提供することである。
【０００６】
【課題を解決するための手段】
上記目的を達成するために、請求項１に記載の発明では、蒸発燃料を含む空気のパージ量を変化させる流量制御弁を備え、内燃機関の運転状態に基づいて前記流量制御弁による蒸発燃料及び空気のパージを実行する（パージ実行手段）。また、パージ実行に伴い前記キャニスタより前記内燃機関の吸気系に放出される蒸発燃料及び空気のパージ量がその時々に許容される限界値としてのパージ限界になるか否かを、パージ非実行時に予測し（パージ限界予測手段）、該パージ限界予測結果からパージ限界となる旨が予測された場合、前記パージ実行を禁止するようにしている（パージ禁止手段）。
【０００７】
要するに、パージ非実行状態からパージを実行する場合において、そのパージ実行後に当該パージが直ちに中断され、さらにその直後にパージが再開されるようなことがあると（パージの開始と停止とが繰り返されるような場合）、空燃比が乱れるという問題が生じる。こうした問題は、パージ実行に伴いキャニスタより内燃機関の吸気系に放出される蒸発燃料及び空気のパージ量が、その時々に許容される限界値に達する際に生じると考えられる。そこで、本発明では、上記許容されるパージ量の限界値を「パージ限界」と定義すると共に、そのパージ限界であるか否かをパージ非実行時に予測し、当該パージ限界と予測される際にはパージ実行を禁止するようにした。かかる場合、パージの開始と停止とが繰り返されるような事態が回避できる。その結果、空燃比を変動なく制御し、ひいては排気エミッションの悪化を抑制することができるという優れた効果が得られる。
【０００８】
ここで、請求項２に記載したように、パージ限界を予測するタイミングは、少なくともパージ非実行状態からパージ実行状態への移行時とするのが望ましい。こうした非パージからパージへの移行時を判定する方法としては、パージ実行を許可するフラグが操作されたことから判定すればよい。
【０００９】
また、請求項３に記載の発明において、前記パージ限界予測手段は、前記蒸発燃料及び空気のパージ量が前記内燃機関の運転に支障を来す領域をパージ限界として予測するようにしている。このとき、内燃機関の運転に支障を来す領域とは、既述したようにパージの実行と停止とが繰り返され、結果として空燃比の乱れを生じさせるような領域を指し、かかる場合においてパージ限界は請求項４〜６に記載したように設定される。
【００１０】
即ち、
・請求項４において、パージ限界は、電磁式インジェクタが燃料噴射量を制御できる最小通電時間を基に設定される。
・請求項５において、パージ限界は、アイドル回転数制御弁が負担空気量を制御できる最小制御値を基に設定される。
・請求項６において、パージ限界は、前記蒸発燃料及び空気のパージ量が気筒分配時に気筒間誤差を生じうるレベルを基に設定される。
【００１１】
【発明の実施の形態】
以下、本発明を具体化した一実施の形態を図面を用いて説明する。
本実施の形態における空燃比制御装置は、周知のマイクロコンピュータ等からなる電子制御装置（Electronic Control Unit ；以下、ＥＣＵという）を中心に構成され、排気ガス中の残存酸素濃度に応じて空燃比フィードバック制御を実施する。またそれに加え、アイドル時における内燃機関への補助空気量を調整するアイドル回転数制御（ＩＳＣ）機能や、燃料タンクにて発生する蒸発燃料（エバポガス）を機関吸気系に放出するパージ制御機能を備える。以下には、その空燃比制御装置の構成及び作用について順を追って説明する。
【００１２】
図１において、内燃機関（エンジン）１は複数のシリンダを有する多気筒ガソリン噴射式内燃機関として構成され、同内燃機関には吸気管２が接続されている。吸気管２にはサージタンク２ａが設けられており、同サージタンク２ａには吸気管２内の圧力である吸気圧ＰＭを検出する吸気圧センサ３が配設されている。サージタンク２ａの上流側には、運転者によるアクセルペダル（図示略）の踏み込み操作に連動するスロットル弁４が配設され、同スロットル弁４の開度（スロットル開度ＴＡ）はスロットルセンサ５により検出されるようになっている。また、吸気管２の最下流側には電磁駆動式のインジェクタ６が配設され、同吸気管２の最上流側にはエアフィルタ７が配設されている。インジェクタ６はＥＣＵ１６からの出力信号に基づき駆動され、内燃機関１に供給される燃料量を調整する。
【００１３】
また、吸気管２には、スロットル弁５を迂回する空気の量を調節するためのＩＳＣ弁２０が配設されている。このＩＳＣ弁２０はＥＣＵ１６からの指令信号により開閉制御され、アイドル時のエンジン回転数が所望の目標値（例えば、７００ｒｐｍ程度）となるようにＩＳＣ弁２０による負担空気量（ＩＳＣ流量）を調整する。
【００１４】
さらに、前記内燃機関１に接続された排気管８には、排気ガス中の酸素濃度に応じてリッチ／リーンのいずれかの電圧信号を出力する酸素（Ｏ２）センサ９が配設されている。
【００１５】
一方、内燃機関１に供給される燃料（ガソリン）を保有する燃料タンク１０には、タンクポート通路１２を介してキャニスタ１１が接続されており、このキャニスタ１１には燃料タンク１０内で発生する蒸発燃料を吸着する吸着剤としての活性炭が収納されている。キャニスタ１１には外気を導入するための大気ポート通路１３が設けられている。また、キャニスタ１１と吸気管２内のスロットル弁４の下流側（サージタンク２ａ）との間は、放出通路１５により接続されており、キャニスタ１１から給送される蒸発ガスは吸気管２の集合部（図示しないインテークマニホールドの上流部）に放出されるようになっている。
【００１６】
放出通路１５の途中には、パージ流量を制御する流量制御弁としてのパージデューティＶＳＶ（Vacuum Switching Valve）１４が配設されている。パージデューティＶＳＶ１４は、スプリング１４ｂにより常時閉側に付勢された可動鉄片からなる弁体１４ｃを有し、コイル１４ａが励磁されることによりスプリング１４ｂの付勢力に抗して弁体１４ｃが開位置に移動する。つまり、コイル励磁時には、弁体１４ｃが当該コイル１４ａに吸引され、負圧通路である放出通路１５が閉から開へと切換えられるようになっている。
【００１７】
かかる場合、上記パージデューティＶＳＶ１４にＥＣＵ１６から制御信号を供給し、それによりキャニスタ１１と前記吸気管２とが放出通路１５を介して連通されるようにしてやれば、大気中から新気が大気ポート通路１３を介して導入される。こうして新気がキャニスタ１１内を換気して内燃機関１の吸気管２内に送り込まれることにより、キャニスタ１１の吸着機能が回復されることになる。このときの新気の導入量に基づくパージ流量Ｑｐ（リットル／ｍｉｎ）は、ＥＣＵ１６からパージデューティＶＳＶ１４に供給されるパルス信号のデューティ比を変えることにより調整される。つまり、パージデューティＶＳＶ１４は、ＥＣＵ１６からのパルス幅変調に基づくデューティ比信号により開度調整され、キャニスタ１１からの蒸発燃料を含む空気のパージ流量Ｑｐを調整する。
【００１８】
ＥＣＵ１６は、入力信号処理回路、演算回路、出力信号回路（駆動回路）及び電源回路等から構成されており、前記吸気圧センサ３、スロットルセンサ５及び酸素センサ９による検出信号を逐次入力する。また、内燃機関１の運転状態を検出するその他のセンサ群として、ディストリビュータ１７には内燃機関１のエンジン回転数ＮＥを検出するための回転数センサ１８が設けられると共に、シリンダブロック１ａにはシリンダ内を循環する冷却水の温度（冷却水温ＴＨＷ）を検出するための水温センサ１９が設けられ、これらセンサ１８，１９による検出信号も併せてＥＣＵ１６に入力される。なお、この他にもエアフィルタ７から吸入される空気温度を検出する吸気温センサ等を必要に応じて配設し、その検出信号をＥＣＵ１６に入力するようにしてもよい。
【００１９】
また、ＥＣＵ１６は、前記酸素センサ９からの電圧信号に基づく混合気のリッチ／リーン判定によりフィードバック補正係数ＦＡＦを求め、このＦＡＦ値を用いて空燃比フィードバック制御を実施する。さらに、ＥＣＵ１６は機関運転状態により基本噴射時間ＴＰを求めると共に、基本噴射時間ＴＰに対しＦＡＦ値等による補正を行って最終噴射時間ＴＡＵを求め、このＴＡＵ値を用い所定の噴射タイミングで前記インジェクタ６により燃料噴射を行わせる。
【００２０】
次に、本実施の形態にかかる内燃機関の空燃比制御装置の作用を説明する。本実施の形態では、ＥＣＵ１６内の複数の演算プログラムを操作することにより、下記の空燃比フィードバック制御、エバポ濃度演算、燃料噴射制御、アイドル回転数制御及びデューティＶＳＶ制御を行うこととしており、以下、各演算プログラムを個々に説明する。
【００２１】
［空燃比フィードバック制御］
空燃比フィードバック制御ルーチンを図２のフローチャートに基づいて説明する。なお、この空燃比フィードバック制御ルーチンは約４ｍｓｅｃ（ミリ秒）毎にＥＣＵ１６にて実行される。
【００２２】
さて、図２のルーチンがスタートすると、ＥＣＵ１６は、先ずステップ１０１でフィードバック（Ｆ／Ｂ）制御条件が成立しているか否かを判定する。このＦ／Ｂ条件は、主として以下に示す条件を全て満足した場合に成立する。
（１）始動時でないこと。
（２）燃料カット中でないこと。
（３）冷却水温ＴＨＷが所定温度以上であること。
（４）酸素センサ９が活性状態であること。
（５）高負荷及び高回転状態でないこと。
【００２３】
上記Ｆ／Ｂ条件が成立すると、ＥＣＵ１６はステップ１０２に進む。そして、ＥＣＵ１６は、ステップ１０２で酸素センサ９の電圧信号と所定判定レベルとを比較し、それぞれ遅延時間Ｈ，Ｉ（ｍｓｅｃ）にて空燃比フラグＸＯＸＲを操作する。このとき、酸素センサ９の電圧信号が空燃比リッチを示すものであれば、空燃比フラグＸＯＸＲに「１」がセットされ、電圧信号が空燃比リーンを示すものであれば空燃比フラグＸＯＸＲが「０」にクリアされる。
【００２４】
次に、ＥＣＵ１６はステップ１０３に進み、空燃比フラグＸＯＸＲに基づいてＦＡＦ値を操作する。即ち、空燃比フラグＸＯＸＲが「０」→「１」、又は「１」→「０」に変化した時、ＦＡＦ値を所定量変化（スキップ）させ、空燃比フラグＸＯＸＲが「１」又は「０」を継続中の時はＦＡＦ値を積分制御する（徐々に増減させる）。その後、ＥＣＵ１６は、ステップ１０４でＦＡＦ値の上下限チェックをした後ステップ１０５に進み、上記の如く決定したＦＡＦ値を基にしてスキップ毎又は所定時間毎になまし（平滑化）処理を行い、ＦＡＦ値の平均値であるＦＡＦＡＶ値を算出する。その一例として、ＦＡＦＡＶ値は次の（１）式にて算出される。
【００２５】
ＦＡＦＡＶ＝｛ＦＡＦＡＶi-1 ・ｎ＋ＦＡＦ・（２５６−ｎ）｝／２５６・・・（１）
ここで、添字「ｉ−１」は、なまし値ＦＡＦＡＶの前回値であることを示し、符号ｎは、なまし度合を決定するための定数である。
【００２６】
なお、ステップ１０１でＦ／Ｂ制御が成立しない時には、ＥＣＵ１６はステップ１０６に進み、ＦＡＦ値を「１．０」とする。そして、ＥＣＵ１６は、本ルーチンを一旦終了する。因みに、ＦＡＦ値は理論空燃比（１４．７）からどれだけ外れているかを示す指標である。
【００２７】
［エバポ濃度演算］
次に、蒸発ガスの濃度であるエバポ濃度を算出するためのエバポ濃度演算ルーチンを図３のフローチャートに基づいて説明する。なお、このエバポ濃度演算ルーチンは約１６ｍｓｅｃ毎にＥＣＵ１６にて実行される。
【００２８】
本ルーチンがスタートすると、ＥＣＵ１６は、先ずステップ２０１でパージ実行許可フラグＸＰＲＧ＝「１」であり、且つパージデューティＶＳＶ１４が出力中であるか否かを判別する。ここで、パージ実行許可フラグＸＰＲＧは、パージ実行を許可するか否かを表すフラグであって、ＸＰＲＧ＝１はパージ実行を許可する旨を表し、ＸＰＲＧ＝０はパージ実行を許可しない旨を表す。また、パージデューティＶＳＶ１４が出力中であることは、当該ＶＳＶ１４が駆動され、蒸発ガスのパージが実際に行われていることを意味する。
【００２９】
そして、ステップ２０１が否定判別されれば、即ちパージ非実行中であれば、ＥＣＵ１６はそのままルーチンを終了する。また、ステップ２０１が肯定判別されれば、即ちパージ実行中であれば、ＥＣＵ１６は、ステップ２０２でＦＡＦ値のなまし値であるＦＡＦＡＶ値からＦＡＦ値の基準値（＝１）を減算した値（＝ＦＡＦＡＶ−１）を求め、その後、ステップ２０３〜２０７で（ＦＡＦＡＶ−１）値に応じてエバポ濃度ＦＧＰＧを演算する。
【００３０】
詳しくは、ＥＣＵ１６は、ステップ２０３で（ＦＡＦＡＶ−１）＞２％であるか否か、即ち空燃比がリーン寄りであるか否かを判別する。また、ＥＣＵ１６は、ステップ２０４で（ＦＡＦＡＶ−１）＜−２％であるか否か、即ち空燃比がリッチ寄りであるか否かを判別する。そして、（ＦＡＦＡＶ−１）＞２％であれば、即ち空燃比がリーン寄りであれば、ＥＣＵ１６は現在のＦＧＰＧ値よりも実際のＦＧＰＧ値は薄いと判断し、ステップ２０５でＦＧＰＧ値を所定量（本実施の形態では、０．２％）だけ減少させる。（ＦＡＦＡＶ−１）＜−２％であれば、即ち空燃比がリッチ寄りであれば、ＥＣＵ１６は現在のＦＧＰＧ値よりも実際のＦＧＰＧ値は濃いと判断し、ステップ２０６でＦＧＰＧ値を所定量（本実施の形態では、０．２％）だけ増加させる。また、−２％≦（ＦＡＦＡＶ−１）≦２％であれば、ＥＣＵ１６は現在のＦＧＰＧ値がほぼ実際値に合致すると判断し、ステップ２０７でＦＧＰＧ値をその時の値にホールドする。
【００３１】
ＦＧＰＧ値の演算後、ＥＣＵ１６は、ステップ２０８にてＦＧＰＧ値が上下限値である０〜２５％以内であるか否かをチェックし、その後本ルーチンを終了する。
【００３２】
［燃料噴射制御］
燃料噴射制御ルーチンを図４のフローチャートに基づいて説明する。なお、この燃料噴射制御ルーチンも約４ｍｓｅｃ毎にＥＣＵ１６にて実行され、本ルーチンによりインジェクタ６の駆動が制御されるようになっている。
【００３３】
図４において、ＥＣＵ１６は、先ずステップ３０１でエンジン回転数ＮＥとエンジン負荷（例えば、吸気圧ＰＭ）とに基づいてインジェクタ６の基本噴射時間ＴＰを算出し、続くステップ３０２で前記ＴＰ値に対して各種基本補正（冷却水温、始動後吸気温等による補正）を実行する。次に、ＥＣＵ１６はステップ３０３に進み、前記図３のルーチンにより算出したエバポ濃度ＦＧＰＧと、後述する図６及び図７のルーチンにより算出されるパージ率ＰＧＲとを乗算してパージ補正係数ＦＰＧを算出する（ＦＰＧ＝ＦＧＰＧ・ＰＧＲ）。
【００３４】
最後に、ＥＣＵ１６は、ステップ３０４でＦＡＦ値、ＦＰＧ値を次式、
１＋（ＦＡＦ−１）＋ＦＰＧ
に代入してこれを補正係数として算出し、その補正係数をインジェクタ６による最終噴射時間ＴＡＵに反映させる。
【００３５】
［アイドル回転数制御］
次に、アイドル回転数制御ルーチンを図５のフローチャートに基づいて説明する。なお、このアイドル回転数制御ルーチンは約１６ｍｓｅｃ毎にＥＣＵ１６にて実行され、本ルーチンによりＩＳＣ弁２０の駆動が制御されるようになっている。
【００３６】
図５において、ＥＣＵ１６は、先ずステップ４０１でエンジン冷却水温ＴＨＷに基づいてＩＳＣ弁２０による基本負担空気量ＧＢを算出する。さらに、ＥＣＵ１６は、続くステップ４０２で前記ＧＢ値に対して各種基本補正（大気圧、吸気温等による補正）を行い、それにより基本ＩＳＣ流量ＱＢを算出する。このとき、前記ステップ４０１にて算出されたＧＢ値が質量空気量であるのに対し、前記ステップ４０２にて算出されたＱＢ値は体積空気量となっている。
【００３７】
次に、ＥＣＵ１６は、ステップ４０３で前記算出した基本ＩＳＣ流量ＱＢから後述する図６及び図７のルーチンにより算出されるパージ流量Ｑｐを減算してＩＳＣ流量Ｑｉを算出する（Ｑｉ＝ＱＢ−Ｑｐ）。
【００３８】
最後に、ＥＣＵ１６は、ステップ４０４で図１１に示すＩＳＣ流量Ｑｉ（リットル／ｍｉｎ）とデューティ比（％）との関係を示すマップを用い、ＩＳＣ負担空気量（Ｑｉ）に対応するＩＳＣ弁２０の駆動デューティ比を算出する。そして、該算出した駆動デューティ比によりＩＳＣ弁２０を駆動させる。なお、図１１のマップは、ＩＳＣ弁２０の駆動周波数（Ｈｚ）をパラメータとして実験的に求められるようになっている。
【００３９】
図１１のマップにおいて、横軸のデューティ比が略２０％を越える領域では、ＩＳＣ流量Ｑｉとデューティ比との関係が直線的に安定するが、デューティ比が略２０％以下の領域ではＱｉ値とデューティ比との関係が安定しないことが分かる。ここで、Ｑｉ値とデューティ比との関係が不安定になるＱｉ値が、ＩＳＣ弁２０の最小制御流量ＱｉMIN となっている。
【００４０】
［デューティＶＳＶ制御］
次に、本実施の形態の要旨であるところのデューティＶＳＶ制御ルーチンについて、図６及び図７を用いて説明する。このデューティＶＳＶ制御ルーチンは１００ｍｓｅｃ毎の時間割り込みによりＥＣＵ１６にて実行され、本ルーチンにより蒸発ガスのパージを実行するためのパージデューティＶＳＶ１４が駆動されるようになっている。
【００４１】
先ず最初にその概要を説明すれば、本ルーチンでは、蒸発ガスのパージが機関運転状態に支障を来たすおそれのある領域を「パージ限界」として定義し、機関運転状態がパージ限界に突入しているか否かに応じてパージ率ＰＧＲを増減させるようにしている。また他方で、非パージ状態からパージ状態に移行するに際し、パージ以後の運転状態がパージ限界に突入する可能性があるのか否かを予測するようにしている。
【００４２】
つまり、本実施の形態では、パージ実行時においてパージ限界である旨が判断されると、パージ率ＰＧＲを減少させてパージ量を少なくし、パージ限界でない旨が判断されると、パージ率ＰＧＲを増加させてパージ量を多くする。また、非パージ状態において、パージの実行に伴いパージ限界に突入する可能性があるかどうかを予測し、パージ限界に突入する可能性がある場合には直ぐにパージを実行せず、パージ限界に突入する可能性が無くなってからパージを実行するようにしている。
【００４３】
ここで、本実施の形態では、次の（イ）〜（ハ）に示す各点をパージ限界点として予め設定しており、それを基準にパージ限界か否かを判断、又は予測するようにしている。
【００４４】
（イ）インジェクタ６の噴射パルス幅が制御可能な最小値（ＴＡＵMIN ）付近となる点をパージ限界点とする。そして、インジェクタ６による燃料噴射時間ＴＡＵが最小通電時間ＴＡＵMIN に安全を見込んだ余裕値Ａ（但し、Ａ＝０でも可）を加算した値よりも小さい場合に（即ち、ＴＡＵ＜ＴＡＵMIN ＋Ａの場合に）、パージ限界であると判断、又は予測する。この場合、最小通電時間ＴＡＵMIN は、図１０に示すインジェクタ６の燃料噴射特性を基に設定されるものであって、同図によればＴＡＵ値がＴＡＵMIN 以下では直線的な噴射量特性が得られないことが分かる。なお、図中のＴｖはインジェクタ６の無効噴射時間を表す。
【００４５】
（ロ）ＩＳＣ弁２０の負担空気量が制御可能な最小値（ＱｉMIN ）付近となる点をパージ限界点とする。そして、ＩＳＣ流量Ｑｉが最小制御流量ＱｉMIN に安全を見込んだ余裕値Ｂ（但し、Ｂ＝０でも可）を加算した値よりも小さい場合に（即ち、Ｑｉ＜ＱｉMIN ＋Ｂの場合に）、パージ限界であると判断、又は予測する。この場合、最小制御流量ＱｉMIN は、図１１に示すＩＳＣ特性を基に設定されるものであって、同図によればＱｉ値がＱｉMIN 以下では直線的なＩＳＣ特性が得られないことが分かる。
【００４６】
（ハ）内燃機関１に吸入される蒸発ガスの量が多いか又はその濃度が濃く、蒸発ガスの気筒分配時に気筒間でパージ量（蒸発ガスの吸入量）に差を生じるレベルとなる点をパージ限界点とする。具体的には、パージ補正係数ＦＰＧが所定値ＦＫ以上となる場合に、パージ限界であると判断する。この場合、所定値ＦＫは、事前に実施された実験値により設定されるようになっている。
【００４７】
以下、上記（イ）〜（ハ）の判断基準を前提として、図６及び図７のデューティＶＳＶ制御ルーチンを説明する。
さて、本ルーチンがスタートすると、ＥＣＵ１６は、先ずステップ５０１で機関運転状態を読み込む。具体的には回転数センサ１８による検出信号から求められたエンジン回転数ＮＥと、吸気圧センサ３による検出信号から求められた吸気圧ＰＭとを読み込む。併せて、ＥＣＵ１６は、前記吸気圧ＰＭに基づいて内燃機関１に供給される混合気の吸入空気量Ｑａを算出する。
【００４８】
次に、ＥＣＵ１６は、ステップ５０２で所定のパージ条件が成立しているか否かを判別する。ここで、パージ条件とは、前記ステップ５０１で読み込んだ機関運転状態に基づくものであって、エンジン回転数ＮＥが所定回転数以上であり、且つ吸入空気量Ｑａが所定量以上である場合に成立する。そして、パージ条件が不成立であれば、ＥＣＵ１６はステップ５０１に戻り同様の処理を繰り返し実行する。
【００４９】
前記ステップ５０２のパージ条件が成立していれば、ＥＣＵ１６はステップ５０３に進み、パージ率ＰＧＲの初期値として運転状態に影響を与えない程度のα％（例えば、０又は０．１％程度の小さな値）を設定する。次に、ＥＣＵ１６は、ステップ５０４でその時のパージ率ＰＧＲと吸入空気量Ｑａとを乗算してパージ流量Ｑｐを算出する（Ｑｐ＝ＰＧＲ・Ｑａ）。さらに、ＥＣＵ１６は、続くステップ５０５で図９に示すマップに基づいてパージ流量Ｑｐ（リットル／ｍｉｎ）に対応するパージデューティＶＳＶ１４の駆動デューテイ比（％）を算出する。なお、図９のマップは、パージデューティＶＳＶ１４が配設された放出通路１５の前後差圧（ｍｍＨｇ）とパージデューティＶＳＶ１４の駆動周波数（Ｈｚ）とをパラメータとして実験的に求められるようになている。また、図９のマップにおいて、デューティ比が略１５〜２０％内の所定値γ（ＱｐMIN に相当するデューティ比）を越えると、パージ流量Ｑｐとデューティ比とが直線的に安定して増減するが、デューティ比が所定値γ以下の領域ではＱｐ値とデューティ比との関係が安定しないことが分かる。
【００５０】
次に、ＥＣＵ１６は、ステップ５０６で前記算出したデューティ比が前記図９の所定値γに対してそれ以上であるか否かを判別する。そして、ステップ５０６でデューティ比がγ未満である旨が判別されると、ＥＣＵ１６は、蒸発ガスのパージ実行が内燃機関１の挙動に悪影響を与える可能性があるとみなし、ステップ５０７に進んでパージ実行許可フラグＸＰＲＧを「０」にクリアする。例えば機関始動当初には、ステップ５０６が否定判別され、かかる場合には、パージデューテイＶＳＶ１４を駆動するためのパルス信号を出力することなく、後述するステップ５０８〜５１０が実行される。
【００５１】
つまり、ステップ５０８に進んだＥＣＵ１６は、今現在のパージ実行状態がパージ限界に達しているか否かを判別する。ここで、ステップ５０８では前記（イ）〜（ハ）の判断基準を用い、その時の最終噴射時間ＴＡＵ、ＩＳＣ流量Ｑｉ、パージ補正係数ＦＰＧ等に基づいてパージ限界に突入しているか否かが判断される。即ち、
・ＴＡＵ＜ＴＡＵMIN ＋Ａ、又は、
・Ｑｉ＜ＱｉMIN ＋Ｂ、又は、
・ＦＰＧ＞ＦＫ
が成立する場合には、パージ限界である旨が判別される。
【００５２】
このとき、パージ限界に突入していない旨が判別されれば、ＥＣＵ１６はステップ５０９に進み、パージ率ＰＧＲを所定値β（％）だけ増加させる。また、パージ限界に突入している旨が判別されれば、ＥＣＵ１６はステップ５１０に進み、パージ率ＰＧＲを所定値β（％）だけ減少させる。こうしたパージ率ＰＧＲの設定後、ＥＣＵ１６はステップ５０４に戻り、既述した処理を繰り返し実行する。
【００５３】
一方、前記ステップ５０６でパージデューティ比がγ以上である旨が判別されると、ＥＣＵ１６は、ステップ５１１でパージ実行許可フラグＸＰＲＧに「１」をセットする。即ち、パージデューティＶＳＶ１４が安定領域で駆動できるとみなされ（図９において、デューティ比≧γの領域）、パージ実行が許可される。パージ実行許可フラグＸＰＲＧのセット後、ＥＣＵ１６は図７のステップ５１２に進む。
【００５４】
その後、ＥＣＵ１６は、ステップ５１２でそのタイミングが非パージ状態からパージ状態に移行したタイミングであるか否かを判別する。具体的には、パージ実行許可フラグＸＰＲＧが「０」から「１」に操作されたタイミングであるか否かを判別する。
【００５５】
このとき、非パージ→パージの移行時でなければ（パージ継続状態である場合）、ＥＣＵ１６はステップ５１２を否定判別してステップ５１３に進み、前記算出したデューティ比に対応するパルス信号をパージデューティＶＳＶ１４に出力する。即ち、蒸発ガスのパージを実行する。
【００５６】
その後、ＥＣＵ１６は、ステップ５１４で機関運転条件（ＮＥ，ＰＭ等）を読み込むと共に、続くステップ５１５でパージ実行条件が成立しているか否かを判別する（既述したステップ５０１，５０２に同じ）。この場合、パージ条件が成立していれば、ＥＣＵ１６は前記図６のステップ５０８に進み、ステップ５０８〜５１０でパージ率ＰＧＲの更新処理を実行する。つまり、ステップ５０６が肯定判別され、且つステップ５１２が否定判別された場合には、通常のパージ処理が実行されることとなる。なお、パージ条件が不成立であれば図６のステップ５０１に戻る。
【００５７】
また、前記ステップ５１２で非パージ→パージの移行時である旨が判別されれば（ＸＰＲＧの「０」→「１」操作時）、ＥＣＵ１６はステップ５３０に進んで後述するパージ限界予測ルーチン（図８）を実行し、パージデューティＶＳＶ１４によるパージ実行に伴い、パージ限界に突入する可能性があるか否かを予測する。このステップ５３０の処理では、パージ補正係数ＦＰＧ、燃料噴射時間ＴＡＵ、パージ流量Ｑｐ、ＩＳＣ流量Ｑｉ等の情報からパージ限界に突入しているか否かが予測され、その予測結果に応じてパージ限界予測フラグが「１」又は「０」に操作される。このとき、パージ限界予測フラグ＝１は、パージを実行するとパージ限界に突入する可能性があることを表し、パージ限界予測フラグ＝０は、パージを実行してもパージ限界に突入する可能性がないことを表す。
【００５８】
その後、ＥＣＵ１６は、ステップ５２１でパージ限界予測フラグが「１」であるか否かを判別する。そして、パージ限界予測フラグ＝０の時、即ちパージを実行してもパージ限界に突入する可能性がないと予測された場合、ＥＣＵ１６はステップ５１３に進み、パージデューティＶＳＶ１４を駆動させてパージを実行する。そして、それ以降、既述したようにステップ５１３→５１４→５１５→…の処理を実行する。
【００５９】
また、前記ステップ５２１でパージ限界予測フラグ＝１の時、即ちパージを実行するとパージ限界に突入する可能性があると予測された場合、ＥＣＵ１６はステップ５２２に進み、機関運転状態（ＮＥ，ＰＭ等）を読み込む。さらに、ＥＣＵ１６は、ステップ５２３でパージ流量Ｑｐを算出すると共に、続くステップ５２４でパージ流量Ｑｐに応じたデューティ比を算出する（前記図６のステップ５０４，５０５に同じ）。その後、ＥＣＵ１６は、ステップ５２５でパージ条件が成立しているか否かを判別する（前記ステップ５０２に同じ）。このとき、パージ条件が成立していれば、ＥＣＵ１６はステップ５３０に戻って前述と同様の処理を繰り返し、パージ条件が不成立であれば図６のステップ５０１に戻って前述と同様の処理を繰り返す。
【００６０】
つまり、上記ステップ５１２が肯定判別された際には、パージ限界突入の可能性がないと予測された場合にのみ（パージ限界予測フラグ＝０の場合のみ）、ステップ５１３に進むことが許可され、パージデューティＶＳＶ１４による蒸発ガスのパージが実行されることになる。
【００６１】
次に、非パージ→パージの移行時に（前記図６のステップ５１２が肯定判別された以降に）実行されるステップ５３０のパージ限界予測ルーチンについて、図８のフローチャートを用いて説明する。なお、本予測ルーチンにおいて、パージ限界の予測は、既述した（イ）〜（ハ）の判断基準を前提として行われる。図８のルーチンを大別すると、同図のステップ５３１〜５３５がインジェクタ６の燃料噴射特性に基づくパージ限界の予測処理に、ステップ５３６〜５３８がＩＳＣ特性に基づくパージ限界の予測処理に、また、ステップ５３９〜５４１がパージ補正係数ＦＰＧに基づくパージ限界の予測処理に、それぞれ相当する。
【００６２】
さて、図８のルーチンがスタートすると、ＥＣＵ１６は、先ずステップ５３１で非パージ時におけるインジェクタ６の燃料噴射時間ＴＡＵ１を算出する。このとき、燃料噴射時間ＴＡＵ１は、パージによる燃料加増がないものとして（パージ補正係数ＦＰＧ＝０として）、
ＴＡＵ１＝ＴＰ・｛１＋（ＦＡＦ−１）｝
といった形態で算出される。
【００６３】
また、ＥＣＵ１６は、続くステップ５３２でパージ実行により加増されると想定される燃料噴射時間ＴＡＵ２を算出する。このとき、燃料噴射時間ＴＡＵ２は、その時のエパポ濃度ＦＧＰＧとパージ率ＰＧＲとに応じたパージ補正係数ＦＰＧ（＝ＦＧＰＧ・ＰＧＲ）を用いて、
ＴＡＵ２＝ＴＰ・ＦＰＧ
といった形態で算出される。
【００６４】
その後、ＥＣＵ１６は、ステップ５３３で前記燃料噴射時間ＴＡＵ１とＴＡＵ２との加算値がインジェクタ６による制御可能な最小値（図１０の最小通電時間ＴＡＵMIN ）よりも大きいか否かを判別する。そして、ＴＡＵ１＋ＴＡＵ２＞ＴＡＵMIN であれば、ＥＣＵ１６は、ステップ５３４でパージ限界予測フラグとしてのＸＦ１を「０」とし、ＴＡＵ１＋ＴＡＵ２≦ＴＡＵMIN であれば、ステップ５３５で同ＸＦ１を「１」とする。
【００６５】
また、ＥＣＵ１６は、ステップ５３６でその時のＩＳＣ流量Ｑｉからパージ流量Ｑｐを減算した値が、ＩＳＣ特性上で安定した直線域を確保するための最小値（図１１の最小制御流量ＱｉMIN ）よりも大きいか否かを判別する。そして、Ｑｉ−Ｑｐ＞ＱｉMIN であれば、ＥＣＵ１６は、ステップ５３７でパージ限界予測フラグとしてのＸＦ２を「０」とし、Ｑｉ−Ｑｐ≦ＱｉMIN であれば、ステップ５３８で同ＸＦ２を「１」とする。
【００６６】
さらに、ＥＣＵ１６は、ステップ５３９でその時のパージ補正係数ＦＰＧが所定値ＦＫよりも大きいか否かを判別する。そして、ＦＰＧ＞ＦＫであれば、ＥＣＵ１６は、ステップ５４０でパージ限界予測フラグとしてのＸＦ３を「０」とし、ＦＰＧ≦ＦＫであれば、ステップ５４１で同ＸＦ３を「１」とする。
【００６７】
以上のように本図８のルーチンでは、ステップ５３３，５３６，５３９の３つの判定条件によりパージ限界予測フラグＸＦ１，ＸＦ２，ＸＦ３が操作される。そして、既述した図７のステップ５２１では、上記パージ限界予測フラグＸＦ１，ＸＦ２，ＸＦ３の全てが「０」であれば同ステップが否定判別され、パージデューティＶＳＶ１４によるパージが開始される。また、ＸＦ１，ＸＦ２，ＸＦ３が１つでも「１」であれば同ステップが肯定判別され、パージ限界の予測処理が繰り返し実行される。つまり、パージ非実行状態になり、且つパージ限界である旨が予測される場合には、パージ限界になる条件が解除されるまで（ＸＦ１，ＸＦ２，ＸＦ３が全てクリアされるまで）、パージの実行が禁止されることになる。
【００６８】
なお、本実施の形態では、図６，７のステップ５０８→５０９（５１０）→５１２→５１３→…といった流れの処理が請求項記載のパージ実行手段に、ステップ５３０（図８）が請求項記載のパージ限界予測手段に、また、ステップ５３０→５２１→５２２→…といった流れの処理（フラグＸＦ１，ＸＦ２，ＸＦ３の操作及び判定処理）が請求項記載のパージ禁止手段に、それぞれ相当する。
【００６９】
図１２は、上記の如く実施される制御概要を示すタイムチャートであり、以下には同図を用いて本実施の形態のパージ制御をより具体的に説明する。なお、図１２中、時間ｔ１，ｔ４が前記図７のステップ５１２で言う「非パージ→パージの移行時」に相当し、この時間ｔ１，ｔ４で前記図７のステップ５３０の処理（図８のパージ限界予測処理）が実行される。但し、時間ｔ１ではパージ限界突入の条件が成立しないが、時間ｔ４ではパージ限界突入の条件が成立するものとしている（即ち、時間ｔ１ではパージ限界予測フラグ＝０、時間ｔ４ではパージ限界予測フラグ＝１となっている）。
【００７０】
図１２において、パージ率ＰＧＲが徐々に上昇し（前記図６のステップ５０９による）、それに伴なって時間ｔ１でパージデューティ比が所定値γを超えると、パージ実行許可フラグＸＰＲＧに「１」がセットされ（図６のステップ５１１）、パージ実行が許可される。この時間ｔ１では、パージを実行してもパージ限界突入の可能性がないため（パージ限界予測フラグ＝０）、ＸＰＲＧの操作に従いパージデューティＶＳＶ１４の駆動によるパージ処理が実行される。
【００７１】
その後、時間ｔ２になると、パージ限界に突入した旨が判別され、パージ率ＰＧＲが徐々に下降する（図６のステップ５１０）。パージ率ＰＧＲの下降に伴ってパージデューティ比が所定値γを下回る時間ｔ３では、パージ実行許可フラグＸＰＲＧが「０」にクリアされ、パージデューティＶＳＶ１４によるパージ処理が非実行状態となる。
【００７２】
さらにその後、再びパージデューティ比が上昇し、同デューティ比が所定値γを超える時間ｔ４では、パージ実行許可フラグＸＰＲＧに「１」がセットされる。また、この時間ｔ４では、パージの実行によりパージ限界突入の可能性があるため（パージ限界予測フラグに「１」がセットされる）、パージ実行許可フラグＸＰＲＧがセット状態のまま保持されていてもパージデューティＶＳＶ１４によるパージ処理が実施されることはない。
【００７３】
そして、パージを実行してもパージ限界突入の可能性がないと判断される時間ｔ５になると、パージ限界予測フラグが「０」にクリアされ、この時間ｔ５でパージデューティＶＳＶ１４の駆動によるパージ処理が再開される。
【００７４】
以上のように本実施の形態の制御装置では、パージ限界予測処理を実施するようにしたため、パージ実行によりパージ限界に突入すると思われる状態下（時間ｔ４〜ｔ５）ではパージ処理が中断される。これに対して、パージ限界予測処理を実施しない従来装置の推移を比較例として破線で示すと、時間ｔ４〜ｔ５においてパージの実行と非実行とが繰り返され、結果として空燃比が乱れ、排気エミッションの悪化を招くことが分かる。
【００７５】
以上詳述したように本実施の形態によれば、以下に示す効果が得られる。
（ａ）本実施の形態では、パージ実行に伴いキャニスタ１１より機関吸気系に放出される蒸発ガスのパージ量がパージ限界（その時々に許容される限界値）に突入するか否かを、パージ非実行時に予測し、該パージ限界に突入する旨が予測された場合、前記パージ実行を禁止するようにした。上記構成によれば、非パージからパージへの移行後において、パージの開始と停止とが繰り返されるような従来装置の問題が回避できる。その結果、空燃比を変動なく制御し、ひいては排気エミッションの悪化を抑制することができる。
【００７６】
（ｂ）特に、本実施の形態の構成によれば、非パージの場合にのみ、パージ限界を予測し、その予測に応じてパージ実行の可否を判断するようにした。そのため、パージ実行中に強制的にパージを中断する等、パージを必要以上に停止させることはなく、パージ効率を悪化させるなどの問題を招くこともない。
【００７７】
（ｃ）さらに、本実施の形態では、蒸発ガスのパージ量が内燃機関１の運転に支障を来す領域をパージ限界として予測するようにした。このとき、内燃機関の運転に支障を来す領域とは、既述したようにパージの実行と停止とが繰り返され、結果として空燃比の乱れを生じさせるような以下の領域に設定した。
【００７８】
即ち、
・インジェクタ６が燃料噴射量を制御できる最小通電時間ＴＡＵMIN に基づいて設定された領域（図１０参照）、
・ＩＳＣ弁２０が負担空気量を制御できる最小制御値ＱｉMIN を基に設定された領域（図１１参照）、
・蒸発ガスのパージ量が気筒分配時に気筒間の差を生じうるレベルに基づいて設定された領域（ＦＰＧ＞ＦＫの領域）、
をパージ限界とした。
【００７９】
こうしたパージ限界の設定によれば、当該パージ限界である旨を好適且つ容易に予測することができ、上記パージ限界領域でのパージ開始を確実に行なわせないようにすることができる。
【００８０】
なお、本発明は、上記実施の形態の他に次の形態にて実現できる。
（１）上記実施の形態では、パージ流量を調整するためにパージデューティＶＳＶ１４を用いたが、本発明を実施する場合には上記構成に限定されるものではなく、他の構成に変更してもよい。要は、放出通路１５の経路中に配設され、蒸発燃料を含む空気のパージ率を変化させるものであれば任意に用いることができる。
【００８１】
（２）上記実施の形態では、酸素センサ９の電圧出力と所定判定レベルとを比較することにより空燃比フィードバック制御を実現したが、本発明を実施する場合には上記構成に限定されるものではなく、他の構成に変更してもよい。例えば排気ガス中の酸素濃度をリニアに検出する空燃比センサを用いたり、排気ガス中のＨＣ濃度を検出するＨＣセンサを用いたりして、空燃比フィードバック制御を実現してもよい。要は、内燃機関より排出される排気エミッションを抑制すべく空燃比をフィードバック制御するものであれば任意に具体化できる。
【００８２】
（３）上記実施の形態では、図８のパージ限界予測ルーチンにおいて、インジェクタ６の燃料噴射特性に応じたパージ限界の予測処理（ステップ５３１〜５３５）と、ＩＳＣ特性に応じたパージ限界の予測処理（ステップ５３６〜５３８）と、パージ補正係数ＦＰＧに応じたパージ限界の予測処理（ステップ５３９〜５４１）とを、各々独立して実施し、個々にパージ限界予測フラグを操作していたがこれを変更してもよい。例えば、上記３通りの処理のうち、いずれかでパージ限界予測フラグに「１」がセットされると、他の処理は実施しないようにしてもよい。
【００８３】
また、パージ補正係数ＦＰＧに応じたパージ限界の予測処理は、蒸発ガスのパージ量が気筒分配時に誤差を生じうるレベルを基に設定されるものであるが、その判定にはパージ補正係数ＦＰＧに代えてパージ流量Ｑｐやパージ率ＰＧＲを用いるようにしてもよい。
【００８４】
（４）上記実施の形態では、パージ実行許可フラグＸＰＲＧが「０」から「１」へ操作されるタイミングでパージ限界予測処理（図８の処理）を実行し、パージ限界突入の旨が予測された後には、パージ限界予測フラグが「０」にクリアされるまで同処理を繰り返し実行するようにしていたが、この構成を変更してもよい。例えばパージ限界突入の旨が予測された後に、所定時間間隔で図８の処理を実施してパージ限界突入の状態が回避されたかどうかを判別するようにしてもよい。要は、パージ非実行時にパージ限界を予測する構成であれば任意に変更できる。
【００８５】
（５）上記実施の形態では、図６に示すステップ５０８〜５１０のパージ率演算処理において、パージ限界に突入したか否かに応じてパージ率ＰＧＲを増減させていたが、これを変更してもよい。例えばＦＡＦＡＶ値の偏差ΔＦＡＦ（＝｜ＦＡＦＡＶ−１｜）をパラメータとし、ΔＦＡＦ値が所定値（例えば５％）以下の場合にＰＧＲ値を増加させると共に、ΔＦＡＦ値が所定値（例えば１０％）以上の場合にＰＧＲ値を減少させるように構成してもよい。要は、機関運転状態に応じてＰＧＲ値を設定する構成であれば任意に変更できる。
【００８６】
（７）上記実施の形態では、図２に示すルーチンのステップ１０５でＦＡＦ値の平均値ＦＡＦＡＶをなまし処理（平滑化処理）により算出したが、これに代えて、他の平均化処理によりＦＡＦＡＶ値を算出するようにしてもよい。
【図面の簡単な説明】
【図１】発明の実施の形態における内燃機関の空燃比制御装置の概要を示す構成図。
【図２】空燃比フィードバック制御ルーチンを示すフローチャート。
【図３】エバポ濃度演算ルーチンを示すフローチャート。
【図４】燃料噴射制御ルーチンを示すフローチャート。
【図５】アイドル回転数制御ルーチンを示すフローチャート。
【図６】デューティＶＳＶ制御ルーチンを示すフローチャート。
【図７】図６に続き、デューティＶＳＶ制御ルーチンを示すフローチャート。
【図８】パージ限界予測ルーチンを示すフローチャート。
【図９】パージ流量とデューティ比との関係を示すマップ。
【図１０】インジェクタによる燃料噴射特性を示すマップ。
【図１１】ＩＳＣ流量とデューティ比との関係を示すマップ。
【図１２】実施の形態における作用をより具体的に説明するためのタイムチャート。
【符号の説明】
１…内燃機関、２…吸気管、６…インジェクタ、１０…燃料タンク、１１…キャニスタ、１４…流量制御弁としてのパージデューティＶＳＶ、１５…放出通路、１６…パージ実行手段，パージ限界予測手段，パージ禁止手段を構成するＥＣＵ（電子制御装置）、２０…ＩＳＣ弁。
【図１】
0003651133.tif 000002
【図２】
0003651133.tif 000003
【図３】
0003651133.tif 000004
【図４】
0003651133.tif 000005
【図５】
0003651133.tif 000006
【図６】
0003651133.tif 000007
【図７】
0003651133.tif 000008
【図８】
0003651133.tif 000009
【図９】
0003651133.tif 000010
【図１０】
0003651133.tif 000011
【図１１】
0003651133.tif 000012
【図１２】
0003651133.tif 000013