(57)【要約】【課題】過渡運転状態でもその運転状態を良好な状態に維持できるように点火時期等の機関制御パラメータを制御可能な内燃機関の制御装置を提供すること。【解決手段】この内燃機関の制御装置は、吸気管４１等の吸気通路構成部材に対する付着挙動が互いに相違する高沸点成分及び低沸点成分の各々の燃料付着量を区別して推定するとともに、前記推定された各燃料成分の燃料付着量に基いて燃焼室２５内に流入する各燃料成分の筒内流入燃料量を区別して推定する。そして、過渡運転状態にあるとき時々刻々と変化し得る筒内に流入する各燃料成分間の成分比率を求め、同成分比率の定常運転時での成分比率（一定値）に対する比率に基き点火時期補正量を求め、定常運転時にて予め適合されている基本点火時期を前記点火時期補正量だけ補正することで点火時期を求める。これにより、過渡運転状態でも、燃焼状態が最適に維持できるように点火時期が決定される。【選択図】　図１
2004144031.tif 000002
【特許請求の範囲】
【請求項１】
蒸留性状の異なる複数の燃料成分を有する燃料を燃焼室に接続された吸気通路内に噴射する燃料噴射手段を備えた内燃機関に適用される内燃機関の制御装置であって、
前記吸気通路を構成する部材に付着する各燃料成分の燃料付着量を同燃料成分毎に推定する燃料付着量推定手段と、
前記推定された各燃料成分の燃料付着量に基いて前記燃焼室内に流入する同各燃料成分の筒内流入燃料量を同燃料成分毎に推定する筒内流入燃料量推定手段と、
前記推定された各燃料成分の筒内流入燃料量に基いて前記内燃機関の運転状態を良好な状態になるように調整するための機関制御パラメータを決定する機関制御パラメータ決定手段と、
を備えた内燃機関の制御装置。
【請求項２】
請求項１に記載の内燃機関の制御装置において、
前記燃料付着量推定手段は、
前記燃料の燃料成分毎に決定されるとともに前記燃料噴射手段から噴射される前記各燃料成分のうち前記吸気通路構成部材に付着する同各燃料成分の割合を表す各燃料成分の吸気通路付着率と、
前記燃料の燃料成分毎に決定されるとともに前記吸気通路構成部材に付着している前記各燃料成分のうち同吸気通路構成部材に付着したまま残留する同各燃料成分の割合を表す各燃料成分の吸気通路残留率と、
を用いて前記吸気通路構成部材に付着する各燃料成分の燃料付着量を推定するように構成された内燃機関の制御装置。
【請求項３】
請求項１又は請求項２に記載の内燃機関の制御装置であって、
前記内燃機関の運転状態が定常状態にあるときに同運転状態が良好な状態になる同内燃機関の基本点火時期を記憶した基本点火時期記憶手段を備え、
前記機関制御パラメータ決定手段は、前記機関制御パラメータとしての前記内燃機関の点火時期を、前記推定された各燃料成分の筒内流入燃料量に基いて前記基本点火時期から補正することで決定するように構成された内燃機関の制御装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、内燃機関の制御装置に係り、特に内燃機関の運転状態が良好な状態になるように点火時期等の機関制御パラメータを制御する制御装置に関する。
【０００２】
【従来の技術】
内燃機関のノック限界は同内燃機関に使用される燃料（ガソリン）のオクタン価に応じて変化し、使用される燃料のオクタン価が高くなるほどノック限界は進角側に移動する。従って、内燃機関において最適な燃焼状態を得るためには、使用される燃料のオクタン価に応じて点火時期を変更する必要がある。このため、下記特許文献１に開示された内燃機関の点火時期制御装置は、レギュラガソリンと同レギュラガソリンよりもオクタン価が高いプレミアムガソリンの何れが内燃機関の燃料として使用されているかを検出する燃料性状センサを備え、同燃料性状センサの検出結果に応じて内燃機関の点火時期を変更するようになっている。
【０００３】
【特許文献１】
特開平５−７９４３９号公報
【０００４】
【発明が解決しようとする課題】
ところで、内燃機関の燃料として広く使用されている燃料（ガソリン）は、炭素数が少なくて沸点が低い低沸点成分（例えば、メタン、エタン）や炭素数が比較的多くて沸点が比較的高い高沸点成分（例えば、ベンゼン）等のように、蒸留性状（沸点、飽和蒸気圧）の異なる複数の燃料成分から構成される混合物である。また、通常、前記低沸点成分はオクタン価が低く、前記高沸点成分は同低沸点成分よりもオクタン価が高い。換言すれば、蒸留性状の異なる各燃料成分のオクタン価は燃料成分毎に異なる。従って、混合物である燃料（ガソリン）自体（全体）のオクタン価は、これら複数の燃料成分の各々のオクタン価と、それらの成分比率とに応じて変化する。
【０００５】
一方、燃料を燃焼室に接続された吸気通路内に噴射する燃料噴射手段（インジェクタ）を備えた内燃機関においては、燃料噴射手段から噴射された燃料は、その一部が吸気通路壁面や吸気弁等の吸気通路構成部材に付着する。その際の吸気通路構成部材に対する燃料の付着挙動は、同燃料の各燃料成分の蒸留性状が異なることから同燃料成分毎に異なる。特に、内燃機関が過渡運転状態にあって吸気通路内の圧力や温度が急激に変化すると、各燃料成分の間で気化又は液化の程度が急激に異なるようになり、吸気通路構成部材に対する各燃料成分の付着挙動が各燃料成分の間で急激に異なるようになる。その結果、シリンダ内（筒内）に流入する各燃料成分の間の成分比率が一時的に変化し、筒内に流入する（筒内で実際に燃焼される）燃料全体としてのオクタン価（以下、「筒内オクタン価」と称呼する。）も使用される燃料自体のオクタン価から一時的に変化する。
【０００６】
従って、上記開示された装置においては、使用される燃料（自体）のオクタン価に応じて点火時期を決定しても、特に、過渡運転状態において、筒内に流入する燃料の筒内オクタン価が燃料自体のオクタン価から一時的に変化することから最適な燃焼状態を得ることができない場合（決定された点火時期が最適な点火時期とはならない場合）があるという問題がある。
【０００７】
従って、本発明の目的は、内燃機関が過渡運転状態にあるときでも、その運転状態を良好な状態に維持できるように点火時期等の機関制御パラメータを制御することができる内燃機関の制御装置を提供することにある。
【０００８】
【本発明の概要】
本発明の特徴は、蒸留性状の異なる複数の燃料成分を有する燃料を燃焼室に接続された吸気通路内に噴射する燃料噴射手段を備えた内燃機関に適用される内燃機関の制御装置が、前記吸気通路を構成する部材に付着する各燃料成分の燃料付着量を同燃料成分毎に推定する燃料付着量推定手段と、前記推定された各燃料成分の燃料付着量に基いて前記燃焼室内に流入する同各燃料成分の筒内流入燃料量を同燃料成分毎に推定する筒内流入燃料量推定手段と、前記推定された各燃料成分の筒内流入燃料量に基いて前記内燃機関の運転状態を良好な状態になるように調整するための機関制御パラメータを決定する機関制御パラメータ決定手段とを備えたことにある。
【０００９】
この場合、前記吸気通路を構成する部材は、吸気管（インテークマニホールド）、吸気弁（特に、吸気弁の背面）、吸気通路に設けられるＳＣＶ等の吸気制御弁を含んでよい。また、前記機関制御パラメータは、例えば、点火時期、目標空燃比、吸気弁及び排気弁の開閉タイミング、燃料の噴射時期であって、これらに限定されない。
【００１０】
また、上記内燃機関の制御装置においては、前記燃料付着量推定手段は、前記燃料の燃料成分毎に決定されるとともに前記燃料噴射手段から噴射される前記各燃料成分のうち前記吸気通路構成部材に付着する同各燃料成分の割合を表す各燃料成分の吸気通路付着率と、前記燃料の燃料成分毎に決定されるとともに前記吸気通路構成部材に付着している前記各燃料成分のうち同吸気通路構成部材に付着したまま残留する同各燃料成分の割合を表す各燃料成分の吸気通路残留率と、を用いて前記吸気通路構成部材に付着する各燃料成分の燃料付着量を推定するように構成されることが好適である。
【００１１】
吸気通路構成部材に対する燃料の付着挙動（例えば、吸気通路付着率、吸気通路残留率等を用いる燃料付着モデルで表される。）は、燃料を構成する蒸留性状の異なる複数の燃料成分毎に異なる。従って、上記構成のように、吸気通路構成部材に付着する各燃料成分の燃料付着量を同燃料成分毎に区別して推定することにより、各燃料成分の吸気通路構成部材への燃料付着量がそれぞれ精度良く推定される。
【００１２】
また、上記構成によれば、前記推定された各燃料成分の燃料付着量に基いて燃焼室内に流入する同各燃料成分の筒内流入燃料量が同燃料成分毎に区別して推定されるので、各燃料成分の筒内流入燃料量もそれぞれ精度良く推定される。従って、内燃機関が過渡運転状態にあるときでも、例えば、時々刻々と変化するシリンダ内（筒内）に流入する燃料の各燃料成分の間の成分比率を精度良く求めることができ、この結果、筒内に流入する燃料の筒内オクタン価を精度良く求めることができる。
【００１３】
よって、上記のように、前記推定された各燃料成分の筒内流入燃料量に基いて前記内燃機関の運転状態を良好な状態になるように調整するための機関制御パラメータを決定することにより、例えば、時々刻々と変化する筒内に流入する各燃料成分の間の成分比率（筒内に流入する燃料の筒内オクタン価）に基いて点火時期等の機関制御パラメータを決定することができる。従って、内燃機関が過渡運転状態にあるときでも、その運転状態を良好な状態に維持できるように点火時期等の機関制御パラメータを制御することができる。
【００１４】
また、上記内燃機関の制御装置が、前記内燃機関の運転状態が定常状態にあるときに同運転状態が良好な状態になる同内燃機関の基本点火時期を記憶した基本点火時期記憶手段を備える場合には、前記機関制御パラメータ決定手段は、前記機関制御パラメータとしての前記内燃機関の点火時期を、前記推定された各燃料成分の筒内流入燃料量に基いて前記基本点火時期から補正することで決定するように構成されることが望ましい。
【００１５】
内燃機関における最適な点火時期は、機関の回転速度、吸入空気量、機関の温度（冷却水温）等のノック限界に影響を与える複数の因子の組み合わせに応じて変化する。従って、一般には、これら複数の因子の組み合わせを種々変更した各々の場合において内燃機関が定常状態にあるときに最適な燃焼状態（運転状態）が得られるような基本点火時期を予め実験等により求めておき、内燃機関の制御装置は、それら種々の場合における基本点火時期が格納された基本点火時期記憶手段（ＲＯＭ内に記憶されたテーブル）を搭載することで、同基本点火時期記憶手段とその時点での上記複数の因子とに基き同時点での点火時期を決定するようになっている。
【００１６】
よって、内燃機関の制御装置がかかる基本点火時期記憶手段を備えている場合には、上記のように、内燃機関の点火時期を、前記推定された各燃料成分の筒内流入燃料量に基いて（即ち、例えば、各燃料成分の筒内流入燃料量に基いて得られる筒内に流入する各燃料成分の間の成分比率、筒内に流入する燃料の筒内オクタン価に基いて）、基本点火時期から補正することで決定することにより、基本点火時期を決定するための既存のテーブルをそのまま利用しつつ容易に最適な点火時期を計算することが可能となる。
【００１７】
【発明の実施の形態】
以下、本発明による内燃機関の制御装置の実施形態について図面を参照しつつ説明すると、図１は、同制御装置を火花点火式多気筒（４気筒）内燃機関１０に適用したシステムの概略構成を示している。
【００１８】
この内燃機関１０は、シリンダブロック、シリンダブロックロワーケース、及びオイルパン等を含むシリンダブロック部２０と、シリンダブロック部２０の上に固定されるシリンダヘッド部３０と、シリンダブロック部２０にガソリン混合気を供給するための吸気系統４０と、シリンダブロック部２０からの排ガスを外部に放出するための排気系統５０とを含んでいる。
【００１９】
シリンダブロック部２０は、シリンダ２１、ピストン２２、コンロッド２３、及びクランク軸２４を含んでいる。ピストン２２はシリンダ２１内を往復動し、ピストン２２の往復動がコンロッド２３を介してクランク軸２４に伝達され、これにより同クランク軸２４が回転するようになっている。シリンダ２１とピストン２２のヘッドは、シリンダヘッド部３０とともに燃焼室２５を形成している。
【００２０】
シリンダヘッド部３０は、燃焼室２５に連通した吸気ポート３１、吸気ポート３１を開閉する吸気弁３２、吸気弁３２を駆動するインテークカムシャフトを含むとともに同インテークカムシャフトの位相角を連続的に変更する可変吸気タイミング装置３３、可変吸気タイミング装置３３のアクチュエータ３３ａ、燃焼室２５に連通した排気ポート３４、排気ポート３４を開閉する排気弁３５、排気弁３５を駆動するエキゾーストカムシャフト３６、点火プラグ３７、点火プラグ３７に与える高電圧を発生するイグニッションコイルを含むイグナイタ３８、及び燃料を吸気ポート３１内に噴射するインジェクタ（燃料噴射手段）３９を備えている。
【００２１】
吸気系統４０は、吸気ポート３１に連通し同吸気ポート３１とともに吸気通路を形成するインテークマニホールドを含む吸気管４１、吸気管４１の端部に設けられたエアフィルタ４２、吸気管４１内にあって吸気通路の開口断面積を可変とするスロットルバルブ４３、及びスワールコントロールバルブ（以下、「ＳＣＶ」と称呼する。）４４を備えている。スロットルバルブ４３は、ＤＣモータからなるスロットルバルブアクチュエータ４３ａにより吸気管４１内で回転駆動されるようになっている。ＳＣＶ４４は、前記スロットルバルブ４３よりも下流で前記インジェクタ３９よりも上流の位置にて前記吸気管４１に対し回動可能に支持されるとともに、ＤＣモータからなるＳＣＶアクチュエータ４４ａにより回転駆動されるようになっている。なお、本明細書においては、インテークマニホールドを含む吸気管４１、吸気ポート３１、吸気弁３２、及びＳＣＶ４４等を吸気通路を構成する部材（吸気通路構成部材）と称呼する。
【００２２】
図２は、一つの気筒（特定の気筒）の燃焼室２５、及び同燃焼室２５の近傍部分の概略平面図である。図２に示したように、前記吸気ポート３１は、実際には各気筒に一対ずつ設けられた吸気ポート３１ａ，３１ｂからなっている。吸気ポート３１ａは、燃焼室２５内にスワール（旋回流）を発生させるようにヘリカル状に形成され所謂スワールポートを構成し、吸気ポート３１ｂは所謂ストレートポートを構成している。吸気管４１のサージタンクから各燃焼室２５に至る部分（即ち、インテークマニホールドの一部）には、吸気管４１の長手方向に沿って伸びる隔壁４１ａが形成されていて、これにより吸気管４１は吸気ポート３１ａに連通する第１インテークマニホールド４５と、吸気ポート３１ｂに連通する第２インテークマニホールド４６とに区画されている。前記ＳＣＶ４４は、第２インテークマニホールド４６内において回動可能に支持され、第２インテークマニホールド４６の開口断面積を変更し得るようになっている。
【００２３】
隔壁４１ａの適宜個所には第１，第２インテークマニホールド４５，４６を連通する連通路４１ｂが形成されていて、前記インジェクタ３９は同連通路４１ｂの近傍位置に固定され、吸気ポート３１ａ，３１ｂに向けて燃料を噴射するようになっている。これにより、インジェクタ３９から噴射された燃料は、その一部が前記吸気通路構成部材に付着する。
【００２４】
再び図１を参照すると、排気系統５０は、排気ポート３４に連通したエキゾーストマニホールド５１、エキゾーストマニホールド５１に接続されたエキゾーストパイプ５２、及びエキゾーストパイプ５２に介装された触媒コンバータ（三元触媒装置）５３を備えている。
【００２５】
一方、このシステムは、熱線式エアフローメータ６１、吸気温センサ６２、大気圧センサ（スロットルバルブ上流圧力センサ）６３、スロットルポジションセンサ６４、ＳＣＶ開度センサ６５、カムポジションセンサ６６、クランクポジションセンサ６７、水温センサ６８、Ｏ２センサ６９、及びアクセル開度センサ８１を備えている。
【００２６】
エアフローメータ６１は、概略斜視図である図３に示したように、吸気管４１内を流れる吸入空気の一部をバイパスさせるバイパス通路と、このバイパス通路にバイパスされた吸入空気の質量流量を計測する熱線計量部６１ａと、計測された質量流量に応じた電圧Ｖｇを出力する信号処理部６１ｂとからなっている。熱線計量部６１ａは、その拡大斜視図である図４に示したように、白金熱線からなる吸気温計測用抵抗（ボビン部）６１ａ１と、同吸気温計測用抵抗６１ａ１を前記信号処理部６１ｂに連結して保持するサポート部６１ａ２と、加熱用抵抗（ヒータ）６１ａ３と、同加熱用抵抗６１ａ３を前記信号処理部６１ｂに連結して保持するサポート部６１ａ４とを備えている。信号処理部６１ｂは、吸気温計測用抵抗６１ａ１と加熱用抵抗６１ａ３とで構成されたブリッジ回路を備え、このブリッジ回路により吸気温計測用抵抗６１ａ１と加熱用抵抗６１ａ３との温度差を常に一定に維持するように同加熱用抵抗６１ａ３に供給する電力を調整するとともに、この供給する電力を前記電圧Ｖｇに変換して出力するようになっている。
【００２７】
吸気温センサ６２は、エアフローメータ６１内に備えられていて、吸入空気の温度を検出し、吸気温度ＴＨＡを表す信号を出力するようになっている。大気圧センサ６３は、スロットルバルブ４３の上流の圧力（即ち、大気圧）を検出し、スロットルバルブ上流圧力Ｐａを表す信号を出力するようになっている。スロットルポジションセンサ６４は、スロットルバルブ４３の開度を検出し、スロットルバルブ開度ＴＡを表す信号を出力するようになっている。ＳＣＶ開度センサ６５は、ＳＣＶ４４の開度を検出し、ＳＣＶ開度θｉｖを表す信号を出力するようになっている。カムポジションセンサ６６は、インテークカムシャフトが９０°回転する毎に（即ち、クランク軸２４が１８０°回転する毎に）一つのパルスを有する信号（Ｇ２信号）を発生するようになっている。クランクポジションセンサ６７は、クランク軸２４が１０°回転する毎に幅狭のパルスを有するとともに同クランク軸２４が３６０°回転する毎に幅広のパルスを有する信号を出力するようになっている。この信号は、エンジン回転速度Ｎｅを表す。水温センサ６８は、内燃機関１０の冷却水の温度を検出し、冷却水温ＴＨＷを表す信号を出力するようになっている。Ｏ２センサ６９は、触媒コンバータ５３に流入する排ガス中の酸素濃度に応じた信号（空燃比が理論空燃比よりもリッチであるか、又はリーンであるかを示す信号）を出力するようになっている。アクセル開度センサ８１は、運転者によって操作されるアクセルペダル８２の操作量Ａｃｃｐを表す信号を出力するようになっている。
【００２８】
電気制御装置７０は、互いにバスで接続されたＣＰＵ７１、ＣＰＵ７１が実行するプログラム、テーブル（ルックアップテーブル、マップ）、定数等を予め記憶したＲＯＭ７２、ＣＰＵ７１が必要に応じてデータを一時的に格納するＲＡＭ７３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ７４、及びＡＤコンバータを含むインターフェース７５等からなるマイクロコンピュータである。インターフェース７５は、前記センサ６１〜６９，８１と接続され、ＣＰＵ７１にセンサ６１〜６９，８１からの信号を供給するとともに、同ＣＰＵ７１の指示に応じて可変吸気タイミング装置３３のアクチュエータ３３ａ、イグナイタ３８、インジェクタ３９、スロットルバルブアクチュエータ４３ａ、及びＳＣＶアクチュエータ４４ａに駆動信号を送出するようになっている。
【００２９】
次に、上記のように構成された制御装置による点火プラグ３７の点火時期の決定方法について説明する。以下に述べる処理は、ＣＰＵ７１がプログラムを実行することによりなされる。
【００３０】
（吸入空気量相当値ＫＬＦＷＤの決定方法の概要）
この制御装置は、ノック限界に影響を与える因子の一部であるエンジン回転速度Ｎｅ、冷却水温ＴＨＷ、及び一気筒当たりの吸入空気量相当値ＫＬＦＷＤの組み合わせを種々変更した各々の場合において内燃機関１０が定常状態にあるときに最適な燃焼状態を得るための予め実験にて求められている基本点火時期ＩＴｓｔｅａｄｙと、エンジン回転速度Ｎｅ、冷却水温ＴＨＷ、及び吸入空気量相当値ＫＬＦＷＤとの関係を規定する基本点火時期テーブルをＲＯＭ７２に格納していて、同基本点火時期テーブルと、エンジン回転速度Ｎｅ、冷却水温ＴＨＷ、及び一気筒当たりの吸入空気量相当値ＫＬＦＷＤとに基き基本点火時期ＩＴｓｔｅａｄｙを決定する。ここで、基本点火時期テーブルは、基本点火時期記憶手段に相当している。
【００３１】
ここで、本制御装置は、実際のエンジン回転速度Ｎｅをクランクポジションセンサ６７の出力に基いて取得し、また、実際の冷却水温ＴＨＷを水温センサ６８から取得する。しかしながら、実際の一気筒当たりの吸入空気量（吸入空気量相当値ＫＬＦＷＤ）を直接検出することはできないので、本制御装置は、吸気行程にある気筒の吸気弁３２が閉じた時点で（即ち、吸気弁閉時に）同気筒内に吸入されているであろう吸入空気量（一気筒当たりの吸入空気量、筒内吸入空気量）を予測する必要がある。
【００３２】
一方、吸気弁閉時の吸気管圧力ＰＭＦＷＤは、燃焼室２５に吸入されている空気量と比例関係にある。従って、吸気管圧力ＰＭＦＷＤを予測することができれば、実際の筒内吸入空気量を推定することができる。そこで、本制御装置は、吸気弁閉時の吸気管圧力ＰＭＦＷＤを予測・推定し、推定した吸気管圧力ＰＭＦＷＤを一気筒の排気量と空気密度の積で除することにより一気筒当たりの吸入空気量に相当する値ＫＬＦＷＤを求める。
【００３３】
なお、エアフローメータ６１の出力電圧Ｖｇと吸入空気量ｍｔＡＦＭとの関係を規定した図５に示したＶｇ−ｍｔＡＦＭ変換テーブルと、エアフローメータ６１の実際の出力電圧Ｖｇとに基づいて現時点での吸入空気量ｍｔＡＦＭを求め、この吸入空気量ｍｔＡＦＭが一気筒当たりの吸入空気量相当値ＫＬＦＷＤと等しいものとして、吸入空気量相当値ＫＬＦＷＤを簡易的に求めてもよい。
【００３４】
上記のように、吸気弁閉時の吸気管圧力ＰＭＦＷＤから一気筒当たりの吸入空気量相当値ＫＬＦＷＤを求める本実施形態の制御装置は、以下のようにして吸気弁閉時の吸気管圧力ＰＭＦＷＤを予測する。即ち、図６に示したように、吸気弁閉時のスロットルバルブ開度ＴＡＳを予測するとともに、同予測したスロットルバルブ開度ＴＡＳ及びエンジン回転速度Ｎｅ等から吸気弁閉時の吸気管圧力Ｐｍ１を所定のモデルを用いて推定する。また、現時点においてスロットルポジションセンサ６４が検出する実際のスロットルバルブ開度ＴＡＲとエンジン回転速度Ｎｅ等に基づき、エアフローメータ６１が現時点で出力するであろう値を推定し、この推定値に基づいて現時点の吸気管圧力Ｐｍ２を推定する。同時に、現時点においてエアフローメータ６１の実際の出力電圧Ｖｇに基づいて現時点の吸気管圧力Ｐｍ３を推定する。最後に、下記数１にしたがって吸気弁閉時の吸気管圧力ＰＭＦＷＤを求める。これにより、スロットルバルブ開度の予測値ＴＡＳに基づく推定値である吸気管圧力Ｐｍ１に含まれる定常的な誤差を、エアフローメータ６１の実際の出力電圧Ｖｇにより補正し、吸気弁閉時の吸気管圧力ＰＭＦＷＤを精度良く推定する。
【００３５】
【数１】
ＰＭＦＷＤ＝Ｐｍ３＋（Ｐｍ１−Ｐｍ２）
【００３６】
なお、スロットルバルブ開度が一定に維持されていて内燃機関１０が定常状態にある場合、吸気管圧力Ｐｍ１と吸気管圧力Ｐｍ２は等しくなるので、上記数１から理解されるように、吸気管圧力ＰＭＦＷＤは吸気管圧力Ｐｍ３と等しくなる。換言すると、定常運転状態では、実質的にエアフローメータ６１の出力電圧Ｖｇに基づいて吸気弁閉時の吸気管圧力ＰＭＦＷＤが決定されることになる。
【００３７】
以下、各吸気管圧力Ｐｍ１，Ｐｍ２，Ｐｍ３の推定方法について、同推定に使用するモデルとともに説明する。
【００３８】
（Ｐｍ１の求め方）
図７に示したように、吸気管圧力Ｐｍ１は電子制御スロットルモデルＭ１、スロットルモデルＭ２、吸気弁モデルＭ３、及びインテークマニホールドモデルＭ４により推定される。
【００３９】
（１）電子制御スロットルモデルＭ１
電子制御スロットルモデルＭ１は、現時点までのアクセルペダル操作量Ａｃｃｐに基づいて吸気弁閉時のスロットルバルブ開度ＴＡＳを推定するモデルである。本実施形態においては、スロットルバルブ電子制御ロジックＡ１にて、アクセル開度センサ８１により検出されたアクセルペダル操作量Ａｃｃｐと、図８に示したアクセルペダル操作量Ａｃｃｐと目標スロットルバルブ開度θｒとの関係を規定するテーブルとに基づいて暫定的な目標スロットルバルブ開度θｒ１が求められ、この暫定的な目標スロットルバルブ開度θｒ１を所定時間Ｔ（例えば、６４ｍｓｅｃ）だけ遅延させた値が最終的な目標スロットルバルブ開度θｒとして決定される。そして、スロットルバルブ電子制御ロジックＡ１（電気制御装置７０）は、実際のスロットルバルブ開度ＴＡが目標スロットルバルブ開度θｒとなるようにスロットルバルブアクチュエータ４３ａに対して駆動信号を送出する。
【００４０】
このように、目標スロットルバルブ開度θｒは、現時点から所定時間Ｔだけ前の時点におけるアクセルペダル操作量Ａｃｃｐに応じて決定されるから、現時点から吸気弁閉時までの時間をｔとすると、吸気弁閉時の目標スロットルバルブ開度θｒは、現時点から時間（Ｔ−ｔ）前における暫定的な目標スロットルバルブ開度θｒ１と等しい。また、目標スロットルバルブ開度θｒは、スロットルバルブアクチュエータ４３ａの作動遅れ時間を無視すれば、スロットルバルブ開度ＴＡＳと等しい。このような考えに基づき、電子制御スロットルモデルＭ１は、検出されるエンジン回転速度Ｎｅと、内燃機関１０の運転状態に応じて別途定められる吸気弁の開閉タイミング（進角量）ＶＴ（上記信号Ｎｅと上記Ｇ２信号とにより求めた実際の開閉タイミングＶＴでも良い。）と等に基づいて現時点から吸気弁閉時までの時間ｔを求め、同時間ｔと、現時点から所定時間Ｔだけ前の時点から現時点までのアクセルペダル操作量Ａｃｃｐ（又は、暫定的な目標スロットルバルブ開度θｒ１）の変化の経緯とに基づいて吸気弁閉時のスロットルバルブ開度ＴＡＳを推定する。なお、スロットルバルブアクチュエータ４３ａの作動遅れ時間を考慮に加えて、吸気弁閉時のスロットルバルブ開度ＴＡＳを推定してもよい。
【００４１】
（２）スロットルモデルＭ２
スロットルモデルＭ２は、スロットルバルブ４３を通過する空気量（スロットル通過空気量）ｍｔを、エネルギー保存則、運動量保存則、質量保存則、及び状態方程式に基づいて得られた下記数２及び下記数３に基づいて推定するモデルである。下記数２及び下記数３において、μは流量係数、Ａｔはスロットル開口面積、νはスロットルバルブ４３を通過する空気の流速、Ｐａはスロットルバルブ上流圧力、Ｐｍは吸気管圧力、Ｔａは吸気温度、ρｍは吸気密度、Ｒは気体定数、及びκは比熱比（以下、κを一定値として扱う。）である。
【００４２】
【数２】
ｍｔ＝μ・Ａｔ・ν・ρｍ＝μ・Ａｔ・｛Ｐａ／（Ｒ・Ｔａ）１／２｝・Φ（Ｐｍ／Ｐａ）
【００４３】
【数３】
2004144031.tif 000003
【００４４】
ここで、上記数２は、ｋ１を所定の係数（＝μ・Ａｔ・｛Ｐａ／（Ｒ・Ｔａ）１／２｝）、ｍｔｓを吸気弁閉時のスロットル通過空気量とするとき下記数４に書き換えられる。また、数４において、内燃機関１０が定常状態にある場合（スロットルバルブ開度が一定である場合）のスロットル通過空気量をｍｔｓＴＡ、及び吸気管圧力をＰｍＴＡとすると、下記数５が得られるので、数４及び数５から係数ｋ１を消去して下記数６を得ることができる。
【００４５】
【数４】
ｍｔｓ＝ｋ１・Φ（Ｐｍ／Ｐａ）
【００４６】
【数５】
ｍｔｓＴＡ＝ｋ１・Φ（ＰｍＴＡ／Ｐａ）
【００４７】
【数６】
ｍｔｓ＝｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝・Φ（Ｐｍ／Ｐａ）
【００４８】
上記数６の右辺における値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝は、スロットルバルブ開度ＴＡが一定であるときの吸入空気流量（スロットル通過空気量）に関する値であり、スロットルバルブ開度ＴＡ、エンジン回転速度Ｎｅ、吸気弁の開閉タイミングＶＴ、及びスロットルバルブ上流圧力Ｐａが決定されると、実質的に一意に定まる値である。スロットルモデルＭ２は、スロットルバルブ開度ＴＡ、エンジン回転速度Ｎｅ、吸気弁の開閉タイミングＶＴ、及びスロットルバルブ上流圧力Ｐａと、値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝との関係を規定したテーブルをＲＯＭ７２内に記憶していて、このテーブルと吸気弁閉時の推定スロットルバルブ開度ＴＡＳ、実際のエンジン回転速度Ｎｅ、実際の吸気弁の開閉タイミングＶＴ、及び実際のスロットルバルブ上流圧力Ｐａとに基づいて値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝を求める。
【００４９】
また、数６の右辺における値Φ（Ｐｍ／Ｐａ）は、上記数３から理解されるように、比熱比κが一定であるとき、吸気管圧力Ｐｍとスロットルバルブ上流圧力Ｐａにより決定される値である。スロットルモデルＭ２は、吸気管圧力Ｐｍ及びスロットルバルブ上流圧力Ｐａと、値Φ（Ｐｍ／Ｐａ）との関係を規定したテーブルをＲＯＭ７２内に記憶していて、このテーブルと、後述するインテークマニホールドモデルＭ４が現時点で既に演算している最新の吸気管圧力Ｐｍ、及び実際のスロットルバルブ上流圧力Ｐａとに基づいて値Φ（Ｐｍ／Ｐａ）を求める。以上により、吸気弁閉時のスロットル通過空気量ｍｔｓが求められる。
【００５０】
（３）吸気弁モデルＭ３
吸気弁モデルＭ３は、吸気管圧力Ｐｍ、吸気管内温度Ｔｍ、及び吸気温度ＴＨＡ等から筒内吸入空気流量ｍｃを推定するモデルである。吸気弁閉弁時の気筒内圧力は吸気弁３２の上流の圧力、即ち吸気弁閉時の吸気管圧力Ｐｍとみなすことができるので、筒内吸入空気流量ｍｃは吸気弁閉時の吸気管圧力Ｐｍに比例する。そこで、吸気弁モデルＭ３は筒内吸入空気流量ｍｃを、経験則に基づく下記数７にしたがって求める。
【００５１】
【数７】
ｍｃ＝（ＴＨＡ／Ｔｍ）・（ｃ・Ｐｍ−ｄ）
【００５２】
数７において、値ｃは比例係数、値ｄは筒内に残存していた既燃ガス量である。吸気弁モデルＭ３は、エンジン回転速度Ｎｅ、及び吸気弁の開閉タイミングＶＴと、比例係数ｃ、及び既燃ガス量ｄとの関係をそれぞれ規定するテーブルをＲＯＭ７２内に格納していて、実際のエンジン回転速度Ｎｅと、実際の吸気弁の開閉タイミングＶＴと前記格納しているテーブルとから比例係数ｃ、及び既燃ガス量ｄを求める。また、吸気弁モデルＭ３は、演算時点において、後述するインテークマニホールドモデルＭ４により既に推定されている直前（最新）の吸気弁閉時の吸気管圧力Ｐｍと直前の吸気管内空気温度Ｔｍとを上記数７に適用し、吸気弁閉時の筒内吸入空気流量ｍｃを推定する。
【００５３】
（４）インテークマニホールドモデルＭ４
インテークマニホールドモデルＭ４は、質量保存則とエネルギー保存則とにそれぞれ基づいた下記数８及び下記数９にしたがって、吸気弁閉時の吸気管圧力Ｐｍと、吸気弁閉時の吸気管内温度Ｔｍとを求める。なお、Ｖは吸気管の容積、Ｒは気体定数、ｍｔはスロットル通過空気量、Ｔａはスロットルバルブ通過空気温度（即ち、吸気温度ＴＨＡ）である。
【００５４】
【数８】
ｄＰｍ／ｄｔ＝κ・（Ｒ／Ｖ）・（ｍｔ・Ｔａ−ｍｃ・Ｔｍ）
【００５５】
【数９】
ｄ（Ｐｍ／Ｔｍ）／ｄｔ＝（Ｒ／Ｖ）・（ｍｔ−ｍｃ）
【００５６】
図７に示したように、インテークマニホールドモデルＭ４は、スロットルモデルＭ２により推定されたスロットル通過空気量ｍｔｓを上記数８，数９におけるスロットル通過空気量ｍｔとして使用し、吸気弁モデルＭ３により推定された吸気弁閉時の筒内吸入空気流量ｍｃを上記数８，数９の筒内吸入空気流量ｍｃとして使用する。このインテークマニホールドモデルＭ４により推定された吸気管圧力Ｐｍが、前記吸気弁閉時の推定吸気管圧力Ｐｍ１となる。
【００５７】
（Ｐｍ２の求め方）
上記エアフローメータ６１が現時点で出力するであろう値に基づく吸気管圧力Ｐｍ２は、上記スロットルモデルＭ２と同じモデルであるスロットルモデルＭ５、エアフローメータモデルＭ６、上記吸気弁モデルＭ３と同じ吸気弁モデルＭ７、及び上記インテークマニホールドモデルＭ４と同じインテークマニホールドモデルＭ８により求められる。
【００５８】
（５）スロットルモデルＭ５
具体的に述べると、スロットルモデルＭ５は、上記数６を書換えた下記数１０に従って、現時点におけるスロットル通過空気量ｍｔＴＨＲを推定する。
【００５９】
【数１０】
ｍｔＴＨＲ＝｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝・Φ（Ｐｍ／Ｐａ）
【００６０】
スロットルモデルＭ５は、上記数１０の右辺における値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝を、スロットルバルブ開度ＴＡ、エンジン回転速度Ｎｅ、吸気弁の開閉タイミングＶＴ、及びスロットルバルブ上流圧力Ｐａと、値｛ｍｔｓＴＡ／Φ（ＰｍＴＡ／Ｐａ）｝との関係を規定した前記テーブルと、スロットルポジションセンサ６４が実際に検出したスロットルバルブ開度ＴＡ（以下、「実スロットルバルブ開度ＴＡＲ」と称呼する。）、実際のエンジン回転速度Ｎｅ、実際の又は計算された吸気弁の開閉タイミングＶＴ、及び実際のスロットルバルブ上流圧力Ｐａとに基づいて求める。
【００６１】
また、スロットルモデルＭ５は、数１０の右辺における値Φ（Ｐｍ／Ｐａ）を、吸気管圧力Ｐｍ及びスロットルバルブ上流圧力Ｐａと値Φ（Ｐｍ／Ｐａ）との関係を規定した前記テーブルと、後述するインテークマニホールドモデルＭ８が既に計算している最新の吸気管圧力ＰｍＲ、及び実際のスロットルバルブ上流圧力Ｐａとに基づいて求める。以上により、現時点におけるスロットル通過空気量ｍｔＴＨＲが求められる。
【００６２】
（６）エアフローメータモデルＭ６
エアフローメータモデルＭ６は、スロットル通過空気量が所定の量αである場合に、エアフローメータ６１が出力するであろう値を推定し、この推定値に基づいてスロットル通過空気量ｍｔＲを推定するモデルである。この場合、上記所定の量αは、スロットルモデルＭ５が推定したスロットル通過空気量ｍｔＴＨＲである。
【００６３】
エアフローメータモデルＭ６は、先ず、スロットル通過空気量ｍｔＴＨＲに対する完全放熱量Ｗ１，Ｗ２を、同完全放熱量Ｗ１，Ｗ２とスロットル通過空気量ｍｔとの関係を規定するテーブルと、前記求められたスロットル通過空気量ｍｔＴＨＲとに基づいて求める。完全放熱量Ｗ１、及び完全放熱量Ｗ２は、図４に示した熱線計量部６１ａのボビン部６１ａ１、及び同熱線計量部６１ａのサポート部６１ａ２にそれぞれ対応した放熱遅れを含まない放熱量である。
【００６４】
次に、エアフローメータモデルＭ６は、ボビン部６１ａ１、及びサポート部６１ａ２にそれぞれ対応する放熱量であり、完全放熱量Ｗ１，Ｗ２に対してそれぞれ一次遅れの特性を有する応答遅れを含む放熱量（応答放熱量）ｗ１，ｗ２を下記数１１及び下記数１２にしたがって求める。数１１，数１２における添え字ｉは今回の演算値、添え字ｉ−１は前回の演算値を表し、Δｔは前回の演算値を求めてから今回の演算値を求めるまでの時間である。
【００６５】
【数１１】
ｗ１ｉ＝Δｔ・（Ｗ１ｉ−ｗ１ｉ−１）／τ１＋ｗ１ｉ−１
【００６６】
【数１２】
ｗ２ｉ＝Δｔ・（Ｗ２ｉ−ｗ２ｉ−１）／τ２＋ｗ２ｉ−１
【００６７】
上記数１１，数１２において、τ１、及びτ２は、ボビン部６１ａ１、及びサポート部６１ａ２にそれぞれ対応する上記一次遅れ特性の時定数であり、下記数１３及び下記数１４により求められる。数１３，数１４中の値ｋ１０，ｋ２０、及び値ｍ１，ｍ２は、実験的に求められた値である。また、値ｕはエアフローメータ６１の熱線計量部６１ａにバイパスされた単位断面積当たりの通過空気量であり、図５に示したエアフローメータ６１の出力電圧Ｖｇと実測された吸入空気量ｍｔＡＦＭとの関係を規定するＶｇ−ｍｔＡＦＭ変換テーブルと、エアフローメータ６１の実際の出力電圧Ｖｇとに基づいて求められた吸入空気量ｍｔＡＦＭを、前記熱線計量部６１ａのバイパス流路断面積Ｓで除した値（ｍｔＡＦＭ／Ｓ）である。
【００６８】
【数１３】
τ１＝ｋ１０・ｕｍ１
【００６９】
【数１４】
τ２＝ｋ２０・ｕｍ２
【００７０】
そして、エアフローメータモデルＭ６は、応答放熱量ｗ１，ｗ２の和（ｗ１＋ｗ２）とエアフローメータ６１が出力するであろう値に基づくスロットル通過空気量ｍｔＲとの関係を規定した図９に示したテーブルと、上記数１１〜数１４により求められた応答放熱量ｗ１，ｗ２の和（ｗ１＋ｗ２）とに基づいて、現時点でエアフローメータ６１が出力するであろう値に基づくスロットル通過空気量ｍｔＲを求める。
【００７１】
（７）吸気弁モデルＭ７
吸気弁モデルＭ７は、上記吸気弁モデルＭ３と同様に、上記数７にしたがって現時点における筒内吸入空気流量ｍｃＲを求める。但し、吸気弁モデルＭ７は、後述するインテークマニホールドモデルＭ８により既に求めらている現時点の吸気管圧力ＰｍＲ、及び現時点の吸気管内温度ＴｍＲを、上記数７における吸気管圧力Ｐｍ、及び吸気管内温度Ｔｍに適用する等、必要なパラメータを全て現時点のものとして数７の計算を行う。
【００７２】
（８）インテークマニホールドモデルＭ８
インテークマニホールドモデルＭ８は、インテークマニホールドモデルＭ４と同様に、上記数８，数９を用いて現時点における吸気管圧力Ｐｍを求める。但し、インテークマニホールドモデルＭ８は、上記エアフローメータモデルＭ６により求められたスロットル通過空気量ｍｔＲ、及び上記吸気弁モデルＭ７により求められた現時点における筒内吸入空気流量ｍｃＲを、それぞれ数８，数９におけるスロットル通過空気量ｍｔ、及び筒内吸入空気流量ｍｃとして使用する。このインテークマニホールドモデルＭ８により推定されたＰｍが、前記エアフローメータ６１が現時点で出力するであろう値に基づく吸気管圧力Ｐｍ２となる。
【００７３】
（Ｐｍ３の求め方）
上記エアフローメータ６１の現時点における実際の出力電圧Ｖｇに基づく吸気管圧力Ｐｍ３は、上記インテークマニホールドモデルＭ４，Ｍ８と同じモデルであるインテークマニホールドモデルＭ９により求められる。
【００７４】
（９）インテークマニホールドモデルＭ９
具体的に述べると、インテークマニホールドモデルＭ９は、エアフローメータ６１の出力電圧Ｖｇと図５に示したＶｇ−ｍｔＡＦＭ変換テーブルとにより求められる現時点の実測された吸入空気量ｍｔＡＦＭを上記数８，数９におけるスロットル通過空気量ｍｔとして使用するとともに、上記吸気弁モデルＭ７により求められた現時点での筒内吸入空気流量ｍｃＲを同数８，数９の筒内吸入空気流量ｍｃとして使用し、吸気管圧力Ｐｍを求める。このインテークマニホールドモデルＭ９により推定された吸気管圧力Ｐｍが、エアフローメータ６１の現時点における実際の出力電圧Ｖｇに基づく吸気管圧力Ｐｍ３となる。以上により、吸気管圧力Ｐｍ１〜Ｐｍ３が求められ、本制御装置は、上記数１にしたがって吸気弁閉時の吸気管圧力ＰＭＦＷＤを予測・推定するとともに、推定した吸気管圧力ＰＭＦＷＤに基いて一気筒当たりの吸入空気量相当値ＫＬＦＷＤを求める。
【００７５】
（各燃料成分の燃料付着量及び各燃料成分の筒内流入燃料量の推定方法の概要）次に、本制御装置が行う、燃料を構成する燃料成分毎の燃料付着量の推定方法、同燃料成分毎のシリンダ内（筒内）流入燃料量の推定方法の概要について説明する。図１０に概念的に示したように、インジェクタ３９から噴射された燃料は、その一部が吸気管４１の壁面部、及び図１０において図示を省略した吸気弁等からなる吸気通路構成部材に付着する。
【００７６】
このように吸気通路構成部材に付着した燃料は、先に説明したように、蒸留性状（沸点、飽和蒸気圧）の異なる多数の燃料成分から構成される混合物である。しかし、本実施形態の制御装置は、計算の簡便化によるＣＰＵ７１の計算負担の軽減等の観点から、この燃料を、以下に説明する低沸点成分と高沸点成分との２つの燃料成分に二分し、これら２つの燃料成分から構成されているものとして扱う。
【００７７】
図１１は、吸気管４１の温度（吸気管４１の壁面温度、吸気通路内の温度。この温度は冷却水温ＴＨＷで代用することができる。）が内燃機関１０が完全に暖機した状態において発生し得る所定温度Ｔ０に維持されていて、且つ内燃機関１０が定常状態にあるものとした場合における、吸気管圧力ＰＭＦＷＤと、吸気通路構成部材への低沸点成分の燃料付着量ｆｗｖ及び吸気通路構成部材への高沸点成分の燃料付着量ｆｗｐとの関係をそれぞれ示したグラフである。
【００７８】
図１１から理解できるように、低沸点成分は、飽和蒸気圧が比較的高いことから、前記所定温度Ｔ０において吸気管圧力ＰＭＦＷＤが内燃機関１０の運転中において発生し得る範囲内における或る所定圧力Ｐ０を下回ると吸気通路構成部材に付着していた燃料付着量ｆｗｖの燃料の全量が蒸発して燃料付着量ｆｗｖが「０」になるような、炭素数が少なくて沸点が低い燃料成分の群を総称したものである。一方、高沸点成分は、飽和蒸気圧が比較的低いことから、前記所定温度Ｔ０において吸気管圧力ＰＭＦＷＤが内燃機関１０の運転中において発生し得る範囲内にあるときには常に吸気通路構成部材に付着する燃料付着量ｆｗｐの燃料が存在し得るような、炭素数が比較的多くて沸点が比較的高い燃料成分の群を総称したものである。
【００７９】
このように、吸気通路構成部材に付着する低沸点成分の挙動は、吸気通路構成部材に付着する高沸点成分の挙動と大きく相違する。換言すると、一般に、燃料付着モデルで使用される燃料の付着率と残留率は、低沸点成分に対するものと高沸点成分に対するものとで大きく相違する。従って、本制御装置は、燃料付着モデルを、吸気通路構成部材に付着する低沸点成分についてのモデルと、吸気通路構成部材に付着する高沸点成分についてのモデルとで独立させ、付着率及び残留率をそれぞれに対して設定することにより、吸気通路構成部材への低沸点成分の燃料付着量ｆｗｖと吸気通路構成部材への高沸点成分の燃料付着量ｆｗｐとを分けて求める。
【００８０】
より具体的に述べると、特定の気筒に着目した図１２に示したように、ｆｉをインジェクタ３９から同特定気筒の一吸気行程に対して噴射される燃料噴射量、ｆｗｐを吸気通路構成部材にすでに付着している高沸点成分の燃料付着量、Ｐｐを吸気通路構成部材にすでに付着している高沸点成分のうち一吸気行程を経た後に同吸気通路構成部材に付着したまま残留している高沸点成分の割合（高沸点成分の残留率Ｐｐ）、Ｒｐをインジェクタ３９から噴射された前記燃料のうち吸気通路構成部材へ付着する高沸点成分の割合（高沸点成分の付着率Ｒｐ）、添え字ｋを今回の演算値（今回の吸気行程に対する値）、添え字ｋ＋１は次回の演算値（次回の吸気行程に対する値）とすると、今回噴射された燃料のうち吸気通路構成部材に新たに付着する高沸点成分の量はＲｐ・ｆｉ（ｋ）であり、吸気通路構成部材にすでに付着していた高沸点成分のうち同吸気通路構成部材に残留する高沸点成分の量はＰｐ・ｆｗ（ｋ）であるから、高沸点成分の燃料付着量ｆｗｐ（ｋ＋１）について下記数１５が成立する。下記数１５は、高沸点成分についての燃料付着モデルを記述したものである。
【００８１】
【数１５】
ｆｗｐ（ｋ＋１）＝Ｒｐ・ｆｉ（ｋ）＋Ｐｐ・ｆｗｐ（ｋ）
【００８２】
同様に、ｆｗｖを吸気通路構成部材にすでに付着している低沸点成分の燃料付着量、Ｐｖを吸気通路構成部材にすでに付着している低沸点成分のうち一吸気行程を経た後に同吸気通路構成部材に付着したまま残留している低沸点成分の割合（低沸点成分の残留率Ｐｖ）、Ｒｖをインジェクタ３９から噴射された前記燃料のうち吸気通路構成部材へ付着する低沸点成分の割合（低沸点成分の付着率Ｒｖ）とすると、今回噴射された燃料のうち吸気通路構成部材に新たに付着する低沸点成分の量はＲｖ・ｆｉ（ｋ）であり、吸気通路構成部材にすでに付着していた低沸点成分のうち同吸気通路構成部材に残留する低沸点成分の量はＰｖ・ｆｗ（ｋ）であるから、低沸点成分の燃料付着量ｆｗｖ（ｋ＋１）について下記数１６が成立する。下記数１６は、低沸点成分についての燃料付着モデルを記述したものである。上記数１５及び下記数１６の演算を行う手段が燃料付着量推定手段に相当する。
【００８３】
【数１６】
ｆｗｖ（ｋ＋１）＝Ｒｖ・ｆｉ（ｋ）＋Ｐｖ・ｆｗｖ（ｋ）
【００８４】
ここで、上記高沸点成分の付着率Ｒｐは、噴射された前記燃料中の高沸点成分の割合Ｋｐと噴射された高沸点成分のうち吸気通路構成部材へ付着する高沸点成分の割合Ｒ’ｐとの積であり、上記低沸点成分の付着率Ｒｖは、噴射された前記燃料中の低沸点成分の割合Ｋｖと噴射された低沸点成分のうち吸気通路構成部材へ付着する低沸点成分の割合Ｒ’ｖとの積である。前記燃料中の高沸点成分の割合Ｋｐと前記燃料中の低沸点成分の割合Ｋｖ（Ｋｐ＋Ｋｖ＝１）は使用される燃料自体に依存し、本実施例では、前記燃料中の高沸点成分の割合Ｋｐと前記燃料中の低沸点成分の割合Ｋｖとをそれぞれ一定値として扱う。
【００８５】
また、上記数１５又は数１６にて使用される上記高沸点成分の残留率Ｐｐ、低沸点成分の残留率Ｐｖ、高沸点成分の付着率Ｒｐ、及び低沸点成分の付着率Ｒｖは、ＳＣＶ開度センサ６５により検出されるＳＣＶ開度θｉｖ、クランクポジションセンサ６７の出力に基くエンジン回転速度Ｎｅ、水温センサ６８により検出される冷却水温ＴＨＷ、及び上記数１により計算されている最新の吸気弁閉時の吸気管圧力ＰＭＦＷＤの関数Ｆｐｐ、関数Ｆｐｖ、関数Ｆｒｐ、及び関数Ｆｒｖにそれぞれ基く所定のテーブルにより決定される。
【００８６】
図１３（ａ）は、前記所定のテーブルにおいて、ＳＣＶ開度θｉｖ、及びエンジン回転速度Ｎｅがそれぞれ所定の一定値に維持されていて、且つ冷却水温ＴＨＷが前記所定温度Ｔ０に維持されているものと仮定した場合における、吸気管圧力ＰＭＦＷＤと、高沸点成分の残留率Ｐｐ及び低沸点成分の残留率Ｐｖとの関係をそれぞれ示したグラフである。また、図１３（ｂ）は、図１３（ａ）において仮定した場合（条件）と同一の場合における、吸気管圧力ＰＭＦＷＤと、高沸点成分の付着率Ｒｐ及び低沸点成分の付着率Ｒｖとの関係をそれぞれ示したグラフである。
【００８７】
図１３（ａ）及び図１３（ｂ）に示したように、上記高沸点成分の残留率Ｐｐ、低沸点成分の残留率Ｐｖ、高沸点成分の付着率Ｒｐ、及び低沸点成分の付着率Ｒｖは全て、吸気管圧力ＰＭＦＷＤが高くなるほど大きくなる。これは、吸気管圧力ＰＭＦＷＤが高くなるほど高沸点成分及び低沸点成分の各々の沸点が高くなり、その結果、吸気通路構成部材に付着していた高沸点成分及び低沸点成分の蒸発が抑制され、且つ、インジェクタ３９から噴射された高沸点成分及び低沸点成分の霧化が抑制されることに基く。
【００８８】
また、吸気管圧力ＰＭＦＷＤが内燃機関１０の運転中において発生し得る範囲内にある状態においては、高沸点成分の残留率Ｐｐは低沸点成分の残留率Ｐｖよりも大きく、且つ、高沸点成分の付着率Ｒｐは低沸点成分の付着率Ｒｖよりも大きい。これは、同一条件下においては高沸点成分の沸点が低沸点成分の沸点よりも高く、その結果、高沸点成分の方が低沸点成分よりも、付着していた吸気通路構成部材からの蒸発が抑制され、且つ、インジェクタ３９から噴射されたときの霧化が抑制されることに基く。
【００８９】
また、低沸点成分の残留率Ｐｖ及び低沸点成分の付着率Ｒｖは共に、吸気管圧力ＰＭＦＷＤが前記所定圧力Ｐ０を下回ると「０」になる。これは、先に図１１を参照しつつ説明したように、冷却水温ＴＨＷが前記所定温度Ｔ０に維持されている場合、吸気管圧力ＰＭＦＷＤが前記所定圧力Ｐ０を下回ると低沸点成分が吸気通路構成部材に付着し得なくなることに基く。
【００９０】
次に、本制御装置は、上述したように高沸点成分についての燃料付着モデルを利用して求めた吸気通路構成部材への高沸点成分の燃料付着量ｆｗｐと低沸点成分についての燃料付着モデルを利用して求めた吸気通路構成部材への低沸点成分の燃料付着量ｆｗｖとに基いて、高沸点成分の筒内流入燃料量ｆｃｐと低沸点成分の筒内流入燃料量ｆｃｖとを分けて求める。
【００９１】
即ち、一吸気行程において、今回噴射された燃料噴射量ｆｉ（ｋ）の燃料のうち吸気通路構成部材に付着せず筒内に直接吸入される高沸点成分の量はＫｐ・（１−Ｒ’ｐ）・ｆｉ（ｋ）＝（Ｋｐ−Ｒｐ）・ｆｉ（ｋ）となり、吸気通路構成部材に付着していた高沸点成分のうち同吸気通路構成部材から離脱して筒内に吸入される高沸点成分の量は（１−Ｐｐ）・ｆｗｐ（ｋ）となるから、筒内に吸入される高沸点成分の筒内流入燃料量ｆｃｐ（ｋ）は下記数１７により表すことができる。
【００９２】
【数１７】
ｆｃｐ（ｋ）＝（１−Ｐｐ）・ｆｗｐ（ｋ）＋（Ｋｐ−Ｒｐ）・ｆｉ（ｋ）
【００９３】
同様に、一吸気行程において、今回噴射された燃料噴射量ｆｉ（ｋ）の燃料のうち吸気通路構成部材に付着せず筒内に直接吸入される低沸点成分の量はＫｖ・（１−Ｒ’ｖ）・ｆｉ（ｋ）＝（Ｋｖ−Ｒｖ）・ｆｉ（ｋ）となり、吸気通路構成部材に付着していた低沸点成分のうち同吸気通路構成部材から離脱して筒内に吸入される低沸点成分の量は（１−Ｐｖ）・ｆｗｖ（ｋ）となるから、筒内に吸入される低沸点成分の筒内流入燃料量ｆｃｖ（ｋ）は下記数１８により表すことができる。上記数１７及び下記数１８の演算を行う手段が筒内流入燃料量推定手段に相当する。
【００９４】
【数１８】
ｆｃｖ（ｋ）＝（１−Ｐｖ）・ｆｗｖ（ｋ）＋（Ｋｖ−Ｒｖ）・ｆｉ（ｋ）
【００９５】
以上のように、本実施形態の制御装置は、高沸点成分についての燃料付着モデルと低沸点成分についての燃料付着モデルを利用して、吸気通路構成部材への高沸点成分の燃料付着量ｆｗｐと吸気通路構成部材への低沸点成分の燃料付着量ｆｗｖとを分けて求めるとともに、吸気通路構成部材への高沸点成分の燃料付着量ｆｗｐと吸気通路構成部材への低沸点成分の燃料付着量ｆｗｖとに基いて、高沸点成分の筒内流入燃料量ｆｃｐと低沸点成分の筒内流入燃料量ｆｃｖとを分けて求める。
【００９６】
（点火時期の決定方法の概要）
一般に、上記低沸点成分はオクタン価が低く、上記高沸点成分は同低沸点成分よりもオクタン価が高い。従って、筒内に流入する（筒内で実際に燃焼される）燃料全体としてのオクタン価（筒内オクタン価）は、筒内に流入する低沸点成分に対する筒内に流入する高沸点成分の比率に応じて変化し、同比率が大きくなるほど筒内オクタン価は高くなる。よって、前記比率の増加に応じて内燃機関１０の点火時期を進角側に移動することが好ましい。ここで、前記筒内に流入する低沸点成分に対する筒内に流入する高沸点成分の比率（以下、「成分比率」と称呼する。）ＲＡＴＩＯは、上記数１７及び数１８に基き、下記数１９により表すことができる。
【００９７】
【数１９】
ＲＡＴＩＯ＝ｆｃｐ（ｋ）／ｆｃｖ（ｋ）＝（（１−Ｐｐ）・ｆｗｐ（ｋ）＋（Ｋｐ−Ｒｐ）・ｆｉ（ｋ））／（（１−Ｐｖ）・ｆｗｖ（ｋ）＋（Ｋｖ−Ｒｖ）・ｆｉ（ｋ））
【００９８】
上記数１９において使用されている高沸点成分の残留率Ｐｐ、低沸点成分の残留率Ｐｖ、高沸点成分の付着率Ｒｐ、及び低沸点成分の付着率Ｒｖは、先に図１３等を参照しつつ説明したように吸気弁閉時の吸気管圧力ＰＭＦＷＤ等が変化すると互いに独立して変化する。従って、内燃機関１０の運転状態が過渡運転状態にあるときには、成分比率ＲＡＴＩＯの値は時々刻々と変化する。
【００９９】
一方、内燃機関１０の運転状態が定常状態にあるときには、吸気通路構成部材への高沸点成分の燃料付着量ｆｗｐも吸気通路構成部材への低沸点成分の燃料付着量ｆｗｖも時間的に変化しない状態となっている。従って、上記数１５においてｆｗｐ（ｋ＋１）をｆｗｐ（ｋ）と置き、上記数１６においてｆｗｖ（ｋ＋１）をｆｗｖ（ｋ）と置くと、それぞれ、下記数２０及び数２１を得ることができる。
【０１００】
【数２０】
ｆｗｐ（ｋ）＝（Ｒｐ／（１−Ｐｐ））・ｆｉ（ｋ）
【０１０１】
【数２１】
ｆｗｖ（ｋ）＝（Ｒｖ／（１−Ｐｖ））・ｆｉ（ｋ）
【０１０２】
上記数１９〜数２１から高沸点成分の燃料付着量ｆｗｐ（ｋ）及び低沸点成分の燃料付着量ｆｗｖ（ｋ）を消去して数１９を整理すると、内燃機関１０が定常状態にあるときにおける筒内に流入する低沸点成分に対する筒内に流入する高沸点成分の比率（以下、「定常状態時成分比率」と称呼する。）ＲＡＴＩＯｓｔｅａｄｙが得られ、定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙは、下記数２２により表すことができる。
【０１０３】
【数２２】
ＲＡＴＩＯｓｔｅａｄｙ＝Ｋｐ／Ｋｖ
【０１０４】
先に述べたように、前記燃料中の高沸点成分の割合Ｋｐと前記燃料中の低沸点成分の割合Ｋｖはそれぞれ一定値である。従って、前記成分比率ＲＡＴＩＯの値は、内燃機関１０が定常状態にある限りにおいて一定値である定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙとなる。なお、上記数２２にて示される定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙの値は、内燃機関１０に使用されている燃料自体の低沸点成分に対する高沸点成分の比率と等しい。従って、内燃機関１０が定常状態にあるときの筒内オクタン価は、燃料自体のオクタン価と等しいことになる。
【０１０５】
ここで、先に説明したように、本制御装置がＲＯＭ７２に格納している上述した基本点火時期テーブルにより求められる基本点火時期ＩＴｓｔｅａｄｙは、内燃機関１０が定常状態にあるときに最適な燃焼状態が得られるように決定されている。従って、基本点火時期ＩＴｓｔｅａｄｙは、筒内に流入する低沸点成分に対する筒内に流入する高沸点成分の比率が上記定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙと等しいときに最適な燃焼状態が得られるように決定されていることになる。
【０１０６】
よって、内燃機関１０が過渡状態にあって、上記数１９に基く成分比率ＲＡＴＩＯの値が上記定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙより大きいときは、内燃機関１０の点火時期を基本点火時期ＩＴｓｔｅａｄｙよりも進角側に補正し、成分比率ＲＡＴＩＯの値が定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙより小さいときは、同点火時期を基本点火時期ＩＴｓｔｅａｄｙよりも遅角側に補正することが好ましい。
【０１０７】
そこで、本実施形態の制御装置は、内燃機関１０の点火時期ＩＴを求めるため、先ず、成分比率ＲＡＴＩＯと定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙとの大小関係の程度を示す指標値として、下記数２３により比率比ＸＲを求める。
【０１０８】
【数２３】
ＸＲ＝ＲＡＴＩＯ／ＲＡＴＩＯｓｔｅａｄｙ
【０１０９】
次に、本制御装置は、比率比ＸＲと点火時期補正量ΔＩＴとの関係を規定する図１４にグラフにより示したテーブルと、同比率比ＸＲとに基いて点火時期補正量ΔＩＴを決定する。これにより、点火時期補正量ΔＩＴは、比率比ＸＲが「１」より大きいとき、比率比ＸＲの増加に応じて増大し（より進角側に補正する正の量となり）、比率比ＸＲが「１」より小さいとき、比率比ＸＲの減少に応じて減少する（より遅角側に補正する負の量となる）。換言すれば、点火時期補正量ΔＩＴは、成分比率ＲＡＴＩＯと定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙとの相違の程度に応じて決定される。
【０１１０】
そして、本制御装置は、下記数２４に基いて機関制御パラメータとしての内燃機関１０の点火時期ＩＴを求める。下記数２４により点火時期ＩＴを決定する手段が機関制御パラメータ決定手段を構成する。
【０１１１】
【数２４】
ＩＴ＝ＩＴｓｔｅａｄｙ＋ΔＩＴ
【０１１２】
このようにして、本制御装置は、上記数１７及び数１８により求められた高沸点成分の筒内流入燃料量ｆｃｐと低沸点成分の筒内流入燃料量ｆｃｖとに基いて基本点火時期ＩＴｓｔｅａｄｙを点火時期補正量ΔＩＴだけ補正することで、内燃機関１０の点火時期ＩＴを求める。
【０１１３】
（点火時期制御の一例）
図１５は、内燃機関１０が定常状態から急減速運転された場合における上記数１により計算される吸気管圧力ＰＭＦＷＤ、比率比ＸＲ、及び点火時期補正量ΔＩＴの変化の一例を示したタイムチャートである。先ず、このタイムチャートにおいては、時刻ｔ１になるまで、内燃機関１０は、吸気管４１の温度（冷却水温ＴＨＷ）が前記所定温度Ｔ０に近い温度になっていて、且つ、スロットルバルブ開度ＴＡが比較的大きい一定値に維持されていることで吸気管圧力ＰＭＦＷＤが前記所定圧力Ｐ０よりも大きい一定値となる、定常状態にある。
【０１１４】
この定常状態では、上述したように成分比率ＲＡＴＩＯの値は定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙと等しいので、図１５（ｂ）に示すように、比率比ＸＲは「１」に維持され、この結果、図１５（ｃ）に示すように、点火時期補正量ΔＩＴは「０」に維持されている。
【０１１５】
そして、時刻ｔ１になると、運転者が、アクセルペダル８２の操作量Ａｃｃｐが急激に所定量だけ減少した後その時点での同操作量Ａｃｃｐに維持されるようにアクセルペダル８２を操作したものとする。これにより、スロットルバルブ開度ＴＡも時刻ｔ１以降、急減した後に一定値に維持される。その結果、図１５（ａ）に示すように、吸気管圧力ＰＭＦＷＤも、時刻ｔ１以降、前記所定圧力よりも大きい一定値から急減を開始するとともに、時刻ｔ２になると、前記所定圧力Ｐ０よりも小さい一定値に維持されて、内燃機関１０は再び定常状態となる。
【０１１６】
このように、時刻ｔ１から時刻ｔ２までの過渡運転状態の間に吸気管圧力ＰＭＦＷＤが急激に減少すると、先に説明した図１３のグラフ、及び成分比率ＲＡＴＩＯを計算する上記数１９から理解できるように、上記低沸点成分の残留率Ｐｖ、及び低沸点成分の付着率Ｒｖが共に直ちに「０」になって、成分比率ＲＡＴＩＯが一時的に定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙよりも小さくなる。
【０１１７】
その結果、図１５（ｂ）に示すように、時刻ｔ１から時刻ｔ２の間、比率比ＸＲが一時的に「１」よりも小さくなるとともに、図１５（ｃ）に示すように、点火時期補正量ΔＩＴも一時的に負の値（点火時期を遅角側に補正する量）になる。そして、時刻ｔ２以降、内燃機関１０が定常状態に復帰したことにより、比率比ＸＲは「１」に維持されるとともに、点火時期補正量ΔＩＴは「０」に維持される。
【０１１８】
従って、図１５に示すタイムチャートにおいては、本制御装置は、点火時期ＩＴを、時刻ｔ１になるまでは基本点火時期ＩＴｓｔｅａｄｙになるように、時刻ｔ１以降時刻ｔ２までの過渡運転状態においては基本点火時期ＩＴｓｔｅａｄｙよりも図１５（ｃ）に示す点火時期補正量ΔＩＴだけ遅角側に補正した時期になるように、時刻ｔ２以降は再び基本点火時期ＩＴｓｔｅａｄｙになるように制御する。
【０１１９】
（実際の作動）
以下、上記内燃機関の制御装置の作動について、ＣＰＵ７１が実行するルーチン（プログラム）をフローチャートにより示した図１６〜図１９を参照しながら説明する。
【０１２０】
（スロットルバルブ制御）
ＣＰＵ７１は、図１６のスロットルバルブ制御ルーチンの処理を所定時間（２ｍｓｅｃ）の経過毎にステップ１６００から開始し、ステップ１６０５に進んでアクセルペダル操作量Ａｃｃｐを読み込む。次いで、ＣＰＵ７１はステップ１６１０に進み、同ステップ１６１０にて図８と同じテーブルを用いることにより上記読み込んだアクセルペダル操作量Ａｃｃｐに基づく暫定的な目標スロットルバルブ開度θｒ１を求める。
【０１２１】
次に、ＣＰＵ７１はステップ１６１５に進んで変数Ｉを「６４」に設定し、続くステップ１６２０にて記憶値θｒ（Ｉ）にθｒ（Ｉ−２）の値を格納する。現時点では、変数Ｉは「６４」であるから、記憶値θｒ（６４）に記憶値θｒ（６２）の値が格納される。次いで、ＣＰＵ７１はステップ１６２５に進み、変数Ｉが「２」と等しくなったか否かを判定する。この場合、変数Ｉの値は「６４」であるから、ＣＰＵ７１はステップ１６２５にて「Ｎｏ」と判定してステップ１６３０に進み、同ステップ１６３０にて変数Ｉの値を「２」だけ減少し、その後上記ステップ１６２０に戻る。この結果、ステップ１６２０が実行されると、記憶値θｒ（６２）に記憶値θｒ（６０）の値が格納される。このような処理は、変数Ｉの値が「２」となるまで繰り返し実行される。
【０１２２】
その後、ステップ１６３０の処理が繰り返されて変数Ｉの値が「２」となると、ＣＰＵ７１はステップ１６２５にて「Ｙｅｓ」と判定してステップ１６３５に進み、同ステップ１６３５にて前記ステップ１６１０にて求めた現時点における暫定的な目標スロットルバルブ開度θｒ１を記憶値θｒ（０）に格納する。以上により、現時点からＩｍｓｅｃ前（０ｍｓｅｃ≦Ｉｍｓｅｃ≦６４ｍｓｅｃ）の暫定的な目標スロットルバルブ開度θｒ（Ｉ）（Ｉ＝６４，６２，・・・，４，２，０）がＲＡＭ７３内に記憶されることになる。
【０１２３】
次に、ＣＰＵ７１はステップ１６４０に進み、同ステップ１６４０にて記憶値θｒ（６４）を最終的な目標スロットルバルブ開度θｒとして設定し、続くステップ１６４５にて実際のスロットルバルブ開度が目標スロットルバルブ開度θｒと等しくなるように、スロットルバルブアクチュエータ４３ａに対し駆動信号を出力し、その後ステップ１６９５にて本ルーチンを一旦終了する。
【０１２４】
以降においても、上記ルーチンの処理は２ｍｓｅｃの経過毎に実行される。この結果、実際のスロットルバルブ開度が、６４ｍｓｅｃ前のアクセルペダル操作量Ａｃｃｐに基づく目標スロットルバルブ開度θｒと等しくなるように制御される。これにより、上記電子制御スロットルモデルＭ１による吸気弁閉時のスロットルバルブ開度ＴＡＳの推定が可能となる。
【０１２５】
（吸気弁開閉タイミング制御、及びＳＣＶ開度制御）
ＣＰＵ７１は、図１７の吸気弁開閉タイミング・ＳＣＶ開度制御ルーチンを所定時間（例えば、２ｍｓｅｃ）の経過毎にステップ１７００から開始し、ステップ１７０５に進んクランクポジションセンサ６７の出力に基づくエンジン回転速度Ｎｅを読み込むとともに、ステップ１７１０にて前述した筒内吸入空気量に相当する値ＫＬＦＷＤ（即ち、エンジン負荷）を読込む。なお、筒内吸入空気量に相当する値ＫＬＦＷＤは、所定時間毎に繰り返し実行される前述した吸入空気量相当値の決定方法（図７に示したモデル）に従う図示しないルーチンにより求められている。
【０１２６】
次に、ＣＰＵ７１はステップ１７１５に進み、同ステップ１７１５内に示したテーブルと上記読み込んだエンジン回転速度Ｎｅ及び筒内吸入空気量相当値ＫＬＦＷＤとに基づいて吸気弁の開閉タイミング（進角量）ＶＴを決定し、続くステップ１７２０にて実際の進角量が前記決定した進角量ＶＴとなるように、アクチュエータ３３ａに駆動信号を出力する。なお、ステップ１７１５に示したテーブルにおいては、ＶＴ１，ＶＴ２，ＶＴ３の順に進角量が大きくなるように設定されている。
【０１２７】
次いで、ＣＰＵ７１はステップ１７２５に進み、同ステップ１７２５内に示したテーブルと上記読み込んだエンジン回転速度Ｎｅ及び筒内吸入空気量相当値ＫＬＦＷＤとに基づいて目標ＳＣＶ開度θｉｖｒを決定し、続くステップ１７３０にて実際のＳＣＶ開度が前記決定した目標ＳＣＶ開度θｉｖｒとなるように、アクチュエータ４４ａに駆動信号を出力した後、ステップ１７９５に進んで本ルーチンを一旦終了する。なお、ステップ１７２５に示したテーブルにおいては、θ１，θ２，θ３の順に値が大きくなるように設定されている。
【０１２８】
以降においても、上記処理は２ｍｓｅｃの経過毎に実行される。この結果、実際の吸気弁開閉タイミングの進角量と実際のＳＣＶ開度が、エンジン回転速度Ｎｅと筒内吸入空気量相当値ＫＬＦＷＤに応じた値に変更される。
【０１２９】
（燃料付着量の推定、及び燃料噴射制御）
ＣＰＵ７１は、特定気筒のクランク角が、その気筒の吸気上死点から所定クランク角度だけ前の角度（例えば、ＢＴＤＣ９０°）になると、図１８の燃料噴射制御ルーチンの処理をステップ１８００から開始してステップ１８０５に進み、同ステップ１８０５にて、図７に示したモデルに従って別途計算されている上記吸気弁閉時の筒内吸入空気量に相当する値ＫＬＦＷＤと、下記数２５とに基いて今回の要求燃料噴射量ｆｃ（ｋ）を算出する。下記数２５において、Ｋは設定空燃比に応じて変化する係数である。
【０１３０】
【数２５】
ｆｃ（ｋ）＝Ｋ・ＫＬＦＷＤ
【０１３１】
次いで、ＣＰＵ７１はステップ１８１０に進んで、ＳＣＶ開度センサ６５が検出するＳＣＶ開度θｉｖ、クランクポジションセンサ６７の出力に基づくエンジン回転速度Ｎｅ、水温センサ６８が検出する冷却水温ＴＨＷ、及び上記筒内吸入空気量相当値ＫＬＦＷＤを推定する際に求めた吸気管圧力ＰＭＦＷＤ等のパラメータ（以下、このパラメータを「引数パラメータ」と称呼する。）を読み込む。
【０１３２】
次に、ＣＰＵ７１はステップ１８１５に進んで、上記引数パラメータと、高沸点成分の付着率Ｒｐ、低沸点成分の付着率Ｒｖ、高沸点成分の残留率Ｐｐ、及び低沸点成分の残留率Ｐｖとの関係をそれぞれ規定する予めＲＯＭ７２に記憶したテーブルと、上記ステップ１８１０にて読み込んだ引数パラメータとに基づき、現時点での高沸点成分の付着率Ｒｐ、低沸点成分の付着率Ｒｖ、高沸点成分の残留率Ｐｐ、及び低沸点成分の残留率Ｐｖをそれぞれ決定する。
【０１３３】
次いで、ＣＰＵ７１はステップ１８２０に進み、下記数２６を変形した得た同ステップ１８２０内に記載した式、上記ステップ１８０５にて求めた要求燃料噴射量ｆｃ（ｋ）、及び上記ステップ１８１５にて決定した付着率Ｒｐ、付着率Ｒｖ、残留率Ｐｐ、及び残留率Ｐｖに基づいて今回の燃料噴射量ｆｉ（ｋ）を算出する。
【０１３４】
【数２６】
ｆｃ（ｋ）＝（１−Ｐｐ）・ｆｗｐ（ｋ）＋（１−Ｐｖ）・ｆｗｖ（ｋ）＋（１−Ｒｐ−Ｒｖ）・ｆｉ（ｋ）
【０１３５】
上記数２６は、筒内流入燃料量（要求燃料噴射量）ｆｃ（ｋ）が高沸点成分の筒内流入燃料量ｆｃｐ（ｋ）と低沸点成分の筒内流入燃料量ｆｃｖ（ｋ）の和であることに基く下記数２７と、上記数１７及び数１８とから、高沸点成分の筒内流入燃料量ｆｃｐ（ｋ）と低沸点成分の筒内流入燃料量ｆｃｖ（ｋ）を消去して数２７を整理し、同整理した後の式に、前記燃料中の高沸点成分の割合Ｋｐと前記燃料中の低沸点成分の割合Ｋｖとの関係を示す下記数２８を適用することにより得ることができる。
【０１３６】
【数２７】
ｆｃ（ｋ）＝ｆｃｐ（ｋ）＋ｆｃｖ（ｋ）
【０１３７】
【数２８】
Ｋｐ＋Ｋｖ＝１
【０１３８】
次に、ＣＰＵ７１はステップ１８２５に進み、上記数１５にしたがって高沸点成分の燃料付着量ｆｗｐ（ｋ＋１）を求めるとともに、続くステップ１８３０にて上記数１６にしたがって低沸点成分の燃料付着量ｆｗｖ（ｋ＋１）を求める。
【０１３９】
そして、ＣＰＵ７１はステップ１８３５に進んで、上記ステップ１８２０にて決定した今回の燃料噴射量ｆｉ（ｋ）だけ燃料を噴射するように前記特定気筒に対するインジェクタ３９に駆動信号を送出し、ステップ１８９５に進んで本ルーチンを一旦終了する。
【０１４０】
以上により、上記特定の気筒に対する燃料噴射量が、高沸点成分及び低沸点成分の各々の燃料付着量、従って高沸点成分及び低沸点成分の各々の筒内流入燃料量に基づいて決定され、同燃料噴射量の燃料が同特定気筒に対するインジェクタ３９から噴射される。なお、ＣＰＵ７１は、他の気筒に対しても、図１８のルーチンと同様なルーチンを同様なタイミングで実行する。
【０１４１】
（点火時期制御）
ＣＰＵ７１は、特定気筒のクランク角が、その気筒の圧縮上死点（爆発上死点）から所定クランク角度だけ前の角度（例えば、ＢＴＤＣ９０°）になると、図１９の点火時期制御ルーチンの処理をステップ１９００から開始してステップ１９０５に進み、同ステップ１９０５にて、上記数１９に基く同ステップ１９０５内に記載した式、図１８のステップ１８２０にて求めた今回の燃料噴射量ｆｉ（ｋ）、図１８のステップ１８１５にて決定した高沸点成分の付着率Ｒｐ、低沸点成分の付着率Ｒｖ、高沸点成分の残留率Ｐｐ、及び低沸点成分の残留率Ｐｖ、並びに、前回本ルーチン実行時において後述するステップ１９３５にて更新（格納）されている高沸点成分の燃料付着量ｆｗｐ（ｋ）及び低沸点成分の燃料付着量ｆｗｖ（ｋ）に基づいて現時点での成分比率ＲＡＴＩＯを算出する。
【０１４２】
次に、ＣＰＵ７１はステップ１９１０に進んで、前記ステップ１９０５にて算出した成分比率ＲＡＴＩＯと前記定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙ（一定値）と、上記数２３に相当するステップ１９１０内に記載の式とに基き比率比ＸＲを求める。次いで、ＣＰＵ７１はステップ１９１５に進んで、図１４のテーブルと、前記ステップ１９１０にて求めた比率比ＸＲとに基き点火時期補正量ΔＩＴを算出する。
【０１４３】
続いて、ＣＰＵ７１はステップ１９２０に進み、上記引数パラメータと、基本点火時期ＩＴｓｔｅａｄｙとの関係を規定する予めＲＯＭ７２に記憶した基本点火時期テーブルと、図１８のステップ１８１０にて読み込んだ引数パラメータとに基づき、現時点での基本点火時期ＩＴｓｔｅａｄｙを決定する。
【０１４４】
次に、ＣＰＵ７１はステップ１９２５に進んで、前記ステップ１９１５にて算出した点火時期補正量ΔＩＴと、前記ステップ１９２０にて決定した基本点火時期ＩＴｓｔｅａｄｙと、上記数２４に相当するステップ１９２５内に記載の式とに基き現時点で最適な点火時期ＩＴを算出する。
【０１４５】
そして、ＣＰＵ７１はステップ１９３０に進み、前記ステップ１９２５にて算出した点火時期ＩＴにて上記特定気筒に対する点火プラグ３７が点火するようにイグナイタ３８に駆動信号を出力した後、続くステップ１９３５にて、図１８のステップ１８２５にて求められている高沸点成分の燃料付着量ｆｗｐ（ｋ＋１）、及び図１８のステップ１８３０にて求められている低沸点成分の燃料付着量ｆｗｖ（ｋ＋１）を、次回の演算のために、高沸点成分の燃料付着量ｆｗｐ（ｋ）、及び低沸点成分の燃料付着量ｆｗｖ（ｋ）にそれぞれ置き換え、ステップ１９９５に進んで、本ルーチンを一旦終了する。
【０１４６】
以上により、上記特定の気筒に対する最適な点火時期が、高沸点成分及び低沸点成分の各々の筒内流入燃料量に基づいて決定され、同特定気筒に対する点火プラグ３７が前記決定された点火時期に点火する。なお、ＣＰＵ７１は、他の気筒に対しても、図１９のルーチンと同様なルーチンを同様なタイミングで実行する。
【０１４７】
以上説明したように、本発明による内燃機関の制御装置の実施形態によれば、吸気通路構成部材に対する付着挙動が互いに相違する高沸点成分及び低沸点成分の各々の燃料付着量ｆｗｐ（ｋ），ｆｗｖ（ｋ）を区別して推定することにより、各燃料成分の吸気通路構成部材への燃料付着量ｆｗｐ（ｋ），ｆｗｖ（ｋ）がそれぞれ精度良く推定される。従って、前記推定された各燃料成分の燃料付着量ｆｗｐ（ｋ），ｆｗｖ（ｋ）に基いて燃焼室２５内に流入する同各燃料成分の筒内流入燃料量ｆｃｐ（ｋ），ｆｃｖ（ｋ）が同燃料成分毎に区別して推定できるので、各燃料成分の筒内流入燃料量ｆｃｐ（ｋ），ｆｃｖ（ｋ）もそれぞれ精度良く推定される。この結果、内燃機関１０が過渡運転状態にあるときでも、時々刻々と変化する筒内に流入する燃料の各燃料成分間の成分比率ＲＡＴＩＯを精度良く求めることができる。そして、定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙ（一定値）に対する成分比率ＲＡＴＩＯの比率である比率比ＸＲに基き点火時期補正量ΔＩＴを求め、定常運転時にて予め適合されている基本点火時期ＩＴｓｔｅａｄｙを点火時期補正量ΔＩＴだけ補正することで点火時期ＩＴが求められる。従って、内燃機関１０が過渡運転状態にあるときでも、その運転状態を良好な状態に維持できるように点火時期ＩＴを決定することができた。
【０１４８】
本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記実施形態では、使用される燃料を高沸点成分と低沸点成分の２つの燃料成分に分けているが、使用される燃料を蒸留性状（沸点、飽和蒸気圧）が異なる３つ以上の燃料成分に分けてもよい。
【０１４９】
また、上記実施形態では、成分比率ＲＡＴＩＯに対する定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙの比率である比率比ＸＲに基き点火時期補正量ΔＩＴを求めているが、成分比率ＲＡＴＩＯと定常状態時成分比率ＲＡＴＩＯｓｔｅａｄｙとの差に基き点火時期補正量ΔＩＴを求めてもよい。
【図面の簡単な説明】
【図１】本発明による内燃機関の制御装置を火花点火式多気筒内燃機関に適用したシステムの概略構成図である。
【図２】図１に示した特定の気筒の燃焼室、及び同燃焼室の近傍部分を概念的に示した平面図である。
【図３】図１に示したエアフローメータの概略斜視図である。
【図４】図３に示したエアフローメータの熱線計量部の拡大斜視図である。
【図５】図１に示したＣＰＵが参照するエアフローメータの出力と吸入空気量（吸入空気流量）との関係を規定したテーブルを表したグラフである。
【図６】吸気弁閉時の吸気管圧力を予測する方法を説明するために、スロットルバルブ開度の変化と各種のモデルにより計算される吸気管圧力の変化を示したタイムチャートである。
【図７】図１に示した制御装置が吸気弁閉時の筒内吸入空気量に相当する値を推定するために採用した各種モデルの接続関係を示す機能ブロック図である。
【図８】図１に示したＣＰＵが参照するアクセルペダル操作量と目標スロットルバルブ開度との関係を規定したテーブルを表したグラフである。
【図９】図１に示したＣＰＵが参照する応答放熱量の和とエアフローメータが出力するであろう値に基づくスロットル通過空気量との関係を規定したテーブルを表したグラフである。
【図１０】図１に示した制御装置による燃料付着量の推定方法を説明するために、インジェクタから噴射された燃料が吸気通路に付着する様子を概念的に示した図である。
【図１１】図１に示した吸気管内の圧力と、高沸点成分と低沸点成分の各々の吸気通路構成部材への燃料付着量との関係を説明するための図である。
【図１２】図１に示したインジェクタから噴射された燃料量と、高沸点成分と低沸点成分の各々の吸気通路構成部材への燃料付着量、及び高沸点成分と低沸点成分の各々の筒内に流入する燃料量の関係を説明するための図である。
【図１３】図１に示したＣＰＵが参照する、吸気管圧力と、高沸点成分及び低沸点成分の各々の残留率及び付着率との関係を規定したテーブルを表したグラフである。
【図１４】図１に示したＣＰＵが参照する比率比と点火時期補正量との関係を規定したテーブルを表したグラフである。
【図１５】内燃機関が定常状態から急減速運転された場合における吸気管圧力、比率比、及び点火時期補正量の変化の一例を示したタイムチャートである。
【図１６】図１に示したＣＰＵが実行するスロットルバルブ開度を制御するためのルーチンを示したフローチャートである。
【図１７】図１に示したＣＰＵが実行する吸気弁開閉タイミング、及びＳＣＶ開度を制御するためのルーチンを示したフローチャートである。
【図１８】図１に示したＣＰＵが実行する燃料噴射量を制御するためのルーチンを示したフローチャートである。
【図１９】図１に示したＣＰＵが実行する点火時期を制御するためのルーチンを示したフローチャートである。
【符号の説明】
１０…火花点火式多気筒内燃機関、２０…シリンダブロック部（エンジン本体部）、２５…燃焼室、３１…吸気ポート、３２…吸気弁、３９…インジェクタ、４１…吸気管、４３…スロットルバルブ、７０…電気制御装置、７１…ＣＰＵ。
【図１】
2004144031.tif 000004
【図２】
2004144031.tif 000005
【図３】
2004144031.tif 000006
【図４】
2004144031.tif 000007
【図５】
2004144031.tif 000008
【図６】
2004144031.tif 000009
【図７】
2004144031.tif 000010
【図８】
2004144031.tif 000011
【図９】
2004144031.tif 000012
【図１０】
2004144031.tif 000013
【図１１】
2004144031.tif 000014
【図１２】
2004144031.tif 000015
【図１３】
2004144031.tif 000016
【図１４】
2004144031.tif 000017
【図１５】
2004144031.tif 000018
【図１６】
2004144031.tif 000019
【図１７】
2004144031.tif 000020
【図１８】
2004144031.tif 000021
【図１９】
2004144031.tif 000022