(57)【特許請求の範囲】
【請求項１】
燃料タンクで発生したベーパの吸着及びその再脱離を許容する吸着剤と、エンジン吸気系への連通路であるパージラインとこの吸着剤との間に介在された空気層と、前記吸着剤を通過して前記パージラインに流入される空気の導入を許容する大気孔とを有するキャニスタを備え、前記燃料タンクで発生したベーパを前記キャニスタから前記パージラインを経由してエンジン吸気系にパージしつつ、燃焼室内で燃焼される混合気の空燃比を目標空燃比とするように燃料供給装置からエンジンへの燃料供給量をフィードバック補正することでエンジンの空燃比を制御する方法であって、
前記キャニスタ空気層でのベーパの蓄積状況を示す物理状態量（Mgair）、前記吸着剤でのベーパの蓄積状況を示す物理状態量（Mgcan）、前記燃料タンクでのベーパの発生状況を示す物理状態量（Fvptnk）、及び前記大気孔より導入される空気の流勢によって前記吸着剤から脱離して前記エンジン吸気系にパージされるベーパの流量である吸着剤脱離ベーパ流量（Fvpcan）に基づいたベーパ挙動についての物理モデルを用いて、前記パージラインを通じて前記エンジン吸気系にパージされるガスの総流量である総パージ流量（Fpgall）に応じた同ガス中のベーパの流量である総ベーパ流量（Fvpall）を推定し、
その推定した総ベーパ流量（Fvpall）に応じて前記フィードバック補正される燃料供給量を補正する
ことを特徴とするエンジンの空燃比制御方法。
【請求項２】
請求項１に記載のエンジンの空燃比制御方法において、
前記キャニスタ空気層に蓄積されたベーパの量である空気層蓄積ベーパ量（Mgair ）、及び前記吸着剤に蓄積されたベーパの量である吸着剤蓄積ベーパ量（Mgcan ）の算出、及びそれら算出される空気層蓄積ベーパ量（Mgair ）及び吸着剤蓄積ベーパ量（Mgcan ）に基づいて前記総パージ流量（Fpgall）に応じた前記総ベーパ流量（Fvpall）を推定する  ことを特徴とするエンジンの空燃比制御方法。
【請求項３】
請求項２に記載のエンジンの空燃比制御方法において、
前記空気層蓄積ベーパ量（Mgair ）と前記総パージ流量（Fpgall）とに応じて、前記キャニスタ空気層から前記パージラインに直接的に吸引されて前記エンジン吸気系にパージされるベーパの流量である空気層ベーパ流量（Fvpair）を算出するとともに、  前記吸着剤蓄積ベーパ量（Mgcan ）と前記総パージ流量（Fpgall）とに応じて、前記吸着剤脱離ベーパ流量（Fvpcan）を算出し、
それら算出される空気層ベーパ流量（Fvpair）と吸着剤脱離ベーパ流量（Fvpcan）との和として前記総ベーパ流量（Fvpall）を推定する  ことを特徴とするエンジンの空燃比制御方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、エンジンの空燃比制御方法、特に燃料タンクで発生したベーパ（燃料蒸気）をエンジン吸気系にパージ（放出）して処理するベーパパージシステムを備えるエンジンに適用されて好適な空燃比制御方法に関するものである。
【０００２】
【従来の技術】
一般に、揮発性液体燃料タンクを備える車両には、上記のようなベーパパージシステムが採用されている。代表的なチャコール・キャニスタ方式のパージシステムでは、燃料タンクで発生したベーパをキャニスタに一旦捕集するようにしている。キャニスタには、活性炭等の吸着剤が内蔵されており、ベーパを吸着させて一旦捕集した後、大気圧よりも低い圧力下におかれることでその吸着剤に蓄積させたベーパを再脱離させることが可能な構成となっている。こうしてキャニスタに捕集されたベーパは適宜、そのキャニスタからパージラインを通じてエンジン吸気系へとパージされ、エンジンに吸入された空気中に混入されている。そしてインジェクタから噴射された燃料とともに、エンジンの燃焼室内で燃焼させることで、上記燃料タンクで発生したベーパを処理するようにしている。
【０００３】
一方、エンジンの燃焼室に供給される可燃混合気の空燃比、すなわち吸入空気に占める噴射燃料の量（燃料供給装置からの燃料供給量）の比率を制御するエンジンの空燃比制御装置が知られている。こうした制御装置では、センサによって検出される実際の空燃比が目標空燃比に合致するように、インジェクタから噴射供給される燃料供給量をフィードバック補正している。
【０００４】
ところが、上記パージシステムを備えるエンジンでは、上記ベーパを含んだパージガスが、燃焼室に供給される本来の混合気に加えられる。したがって、こうしたパージシステムを備えるエンジンにおいて、空燃比制御のような燃焼室内で燃焼させる燃料供給量の厳密なコントロールの要求される制御を適合させるには、そうした制御にパージガスの影響分を見込んで燃料供給量を調整することが必要となる。
【０００５】
そこで従来は、次のようにしてパージガスの影響分を見込んだ空燃比制御を実現している。すなわち、上記空燃比のフィードバックに係る燃料供給量の補正値（空燃比フィードバック補正値）について、そのパージガスのパージガス流量変化時の値の変化から、そのパージガス中の燃料成分の濃度（ベーパ濃度）を推定する。以降、その推定された燃料成分の濃度とパージガスの流量とから、パージを通じてエンジンに供給されるベーパ流量を求め、その分、インジェクタからの燃料の噴射量を減量補正する。そしてエンジンの運転条件が所定条件を満たす都度、上記ベーパ濃度を同様に求め、その推定値を修正していくことで、制御の適合を図るようにしている。
【０００６】
こうした態様での空燃比制御は、ベーパ濃度がパージ流量に拘わらず一定であって、且つパージガス中のベーパ成分の濃度の変化が充分に緩やかである場合には、充分に有効に機能する。すなわち、エンジン吸気通路へのパージ流量とその中に含まれるベーパの流量とが、線形関係に有ることを前提条件として、空燃比制御の適合がなされている。
【０００７】
ところで、給油時等のベーパの大量発生時には、一時に過剰なベーパが吸着剤に吸着されてその劣化を招くことがある。このため従来は、例えば特開平９−１８４４４４号公報にみられるように、キャニスタ内に吸着剤の充填されていない空間を設け、その空間内の空気の層（キャニスタ空気層）を緩衝帯とすることで、吸着剤の劣化を抑制する構成のパージシステムも提案されている。
【０００８】
こうしたパージシステムでは、状況によっては、燃料タンクで発生したベーパの一部が、上記キャニスタ空気層を通過して、吸着剤に捕集されることなく直接にエンジンの吸気通路にパージされることが生じ得る。
【０００９】
そのため、上記公報記載のエンジンの空燃比制御装置では、
（ａ）吸着剤に吸着されることなく燃料タンクから吸気通路に直接パージされる態様、及び、
（ｂ）吸着剤に一旦吸着された後にその吸着剤から再脱離されて吸気通路にパージされる態様、
の２つの態様で、エンジンにベーパが流入すると想定し、パージの影響分を見込んだ制御の適合を図っている。以下では、前者の態様（ａ）でのパージを「タンク流入パージ」と云い、後者の態様（ｂ）でのパージを「吸着剤脱離パージ」と云うこととする。これら「タンク流入パージ」と「吸着剤脱離パージ」とでは、パージ中のベーパ挙動は自ずと異なったものとなる。そしてその結果、上記制御の前提条件の一つであるパージ流量とベーパ流量との線形関係は、常には成立しなくなる。例えば、吸気通路へのパージ流量が同一であれ、燃料タンクからのベーパ流入が有るときと無いときとでは、パージ中のベーパ挙動は全く異なったものとなる。
【００１０】
そこで上記公報記載のエンジンの空燃比制御装置では、上記「タンク流入パージ」分の吸気通路へのベーパ流量Fvptnkと、上記「吸着剤脱離パージ」分の吸気通路へのベーパ流量Fvpcanとをそれぞれ別個に求めている。そしてそれら２つのベーパ流量を各別の算出態様で算出し、エンジン吸気系にパージされるベーパの総流量（総ベーパ流量）の推定値Fvpallを求めるようにしている。
【００１１】
具体的には、次の数式によって上記各ベーパ流量Fvptnk、Fvpcanを各算出し、それらの総和（Fvptnk＋Fvpcan）を上記総ベーパ流量Fvpallとして推定し、その推定値に基づいてインジェクタからの燃料噴射量を補正するようにしている。
<<参考式>>
・  Fvptnk  ←  rvptnk／（Ｑ・Fpgall）
・  Fvpcan  ←  rvpcan・Fpgall
・  Fvpall  ←  Fvptnk＋Fvpcan
ここで、「Ｑ」：吸入空気量、「rvptnk」：タンク流入パージのベーパ濃度（パージガス中のベーパの含有率）、「rvpcan」：吸着剤脱離パージのベーパ濃度、を各示している。
【００１２】
すなわち、上記公報記載の空燃比制御装置では、タンク流入パージのベーパ流量Fvptnkと吸着剤脱離パージのベーパ流量Fvpcanとをそれぞれ別途に算出し、それらのベーパ流量の総和として総ベーパ流量Fvpallを算出する構成となっている。
【００１３】
【発明が解決しようとする課題】
以上の態様でベーパ流量を推定すれば、燃料タンクからキャニスタへのベーパ流入の有無によるベーパ濃度条件の差異に応じたベーパ流量の推定が可能となって、ある程度は空燃比制御等の精度向上を期することができる。
【００１４】
しかしながら、実際のパージシステムでのベーパ挙動は、上記制御装置でのベーパ流量の推定ロジックの設定に際して想定されているものよりも遙かに複雑であることが、発明者らの試験等の結果によって確かめられている。このため、上記公報記載の態様でのベーパ流量の算出ロジックを採用した場合であれ、その算出精度を充分に高めることはできず、空燃比制御等に対するパージの影響緩和にも自ずと限界があった。
【００１５】
本発明は、こうした実情に鑑みてなされたものであって、その目的は、燃料タンクで発生したベーパをパージして処理するベーパパージシステムを備えるエンジンにあって、そのパージによるエンジンへのベーパ流量をより正確に推定して、空燃比制御等に与えるパージの影響を好適に抑制することのできるエンジンの空燃比制御方法を提供することにある。
【００１６】
【課題を解決するための手段】
以下、上記目的を達成するための構成及びその作用効果について記載する。
［請求項１］
エンジンの空燃比制御方法としての請求項１に記載の発明は、燃料タンクで発生したベーパの吸着及びその再脱離を許容する吸着剤と、エンジン吸気系への連通路であるパージラインとこの吸着剤との間に介在された空気層と、前記吸着剤を通過して前記パージラインに流入される空気の導入を許容する大気孔とを有するキャニスタを備え、前記燃料タンクで発生したベーパを前記キャニスタから前記パージラインを経由してエンジン吸気系にパージしつつ、燃焼室内で燃焼される混合気の空燃比を目標空燃比とするように燃料供給装置からエンジンへの燃料供給量をフィードバック補正することでエンジンの空燃比を制御する方法であって、前記キャニスタ空気層でのベーパの蓄積状況を示す物理状態量（Mgair）、前記吸着剤でのベーパの蓄積状況を示す物理状態量（Mgcan）、前記燃料タンクでのベーパの発生状況を示す物理状態量（Fvptnk）、及び前記大気孔より導入される空気の流勢によって前記吸着剤から脱離して前記エンジン吸気系にパージされるベーパの流量である吸着剤脱離ベーパ流量（Fvpcan）に基づいたベーパ挙動についての物理モデルを用いて、前記パージラインを通じて前記エンジン吸気系にパージされるガスの総流量である総パージ流量（Fpgall）に応じた同ガス中のベーパの流量である総ベーパ流量（Fvpall）を推定し、その推定した総ベーパ流量（Fvpall）に応じて前記フィードバック補正される燃料供給量を更に補正するようにしたものである。
【００１７】
ベーパパージシステムでのベーパ挙動は、キャニスタにおけるキャニスタ空気層及び吸着剤でのベーパの各蓄積状況、燃料タンクでのベーパの発生状況の３つの物理状態量（Mgair、Mgcan、Fvptnk）、及び吸着剤脱離ベーパ流量（Fvpcan）に基づいた物理モデルによって説明できる（図１３、図４６等参照）。また、パージシステムでのベーパ挙動は、パージ状況や燃料タンクでのベーパの発生状況等に応じて、刻一刻と変動していく。しかしながら、上記物理モデルによれば、上記各物理状態量（Mgair、Mgcan、Fvptnk）、及び吸着剤脱離ベーパ流量（Fvpcan）に基づくことで、かかるベーパ挙動の変化に合わせて、パージラインからエンジン吸気系にパージされるベーパの流量（総ベーパ流量Fvpall）を正確に推定可能である。したがって、上記構成によれば、パージシステムでのベーパ挙動の変化に拘わらず、パージラインからエンジン吸気系にパージされるベーパの流量を正確に予測し、パージ実施中の空燃比を高精度に制御することができるようになる。
【００１８】
なお、上記各物理状態量は、パージ実施中の空燃比フィードバックのずれや燃料タンクの内圧などから、直接的、或いは間接的に検出することができる。ただし、上記物理モデルを用いれば、パージ状況や各物理状態量自体の現状値に基づいて、それら各物理状態量の推移を推定することができる。
【００１９】
そこで上記構成において、上記物理モデルを用い、パージ状況及び各物理状態量の現状値に基づいて、各物理状態量の値を周期的に更新するようにしても良い。この場合、オープンループの演算処理のみによる、すなわちフィードフォワードによる総ベーパ流量Fvpallの推定も可能となる。そしてこれにより、センサ等の検出結果に基づくフィードバックに依らずとも、上記ベーパ挙動の変化に対応したパージ実施中の空燃比の高精度な制御が可能となる。
【００２０】
一方、上記各物理状態量の値に誤差が有れば、空燃比制御等に好ましからぬ影響を与えるおそれがある。特にそれら各物理状態量は、パージ状況等に応じて推移していくため、その変化に併せてその値を修正しなければ、たとえある時点において適正な値となっていても、やがては誤差が生じてしまう。また、上記のように物理モデルに基づいて各物理状態量の値を周期的に更新する場合であれ、更新毎の各物理状態量の値のわずかなエラーが更新を重ねる毎に積み重なり、やがては大きな誤差を生むおそれがある。更には、上記各物理状態量を直接的、或いは間接的に測定する場合にせよ、センサやシステムの劣化等に伴う経時変化や個体差などのため、各物理状態量の値を誤る可能性は十分に考えられる。
【００２１】
そこで上記構成において、パージ実施に伴う空燃比フィードバック補正値の推移を監視し、その値のずれに応じて各物理状態量の値を修正するようにしても良い。この場合には、更新時の誤差の積み重ね等によって各物理状態量の値にずれが生じようとも、そのずれを補正可能となり、それら各物理状態量を正確な値に保持できる。
【００２２】
［請求項２］
請求項２に記載の発明は、請求項１に記載のエンジンの空燃比制御方法において、前記キャニスタ空気層に蓄積されたベーパの量である空気層蓄積ベーパ量（Mgair ）、及び前記吸着剤に蓄積されたベーパの量である吸着剤蓄積ベーパ量（Mgcan ）の算出、及びそれら算出される空気層蓄積ベーパ量（Mgair ）及び吸着剤蓄積ベーパ量（Mgcan ）に基づいて前記総パージ流量（Fpgall）に応じた前記総ベーパ流量（Fvpall）を推定するようにしたものである。
【００２３】
上記物理モデル（図１３等参照）によれば、エンジン吸気系にパージされるベーパの流量、すなわち総ベーパ流量Fvpallは、大きくは空気層蓄積ベーパ量Mgair 、吸着剤蓄積ベーパ量Mgcan 、及び総パージ流量Fpgallの関数として求めることができる。このため上記構成によれば、それら空気層蓄積ベーパ量Mgair 、吸着剤蓄積ベーパ量Mgcan 、及び総パージ流量Fpgallさえ的確に把握すれば、パージシステムでのベーパ挙動の変化に拘わらず、上記総ベーパ流量Fvpallを正確に推定可能となる。
<<参考式>>
・  Fvpall  ←  Fnc.｛Fpgall、Mgair 、Mgcan ｝
なお本明細書の説明では、「Fnc.｛ｘ１、ｘ２、・・・｝」は、パラメータｘ１、ｘ２・・・による所定の関数を表すものとする。
【００２４】
［請求項３］
請求項３に記載の発明は、請求項２に記載のエンジンの空燃比制御方法において、前記空気層蓄積ベーパ量（Mgair ）と前記総パージ流量（Fpgall）とに応じて、前記キャニスタ空気層から前記パージラインに直接的に吸引されて前記エンジン吸気系にパージされるベーパの流量である空気層ベーパ流量（Fvpair）を算出するとともに、前記吸着剤蓄積ベーパ量（Mgcan ）と前記総パージ流量（Fpgall）とに応じて、前記吸着剤脱離ベーパ流量（Fvpcan）を算出し、それら算出される空気層ベーパ流量（Fvpair）と吸着剤脱離ベーパ流量（Fvpcan）との和として前記総ベーパ流量（Fvpall）を推定するようにしたものである。
【００２５】
上記物理モデル（図１３等参照）によれば、パージ実施中にパージラインからエンジン吸気系へのパージの成分は、次の２つの成分に大別できる。すなわち、キャニスタ空気層に蓄積された状態のベーパがパージラインに直接的に吸引されてパージされる態様（空気層パージ）と、吸着剤に蓄積された状態から大気孔より導入された空気の流勢によって脱離してその空気と共にパージされる態様（吸着剤脱離パージ）とである（［１−３−１］及び［１−３−３］、図１３等参照）。パージ実施中のそれら両パージ成分それぞれでのベーパの挙動は、自ずと異なったものとなる。
【００２６】
一方、上記物理モデルによれば、空気層パージ中のベーパの流量（空気層ベーパ流量Fvpair）は、上記空気層蓄積ベーパ量Mgair と総パージ流量Fpgallとの関数として求めることができる。また、吸着剤脱離パージ中のベーパの流量（吸着剤脱離ベーパ流量Fvpcan）は、上記吸着剤蓄積ベーパ量Mgcan と総パージ流量Fpgallとの関数として求めることができる（［１−３−３］、図１７等参照）。またそのときの総ベーパ量は、それら空気層ベーパ流量Fvpairと吸着剤脱離ベーパ流量Fvpcanとの和として求められる。
【００２７】
このため上記構成では、そうしたベーパ挙動の異なる２つのパージ態様でのベーパ流量を別途求めた後、総ベーパ流量Fvpallを算出するようにしている。したがって、ベーパ挙動の異なる、すなわち総パージ流量Fpgallに対するベーパ流量の相関関係の異なる両成分でのベーパの各流量を的確に求めて総ベーパ流量Fvpallをより正確に推定可能となる。
<<参考式>>
・  Fvpair  ←  Fnc.｛Fpgall、Mgair ｝
・  Fvpcan  ←  Fnc.｛Fpgall、Mgcan ｝
・  Fvpall  ←  Fvpair＋Fvpcan。
【００４９】
また上記の如く、パージ実施中のタンク流入ベーパ流量Fvptnkが多い場合には、更にそのタンク流入ベーパ流量Fvptnkも併せ加味して各流量を算出するようにすれば、吸着剤内空気流量Fpgcanの値をより的確に求めることができるようにもなる。
<<参考式>>
・  Fvpcan  ←  Fnc.｛Fpgcan、Mgcan ｝。
【０１４３】
【発明の実施の形態】
以下、本発明の実施の形態を、図面を参照して説明する。
発明者らは、下記の如く構成されたベーパパージシステム（以下「パージシステム」という）でのベーパの挙動を試験等によって詳細に検証した。そしてその結果に基づいて、後述するパージシステムのベーパ挙動の物理モデルを提唱した（図１３等参照）。
【０１４４】
この物理モデルによれば、パージシステムでのベーパ挙動についての下記各特性が導き出される。
・パージシステムのベーパ挙動は、キャニスタのキャニスタ空気層及び吸着剤でのベーパの各蓄積状態、及び燃料タンクでのベーパの発生状態をそれぞれ示す３つの物理状態量の相関関係によって表される。
・その物理モデルに従えば、キャニスタからエンジン吸気系へとパージされるベーパの流量（総ベーパ流量）は、エンジン吸気系にパージされるガスの流量（総パージ流量）、及び上記各物理状態量の関数として表すことができる。
・また上記物理モデルに従えば、パージシステム内での上記各物理状態量の推移も、エンジン吸気系へのパージ状況、及びそれら各物理状態量の現状値から、詳細に把握可能となる。
【０１４５】
以下、こうした物理モデルの詳細、及びそれを適用したエンジンの空燃比制御方法について、詳細に説明する。
まず始めに、［１］として、上記パージシステムにおけるベーパ挙動の物理モデルの詳細を説明する。ここで説明の概要は次の通りである。
【０１４６】
［１−１］  パージシステムの基本構成
ここでは、かかる物理モデルの適用対象となるパージシステムの基本構成を、図１を参照して説明する。
【０１４７】
［１−２］  パージシステムのベーパ挙動の検証
ここでは、上記［１−１］で説明するパージシステムを用いて行われたベーパ挙動の検証試験の結果、及びその結果から導出されたベーパ挙動の特性について、図２〜図６を併せ参照して説明する。なお、［１−２−１］ではタンク発生ベーパの挙動について、［１−２−２］では空気層蓄積ベーパの挙動について、［１−２−３］では吸着剤蓄積ベーパの挙動について、それぞれ説明する。
【０１４８】
［１−３］  パージシステムでのベーパ挙動の物理モデル
ここでは、上記［１−２］の検証試験結果に基づき提起された物理モデルの詳細を、図７〜図１３を併せ参照して説明する。なお、［１−３−１］ではキャニスタ空気層内でのベーパ挙動の物理モデルについて、［１−３−２］ではキャニスタ内でのベーパ挙動の物理モデルについて、［１−３−３］ではパージ中のベーパ挙動の物理モデルについて、それぞれ説明する。また［１−３−４］では、それらを総括したパージシステム全体でのベーパ挙動の物理モデルの全体像を説明する。
【０１４９】
これに続き、［２］として上記物理モデルを適用したエンジンの空燃比制御装置の具体例を説明する。ここでの説明の概要は次の通りである。
［２−１］  空燃比制御装置の全体構成
ここでは、上記物理モデルに基づく制御の適用される空燃比制御装置について、その全体構成を、図１４を参照して説明する。
【０１５０】
［２−２］  パージ制御の概要
ここでは、上記物理モデルに基づくパージに係る制御の全体像を、図１５を併せ参照して概略的に説明する。
【０１５１】
またこれに引き続き、上記［２−２］でその概要の説明される制御の詳細を説明する。その説明内容は次の通りである。
［２−３］  物理モデルに基づく各物理状態量の定期更新処理
ここでは、上記物理モデルに基づいて行われる各物理状態量の定期更新処理について、同図１５を参照してその詳細を説明する。
【０１５２】
［２−４］  パージ補正値の算出処理
ここでは、上記物理モデルに基づいて行われる空燃比制御にかかるパージ補正値の算出にかかる処理について、図１６〜図２２を併せ参照して、詳細に説明する。
【０１５３】
［２−５］  燃料噴射量の算出処理
ここでは、上記パージ補正量に応じた燃料噴射量の算出にかかる処理について説明する。なおここでは、空燃比フィードバック制御の概要についても、図２３を参照して併せ説明する。
【０１５４】
［２−６］  各物理状態量の初期化処理
ここでは、各物理状態量の初期化にかかる処理の詳細について、図２４〜図２６を併せ参照して説明する。
【０１５５】
［２−７］  各物理状態量の修正処理
ここでは、各物理状態量の修正にかかる処理の詳細について、図２７〜図３４を併せ参照して説明する。
【０１５６】
［２−８］  ＶＳＶ開度の算出処理
ここでは、ＶＳＶ開度の算出に係る処理の詳細について、図３５〜図３７を併せ参照して説明する。
【０１５７】
更に引き続いて、［２−９］では、本空燃比制御装置に加え得る更なる改良点を説明する。そうした各改良点の説明の概要は次の通りである。
［２−９−１］では、低開度でのＶＳＶ制御に係る改良点について図３８〜図４０を併せ参照して説明する。
【０１５８】
［２−９−２］では、空燃比フィードバック補正値の中心値の算出処理に係る改良点について、図４１を参照して説明する。
［２−９−３］ではパージ流量の濃度補正処理に係る改良点について、図４２を参照して説明する。
【０１５９】
［２−９−４］では、各物理量の更新エラーの低減処理に係る改良点について、図４３〜図４５を参照して説明する。
［２−９−５］では、タンク発生ベーパの直接流入への対策処理について、図４６に基づき説明する。
【０１６０】
以上が、本明細書及び図面における「発明の実施の形態」の説明内容の概要である。
なお本明細書及び図面の説明では、「ベーパ」とは燃料タンクで発生した燃料蒸気を、「パージガス」とはそのベーパと空気との混合気を、各々示している。そして「ベーパ量」は、パージガス中等に含まれるベーパ成分の質量を示し、また「パージ流量」及び「ベーパ流量」はそれぞれパージガス、ベーパについて、それらが単位時間当たりに移動した質量を示している。
【０１６１】
［１］  パージシステムのベーパ挙動の物理モデル
［１−１］  パージシステムの基本構成
まずは、上記物理モデルの適用対象となるパージシステムの基本構成を、図１を参照して説明する。
【０１６２】
同図１に示すように、このパージシステムは、ベーパを捕集するキャニスタ３、そのキャニスタ３と燃料タンク１とを連結するエバポライン２、及びエンジン５の吸気系を構成する吸気通路６とキャニスタ３とを連結するパージライン４とを備えて構成されている。パージライン４の途中には、パージ調整バルブ（ＶＳＶ）７が設けられており、吸気通路６に導入されるパージガスの流量がその開度制御を通じて調整可能となっている。
【０１６３】
キャニスタ３の下部には、外部の空気（大気）を導入するための大気孔８が形成され、そのキャニスタ３の内部には、上記エバポライン２を通じて導入されたベーパを吸着して捕集するための吸着剤３ａが内蔵されている。ただしこのキャニスタ３の内部には、吸着剤３ａの上部に空間が残されており、その空間を満たす空気の層（キャニスタ空気層）３ｂが形成されている。このキャニスタ３にあって上記エバポライン２及びパージライン４は共に、このキャニスタ空気層３ｂに開口されている。
【０１６４】
さて、このパージシステムでは、燃料タンク１で発生したベーパは、エバポライン２を通じてキャニスタ３に送られ、キャニスタ空気層３ｂのパージガス中に一旦混入された後、徐々に吸着剤３ａに吸着されるようになっている。
【０１６５】
ここでエンジンの運転時にＶＳＶ７を開弁すると、吸気通路６内の圧力（吸気通路内圧ＰＭ）との差圧によって、キャニスタ３の内部からパージライン４にガスが引き込まれ、その引き込まれたガスが吸気通路６にパージされる。このとき、そうしたパージガスの流量が十分に多ければ、大気孔８から外部の空気が導入され、吸着剤３ａ内を通過してパージライン４に流入する。そしてこうした空気の流れによって、吸着剤３ａに吸着されたベーパが再脱離され、パージライン４を通じて吸気通路６にパージされる。以上が本パージシステムでのベーパ挙動の概要である。
【０１６６】
［１−２］  パージシステムのベーパ挙動の検証
発明者らは、以上のように構成されたパージシステムにおいて、下記の各種試験を行い、そのシステムでのベーパ挙動を更に詳細に検証した。以下に、その結果を記載する。
【０１６７】
［１−２−１］  タンク発生ベーパの挙動
ここでは、燃料タンク１でのベーパの発生状況、すなわち燃料タンク１からキャニスタ３に送られるベーパの流量（タンク発生ベーパ流量Fvptnk）が、吸気通路６へのパージに与える影響について説明する。発明者らは、そうした影響を検証すべく、下記（Ｉ）（II）の試験を行った。
【０１６８】
（Ｉ）  定常状態でのタンク発生ベーパ流量Fvptnkの影響の検証試験
まずは、定常状態、すなわちタンク発生ベーパ流量Fvptnkを一定に保持し、且つ吸気通路６の内圧（吸気通路内圧ＰＭ）及びＶＳＶ７の開度（ＶＳＶ開度）を一定に、すなわち吸気通路６へのパージ流量を一定に保持した状態で、吸気通路６にパージされるベーパの流量を測定した。なお本試験では、吸着剤３ａから脱離したベーパの影響を排除すべく、キャニスタ３内に蓄積されたベーパをほぼ完全に空とした状態で上記測定を開始した。そして、タンク発生ベーパ流量FvptnkやＶＳＶ開度の異なる複数の条件で測定を行った。
【０１６９】
その結果は、図２に示す通りである。すなわち、
（イ）燃料タンク１でのベーパの発生量、すなわち上記タンク発生ベーパ流量Fvptnkが一定であれば、ＶＳＶ開度が十分に小さな領域では、すなわち総パージ流量Fpgallが十分に少ない領域では、その開度の増大に応じて吸気通路６へのベーパ流量も増加する。ただし、総パージ流量Fpgallがある程度を上回ると、吸気通路６へのベーパ流量は一定値に飽和する。
（ロ）飽和して一定値となったときの吸気通路６へのベーパ流量は、タンク発生ベーパ流量Fvptnkによって決まり、その吸気通路６へのベーパ流量は、タンク発生ベーパ流量Fvptnkの増大に応じて多くなる。
【０１７０】
（II）  過渡状態でのタンク発生ベーパ流量Fvptnkの影響の検証試験
続いて、タンク発生ベーパ流量Fvptnkを一定に保持した状態で、パージカット（ＶＳＶ８を全閉）を所定時間行ってからパージを開始した後のベーパ流量の推移を測定した。なお、本試験においても、吸着剤３ａから脱離したベーパの影響を排除すべく、キャニスタ３内に蓄積されたベーパをほぼ完全に空とした状態としてから、パージカットを開始している。また、パージ開始後の吸気通路内圧ＰＭ、及びＶＳＶ開度を一定として、パージ流量を一定に保持した状態でベーパ流量の推移を測定した。そして、タンク発生ベーパ流量Fvptnkやパージカット時間の異なる複数の条件で測定を行った。
【０１７１】
その結果は、図３に示す通りである。すなわち、
（ハ）パージ開始直後のベーパ流量は、パージカットの時間の長さ、及びタンク発生ベーパ流量Fvptnkの多さに応じて増加する。ただし、ベーパ流量は、パージカット時間、或いは燃料タンク１でのベーパ発生量（タンク発生ベーパ流量Fvptnk）に対して単純な比例関係とはならない。
（ニ）パージ開始後、ベーパ流量は時間と共に徐々に減少し、やがては一定値に飽和する。そして、タンク発生ベーパ流量Fvptnkが一定であれば、飽和して一定となった後のベーパ流量は同じ値となる。なお、飽和して一定となった後のベーパ流量は、タンク発生ベーパ流量Fvptnkが同一であれば、上記定常状態において、やはり飽和して一定となった後のベーパ流量と同じである。
【０１７２】
ちなみに、同図には図示しないが、ＶＳＶ開度や吸気通路内圧ＰＭの異なる複数の条件においても、同様の測定を行った。その結果によれば、それらを変化させようとも、すなわち上記総パージ流量Fpgall（吸気通路６に流入されるパージガスの流量）を変化させても、上記試験結果（ハ）、（ニ）に示される傾向は変化しないことが確認されている。
【０１７３】
なおそのとき、総パージ流量Fpgallがある程度よりも少ない領域では、上記傾向は変化しないものの、パージ開始時のベーパ流量の値やその後のベーパ流量の減少率は、総パージ流量Fpgallの増大に応じて大きくなる。ただし、総パージ流量Fpgallがある程度よりも多くなると、総パージ流量Fpgallをそれ以上に幾ら増大しようとも、それらパージ開始時のベーパ流量の値やその後のベーパ流量の減少率はほとんど変化しなくなることが確認されている。
【０１７４】
［１−２−２］  空気層蓄積ベーパの挙動
ここでは、キャニスタ空気層３ｂに混入した状態で蓄積されているベーパ（空気層蓄積ベーパ）が吸気通路６に流入するベーパ流量に与える影響について説明する。発明者らは、そうした空気層蓄積ベーパの影響を、上記過渡状態でのタンク発生ベーパ流量Fvptnkの影響の検証試験の結果に基づいて検証した。
【０１７５】
発明者らは、上記パージカット期間に燃料タンク１からキャニスタ３に導入されたベーパの総量、すなわちパージカット時間中のタンク発生ベーパ流量Fvptnkの積算値と、パージカット復帰直後、すなわちパージカットを中止してパージを開始した直後における吸気通路６へのベーパ流量との関係を調査した。なお、ここでの調査にはパージ開始直後のベーパ流量として、総パージ流量Fpgallが十分に大きく、その変化に依存しない領域での測定値を、すなわちその条件下での上記ベーパ流量の最大値を、その調査に用いている。
【０１７６】
このとき、パージカット期間にキャニスタ３に導入されたベーパの全てがキャニスタ空気層３ｂに蓄積されるものと仮定すれば、上記空気層蓄積ベーパ量Mgair と吸気通路６へのベーパ流量の最大値との間には、特別な因果関係はみられない。しかしながら、実際には、燃料タンク１からキャニスタ空気層３ｂに導入されたベーパは、吸着剤３ａに徐々に吸着されていく。
【０１７７】
ここで後述するキャニスタ３内でのベーパ挙動の物理モデル（本章２−５参照）に従い、キャニスタ空気層３ｂから吸着剤３ａへのベーパの吸着速度が空気層蓄積ベーパ量Mgair に比例すると仮定する。この仮定によれば、パージ開始時の空気層蓄積ベーパ量Mgair は、パージカット期間のタンク発生ベーパ流量Fvptnkの積算値から、上記関係に基づいて推定される同期間のキャニスタ空気層３ｂから吸着剤３ａへのベーパの吸着分の総量を差し引いたものとなる。この仮定に基づき推定される空気層蓄積ベーパ量Mgair は、パージ開始直後のベーパ流量と非常に高い相関性を持つことが上記試験結果の検証から確認された。
【０１７８】
そして、そうした検証結果から、空気層蓄積ベーパの挙動に関する下記の２つの傾向が確認された。
（ホ）吸着剤３ａから脱離してパージされるベーパ流量（吸着剤脱離ベーパ流量Fvpcan）が無いときには、吸気通路６への最大ベーパ流量は、空気層蓄積ベーパ量Mgair によってほぼ一義的に求められる。
（ヘ）またそのときの最大ベーパ流量は、空気層蓄積ベーパ量Mgair の増大に応じて増加するものの、やがては飽和する。
【０１７９】
［１−２−３］  吸着剤蓄積ベーパの挙動
ここでは、吸着剤３ａに吸着された状態で蓄積されたベーパの量（吸着剤蓄積ベーパ量Mgcan ）の挙動について説明する。発明者らは、そうした挙動を検証すべく、下記（Ｉ）（II）の試験を行った。
【０１８０】
（Ｉ）  パージ実施中の吸着剤蓄積ベーパの影響の検証試験
まずは、吸着剤３ａに所定量のベーパを吸着させた状態でパージを開始し、その後の吸気通路６へのベーパ流量の推移を測定した。またこれと同時に、吸着剤３ａに残存するベーパの吸着量（吸着剤蓄積ベーパ量Mgcan ）の測定も並行して行った。こうした測定を、吸着剤蓄積ベーパ量Mgcan の初期条件を変えて複数回行った。ちなみに本試験では、タンク発生ベーパの影響を排除すべく、燃料タンク１からのベーパの流入を遮断した状態で上記測定を行っている。
【０１８１】
その結果は、図４に示す通りである。同図４には、吸着剤蓄積ベーパ量Mgcan の初期条件や吸気通路６へのパージ流量を変化させて行った複数回の測定について、その結果から求められたベーパ濃度（吸気通路６にパージされるガス中のベーパ成分の濃度）と吸着量（吸着剤蓄積ベーパ量Mgcan ）との関係が示されている。この図４に示されるように、吸着剤蓄積ベーパ量Mgcan の初期条件や吸気通路６へのパージ流量を変化させても、上記関係は一定である。
【０１８２】
この結果より、パージ実施中の吸着剤蓄積ベーパの挙動について、以下の傾向が確認された。
（ト）キャニスタ空気層３ｂからパージライン４へのベーパ流入が無いときには、吸着剤蓄積ベーパ量Mgcan が一定であれば、吸気通路６へのパージ流量に拘わらず、そのベーパ濃度は一定である。よって吸着剤蓄積ベーパ量Mgcan が一定であれば、パージ実施中に、大気孔８から導入された空気の流勢によって吸着剤３ａから脱離してパージされるベーパ、すなわち吸着剤脱離ベーパの流量（吸着剤脱離ベーパ流量Fvpcan）は、図５に示すように、吸気通路６へのパージ流量に比例する。
（チ）一方、吸着剤蓄積ベーパ量Mgcan は、吸着剤３ａに蓄積されたベーパが脱離してパージされるにつれて、その分ずつ減少していくのは自明である。よって、吸着剤蓄積ベーパ量Mgcan は、吸着剤３ａから脱離してパージされたベーパ流量の積算値から相対的に求めることができる。
【０１８３】
（II）  パージカット中の吸着剤蓄積ベーパの影響の検証試験
一方、吸着剤蓄積ベーパは、大気孔８からの空気の流勢に依らずとも、自然にその一部が、キャニスタ空気層３ｂへと少しずつ脱離していくものと考えられる。そこで発明者らは、上記パージ実施中の吸着剤蓄積ベーパの影響の検証にかかる測定の途中、所定時間ずつパージカットを実施し、そのパージカット前後のベーパ流量の推移から、そうした吸着剤蓄積ベーパの挙動を検証した。
【０１８４】
図６は、そうした検証の結果から求められた、パージカット中の吸着剤３ａからのベーパの脱離速度と吸着剤蓄積ベーパ量Mgcan との関係を示している。ここでのベーパの脱離速度は、パージカット前後のベーパ流量の差とそのパージカットの実施時間とから求められている。
【０１８５】
これらの関係からは、同図６に示されるように、吸着剤３ａから自然脱離するベーパの挙動についての以下の傾向が確認された。
（リ）パージカット中に吸着剤３ａからキャニスタ空気層３ｂへと自然脱離するベーパの流量、すなわち自然脱離速度Fvpctaと上記吸着剤蓄積ベーパ量Mgcan とは、ほぼ線形関係にある。
（ヌ）ただし、そうした自然に脱離するベーパの流量は、パージ実施中に大気孔８から導入された空気の流勢によって吸着剤３ａから脱離してパージされるベーパの流量（吸着剤脱離ベーパ流量Fvpcan）に比して、非常に少ない。
【０１８６】
［１−３］  パージシステムでのベーパ挙動の物理モデル
続いて以上のベーパ挙動の検証結果に鑑みて、発明者らにより提唱された物理モデルの詳細を説明する。
【０１８７】
［１−３−１］  キャニスタ空気層内でのベーパ挙動の物理モデル
まず、パージ実施時におけるキャニスタ空気層３ｂに蓄積されたベーパの挙動についての物理モデルについて、図７〜図１０を併せ参照して説明する。この物理モデルによれば、パージ実施中のキャニスタ空気層３ｂに蓄積されたベーパの挙動は、以下の通りとなる。
（ａ）パージ実施中に、キャニスタ空気層３ｂのベーパ分を含んだパージガスは、大気孔から導入されて吸着剤内を通過する空気よりも、優先的にパージライン４に吸入されてパージされる。すなわち上記空気層パージは、吸着剤脱離パージよりも優先的にエンジンの吸気通路６にパージされる。
（ｂ）パージ実施中に、キャニスタ空気層３ｂから吸気通路６にパージされるガスの流量（空気層パージ流量Fpgair）の最大値である最大空気層パージ流量Fpgairmxは、キャニスタ空気層３ｂに蓄積されたベーパの量、すなわち空気層蓄積ベーパ量Mgair の値によって一義的に導き出される。また同様に、キャニスタ空気層３ｂから吸気通路６にパージされるガス中のベーパの流量（空気層ベーパ流量Fvpair）の最大値である最大空気層ベーパ流量Fvpairmxも、空気層蓄積ベーパ量Mgair の値によって一義的に導き出される。
【０１８８】
以下に、上記仮定事項（ａ）及び（ｂ）について、その理論的根拠、及びその詳細を説明する。
上記［１−２−１］及び［１−２−２］で説明したように、吸着剤蓄積ベーパ量Mgcan が「０」であるときのパージ実施中に、総パージ流量Fpgall、及びタンク発生ベーパ流量Fvptnkが所定の限界値を超えると、総ベーパ流量Fvpallは一定値となることが確認されている（図２、図３等参照）。こうした測定結果に鑑みて発明者らは、パージ実施中のキャニスタ３内でのパージガスの挙動について、図７に示されるような物理モデルを想定した。
【０１８９】
パージ実施中には、キャニスタ空気層３ｂに蓄積されたベーパ分を含んだガス（パージガス）がパージライン４に吸入されるのと同時に、大気孔８を通じて外部から導入された空気（大気）もパージライン４に吸入される。よってキャニスタ空気層３ｂのパージガスは、大気孔８から導入された空気の干渉を受けつつ、パージライン４に吸入されることとなる。そこで上記物理モデルでは、同図７に示されるように、「キャニスタ空気層３ｂのパージガスは、パージ実施中に、大気孔８から導入された空気を介してパージライン４に吸入される」との想定のもとにパージガス挙動のモデル化を図っている。
【０１９０】
ここでキャニスタ空気層３ｂのパージガスは、その内部に含まれるベーパの分圧だけ大気圧よりも高圧となっている。なお本明細書では、大気圧を基準としてそれよりも低い圧力を「負圧」といい、大気圧を基準としてそれよりも高い圧力を「正圧」という。よってキャニスタ空気層３ｂのパージガスの圧力は、正圧となっている。これに対して、大気孔３からの空気の圧力は大気圧であり、パージ実施中のパージライン４の内圧は負圧となっている。
【０１９１】
かかる圧力関係にあっては、正圧であって最も高圧なキャニスタ空気層３ｂのパージガスが、大気圧である大気孔８からの空気を押しのけて、負圧となったパージライン４に優先的に吸入されることとなる。これより、上記（ａ）記載の仮定が導き出されている。ちなみに、この（ａ）記載の仮定事項は、上述の検証試験の結果に示されるように（図２、図３等参照）、実証によって裏付けられている。
【０１９２】
一方、たとえパージライン４に吸入されるパージガスの総流量が無制限であるとしても、キャニスタ空気層３ｂからパージライン４に吸入されるパージガスの流量には、すなわち上記空気層パージ流量Fpgairには自ずと限界がある。上記物理モデルでは、吸気通路６にパージされるガスの流量（総パージ流量Fpgall）のうち、上記空気層パージ流量Fpgairの限界、すなわち上記最大空気層パージ流量Fpgairmxを超えた不足分が、大気孔８からの空気によって補われると想定している。この最大空気層パージ流量Fpgairmxは、大気孔８からの空気を押しのけてキャニスタ空気層３ｂから流出可能なパージガスの流量の限界によって決まる。そしてその最大空気層パージ流量Fpgairmxの値は、図８に示されるようなモデルを想定することで、理論的に求めることができる。
【０１９３】
同図８のモデルでは、大気中に置かれ、開口部を有した容器に上記キャニスタ空気層３ｂを見立てている。上記最大空気層パージ流量Fpgairmxは、キャニスタ空気層３ｂに見立てられた上記容器から噴出されるパージガスの流量として求めることができる。ここで同図８に示されるように、その容器の内圧、すなわちキャニスタ空気層３ｂの内圧を記号「Ｐ」、その容器の外気圧、すなわち大気圧を記号「Ｐ０」、その容器から噴出されるパージガスの流量、すなわち最大空気層パージ流量Fpgairmxを記号「ｑ」で、それぞれ各示すものとする。そして更に、その容器内（キャニスタ空気層３ｂ）のパージガスの密度を記号「ρ」で表すこととすれば、上記流量ｑは、ベルヌーイの定理に基づき、次式によって求められる。
【０１９４】
【数１】
0004400003.tif 000002
ここで同図８のモデルにおける容器内の圧力Ｐは、キャニスタ空気層３ｂのパージガスにおけるベーパ成分の分圧Ｐｘと空気成分の分圧Ｐ０との和として表すことができる。またここで、キャニスタ空気層３ｂに蓄積されているベーパの量（空気層蓄積ベーパ量Mgair ）を記号「Ｇ」で表すこととする。そして更に、キャニスタ空気層３ｂの容積を記号「Ｖ」、同空気層３ｂのパージガスの絶対温度を記号「Ｔ」、同パージガスの質量を記号「Ｍ」、ベーパの分子量を記号「mx」、更に気体定数を記号「Ｒ」で各々示すものとすれば、上記流量ｑ（＝最大空気層パージ流量Fpgairmx）は、更に次式となる。
【０１９５】
【数２】
0004400003.tif 000003
ここでキャニスタ空気層３ｂのパージガス中の空気成分の分圧Ｐ０が常に大気圧であると仮定し、値αを「α＝１／Ｍ」、値βを「β＾２＝２ＲＴ／（mx・Ｍ）」とすると、次式が求められる。
【０１９６】
【数３】
0004400003.tif 000004
なお、本明細書では、任意のパラメータを「Ｘ」、任意の数を「ｎ」としたとき、そのパラメータＸのｎ階のべき乗を「Ｘ＾ｎ」と表記するようにしている。したがって上記「β＾２」は、βの２乗を示している。
【０１９７】
更に、上記流量ｑのうちのベーパ成分の流量、すなわち空気層ベーパ流量Fvpairを記号「ｖ」で表すものとする。このベーパ成分の流量ｖは、パージガス中のベーパ成分の濃度と上記流量ｑとに比例する。そこで、値γ＾２＝２ＲＴ／（mx・Ｍ＾３）とすると、次式が求められる。
【０１９８】
【数４】
0004400003.tif 000005
パージシステム使用時のキャニスタ空気層３ｂの温度変化が十分に小さく、上記絶対温度Ｔが一定値であると仮定すれば、上記各値α、β、γはいずれも、パージシステムに固有の定数と見なすことができる。また上記各値α、β、γの適正値は、試験等により求めることができる。
【０１９９】
なお一般のパージシステムの通常の使用条件では、上記絶対温度Ｔの変化は、上記各流量ｑ、ｖの算出精度に影響を与えるほど大きくはなく、上記の仮定は十分に成立する。ちなみに、そうした絶対温度Ｔの変化の影響が無視できない場合の処置については、後述する（［２−４］、図２０等参照）。
【０２００】
したがって、上記流量ｑ及び流量ｖは、すなわち最大空気層パージ流量Fpgairmx並びに最大空気層ベーパ流量Fvpairmxは、キャニスタ空気層３ｂに蓄積されたベーパの量Ｇ、すなわち空気層蓄積ベーパ量Mgair の関数として表される。これにより、上記（ｂ）の仮定が導き出されている。
【０２０１】
以上説明したように、上記想定された物理モデルによれば、空気層蓄積ベーパ量Mgair が一定であれば、パージ実施中に吸気通路６に放出されるパージガスの各成分と総パージ流量Fpgallとの関係は、図９に示す通りとなる。
【０２０２】
すなわち、上記空気層蓄積ベーパ量Mgair に応じて決まった最大空気層パージ流量Fpgairmxに総パージ流量Fpgallが達するまでは（Fpgall＜Fpgairmx）、吸気通路６へのパージガスの全てがキャニスタ空気層３ｂからのパージガスによって占められる。したがって同図９に示されるように、そのときの空気層パージ流量Fpgairは総パージ流量Fpgallと同一の値となる（Fpgair＝Fpgall）。一方、総パージ流量Fpgallが最大空気層パージ流量Fpgairmxを超えると（Fpgall≧Fpgairmx）、空気層パージ流量Fpgairはその最大空気層パージ流量Fpgairmxにて飽和する（Fpgair＝Fpgairmx）。またそのときのパージガスの流量の不足分（Fpgall−Fpgairmx）は、大気孔８から導入された空気の流量によって補われる。
【０２０３】
なお、空気層ベーパ流量Fvpairは、キャニスタ空気層３ｂのパージガスのベーパ濃度と空気層パージ流量Fpgairとによって求まり、その濃度は、空気層蓄積ベーパ量Mgair によって決まる。したがって空気層蓄積ベーパ量Mgair が一定であれば、同図９に示されるように、空気層ベーパ流量Fvpairは空気層パージ流量Fpgairに比例した値となる。そして、空気層パージ流量Fpgairがその最大流量Fpgairmxに飽和すれば、空気層ベーパ流量Fvpairも自ずと最大空気層ベーパ流量Fvpairmxに飽和する。ちなみに空気層パージのベーパ濃度rvpair、すなわち空気層パージ流量Fpgairに占める空気層ベーパ流量Fvpairの比率は、上記［数３］［数４］の式に基づき算出される最大空気層パージ流量Fpgairmx、及び最大空気層ベーパ流量Fvpairmxの比率として求められる（rvpair←Fvpairmx／Fpgairmx）。
【０２０４】
また上記物理モデルによれば、総パージ流量Fpgallを一定としたときの空気層蓄積ベーパ量Mgair と空気層ベーパ流量Fvpairとの対応関係は、図１０に示す通りとなる。
【０２０５】
すなわち、総パージ流量Fpgallを一定としたときには、同図１０に示されるように、上記［数３］の式に従い、空気層蓄積ベーパ量Mgair の増加に伴って空気層ベーパ流量Fvpairが増大する。ただしその空気層ベーパ流量Fvpairの増大率は、空気層蓄積ベーパ量Mgair の増加に応じて徐々に減少していく傾向にある。
【０２０６】
ちなみに以上説明した物理モデルに基づき求められた空気層パージ流量Fpgair及び空気層ベーパ流量Fvpairの理論値は、発明者らによる実機での検証試験の結果にほぼ合致しており、上記（ｂ）記載の仮定事項についても、実証によって裏付けられている。
【０２０７】
［１−３−２］  定常時のキャニスタ内でのベーパ挙動の物理モデル
続いて、定常時のキャニスタ３内でのベーパの挙動についての物理モデルを、図１１を併せ参照して説明する。この物理モデルは、定常時、すなわち燃料タンク１からのベーパ流入やパージ実施による吸気通路６へのパージガスの流出のないときの、キャニスタ３内でのベーパの挙動を説明するものである。そして本モデルによれば、定常時にキャニスタ空気層３ｂと吸着剤３ａとの間で授受されるベーパの挙動は次の通りとなる。
（ｃ）定常時にキャニスタ空気層３ｂに蓄積されたパージガス中から吸着剤３ａに吸着されるベーパの流量、すなわちベーパ吸着速度Fvpatcは、空気層蓄積ベーパ量Mgair に応じて増大する。
（ｄ）また、ベーパ吸着速度Fvpatcは、吸着剤３ａのベーパの未吸着面積が大きなほど増大する。
（ｅ）定常時に吸着剤３ａから自然に脱離してキャニスタ空気層３ｂのパージガス中に放出されるベーパの流量、すなわち自然脱離速度Fvpctaは、吸着剤蓄積ベーパ量Mgcan に応じて増大する。
【０２０８】
以下に、上記仮定事項（ｃ）〜（ｅ）について、その理論的根拠、及びその詳細を説明する。
吸着剤３ａは、活性炭のような体積比の表面積の大きな多数の粒子の表面にベーパを付着させることでベーパを吸蔵する構成となっている。吸着剤３ａ全体のベーパ吸着可能な表面積は膨大ではあるものの、その吸蔵能力は有限である。そこである程度のベーパが吸着した状態では、吸着剤３ａ全体の表面に、ベーパが既に吸着した部分（ベーパ吸着部）と未だに吸着していない部分（ベーパ未吸着部）とが存在するものとして図１１に示されるようなモデルを提起する。
【０２０９】
このモデルによれば、定常時には、吸着剤３ａのベーパ吸着部からキャニスタ空気層３ｂのパージガスにベーパが徐々に移動し、そのパージガスから吸着剤３ａのベーパ未吸着部にベーパが徐々に移動すると想定される。
【０２１０】
そしてキャニスタ空気層３ｂのパージガス中のベーパ分圧が高ければ、定常時に吸着剤３ａのベーパ未吸着部に移動するベーパ量が増大することは容易に推測できる。上記ベーパ分圧は、吸着剤蓄積ベーパ量Mgcan の増大に伴いそれにほぼ比例して上昇する。このため、上記仮定事項（ｃ）のように、空気層蓄積ベーパ量Mgair の増加に応じてベーパ吸着速度Fvpatcも増大することが推測できる。ちなみに本実施形態では、ベーパ吸着速度Fvpatcは空気層蓄積ベーパ量Mgair に単純に比例するものとして取り扱っている（Fvpatc∝Mgair ）。
【０２１１】
なお厳密には、ベーパ吸着速度Fvpatcと空気層蓄積ベーパ量Mgair とが単純な比例関係となることは実証されていない。しかしながら通常、上記タンク発生ベーパ流量Fvptnk、或いはキャニスタ空気層３ｂや吸着剤３ａから吸気通路６へのパージ実施中のベーパ流量Fvpair、Fvpcan等に比して、上記ベーパ吸着速度Fvpatcは非常に小さな値となる。よって上記比例関係にあるとの仮定に従いベーパ吸着速度Fvpatcを算出しても、実用上は十分である。勿論、更なる検証試験を行い、ベーパ吸着速度Fvpatcと空気層蓄積ベーパ量Mgair と詳細な相関関係を求めて算出に用いれば、ベーパ吸着速度Fvpatcのより厳密な推定が可能となる。
【０２１２】
一方、吸着剤３ａのベーパ未吸着部の表面積が減少すれば、ベーパの吸着能力は一時的に低下する。よって上記仮定事項（ｄ）のように、吸着剤３ａのベーパ未吸着面積が大きなほどベーパ吸着速度Fvpatcが増大することは容易に予測できる。また吸着剤３ａのベーパの最大吸着量VPCANMX 、すなわち吸着剤３ａの吸着表面の全てがベーパによって満たされ、それ以上のベーパの吸着を許容できなくなった飽和時の吸着剤蓄積ベーパ量Mgcan は、試験等により求めることができる。そして上記ベーパの未吸着面積は、その最大吸着量VPCANMX から現状の吸着剤蓄積ベーパ量Mgcan を減算した値に比例する。ちなみに本実施形態では、ベーパ吸着速度Fvpatcは空気層蓄積ベーパ量Mgair に単純に比例するものとして取り扱っている。すなわち、ベーパ吸着速度Fvpatcは、上記最大吸着量VPCANMX から現状の吸着剤蓄積ベーパ量Mgcan を減算した値に比例するものとしている（Fvpatc∝｜VPCANMX −Mgcan ｜）。これについても、上記比例関係は実証されていないものの、上記仮定事項（ｃ）の場合と同様に実用上は十分である。また更なる検証試験によって、ベーパ吸着速度Fvpatcと空気層蓄積ベーパ量Mgair との相関関係を詳細に求めて算出に用いれば、ベーパ吸着速度Fvpatcのより厳密な推定が可能となることは勿論である。
【０２１３】
また更に、定常時に吸着剤３ａからのベーパの自然脱離は、吸着されたベーパに対して一定の確率で生じることが確認されている。このため、上記仮定事項（ｅ）のように、上記自然脱離速度Fvpctaは、吸着剤３ａに吸着されたベーパ量、すなわち吸着剤蓄積ベーパ量Mgcan の増加に応じて増大する。また上記のようにベーパの自然脱離の確率が一定なため、自然脱離速度Fvpctaは吸着剤蓄積ベーパ量Mgcan に比例する（Fvpcta∝Mgcan ）。
【０２１４】
以上のように、定常時のキャニスタ３内でのベーパ挙動を物理モデルに基づいて推定することができる。また非定常時においても、燃料タンク１からのベーパ流入やパージ実施による吸気通路６へのベーパ流出が更に加わるだけで、上記定常時のベーパ挙動は保持されると考えられる。
【０２１５】
ところで、吸着剤３ａは、ベーパの吸着・脱離を繰り返す毎に徐々に劣化し、そのベーパ吸着能力が低下する。こうした劣化は、上記最大吸着量VPCANMX の減少として説明できる。このため、かかる劣化により、上記ベーパ吸着速度Fvpatcの推定値に多少の誤差が生じるおそれはある。そうした場合であれ、吸着剤３ａの劣化の度合いに応じて最大吸着量VPCANMX の値を適宜更新していくようにすれば、そうした劣化に拘わらず、ベーパ吸着速度Fvpatcを正確に推定できる。なお、実機においては劣化による最大吸着量VPCANMX の減少は微量であり、各種エンジン制御にほとんど影響を与えることはなく、上記劣化に対応した処置を行わずとも、実用上ほとんど問題にはならない。
【０２１６】
［１−３−３］  パージ実施中のベーパ挙動の物理モデル
続いて、パージ実施中のベーパ挙動についての物理モデルを説明する。なお、キャニスタ空気層３ｂのパージガスから吸気通路６にパージされるベーパの挙動については、上記［１−３−１］で説明した通りであるため、ここではパージ実施中に吸着剤３から脱離してパージされるベーパの挙動について考察する。
【０２１７】
パージ実施中には、大気孔８から導入された空気の流勢によって、吸着剤３ａに吸着されたベーパが脱離して、吸気通路６にパージされる。このため、パージ実施中に吸着剤３ａから脱離してパージされるベーパの流量、すなわち吸着剤脱離ベーパ流量Fvpcanは、吸着剤３ａ内を流過する空気の流量、すなわち吸着剤内空気流量Fpgcanにほぼ比例する（Fvpcan∝Fpgcan）。
【０２１８】
また吸着剤３ａに吸着されたベーパが多量であるほど、吸着剤３ａから脱離するベーパの流量が増大することは容易に推定できる。また更に、大気孔８から導入された空気と共に吸気通路６にパージされるパージガス（吸着剤脱離パージガス）中のベーパ濃度rvpcanは、吸着剤蓄積ベーパ量Mgcan に応じて一義的に求まることが従来より確認されている（rvpcan←Fnc.｛Mgcan ｝）。
【０２１９】
以上により、以下の結論が導き出される。
（ｆ）吸着剤脱離ベーパ流量Fvpcanは、パージ実施中に大気孔８から導入される空気の流量、すなわち吸着剤内空気流量Fpgcanに比例する。
（ｇ）また吸着剤脱離パージガスのベーパ濃度rvpcanは、吸着剤蓄積ベーパ量Mgcan によって一義的に求められる。すなわち、吸着剤内空気流量Fpgcanを一定としたときの吸着剤脱離ベーパ流量Fvpcanは、吸着剤蓄積ベーパ量Mgcan によって一義的に決まる。
【０２２０】
また更に、上記［１−３−１］にて導出した空気層パージのベーパ挙動（図９参照）を併せ考慮すれば、パージ実施中に吸気通路６に放出されるパージガス中の各成分を推定可能となる。そして、空気層蓄積ベーパ量Mgair 及び吸着剤蓄積ベーパ量Mgcan が各一定であれば、吸気通路６へのパージガスの各成分と総パージ流量Fpgallとの関係は、図１２に示す通りとなる。
【０２２１】
すなわち、上述のように総パージ流量Fpgallが最大空気層パージ流量Fpgairmxを超えると、キャニスタ空気層３ｂａ内からのパージガスの流量（空気層パージ流量Fpgair）はそこで頭打ちとなり、その不足分が大気孔から導入された空気の流量によって補われる。よって、その不足分の流量、すなわち総パージ流量Fpgallと最大空気層パージ流量Fpgairmxとの差分（Fpgall−Fpgair）分の流量が、大気孔８から導入される吸着剤内空気流量Fpgcanとなる。
【０２２２】
このとき吸着剤内空気流量Fpgcanに占めるベーパの濃度rvpcanは、吸着剤蓄積ベーパ量Mgcan が変化しない限り一定である。このため、吸着剤蓄積ベーパ量Mgcan が一定であれば、総パージ流量Fpgallが最大空気層パージ流量Fpgairmxを超えた領域での吸着剤脱離ベーパ流量Fvpcanは、吸着剤内空気流量Fpgcanに比例する。よって上記領域での吸着剤脱離ベーパ流量Fvpcanは、総パージ流量Fpgallの増大に伴い単調増加する。
【０２２３】
［１−３−４］  パージシステム全体でのベーパ挙動の物理モデル
以上を総括すれば、図１３に示されるような、パージシステム全体でのベーパ挙動を表す物理モデルを導き出すことができる。ここで、その物理モデルにおいて、同図１３に示される各パラメータの説明、及びその値の算出にかかる関係式を以下に列挙する。
【０２２４】
（Ａ）タンク発生ベーパ流量Fvptnk
燃料タンク１で発生したベーパのキャニスタ空気層３ｂへの流量（流入率）［ｇ毎秒］。燃料タンク１の内圧の推移等から測定によって求めることもできるが、空気層蓄積ベーパ量Mgair の推定値のずれ率（ずれ量の時間的変化）に応じて推定することもできる。
【０２２５】
（Ｂ）空気層蓄積ベーパ量Mgair
キャニスタ空気層３ｂに蓄積されているベーパの量［ｇ］。タンク発生ベーパ流量Fvptnk、ベーパ吸着速度Fvpatc、自然脱離速度Fvpcta、空気層ベーパ流量Fvpairに応じて所定時間毎にその値を更新する。また空燃比フィードバック補正値の監視により検出される空気層ベーパ流量Fvpairの推定値のずれ量に応じてその値を修正する。
<<関係式>>
・  ΔMgair ←  Fvptnk−Fvpatc＋Fvpcta−Fvpair
ここで「ΔMgair 」は、単位時間（１秒）当たりの空気層蓄積ベーパ量Mgair の更新量を示す。
【０２２６】
（Ｃ）吸着剤蓄積ベーパ量Mgcan
キャニスタ３の吸着剤３ａ内に蓄積されているベーパの量［ｇ］。ベーパ吸着速度Fvpatc、自然脱離速度Fvpcta、吸着剤脱離ベーパ流量Fvpcanに応じて所定時間毎にその値を更新する。
<<関係式>>
・  ΔMgcan ←  Fvpatc−Fvpcta−Fvpcan
ここで「ΔMgcan 」は、単位時間（１秒）当たりの吸着剤蓄積ベーパ量Mgcan の更新量を示す。
【０２２７】
（Ｄ）ベーパ吸着速度Fvpatc
定常時のキャニスタ空気層３ｂから吸着剤３ａに吸着するベーパの流量（単位時間当たりの吸着量）［ｇ毎秒］。空気層蓄積ベーパ量Mgair 及び吸着剤３ａのベーパ未吸着面積（VPCANMX −Mgcan ）に比例する。
<<関係式>>
・  Fvpatc  ←  k1・Mgair ・（VPCANMX −Mgcan ）
ここで「k1」は、所定の定数を示す。
【０２２８】
（Ｅ）自然脱離速度Fvpcta
大気孔８からの空気流入に伴うことなく吸着剤３ａからキャニスタ空気層３ｂへと自然に脱離するベーパの流量［ｇ毎秒］。その値は吸着剤蓄積ベーパ量Mgcan に比例する。
<<関係式>>
・  Fvpcta  ←  k2・Mgcan
ここで「k2」は、所定の定数を示す。
【０２２９】
（Ｆ）空気層ベーパ流量Fvpair
パージ実施中に、キャニスタ空気層３ｂから吸気通路６にパージされるベーパの流量［ｇ毎秒］。その値は、空気層蓄積ベーパ量Mgair と総パージ流量Fpgallとの関数として求められる。
<<関係式>>
・  Fvpair  ←  rvpair・Fpgair
・  rvpair  ←  Fvpairmx／Fpgairmx（＝Fnc.｛Mgair ｝）
・  Fpgair  ←  Fpgall  （ただし、Fpgair≦Fpgairmx）。
【０２３０】
（Ｇ）吸着剤脱離ベーパ流量Fvpcan
パージ実施中に、大気孔８からの空気流入に伴って吸着剤３ａから脱離して吸気通路６にパージされるベーパ流量［ｇ毎秒］。その値は吸着剤内空気流量Fpgcanに比例する。またその比例定数（吸着剤脱離パージのベーパ濃度rvpcanに相当）は、吸着剤蓄積ベーパ量Mgcan によって一義的に求まる。
<<関係式>>
・  Fvpcan  ←  rvpcan・Fpgcan
・  rvpcan  ←  Fnc.｛Mgcan ｝
・  Fpgcan  ←  Fpgall−Fpgairmx  （ただし、Fpgcan≧０）
以上、［１−３−１］、［２−４−２］等を参照。
【０２３１】
（Ｈ）総ベーパ流量Fvpall
パージ実施中に吸気通路６に放出されるベーパの総流量［ｇ毎秒］。その値は、空気層ベーパ流量Fvpairと吸着剤脱離ベーパ流量Fvpcanとの和となる。
<<関係式>>
・  Fvpall  ←  Fvpair＋Fvpcan
以上、［２−４−２］等を参照。
【０２３２】
以上のように、この物理モデルによれば、センサ等による実測結果に頼らずとも、パージシステム内でのベーパ挙動の変化を的確に把握して、パージ実施中のエンジンへの総ベーパ流量Fvpallを正確に推定できる。そして、こうして推定された総ベーパ流量Fvpallを用いることで、空燃比フィードバック制御の更なる精度向上を図ることが可能となる。
【０２３３】
また同物理モデルによれば、上記ベーパ挙動にかかる各パラメータの推移を詳細に常時把握できるため、上記空燃比フィードバック制御以外の各種エンジン制御についても、それら各パラメータの推移の監視のもとに、より細密な制御を行うことができるようになる。
【０２３４】
［２］  物理モデル適用の具体例
［２−１］  空燃比制御装置の全体構成
続いて、上記物理モデルに基づく制御を適用したエンジンの空燃比制御装置の具体例について、図１４を参照してその全体構成を説明する。
【０２３５】
同図１４に示すようにエンジン１０は、燃焼室１１、吸気通路１２及び排気通路１３を備えている。エンジン１０の運転にあたっては燃料タンク３０内に備蓄された燃料（例えばガソリン）が燃料ポンプ３１によって汲み出され、燃料供給通路を通じてデリバリパイプ１２ａに送られた後、インジェクタ１２ｂによって吸気通路１２内に噴射供給される。この吸気通路１２の上流には、アクセルペダル（図示略）の踏み込み操作に基づいて吸気通路１２の流路面積を可変とするスロットルバルブ１２ｃが設けられている。更に吸気通路１２には、吸入空気を浄化するためのエアクリーナ１２ｄ及び吸気通路１２の内圧（吸気通路内圧ＰＭ）を検出するための吸気圧センサ１２ｅが設けられている。
【０２３６】
一方、排気通路１３には、エンジン１０からの排気ガスを浄化するための触媒コンバータ１３ａが設けられ、その上流には、排気ガス中の酸素濃度を検出するための空燃比センサ１３ｂが配設されている。そしてこの空燃比センサ１３ｂの検出信号に応じて、燃焼室１１で燃焼される混合気の空燃比が求められるようになっている。
【０２３７】
ベーパパージシステム２０は、燃料タンク３０から発生するベーパを捕集するキャニスタ４０や、その捕集されたベーパをエンジン１０の吸気通路１２にパージするパージライン７１を備えている。
【０２３８】
このパージシステム２０にあって、燃料タンク３０の天井部分には、燃料タンク３０内の圧力を検知するためのタンク内圧センサ３２と、ブリーザ制御バルブ３３とが設けられている。タンク内圧センサ３２は、燃料タンク３０及び同タンク３０と連通する領域の圧力を検出する。ブリーザ制御バルブ３３はダイアフラム式の差圧弁であり、給油時等、燃料タンクの内圧がブリーザライン３４内の圧力より所定圧以上高くなるときに自律的に開弁してブリーザライン３４を介してベーパをキャニスタ４０に逃がす。
【０２３９】
更に燃料タンク３０は、ブリーザライン３４よりも通路内径の小さなベーパライン３５を介してキャニスタ４０に連通可能となっている。ベーパライン３５とキャニスタ４０との間に設けられたタンク内圧制御バルブ６０は、先のブリーザ制御バルブ３３と同様の機能を有するダイアフラム式差圧弁である。タンク内圧制御バルブ６０内のダイアフラム弁体６１は、燃料タンク３０内の圧力がキャニスタ４０内の圧力より所定圧以上高くなるときのみ、タンク内圧制御バルブ６０を開弁させる。
【０２４０】
キャニスタ４０は、その内部に吸着剤（例えば活性炭）を備え、ベーパを該吸着剤に吸着させて一時的に蓄えた後、大気圧よりも低い圧力下におかれることによって、すなわち負圧状態となることによって、この吸着剤に吸着させたベーパを再離脱させることが可能な構成となっている。キャニスタ４０は、ブリーザライン３４及びベーパライン３５を介して燃料タンク３０と連通可能である他、パージライン７１を経由して吸気通路１２に連通可能であり、更に大気バルブ７０を介して大気導入ライン７２及び大気排出ライン７３にも連通している。
【０２４１】
なお、パージライン７１の途中にはパージ調整バルブ（ＶＳＶ）７１ａが設けられている。このＶＳＶ７１ａは単なる開閉弁ではなく、全閉状態（開度０％）から全開状態（開度１００％）まで任意に開度調節可能なタイプであり、外部からのデューティ制御によって駆動する。
【０２４２】
また、エアクリーナ１２ｄに連通する大気導入ライン７２の途中には、大気導入バルブ７２ａが設けられている。
大気バルブ７０内には、各々が異なる機能を有する二つのダイアフラム弁体７４，７５が設けられている。第１のダイアフラム弁体７４は、その背面側の空間７４ａがパージライン７１と連通しており、パージライン７１が所定圧以下の負圧状態になると開弁して大気導入ライン７２からキャニスタ４０内への外気の流入を許容する。一方、第２のダイアフラム弁体７５は、キャニスタ４０内が所定圧以上の正圧に達すると開弁してキャニスタ４０から大気排出ライン７３へ余分な空気を排出させる。
【０２４３】
キャニスタ４０の内部は仕切板４１によって、第１吸着剤室４２と第２吸着剤室４３とに区画されている。両吸着剤室４２、４３は吸着剤（活性炭）で満たされるも、両室はキャニスタ底部（図１４では右側）において通気性フィルタ４４を介して連通している。燃料タンク３０は、一方ではベーパライン３５及びタンク内圧制御バルブ６０を介して、他方ではブリーザ制御バルブ３３及びブリーザライン３４を介して第１吸着剤室４２に連通可能となっている。また、大気導入ライン７２及び大気排出ライン７３は大気バルブ７０を介して第２吸着剤室４３に連通可能となっている。そして、ＶＳＶ７１ａを備えたパージライン７１は、キャニスタ４０の第１吸着剤室４２と、吸気通路１２のスロットルバルブ１２ｃ下流位置とを連結しており、ＶＳＶ７１ａの開弁動作に応じて第１吸着剤室４２とスロットルバルブ１２ｃ下流位置とを連通する。
【０２４４】
なお第１吸着剤室４２には、上記ベーパライン３３、ブリーザライン３４、及びパージライン７１の各開口するキャニスタ４０の天井部分と吸着剤とを隔てるキャニスタ空気層４５が形成されている。よってベーパライン３５やブリーザライン３４から導入されたベーパは、一旦はキャニスタ空気層４５のパージガスに混入された後、徐々に第１吸着剤室４２の内部の吸着剤に吸着される。このため給油時等の燃料タンク３０からの多量のベーパ流入時であれ、そのキャニスタ空気層４５が緩衝となり、吸着剤の劣化が抑制される。
【０２４５】
また、大気バルブ７０を構成する第２のダイアフラム弁体７５が開弁してキャニスタ４０内の余分な空気が大気排出ライン７３から排出される場合でも、キャニスタ空気層４５のパージガス中に蓄積されたベーパは、第２吸着剤室４３を通過する際にその内部の吸着剤に吸着される。
【０２４６】
加えてベーパパージシステム２０には、タンク内圧制御バルブ６０（またはべーパ通路３５の一端部）とキャニスタの第２吸着剤室４３とを連絡するように負圧導入用のバイパスライン８０が設けられている。このバイパスライン８０の途中には、バイパス制御バルブ８０ａが設けられている。このバイパス制御バルブ８０ａの開弁時には、バイパスライン８０及びベーパライン３５を介して第２吸着剤室４３と燃料タンク３０とが直接連通する。
【０２４７】
更に、エンジン１０及びベーパパージシステム２０は、エンジン制御系及びパージ制御系としての電子制御装置（ＥＣＵ）５０を備えている。ＥＣＵ５０には、上記吸気圧センサ１２ｅ、タンク内圧センサ３２のほかに、エンジン回転速度（ＮＥ）センサや気筒判別センサ等のエンジン１０の運転制御に必要な各種センサが直接的または間接的に接続されている。またＥＣＵ５０には、インジェクタ１２ｂ、燃料ポンプ３１、ＶＳＶ７１ａ、大気導入バルブ７２ａ及びバイパス制御バルブ８０ａが、それぞれの駆動回路を介して接続されている。
【０２４８】
ＥＣＵ５０は、各センサから提供される各種情報に基づき、空燃比のフィードバック制御、燃料噴射量の制御、点火タイミング制御等のエンジン制御を実行する。また、ＥＣＵ５０は、タンク内圧センサ３２からの出力信号を認識しつつ、ＶＳＶ７１ａ、大気導入バルブ７２ａ及びバイパス制御バルブ８０ａを適宜開閉制御することにより、ベーパパージ制御やパージシステムの自己診断（即ち上記パージ経路の漏れ診断等）を行う。
【０２４９】
なお、ＶＳＶ７１ａの開度調節は、駆動回路からＶＳＶ７１ａに出力される駆動信号のデューティ比を制御することにより行われる。具体的には、前記デューティ比が０％のときにＶＳＶ７１ａは全閉となり、前記デューティ比が１００％のときにＶＳＶ７１ａは全開となる。なお、本パージシステム２０のＶＳＶ７１ａは、吸気通路内圧ＰＭが一定の条件下では、キャニスタ４０から吸気通路１２にパージされるガスの流量（総パージ流量Fpgall）と上記デューティ比とが比例関係となるように構成されている。ちなみに、上記デューティ比は、ＶＳＶ７１ａの実開度と一義的に対応する制御パラメータであるため、以下では上記デューティ比を「ＶＳＶ開度Dvsv」と記して説明する。
【０２５０】
（ベーパパージシステムにおけるベーパパージ処理の概要）
燃料タンク３０内の燃料が蒸気化しその蒸気圧が所定圧以上に達すると、タンク内圧制御バルブ６０が自律開弁して燃料タンク３０からキャニスタ４０内へベーパの流入が許容される。また、例えば燃料給油時のように、ベーパの蒸気圧が燃料タンク３０内で急激に高まるような場合には、ブリーザ制御バルブ３３が自律開弁して、燃料タンク３０からキャニスタ４０内へのより大量のベーパの流入が許容される。キャニスタ４０内に流入したベーパは、一旦はキャニスタ空気層４５のパージガスに混入された後、キャニスタ４０内の吸着剤に徐々に吸着される。
【０２５１】
その後、エンジン１０の冷却水温が所定のパージ開始水温に達するなど、エンジン運転条件が所定の条件を満たすと、ＥＣＵ５０からの制御信号に基づき、閉じていたＶＳＶ７１ａが開弁される。パージライン７１を介して吸気通路１２からキャニスタ４０内に吸引負圧が導かれ、キャニスタ４０内に蓄積されたベーパを含んだパージガスが吸気通路１２にパージされる。
【０２５２】
更にパージされるガスの流量（総パージ流量Fpgall）が一定の流量以上となると、大気導入バルブ７２ａの開弁状態が維持され、大気導入ライン７２を通じてエアクリーナ１２ｄからキャニスタ４０内に新気が導入される。この負圧及び新気の導入によって吸着剤からベーパが離脱し、パージライン７１を介して吸気通路１２にパージされる。よって本パージシステム２０では、大気導入ライン７２、大気導入バルブ７２ａ、大気バルブ７０等が上記「大気孔」に相当する構成となっている。
【０２５３】
［２−２］  パージ制御の概要
続いて、本制御装置でのパージ制御の概要について、図１５を併せ参照して説明する。
【０２５４】
本制御装置では、ＥＣＵ５０は、上記のベーパパージ処理を行いつつ、インジェクタ１２ｂからの燃料噴射量（噴射時間）ＴＡＵの調整に基づいて、燃焼室１１で燃焼される混合気の空燃比を所望とする目標値（例えば理論空燃比）に保持する処理を行っている。ここでＥＣＵ５０は、上記物理モデルに基づき推定される総ベーパ流量Fvpallに応じて上記燃料噴射量を補正することで、ベーパパージの影響分を見込んだ空燃比制御の適合を図っている。更にＥＣＵ５０は、上記総ベーパ流量Fvpallの推定精度の保持や、空燃比制御へのベーパパージの影響緩和のための各種処理を併せ実施することで、上記空燃比制御の適合を更に充実させている。
【０２５５】
図１５は、こうした空燃比制御へのベーパパージの適合にかかる処理内容の概要を示す「基本ルーチン」である。本ルーチンの処理は、エンジン１０の運転中、ＥＣＵ５０によって繰り返し実行される。なお本ルーチンは、かかる処理の全体像を理解し易い態様で示したものであり、実際のＥＣＵ５０の処理手順に完全に一致するものではない。
【０２５６】
まずＥＣＵ５０は、同図１５のステップ１００に示すように、ＶＳＶ７１ａの開度（デューティ比）Dvsvの算出処理を行っている。ここでは、上記物理モデルに基づいて、空燃比制御に与える影響を緩和可能な範囲内に上記総ベーパ流量Fvpallを調整すべく、ＶＳＶ開度Dvsvを設定するようにしている。この処理の詳細については、後記［２−８］で説明する。
【０２５７】
そしてＥＣＵ５０は、続くステップ２００において、上記物理モデルに基づいて、現状の総ベーパ流量Fvpallを推定し、その推定値に応じてパージ補正量を算出する。このときＥＣＵ５０は、上記ステップ１００で算出されたＶＳＶ開度Dvsvに基づき把握される上記総パージ流量Fpgallと、上述の各物理状態量（Mgair 、Mgcan 等）とに基づいて、総ベーパ流量Fvpallを推定する。この処理の詳細については、後記［２−４］で説明する。
【０２５８】
更に続くステップ３００においてＥＣＵ５０は、ここで算出されたパージ補正量に応じて、インジェクタ１２ｂからの燃料噴射量ＴＡＵを算出する。そしてステップ４００において、その算出された燃料噴射量ＴＡＵに応じてインジェクタ１２ｂを駆動制御して燃料噴射を実行する。ここでの燃料噴射量ＴＡＵの算出にかかる処理の詳細は、後記［２−５］で説明する。
【０２５９】
またＥＣＵ５０は、ステップ５００に示すように、上記物理モデルに従って、上記各物理状態量の値の定期更新にかかる処理を行う。この定期更新処理により、上記各物理状態量は、パージシステム２０内でのベーパ挙動の変動に応じた適正な値に保持される。なおこの定期更新処理の詳細については、下記［２−３］で説明する。
【０２６０】
更にＥＣＵ５０は、ステップ６００に示すように、パージ実施中の空燃比フィードバック補正項（以下「空燃比Ｆ／Ｂ補正項」）のずれに応じて、上記各物理状態量の誤差を把握し、それらの値を修正する処理を併せ行っている。この修正処理によって、上記各物理状態量は適正な値に保たれる。この処理の詳細は、後記［２−７］で説明する。
【０２６１】
［２−３］  物理モデルに基づく各物理状態量の定期更新処理（図１５のＳ５００）
続いて、本制御装置における上記各物理状態量の定期更新にかかるＥＣＵ５０の処理の詳細を説明する。
【０２６２】
上述したように上記物理モデルによれば、空気層蓄積ベーパ量Mgair は単位時間に、燃料タンク３０から流入したベーパの流量（タンク発生ベーパ流量Fvptnk）分増加する。また空気層蓄積ベーパ量Mgair は単位時間に、キャニスタ空気層４５と吸着剤４２との間で授受されるベーパの流量分、増減する。詳しくは、ベーパ吸着速度Fvpatc分減少し、自然脱離速度Fvpcta分増加する。またパージ実施中であれば、空気層蓄積ベーパ量Mgair は単位時間に、空気層ベーパ流量Fvpair分減少する。
【０２６３】
更に上記物理モデルによれば、吸着剤蓄積ベーパ量Mgcan は単位時間に、ベーパ吸着速度Fvpatc分増加し、自然脱離速度Fvpcta分減少する。またパージ実施中であれば、吸着剤蓄積ベーパ量Mgcan はその単位時間に、吸着剤脱離ベーパ流量Fvpcan分減少する。
【０２６４】
このため、上記各蓄積ベーパ量の単位時間の変化量ΔMgair 、ΔMgcan は、同図１５に併せ示される数式によって表される。ちなみに、ベーパ吸着速度Fvpatcは、空気層蓄積ベーパ量Mgair 及び吸着剤のベーパ未吸着面積に比例するパラメータとして、自然脱離速度Fvpctaは、吸着剤蓄積ベーパ量Mgcan に比例するパラメータとして、各求められることも上述した通りである（以上、［１−３−２］、［１−３−４］及び図１３等を参照）。よって、所定時間Ｔｓ［秒］毎に定期更新処理を実施するとすれば、処理毎の各蓄積ベーパ量Mgair 、Mgcan の更新量はそれぞれ、上記単位時間の変化量ΔMgair 、ΔMgcan の所定時間Ｔｓにおける時間積分値となる。
【０２６５】
なお、本制御装置では、ＥＣＵ５０は、こうした定期更新処理を単位時間（１秒）毎に実施して、各蓄積ベーパ量Mgair 、Mgcan の値を更新している。このため、本制御装置での本処理時の各蓄積ベーパ量Mgair 、Mgcan の更新量は、上記単位時間の変化量ΔMgair 、ΔMgcan と同一の値となる。
【０２６６】
［２−４］  パージ補正量の算出処理（図１５のＳ２００）
次に、本制御装置でのパージ補正量の算出処理の詳細について、図１６〜図２２を併せ参照して説明する。
【０２６７】
本制御装置では上述したように、上記物理モデルに従い、上記総パージ流量Fpgallと上記各物理状態量とに基づいて総ベーパ流量Fvpallを推定し、その推定値からパージ補正量を求めている。図１６は、こうした総ベーパ流量Fvpallの推定にかかる各パージ流量の算出ロジックを、また図１７は、その推定にかかる各ベーパ流量の算出ロジックを、それぞれ示している。以下、本制御装置でのＥＣＵ５０による総ベーパ流量Fvpallの算出処理を、それら図１６、１７を参照して説明する。
【０２６８】
［２−４−１］  各パージ流量の算出処理（図１６）
まずＥＣＵ５０は、吸気圧センサ１２ｅの検出する吸気通路内圧ＰＭと、上記ＶＳＶ７１ａへの指令信号に基づき把握されるＶＳＶ開度Dvsvとに基づいて、総パージ流量Fpgallを算出する。ここでは詳しくは、以下の演算処理にて総パージ流量Fpgallが算出されている。
【０２６９】
ＶＳＶ７１ａを全開（ＶＳＶ開度Dvsvを１００％）としたときの所定の吸気通路内圧ＰＭでの総パージ流量Fpgall、すなわち所定の吸気通路内圧ＰＭにおける総パージ流量の最大値（最大総パージ流量）Fpgmx は一義的に求めることができる。また上述のように本パージシステムでは、吸気通路内圧ＰＭ一定の条件下では、ＶＳＶ開度Dvsvと総パージ流量Fpgallとが比例する構成となっている。
【０２７０】
本制御装置では、試験等により求められた吸気通路内圧ＰＭと最大総パージ流量Fpgmx との関係が、図１８に例示するような演算マップとして、ＥＣＵ５０のメモリ内に予め記憶されている。そこでＥＣＵ５０は、吸気通路内圧ＰＭの検出値からその演算マップを用いて上記最大総パージ流量Fpgmx を求め、その最大総パージ流量Fpgmx とＶＳＶ開度（デューティ比）Dvsvとの積算によって総パージ流量Fpgallを算出している。
【０２７１】
続いてＥＣＵ５０は、その総パージ流量Fpgallに占める各パージ成分の流量、すなわち空気層パージ流量Fpgairと吸着剤内空気流量Fpgcanとをそれぞれ算出する。それらの算出は、詳しくは以下の態様で行われる。
【０２７２】
上述したように、総パージ流量Fpgallが最大空気層パージ流量Fpgairmxに達するまでは、総パージ流量Fpgallのほぼ全てが空気層パージ流量Fpgairによって占められる。またその最大空気層パージ流量Fpgairmxが、空気層蓄積ベーパ量Mgair によって一義的に決まることも上述した通りである（以上、［１−３−１］、図９等を参照）。
【０２７３】
ＥＣＵ５０のメモリ内には、図１９に例示するような、試験等によって求められた空気層蓄積ベーパ量Mgair と最大空気層パージ流量Fpgairmxとの対応関係を示す演算マップが、予め記憶されている。そこでＥＣＵ５０は、まずその演算マップを用いて最大空気層パージ流量Fpgairmxを算出し、その算出値と上記求められた総パージ流量Fpgallとの対比のもとに、上記各パージ流量Fpgair、Fpgcanを求めている。すなわち、総パージ流量Fpgallが最大空気層パージ流量Fpgairmx未満であれば、空気層パージ流量Fpgairを総パージ流量Fpgallと同値に設定するとともに、吸着剤内空気流量Fpgcanを「０」とする。また総パージ流量Fpgallが最大空気層パージ流量Fpgairmx以上であれば、空気層パージ流量Fpgallを最大空気層パージ流量Fpgairmxと同値に設定する。それと共に、総パージ流量Fpgallから最大空気層パージ流量Fpgairmxを減算した値を、吸着剤内空気流量Fpgcanの値として設定する。以上が図１６に示される各パージ流量の演算処理の内容である。
【０２７４】
なお、最大空気層パージ流量Fpgairmxの理論式である上記［数３］に示されるように、同流量Fpgairmxは、キャニスタ空気層４５のパージガスの絶対温度Ｔにも、ある程度依存するパラメータである。本制御装置では、通常の使用条件では、上記絶対温度Ｔの変化は小さく、算出精度にほとんど影響しないものと見なして上記流量Fpgairmxの算出を行っている。しかしながら、パージシステムの構成やその使用状況等によっては、そうした上記絶対温度Ｔの影響を無視できない場合もある。その場合、以下の態様で同流量Fpgairmxを算出することで、その算出精度の低下を好適に回避できる。
【０２７５】
上記理論式［数３］に示されるように、最大空気層パージ流量Fpgairmxは、上記絶対温度Ｔの平方根に比例する。よって、上記試験等による上記演算マップ（図１９）の作成にあたり測定、或いは想定されたキャニスタ空気層４５のパージガスの絶対温度Ｔｓ［Ｋ］と、流量Fpgairmx算出時の同パージガスの絶対温度Ｔｎ［Ｋ］とを各求めておく。そして、それらの絶対温度の比（Ｔｒ／Ｔｓ）の平方根を、上記演算マップを用いて算出された値に乗算すれば、かかる温度の影響を上記流量Fpgairmxの算出値に織り込むことができる。以下にそうした算出処理態様の一例を説明する。
【０２７６】
キャニスタ空気層４５のパージガスの温度は、吸気通路１２に吸入される空気の温度（吸気温度）ｔｈａと、ほぼ同じであると考えられる。また多くの車載エンジンの制御系では、吸気温度ｔｈａを監視しており、その吸気温度ｔｈａの値は摂氏温度［℃］で示されている。ここで、最大空気層パージ流量Fpgairmxの算出用の演算マップの作成時の想定温度をＴｓ［℃］とすれば、上記絶対温度の比ｋｔｈａは、図２０右上に記した式にて表される（ｋｔｈａ＾２←（ｔｈａ＋２７３）／（Ｔｓ＋２７３））。そしてその比ｋｔｈａと吸気温度ｔｈａとの対応関係は、同図２０に併せ示されるグラフの通りとなる。よって、ＥＣＵ５０のメモリ内に予め記憶しておいた上記対応関係を示す演算マップを利用して、吸気温度ｔｈａに応じて上記比ｋｔｈａを流量の温度補正係数として算出する。そしてその温度補正係数ｋｔｈａを、上記図１９に例示された演算マップにより算出された値に乗算して最大空気層パージ流量Fpgairmxを求めるようにする。勿論、その流量Fpgairmxの算出の都度、上記図２０に示した関係式により温度補正係数ｋｔｈａを毎回演算しても、同じ結果が得られる。
【０２７７】
［２−４−２］  各ベーパ流量の算出処理（図１７）
さらにＥＣＵ５０は、上記算出された各パージ流量、すなわち総パージ流量Fpgall、空気層パージ流量Fpgair、吸着剤内空気流量Fpgcan、及び最大空気層パージ流量Fpgairmxを用い、図１７に示される態様で各ベーパ流量の算出処理を行う。以下にその算出処理の詳細を説明する。
【０２７８】
上記［１−３−１］で述べたように、空気層パージのベーパ挙動は下記の特性を有している。
・  空気層パージのベーパ濃度rvpairは、上述の理論式［数４］に基づき求められる最大空気層ベーパ流量Fvpairmxと、上記算出された最大空気層パージ流量Fpgairmxとの比（Fvpairmx／Fpgairmx）として求められる。
・  最大空気層ベーパ流量Fvpairmxは、同理論式［数４］に従い、空気層蓄積ベーパ量Mgair から一義的に求まる。
【０２７９】
よってＥＣＵ５０は、まず空気層蓄積ベーパ量Mgair から最大空気層ベーパ流量Fvpairmxを求め、それと上記算出された最大空気層パージ流量Fpgairmxとの比として空気層パージのベーパ濃度rvpairを算出している。そしてそのベーパ濃度rvpairと上記算出された最大空気層ベーパ流量Fvpairmxとを積算し、空気層ベーパ流量Fvpairを求めている。本制御装置ではＥＣＵ５０のメモリ内には、空気層蓄積ベーパ量Mgair と最大空気層ベーパ流量Fvpairmxとの対応関係を示す演算マップが記憶されており、その演算マップを用いて、最大空気層ベーパ流量Fvpairmxを求めている。図２１は、そうした演算マップの一例を示している。
【０２８０】
また上記［１−３−２］で述べたように、吸着剤脱離パージのベーパ濃度rvpcanは、吸着剤蓄積ベーパ量Mgcan より一義的に求められる。そこでＥＣＵ５０は、まず吸着剤蓄積ベーパ量Mgcan からそのベーパ濃度rvpcanを求めている。本制御装置では、予めＥＣＵ５０のメモリ内に記憶された吸着剤蓄積ベーパ量Mgcan とベーパ濃度rvpcanとの対応関係を示す演算マップを用いて、同ベーパ濃度rvpcanの算出処理を行っている。図２２は、そうした演算マップの一例を示している。そしてＥＣＵ５０は、先に算出された吸着剤内空気流量Fpgcanとそのベーパ濃度rvpcanとを積算して、吸着剤脱離ベーパ流量Fvpcanを求めている。
【０２８１】
更にＥＣＵ５０は、それら求められた空気層ベーパ流量Fvpairと吸着剤脱離ベーパ流量Fvpcanとの和として、総ベーパ流量Fvpallを求めている（Fvpall←Fvpair＋Fvpcan）。以上が図１７に示される各ベーパ流量の算出処理の内容である。
【０２８２】
なお上記理論式［数４］に示されるように、最大空気層ベーパ流量Fvpairmxも、やはりキャニスタ空気層４５のパージガスの温度について上記最大空気層パージ流量Fpgairmxと同様の依存性を持ったパラメータである。そうした温度依存性が問題となる場合には、最大空気層パージ流量Fpgairmxの場合と同様にして求められた上記温度補正項ｋｔｈａを、演算マップより求められた値に乗算して最大空気層ベーパ流量Fvpairmxを求めれば、その問題を回避できる。
【０２８３】
以上説明した算出処理の後、ＥＣＵ５０は、上記求められた総ベーパ流量Fvpallに応じてパージ補正値ｆｐｇを算出する。このパージ補正値ｆｐｇは、インジェクタ１２ｂからの単位時間（例えば１秒）当たりの燃料噴射量Ｑｆｉｎについて、そのベーパパージの影響分に相当する補正項である。よって、上記物理モデル（図１３等参照）に基づくベーパパージ処理が実施されているときのパージ補正値ｆｐｇは、上記そうベーパ流量Fvpallを符号反転した値となる（ｆｐｇ←−Fvpall）。
【０２８４】
［２−５］  燃料噴射量の算出処理（図１５のＳ３００）
次に、本制御装置での燃料噴射量の算出処理の詳細について説明する。
本制御装置では、ＥＣＵ５０は、次式に概ね従って、インジェクタ１２ｂからの単位時間当たりの燃料噴射量Ｑｆｉｎ［ｇ毎秒］を求めている。
<<燃料噴射量の演算式>>
・  Ｑｆｉｎ  ←  Ｑｂａｓｅ＋ｆａｆ＋ＫＧ＋ｆｐｇ
上式において、「Ｑｂａｓｅ」は基本燃料噴射量［ｇ毎秒］であり、ＥＣＵ５０のメモリ内に予め記憶された所定の演算マップを用い、エンジン回転速度ＮＥやエンジン負荷Ｑに応じて算出されている。また「ｆａｆ」は空燃比フィードバック補正値（以下「空燃比Ｆ／Ｂ補正値」と表記）を、「ＫＧ」は空燃比学習値を各示している。これら空燃比Ｆ／Ｂ補正値ｆａｆ及び空燃比学習値ＫＧは、下記の空燃比フィードバック制御の処理においてその値がそれぞれ設定される。
【０２８５】
ここで本制御装置における空燃比フィードバック制御の概要を、図２３を参照して説明する。この空燃比フィードバック制御は、燃焼室１１で燃焼される混合気の空燃比を目標とする空燃比（例えば理論空燃比）とするための制御であり、上記空燃比Ｆ／Ｂ補正値ｆａｆ及び空燃比学習値ＫＧによる燃料噴射量Ｑｆｉｎの補正を通じてＥＣＵ５０により処理される。
【０２８６】
同図２３には、上記空燃比センサ１３ｂの検出結果に応じた空燃比Ｆ／Ｂ補正値ｆａｆの推移が例示されている。同図２３にその推移の示されるパラメータ「ＸＯ」は、上記空燃比センサ１３ｂの検出信号より把握される空燃比の測定値に基づき、その測定値が上記目標値よりも小さいか、大きいかに応じて、二値化したものである。よって「ＸＯ」は、測定結果に基づく現状のエンジン１０の空燃比が、目標空燃比よりもリーンな状況にあるか、リッチな状況にあるかを示す実空燃比の指標値である。
【０２８７】
ＥＣＵ５０は、この実空燃比の指標値ＸＯに応じて空燃比Ｆ／Ｂ補正値ｆａｆの値を操作し、燃料噴射量Ｑｆｉｎを調整することで、エンジン１０の実際の空燃比をその目標値近傍に保持している。より詳細には、そうした空燃比Ｆ／Ｂ補正値ｆａｆの操作は、次の態様で行われる。
【０２８８】
すなわち、同図２５の時刻ｔ１のように、上記指標値Ｘ０により把握される実空燃比がリッチからリーンに転じると、ＥＣＵ５０は、空燃比Ｆ／Ｂ補正値ｆａｆを一時に所定量増大させ、燃料噴射量Ｑｆｉｎをその分増量する。そして上記実空燃比がリーンからリッチに反転するまでの期間（時刻ｔ１〜時刻ｔ２までの期間）、空燃比Ｆ／Ｂ補正値ｆａｆの値を所定率ずつ徐増していく。一方、時刻ｔ２のように上記実空燃比がリーンからリッチに転じると、ＥＣＵ５０は空燃比Ｆ／Ｂ補正値ｆａｆの値を、今度は一時に所定量低減させる。そして更に上記実空燃比がリーンからリッチに再反転するまでの期間（時刻ｔ２〜時刻ｔ３までの期間）、空燃比Ｆ／Ｂ補正値ｆａｆの値を所定率ずつ徐減していく。以上により、空燃比をその目標値近傍に保持すべく、燃料噴射量Ｑｆｉｎのフィードバック補正が行われる。以後、実空燃比のリーン／リッチ間の移行時の空燃比Ｆ／Ｂ補正値ｆａｆの一時増減を空燃比Ｆ／Ｂ補正値ｆａｆの「スキップ」という。また、その移行時から更に実空燃比のリーン／リッチが再反転するまでの空燃比Ｆ／Ｂ補正値ｆａｆの徐減増を空燃比Ｆ／Ｂ補正値ｆａｆの「積分」といい、それがなされている期間を「積分期間」という。
【０２８９】
なおＥＣＵ５０は、かかる空燃比Ｆ／Ｂ補正値ｆａｆ推移より、その増減の中心値（空燃比Ｆ／Ｂ中心値）fafav を求めている。そしてＥＣＵ５０は、所定のエンジン運転条件が満たされたときのその空燃比Ｆ／Ｂ中心値fafav の値に基づいて、その中心値fafav の値がほぼ「０」となるように、空燃比学習値ＫＧを求め、記憶保持している。この空燃比学習値ＫＧは、エンジン回転速度ＮＥやエンジン負荷Ｑなどのエンジン運転状態に応じて区分けされた複数の領域毎にその値が別途求められ、記憶保持されている。そしてこれにより、エンジン運転状況の遷移時であれ、上記空燃比Ｆ／Ｂ補正値ｆａｆの積分による追従を待たずして、早急に所望とする空燃比を確保できる。なお、そうした空燃比学習値ＫＧの設定にかかる上記所定のエンジン運転条件には通常、ベーパパージ処理の実施やエンジン運転状況の変化を始めとする不安定化要素の十分に小さな、空燃比の安定したエンジン運転条件が選定される。
【０２９０】
またＥＣＵ５０は、その空燃比Ｆ／Ｂ中心値fafav の徐変値fafsm 、すなわち空燃比Ｆ／Ｂ中心値fafav の変化に追従して、徐増若しくは徐減する値を求めており、外乱の影響を排した上記空燃比Ｆ／Ｂ補正値ｆａｆの推移を掌握するようにしている。
【０２９１】
さて上述のようにパージ実施中には、ベーパパージ処理に伴い吸気通路１２に放出されるベーパが燃焼室１１内で燃焼される混合気に混入されるため、本来は、ベーパパージ処理により混入されたベーパ分、混合気の空燃比は減少する（リッチとなる）。ただし本制御装置では、上記算出式に示されるようにパージ補正値ｆｐｇによって、燃料噴射量Ｑｆｉｎはその混入されたベーパ分減量されるようになっている。よって、上記総ベーパ流量Fvpallが適正に推定されており、パージ補正値ｆｐｇに適切な値が設定されていれば、パージの実施の有無や総パージ流量Fpgallの変化等のパージ状況に変化が生じようとも、空燃比Ｆ／Ｂ補正値ｆａｆには何ら影響しない。なお、これを逆に取れば、パージ状況の変化によって空燃比Ｆ／Ｂ補正値ｆａｆの値にずれが生じれば、上記パージ補正値ｆｐｇに、そしてひいては総ベーパ流量Fvpallの推定に誤りがあると考えられる。
【０２９２】
なお、ＥＣＵ５０は、こうして算出された燃料噴射量Ｑｆｉｎより、それをエンジン回転速度ＮＥ等に応じて各インジェクタ１２ｂの１噴射あたりの噴射時間ＴＡＵに換算する。そしてＥＣＵ５０は、その噴射時間ＴＡＵに基づいて各インジェクタ１２ｂに指令信号を出力し、エンジン１０に燃料を噴射供給させる。以上により、ベーパパージの影響を加味した上での、燃料噴射量の調整に基づく空燃比フィードバック制御が行われている。
【０２９３】
［２−６］  各物理状態量の初期化処理
以上説明したように、上記物理モデル（図１３等参照）に従えば、上記各物理状態量（タンク発生ベーパ流量Fvptnk、空気層蓄積ベーパ量Mgair 、吸着剤蓄積ベーパ量Mgcan ）より総ベーパ流量Fvpallを推定し、パージ補正値ｆｐｇを適正に求めることができる。また同モデルに従い、上記［２−３］で述べた定期更新処理を行えば、キャニスタ４０でのベーパ挙動の変化に応じて上記各物理状態量を現状に則した適正値に保持できる。しかしながら、エンジン１０始動後の始めてのパージ実施時のように、上記各物理状態量が不明なときには、そうした物理モデルに基づく総パージ流量Fvpallの推定等も不能となる。
【０２９４】
そこで本制御装置では、上記各物理状態量が不明なときには、下記態様でそれらの値を初期化する処理、すなわちそれらの初期値を算出する処理を行っている。ここではまず、そうした初期化にかかる処理の詳細を、図２４〜図２６を併せ参照して説明する。
【０２９５】
［２−６−１］  初期化完了前のベーパパージ処理
さて本制御装置では、ＥＣＵ５０は、上記物理モデルに従い推定される総ベーパ流量Fvpallに加え、パージ状況の変化に応じた上記空燃比Ｆ／Ｂ補正値ｆａｆの推移に基づき算出される総ベーパ流量実測値Fvpsを求めている。そして上記各物理状態量が不明で総ベーパ流量Fvpallの推定が不能な上記初期化の完了前には、総ベーパ流量Fvpallに代えてその総ベーパ流量実測値Fvpsを用いて上記パージ補正値ｆｐｇを求めている。そして初期化完了前には、パージ実施に伴う空燃比Ｆ／Ｂ補正値ｆａｆの推移を監視しつつ、そのずれを所定範囲内に抑えるように総パージ流量Fpgallを調整している。
【０２９６】
図２４は、そうした初期化完了前の制御態様例を示している。以下、同図２４を例として、初期化完了前の総パージ流量Fpgallの調整（ＶＳＶ開度Dvsvの調整）、及び上記総ベーパ流量実測値Fvpsの算出にかかるＥＣＵ５０の各処理について説明する。
【０２９７】
さて同図２４の例ではエンジン始動後、暖機が完了する、或いは空燃比Ｆ／Ｂ補正値ｆａｆの値が安定する（その中心値fafav が「０」近傍に保持されている）等、パージ実施に必要な諸条件が時刻ｔ０に成立したものとする。この時刻ｔ０よりＥＣＵ１０は、それまで全閉に保持されていたＶＳＶ７１ａを徐々に開弁してパージを開始する。これにより、その時刻ｔ０以降、ＶＳＶ７１ａの開弁に応じて総パージ流量Fpgallも徐々に増加する。
【０２９８】
なお、上記総ベーパ流量実測値Fvpsの初期値、すなわちエンジン始動時の値は「０」に設定されており、よってパージ補正値ｆｐｇも「０」となっている。このためその時刻ｔ０以降、総パージ流量Fpgallの増加に伴う吸気通路１２へのベーパ流入量の増加分を補償すべく、空燃比Ｆ／Ｂ補正値ｆａｆはその値が減少する方向にずれて行く。ちなみに同図２４では、総ベーパ流量実測値Fvpsについて、その符号を反転した値の推移を示している。
【０２９９】
本制御装置ではＥＣＵ５０は、次の２つのしきい値α、βを用いて、空燃比Ｆ／Ｂ補正値ｆａｆに、パージ実施に伴う有意なずれが生じているか否かを検知している。まずＥＣＵ５０は、実空燃比ＸＯのリーン／リッチ移行時のスキップ処理後の空燃比Ｆ／Ｂ中心値fafav の絶対値が上記しきい値αを超えていれば（fafav ＜−α、またはfafav ＞α）、上記ずれ発生と判断する。このときＥＣＵ５０は、その時点で、上記総ベーパ流量実測値Fvpsを所定値加増する、若しくは所定値低減することで、そのずれを補正する。
【０３００】
またＥＣＵ５０は、上記積分期間中の空燃比Ｆ／Ｂ補正値ｆａｆの絶対値が上記しきい値βを超えているときにも、上記ずれの発生と判断する（ｆａｆ＜−β、またはｆａｆ＞β）。同図２４に示されるように、しきい値βは、しきい値αに比して大きな値が設定されている。そしてこのときＥＣＵ５０は、ずれ発生と判断されている間、上記総ベーパ流量実測値Fvpsを所定率ずつ加増する、若しくは所定率ずつ低減することで、そのずれを補正する。
【０３０１】
そこで同図２４の例では、時刻ｔ０以降の空燃比Ｆ／Ｂ補正値ｆａｆの負方向のずれにより、時刻ｔ１においてその補正項ｆａｆの絶対値がしきい値βを超えると、ＥＣＵ５０は以降、総ベーパ流量実測値Fvpsを所定率ずつ加増していく。こうした態様での実測値Fvpsの加増は、上記実空燃比ＸＯがリッチからリーンに転じて空燃比Ｆ／Ｂ項ｆａｆの増加方向へのスキップ処理がなされる時刻ｔ２まで継続される。
【０３０２】
更にＥＣＵ５０は、そうしたずれ発生の検知から空燃比Ｆ／Ｂ補正値ｆａｆの安定が確認されるまで、ＶＳＶ７１ａの開弁方向への開度変更を中断してその開度をそのまま保持し、総パージ流量Fpgallを一定の状態に保持している。ちなみに本制御装置ではＥＣＵ５０は、上記スキップ処理後の空燃比Ｆ／Ｂ中心値fafav の絶対値がしきい値α以下となることをもって、空燃比Ｆ／Ｂ補正値ｆａｆの安定を確認している。
【０３０３】
この時刻ｔ２でのスキップ処理後、空燃比Ｆ／Ｂ補正値ｆａｆの負方向へのずれによって、その中心値fafav の絶対値が上記しきい値αを超えていれば、ＥＣＵ５０は、総ベーパ流量実測値Fvpsの値を所定値加増する。またこれと共にＥＣＵ５０は、空燃比Ｆ／Ｂ補正値ｆａｆ及びその中心値fafav の値も、上記実測値Fvpsの加増量に相応する分だけ修正する。なお、同図２４の例では、時刻ｔ２に続くスキップ処理の行われた時刻ｔ３においても、同様のずれ発生が検知され、同様の処理がなされている。
【０３０４】
さて、その時刻ｔ３に続くスキップ処理がなされた時刻ｔ４において、空燃比Ｆ／Ｂ補正値ｆａｆの安定が確認されると、ＥＣＵ５０はその時刻ｔ４より、ＶＳＶ７１ａの開弁方向への開度変更を再開し、総パージ流量Fpgallを再び徐増させていく。
【０３０５】
一方同図２４の例では、その時刻ｔ４に引き続いてのスキップ処理がなされた時刻ｔ５においては、今度は空燃比Ｆ／Ｂ補正値ｆａｆの正方向のずれにより、その中心値fafav の絶対値が上記しきい値αを超えている。このときＥＣＵ５０は、総ベーパ流量実測値Fvpsの見積もりが多すぎたものと見なし、総ベーパ流量実測値Fvpsの値を所定値低減するとともに、空燃比Ｆ／Ｂ補正値ｆａｆ及びその中心値fafav の値もその低減量に相応する分だけ修正する。またＥＣＵ５０は、ずれ発生の検知に応じて上記再開したＶＳＶ７１ａの開弁方向への開度変更を再び中断して、その開度を現状のまま保持する。
【０３０６】
以後ＥＣＵ５０は、時刻ｔ６のように空燃比Ｆ／Ｂ補正値ｆａｆの安定が確認されれば、ＶＳＶａの開弁方向への開度変更を再開し、再びずれ発生が検知されればその開弁方向への開度変更を中断するとともに上記実測値Fvps等を修正する。ＥＣＵ５０は、以上例示したような態様で、総パージ流量Fpgallを徐増しながら上記実測値Fvpsを求めている。以上が、初期化完了前の総パージ流量Fpgallの調整（ＶＳＶ開度Dvsvの調整）、及び上記実測値Fvpsの算出にかかるＥＣＵ５０の処理の詳細である。
【０３０７】
以上説明した処理を通じて求められる総ベーパ流量実測値Fvpsによれば、初期化完了前であれ、吸気通路１２へのパージガスのベーパ濃度rvpsを把握することができる（rvps←Fvps／Fpgall）。本制御装置では、そのベーパ濃度rvpsの推移の監視のもとに、上記各物理状態量（タンク発生ベーパ流量Fvptnk、空気層蓄積ベーパ量Mgair 、吸着剤蓄積ベーパ量Mgcan の各値）を初期化するようにしている。
【０３０８】
［２−６−２］  タンク発生ベーパ流量の初期化
続いてその初期化処理の詳細について、図２５及び図２６を併せ参照して説明する。ここではまず、タンク発生ベーパ流量の初期化にかかる処理の詳細を、図２５を参照して説明する。
【０３０９】
さて上述したように、燃料タンク３０の内圧がキャニスタ空気層４５の内圧よりも所定圧以上高く、上記タンク内圧制御バルブ６０が開弁していれば、ベーパライン３５を通じて燃料タンク３０からキャニスタ４０にベーパが流入する（［２−１］、図１４参照）。このとき、ベーパパージ処理が実施されていれば、燃料タンク３０から流入したより高圧なベーパは、キャニスタ空気層４５のパージガス、及び大気導入ライン７２等から導入された空気よりも、優先的にパージライン７１に吸入される。
【０３１０】
よって、ＶＳＶ７１ａを全閉から徐々に開弁していった場合、その開弁開始直後の吸気通路１２へのパージガスは、そのほとんどが上記燃料タンク１０からキャニスタ空気層４５を素通りしてパージラインに直接流入したベーパによって占められる。以下、こうした態様で吸気通路１２に放出されるベーパを、「タンク流入ベーパ」という。またそうしたタンク流入ベーパは、そのほぼ全てがベーパ分によって占められる、すなわちベーパ濃度rvps１００％である、と考えられる。
【０３１１】
そこで上記初期化完了前において、ＶＳＶ７１ａの開弁開始時に、タンク内圧制御バルブ６０が開弁しており、燃料タンク３０からのベーパ流入が許容されていれば、その直後にはタンク流入ベーパのみが吸気通路１２に流入する。なお、こうしたタンク流入ベーパの流量は、上記タンク発生ベーパ流量Fvptnkに対して一定の割合となると考えられる。よってタンク流入ベーパ流量の上限は、タンク発生ベーパ流量Fvptnkに応じてほぼ一義的に決まり、上記割合を定数rvptnkとすれば次の算出式によって求められる（０≦rvptnk≦１）。
<<算出式>>
・  ［タンク流入ベーパ流量］  ←  rvptnk・Fvptnk
なお、定数rvptnkの値は、ベーパパージシステム構成に固有の定数として、試験等によって求めることができる。
【０３１２】
こうした状況は、総パージ流量Fpgallがある程度よりも増大して、パージライン７１へのキャニスタ空気層４５のパージガスの流入、すなわち空気層パージが許容されるようになるまで続く。その間、吸気通路１２に放出されるパージガス中のベーパ濃度rvpsは上記のようにほぼ１００％であるため、総パージ流量Fpgallと上記総ベーパ流量実測値Fvpsとは、ほぼ同一の値となる（Fvps＝Fpgall；Fvps／Fpgall＝１）。
【０３１３】
一方、空気層パージのベーパ濃度rvpairは空気層蓄積ベーパ量Mgair に依存して増減するものの、１００％でないことだけは確かである（［１−３−１］、図９等参照）。よって、総パージ流量Fpgallが空気層パージを許容するレベルを超えると、図２５に示すように、総パージ流量Fpgallの増加率に比して上記実測値Fvpsの増加率が小さくなり、それまでほぼ同一値だった両流量に差が生じる。よって同図２５に示されるように、初期化完了前のパージ開始後に、総パージ流量Fpgallと上記実測値Fvpsとの間に有意な差Δ１が生じたときのそれら流量より、タンク発生ベーパ流量Fvptnkの初期値が求められる。詳しくは、上記有意な差が生じたときの総ベーパ流量実測値Fvpsが上記タンク流入ベーパ流量と同一値であると仮定し、その値をもとに上記算出式から逆算してタンク発生ベーパ流量Fvptnkの初期値を求めている（Fvptnk［初期値］←［タンク流入ベーパ流量］／rvptnk）。なお本制御装置では、ＥＣＵ５０は、上記差（Fpgall−Fvps）が所定値Δ１以上となるときに、かかるタンク発生ベーパ流量Fvptnkの初期化処理を実施している。
【０３１４】
よって本実施形態では、吸気通路１２へのパージ成分のすべてがタンク流入ベーパ流量によって占められると仮定したときの総ベーパ流量Fvpallの理論値（仮定に従えば、この理論値は総パージ流量Fpgallと同一値となる）と、その実測値Fvpsとの対比のもとに、タンク発生ベーパ流量Fvptnkの初期値を求めている。換言すれば、吸気通路１２へのパージガスのベーパ濃度についての上記仮定に基づく理論値（＝１００％）とその実測値（Fvps／Fpgall）との対比により、タンク発生ベーパ流量Fvptnkの初期値が求められている。
【０３１５】
なお、かかる初期化の処理中に上記タンク内圧制御バルブ６０が閉弁していれば、タンク発生ベーパ流量Fvptnkの初期値は勿論「０」となる。ちなみにタンク内圧制御バルブ６０の開閉は、例えば上記タンク内圧センサ３２により検出される燃料タンク３０の内圧の推移によって確認することも可能である。
【０３１６】
［２−６−３］  空気層蓄積ベーパ量の初期化
続いて、空気層蓄積ベーパ量Mgair の初期化にかかる処理の詳細を、図２６を参照して説明する。
【０３１７】
さて、タンク流入ベーパ流量Fvptnkの初期化後、さらに総パージ流量Fpgallを徐増していくと、その成分は、タンク流入ベーパと空気層パージとによって占められることとなる。
【０３１８】
空気層パージのベーパ濃度rvpairは空気層蓄積ベーパ量Mgair によって一義的に求まり、同蓄積ベーパ量Mgair が一定であれば、そのベーパ濃度rvpairも一定に保持される。また空気層パージ流量Fpgairには上限（最大空気層パージ流量Fpgairmx）があり、その値も空気層蓄積ベーパ量Mgair から一義的に求まる（以上、［１−３−１］、図９等参照）。
【０３１９】
よって先の図２５にも示すように、上記初期化完了前のベーパパージ処理において、許容されるタンク流入ベーパをパージしきれるレベル以上に総パージ流量Fpgallが増加した後の総ベーパ流量実測値Fvpsの増加率は一定となる。またそのときの総ベーパ流量実測値Fvpsの増加率は、上記空気層蓄積ベーパ量Mgair より求まる空気層パージのベーパ濃度rvpairに従うと考えられる。
【０３２０】
ＥＣＵ５０は、上記初期化完了前のベーパパージ処理において、総ベーパ流量実測値Fvpsが更新される都度、その更新値に応じて、空気層パージのベーパ濃度rvpairの仮値rvpsを求めている。この空気層パージのベーパ濃度仮値rvpsの値は、吸気通路１２へのパージ成分がタンク流入ベーパ及び空気層パージのみによって占められている限りは、ほとんど不変であると推定される。ＥＣＵ５０は、そのベーパ濃度仮値rvpsを下記算出式に従って求めている。
<<算出式>>
・  rvps  ←  （Fvps−rvptnk・Fvptnk）／（Fpgall−rvptnk・Fvptnk）
更にＥＣＵ５０は、このベーパ濃度仮値rvpsに対する総ベーパ流量の推定値Fvptを求めている（Fvpt←rvps・Fpgall）。この推定値Fvptは、吸気通路１２へのパージ成分がタンク流入ベーパ及び空気層パージのみで占められていると仮定したときの、総ベーパ流量Fvpallの理論値に相当する。
【０３２１】
さてその後、更に総ベーパ流量Fpgallが増加され、許容される空気層パージのすべてをパージしきれるようになると、即ち空気層パージ流量Fpgairがその最大値Fpgairmxに達すると、吸気通路１２へのパージ成分に吸着剤脱離パージが更に加わるようになる。そしてその結果、吸気通路１２へのパージガスのベーパ濃度が変化し、図２６に示されるように、総ベーパ流量実測値Fvpsの増加傾向が変化して、上記総ベーパ流量Fvpallについての実測値Fvpsと理論値Fvptとの間に有意な差が生じる。よってこれにより、空気層パージ成分の総量を把握でき、それをもとに上記各蓄積ベーパ量Mgair 、Mgcan の初期値を求めることができる。本制御装置ではＥＣＵ５０は、上記実測値Fvpsと理論値Fvptとの差が所定値Δ２となったとき、上記各蓄積ベーパ量Mgair 、Mgcan の初期値算出にかかる初期化処理を実施する。
【０３２２】
上記実測値Fvpsと理論値Fvptとの間に有意な差が認められ、吸着剤脱離パージの合流が確認されれば、そのときの総パージ流量Fpgall及び上記初期化されたタンク流入ベーパ流量（rvptnk・Fvptnk）から、空気層パージ流量Fpgairの最大値、即ち最大空気層パージ流量Fpgairmxを求めることができる（Fpgairmx←Fpgall−rvptnk・Fvptnk）。また上記総ベーパ流量実測値Fvps及びタンク流入ベーパ流量（rvptnk・Fvptnk）からは、最大空気層ベーパ流量Fvpairmxを求めることができる（Fvpairmx←Fvps−rvptnk・Fvptnk）。更には上記のように、それら最大空気層パージ流量Fpgairmx及び最大空気層ベーパ流量Fvpairmxは、空気層蓄積ベーパ量Mgair より一義的に求められる。
【０３２３】
よってそれらの相関関係に基づき、上述の上記両最大流量Fpgairmx、Fvpairmxの演算ロジックを逆算すれば、空気層蓄積ベーパ量Mgair の初期値が求められる。
【０３２４】
一方、本制御装置ではＥＣＵ５０は、総ベーパ流量実測値Fvpsより、総パージ流量Fpgallの推定最大値tFpgmxを求めている。この推定最大値tFpgmxは、吸気通路１２へのパージ成分がタンク流入ベーパと空気層パージのみであるとの仮定に基づいたときの総パージ流量Fpgallの最大値の理論値である。この推定最大値tFpgmxの値は、以下の態様にて求められる。
【０３２５】
上記仮定の成立するときの空気層ベーパ流量Fvpairの値は、上記実測値Fvpsからタンク流入ベーパ流量を減算した値となる。一方、同一の空気層蓄積ベーパ量Mgair での最大空気層ベーパ流量Fvpairmxと最大空気層パージ流量Fpgairmxとの対応関係（図１９、図２１等参照）より把握されるように、空気層パージのベーパ濃度rvpairには上限がある。そこで空気層パージのベーパ濃度rvpairを最大限に見積もったときの上記実測値Fvpsから推定される総パージ流量Fpgallの最大値が、上記総パージ流量Fpgallの推定最大値tFpgmxとなる。そこで上記ベーパ濃度rvpairの最大値を定数RVPAIRMXとおけば、その推定最大値tFpgmxを下記の算出式より求めることができる。
<<算出式>>
・  tFpgmx  ←  rvptnk・Fvptnk＋RVPAIRMX・（Fvps−rvptnk・Fvptnk）
ちなみに本制御装置では、ＥＣＵ５０は、予めそのメモリ内に記憶された上記実測値Fvpsと推定最大値tFpgmxとの対応関係を示した演算マップ（図示略）を用いて、その推定最大値tFpgmxを求めている。
【０３２６】
なお本パージシステムでは、吸着剤だつパージのベーパ濃度rvpcanは通常、空気層パージのベーパ濃度rvpairに比して小さい値となる。このため、吸気通路１２へのパージ成分に吸着剤脱離パージが合流すれば、同図２６に示されるように、総ベーパ流量実測値Fvpsの増加率は減少する傾向にある。その結果、吸着剤脱離パージの流量（吸着剤内空気流量Fpgcan）が増加するほどに、総パージ流量Fpgallと上記推定最大値tFpgmxとの差が開くこととなる。
【０３２７】
したがって総パージ流量Fpgallと上記推定最大値ｔFpgmxとの差が、空気層蓄積ベーパ量Mgair の変化に対する空気層パージのベーパ濃度rvpairの変化量に比して十分に大きな値となれば、上記吸着剤脱離パージの合流を確認できる。そこで本制御装置では、ＥＣＵ５０は、総ベーパ流量Fpgallとその推定最大値tFpgmxとの差が所定値以上となったときにも、そのときの総ベーパ流量実測値Fvpsの値に基づいて空気層蓄積ベーパ量Mgair の初期化を行っている。
【０３２８】
［２−６−４］  吸着剤蓄積ベーパ量の初期化
続いて、残る吸着剤蓄積ベーパ量Mgcan の初期化にかかる処理の詳細を同図２６を参照して説明する。
【０３２９】
吸着剤蓄積ベーパ量Mgcan は、吸着剤脱離パージのベーパ濃度rvpcanから一義的に求められる（［１−３−３］、図２２等参照）。そこで上記空気層蓄積ベーパ量Mgair の初期化完了後も総パージ流量Fpgallの徐増を続け、総ベーパ流量実測値Fvpsの増加率から上記ベーパ濃度rvpcanを求めれば、吸着剤蓄積ベーパ量Mgcan の初期値を求めることができる。
【０３３０】
ただし上記態様で吸着剤蓄積ベーパ量Mgcan の初期化を行えば、空気層蓄積ベーパ量Mgair の初期化完了後もしばらくは、上記初期化完了前のベーパパージ処理を継続する必要があり、上記物理モデルに基づくベーパパージ処理への移行が遅れることとなる。そこで本制御装置では、下記態様にて初期値を求めることで、上記空気層蓄積ベーパ量Mgair の初期化と同時に吸着剤蓄積ベーパ量Mgcan の初期化を行っている。
【０３３１】
こうした初期化のための上記初期化完了前のベーパパージ処理（初期化処理）の開始前には、即ちエンジン始動後の始めてのベーパパージ処理の開始前には、パージシステム２０は長期に亘り、定常状態に保持されている。このため、その間にキャニスタ４０内は平衡状態となっており、上記ベーパ吸着速度Fvpatcと上記自然脱離速度Fvpctaとが釣り合っているものと考えられる（Fvpatc＝Fvpcta）。そこで本制御装置では、上記初期化処理の開始時のキャニスタ４０内が平衡状態にあると仮定して、上記求められた空気層蓄積ベーパ量Mgair の初期値より、吸着剤蓄積ベーパ量Mgcan の初期値を求めている。
【０３３２】
上記［１−３−３］で説明したように、ベーパ吸着速度Fvpatc及び自然脱離速度Fvpctaは、各々下記の算出式で求められる。
・  Fvpatc  ←  k1・Mgair ・（VPCANMX−Mgcan ）
・  Fvpcta  ←  k2・Mgcan
よってそれら速度の釣り合う平衡状態での吸着剤蓄積ベーパ量Mgcan は、下記算出式により求められる。
・  Mgcan ←  k1・VPCANMX ・Mgair／（k1・Mgair ＋k2）
これにより、空気層蓄積ベーパ量Mgair の初期化と同時に吸着剤蓄積ベーパ量Mgcan の初期化も完了し、直ちに上記物理モデルに基づくベーパパージ処理に移行することができる。なお本制御装置を、上記のような総ベーパ流量実測値Fvpsの増加率による吸着剤蓄積ベーパ量Mgcan の初期化を行うよう変更して実施することも勿論可能である。
【０３３３】
［２−７］  各物理状態量の修正処理（図１５のＳ６００）
次に、本制御装置での上記各物理状態量の修正処理の詳細について、図２７〜図３４を併せ参照して説明する。
【０３３４】
本制御装置では、上記定期更新処理（［２−３］参照）を通じて、両蓄積ベーパ量Mgair 、Mgcan の値を適正に保持するようにしているが、それでもそれらの値に誤差が生じることがある。たとえそれらの値の見積もりの誤りがほんの僅かでもあっても、上記定期更新処理に際しての各蓄積ベーパ量Mgair 、Mgcan の値の更新量にエラーが生じることとなる。そして定期更新を重ねる毎にそれら両蓄積ベーパ量Mgair 、Mgcan の値にエラーが積み重なり、やがてはそれらの値に大きなずれを生じさせる。こうしたずれは、それらの値に基づき推定される総ベーパ流量Fvpallの、そしてひいてはパージ補正値ｆｐｇの見積もりを誤らせることとなる。
【０３３５】
一方、空燃比Ｆ／Ｂ制御では、空燃比Ｆ／Ｂ補正値ｆａｆの中心値fafav が「０」近傍に保持されるように空燃比学習値ＫＧが設定されている。なお以下では、特に必要のない限り、その中心値fafav を単に「空燃比Ｆ／Ｂ中心」という。またベーパパージの実施中であれ、吸気通路１２に放出されるベーパパージの影響分を上記パージ補正値ｆｐｇにより吸収することで、表面上はベーパパージの影響が無いかの如く、空燃比Ｆ／Ｂ制御が継続されるようになっている。よって、上記パージ補正値ｆｐｇの見積もりにエラーがあれば、ベーパパージの実施に伴い、上記空燃比Ｆ／Ｂ中心が「０」近傍からずれてしまうようになる（［２−４］参照）。
【０３３６】
よって本制御装置では、パージ実施中の空燃比Ｆ／Ｂ中心の推移を監視し、そのずれの検知に応じて上記各蓄積ベーパ量Mgair 、Mgcan の値を修正する処理を行っている。ちなみに、かかる修正処理について本制御装置では厳密には、そうした空燃比Ｆ／Ｂ中心のずれを、上記中心値fafav の徐変値fafsm に基づき検知している。
【０３３７】
図２７は、こうした修正処理にかかる「物理状態量の修正ルーチン」を示している。本ルーチンの処理は、ＥＣＵ５０によって、上記パージ補正値の算出処理（［２−４］参照）に引き続き実行される。以下、同図２７を併せ参照しつつ、本制御装置での修正処理の詳細な内容を説明する。
【０３３８】
同図２７に示されるように、ＥＣＵ５０は、両蓄積ベーパ量Mgair 、Mgcan の値のうち、空燃比Ｆ／Ｂ中心のずれ態様に応じて修正の必要な値を選択し（同図２７のＳ６１０〜Ｓ６３０）、その選択した値を修正するようにしている。
【０３３９】
［２−７−１］  空燃比Ｆ／Ｂ中心のずれの要因判定（Ｓ６１０〜Ｓ６３０）空気層ベーパ流量Fvpairの算出にかかる空気層蓄積ベーパ量Mgair の値にエラーがあるときと、吸着剤脱離ベーパ流量Fvpcanの算出にかかる吸着剤蓄積ベーパ量Mgcan の値にエラーがあるときとでは、上記空燃比Ｆ／Ｂ中心のずれ態様も自ずと異なったものとなる。
【０３４０】
パージ実施中の吸気通路内圧ＰＭ等のエンジン運転状況の変化に伴う空気層パージ状況の変化に応じて、空気層蓄積ベーパ量Mgair の値は大きく変動する。また燃料タンク３０でのベーパ発生状況、すなわちタンク発生ベーパ流量Fvptnkの変化によっても、その値Mgair は大きく変動する。更に空気層パージの成分は、キャニスタ空気層４５のパージガスそのものであるため、空気層蓄積ベーパ量Mgair のエラーは、空気層ベーパ流量Fvpairの見積もりに鋭敏に反映される。よって空気層蓄積ベーパ量Mgair にエラーがあれば、ベーパパージ実施中の空燃比Ｆ／Ｂ中心に大きいずれが発生する。
【０３４１】
一方、キャニスタ４０の吸着剤中に蓄積されたベーパ量（吸着剤蓄積ベーパ量Mgcan ）の変化は比較的緩やかである。また吸着剤蓄積ベーパ量Mgcan のエラーは、吸着剤脱離パージのベーパ濃度rvpcanの誤差として反映されるのみで、吸着剤脱離ベーパ流量Fvpcanの値自体への影響は比較的小さいものとなる。よって吸着剤蓄積ベーパ量Mgcan にエラーのある場合には、空燃比Ｆ／Ｂ中心（fafsm[can]）のずれは、図２８に例示されるように総パージ流量Fpgallの変化に対応した態様で生じることとなる。
【０３４２】
なお、本制御装置では、修正に際して用いられる空燃比Ｆ／Ｂ中心の指標値として、空気層蓄積ベーパ量Mgair の修正用と吸着剤蓄積ベーパ量Mgcan の修正用とで、それぞれ別の空燃比Ｆ／Ｂ中心値fafav の徐変値fafsm を用いている。空気層蓄積ベーパ量Mgair の修正用の徐変値fafsm[air]は、吸着剤蓄積ベーパ量Mgcan の修正用の徐変値fafsm[can]よりも、空燃比Ｆ／Ｂ中心値fafav の変化への追従性が大きくなるように設定されている。こうした空燃比Ｆ／Ｂ中心値fafav の変化への追従性の度合いは、各徐変値fafsm[air]、fafsm[can]の徐変率、値の更新周期等の調整により、適宜設定することができる。
【０３４３】
こうして本制御装置では、上記各値Mgair 、Mgcan のエラーが上記空燃比Ｆ／Ｂ中心のずれに与える影響の傾向に応じて、追従性の異なる徐変値を使い分けて修正処理を行うようにしている。このため、空燃比Ｆ／Ｂ中心のずれに応じた要修正値の判定やその修正量の設定等をより的確に行える。
【０３４４】
より詳しくは、同図２７にも示されるように、ＥＣＵ５０は、下記態様で要修正値のエラー判定を行っている。
（吸着剤蓄積ベーパ量Mgcan のエラー判定）
ＥＣＵ５０は、下記エラー条件（ａ）（ｂ）のいずれかの成立をもって、総ベーパ流量Fpgallの変化に応じた空燃比Ｆ／Ｂ中心のずれの検知有りと判定する（Ｓ６１０：ＹＥＳ）。そして所定の修正条件（［２−７−２］参照）さえ成立していれば（Ｓ６８０：ＹＥＳ）、吸着剤蓄積ベーパ量Mgcan の値を修正する（Ｓ６９０）。
【０３４５】
（ａ）  空燃比Ｆ／Ｂ中心の安定時とずれ時との総パージ流量Fpgallの差（若しくは吸着剤内空気流量Fpgcanを用いても良い）が一定値以上である。なお本制御装置では、下記判定式(a1)の成立をもって空燃比Ｆ／Ｂ中心の安定時と判定し、下記判定式（a2）の成立をもってそのずれ時と判定している。
<<判定式>>
・  ｜fafsm[can]｜＜SFFAFSMCAN      …（a1）
ここで「SFFAFSMCAN」はfafsm[can]用の安定判定値であり、上式（a1）の成立時には空燃比Ｆ／Ｂ中心が「０」近傍に留まっているように、所定の定数としてその値が設定されている。
・  ｜fafsm[can]｜＞ERFAFSMCAN      …（a2）
ここで「ERFAFSMCAN」はfafsm[can]用のずれ判定値であり、試験等の結果に基づき設定された所定の定数である（SFFAFSMCAN＜ERFAFSMCAN）。
【０３４６】
（ｂ）  総パージ流量Fpgall（若しくは吸着剤内空気流量Fpgcanを用いても良い）の流量変化が所定期間よりも長く継続しており、且つその流量変化に応じた側への噴射補正絶対量のずれも、その所定期間継続している。
【０３４７】
図２８は、こうした総パージ流量Fpgallの変化に伴う空燃比Ｆ／Ｂ中心のずれ発生時の各パラメータの推移を示している。同図２８において、時刻ｔ１以前は、空燃比Ｆ／Ｂ中心は安定した状態となっており（上式（a1）成立）、時刻ｔ２には空燃比Ｆ／Ｂ中心のずれ有りと判定されている（上式（a2）成立）。
【０３４８】
またＥＣＵ５０は、上記エラー判定条件（ａ）（ｂ）が共に不成立であれ（Ｓ６１０：ＮＯ）、下記判定式（c1）（c2）の双方が共に成立することをもって、あまり大きくはないものの無視し得ない空燃比Ｆ／Ｂ中心のずれ有りと判定する（Ｓ６３０：ＹＥＳ）。このときにも、やはり上記所定の修正条件（［２−７−２］参照）さえ成立していれば（Ｓ６８０：ＹＥＳ）、吸着剤蓄積ベーパ量Mgcan の値を修正する（Ｓ６９０）。
<<判定式>>
・  ｜fafsm[can]｜＞ERFAFSMCAN      …（c1）
・  ｜fafsm[air]｜≦ERFAFSMAIR      …（c2）
ここで「ERFAFSMAIR」はfafsm[air]用のずれ判定値である。なお上式（c2）の不成立時には、空燃比Ｆ／Ｂ中心の大きいずれ有りと判定される。
【０３４９】
（空気層蓄積ベーパ量Mgair のエラー判定）
ＥＣＵ５０は、下記判定式（d1）（d2）のいずれかの成立をもって、空燃比Ｆ／Ｂ中心の大きいずれ有りと判定する（Ｓ６２０：ＹＥＳ）。
<<判定式>>
・  ｜fafsm[air]｜＞ERFAFSMAIR      …（d1）
・  ｜ｆａｆ｜＞ERFAFAIR            …（d2）
ここで上記判定式（d1）が成立し、所定の修正条件（［２−７−３］参照）さえ成立すれば（Ｓ６４０：ＹＥＳ）、ＥＣＵ５０は空気層蓄積ベーパ量Mgair の値を修正する（Ｓ６５０）。
【０３５０】
なお、上記各ずれ判定値ERFAFAIR、ERFAFSMAIR、ERFAFSMCANは、それぞれ空燃比Ｆ／Ｂ中心のずれが許容できないレベルに達したときの空燃比Ｆ／Ｂ補正値ｆａｆ、及びその中心値fafav の各徐変値fafsm[air]、fafsm[can]の値に相当する所定の定数として設定されている（ERFAFAIR＞ERFAFSMAIR＞ERFAFSMCAN）。
【０３５１】
［２−７−２］  吸着剤蓄積ベーパ量Mgcan の修正処理（図２４のＳ６８０、Ｓ６９０）
続いて、上記判定処理の結果に応じて実施される吸着剤蓄積ベーパ量Mgcan の修正処理の詳細な内容を説明する。
【０３５２】
さて図２７に示されるように、下記のいずれかのとき、ＥＣＵ５０は吸着剤蓄積ベーパ量Mgcan の修正条件の成立の有無を判定する（Ｓ６８０）。
・  上記判定処理により吸着剤蓄積ベーパ量Mgcan が要修正と判定されたとき（Ｓ６１０：ＹＥＳ、またはＳ６３０：ＹＥＳ）。
・  空気層蓄積ベーパMgair の修正要求がなされるも、その修正条件が満たされていないとき（Ｓ６４０：ＮＯ）。
ここでの吸着剤蓄積ベーパ量Mgcan の修正条件は、同量Mgcan の不適切な修正がなされないよう設定されている。以下に、その修正条件が不成立となる条件の一部を示す。
（１）  吸着剤内空気流量Fpgcanが所定量よりも少ないとき。
（２）  吸着剤蓄積ベーパ量Mgcan の現状値が、その設定許容範囲の上下限に既に到達しており、且つその設定許容範囲から外れる側に対して同蓄積ベーパ量Mgcan の修正要求がなされたとき、など。
【０３５３】
なお上記（１）の条件が満たされるときには、空燃比Ｆ／Ｂに影響を与えるほどの吸着剤脱離パージは実際には行われておらず、空気層蓄積ベーパ量Mgair 側にエラーがあると推定される。
【０３５４】
一方、上記（２）の条件での吸着剤蓄積ベーパ量Mgcan の設定許容範囲は、下記の２つのガードにより規定されている。これにより、如何に修正要求がなされようとも、上記（２）の修正不成立条件により、その設定許容範囲からの吸着剤蓄積ベーパ量Mgcan の逸脱が禁止されている。
・  絶対値ガード  ：  吸着剤蓄積ベーパ量Mgcan の設定許容範囲は、物理的に吸着剤に吸着可能なベーパの量により規定されている。即ち、同吸着剤蓄積ベーパ量Mgcan の値の設定許容範囲は「０」以上で、且つ吸着剤に吸着の許容されるベーパの最大量である上記最大吸着量VPCANMX 以下となる。
・  相対値ガード  ：  エンジン始動後に適宜にベーパパージ処理を実施して、吸着されたベーパを十分に脱離しておけば、燃料タンク３０から多量のベーパが流入しようとも、キャニスタ空気層４５が緩衝となり、吸着剤のベーパ吸着量、すなわち吸着剤蓄積ベーパ量Mgcan の急増は抑制される。このため、エンジン始動後に十分な脱離がなされた後に、吸着剤脱離パージのベーパ濃度rvpcanが急増することは理論上は有り得ても、実際にそれが生じる可能性はほぼ皆無である。よって、エンジン始動後における上記ベーパ濃度rvpcanの最小値を記憶保持しておき、上記物理モデルに従い推定されるベーパ濃度rvpcanがその最小値に所定値を加算した値を超えないように、吸着剤蓄積ベーパ量Mgcan の上限値を規定する。なお、例えばキャニスタ４０内が十分に暖機されている、空燃比Ｆ／Ｂが安定した状態にある、総パージ流量Fpgallが所定値以上であるなど、十分な脱離がなされた後で、且つ上記ベーパ濃度rvpcanの推定値の信頼性が十分な状況に限り、上記最小値を記憶保持することが望ましい。
【０３５５】
ＥＣＵ５０は、上記修正条件が成立していれば（Ｓ６８０：ＹＥＳ）、次の態様で吸着剤蓄積ベーパ量Mgcan の修正処理を実施する（Ｓ６９０）。
このときの空燃比Ｆ／Ｂ中心（fafsm[can]）のずれは、吸着剤脱離ベーパ流量Fvpcanの見積もりの誤りによるパージ補正値ｆｐｇの誤差が原因であると推定できる。一方、吸着剤脱離ベーパFvpcanは、吸着剤脱離パージのベーパ濃度rvpcanと吸着剤内空気流量Fpgcanとの積として求められ、そのベーパ濃度rvpcanは、吸着剤蓄積ベーパ量Mgcan により一義的に求められる（［１−３−３］、［２−４−２］、図２２等参照）。よって、上記空燃比Ｆ／Ｂ中心のずれに相当する吸着剤脱離ベーパ流量Fvpcanの誤差を求め、その誤差より把握される上記ベーパ濃度rvpcanの推定誤差に応じて吸着剤蓄積ベーパ量Mgcan の修正量を求めることができる（図２８参照）。
【０３５６】
本制御装置では、ＥＣＵ５０は、下記算出式に順次従って吸着剤蓄積ベーパ量Mgcan の修正処理を行っている。
<<算出式>>
・  △rvpcan  ←  fafsm[can]／Fpgcan
・  rvpcan［修正値］  ←  rvpcan［現状値］十△rvpcan
・  △Mgcan ←  fnc.｛rvpcan［修正値］｝
・  Mgcan ［修正値］  ←  Mgcan ［現状値］十△Mgcan
なお「△rvpcan」はベーパ濃度rvpcanの推定誤差を、「△Mgcan 」は吸着剤蓄積ベーパ量Mgcan の修正量をそれぞれ示している。また関数fnc.｛rvpcan［修正値］｝は、吸着剤脱離ベーパ流量Fvpcanの算出処理にかかる上記ベーパ濃度rvpcanの演算ロジックの逆算関数であり、図２２に例示するような演算マップにより示される吸着剤蓄積ベーパ量Mgcan と上記ベーパ濃度rvpcanとの相関関係に基づき求められている。
【０３５７】
一方ＥＣＵ５０は、上記修正条件が不成立であれば（Ｓ６８０：ＮＯ）、空気層蓄積ベーパ量Mgair 側の修正処理（Ｓ６４０）に移行する。即ち、吸着剤蓄積ベーパ量Mgcan の修正に不適切な状況に有れば、たとえ同値Mgcan の修正要求がなされていようとも、空気層蓄積ベーパ量Mgair の値を修正することで、現状のパージ補正値ｆｐｇの不適合の解消を図るようにしている。
【０３５８】
［２−７−３］  空気層蓄積ベーパ量Mgair の修正処理（Ｓ６４０、Ｓ６５０）
更に上記判定処理の結果に応じて実施される空気層蓄積ベーパ量Mgair の修正処理の内容の詳細を、図２９〜図３２を併せ参照して説明する。
【０３５９】
図２７に示されるように、下記のいずれかの場合、ＥＣＵ５０は、空気層蓄積ベーパ量Mgair の修正条件の成立の有無を判定する（Ｓ６４０）。
・  上記判定処理により空気層蓄積ベーパ量Mgair が要修正と判定されたとき（Ｓ６２０：ＹＥＳ）。
・  吸着剤蓄積ベーパ量Mgcan の修正要求がなされるも、その修正条件が満たされていないとき（Ｓ６８０：ＮＯ）。
ここでの空気層蓄積ベーパ量Mgair の修正条件は、同量Mgair の不適切な修正がなされないよう設定されている。以下に、その修正条件が不成立となる条件の一部を示す。
（１）  上記燃料噴射量Ｑｆｉｎを減量させる側への空燃比Ｆ／Ｂ中心のずれが検知されており、且つそのずれ分を現状の空気層ベーパ流量Fvpairに加算した値（修正後の空気層ベーパ流量Fvpair）が、現状の最大空気層ベーパ流量Fvpairmxに達していないとき。
（２）  上記燃料噴射量Ｑｆｉｎを減量させる側への空燃比Ｆ／Ｂ中心のずれが検知されており、且つ現状の総パージ流量Fpgallが、最大空気層パージ流量Fpgairmxの仮定値tFpgairmx に達していないとき（図３２の時刻ｔ２）。この仮定値tFpgairmx は、現状の空気層ベーパ流量Fvpairが最大流量、即ち最大空気層ベーパ流量Fvpairmxであると想定したときの最大空気層パージ流量Fpgairmxの理論値である。
（３）  空気層蓄積ベーパ量Mgair の現状値が、その設定許容範囲の上下限に既に到達しており、且つその設定許容範囲から外れる側に対して同蓄積ベーパ量Mgair の修正要求がなされたとき、など。
【０３６０】
上記（３）の不成立条件における空気層蓄積ベーパ量Mgair の設定許容範囲は、物理的にキャニスタ空気層４５に存在可能なベーパの量により規定されている。即ち、同空気層蓄積ベーパ量Mgair の設定許容範囲は「０」以上で、且つキャニスタ空気層４５に存在の許容されるベーパの上限量である空気層飽和ベーパ量VPAIRMX 以下となる。なお空気層飽和ベーパ量VPAIRMX は、キャニスタ空気層４５の容積（同空気層４５に存在する空気の体積）に応じて決まる定数として求められる。
【０３６１】
ＥＣＵ５０は、上記修正条件が成立していれば（Ｓ６４０：ＹＥＳ）、次の態様で空気層蓄積ベーパ量Mgair の修正処理を実施する。
上記判定式（d2）の成立する空燃比Ｆ／Ｂ補正値ｆａｆ自体の大きな逸脱が検知された場合（ｌｆａｆｌ＞ERFAFAIR）、ＥＣＵ５０は、図２９に示されるように、かかる補正項ｆａｆの逸脱が検知されている間、最大空気層ベーパ流量Fvpairmxを所定率で修正している。即ち、その期間、最大空気層ベーパ流量Fvpairmxを、予め定められた周期毎に所定値ずつ修正するようにしている。そして、図２１に例示するような演算マップにより示される相関関係に基づいて、その修正された最大空気層ベーパ流量Fvpairmxに応じて空気層蓄積ベーパ量Mgair を修正する。
【０３６２】
一方、上記判定式（d1）の成立する空燃比Ｆ／Ｂ中心の大きな逸脱が検知された場合（｜fafsm[air]｜＞ERFAFSMAIR）、ＥＣＵ５０は、図３０〜図３２に例示されるような態様で、空気層蓄積ベーパ量Mgair の修正処理を実施する。
【０３６３】
このときの空燃比Ｆ／Ｂ中心（fafsm[air]）のずれは、空気層ベーパ流量Fvpairの見積もりの誤りによるパージ補正値ｆｐｇの誤差が原因であると推定できる。そして上記のように現状の空気層ベーパ流量Fvpairに空燃比Ｆ／Ｂ中心のずれ分を加算した値が現状の最大空気層ベーパ流量Fvpairmxを上回つていれば、現状の最大空気層ベーパ流量Fvpairmxに見積もりの誤りがあると推定できる。
【０３６４】
よってこのとき本制御装置では、ＥＣＵ５０は、下記算出式に順次従って空気層蓄積ベーパ量Mgair の修正処理を行っている。
<<算出式>>
・  Fvpairmx［修正値］  ←  Fvpair［現状値］＋fafsm[air]
・  Mgair ［修正値］  ←  fnc.｛Fvpairmx［修正値］｝
ここで関数fnc.｛Fvpairmx［修正値］｝は、上述の最大空気層ベーパ流量Fvpairmxの演算ロジック（［２−４−２］等参照）の逆算関数であり、図２１に例示するような演算マップにより示される空気層蓄積ベーパ量Mgair と最大空気層ベーパ流量Fvpairmxとの相関関係に基づき求められている。図３０の時刻ｔ１や時刻ｔ３、図３２の時刻ｔ１や時刻ｔ３においては、以上のようにして空気層蓄積ベーパ量Mgair の修正処理がなされている。
【０３６５】
（空気層蓄積ベーパ量Mgair の修正条件不成立時の処理）
なお本制御装置では、上記（１）の条件により修正条件が不成立となったときには、ＥＣＵ５０は下記の処理を実施する。
【０３６６】
図３０の時刻ｔ２のように、上記（１）の条件が成立するときには、許容される空気層パージの全てが吸気通路１２にパージされておらず、そのときの上記空燃比Ｆ／Ｂ中心のずれ分は、最大空気層ベーパ流量Fvpairmxに要する修正分の一部を示すのみである。そのため、かかる状況下では必要な最大空気層ベーパ流量Fvpairmxの修正量を適正に求めることができず、よって空気層蓄積ベーパ量Mgair の適正な修正も不能である。
【０３６７】
そこで上記（１）の条件成立時においてＥＣＵ５０は、空気層蓄積ベーパ量Mgair の修正処理は取り敢えずは禁止する。そして図３１に例示のように、現状の値のまま保持されたその空気層蓄積ベーパ量Mgair に応じた空気層ベーパ流量Fvpairの算出値に、上記空燃比Ｆ／Ｂ中心のずれ分を加算した値を空気層ベーパ流量Fvpairの仮値とするようにしている。そしてこれにより、取り敢えずは、現状のパージ補正値ｆｐｇの不整合だけは解消するようにしている。
【０３６８】
なお同図３０の時刻ｔ３のように、修正後の空気層ベーパ流量Fvpairの値、即ち空燃比Ｆ／Ｂ中心のずれ分を補正した後の空気層ベーパ流量Fvpairの値により見積もられる最大空気層ベーパ流量Fvpairmxを、現状の総パージ流量Fpgallが下回っている。このときにも、やはり許容される空気層パージの全てが吸気通路１２にパージされておらず、本来あるべき最大空気層パージ流量Fpgairmxを適正に把握すること、即ち空気層蓄積ベーパ量Mgair の修正値を厳密に求めることはできない。ただし、このときの本来あるべき最大空気層ベーパ流量Fvpairmxが少なくとも、上記修正後の空気層ベーパ流量Fvpairの値により推定されるその理論値以上であることだけは確かである。よってこの場合には、上記算出式に従い空気層蓄積ベーパ量Mgair の修正処理を行い、予測可能な範囲で同量Mgair を修正する。
【０３６９】
また本制御装置では、上記（２）の条件により修正条件が不成立となったときには、ＥＣＵ５０は下記の処理を実施する。
図３２の時刻ｔ２のように、上記（２）の条件が成立するときには、やはり許容される空気層パージのすべてが吸気通路１２にパージされていないため、最大空気層ベーパ流量Fvpairmxの要修正量を適正に見積ることはできない。そこでＥＣＵ５０は、この場合にも、空気層蓄積ベーパ量Mgair の値は現状のまま保持して、空気層ベーパ流量Fvpairを空燃比Ｆ／Ｂ中心のずれ分補正するに留め、現状のパージ補正値ｆｐｇの不整合だけは解消するようにしている。
【０３７０】
なお上記（３）の条件により修正条件が不成立となったときには、ＥＣＵ５０は、吸着剤蓄積ベーパ量Mgcan 側の修正処理（Ｓ６８０）に移行する。即ち、その設定許容範囲の逸脱を規制するため空気層蓄積ベーパ量Mgair の修正が不能であれば、たとえ同量Mgair の修正要求がなされていようとも、吸着剤蓄積ベーパ量Mgcan の値を修正することで、現状のパージ補正値ｆｐｇの不適合の解消を図るようにしている。
【０３７１】
以上の各場合、適正な空気層蓄積ベーパ量Mgair の修正が許容されるようになった時点で、即ち上記修正後の空気層ベーパ流量Fvpairより推定される最大空気層ベーパ流量Fvpairmxの理論値を総パージ流量Fpgallが上回った時点で、空気層蓄積ベーパ量Mgair は適正な値に修正されることとなる。
【０３７２】
ところで、両蓄積ベーパ量Mgair 、Mgcan が共に設定許容範囲の上下限に達し、いずれの修正も不能となったときには、下記のいずれかの処理によって対応することで、現状のパージ補正値ｆｐｇの不適合の解消することができる。
・  設定許容範囲の上限値VPAIRMX 、VPCANMX に、経時変化や個体差などによる見積もり誤りがあると見なし、それら上限値VPAIRMX 、VPCANMX の少なくとも一方を補正し、修正を許容する。
・  空燃比学習値ＫＧに誤りがあると見なし、その空燃比Ｆ／Ｂ中心のずれに応じて同学習値ＫＧを補正する。
【０３７３】
以上が空気層蓄積ベーパ量Mgair の修正に係る処理内容の詳細である。なお本制御装置では、空気層蓄積ベーパ量Mgair の修正処理に引き続き、ＥＣＵ５０は、タンク発生ベーパ流量Fvptnkの修正処理（Ｓ６６０）、及び吸着剤蓄積ベーパ量Mgcan の反省処理（Ｓ６７０）を併せて実施する。
【０３７４】
［２−７−４］  吸着剤蓄積ベーパ量Mgcan の反省処理（Ｓ６７０）
図３３に示すように、空気層蓄積ベーパ量Mgair の値が修正されると、それに応じて最大空気層パージ流量Fpgairmxの値も変化する（同図３３では空気層蓄積ベーパ量Mgair が減量修正されたときを例示）。ただし修正処理では総パージ流量Fpgallは変化しないため、空気層蓄積ベーパ量Mgair の修正に応じて上記吸着剤内空気流量Fpgcanの値も変化する。即ち、上記修正処理による最大空気層パージ流量Fpgairmxの「増減」分、吸着剤内空気流量Fpgcanが「減増」されることとなる（△Fpgairmx＝−△Fpgcan）。ここで吸着剤蓄積ベーパ量Mgcan の値を現状のまま、即ち空気層蓄積ベーパ量Mgair の修正前の値のまま保持すれば、吸着剤脱離パージのベーパ濃度rvpcanは一定のまま吸着剤内空気流量Fpgcanが変化することとなる。そしてその結果、吸着剤脱離ベーパ流量Fvpcanの推定値も変化してしまい、上記修正処理により適合されるべきパージ補正値ｆｐｇも結局は不適合な値となってしまう。
【０３７５】
よって本制御装置では、修正前の吸着剤蓄積ベーパ量Mgcan が空気層蓄積ベーパ量Mgair のエラー分を吸収していたと見なし、ＥＣＵ５０は上記修正処理に併せて吸着剤蓄積ベーパ量Mgcan の反省処理を実施する（図２７のＳ６７０）。この反省処理とは即ち、上記修正処理前後の吸着剤脱離ベーパ流量Fvpcanの値が一定に保持されるように、吸着剤蓄積ベーパ量Mgcan を修正する処理である。
【０３７６】
この反省処理での吸着剤蓄積ベーパ量Mgcan の修正値（反省処理後の値）は、以下の態様で求められる。即ち、上記空気層蓄積ベーパ量Mgair の修正処理に応じて変更された吸着剤内空気流量Fpgcanの値と上記修正処理前の吸着剤脱離ベーパ流量Fvpcanの値とより、反省処理後の吸着剤脱離パージのベーパ濃度rvpcanをまず求める。そして、図２２に例示されるベーパ濃度rvpcanと吸着剤蓄積ベーパ量Mgcan との相関関係に基づいて、上記求められた反省処理後のベーパ濃度rvpcanより、この反省処理での吸着剤蓄積ベーパ量Mgcan の修正値を算出する。
【０３７７】
［２−７−５］  タンク発生ベーパ流量Fvptnkの修正処理（Ｓ６６０）
ところで、上記修正処理により空気層蓄積ベーパ量Mgair に修正が必要となった背景には、タンク発生ベーパ流量Fvptnkの見積もりの誤りによる定期更新処理時の同蓄積ベーパ量Mgair の更新エラーの積み重なりにその一因があると推定される。よって本制御装置では、ＥＣＵ５０は上記空気層蓄積ベーパ量Mgair の修正処理に併せて、図３４に例示される態様でタンク発生ベーパ流量Fvptnkの修正処理を実施する（図２７のＳ６６０）。
【０３７８】
図３４の時刻ｔ１以前のように、空燃比Ｆ／Ｂ中心が安定している間は、空気層蓄積ベーパ量Mgair は適正な値に保持され、タンク発生ベーパ流量Fvptnkの見積もりも正しいものと考えられる。ちなみに本制御装置では、空気層蓄積ベーパ量Mgair の修正用の空燃比Ｆ／Ｂ中心の徐変値fafsm[air]の絶対値が所定の安定判定値SFFAFSMAIR以下であるとき、上記空燃比Ｆ／Ｂ中心が安定していると判定している。
【０３７９】
一方上記推定によれば、同図３４の時刻ｔ１から時刻ｔ２までの期間のように、空燃比Ｆ／Ｂ中心のずれが生じ始めて（｜fafsm[air]｜＞SFFSFSMAIR）から空気層蓄積ベーパ量Mgair の修正を要するまでの期間には、タンク発生ベーパ流量Fvptnkの見積もりに誤りが生じていると考えられる。そしてその空燃比Ｆ／Ｂ中心のずれは、そのタンク発生ベーパ流量Fvptnkの見積もりの誤りによる空気層蓄積ベーパ量Mgair の更新エラーの積算によるものと考えられる。よって上記空気層蓄積ベーパ量Mgair の修正処理時の空燃比Ｆ／Ｂ中心のずれ量は、即ち同蓄積ベーパ量Mgair の修正量は、上記ずれ発生から修正処理実施までの期間（時刻ｔ１〜時刻ｔ２：時間Ｔ１２）におけるタンク発生ベーパ流量Fvptnkの誤差の積算値と見なすことができる。
【０３８０】
そこで本制御装置では、ＥＣＵ５０は上記空気層蓄積ベーパ量Mgair の修正処理に併せて、以下の態様でタンク発生ベーパ流量Fvptnkの修正処理を実施する。詳しくは、上記修正処理時の空燃比Ｆ／Ｂ中心のずれ量（修正処理時のfafsm[air]の現状値）、即ち修正処理時の空気層蓄積ベーパ量Mgair の修正項△Mgair について、その上記ずれ発生から修正処理実施までの時間（時間Ｔ１２）による時間微分値をタンク発生ベーパ流量Fvptnkの修正項△Fvptnkとし、その修正処理を行う。
<<算出式>>
・  △Fvptnk  ←  △Mgair ／Ｔ１２
・  Fvptnk［修正値］  ←  Fvptnk［現状値］＋△Fvptnk
例えば、空気層蓄積ベーパ量Mgair の修正処理後に、上記算出式に従って演算することで、かかるタンク発生ベーパ流量Fvptnkの修正処理を行うことができる。また上記算出式において、修正項△Mgair に替わりに、空気層蓄積ベーパ量Mgair の修正処理開始時の空燃比Ｆ／Ｂ中心値fafav （より望ましくはその徐変値fafsm[air]）を用いても、同様の上記タンク発生ベーパ流量Fvptnkの修正処理は可能なことは勿論である。
【０３８１】
なお上記本パージシステム２０のように、燃料タンク３０の内圧を検出するタンク内圧センサ３２を備える構成にあっては（図１４参照）、その検出値より燃料タンク３０内でのベーパ発生状況を把握してタンク発生ベーパ流量Fvptnkをある程度推定することは可能である。そこで検出される燃料タンク３０の内圧（タンク内圧）に応じて、上記タンク発生ベーパ流量Fvptnkの設定許容範囲を規定する。そして、たとえ上記修正処理による修正要求が如何になされようとも、その規定された設定許容範囲外への同流量Fvptnkの修正を制限するようにしても良い。例えば、タンク発生ベーパ流量Fvptnkの許容される上限値を、タンク内圧に応じて、そのタンク内圧が高いほど上限値を大きくするように設定する、などの設定許容範囲の規定態様が考えられる。このようにタンク内圧に応じて設定許容範囲を規定することで、検出される実情に反した不適切なタンク発生ベーパ流量Fvptnkの設定を回避することができる。
【０３８２】
［２−８］  ＶＳＶ開度の算出処理（図１５のＳ１００）
本制御装置では、上記物理モデルに基づく総ベーパ流量Fvpallの予測のもと、上記ＶＳＶ７１ａの開度調整による総パージ流量Fpgallの調整を実施することで、パージシステム２０内のベーパ挙動に適合したベーパパージ処理を行つている。これにより、空燃比Ｆ／Ｂ制御に与えるベーパパージ処理の影響が好適に抑制されるようになっている。以下、こうした好適なベーパパージ処理にかかるＶＳＶ開度の算出処理の詳細を、図３５〜図３７を併せ参照して説明する。
【０３８３】
図３６は、そうしたＶＳＶ開度算出にかかる処理ルーチンを示している。本ルーチンの処理は、ベーパパージ処理の実施条件の成立時にＥＣＵ５０によって周期的に実行される。
【０３８４】
本ルーチンでは、ＥＣＵ５０はまず、そのときのエンジン運転状況に応じたＶＳＶ開度（デューティ比）Dvsvの目標値（目標ＶＳＶ開度）tDvsv を求める（Ｓ１１０）。この目標ＶＳＶ開度tDvsv は、例えばエンジン回転速度ＮＥや吸気通路内圧ＰＭ、吸入空気量Ｇａ、エンジン１０及びキャニスタ４０の暖機状態等に応じて、適度な総パージ流量Fpgallを確保するように設定される。
【０３８５】
ただしパージシステム２０内でのベーパ挙動の如何により、総パージ流量Fpgallと総ベーパ流量Fvpallとの相関関係は変化してしまう。よってここで如何様に目標ＶＳＶ開度tDvsv を設定しようとも、それだけでは、空燃比Ｆ／Ｂ制御にベーパパージの与える影響を予測し、その影響を抑制し得るようにＶＳＶ開度Dvsv、及び総パージ流量Fpgallを設定することは困難である。
【０３８６】
そこで本制御装置では、上記物理モデルを用いてＶＳＶ開度Dvsv変更後の総ベーパ流量Fvpallを予測するとともに、下記ガード値により、好適な空燃比Ｆ／Ｂ制御を保証するようにＶＳＶ開度Dvsvを設定する。
【０３８７】
［２−８−１］  ガード値tFvpmxの算出（Ｓ１２０〜Ｓ１２２）
（ａ）  絶対ガード値tFvpmx[AB]の算出（Ｓ１２０）
エンジン１０の吸入空気量Ｇａが少ないときには、たとえ吸気通路１２にパージされるベーパ流量（総ベーパ流量Fvpall）が少量であろうとも、それが空燃比Ｆ／Ｂ制御に与える影響は大となる。そこで吸入空気量Ｇａに応じて許容される総ベーパ流量Fvpallの上限、即ち絶対ガード値tFvpmx[AB]を規定する。ここでは絶対ガード値tFvpmx[AB]は、図３５に例示するように、吸入空気量Ｇａが大きなほど、より多くのベーパの吸気通路１２へのパージを許容するようにその値が設定されている。
【０３８８】
（ｂ）  相対ガード値tFvpmx[RE]の算出（Ｓ１２１）
ＶＳＶ開度Dvsvの変更の結果、総ベーパ流量Fvpallが急変してパージ補正値ｆｐｇの値が一時に大きく変化すると、空燃比Ｆ／Ｂ制御に好ましからぬ影響を与えるおそれがある。特に、パージ補正値ｆｐｇの絶対値が急増する側にＶＳＶ開度Dvsvが変更された場合、パージライン７１でのベーパの輸送遅れの影響が大となり、また上記各物理状態量の見積もりの誤りによるパージ補正値ｆｐｇの誤差も拡大し、空燃比Ｆ／Ｂ制御に悪影響を与える公算がより高くなる。
【０３８９】
そこでここでは、総ベーパ流量Fvpallの現状値に応じて、ＶＳＶ開度Dvsvの変更による総ベーパ流量Fvpallの増大側への変化率を所定値以内とすべく許容される開度変更後の総ベーパ流量Fvpallの上限、即ち相対ガード値tFvpmx[RE]を規定する。ここでの相対ガード値tFvpmx[RE]は、以下の算出式より求められる。
<<算出式>>
・  tFvpmx[RE]  ←  Fvpall［現状値］＋DFVP
ここで「DFVP」は、空燃比Ｆ／Ｂ制御に与える影響を十分に抑制可能な総ベーパ流量Fvpallの増加率の上限値であり、ここでは試験等の結果求められた所定の定数として設定されている。なお、この増加率の上限値DFVPを吸入空気量Ｇａ等に応じて可変設定するようにしても良い。その場合、例えば吸入空気量Ｇａが多いほど増加率の上限値DFVPを大きくするような設定態様が考えられる。また総ベーパ流量Fvpallの値が減少する側についても、同様の相対ガードを設定するようにしても良い。
【０３９０】
（ｃ）  ガード値tFvpmxの算出（Ｓ１２２）
以上によって求められた両ガード値tFvpmx[AB]、tFvpmx[RE]のうちの小さい方の値を最終的なガード値tFvpmxとして設定する。以後、開度変更後の総ベーパ流量Fvpallの予測値が、この最終的なガード値tFvpmxを超えないようにＶＳＶ開度Dvsvの算出が行われる。
【０３９１】
［２−８−２］  ＶＳＶ開度ガード値tDvsvgd の算出（Ｓ１３０〜Ｓ１５０）以上の処理によりガード値tFvpmxを求めた後、ＥＣＵ５０はまず、現状のベーパパージシステム２０のベーパ挙動に則して、総ベーパ流量Fvpallが丁度そのガード値tFvpmxとなるＶＳＶ開度、即ちＶＳＶ開度ガード値tDvsvgd を算出する。このＶＳＶ開度ガード値tDvsvgd の算出処理は、総ベーパ流量Fvpallを上記ガード値tFvpmxとしたときの、上記物理モデルに基づく総ベーパ流量Fvpallの算出ロジックの逆算処理を通じて行われる。その算出処理の詳細は、図３６のステップ１３０〜ステップ１５０に示す通りである。
【０３９２】
上記総ベーパ流量Fvpallをガード値tFvpmxとしたときの吸気通路１２へのパージが空気層パージのみであると予測されるとき、即ちガード値tFvpmxが現状の最大空気層ベーパ流量Fvpairmx未満のときには（Ｓ１３０：ＹＥＳ）、同図３６のステップ１３５に示される算出式によりＶＳＶ開度ガード値tDvsvgd が求められる。
【０３９３】
またそのときの吸気通路１２へのパージが空気層パージと吸着剤脱離パージとの双方が含まれると予測されるときには（Ｓ１４０：ＹＥＳ）、ステップ１４５に示される算出式によりＶＳＶ開度ガード値tDvsvgd が求められる。
【０３９４】
更に、ガード値tFvpmxが現状でパージ可能なベーパ流量の限界を超えている場合には（Ｓ１４０：ＮＯ）、ＶＳＶ開度ガード値tDvsvgd をＶＳＶ開度Dvsvの上限、即ち１００％に設定する（Ｓ１５０）。
【０３９５】
［２−８−３］  ＶＳＶ開度の算出（Ｓ１６０〜Ｓ１８０）
そしてＥＣＵ５０は、こうして求められたＶＳＶ開度ガード値tDvsvgd と目標ＶＳＶ開度tDvsv とを対比する（Ｓ１６０）。そしてＶＳＶ開度ガード値tDvsvgd が目標ＶＳＶ開度tDvsv 未満で有れば（ＹＥＳ）、ガード値tDvsvgd をＶＳＶ開度Dvsvに設定し（Ｓ１７０）、そうでなければ（ＮＯ）目標ＶＳＶ開度tDvsv をそのままＶＳＶ開度Dvsvに設定する（Ｓ１８０）。
【０３９６】
図３７は、以上説明したＶＳＶ開度算出処理に基づく制御態様の例を示している。同図３７では、
（ａ）キャニスタ４０全体に蓄積されたベーパ量が少ないとき、
（ｂ）空気層蓄積ベーパ量Mgair が多く、空気層パージのベーパ濃度rvpairが高いとき、
（ｃ）吸着剤蓄積ベーパ量Mgcan が多く、吸着剤脱離パージのベーパ濃度rvpcanが高いとき、
の３つの状況を例として、ベーパパージ処理開始からのＶＳＶ開度Dvsv、総ベーパ流量Fvpallの推移を示している。
【０３９７】
なお、ここでは総ベーパ流量Fvpallを基本パラメータとして、以上のＶＳＶ開度算出処理を行うようにしているが、パージ補正項ｆｐｇを基本パラメータとして同様のＶＳＶ開度算出処理を行うようにしても良い。ちなみに本制御装置では、パージ補正値ｆｐｇと総ベーパ流量Fvpallとは、符号の正負が異なるのみで一義的な関係に有り、いずれを用いようともその制御結果は同じである。ただし、燃料噴射量Ｑｆｉｎの算出ロジック（［２−５］参照）次第では、それら両パラメータが必ずしも一義的な関係とはならないことがある。その場合には、基本パラメータとしてパージ補正値ｆｐｇを用いることで、空燃比Ｆ／Ｂ制御にベーパパージが与える影響に、より則した態様でかかるＶＳＶ開度算出処理を行うことができる。
【０３９８】
［２−９］  その他の改良点
本発明を具体化したエンジンの空燃比制御装置の一実施形態の詳細は、以上の通りである。続いて本制御装置に加え得る更なる改良点について説明する。
【０３９９】
［２−９−１］  低開度でのＶＳＶ制御
上述したように上記パージシステム２０のＶＳＶ７１ａは、その開度制御にかかる指令値であるＶＳＶ開度（デューティ比）Dvsvと同ＶＳＶ７１ａを通じて吸気通路１２にパージされる総パージ流量Fpgallとが、吸気通路内圧ＰＭ一定の条件下で比例関係（リニアリティ）となるように構成されている。そして上記制御装置では、その関係を用い、吸気通路内圧ＰＭとＶＳＶ開度Dvsvとより、総パージ流量Fpgallを求め、各種処理を実施している（［２−４−１］等参照）。
【０４００】
ところが、ＶＳＶ７１ａの構成部品の寸法公差、温度の影響による寸法変化などのため、図３８に例示するように、ＶＳＶ７１ａは、ある程度よりも小さな開度では、上記比例関係を確保できなくなることがある。以下、こうした比例関係の確保可能なＶＳＶ開度Dvsvの下限値をリニアリティ下限値DVSVL という。ＶＳＶ開度Dvsvがその下限値DVSVL 未満となれば、総パージ流量Fpgallの正確な把握ができなくなり、よって上記パージ補正値の算出処理等、上記物理モデルに基づく各種処理の実施が不能となってしまう。
【０４０１】
こうした場合、かかる低開度へのＶＳＶ開度Dvsvの設定を禁止するのも一つの手ではある。ただし、図３９に示される「低開度でのＶＳＶ制御処理」を実施すれば、かかる状況下であれ、空燃比Ｆ／Ｂ制御に悪影響を与えることなくベーパパージ処理を実施可能となる。
【０４０２】
更に低開度処理中は、正確な総パージ流量Fpgallが不明なことに起因してパージシステム２０でのベーパ挙動を正確には把握できない状態となる。このため、低開度処理中は、上記各物理状態量の更新処理や修正処理（［２−３］、［２−７］参照）についても一時中止して、それら各物理状態量の見積もりの誤りの拡大を防止している。なお低開度処理中に生じた各物理状態量の変化や誤差は、低開度処理終了後の修正処理により補正されることとなる。
【０４０３】
以下、同図３９及び図４０を併せ参照して、その処理の詳細を説明する。なお本ルーチンの処理は、上記ＶＳＶ開度の算出処理（［２−８］、図３６等参照）に引き続きＥＣＵ５０により実行される。そして上記算出処理により求められたＶＳＶ開度ガード値tDvsvgd がリニアリティ下限値DVSVL 未満となると（Ｓ７００：ＹＥＳ）、ＥＣＵ５０は、上述のパージ補正値の算出処理（［２−４］等参照）等の通常のベーパパージ処理を一時中止し、以下の態様で処理を実施する。
【０４０４】
図４０の時刻ｔ０のように、通常の処理から低開度処理に移行したとき（Ｓ７１０：ＮＯ）、ＥＣＵ５０は、ＶＳＶ開度Dvsvを一旦全閉（Dvsv＝「０」％）にするとともに、流入率rvpdtlの値を「０」に設定する（Ｓ７２５）。上記処理の移行は、前回の本ルーチン実施時に低開度での処理が行われたことを示すフラグxDvsvlのオン／オフにより確認されている（Ｓ７０５、Ｓ７６０参照）。
【０４０５】
また流入率rvpdtlは、低開度処理中にのみ用いられる空気層パージのベーパ濃度rvpairの代替値である。流入率rvpdtlは、空燃比Ｆ／Ｂ中心値fafav のずれに基づいて求められている。また低開度処理中は、その流入率rvpdltに応じて、下記算出式に従い総ベーパ流量Fvpallを求めている。
<<算出式>>
・  Fvpall［低開度時］←rvpdtl・Fvpairmx
これにより求められた総ベーパ流量Fvpallに応じてパージ補正値ｆｐｇが求められる。よって、低開度処理中は、空燃比Ｆ／Ｂ中心のずれに応じたフィードバック処理によりパージ補正値ｆｐｇが求められている。
【０４０６】
さてＥＣＵ５０は、空燃比Ｆ／Ｂ中心のずれが検知されるか、総ベーパ流量Fvpallが所定の上限値Fvpmxに達するかのいずれかが成立しない限りは（Ｓ７３０：ＹＥＳ）、ＶＳＶ開度Dvsvを徐々に開弁駆動する（Ｓ７５２）。このときの開弁速度、即ちＶＳＶ開度Dvsvの増加率は、最大空気層ベーパ流量Fvpairmxに応じ、同流量Fvpairmxが多いほど、即ちパージ中のベーパ濃度が高いと推定されるときほど、ＶＳＶ７１ａがゆっくりと開弁駆動させるように設定される。
【０４０７】
そして空燃比Ｆ／Ｂ中心のずれが検知されたときには（Ｓ７３０：ＮＯ）、ＶＳＶ７１ａの開弁方向への駆動を一時中断して開度を保持し、そのずれに応じて上記流入率rvpdtlを更新する。またそれに合わせて上記総ベーパ流量Fvpallも更新される（Ｓ７４０）。なおここでは、空燃比Ｆ／Ｂ中心値fafav の絶対値が所定のずれ判定値FAFAVHを超えたことをもって、空燃比Ｆ／Ｂ中心のずれを検知している。
【０４０８】
このときの流入率rvpdtlの更新は、その空燃比Ｆ／Ｂ中心のずれを補償すべく増減される（Ｓ７４０）。即ち、図４０の時刻ｔ１、ｔ３、ｔ４等のように、空燃比Ｆ／Ｂの目標値に対して中心値fafav のリーン側へのずれが検知されたときには、そのずれ量に相当する値が上記流入率rvpdtlに加算される。また同図４０の時刻ｔ８のように、上記中心値fafav のリッチ側へのずれが検知されたときには、そのずれ量に相当する値が同推定値rvpdtlから減算される。
【０４０９】
そして、同図４０の時刻ｔ２、ｔ９のように、こうした流入率rvpdtl及び総ベーパ流量Fvpallの更新に伴うパージ補正値ｆｐｇの修正により、空燃比Ｆ／Ｂ中心のずれが解消されると、ＶＳＶ７１ａの開弁方向への駆動が再開される。
【０４１０】
また同図４０の時刻ｔ５〜ｔ６までの期間のように、空燃比Ｆ／Ｂ中心のずれは検知されていなくても、総ベーパ流量Fvpallが上限値Fvpmxに達していれば、ＶＳＶ７１ａの開弁方向への駆動を一時中止してその開度を保持する（Ｓ７３０：ＮＯ）。この上限値Fvpmxは、低開度処理中に許容される総ベーパ流量の上限値であり、試験等により求められた所定の定数として設定されている。
【０４１１】
そして同図４０の時刻ｔ６〜ｔ７までの期間のように、総ベーパ流量Fvpallが上限値Fvpmxに達しており、且つ空燃比Ｆ／Ｂ中心がリッチ側にずれているときには、ＶＳＶ７１ａを所定率で閉弁側に駆動する（Ｓ７５４）。
【０４１２】
なお低開度処理中は、上記求められた総ベーパ流量Fvpallに応じて、ＶＳＶ７１ａの実開度の推定値Dvsvl を下記算出式に従って求めている（Ｓ７６０）。
<<算出式>>
・  Dvsvl   ←  （Fvpall／Fvpairmx）・（Fpgairmx／Fpgmx）
この実開度の推定値Dvsvl が上記ＶＳＶ開度ガード値tDvsvgd を上回ってしまったときには、ＶＳＶ開度Dvsvを再び全閉（０％）とし、そこから再びＶＳＶ７１ａの開弁駆動を再開する。
【０４１３】
そして同図４０の時刻ｔ１０のように、ＶＳＶ開度Dvsvが上記リニアリティ下限値をDVSVL を超えれば、通常のベーパパージ処理に復帰する。なおこのとき、パージ補正値ｆｐｇの不連続的な変化を防止するため、空気層蓄積ベーパ量Mgair の値はそのままに保持しつつ、上記低開度処理終了時の総ベーパ流量Fvpallと一致するように空気層ベーパ流量Fvpairのみを修正する。
【０４１４】
以上が低開度時のＶＳＶ開度制御の詳細である。なお、上記低開度処理中の空燃比Ｆ／Ｂ中心のずれに応じて推定される総ベーパ流量Fpgallや実開度の推定値Dvsvl の値に信頼をおけるなら、上記各物理状態量の更新処理や修正処理をそれらの値に基づいて実施するようにしても良い。無論、低開度処理中にそれら更新処理及び修正処理のいずれか一方の処理のみ禁止し、他方の処理については継続して実施するようにすることもできる。
【０４１５】
なお、低開度時の上記リニアリティの破綻は、ＶＳＶを有したベーパパージシステムを備えるエンジンの空燃比制御装置全般に、生じ得る普遍的な問題である。よって上記低開度処理は、上記実施形態の空燃比制御装置に限らず、ＶＳＶを有したベーパパージシステムを備えるエンジンの空燃比制御装置であれば、同様またはそれに準じた態様で適用可能である。
【０４１６】
［２−９−２］  空燃比フィードバック補正値の中心値の算出処理
次に空燃比Ｆ／Ｂ中心値fafav の算出処理について、上記制御装置に加え得る改良点を、図４１を参照して説明する。
【０４１７】
従来より、空燃比Ｆ／Ｂ中心値fafav は、同図４１（ａ）に示すように、空燃比Ｆ／Ｂ補正値ｆａｆのスキップ時にのみ、その値が更新されている。このため、パージ状況やエンジン運転状況が大きく変化して積分期間が長期化した場合などには、空燃比Ｆ／Ｂ中心値fafav の更新もその間途絶え、上記状況の変化前の値が維持されることとなる。そしてその結果、その空燃比Ｆ／Ｂ中心値fafav を参照して行われる各種処理に好ましくない影響を与えることがある。
【０４１８】
特に上記実施形態の空燃比制御装置では、その影響はより深刻である。
上記実施形態の空燃比制御装置では、上記修正処理（［２−７］参照）を始めとする各種処理を空燃比Ｆ／Ｂ中心のずれに基づいて行っている。そしてそうした処理を通じて設定された上記各物理状態量の値より、パージ補正値ｆｐｇを求めている。このため、上記態様で算出された空燃比Ｆ／Ｂ中心値fafav を用いれば、各物理状態量の更新が上記状況の変化に十分に追従できず、上記物理モデルに基づくベーパパージ処理の採用により向上した空燃比Ｆ／Ｂの精度を充分に維持できなくなってしまう。
【０４１９】
その場合であれ、同図４１（ｂ）に示される態様で、空燃比Ｆ／Ｂ中心値fafav の算出処理を採用すれば、そうした間題を解消可能となる。即ち、同図４１（ｂ）の例では、空燃比Ｆ／Ｂ補正値ｆａｆの積分期間中も、その補正値ｆａｆの振幅を監視し、その中心値fafav の値を更新するようにしている。
【０４２０】
ここでは、補正値ｆａｆの積分期間中に、下記算出式で求められる値fafavlが空燃比Ｆ／Ｂ中心値fafav の現状値よりもその補正値ｆａｆの現状値に近いときには（｜ｆａｆ−fafav ｜＞｜ｆａｆ−fafavl｜）、空燃比Ｆ／Ｂ中心値fafav をその値fafavlに更新する（fafav ←fafavl）。
<<算出式>>
・  fafavl←（ｆａｆ０＋ｆａｆ）／２
ここで「ｆａｆ０」は、その積分期間の前の空燃比Ｆ／Ｂ補正値ｆａｆのスキップ処理時のスキップ中心値、即ちスキップ処理前の補正値ｆａｆとスキップ処理後の補正値ｆａｆとの平均値である。
【０４２１】
［２−９−３］  パージ流量の濃度補正処理
パージライン７１から吸気通路１２にパージされるガスの流量は、実験等により予め吸気通路内圧ＰＭ及びＶＳＶ開度Dvsvとの相関関係を求めておけば、エンジン運転中に実測するまでもなく求めることができる。そこで上記実施形態では、図１８に例示されるような演算マップを用いて許容されるパージ流量の最大値Fpgmxを吸気通路内圧ＰＭから算出し、その最大値FpgmxとＶＳＶ開度Dvsvとの対比のもとに総パージ流量Fpgallを求めている（［２−４−１］参照）。
【０４２２】
ただし厳密に云えば、こうして求められた総パージ流量Fpgallは、パージガスの比重を一定値としたときの体積流量でしかない。ちなみに、上記実施形態では、パージライン７１から吸気通路１２へのパージガスの比重を空気の比重（約１．２ｇ毎リットル）であるとの想定のもとに、図１８に例示される演算マップが作成されている。
【０４２３】
一方、実際には、パージガスの比重はそのパージガスのベーパ含有率、即ちパージガスのベーパ濃度rvpt（＝Fvpall／Fpgall）に応じて変化してしまう。そして上記実施形態では、そのように求められた総パージ流量Fpgallを吸気通路１２にパージされるガスの体積流量［ｇ毎秒］と見なして、上記各種処理を行っている。
【０４２４】
勿論、こうした場合であれ、ベーパパージ処理中の上記パージガスの比重が、上記演算マップの作成時に想定されたパージガスの比重（上記実施形態では空気の比重）に対して大差がなければ、それでも必要な総パージ流量Fpgallの算出精度の確保は充分に可能である。即ち、上記実施形態では、上記ベーパ濃度rvptがある程度よりも小さいことを条件に、総パージ流量Fpgallの算出精度が確保されている。よって上記実施形態では、ベーパ濃度rvptが大きいときには、総パージ流量Fpgallの算出精度のある程度の低下は避け難い。
【０４２５】
その場合であれ、上記総パージ流量Fpgallの算出値を上記パージガスの比重、即ちそのベーパ濃度rvptに応じて適宜補正すれば、ベーパ濃度rvptの変化に拘わらず、その算出精度を維持できる。
【０４２６】
例えば、上記パージガスの比重と空気の比重との比（比重比）と、そのパージガスのベーパ含有率（ベーパ濃度rvpt）との相関関係を予め求めておき、図４２に示されるような演算マップを作成する。そして、現状の総パージ流量Fpgallと総ベーパ流量Fvpallとより、現状のパージガスのベーパ濃度rvptを算出し、その演算マップより、上記比重比を流量補正係数として求める。こうして求められた流量補正係数を、上記最大流量Fpgmx演算用のマップ（図１８）より吸気通路内圧ＰＭに応じて算出された最大流量Fpgmxに乗算し、その値を最終的な最大流量Fpgmxとして総パージ流量Fpgallを算出する。或いは上記実施形態の算出ロジックに従い算出された総パージ流量Fpgallに、その流量補正係数を乗算した値を最終的な総パージ流量Fpgallとする。以上により、パージガスの比重変化を加味した精度の高い総パージ流量Fpgallの算出が可能になる。
【０４２７】
［２−９−４］  各物理状態量の修正エラーの低減処理
上記実施形態では、空燃比Ｆ／Ｂ中心のずれに応じて各物理状態量、即ち空気層蓄積ベーパ量Mgair 、吸着剤蓄積ベーパ量Mgcan 及びタンク発生ベーパ流量Fvptnkの値を補正する修正処理を実施している（［２−７］参照）。こうした修正処理に対して、下記の各修正エラーの低減処理を行えば、各物理状態量の値の精度の更なる向上が許容される。
【０４２８】
（ａ）  吸入空気量Ｇａの影響による修正エラーの低減処理
上記空燃比Ｆ／Ｂ制御にかかる空燃比学習値ＫＧなどに誤差があるときには、吸入空気量Ｇａの増大に伴いその誤差による燃料噴射量Ｑｆｉｎの算出エラーが増幅され、空燃比Ｆ／Ｂのずれも拡大される。そうした状況での空燃比Ｆ／Ｂのずれの拡大に従ってそのまま修正処理を行えば、各物理状態量を過剰に補正してしまい、吸入空気量Ｇａが低減されたときにパージ補正が過剰となるおそれがある。
【０４２９】
こうした問題は、吸入空気量Ｇａが多いほど各物理状態量の修正度合いを小さくするようにすれば、即ち吸入空気量Ｇａが多いときほど空燃比Ｆ／Ｂのずれに対する各物理状態量の修正の度合いを小さくすれば、容易に回避できる。具体的には、以下に例示する各処置の少なくとも１つを採用することで、上記問題の回避が可能となる。
【０４３０】
（ａ−Ｉ）  吸入空気量Ｇａに応じた修正反映率の変更
吸入空気量Ｇａが多いときほど、空燃比Ｆ／Ｂのずれに対する各物理状態量の修正量の比率、即ち修正反映率を小さくすれば、上記問題を回避することができる。修正反映率は、例えば図４３に例示する吸入空気量Ｇａによる演算マップなどを用いて求めることができる。
【０４３１】
（ａ−II）  吸入空気量Ｇａに応じた空燃比Ｆ／Ｂのずれ判定値の変更
図４４に例示されるように、吸入空気量Ｇａが多いときほど空燃比Ｆ／Ｂの各ずれ判定値ERFAFAIR、ERFAFSMAIR、ERFAFSMCAN（［２−７−１］参照）を大きくする処置を行う。この処置によっても、吸入空気量Ｇａが多いときほど、空燃比Ｆ／Ｂ中心のずれに対する各物理状態量の修正度合いが低減され、上記問題を回避することができる。なお、上述の各安定判定値SFFSFSMAIR、SFFAFSMCAN（同じく［２−７−１］参照）についても同様に、吸入空気量Ｇａが多いときほどその値を大きくする処置を行えば、より好適に修正処理を行えるようになる。
【０４３２】
（ｂ）  吸着剤内空気流量Fpgcanの影響による修正エラーの低減処理
上記実施形態では、所定条件が満たされたときの空燃比Ｆ／Ｂ中心のずれを吸着剤脱離パージのベーパ濃度rvpcanの誤差によるものと見なし、そのずれに応じて吸着剤蓄積ベーパ量Mgcan の修正を行っている。ただし厳密には、空燃比Ｆ／Ｂ中心のずれの要因には、空燃比学習値ＫＧの誤差等の他の要因も含まれている可能性がある。このため、吸着剤内空気流量Fpgcanが少ないときに、その空燃比Ｆ／Ｂ中心のずれの要因の全てをベーパ濃度rvpcanの誤差に求めてしまえば、同流量Fpgcanが増大されたときに、上記他の要因による見積もりエラーが増幅され、過補正となるおそれがある。そこで上記実施形態では、吸着剤内空気流量Fpgcanが一定値に満たないときには、吸着剤蓄積ベーパ量Mgcan の修正を禁止することで、かかる間題に対応している（以上、［２−７−２］参照）。
【０４３３】
ただし、修正の許容／禁止といった二者択一の処置だけでは、上記問題に十分に対応しきれないことがある。そうした場合であれ、吸着剤内空気流量Fpgcanに応じ、その流量Fpgcanが少ないときほど吸着剤蓄積ベーパ量Mgcan の修正度合いを小さくするように処置すれば、上記間題により好適に対処できる。例えば図４５に例示されるように、吸着剤蓄積ベーパ量Mgcan 修正用の空燃比Ｆ／Ｂ中心の徐変値fafsm[can]の中心値fafav に対する追従性を、吸着剤内空気流量Fpgcanに応じて、その流量Fpgcanが少ないときほど小さくすることによっても、そうした処置が可能である。
【０４３４】
なお、空気層蓄積ベーパ量Mgair の修正についても、同様の傾向がないとは云えない。そこで空気層パージ流量Fpgairが少ないときほど、空気層蓄積ベーパ量Mgair の修正度合いを小さくするようにすることも考えられる。
【０４３５】
［２−９−５］  タンク発生ベーパの直接流入への対策処理
パージ実施中には、空気層パージ、吸着剤脱離パージの他にも、燃料タンク３０からパージライン７１にベーパが直接流入する可能性がある。上記実施形態では、こうしたタンク流入ベーパは、ごく微量であるとして、総ベーパ流量Fvpallの算出にあたっては無いものと見なして、上記物理モデル（図１３参照）を構築している。
【０４３６】
ただし、より厳密に総ベーパ流量Fvpallを求める必要がある場合、パージシステムの構成やその使用条件などのため上記タンク流入ベーパが無視し得ない流量となる場合などには、図４６に例示するような、タンク流入ベーパも考慮した物理モデルの構築が必要となる。
【０４３７】
ちなみに上記初期化処理の説明で述べたように、ベーパパージ中に許容されるタンク流入ベーパの流量の上限は、タンク発生ベーパ流量Fvptnkに対して一定の割合となると推定される（［２−６−２］参照）。またそのほとんどがベーパ分によって占められ、より高圧なタンク流入ベーパは、空気層パージ及び吸着剤脱離パージに優先してパージされると考えられる。
【０４３８】
よって図４６の物理モデルに基づいたときの総ベーパ流量Fvpallは、例えば下記の算出式に従って求めることができる。
<<算出式>>
・  Fvpttp  ←  Fpgall／RVPTNK（ただし、Fvpttp≦Fvptnk／RVPTNK）
・  rvptnk  ←  RVPTNK・Fvpttp／Fvptnk
・  Fpgair  ←  Fpgall−rvptnk・Fvptnk（ただし、０≦Fpgair≦Fpgairmx）
・  Fpgcan  ←  Fpgall−rvptｎｋ・Fvptnk−Fpgair（ただし、Fpgcan≧０）
・  rvpair  ←  Fpgair／Fpgairmx
・  Fvpair  ←  −rvpair・Fvpairmx
・  Fvpcan  ←  −rvpcan・Fpgcan
・  Fvpall  ←  rvptnk・Fvptnk＋Fvpair＋Fvpcan
ここで「Fvpttp」は、ベーパパージ中に許容されるタンク流入ベーパ流量の最大値である。また「RVPTNK」は、ベーパパージ中に許容されるタンク流入ベーパ流量Fvpttpの上限値について、そのタンク発生ベーパ流量Fvptnkに対するの比率を示す。ここでは、その比率RVPTNKを所定の定数としている。更にここでの「rvptnk」は、タンク発生ベーパ流量Fvptnkに対するタンク流入ベーパ流量の比率を示している。以上に記載した以外のパラメータについては、上記実施形態と同様である。
【０４３９】
更に空気層蓄積パージ量の定期更新処理については、以下の算出式によりその更新量△Mgair を求めることができる。
<<算出式>>
・  △Mgair ←  （１−rvptnk）・Fvptnk＋Fvpcta−Fvpatc−Fvpair
このように、タンク流入ベーパ流量Fvpttpを考慮して算出式等を適宜変更すれば、図４６の物理モデルに従って上記実施形態に記載の各種処理を同様に行うことができる。
【０４４０】
以上が上記実施形態の空燃比制御装置に加え得る各改良点の詳細である。
上記実施形態により得られる効果は以上詳述した通りであるが、その効果のうちの主だったものを以下に列記する。
【０４４１】
（１）上記実施形態のエンジンの空燃比制御装置では、空気層蓄積ベーパ量Mgair 、吸着剤蓄積ベーパ量Mgcan 、及びタンク発生ベーパ流量Fvptnkに基づいたベーパ挙動についての物理モデルに従い、総パージ流量Fpgallに応じた総ベーパ流量Fvpallを推定している。そしてその推定値に応じて燃料噴射量を補正するようにしている。したがって、上記実施形態によれば、パージシステム２０でのベーパ挙動の変化に拘わらず、パージライン７１から吸気通路１２にパージされる総ベーパ流量Fvpallを正確に予測し、パージ実施中の空燃比を高精度に制御できる。
【０４４２】
（２）上記実施形態では、上記物理モデルを用い、パージ状況及び各物理状態量の現状値に基づいて、各物理状態量の値を周期的に更新するようにしている。このため、オープンループの演算処理のみによる、すなわちフィードフォワードによる総ベーパ流量Fvpallの推定が可能となり、空燃比Ｆ／Ｂのずれに基づくフィードバックに依らずとも、上記ベーパ挙動の変化に対応したパージ実施中の空燃比の高精度な制御が可能となる。
【０４４３】
（３）上記実施形態では、ベーパパージ実施中の空燃比Ｆ／Ｂの推移を監視し、その値のずれに応じて各物理状態量の値を修正している。このため、各物理状態量を正確な値に保持し、空燃比を高精度に維持することができる。
【０４４４】
（４）上記実施形態では、空燃比Ｆ／Ｂのずれに基づいて総ベーパ流量の仮値Fvpsを求め、ＶＳＶ７１ａを全閉から徐々に開弁駆動して総パージ流量Fpgallを「０」から徐増していったときの、その仮値Fvpsの推移に基づいて上記各物理状態量の初期値を求めている。このため、各物理状態量の値が不明であれば、それらの値を求めだして、上記物理モデルに基づく制御を実施可能とすることができる。
【０４４５】
（５）上記実施形態では、上記物理モデルに基づく総ベーパ流量Fvpallの推定ロジックを用いて任意のＶＳＶ開度Dvsvでの総ベーパ流量Fvpallを予測しつつ、ＶＳＶ７１ａの開度制御を行っている。これにより、所望とする総ベーパ流量Fvpallを的確に確保可能なように、ＶＳＶ７１ａの開度制御に基づいた総パージ流量Fpgallの調整が可能となる。
【０４４６】
（６）また上記実施形態では、［２−９−１］の処理を更に追加することで、吸気通路内圧ＰＭ及びＶＳＶ開度Dvsvと総パージ流量Fpgallとの相関関係が不明なＶＳＶ７１ａの低開度時には、ＶＳＶ７１ａを一旦全閉とした後、空燃比Ｆ／Ｂの変化度合いに応じてＶＳＶ７１ａの開度制御を行う構成となる。この構成では、正確な総パージ流量Fpgallの把握の困難な状況下であれ、空燃比Ｆ／Ｂ制御への影響を抑えて好適にベーパパージ処理を行うことができる。
【０４４７】
ちなみに、上記実施形態の制御の詳細部分については適宜に変更し、例えば上記［課題を解決するための手段］の欄に記載の各態様で実施可能である。またベーパパージシステムの構成についても、それが上記のような吸着剤、キャニスタ空気層及び大気孔を有するキャニスタを備え、燃料タンクで発生したベーパをキャニスタからパージラインを経由してエンジン吸気系にパージするパージシステムであれば、任意に本発明を適用可能である。
【図面の簡単な説明】
【図１】ベーパパージシステムの基本構成を示す模式図。
【図２】ベーパ流量とＶＳＶ開度との関係を示すグラフ。
【図３】パージ開始からのベーパ流量の推移を示すグラフ。
【図４】吸着ベーパ量とベーパ濃度との関係を示すグラフ。
【図５】パージガス各成分の流量と大気孔からの空気流量との関係を示すグラフ。
【図６】吸着ベーパ量と脱離速度との関係を示すグラフ。
【図７】パージ実施時のパージシステム内でのパージガスの挙動を示すモデル図。
【図８】パージ実施時の空気層パージガスの挙動を示すモデル図。
【図９】パージガス各成分の流量と総パージ流量との関係を示すグラフ。
【図１０】空気層蓄積ベーパ量と空気層ベーパ流量との関係を示すグラフ。
【図１１】定常時のキャニスタ内でのベーパ挙動を示すモデル図。
【図１２】パージガス各成分の流量と総パージ流量との関係を示すグラフ。
【図１３】パージシステム全体でのベーパ挙動を示すモデル図。
【図１４】本発明の一実施形態についてそのパージシステムの全体構成を示す模式図。
【図１５】基本ルーチンの処理手順を示すフローチャート。
【図１６】各パージ流量の算出ロジックを示すブロック図。
【図１７】各ベーパ流量の算出ロジックを示すブロック図。
【図１８】吸気通路内圧と最大総パージ流量との関係を示すグラフ。
【図１９】空気層蓄積ベーパ量と最大空気層パージ流量との関係を示すグラフ。
【図２０】流量の温度補正係数と吸気温度との関係を示すグラフ。
【図２１】空気層蓄積ベーパ量と最大空気層ベーパ流量との関係を示すグラフ。
【図２２】吸着剤蓄積ベーパ量と吸着剤脱離ベーパ濃度との関係を示すグラフ。
【図２３】空燃比制御における制御態様例を示すタイムチャート。
【図２４】物理状態量の初期化処理における制御態様例を示すタイムチャート。
【図２５】総パージ流量と各ベーパ成分の流量との関係を示すグラフ。
【図２６】総パージ流量と各ベーパ成分の流量との関係を示すグラフ。
【図２７】物理状態量の修正ルーチンの処理手順を示すフローチャート。
【図２８】吸着剤蓄積ベーパ量の修正にかかる制御態様例を示すタイムチャート。
【図２９】空気層蓄積ベーパ量の修正にかかる制御態様例を示すタイムチャート。
【図３０】同じく修正にかかる制御態様例を示すタイムチャート。
【図３１】同じく修正にかかる制御態様例を示すタイムチャート。
【図３２】同じく修正にかかる制御態様例を示すタイムチャート。
【図３３】反省処理にかかる制御態様例を示すグラフ。
【図３４】タンク発生ベーパ流量の修正にかかる制御態様例を示すタイムチャート。
【図３５】吸入空気量と絶対ガード値との関係を示すグラフ。
【図３６】ＶＳＶ開度の算出ルーチンの処理手順を示すフローチャート。
【図３７】パージ開始後のＶＳＶ開度及び総ベーパ流量の推移を示すタイムチャート。
【図３８】ＶＳＶ開度と総パージ流量との関係を示すグラフ。
【図３９】低開度時ＶＳＶ制御の処理手順を示すフローチャート。
【図４０】低開度時ＶＳＶ制御の制御態様例を示すタイムチャート。
【図４１】空燃比Ｆ／Ｂ補正値及びその中心値の推移を示すタイムチャート。
【図４２】ベーパ濃度と流量補正係数との関係を示すグラフ。
【図４３】吸入空気量と修正量反映係数との関係を示すグラフ。
【図４４】吸入空気量とずれ判定値との関係を示すグラフ。
【図４５】徐変時定数と総パージ流量との関係を示すグラフ。
【図４６】他の実施形態についてそのパージシステム全体でのベーパ挙動を示すモデル図。
【符号の説明】
１，３０…燃料タンク、２，３５…ベーパライン、３，４０…キャニスタ、３ａ，４２，４３…吸着剤、３ｂ，４５…キャニスタ空気層、４，７１…ベーパライン、５，１０…エンジン、１１…燃焼室、６，１２…吸気通路（エンジン吸気系）、１２ａ…デリバリパイプ（燃料供給装置）、１２ｂ…インジェクタ（燃料供給装置）、３１…燃料ポンプ（燃料供給装置）、７，７１ａ…パージ調整バルブ（ＶＳＶ）、８…大気孔、７０…大気バルブ（大気孔）、７２…大気導入ライン（大気孔）、７２ａ…大気導入バルブ（大気孔）、２０…ベーパパージシステム、５０…ＥＣＵ）。
【図１】
0004400003.tif 000006
【図２】
0004400003.tif 000007
【図３】
0004400003.tif 000008
【図４】
0004400003.tif 000009
【図５】
0004400003.tif 000010
【図６】
0004400003.tif 000011
【図７】
0004400003.tif 000012
【図８】
0004400003.tif 000013
【図９】
0004400003.tif 000014
【図１０】
0004400003.tif 000015
【図１１】
0004400003.tif 000016
【図１２】
0004400003.tif 000017
【図１３】
0004400003.tif 000018
【図１４】
0004400003.tif 000019
【図１５】
0004400003.tif 000020
【図１６】
0004400003.tif 000021
【図１７】
0004400003.tif 000022
【図１８】
0004400003.tif 000023
【図１９】
0004400003.tif 000024
【図２０】
0004400003.tif 000025
【図２１】
0004400003.tif 000026
【図２２】
0004400003.tif 000027
【図２３】
0004400003.tif 000028
【図２４】
0004400003.tif 000029
【図２５】
0004400003.tif 000030
【図２６】
0004400003.tif 000031
【図２７】
0004400003.tif 000032
【図２８】
0004400003.tif 000033
【図２９】
0004400003.tif 000034
【図３０】
0004400003.tif 000035
【図３１】
0004400003.tif 000036
【図３２】
0004400003.tif 000037
【図３３】
0004400003.tif 000038
【図３４】
0004400003.tif 000039
【図３５】
0004400003.tif 000040
【図３６】
0004400003.tif 000041
【図３７】
0004400003.tif 000042
【図３８】
0004400003.tif 000043
【図３９】
0004400003.tif 000044
【図４０】
0004400003.tif 000045
【図４１】
0004400003.tif 000046
【図４２】
0004400003.tif 000047
【図４３】
0004400003.tif 000048
【図４４】
0004400003.tif 000049
【図４５】
0004400003.tif 000050
【図４６】
0004400003.tif 000051