JP2004339958A - 内燃機関の筒内吸入新気量推定装置 - Google Patents

Abstract

【課題】バルブオーバーラップ時に吸気通路に吹き返された既燃ガスをシリンダ内に再吸入するときの吸気弁通過ガス流量を精度良く求め、これにより、筒内吸入新気量を精度良く求めること。
【解決手段】電気制御装置８０は、吸気弁３２を通過してシリンダ２１内に吸入される吸気弁通過ガスの流量を絞り部を通過する流体の流量を表す式に基づいて計算する。電気制御装置は、吸気弁通過ガス流量の積算値Ｍｉが０より大きく、シリンダ内に吸入されているガスが新気であると判定すると、同新気の温度を用いて吸気弁通過ガス流量を計算する。また、積算値Ｍｉが０より小さく、シリンダ内に吸入されているガスが吸気通路に吹き返された既燃ガスであると判定すると、排気ガス温度を用いて吸気弁通過ガス流量を計算する。電気制御装置は、このように求められる吸気弁通過ガス流量を積分して筒内吸入新気量を推定する。
【選択図】 図１０

Description

【０００１】
【発明の属する技術分野】
本発明は、内燃機関のシリンダ内に吸入される新気の量（筒内吸入新気量）等を推定する内燃機関の吸入新気量推定装置に関する。
【０００２】
【従来の技術】
内燃機関により燃焼される混合気の空燃比を所望の値に制御したり、点火時期を適切に制御するためには、同内燃機関のシリンダ内に吸入される新気の量（以下、「筒内吸入新気量Ｍａ」と称呼する。）を精度良く求める必要がある。新気とは大気のことであり燃焼ガスを含まない。通常、この筒内吸入新気量Ｍａは、内燃機関の吸気通路に設けられた空気流量センサの出力値により推定される。ところが、スロットルバルブ開度が急激に変化する過渡運転状態においては、空気流量センサの出力値の挙動と実際の吸入新気量の挙動とが一致しなくなるため、空気流量センサの出力値に基いて筒内吸入新気量Ｍａをあらゆる運転条件下で精度良く求めることは一般に困難である。
【０００３】
そこで、近年においては、エネルギー保存則や運動量保存則等の物理法則に基づいて得られた式により吸気通路における空気の挙動を表すモデルを構築し、このモデルを用いることにより、筒内吸入新気量Ｍａに応じた値を精度良く推定する種々の試みがなされている（例えば、特許文献１参照。）。
【０００４】
【特許文献１】
特開２００１−４１０９５号公報（第４−１２頁、図３、図４）
【０００５】
本出願人は、このようなモデルの一つとして、質量保存則及びエネルギー保存則に基いてシリンダ（シリンダ内の空気）に関するモデル（シリンダモデル）を構築し、このモデルを用いて筒内吸入新気量Ｍａを推定する技術を検討している。
【０００６】
より具体的に述べると、上記シリンダモデルは、下記（１）式及び下記（２）式を基本として構築されている。
【０００７】
【数１】

【０００８】
【数２】

【０００９】
（１）式は質量保存則に基づく式である。（１）式は「新気及び既燃ガスからなるシリンダ内のガスの量である筒内ガス量Ｍｃの時間的変化量（ｄＭｃ／ｄｔ）」は、「吸気弁の周囲を通過してシリンダ内に吸入されるガス（以下、「吸気弁通過ガス」と称呼する。）の流量である吸気弁通過ガス流量ｍｉ」と「排気弁の周囲を通過してシリンダ内に吸入されるガス（以下、「排気弁通過ガス」と称呼する。）の流量である排気弁通過ガス流量ｍｅ」との和に等しいことを表している。なお、この場合、吸気弁通過ガス流量ｍｉ及び排気弁通過ガス流量ｍｅの何れも、シリンダ内にガスが流れ込む場合に正の値となるように規定されている。
【００１０】
（２）式は、エネルギー保存則に基づく式である。（２）式は「シリンダ内のガスのエネルギーＥの時間的変化量（ｄＥ／ｄｔ）」が、「吸気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｉ’」、「排気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｅ’」及び「シリンダ壁面からシリンダ内のガスに単位時間あたりに与えられる熱量（以下、「伝達熱流量」と称呼する。）Ｑｗ’」の和から「シリンダ内のガスが単位時間あたりにピストンに対して行う仕事Ｗ’」を減じた値に等しいことを表している。
【００１１】
一方、吸気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｉ’は、下記（３）式に示したように、吸気弁通過ガス流量ｍｉ及び吸気弁通過ガス温度（吸気ガス温度）Ｔｉに比例する。ここで、Ｃｐｉは吸気弁通過ガスの定圧比熱である。また、排気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｅ’は、下記（４）式に示したように、排気弁通過ガス流量ｍｅ及び排気弁通過ガス温度（排気ガス温度）Ｔｅに比例する。ここで、Ｃｐｅは排気弁通過ガスの定圧比熱である。
【００１２】
【数３】
Ｅｉ’＝ｍｉ・Ｃｐｉ・Ｔｉ …（３）
【００１３】
【数４】
Ｅｅ’＝ｍｅ・Ｃｐｅ・Ｔｅ …（４）
【００１４】
そこで、シリンダモデルは、（３）式及び（４）式の関係を（２）式に適用した式をエネルギー保存則に基づく式とし、その式と（１）式とを連立させることにより、シリンダ内のガスの圧力（以下、「筒内ガス圧力」と称呼する。）Ｐｃとシリンダ内のガスの温度（以下、「筒内ガス温度」と称呼する。）Ｔｃとを求めるようになっている。そして、このように求められた筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃは、下記（５）式乃至下記（８）式に適用され、これにより吸気弁通過ガス流量ｍｉ及び排気弁通過ガス流量ｍｅが更新される。なお、（５）式乃至（８）式は、絞り部を通過する流体の流量を求める一般的な式に基づく式である。Ｃｄｉ，Ｃｄｅは流量係数であり、Ａｉ，Ａｅは流路断面積である。
【００１５】
【数５】

【００１６】
【数６】

【００１７】
【数７】

【００１８】
【数８】

【００１９】
（５）式又は（６）式により更新された吸気弁通過ガス流量ｍｉ及び（７）式又は（８）式により更新された排気弁通過ガス流量ｍｅは、上記（１）乃至（４）式に適用され、新たな筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃが求められる。その新たな筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃは、再び上記（５）式乃至（８）式に適用され、更に新たな吸気弁通過ガス流量ｍｉ及び排気弁通過ガス流量ｍｅが求められる。
【００２０】
このように、吸気弁通過ガス流量ｍｉ、排気弁通過ガス流量ｍｅ、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃは順次更新されて行く。また、更新された吸気弁通過ガス流量ｍｉは、吸気弁開弁時から吸気弁閉弁時まで積算され、これにより筒内吸入新気量Ｍａが推定される。更新された排気弁通過ガス流量ｍｅは、吸気弁開弁時から排気弁閉弁時まで積算され、これにより筒内に残留する既燃ガス量（以下、「筒内残留ガス量」と称呼する。）Ｍｅが推定される。
【００２１】
【発明が解決しようとする課題】
一般に、内燃機関の吸気弁と排気弁は同時に開弁されている状態をとるときがある。即ち、図１６の（Ａ）における曲線ＬＩＮ及び曲線ＬＥＸにより吸気弁及び排気弁のリフト量をそれぞれ示したように、吸気弁はクランク角が吸気上死点（ＴＤＣ）前の所定クランク角θ０（時刻ｔ０）となったときに開弁する。排気弁はこの時刻ｔ０において依然として開弁していて、その後、クランク角が吸気上死点（ＴＤＣ）後の所定クランク角θ２（時刻ｔ２）となったときに閉弁する。この吸気弁と排気弁とが共に開弁している期間（時刻ｔ０〜ｔ２）は、バルブオーバーラップ時（ＶＯＬ）と称呼される。
【００２２】
かかるバルブオーバーラップ時の初期は、排気行程の後半であるから、燃焼室（シリンダ）の容積は減少している。従って、吸気弁上流ガス圧力Ｐｉは筒内ガス圧力Ｐｃ及び排気ガス圧力Ｐｅより小さいので、排気系（排気ポート及び排気管）に排出された既燃ガスが排気弁の周囲、シリンダ及び吸気弁の周囲を介して吸気通路（吸気系、即ち、スロットルバルブ下流の吸気通路を構成する吸気ポート及び吸気管）に吹き返される。この場合、吸気弁通過ガス流量ｍｉは（６）式により計算されて負の値となるから、図１６の（Ｂ）に示したように、吸気弁通過ガス流量ｍｉの積算値は時刻ｔ０から減少して負の値となる。
【００２３】
その後、クランク角が吸気上死点（ＴＤＣ）後の所定のクランク角θ１となると、排気弁は実質的に閉弁する（時刻ｔ１）。この時点は吸気行程の初期であり、燃焼室（シリンダ）の容積は増大しているから、吸気弁上流ガス圧力Ｐｉは筒内ガス圧力Ｐｃより大きくなる。従って、吸気通路に吹き返された既燃ガスは、吸気弁の周囲を介してシリンダ内に再度吸入され始める。この場合、吸気弁通過ガス流量ｍｉは（５）式により計算されて正の値となるから、図１６の（Ｂ）に示したように、吸気弁通過ガス流量ｍｉの積算値は時刻ｔ１から増大し始める。そして、時刻ｔ４にて、吸気通路に吹き返された既燃ガスが総べてシリンダ内に吸入されると、新気がシリンダ内に吸入され始める。吸気弁通過ガス流量ｍｉの積算値は、時刻ｔ４にて「０」となり、時刻ｔ４以降は新気がシリンダ内に吸入されるから、最終的には筒内吸入新気量Ｍａを表す値となる。
【００２４】
ところで、吸気通路に吹き返された既燃ガスは燃焼直後のガスであるから、その温度は吸気通路に存在している新気の温度よりも高い。しかしながら、前述した装置は、吸気通路に吹き返された既燃ガスがシリンダ内に再吸入されるときに使用する吸気弁通過ガス流量ｍｉを求めるための（３）式において、吸気弁通過ガス温度Ｔｉを吸気管内空気温度Ｔｍ（スロットルバルブ下流の新気の温度）と等しい温度に設定しているので、計算される吸気弁通過ガス流量ｍｉの絶対値が実際の吸気弁通過ガス流量の絶対値よりも大きくなる。その結果、図１６の（Ｂ）において、破線により示したように、計算上の新気の吸入開始時刻ｔ３（吸気弁通過ガス流量ｍｉの積算値が「０」に復帰する時点）が実際の新気の吸入開始時刻ｔ４よりも早い時刻となり、それ故に筒内吸入新気量Ｍａが実際の筒内吸入新気量と異なってしまうという問題がある。
【００２５】
【発明の概要】
本発明による内燃機関の筒内吸入新気量推定装置は、上記課題に対処するためになされたものであって、吸気弁の周囲を通過して吸気通路からシリンダ内に吸入されているガスが「新気であるか」又は「同シリンダ内から同吸気通路へ吹き返された既燃ガスであるか」を判定し、その判定結果に応じて適切なガスの状態パラメータを選択することにより吸気弁を通過するガスの流量（吸気弁通過ガス流量）を推定する。従って、吸気弁通過ガス流量が精度良く求められるので、同吸気弁通過ガス流量に基づいて推定される筒内吸入新気量も精度良く求められる。
【００２６】
より具体的に述べると、本発明による内燃機関の筒内吸入新気量推定装置は、内燃機関の吸気弁の周囲を通過して吸気通路からシリンダ内に吸入されているガスが新気であるか同シリンダ内から同吸気通路へ吹き返された既燃ガスであるかを判定する吸気弁通過ガス種判定手段と、前記シリンダ内に吸入されているガスが前記新気であると判定されているとき同新気の状態パラメータを使用して前記吸気弁の周囲を通過するガスの流量である吸気弁通過ガス流量を算出するとともに同シリンダ内に吸入されているガスが前記既燃ガスであると判定されているとき同既燃ガスの状態パラメータを使用して同吸気弁通過ガス流量を推定する吸気弁通過ガス流量算出手段と、前記算出された吸気弁通過ガス流量に基づいて前記シリンダ内に吸入される新気の量である筒内吸入新気量を推定する筒内吸入新気量推定手段とを備えている。状態パラメータとはガスの温度や比熱比を含む。
【００２７】
これによれば、吸気弁通過ガス流量算出手段は、吸気弁の周囲を通過するガスの流量（吸気弁通過ガス流量）をそのガスの状態パラメータを使用しながら算出する。このとき、シリンダ内に吸入されているガスが新気であると判定されていれば、そのガスの状態パラメータとして新気の状態パラメータが使用される。一方、シリンダ内に吸入されているガスが既燃ガスであると判定されていれば、そのガスの状態パラメータとして既燃ガスの状態パラメータが使用される。従って、シリンダに吸入されているガスの種類に応じた状態パラメータにより吸気弁通過ガス流量が精度良く求められるので、同吸気弁通過ガス流量に基づいて算出される筒内吸入新気量が精度良く求められる。
【００２８】
この場合、前記吸気弁通過ガス流量算出手段は、前記吸気通路から前記シリンダ内にガスが吸入されているときに前記吸気弁通過ガス流量を正及び負の何れか一方の符号付きガス流量として、前記シリンダ内から前記吸気通路にガスが吹き返されているときに前記吸気弁通過ガス流量を正及び負の何れか他方の符号付きガス流量として、所定時間の経過毎に算出するように構成され、前記吸気弁通過ガス種判定手段は、前記算出された吸気弁通過ガス流量の前記吸気弁の開弁時からの積算値が同積算値の初期値（同積算値の初期値が０である場合は０）より大きいか否かに応じて前記シリンダ内に吸入されているガスが前記新気であるか前記既燃ガスであるかを判定するように構成されることが好適である。
【００２９】
一般の内燃機関においては、吸気弁は排気行程の後半であって排気弁が閉弁する前に開弁する。このとき、吸気通路の圧力はシリンダ内の圧力より小さいので、排気通路に排出された排気ガス（既燃ガス）が排気弁の周囲、シリンダ及び吸気弁の周囲を介して吸気通路に吹き返す（逆流する）。この場合、吸気弁通過ガス流量算出手段は、吸気弁通過ガス流量を正及び負の何れか一方の符号を有する値として算出する。
【００３０】
その後、吸気行程に進んで排気弁が実質的に閉弁するタイミングになると、吸気通路の圧力はシリンダ内の圧力より大きくなる。従って、吸気通路内に吹き返された既燃ガスは吸気弁の周囲を介して再びシリンダ内に吸入される。この場合、吸気弁通過ガス流量算出手段は、吸気弁通過ガス流量を正及び負の何れか他方の符号を有する値として算出する。
【００３１】
従って、吸気弁通過ガス流量が精度良く算出されていれば、吸気弁が開弁されてからの吸気弁通過ガス流量の積算値は、吸気通路内に吹き返された既燃ガスの総べてがシリンダ内に吸入された時点で吸気弁が開弁されたときの積算の初期値に復帰する。かかる知見に基づき、上記筒内吸入新気量推定装置は、吸気弁通過ガス流量の積算値が同積算値の初期値より大きいか否かに応じて前記シリンダ内に吸入されているガスが前記新気であるか前記既燃ガスであるかを判定する。
【００３２】
即ち、説明の便宜上、吸気弁通過ガス流量算出手段は、ガスがシリンダ内から吸気通路に逆流している場合に負の符号を有する吸気弁通過ガス流量、ガスが吸気通路からシリンダ内に流入している順流の場合に正の符号を有する吸気弁通過ガス流量を算出すると仮定し、且つ、吸気弁通過ガス流量の積算開始時（吸気弁が開弁したとき）における初期値を０と仮定すると、吸気弁通過ガス流量の積算値は吸気弁の開弁時から負の値となって減少する。そして、吸気弁通過ガス流量の積算値は、それまでに吸気通路に吹き返された既燃ガスが再びシリンダ内に吸入されはじめたときに増大を開始するとともに、吸気通路に吹き返された既燃ガスの総べてがシリンダ内に吸入された時点で０に復帰し、その後、正の値となって増大する。
【００３３】
従って、吸気弁通過ガス流量の積算値が負の値であるときは既燃ガスが吸気通路へと吹き返されているか又はそれまでに吸気通路に吹き返された既燃ガスがシリンダへと吸入されていると判定することができる。また、吸気弁通過ガス流量の積算値が正の値であるときは、新気がシリンダへと吸入されていると判定することができる。吸気弁通過ガス種判定手段は、かかる知見に基づいてなされていて、吸気弁通過ガス流量の積算値が０（積算値の初期値）より大きいか否かを判定するという簡単な構成で、吸気弁通過ガスの種類を判定するのである。
【００３４】
この場合、前記筒内吸入新気量推定装置は、前記吸気通路内の新気の温度を取得する新気温取得手段と、前記内燃機関の排気通路内の排気ガスの温度を取得する排気ガス温取得手段と、を備えるとともに、前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の温度を使用するとともに前記既燃ガスの状態パラメータとして前記取得された排気ガスの温度を使用するように構成されていることが好適である。
【００３５】
シリンダから吸気通路に吹き返されたガスの温度は、新気の温度よりも排気ガス温度により近いので、上記構成により、吸気通路に吹き返された既燃ガスをシリンダに再吸入するときの吸気弁通過ガス流量がより精度良く算出され得る。
【００３６】
また、前記筒内吸入新気量推定装置は、前記吸気通路内の新気の温度を取得する新気温取得手段と、前記シリンダ内のガスの温度を取得する筒内ガス温取得手段と、を備えるとともに、前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の温度を使用するとともに前記既燃ガスの状態パラメータとして前記取得されたシリンダ内のガスの温度を使用するように構成されていることが好適である。
【００３７】
シリンダから吸気通路に吹き返されたガスは、排気通路からシリンダ内を通過して吸気通路に吹き返されたガスであるから、その温度は新気の温度よりもシリンダ内のガスの温度により近い。従って、上記構成により、吸気通路に吹き返された既燃ガスをシリンダに再吸入するときの吸気弁通過ガス流量がより精度良く算出され得る。
【００３８】
また、前記筒内吸入新気量推定装置は、前記吸気通路内の新気の温度を取得する新気温取得手段と、前記吸気弁の周囲を通過して前記シリンダから前記吸気通路へ吹き返された既燃ガスのエネルギーを取得するとともに同取得した既燃ガスのエネルギーに基づいて同既燃ガスの温度を取得する吹き返しガス温取得手段と、を備えるとともに、前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の温度を使用するとともに前記既燃ガスの状態パラメータとして前記取得された既燃ガスの温度を使用するように構成されていることが好適である。
【００３９】
これによれば、シリンダから吸気通路へ吹き返された既燃ガスのエネルギーが取得され、その取得された既燃ガスのエネルギーに基づいて同既燃ガスの温度が推定される。そして、この既燃ガスの温度は、吸気通路へ吹き返された既燃ガスがシリンダへと再吸気されているときに、吸気弁通過ガス流量を求めるための状態パラメータとして使用される。従って、シリンダから吸気通路に吹き返された既燃ガスがシリンダに再吸入されるときの吸気弁通過ガス流量がより精度良く算出され得る。
【００４０】
また、前記筒内吸入新気量推定装置は、前記吸気通路内の新気の比熱比を取得する新気比熱比取得手段と、前記排気通路内の排気ガスの比熱比を取得する排気ガス比熱比取得手段と、を備えるとともに、前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の比熱比を使用するとともに前記既燃ガスの状態パラメータとして前記取得された排気ガスの比熱比を使用するように構成されることが好適である。
【００４１】
これによれば、吸気通路からシリンダ内に吸入されるガスの比熱比が実際の比熱比により近似した値となるので、吸気弁通過ガス流量が一層精度よく推定され得る。
【００４２】
また、上記何れかの筒内吸入新気量推定装置において、前記吸気弁通過ガス流量算出手段は、前記吸気弁の周囲に形成されるガス通路を絞り部とみなしたときに成立する同絞り部を通過するガス流量を算出する式に基づいて前記吸気弁通過ガス流量を算出するように構成されることが好適である。
【００４３】
吸気弁が開弁した状態においては、吸気弁の開弁により形成されるガス通路は絞り部（オリフィス）とみなすことができるので、上記構成により吸気弁通過ガス流量がより精度良く推定される。
【００４４】
【発明の実施の形態】
以下、本発明による内燃機関の吸入新気量推定装置（この装置は、筒内吸入空気量推定装置、筒内吸入ガス量推定装置、筒内残留ガス量推定装置と称呼することもできる。）の第１実施形態について図面を参照しつつ説明する。この吸入新気量推定装置は燃料噴射量制御装置の一部である。図１は、係る燃料噴射量制御装置を火花点火式多気筒（４気筒）内燃機関１０に適用したシステムの概略構成を示している。なお、図１は、特定気筒の断面のみを示しているが、他の気筒も同様な構成を備えている。
【００４５】
内燃機関１０は、シリンダブロック、シリンダブロックロワーケース、及びオイルパン等を含むシリンダブロック部２０と、シリンダブロック部２０の上に固定されるシリンダヘッド部３０と、シリンダブロック部２０にガソリン混合気を供給するための吸気系統４０と、シリンダブロック部２０からの排ガスを外部に放出するための排気系統５０とを含んでいる。
【００４６】
シリンダブロック部２０は、シリンダ２１、ピストン２２、コンロッド２３、及びクランク軸２４を含んでいる。ピストン２２はシリンダ２１内を往復動し、ピストン２２の往復動がコンロッド２３を介してクランク軸２４に伝達され、これにより同クランク軸２４が回転するようになっている。シリンダ２１とピストン２２のヘッドは、シリンダヘッド部３０とともに燃焼室２５を形成している。
【００４７】
シリンダヘッド部３０は、燃焼室２５に連通した吸気ポート３１、吸気ポート３１を開閉する吸気弁３２、吸気弁３２を駆動するインテークカムシャフトを含むとともに同インテークカムシャフトの位相角及び同吸気弁３２の吸気弁リフト量（最大吸気弁リフト量）を連続的に変更し得る吸気弁制御装置３３、吸気弁制御装置３３のアクチュエータ３３ａ、燃焼室２５に連通した排気ポート３４、排気ポート３４を開閉する排気弁３５、排気弁３５を駆動するエキゾーストカムシャフト３６、点火プラグ３７、点火プラグ３７に与える高電圧を発生するイグニッションコイルを含むイグナイタ３８及び燃料を吸気ポート３１内に噴射するインジェクタ（燃料噴射手段）３９を備えている。
【００４８】
吸気系統４０は、吸気ポート３１に連通し同吸気ポート３１とともに吸気通路を形成するインテークマニホールドを含む吸気管４１、吸気管４１の端部に設けられたエアフィルタ４２、吸気管４１内にあって吸気通路の開口断面積を可変とするスロットルバルブ４３及びスワールコントロールバルブ（以下、「ＳＣＶ」と称呼する。）４４を備えている。スロットルバルブ４３は、ＤＣモータからなるスロットルバルブアクチュエータ４３ａにより吸気管４１内で回転駆動されるようになっている。ＳＣＶ４４は、前記スロットルバルブ４３よりも下流で前記インジェクタ３９よりも上流の位置にて前記吸気管４１に対し回動可能に支持されるとともに、ＤＣモータからなるＳＣＶアクチュエータ４４ａにより回転駆動されるようになっている。
【００４９】
排気系統５０は、排気ポート３４に連通したエキゾーストマニホールド５１、エキゾーストマニホールド５１に接続された排気管５２及び排気管５２に介装された触媒コンバータ（三元触媒装置）５３を備えている。
【００５０】
一方、このシステムは、熱線式エアフローメータ６１、吸気温センサ６２、大気圧センサ（スロットルバルブ上流圧力センサ）６３、スロットルポジションセンサ６４、ＳＣＶ開度センサ６５、カムポジションセンサ６６、吸気弁リフト量センサ６７、クランクポジションセンサ６８、水温センサ６９、空燃比センサ（Ｏ_２センサ）７０及びアクセル開度センサ７１を備えている。
【００５１】
エアフローメータ６１は、内燃機関１０の吸気通路内を流れる空気の流量（吸入空気流量）を実際に測定し、測定した吸入空気流量ＡＦＭを表す信号を出力するようになっている。吸気温センサ６２は、エアフローメータ６１内に備えられていて、吸入空気の温度を検出し、吸気温度Ｔａを表す信号を出力するようになっている。大気圧センサ６３は、スロットルバルブ４３の上流の圧力（即ち、大気圧）を検出し、スロットルバルブ上流圧力Ｐａを表す信号を出力するようになっている。スロットルポジションセンサ６４は、スロットルバルブ４３の開度（スロットルバルブ開度）を検出し、スロットルバルブ開度ＴＡを表す信号を出力するようになっている。ＳＣＶ開度センサ６５は、ＳＣＶ４４の開度を検出し、ＳＣＶ開度θｉｖを表す信号を出力するようになっている。
【００５２】
カムポジションセンサ６６は、インテークカムシャフトが９０°回転する毎に（即ち、クランク軸２４が１８０°回転する毎に）一つのパルスを有する信号（Ｇ２信号）を発生するようになっている。吸気弁リフト量センサ６７は、吸気弁３２のリフト量を検出し、吸気弁３２が全閉のとき「０」の値をとる吸気弁リフト量Ｌｉを表す信号を出力するようになっている。クランクポジションセンサ（エンジン回転速度センサ）６８は、クランク軸２４が１０°回転する毎に幅狭のパルスを有するとともに同クランク軸２４が３６０°回転する毎に幅広のパルスを有する信号を出力するようになっている。この信号は、エンジン回転速度ＮＥを表す。
【００５３】
水温センサ６９は、内燃機関１０の冷却水の温度を検出し、冷却水温ＴＨＷを表す信号を出力するようになっている。Ｏ_２センサ７０は、触媒コンバータ５３に流入する排ガス中の酸素濃度に応じた信号（排ガスの空燃比に応じた値）を出力するようになっている。アクセル開度センサ７１は、運転者によって操作されるアクセルペダルの操作量Ａｃｃｐを表す信号を出力するようになっている。
【００５４】
電気制御装置８０は、互いにバスで接続された、ＣＰＵ８１、ＲＯＭ８２、ＲＡＭ８３、バックアップＲＡＭ８４、及びインターフェース８５等からなるマイクロコンピュータである。ＲＯＭ８２は、テーブル（マップ）及び定数等を予め記憶するようになっている。ＲＡＭ８３は、ＣＰＵ８１の必要に応じてデータを一時的に格納するようになっている。バックアップＲＡＭ８４は、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するようになっている。インターフェース８５は、ＡＤコンバータを含み、前記センサ６１〜７１と接続され、ＣＰＵ８１にセンサ６１〜７１からの信号を供給するとともに、ＣＰＵ８１の指示に応じて吸気弁制御装置３３のアクチュエータ３３ａ、イグナイタ３８、インジェクタ３９、スロットルバルブアクチュエータ４３ａ及びＳＣＶアクチュエータ４４ａ等に駆動信号を送出するようになっている。
【００５５】
次に、上記のように構成された燃料噴射量制御装置による燃料噴射量の決定方法（筒内吸入新気量Ｍａの推定方法）について説明する。
【００５６】
（燃料噴射量ｆｃの決定方法・筒内吸入新気量の推定方法）
燃料噴射量制御装置は、吸気行程にある気筒の吸気弁３２が閉じる前に同気筒に対して燃料を噴射しなければならない。このため、燃料噴射量制御装置は、吸気弁３２が閉じた時点で（即ち、吸気弁閉時に）同気筒内に吸入されているであろう筒内吸入新気量Ｍａを吸気弁３２が閉弁する前に予測し、ｆｃ＝ｋλ・Ｍａなる式に基いて燃料噴射量（基本噴射量）ｆｃを決定する。ここで、ｋλは運転状態に応じて変化する設定空燃比に基づく係数である。そして、燃料噴射量制御装置は、吸気行程にある気筒のインジェクタ３９から、その気筒の吸気弁３２が閉弁するよりも前の時点で前記決定された燃料噴射量ｆｃの燃料を噴射する。
【００５７】
以下、燃料噴射量制御装置について具体的に説明する。図２は、この燃料噴射量制御装置の機能ブロック図である。燃料噴射量制御装置は、電子制御スロットルモデルＭ１、スロットルモデルＭ２、吸気管モデルＭ３、吸気弁モデルＭ４、シリンダモデルＭ５及び排気弁モデルＭ６からなるシミュレーションモデルを用いて筒内吸入新気量Ｍａを推定するようになっている。
【００５８】
これらのモデルＭ１〜Ｍ６における計算は、ＣＰＵ８１が所定時間（プログラム実行間隔時間Δｔ）の経過毎に各モデルに対応するプログラムを実行することにより達成される。また、以下に述べる微分方程式は、実際には離散化されることによりその解が求められる。更に、燃料噴射量制御装置は、スロットルバルブ電子制御ロジックＡ１を備え、スロットルバルブアクチュエータ４３ａを介してスロットルバルブ４３の開度を制御するようになっている。以下、各モデルについて順に説明する。
【００５９】
（電子制御スロットルモデルＭ１及びスロットルバルブ電子制御ロジックＡ１）
電子制御スロットルモデルＭ１は、現時点までのアクセルペダル操作量Ａｃｃｐに基づいて現時点から所定時間Ｔ０先の時刻ｔにおけるスロットルバルブ開度θｔを推定するモデルである。
【００６０】
具体的に述べると、スロットルバルブ電子制御ロジックＡ１は、アクセルペダル操作量Ａｃｃｐと目標スロットルバルブ開度θｒとの関係を規定する図３に示したテーブル及びアクセル開度センサ７１により検出された実際の（現時点の）アクセルペダル操作量Ａｃｃｐに基づいて暫定的な目標スロットルバルブ開度θｒ１を所定時間（例えば、２ｍｓｅｃ）の経過毎に決定する。アクセルペダル操作量Ａｃｃｐと目標スロットルバルブ開度θｒとの関係を規定するテーブル（マップ）は、アクセルペダル操作量Ａｃｃｐが増大するとともに目標スロットルバルブ開度θｒが増大するようにこれらの関係を規定している。
【００６１】
また、スロットルバルブ電子制御ロジックＡ１は、暫定的な目標スロットルバルブ開度θｒ１を所定時間Ｔ（例えば、６４ｍｓｅｃ）だけ遅延させた値、即ち、現時点より所定時間Ｔだけ前の時点にて決定された暫定的な目標スロットルバルブ開度θｒ１を、現時点での最終的な目標スロットルバルブ開度θｒとして決定する。そして、スロットルバルブ電子制御ロジックＡ１は、実際のスロットルバルブ開度ＴＡが現時点の目標スロットルバルブ開度θｒとなるようにスロットルバルブアクチュエータ４３ａに対して駆動信号を送出する。
【００６２】
このように、目標スロットルバルブ開度θｒは、現時点から所定時間Ｔだけ前の時点におけるアクセルペダル操作量Ａｃｃｐに応じて決定された暫定的な目標スロットルバルブ開度θｒ１と等しいから、現時点から所定時間Ｔ０だけ先の時刻ｔにおける目標スロットルバルブ開度θｒは現時点から時間（Ｔ−Ｔ０）前における暫定的な目標スロットルバルブ開度θｒ１と等しい。また、現時点から時間（Ｔ−Ｔ０）前における暫定的な目標スロットルバルブ開度θｒ１は、スロットルバルブアクチュエータ４３ａの作動遅れ時間を無視すれば、時刻ｔにおけるスロットルバルブ開度θｔと等しい。
【００６３】
このような考えに基づき、電子制御スロットルモデルＭ１は、現時点から所定時間Ｔ０だけ先の時刻ｔにおけるスロットルバルブ開度θｔを推定する。即ち、現時点から時間（Ｔ−Ｔ０）前における暫定的な目標スロットルバルブ開度θｒ１を現時点から所定時間Ｔ０だけ先の時刻ｔにおけるスロットルバルブ開度θｔとして推定する。なお、スロットルバルブアクチュエータ４３ａの作動遅れ時間を考慮に加えて、スロットルバルブ開度θｔを推定してもよい。
【００６４】
（スロットルモデルＭ２）
スロットルモデルＭ２は、スロットルバルブ４３を通過する空気流量（（以下、「スロットル通過空気流量」と称呼する。）ｍｔを、エネルギー保存則、運動量保存則、質量保存則及び状態方程式等の物理法則に基づいて得られた下記（９）式及び下記（１０）式に基づいて推定するモデルである。
【００６５】
下記（９）式及び下記（１０）式において、Ｃｔ（θｔ）はスロットルバルブ開度θｔに応じて変化する流量係数、Ａｔ（θｔ）はスロットルバルブ開度θｔに応じて変化するスロットル開口面積（吸気管４１の開口面積）、Ｐａはスロットルバルブ上流圧力（即ち、大気圧）、Ｐｍは吸気管内空気圧力（吸気管圧力）、Ｔａは吸気温度（大気温度）、Ｔｍはスロットルバルブ４３の下流の吸気管内空気温度、Ｒは気体定数、及びκは比熱比（ここではκを一定値として扱う。）である。スロットルモデルＭ２は、スロットルバルブ上流圧力Ｐａが吸気管内空気圧力Ｐｍより大きい順流の場合（Ｐａ≧Ｐｍ）に（９）式を使用し、スロットルバルブ上流圧力Ｐａが吸気管内空気圧力Ｐｍより小さい逆流の場合（Ｐａ＜Ｐｍ）に（１０）式を使用する。
【００６６】
【数９】

【００６７】
【数１０】

【００６８】
上記（９）式及び（１０）式において、θｔは電子制御スロットルモデルＭ１により推定された現時点から所定時間Ｔ０だけ先の時刻ｔにおける推定スロットルバルブ開度である。スロットルモデルＭ２は、スロットルバルブ開度θｔと流量係数Ｃｔ（θｔ）との関係を規定した図４に示すテーブル及び前記推定したスロットルバルブ開度θｔを用いて流量係数Ｃｔ（θｔ）を求めるとともに、スロットルバルブ開度θｔと開口面積Ａｔ（θｔ）との関係を規定した図５に示すテーブル及び前記推定したスロットルバルブ開度θｔとを用いて開口面積Ａｔ（θｔ）を求める。
【００６９】
なお、スロットルモデルＭ２は、スロットルバルブ開度θｔと積値Ｃｔ（θｔ）・Ａｔ（θｔ）との関係を規定した図６に示すテーブル及び前記推定したスロットルバルブ開度θｔを用いて積値Ｃｔ（θｔ）・Ａｔ（θｔ）を一時に求めるように構成してもよい。また、スロットルモデルＭ２は、スロットルバルブ開度θｔ及び吸気管内空気圧力Ｐｍと流量係数Ｃｔ（θｔ，Ｐｍ）との関係を規定したテーブルＭａｐＣｔ（θｔ，Ｐｍ）と、前記推定したスロットルバルブ開度θｔ及び後述する吸気管モデルＭ３から取得される吸気管内空気圧力Ｐｍと、を用いて、流量係数Ｃｔ（θｔ）に代わる流量係数Ｃｔ（θｔ，Ｐｍ）を求めるように構成されていてもよい。
【００７０】
スロットルモデルＭ２は、スロットルバルブ上流圧力Ｐａ及び吸気温度Ｔａを大気圧センサ６３及び吸気温センサ６２からそれぞれ取得するとともに、吸気管内空気圧力Ｐｍと吸気管内空気温度Ｔｍとを後述する吸気管モデルＭ３から取得し、これらの値を用いて上記（９）式又は（１０）式の計算を行い、時刻ｔにおけるスロットル通過空気流量ｍｔを推定する。なお、吸気管モデルＭ３は、後述するように、スロットルモデルＭ２の出力であるスロットル通過空気流量ｍｔを用いて計算を行う。従って、スロットルモデルＭ２が吸気管モデルＭ３から取得する吸気管内空気圧力Ｐｍと吸気管内空気温度Ｔｍは、前回の（現時点からプログラム実行間隔時間Δｔだけ前の）計算タイミングにて同吸気管モデルＭ３が計算していた値である。かかる計算手法は、他のモデルにおいても同様に使用される。
【００７１】
ここで、上記スロットルモデルＭ２を記述した（９）式及び（１０）式の導出過程について説明する。いま、スロットルバルブ４３の上流の開口断面積をＡｕ、空気密度をρｕ、空気の流速をｖｕとし、スロットルバルブ４３による吸気管４１の開口断面積をＡｄ、そこでの空気密度をρｄ、スロットルバルブ４３を通過する空気の流速をｖｄとすると、スロットル通過空気流量ｍｔは、下記（１１）式で表される。（１１）式は質量保存則を記述した式と言える。
【００７２】
【数１１】
ｍｔ＝Ａｄ・ρｄ・ｖｄ＝Ａｕ・ρｕ・ｖｕ …（１１）
【００７３】
一方、運動エネルギーは、空気の質量をｍとすると、スロットルバルブ４３の上流でｍ・ｖｕ^２／２であり、スロットルバルブ４３を通過する場所でｍ・ｖｄ^２／２である。他方、熱エネルギーは、スロットルバルブ４３の上流でｍ・Ｃｐ・Ｔｕであり、スロットルバルブ４３を通過する場所でｍ・Ｃｐ・Ｔｄである。従って、エネルギー保存則により、下記（１２）式が得られる。なお、Ｔｕはスロットルバルブ上流の空気温度、Ｔｄはスロットルバルブ下流の空気温度、Ｃｐは定圧比熱である。
【００７４】
【数１２】
ｍ・ｖｕ^２／２＋ｍ・Ｃｐ・Ｔｕ＝ｍ・ｖｄ^２／２＋ｍ・Ｃｐ・Ｔｄ …（１２）
【００７５】
ところで、状態方程式は下記（１３）式、比熱比κは下記（１４）式、マイヤーの関係は下記（１５）式で示されるから、（１３）式〜（１５）式よりＣｐ・Ｔは下記（１６）式のように表される。なお、Ｐは気体の圧力、ρは気体の密度、Ｔは気体の温度、Ｒは気体定数、Ｃｖは定容比熱である。
【００７６】
【数１３】
Ｐ＝ρ・Ｒ・Ｔ …（１３）
【００７７】
【数１４】
κ＝Ｃｐ／Ｃｖ …（１４）
【００７８】
【数１５】
Ｃｐ＝Ｃｖ＋Ｒ …（１５）
【００７９】
【数１６】
Ｃｐ・Ｔ＝｛κ／（κ−１）｝・（Ｐ／ρ） …（１６）
【００８０】
上記（１６）式の関係を用いて上記エネルギー保存則に基づく（１２）式を書換えると、下記（１７）式が得られる。
【００８１】
【数１７】
ｖｕ^２／２＋｛κ／（κ−１）｝・（Ｐｕ／ρｕ）＝ｖｄ^２／２＋｛κ／（κ−１）｝・（Ｐｄ／ρｄ） …（１７）
【００８２】
そして、スロットルバルブ４３の無限上流を考えると、Ａｕ＝∞、ｖｕ＝０であるから、エネルギー保存則に基づく上記（１７）式は下記（１８）式に書き換えられる。
【００８３】
【数１８】
｛κ／（κ−１）｝・（Ｐｕ／ρｕ）＝ｖｄ^２／２＋｛κ／（κ−１）｝・（Ｐｄ／ρｄ） …（１８）
【００８４】
次に、運動量について記述する。断面積Ａｕの部分に加わる圧力をＰｕ、断面積Ａｄの部分に加わる圧力をＰｄ、断面積Ａｕの部分と断面積Ａｄの部分との間をつなぐ固定された空間の平均圧力をＰｍｅａｎとすると、下記（１９）式が得られる。
【００８５】
【数１９】
ρｄ・ｖｄ^２・Ａｄ−ρｕ・ｖｕ^２・Ａｕ＝Ｐｕ・Ａｕ−Ｐｄ・Ａｄ＋Ｐｍｅａｎ・（Ａｄ−Ａｕ） …（１９）
【００８６】
上記（１９）式で、Ａｕ＝∞、ｖｕ＝０を考慮すると、下記（２０）式が得られるので、同（２０）式と上記（１９）式とから下記（２１）式の運動量に関する関係（運動量保存則に基づく関係）が得られる。
【００８７】
【数２０】
Ｐｍｅａｎ＝Ｐｕ …（２０）
【００８８】
【数２１】
ρｄ・ｖｄ^２＝Ｐｕ−Ｐｄ …（２１）
【００８９】
従って、上記（１１）式、上記（１８）式、及び上記（２１）式から、下記（２２）式が得られる。
【００９０】
【数２２】

【００９１】
上記（２２）式において、Ｐｕはスロットルバルブ上流圧力Ｐａであり、Ｐｄは吸気管内空気圧力Ｐｍである。また、ρｕ＝Ｐｕ／（Ｒ・Ｔｕ）である。これらの関係を使用するとともに、上記（２２）式に流量係数Ｃｔ（θｔ）を適合のための係数として導入し、開口断面積Ａｄを開口面積Ａｔ（θｔ）とおきなおして整理すると、上記（９）式が得られる。上記（１０）式の導出過程は、上記（９）式の導出過程と同様であるので省略する。
【００９２】
なお、（９）式及び（１０）式は絞り部（オリフィス）を通過する流体の流量ｍｓを表す一般的な下記（２３）式乃至（２５）式から求めることができる。（２３）式及び（２４）式において、Ｃは流量係数、Ａは絞り部の開口面積、Ｐｕは絞り部の上流の圧力、Ｔｕは絞り部の上流のガス温度、Ｐｄは絞り部の下流の圧力、Ｔｄは絞り部の下流のガス温度である。（２５）式は、関数Φを定義する式であり、（２５）式の上段は流速が音速以下の場合、下段は流速が音速になる場合に使用される。
【００９３】
【数２３】

【００９４】
【数２４】

【００９５】
【数２５】

【００９６】
（吸気管モデルＭ３）
吸気管モデルＭ３は、質量保存則とエネルギー保存則とにそれぞれ基づいた下記（２６）式及び下記（２７）式、スロットル通過空気流量ｍｔ、スロットル通過空気温度（即ち、吸気温度）Ｔａ及び吸気管から流出する空気流量である吸気弁通過ガス流量ｍｉから、吸気管内空気圧力Ｐｍ及び吸気管内空気温度Ｔｍを求めるモデルである。なお、下記（２６）式及び下記（２７）式において、Ｖｍはスロットルバルブ４３から吸気弁３２までの吸気管４１（以下、単に「吸気管部」と称呼する。）の容積である。
【００９７】
【数２６】
ｄ（Ｐｍ／Ｔｍ）／ｄｔ＝（Ｒ／Ｖｍ）・（ｍｔ−ｍｉ） …（２６）
【００９８】
【数２７】
ｄＰｍ／ｄｔ＝κ・（Ｒ／Ｖｍ）・（ｍｔ・Ｔａ−ｍｉ・Ｔｍ） …（２７）
【００９９】
吸気管モデルＭ３は、上記（２６）式及び上記（２７）式におけるスロットル通過空気流量ｍｔをスロットルモデルＭ２から取得し、吸気弁通過ガス流量ｍｉを後述する吸気弁モデルＭ４から取得する。そして、（２６）式及び（２７）式に基づく計算を行って、現時点から所定時間Ｔ０だけ先の時刻ｔにおける吸気管内空気圧力Ｐｍ及び吸気管内空気温度Ｔｍを求める。
【０１００】
ここで、上記吸気管モデルＭ３を記述した（２６）式及び（２７）式の導出過程について説明する。いま、吸気管部の総空気量をＭとすると、総空気量Ｍの時間的変化は、吸気管部に流入する空気量に相当するスロットル通過空気流量ｍｔと同吸気管部から流出する空気量に相当する吸気弁通過ガス流量ｍｉの差であるから、質量保存則に基づく下記（２８）式が得られる。
【０１０１】
【数２８】
ｄＭ／ｄｔ＝ｍｔ−ｍｉ …（２８）
【０１０２】
また、状態方程式は下記（２９）式となるから、上記（２８）式と下記（２９）式とから総空気量Ｍを消去することにより、質量保存則に基づく上記（２６）式が得られる。
【０１０３】
【数２９】
Ｐｍ・Ｖｍ＝Ｍ・Ｒ・Ｔｍ …（２９）
【０１０４】
次に、吸気管部に関するエネルギー保存則について検討する。この場合、吸気管部の容積Ｖｍは変化せず、また、エネルギーの殆どが温度上昇に寄与する（運動エネルギーは無視し得る）と考えられる。従って、吸気管部の空気のエネルギーＭ・Ｃｖ・Ｔｍの時間的変化量は、同吸気管部に流入する空気のエネルギーＣｐ・ｍｔ・Ｔａと同吸気管部から流出する空気のエネルギーＣｐ・ｍｉ・Ｔｍとの差に等しいので、下記（３０）式が得られる。
【０１０５】
【数３０】
ｄ（Ｍ・Ｃｖ・Ｔｍ）／ｄｔ＝Ｃｐ・ｍｔ・Ｔａ−Ｃｐ・ｍｉ・Ｔｍ …（３０）
【０１０６】
この（３０）式を、上記（１４）式（κ＝Ｃｐ／Ｃｖ）と、上記（２９）式（Ｐｍ・Ｖｍ＝Ｍ・Ｒ・Ｔｍ）とを用いて変形することにより、上記（２７）式が得られる。
【０１０７】
（吸気弁モデルＭ４）
吸気弁モデルＭ４は、図７に示した吸気弁３２の周囲を通過する空気流量である吸気弁通過ガス流量ｍｉを、エネルギー保存則、運動量保存則、質量保存則及び状態方程式等の物理法則に基づいて得られた下記（３１）式及び下記（３２）式にしたがって推定するモデルである。
【０１０８】
（３１），（３２）式の導出過程は、上記スロットルモデルＭ２の場合と同様である。（３１）式及び（３２）式において、Ｃｄｉ（Ｌｉ）は吸気弁３２のリフト量Ｌｉに応じて変化する流量係数、Ａｉ（Ｌｉ）は同リフト量Ｌｉに応じて変化する吸気弁３２の周囲に形成される開口の面積、Ｔｉは吸気弁３２の上流部のガス温度（吸気弁を通過するガスの温度であり、（以下、「吸気弁通過ガス温度Ｔｉ」と称呼する。））、及びＰｃは筒内ガス圧力（シリンダ２１内の圧力Ｐｃ）である。
【０１０９】
吸気弁モデルＭ４は、吸気弁上流ガス圧力Ｐｉが筒内ガス圧力Ｐｃより大きい順流の場合に（３１）式を使用し、吸気弁上流ガス圧力Ｐｉが筒内ガス圧力Ｐｃより小さい逆流の場合に（３２）式を使用する。このように、吸気弁通過ガス流量ｍｉは、シリンダ２１内にガスが吸入されている場合に正、シリンダ２１内からガスが吸気通路に吹き返されている場合に負の値をとるように規定されている。
【０１１０】
【数３１】

【０１１１】
【数３２】

【０１１２】
吸気弁モデルＭ４は、現時点から所定時間Ｔ０だけ先の時刻ｔにおける吸気弁リフト量Ｌｉ（ｔ）を推定する。この推定は、吸気弁リフト量センサ６７が検出している現時点のリフト量Ｌｉとエンジン回転速度ＮＥとに基づいてなされてもよく、現時点のクランク角とエンジン回転速度ＮＥを含む運転状態に応じて予め設定されるリフト量マップとに基いてなされてもよい。そして、吸気弁モデルＭ４は、吸気弁リフト量Ｌｉと積値Ｃｄｉ（Ｌｉ）・Ａｉ（Ｌｉ）との関係を規定した図８に示したテーブルと、前記推定した吸気弁リフト量Ｌｉ（ｔ）とに基づいて、上記（３１）式及び上記（３２）式にて使用する積値Ｃｄｉ（Ｌｉ）・Ａｉ（Ｌｉ）を求める。
【０１１３】
吸気弁モデルＭ４は、後にフローチャートを参照しながら詳述するように、吸気弁上流ガス圧力Ｐｉを吸気管モデルＭ３から取得される吸気管内空気圧力Ｐｍと等しい値に設定する。吸気弁モデルＭ４は、上記（３１）式における吸気弁通過ガス温度Ｔｉを、シリンダ２１が新気を吸入している場合には吸気管モデルＭ３から取得される吸気管内空気温度Ｔｍと等しい値に設定し、シリンダ２１が吸気通路に吹き返された既燃ガスを吸入している場合には排気ガス温度Ｔｅに設定する。吸気弁モデルＭ４は、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃを後述するシリンダモデルＭ５から取得する。吸気弁モデルＭ４は、これらの変数を用いて上記（３１）式又は上記（３２）式に基づく計算を行うことで、時刻ｔにおける吸気弁通過ガス流量ｍｉを推定する。
【０１１４】
なお、上記（３１）式及び上記（３２）式に代え、これらと等価な（３３）式乃至下記（３５）式により吸気弁通過ガス流量ｍｉを求めても良い。これらの式において、Ｃｄｉ＝Ｃｄｉ（Ｌｉ）、Ａｉ＝Ａｉ（Ｌｉ）である。（３５）式は、（２５）式で説明した関数Φを定義する式であり、（３５）式の上段は流速が音速以下の場合、下段は流速が音速になる場合に使用される。なお、（３１）式〜（３５）式は、吸気弁３２の周囲に形成されるガス通路を絞り部とみなしたときに成立する同絞り部を通過するガス流量を算出する一般的な式である。
【０１１５】
【数３３】

【０１１６】
【数３４】

【０１１７】
【数３５】

【０１１８】
（シリンダモデルＭ５）
シリンダモデルＭ５は、シリンダ２１（燃焼室２５）内の空気についての質量保存則及びエネルギー保存則に基づいた下記（３６）式及び下記（３７）式にしたがって、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃを求めるモデルである。図７に示したように、Ｖｃはシリンダ２１（燃焼室２５）の容積、Ｔｉは吸気弁通過ガス温度（吸気ガス温度）、ｍｉは吸気弁通過ガス流量、ｍｅは排気弁通過ガス流量、Ｔｅは排気ガス温度（排気弁通過ガス温度）である。
【０１１９】
κｃ，κｉ及びκｅは、それぞれシリンダ２１内のガスの比熱比（筒内ガス比熱比）、吸気弁通過ガスの比熱比及び排気ガス（排気弁通過ガス）の比熱比であり、ここでは一定値κとして扱う。Ｒｃ，Ｒｉ及びＲｅは、それぞれシリンダ２１内のガスの気体定数、吸気弁通過ガスの気体定数及び排気ガスの気体定数であり、ここでは一定値Ｒとして扱う。Ｑｗ’はシリンダ２１内のガスに同シリンダ２１（シリンダ壁面）から伝達される単位時間あたりの熱量（伝達熱流量、熱流）である。
【０１２０】
【数３６】

【０１２１】
【数３７】

【０１２２】
（３６）式及び（３７）式における時刻ｔの吸気弁通過ガス流量ｍｉは吸気弁モデルＭ４により与えられ、排気弁通過ガス流量ｍｅは後述する排気弁モデルＭ６により与えられる。時刻ｔのシリンダ容積Ｖｃはクランクポジションセンサ６８が検出している実際のクランク角とエンジン回転速度ＮＥとから求めることができる。時刻ｔの吸気弁通過ガス温度Ｔｉは前述した吸気弁モデルＭ４と同様に取得される。排気ガス温度Ｔｅは、総べての気筒の排気弁通過ガス流量ｍｅを積算（積分、時間積分）して得られる単位時間あたりの排気ガス量Ｇｅ（又は、総べての気筒の吸気弁通過ガス流量ｍｉを積分（時間積分）して得られる単位時間あたりの吸入ガス量Ｇａ）、総べて気筒の燃料噴射量ｆｃを積分（時間積分）して得られる単位時間あたりの燃料量Ｇｆ及びエンジン回転速度ＮＥの関数（マップ値）として求められる。なお、排気ポート３４に排気ガス温度センサを設置し、この排気ガス温度センサの出力から排気ガス温度Ｔｅを求めてもよい。
【０１２３】
伝達熱流量Ｑｗ’は比較的小さいので、この例においては無視する。以上のことから、（３６）式及び（３７）式に基づいて筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃが求められる。なお、実際には、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃは、（３６）式及び（３７）式が離散化された式により求められる。
【０１２４】
ここで、上記（３６）式の導出過程について説明する。先ず、質量保存則に基づけば、下記（３８）式を得ることができる。Ｍｃはシリンダ２１内のガス量（筒内ガス量）である。（３８）式は「新気及び既燃ガスからなるシリンダ内のガスの量である筒内ガス量Ｍｃの時間的変化量（ｄＭｃ／ｄｔ）」は、「吸気弁通過ガス流量ｍｉ」と「排気弁通過ガス流量ｍｅ」との和に等しいことを表している。
【０１２５】
【数３８】

【０１２６】
一方、気体の状態方程式は下記の（３９）式のとおりである。
【０１２７】
【数３９】

【０１２８】
（３８）式に（３９）式を適用して筒内ガス量Ｍｃを消去することにより下記（４０）式が得られる。（４０）式を変形すると下記（４１）式が得られ、（４１）式を整理することにより上記（３６）式が得られる。
【０１２９】
【数４０】

【０１３０】
【数４１】

【０１３１】
次に、上記（３７）式の導出過程について説明する。先ず、エネルギー保存則に基づけば、下記（４２）式を得ることができる。
【０１３２】
【数４２】

【０１３３】
（４２）式において、Ｅはシリンダ内のガスのエネルギー、Ｅｉ’は吸気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギー、Ｅｅ’は排気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギー、Ｗ’はシリンダ内のガスが単位時間あたりにピストンに対して行う仕事、及びＱｗ’はシリンダ壁面からシリンダ内のガスに単位時間あたりに与えられる熱量（伝達熱流量）である。
【０１３４】
（４２）式は、「シリンダ内のガスのエネルギーＥの時間的変化量（ｄＥ／ｄｔ）」が、「吸気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｉ’」と、「排気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｅ’」と、「伝達熱流量Ｑｗ’」との和から、「シリンダ内のガスが単位時間あたりにピストンに対して行う仕事Ｗ’」を減じた値に等しいというエネルギー保存則を表している。
【０１３５】
ここで（４２）式の各変数について検討すると、シリンダ内のガスのエネルギーＥの時間的変化量（ｄＥ／ｄｔ）について、ｕを内部エネルギー、Ｃｖｃをシリンダ内のガスの定容比熱とするとき、下記（４３）式が成立する。
【０１３６】
【数４３】

【０１３７】
また、吸気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｉ’について下記（４４）式が成立する。（４４）式において、ｈｉは吸気弁通過ガスのエンタルピーであり、Ｃｐｉは吸気弁通過ガスの定圧比熱である。なお、エンタルピーｈは、ｈ＝Ｃｐ・Ｔ（Ｃｐは対象とするガスの定圧比熱、Ｔはそのガスの温度）と定義される。
【０１３８】
【数４４】

【０１３９】
同様に、排気弁通過ガスにより単位時間あたりにシリンダ内のガスに与えられるエネルギーＥｅ’について下記（４５）式が成立する。（４５）式において、ｈｅは排気弁通過ガスのエンタルピーであり、Ｃｐｅは排気弁通過ガスの定圧比熱である。
【０１４０】
【数４５】

【０１４１】
更に、シリンダ内のガスが単位時間あたりにピストンに対して行う仕事Ｗ’については、下記（４６）式が成立する。
【０１４２】
【数４６】

【０１４３】
これら（４３）式乃至（４６）式を（４２）式に適用すると、下記（４７）式が得られる。
【０１４４】
【数４７】

【０１４５】
この（４７）式に下記（４８）式の気体の状態方程式を適用して筒内ガス量Ｍｃを消去すると、下記（４９）式が得られる。
【０１４６】
【数４８】

【０１４７】
【数４９】

【０１４８】
一方、上記（１４）式（κ＝Ｃｐ／Ｃｖ）及び上記（１５）式（Ｃｐ＝Ｃｖ＋Ｒ）から、ガスｊ（ｊ＝ｉ，ｃ，ｅ）について下記（５０）式及び下記（５１）式が成立する。この（５０）式及び（５１）式を（４９）式に適用することにより、上記（３７）式が得られる。
【０１４９】
【数５０】

【０１５０】
【数５１】

【０１５１】
なお、本例のシリンダモデルＭ６は、比熱比κｉ，κｅ，κｃを一定の比熱比κ、気体定数Ｒｉ，Ｒｅ，Ｒｃを一定の気体定数Ｒとして扱って、計算を行うようになっている。
【０１５２】
（排気弁モデルＭ６）
排気弁モデルＭ６は、図７に示した排気弁３５の周囲を通過する空気流量（即ち、排気弁通過ガス流量）ｍｅを、エネルギー保存則、運動量保存則、質量保存則、及び状態方程式等の物理法則に基づいて得られた下記（５２）式及び下記（５３）式にしたがって推定するモデルである。
【０１５３】
（５２），（５３）式の導出過程は、上記スロットルモデルＭ２の場合と同様である。（５２）式及び（５３）式において、Ｃｄｅ（Ｌｅ）は排気弁３５のリフト量Ｌｅに応じて変化する流量係数、Ａｅ（Ｌｅ）は同リフト量Ｌｅに応じて変化する排気弁３５の周囲に形成される開口の面積、Ｔｅは排気ガス温度、Ｐｃは筒内ガス圧力及びＰｅは排気ガス圧力（排気弁３５の下流側のガスの圧力）である。
【０１５４】
排気弁モデルＭ６は、排気ガス圧力Ｐｅが筒内ガス圧力Ｐｃより大きい場合に（５２）式を使用し、排気ガス圧力Ｐｅが筒内ガス圧力Ｐｃより小さい場合に（５３）式を使用する。このように、排気弁通過ガス流量ｍｅは、排気系からシリンダ２１内にガスが吸入されている場合に正、シリンダ２１内からガスが排気系に排出されている場合に負の値をとるように規定されている。
【０１５５】
【数５２】

【０１５６】
【数５３】

【０１５７】
排気弁モデルＭ６は、現時点から所定時間Ｔ０だけ先の時刻ｔにおける排気弁リフト量Ｌｅ（ｔ）を推定する。この推定は、現時点のクランク角とエンジン回転速度ＮＥとに基づいて行われ得る。そして、排気弁モデルＭ６は、排気弁リフト量Ｌｅと積値Ｃｄｅ（Ｌｅ）・Ａｅ（Ｌｅ）との関係を規定した図９に示したテーブルと、前記推定した排気弁リフト量Ｌｅ（ｔ）とに基づいて、上記（５２）式及び上記（５３）式にて使用する積値Ｃｄｅ（Ｌｅ）・Ａｅ（Ｌｅ）を求める。
【０１５８】
排気弁モデルＭ６は、排気ガス圧力Ｐｅを、総べての気筒の排気弁通過ガス流量ｍｅを積算（時間積分）して得られる単位時間あたりの排気ガス量Ｇｅ（又は、総べての気筒の吸気弁通過ガス流量ｍｉを積算（時間積分）して得られる単位時間あたりの吸入ガス量Ｇａ）の関数として求める。なお、排気ポート３４に排気圧力センサを設置し、この排気圧力センサの出力から排気ガス圧力Ｐｅを求めてもよい。
【０１５９】
また、排気弁モデルＭ６は、排気ガス温度Ｔｅを、排気弁通過ガス流量ｍｅを積算（時間積分）して得られる単位時間あたりの排気ガス量Ｇｅ、総べての気筒の燃料噴射量ｆｃを積算（時間積分）して得られる単位時間あたりの燃料量Ｇｆ及びエンジン回転速度ＮＥの関数（マップ値）として求める。排気弁モデルＭ６は筒内ガス圧力Ｐｃ及び筒内ガス温度ＴｃをシリンダモデルＭ５から取得する。排気弁モデルＭ６は、これらの変数を用いて上記（５２）式又は上記（５３）式を計算することで、時刻ｔにおける排気弁通過ガス流量ｍｅを推定する。
【０１６０】
なお、上記（５２）式及び上記（５３）式に代え、これらと等価な下記（５４）式、下記（５５）式及びΦの定義式である上記（２５）式により排気弁通過ガス流量ｍｅを求めても良い。これらの式において、Ｃｄｅ＝Ｃｄｅ（Ｌｅ）、Ａｅ＝Ａｅ（Ｌｅ）である。
【０１６１】
【数５４】

【０１６２】
【数５５】

【０１６３】
本燃料噴射量制御装置は、このような一連の計算により求められる吸気弁通過ガス流量ｍｉを吸気弁３２が開弁する時刻ｔｏｐｅｎから同吸気弁３２が閉弁する時刻ｔｃｌｓｅまで時間積分することにより、一吸気行程にてシリンダ２１内に吸入される筒内吸入新気量Ｍａを推定し、この値Ｍａに基づいて燃料噴射量ｆｃを決定する。なお、吸気弁通過ガス流量ｍｉの積算値Ｍｉは、吸気弁３２が開弁する時刻ｔｏｐｅｎにて初期値０を与えれば、吸気通路に吹き返されたガスを総べて再吸入した時点で「０」になるので、結果として、筒内吸入新気量Ｍａを示す値となる。以上が、筒内吸入新気量Ｍａ及び燃料噴射量ｆｃを決定する原理である。
【０１６４】
次に、本発明の第１実施形態に係る燃料噴射量制御装置の作動について説明する。第１実施形態の燃料噴射量制御装置は、吸気弁３２の開弁にともない既燃ガスが吸気通路に吹き返され始めた時点から、その吹き返された既燃ガスがシリンダ２１内に完全に吸入されるまでの期間（実際には、吸気通路に吹き返されたガスがシリンダ２１内に再吸気される期間）において、吸気弁モデルＭ４にて使用する吸気弁通過ガス温度Ｔｉを排気ガス温度Ｔｅと等しい値に設定し、新気がシリンダ２１内に吸入され始めてからは吸気弁モデルＭ４にて使用する吸気弁通過ガス温度Ｔｉを吸気管モデルＭ３が求めた吸気管内空気温度Ｔｍと等しい温度に設定する。
【０１６５】
以下、ＣＰＵ８１が実行するルーチン（プログラム）を示した図１０のフローチャート及び図１１のタイムチャートを参照して、燃料噴射量制御装置の作動について説明する。ＣＰＵ８１は、図１０に示したルーチンを所定時間（プログラム実行間隔時間Δｔ）の経過毎に、且つ、特定の気筒（ある気筒）に対して実行するようになっている。このルーチン上での時刻は現在より所定時間Ｔ０だけ先の時刻ｔである。即ち、例えば、このルーチンの実行により吸気弁３２が閉弁状態から開弁状態に移行したと判定された時点から所定時間Ｔ０が経過した後に実際に吸気弁３２が閉弁状態から開弁状態に移行する。なお、ＣＰＵ８１は、同様のルーチンを他の気筒に対しても独立して実行している。
【０１６６】
所定のタイミングになると、ＣＰＵはステップ１０００から本ルーチンの処理を開始し、ステップ１００２に進んで特定の気筒の吸気弁３２が閉弁状態から開弁状態へと移行した直後であるか否かを判定する。いま、特定の気筒の吸気弁３２が閉弁状態にあるものとして説明を続けると、ＣＰＵ８１はステップ１００２にて「Ｎｏ」と判定してステップ１００４に進み、その吸気弁３２が開弁状態にあるか否かを判定する。
【０１６７】
この場合、特定の気筒の吸気弁３２は閉弁状態にあるから、ＣＰＵ８１はステップ１００４にて「Ｎｏ」と判定してステップ１００６に進み、吸気弁３２が開弁状態から閉弁状態へと移行した直後であるか否かを判定する。この場合においても、吸気弁３２は閉弁状態にあるから、ＣＰＵ８１はステップ１００６にて「Ｎｏ」と判定してステップ１０９５に進み、本ルーチンを一旦終了する。このような処理は、吸気弁３２が閉弁状態から開弁状態へと移行する時点まで繰り返し行われる。
【０１６８】
その後、吸気弁３２は、クランク角が吸気上死点ＴＤＣより前の開弁クランク角θ０となると、閉弁状態から開弁状態へと移行する（図１１の時刻ｔ０を参照）。この場合、ＣＰＵ８１はステップ１００２に進んだとき、同ステップ１００２にて「Ｙｅｓ」と判定してステップ１００８に進み、吸気管モデルＭ３が推定している吸気管内空気圧力Ｐｍ及び吸気管内空気温度Ｔｍを取得する。
【０１６９】
次いで、ＣＰＵ８１はステップ１０１０にて、吸気弁上流ガス圧力Ｐｉを前記取得した吸気管内空気圧力Ｐｍと等しい値に設定し、ステップ１０１２にて同ステップ１０１２中に記載した式（Ｐｅ＝ａ・Ｇｅ＋ｂ・Ｇｅ^２， ａ，ｂは定数）にしたがって排気ガス圧力Ｐｅを推定する。ここで、Ｇｅは総べての気筒の排気弁通過ガス流量ｍｅを積算（時間積分）して得られる単位時間あたりの排気ガス量であり、図示しないルーチンにより別途計算されている。
【０１７０】
次に、ＣＰＵ８１はステップ１０１４に進み、総べての気筒の吸気弁通過ガス流量ｍｉを積算（時間積分）して得られる単位時間あたりの吸入ガス量Ｇａ、総べての気筒の燃料噴射量ｆｃを積算（時間積分）して得られる単位時間あたりの燃料量Ｇｆ及びエンジン回転速度ＮＥと、これらの値と排気ガス温度Ｔｅとの関係を規定したテーブルＭａｐＴｅとから排気ガス温度Ｔｅ（＝ＭａｐＴｅ（Ｇａ，Ｇｆ，ＮＥ））を求める。
【０１７１】
次いで、ＣＰＵ８１はステップ１０１６にて各値の初期設定を行う。具体的に述べると、ＣＰＵ８１は、筒内ガス圧力Ｐｃに排気ガス圧力Ｐｅを設定するとともに、筒内ガス温度Ｔｃに排気ガス温度Ｔｅを設定する。これは、吸気弁３２が開弁した時刻ｔ０においては排気弁３５が開弁しているから（図１１の破線ＬＥＸにて示した排気弁３５のリフト量を参照。）、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃが排気ガス圧力Ｐｅ及び排気ガス温度Ｔｅにそれぞれ近似した値となるであろうとの知見に基づく。
【０１７２】
更に、ＣＰＵ８１は、吸気弁通過ガス温度Ｔｉを排気ガス温度Ｔｅと等しい値に設定する。これは、吸気弁３２と排気弁３５が共に開弁しているので、排気ポート３４内のガスがシリンダ２１を介して圧力の低い吸気ポート３１に吹き返されるため、吸気弁通過ガス温度Ｔｉが排気ガス温度Ｔｅに近似した値となるであろうとの知見に基づく。また、ＣＰＵ８１は、吸気弁３２が開弁したときからの吸気弁通過ガス流量ｍｉの積算値Ｍｉを初期値「０」に設定するとともに、吸気弁３２が開弁したときからの排気弁通過ガス流量ｍｅの積算値Ｍｅを初期値「０」に設定する
【０１７３】
次いで、ＣＰＵ８１はステップ１０１８に進み、吸気弁モデルＭ４によって吸気弁通過ガス流量ｍｉを算出する。この場合、吸気弁上流ガス圧力Ｐｉ（＝吸気管内空気圧力Ｐｍ）は筒内ガス圧力Ｐｃ（＝排気ガス圧力Ｐｅ）より小さいから、シリンダ２１内から吸気ポート３１へと既燃ガスが逆流する。即ち、（３２）式又は（３４）式が使用されて吸気弁通過ガス流量ｍｉが求められる。従って、吸気弁通過ガス流量ｍｉは負の値となる。
【０１７４】
次に、ＣＰＵ８１はステップ１０２０に進み、その時点の吸気弁通過ガス流量ｍｉの積算値Ｍｉにプログラム実行間隔時間Δｔと吸気弁通過ガス流量ｍｉの積値（Δｔ・ｍｉ）を加えて新たな吸気弁通過ガス流量ｍｉの積算値Ｍｉを求め、ステップ１０２２にて排気弁モデルＭ６によって排気弁通過ガス流量ｍｅを算出する。この場合、吸気ガス圧力Ｐｅは筒内ガス圧力Ｐｃ以上であるから、排気ポート３４からシリンダ２１へとガスが流れ込む。即ち、（５２）式又は（５４）式が使用されて排気弁通過ガス流量ｍｅが求められる。従って、排気弁通過ガス流量ｍｅは正の値となる。
【０１７５】
次いで、ＣＰＵ８１はステップ１０２４にて吸気弁通過ガス流量ｍｉの積算値Ｍｉが初期値「０」より大きいか否か（正の値であるか否か）判定する。この積算値Ｍｉが負の値であれば、シリンダ２１から吸気通路にガスが吹き返されているか、又は、その吹き返されたガスがシリンダ２１内に再吸入されている段階にあることを示す。他方、積算値Ｍｉが正の値であれば、吸気通路に吹き返されたガスが総べてシリンダ２１内に再吸入され、新気がシリンダ２１内に吸入されている段階にあることを示す。即ち、ステップ１０２４は、シリンダ２１内に新気が吸入されているか否かを判定する吸気弁通過ガス種判定手段を構成している。
【０１７６】
現時点では、吸気弁通過ガス流量ｍｉは負の値であり、従って、吸気弁通過ガス流量ｍｉの積算値Ｍｉは負の値となっているから、ＣＰＵ８１はステップ１０２４にて「Ｎｏ」と判定し、ステップ１０２６にて吸気弁通過ガス温度Ｔｉを排気ガス温度Ｔｅと等しい値に設定する。
【０１７７】
その後、ＣＰＵ８１はステップ１０２８にてシリンダモデルＭ５により筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃを算出し、ステップ１０３０にてその時点の排気弁通過ガス流量ｍｅの積算値Ｍｅにプログラム実行間隔時間Δｔと排気弁通過ガス流量ｍｅの積値（Δｔ・ｍｅ）を加えて新たな排気弁通過ガス流量ｍｅの積算値Ｍｅを求める。なお、排気弁通過ガス流量ｍｅの積算値Ｍｅを、以下、「筒内残留ガス量」と称呼する。
【０１７８】
そして、ＣＰＵ８１はステップ１０３２にて吸気弁上流ガス圧力Ｐｉをその時点で吸気管モデルＭ３により求められている吸気管内空気圧力Ｐｍと等しい値に設定し、ステップ１０３４にてその時点の排気ガス量Ｇｅに基づいて上記ステップ１０１２と同様に排気ガス圧力Ｐｅを求める。次いで、ＣＰＵ８１は、ステップ１０３６にてその時点の吸入ガス量Ｇａと燃料量Ｇｆとエンジン回転速度ＮＥとテーブルＭａｐＴｅとからステップ１０１４と同様に排気ガス温度Ｔｅ（＝ＭａｐＴｅ（Ｇａ，Ｇｆ，ＮＥ））を求め、ステップ１０９５に進んで本ルーチンを一旦終了する。
【０１７９】
その後、プログラム実行間隔時間Δｔが経過すると、ＣＰＵ８１は再びステップ１０００から本ルーチンの処理を開始する。このとき、吸気弁３２は開弁している。従って、ＣＰＵ８１はステップ１００２にて「Ｎｏ」と判定するとともに、吸気弁３２が開弁状態にあるか否か（開弁しているか否か）を判定するステップ１００４にて「Ｙｅｓ」と判定し、ステップ１０１８〜ステップ１０２４に進む。
【０１８０】
この時点では、吸気弁３２が開弁してから十分な時間が経過していないから、シリンダ２１から吸気ポート３１へガスが依然として逆流している。従って、吸気弁通過ガス流量ｍｉは負であって、吸気弁通過ガス流量ｍｉの積算値Ｍｉは負であるから、ＣＰＵ８１はステップ１０２４にて「Ｎｏ」と判定してステップ１０２６を実行し、その後、ステップ１０２８〜ステップ１０３６を実行して本ルーチンを一旦終了する。以降、以上に述べた処理が繰り返し実行される。
【０１８１】
その後、時間が経過して図１１に示した時刻ｔ１（クランク角が吸気上死点ＴＤＣよりも僅かだけ後のクランク角θ１となる時刻）になると、排気弁３５は実質的に閉弁する。また、この時点では、ピストン２２が下降して燃焼室２５の容積が増大している。これにより、吸気ポート３１内に吹き返された既燃ガスは吸気弁３２の周囲を介してシリンダ２１内に吸入され始める。即ち、吸気弁上流ガス圧力Ｐｉは筒内ガス圧力Ｐｃより大きくなるので、ステップ１０１８において上記（３１）式又は上記（３３）式が使用されて吸気弁通過ガス流量ｍｉが求められる。
【０１８２】
このとき、前回の本ルーチン実行時におけるステップ１０２６により、吸気弁通過ガス温度Ｔｉは排気ガス温度Ｔｅと等しい値に設定されているから、（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉには排気ガス温度Ｔｅと等しい値が使用される。この結果、吸気弁通過ガス流量ｍｉは実際の値により近い値となる。また、吸気弁通過ガス流量ｍｉは正の値となるから、吸気弁通過ガス流量ｍｉの積算値Ｍｉは時刻ｔ１以降において増大する。
【０１８３】
しかしながら、吸気弁通過ガス流量ｍｉの積算値Ｍｉは依然として負の値であるから、ＣＰＵ８１はステップ１０２４にて「Ｎｏ」と判定してステップ１０２６を実行する。従って、次回の本ルーチン実行時においても上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉは排気ガス温度Ｔｅとなる。
【０１８４】
そして、かかる状態が継続すると、吸気弁通過ガス流量ｍｉの積算値Ｍｉは時刻ｔ４’にて「０」となる。本実施形態では、この時刻ｔ４’が、吸気通路に吹き返された既燃ガスの総べてがシリンダ２１内に吸入され、新気がシリンダ２１内に吸入され始める時点に相当する。
【０１８５】
この時刻ｔ４’の直後に本ルーチンが実行されると、ＣＰＵ８１はステップ１０２４に進んだとき同ステップ１０２４にて「Ｙｅｓ」と判定してステップ１０３８に進み、吸気弁通過ガス温度Ｔｉを吸気管モデルＭ３が推定している吸気管内空気温度Ｔｍと等しい値に設定するようになる。
【０１８６】
これにより、ＣＰＵ８１が次回の本ルーチンの実行時にステップ１０１８に進んだとき、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉが吸気管内空気温度Ｔｍとされて吸気弁通過ガス流量ｍｉが計算される。このように、時刻ｔ４’以降においては新気がシリンダ２１内に吸入されているはずであるから、吸気弁通過ガス温度Ｔｉを吸気管内空気温度Ｔｍと設定することにより吸気弁通過ガス流量ｍｉがより精度良く計算される。
【０１８７】
その後、吸気弁３２は、クランク角が閉弁クランク角θ５となると、開弁状態から閉弁状態へと移行する（図１１の時刻ｔ５を参照。）。このとき、ＣＰＵ８１は、ステップ１００２及びステップ１００４にて「Ｎｏ」と判定し、ステップ１００６にて吸気弁３２が開弁状態から閉弁状態へと移行したか否かを判定する。この判定結果は「Ｙｅｓ」となるので、ＣＰＵ８１はステップ１００６からステップ１０４０へと進み、その時点までに計算されている吸気弁３２が開弁したときからの吸気弁通過ガス流量ｍｉの積算値Ｍｉを筒内吸入新気量Ｍａとして格納する。
【０１８８】
そして、ＣＰＵ８１は、ステップ１０４２にて燃料噴射量ｆｃを筒内吸入新気量Ｍａに基づいて求め、続くステップ１０４４にて筒内吸入新気量Ｍａ、筒内残留ガス量Ｍｅ及びエンジン回転速度ＮＥと点火時期マップＭａｐθｉｇとから、点火時期θｉｇ（＝Ｍａｐθｉｇ（Ｍａ，Ｍｅ，ＮＥ））を決定し、ステップ１０９５に進んで本ルーチンを一旦終了する。このように、本燃料噴射量制御装置は、筒内吸入新気量Ｍａ及び筒内残留ガス量Ｍｅに基づいて燃料噴射量ｆｃや点火時期θｉｇ等のエンジン制御パラメータを決定する。
【０１８９】
以上、説明したように、第１実施形態の筒内新気量推定装置は、吸気弁３２が排気行程後半において（クランク角が吸気上死点ＴＤＣより前のクランク角θ０となったとき）開弁し、その時点では排気弁３５が開弁しているバルブオーバーラップを有する内燃機関１０に適用されている。なお、排気弁３５は、吸気行程開始直後において（クランク角が吸気上死点ＴＤＣの直後のクランク角θ２となったとき）閉弁する。
【０１９０】
また、第１実施形態の筒内吸入新気量推定装置は、内燃機関１０の吸気弁３２の周囲を通過して吸気通路（吸気ポート３１及び吸気管４１）からシリンダ２１内に吸入されているガスが新気であるかシリンダ２１内から吸気通路へ吹き返された既燃ガスであるかを判定する吸気弁通過ガス種判定手段（ステップ１０２４）を備えている。
【０１９１】
更に、この筒内吸入新気量推定装置は、シリンダ２１内に吸入されているガスが新気であると判定されているとき（積算値Ｍｉ＞０）、新気の状態パラメータを使用して吸気弁３２の周囲を通過するガスの流量である吸気弁通過ガス流量ｍｉを算出する。即ち、この筒内吸入新気量推定装置は、シリンダ２１内に吸入されているガスが新気であると判定されているとき、新気の状態パラメータである吸気管内空気温度Ｔｍを吸気弁通過ガス温度Ｔｉの値として設定し、この吸気弁通過ガス温度Ｔｉを使用して吸気弁通過ガス流量ｍｉを算出する吸気弁通過ガス流量算出手段を備える（ステップ１０３８、ステップ１０１８及びステップ１０２８）。なお、新気の温度Ｔｍは新気温取得手段として機能する吸気管モデルＭ３によって取得される。
【０１９２】
この吸気弁通過ガス流量算出手段は、シリンダ２１内に吸入されているガスが吸気通路へ吹き返された既燃ガスであると判定されているとき（積算値Ｍｉ＜０又はＭｉ≦０）、その既燃ガスの状態パラメータを使用して吸気弁通過ガス流量ｍｉを推定する。即ち、吸気弁通過ガス流量算出手段は、既燃ガスの状態パラメータである排気ガス温度Ｔｅを使用して吸気弁通過ガス流量ｍｉを算出する（ステップ１０２６、ステップ１０１８及びステップ１０２８）。なお、排気ガスの温度は排気ガス温取得手段として機能するステップ１０１４及びステップ１０３６又は排気温度センサにより取得される。
【０１９３】
更に、この筒内吸入新気量推定装置は、算出された吸気弁通過ガス流量ｍｉに基づいてシリンダ２１内に吸入される新気の量である筒内吸入新気量Ｍａを推定する筒内吸入新気量推定手段（ステップ１０２０及びステップ１０４０）を備えている。
【０１９４】
このように、上記筒内吸入新気量推定装置においては、シリンダ２１に吸入されるガスの種類（新気か吹き返された既燃ガスか）に応じて適切な吸気弁通過ガス温度Ｔｉが使用されて吸気弁通過ガス流量ｍｉが算出されるので、吸気弁通過ガス流量ｍｉに基づいて算出される筒内吸入新気量Ｍａの推定精度を向上することができる。
【０１９５】
また、前記吸気弁通過ガス流量算出手段は、吸気通路からシリンダ２１内にガスが吸入されているときに正及び負の何れか一方の符号付きガス流量（本例においては、正の符号を有するガス流量）を吸気弁通過ガス流量ｍｉとして所定時間Δｔの経過毎に算出するように構成されている（ステップ１０１８、（３１）式及び（３３）式）。更に、吸気弁通過ガス流量算出手段は、シリンダ２１内から吸気通路にガスが吹き返されているときに正及び負の何れか他方の符号付きガス流量（本例においては負の符号を有するガス流量）を、吸気弁通過ガス流量ｍｉとして所定時間Δｔの経過毎に算出するように構成されている（ステップ１０１８、（３２）式及び（３４）式）。
【０１９６】
一方、前述した吸気弁通過ガス種判定手段は、算出された吸気弁通過ガス流量ｍｉを吸気弁３２の開弁時から積算して積算値Ｍｉを求める（ステップ１０２０）。吸気弁３２の開弁時において設定される積算値Ｍｉの初期値は、本例では０である。そして、吸気弁通過ガス種判定手段は、この積算値Ｍｉの初期値（即ち「０」）より大きいか否かに応じてシリンダ２１内に吸入されているガスが「新気」であるか「それまでに吸気通路に吹き返された既燃ガス」であるかを判定するように構成されている（ステップ１０２４）。
【０１９７】
従って、この筒内吸入新気量推定装置は、簡単な構成により吸気弁３２の周囲を通過するガスの種類を判別することができる。
【０１９８】
次に、本発明の第２実施形態に係る燃料噴射量制御装置（筒内吸入新気量推定装置）について説明する。この燃料噴射量制御装置は、第１実施形態に係る燃料噴射量制御装置のＣＰＵ８１が実行する図１０のステップ１０２４、ステップ１０２６及びステップ１０３８を図１２に示したステップ１２０２〜ステップ１２０６に置換した点のみにおいて、同第１実施形態に係る燃料噴射量制御装置と相違している。従って、以下、係る相違点を主として説明する。
【０１９９】
第２実施形態に係る燃料噴射量制御装置のＣＰＵ８１は、吸気弁３２が開弁しているときにステップ１０２２の処理を実行すると、図１２のステップ１２０２に進んで吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」より大きいか否かを判定する。そして、吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」より小さい場合、ＣＰＵ８１はステップ１２０２にて「Ｎｏ」と判定してステップ１２０４に進み、吸気弁通過ガス温度ＴｉをシリンダモデルＭ５が推定している筒内ガス温度Ｔｃと等しい値に設定する。これは、吸気ポート３１に吹き返された既燃ガスの温度はシリンダ２１内のガスの温度と近似した温度になっているであろうという知見に基づいている。
【０２００】
そして、ＣＰＵ８１はステップ１２０４から図１０のステップ１０２８以降のステップに進む。この結果、吸気ポート３１に吹き返された既燃ガスがシリンダ２１内に吸入されるとき（図１１の時刻ｔ１〜ｔ４’に相当する期間）、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉに筒内ガス温度Ｔｃが適用されて吸気弁通過ガス流量ｍｉが算出される。これにより、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉが実際の吸気弁通過ガス温度に近い値となるので、吸気弁通過ガス流量ｍｉが精度良く計算される。
【０２０１】
一方、吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」より大きくなると、ＣＰＵ８１はステップ１２０２に進んだとき、同ステップ１２０２にて「Ｙｅｓ」と判定してステップ１２０６に進み、吸気弁通過ガス温度Ｔｉを吸気管モデルＭ３が推定している吸気管内空気温度Ｔｍに設定し、ステップ１０２８以降のステップを実行するようになる。
【０２０２】
これにより、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉには吸気管内空気温度Ｔｍが適用されるようになる。吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」より大きくなった時点以降ではシリンダ２１内に新気が吸入されているはずであるから、吸気弁通過ガス温度Ｔｉを吸気管内空気温度Ｔｍと等しい値に設定することにより吸気弁通過ガス流量ｍｉがより精度良く計算される。
【０２０３】
以上、説明したように、第２実施形態の筒内新気量推定装置は、第１実施形態と同様にバルブオーバーラップを有する内燃機関１０に適用されている。また、第２実施形態の筒内吸入新気量推定装置は、第１実施形態と同様に、吸気弁通過ガス種判定手段（ステップ１２０２）を備えている。
【０２０４】
更に、この筒内吸入新気量推定装置は、シリンダ２１内に吸入されているガスが新気であると判定されているとき（積算値Ｍｉ＞０）、新気の状態パラメータを使用して吸気弁３２の周囲を通過するガスの流量である吸気弁通過ガス流量ｍｉを算出する。即ち、この筒内吸入新気量推定装置は、シリンダ２１内に吸入されているガスが新気であると判定されているとき、新気の状態パラメータである吸気管内空気温度Ｔｍを吸気弁通過ガス温度Ｔｉの値として設定し、この吸気弁通過ガス温度Ｔｉを使用して吸気弁通過ガス流量ｍｉを算出する吸気弁通過ガス流量算出手段（ステップ１２０６、ステップ１０１８及びステップ１０２８）を備える。なお、新気の温度は新気温取得手段として機能する吸気管モデルＭ３によって取得される。
【０２０５】
この吸気弁通過ガス流量算出手段は、シリンダ２１内に吸入されているガスが吸気通路へ吹き返された既燃ガスであると判定されているとき（積算値Ｍｉ＜０又はＭｉ≦０）、その既燃ガスの状態パラメータを使用して吸気弁通過ガス流量ｍｉを推定する。即ち、この吸気弁通過ガス流量算出手段は、既燃ガスの状態パラメータである筒内ガス温度Ｔｃを使用して吸気弁通過ガス流量ｍｉを算出する（ステップ１２０４、ステップ１０１８及びステップ１０２８）。なお、筒内ガス温度Ｔｃは筒内ガス温取得手段として機能するシリンダモデルＭ５（ステップ１０２８）によって取得される。
【０２０６】
更に、この筒内吸入新気量推定装置は、算出された吸気弁通過ガス流量ｍｉに基づいてシリンダ２１内に吸入される新気の量である筒内吸入新気量Ｍａを推定する筒内吸入新気量推定手段（ステップ１０２０及びステップ１０４０）を備えている。
【０２０７】
このように、シリンダ２１に吸入されるガスの種類（新気か吹き返された既燃ガスか）に応じて適切な吸気弁通過ガス温度Ｔｉが使用されて吸気弁通過ガス流量ｍｉが算出されるので、吸気弁通過ガス流量ｍｉが精度良く求められる。従って、吸気弁通過ガス流量ｍｉに基づいて算出される筒内吸入新気量Ｍａの推定精度を向上することができる。
【０２０８】
次に、本発明の第３実施形態に係る燃料噴射量制御装置（筒内吸入新気量推定装置）について説明する。第３実施形態に係る燃料噴射量制御装置は、吸気ポート３１内に吹き返されたガスの温度（（以下、「吹き返しガス温度」と称呼する。）Ｔｐを、同吹き返されたガスのエネルギー（エネルギーの積算値）と同吹き返されたガスの量とから求めておき、同吹き返されたガスがシリンダ２１内に再吸入されるとき、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉに計算された吹き返しガス温度Ｔｐを適用して吸気弁通過ガス流量ｍｉを算出する点で、第１及び第２実施形態の燃料噴射量制御装置と相違している。
【０２０９】
より具体的に述べると、第３実施形態のＣＰＵ８１は、シリンダ２１から吸気ポート３１へガスが吹き返されている期間、その吹き返しガスのエネルギーを積分（積算）した値ＳＥを得るとともに、吹き返しガスの総量ＳＣを求めておく。そして、ＣＰＵ８１は、吹き返し期間が終了したとき（吹き返されたガスがシリンダ２１内に吸入され始めるとき）、下記（５６）式に示した式に従って吹き返しガス温度Ｔｐを求める。そして、ＣＰＵ８１は、新気がシリンダ２１内に吸入され始めるまで（吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」となるまで）、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉに吹き返しガス温度Ｔｐを適用して吸気弁通過ガス流量ｍｉを求める。
【０２１０】
【数５６】

【０２１１】
次に、第３実施形態に係る燃料噴射量制御装置の作動について説明する。この装置は、図１３のフローチャートに示したように、図１０に示したステップ１０１６とステップ１０１８の間に追加されたステップ１３０２を実行する点と、図１０のステップ１０２４、ステップ１０２６及びステップ１０３８に代わるステップ１３０４〜ステップ１３１６を実行する点のみにおいて、第１実施形態に係る燃料噴射量制御装置と相違している。従って、以下、係る相違点を主として説明する。
【０２１２】
第３実施形態に係る燃料噴射量制御装置のＣＰＵ８１は、吸気弁３２が閉弁状態から開弁状態へと変化したとき、ステップ１３００に続くステップ１００２〜１０１６を実行した後、ステップ１３０２にて吹き返しガスのエネルギー積分値ＳＥ及び吹き返しガスの総量ＳＣを共に「０」に設定する。そして、ステップ１０１８〜ステップ１０２２を実行した後、ステップ１３０４にて吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」より大きいか否かを判定する。
【０２１３】
この場合、吸気弁通過ガス流量ｍｉの積算値Ｍｉは「０」より小さいから、ＣＰＵ８１はステップ１３０４にて「Ｎｏ」と判定してステップ１３０６に進み、ステップ１０１０にて吸気管内空気圧力Ｐｍと等しい値に設定されている吸気弁上流ガス圧力Ｐｉが筒内ガス圧力Ｐｃよりも大きいか否かを判定する。即ち、ステップ１３０６では、吸気通路からシリンダ２１内にガスが流れ込んでいるか、シリンダ２１から吸気通路にガスが吹き返されているかが判定される。
【０２１４】
現時点は、吸気弁３２が開弁した直後であるから、吸気弁上流ガス圧力Ｐｉは筒内ガス圧力Ｐｃよりも小さいので、ＣＰＵ８１はステップ１３０６にて「Ｙｅｓ」と判定してステップ１３０８に進み、その時点の吹き返しガスのエネルギーの積分値ＳＥに、前回の本ルーチン実行時から今回の本ルーチン実行時までの吹き返しガスのエネルギーに相当する値Δｔ・ｍｉ・Ｃｐ・Ｔｃを加えて新たな吹き返しガスのエネルギーの積分値ＳＥを得る。
【０２１５】
次いで、ＣＰＵ８１は、ステップ１３１０にてその時点の吹き返しガスの総量ＳＣに、前回の本ルーチン実行時から今回の本ルーチン実行時までの吹き返しガス量に相当する値Δｔ・ｍｉを加えて新たな吹き返しガスの総量ＳＣを得る。そして、ＣＰＵ８１は、続くステップ１３１２にて同ステップ１３１２に示した式（Ｔｐ＝ＳＥ／（Ｃｐ・ＳＣ））にしたがって吹き返しガス温度Ｔｐを求めた後、ステップ１０２８以降に進む。
【０２１６】
次に、ＣＰＵ８１が本ルーチンの処理を開始すると、ＣＰＵ８１はステップ１００２及びステップ１００４を経てステップ１０１８〜ステップ１０２２及びステップ１３０４へ進む。この時点は、吸気弁３２が開弁してから僅かな時間しか経過していないので、吸気弁通過ガス流量ｍｉの積算値Ｍｉは負の値であり、且つ、吸気弁上流ガス圧力Ｐｉは筒内ガス圧力Ｐｃよりも小さい。従って、ＣＰＵ８１はステップ１３０４及びステップ１３０６〜ステップ１３１２へと進み吹き返しガス温度Ｔｐを更新する。
【０２１７】
このように、吸気弁通過ガス流量ｍｉの積算値Ｍｉが負の値であり、且つ、吸気弁上流ガス圧力Ｐｉが筒内ガス圧力Ｐｃよりも小さいとき、吹き返しガス温度Ｔｐが更新されて行く。なお、この場合の吸気弁通過ガス流量ｍｉは、（３２）式又は（３４）式により求められる。
【０２１８】
その後、排気弁３５が閉弁するとともに燃焼室２５の容積が増大するようになるので、吸気弁上流ガス圧力Ｐｉは筒内ガス圧力Ｐｃよりも大きくなり、吹き返しガスがシリンダ２１内に吸入され始める（図１１の時刻ｔ１を参照。）。このとき、ＣＰＵ８１はステップ１３０４及びステップ１３０６の両ステップにて「Ｎｏ」と判定してステップ１３１４に進み、吸気弁通過ガス温度Ｔｉをステップ１３１２の計算により求めた吹き返しガス温度Ｔｐに設定する。
【０２１９】
この結果、吸気ポート３１に吹き返された既燃ガスがシリンダ２１内に吸入されるとき（図１１の時刻ｔ１〜ｔ４’に相当する期間）、ＣＰＵ８１はステップ１０１８において上記（３１）式又は上記（３３）式の吸気弁通過ガス温度Ｔｉに吹き返しガス温度Ｔｐを適用し吸気弁通過ガス流量ｍｉを算出する。これにより、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉが実際の吸気弁通過ガス温度に極めて近い値となるので、吸気弁通過ガス流量ｍｉが精度良く計算される。
【０２２０】
その後、吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」より大きくなると、ＣＰＵ８１はステップ１３０４に進んだとき、同ステップ１３０４にて「Ｙｅｓ」と判定してステップ１３１６に進み、吸気弁通過ガス温度Ｔｉを吸気管モデルＭ３が推定している吸気管内空気温度Ｔｍに設定し、ステップ１０２８以降のステップを実行するようになる。
【０２２１】
これにより、上記（３１）式又は上記（３３）式における吸気弁通過ガス温度Ｔｉには吸気管内空気温度Ｔｍが適用されるようになる。吸気弁通過ガス流量ｍｉの積算値Ｍｉが「０」より大きくなった時点以降ではシリンダ２１内に新気が吸入されているはずであるから、吸気弁通過ガス温度Ｔｉを吸気管内空気温度Ｔｍと等しい値に設定することにより吸気弁通過ガス流量ｍｉがより精度良く計算される。
【０２２２】
以上、説明したように、第３実施形態の筒内新気量推定装置は、第１実施形態と同様にバルブオーバーラップを有する内燃機関１０に適用されている。また、第３実施形態の筒内吸入新気量推定装置は、第１実施形態のステップ１０２４と同様な吸気弁通過ガス種判定手段としてのステップ１３０４を備えている。
【０２２３】
更に、この筒内吸入新気量推定装置は、シリンダ２１内に吸入されているガスが新気であると判定されているとき（積算値Ｍｉ＞０）、新気の状態パラメータを使用して吸気弁３２の周囲を通過するガスの流量である吸気弁通過ガス流量ｍｉを算出する。即ち、この筒内吸入新気量推定装置は、シリンダ２１内に吸入されているガスが新気であると判定されているとき、新気の状態パラメータである吸気管内空気温度Ｔｍを吸気弁通過ガス温度Ｔｉの値として設定し、この吸気弁通過ガス温度Ｔｉを使用して吸気弁通過ガス流量ｍｉを算出する吸気弁通過ガス流量算出手段（ステップ１３１６、ステップ１０１８及びステップ１０２８）を備える。なお、新気の温度Ｔｍは新気温取得手段として機能する吸気管モデルＭ３によって取得される。
【０２２４】
一方、この筒内吸入新気量推定装置は、吸気弁３２の周囲を通過してシリンダ２１から吸気通路へ吹き返された既燃ガスのエネルギーの総量ＳＥを取得するとともに同取得した既燃ガスのエネルギーに基づいて同既燃ガスの温度（再吸入時における吹き返しガス温度）Ｔｐを取得する吹き返しガス温取得手段（ステップ１３０６〜ステップ１３１２）を備える。実際には、吹き返しガス温取得手段は、既燃ガスが吸気通路へ吹き返されている期間において積算された既燃ガスのエネルギーＳＥを、同期間において積算された既燃ガスの流量である既燃ガス総量ＳＣと定圧比熱Ｃｐの積で除することにより既燃ガス温度Ｔｐを算出する（（５６）式を参照。）。
【０２２５】
そして、前述の吸気弁通過ガス流量算出手段は、シリンダ２１内に吸入されているガスが吸気通路へ吹き返された既燃ガスであると判定されているとき（積算値Ｍｉ＜０又はＭｉ≦０）、その既燃ガスの状態パラメータを使用して吸気弁通過ガス流量ｍｉを推定する。即ち、この吸気弁通過ガス流量算出手段は、既燃ガスの状態パラメータである前記算出された既燃ガス温度Ｔｐを使用して吸気弁通過ガス流量ｍｉを算出する（ステップ１３１４、ステップ１０１８及びステップ１０２８）。
【０２２６】
更に、この筒内吸入新気量推定装置は、算出された吸気弁通過ガス流量ｍｉに基づいてシリンダ２１内に吸入される新気の量である筒内吸入新気量Ｍａを推定する筒内吸入新気量推定手段（ステップ１０２０及びステップ１０４０）を備えている。
【０２２７】
従って、シリンダ２１に吸入されるガスの種類（新気か吹き返された既燃ガスか）に応じて適切な吸気弁通過ガス温度Ｔｉが使用されて吸気弁通過ガス流量ｍｉが算出されるので、吸気弁通過ガス流量ｍｉが精度良く求められる。その結果、吸気弁通過ガス流量ｍｉに基づいて算出される筒内吸入新気量Ｍａの推定精度を向上することができる。
【０２２８】
なお、第３実施形態においては、ステップ１３０８にて使用される定圧比熱Ｃｐが一定であるとして吹き返しガス温度Ｔｐを求めていたが、この定圧比熱Ｃｐを、例えば、吹き返しガス温度Ｔｐに基づいて定め（Ｃｐ＝ｆｕｎｃ（Ｔｐ））、このＣｐを使用しながら吹き返し温度Ｔｐを更新するように構成してもよい。この場合、ステップ１３１０をＳＣ＝ＳＣ＋Δｔ・ｍｉ・Ｃｐとして、ステップ１３１２をＴｐ＝ＳＥ／ＳＣとすればよい。
【０２２９】
次に、本発明による燃料噴射量制御装置の第１変形例について説明する。上述した各実施形態の燃料噴射量制御装置は、例えば吸気弁モデルＭ４を表す（３１）式〜（３５）式、シリンダモデルＭ５を表す（３７）式及び排気弁モデルＭ６を表す（５２）式〜（５５）式等において、比熱比κｃ，κｉ及びκｅを一定の値κとして取り扱っていた。
【０２３０】
しかしながら、実際の比熱比はガス組成及びガス温度等により変化する。特に内燃機関においては、数百度の既燃ガスで満たされているシリンダ２１内に数十度の新気が導入されるので、その前後でシリンダ２１内のガスの比熱比（筒内ガス比熱比）κｃは大きく変化する。更に、シリンダ２１内には残留ガス（既燃ガス）が存在した状態で新気が吸入されるから、残留ガスと新気の比率によってシリンダ２１内のガス組成が変動し、従って、筒内ガス比熱比κｃは大きく変動する。
【０２３１】
そこで、第１変形例に係る燃料噴射量制御装置は、これらの比熱比κｃ，κｉ及びκｅを推定し、例えば、推定した比熱比κｃ，κｉ及びκｅを（３１）式〜（３５）式、（３７）式及び（５２）式〜（５５）式等に適用することで吸気弁通過ガス流量ｍｉ、排気弁通過ガス流量ｍｅ、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃ等の各値の計算の精度を一層向上させる。この比熱比κｉは新気の状態パラメータであり、比熱比κｅは既燃ガスの状態パラメータと云うことができる。
【０２３２】
（第１変形例の比熱比算出原理）
先ず、第１変形例における筒内ガスの比熱比κｃを算出する原理について説明する。筒内ガスの比熱比κｃは、下記（５７）式により計算される。
【０２３３】
【数５７】

【０２３４】
（５７）式において、Ｍｉは吸気弁通過ガス流量ｍｉを積分することにより得られる。新気の比熱比κｉはガス温度Ｔと新気の比熱比κｉとの関係を規定するテーブルＭａｐκｉ及び筒内ガス温度Ｔｃから求めることができる（κｉ＝Ｍａｐκｉ（Ｔｃ））。Ｒｉは新気の気体定数であり既知である。
【０２３５】
Ｍｅは排気弁通過ガス流量ｍｅを積分することにより得られる。既燃ガスの比熱比κｅはガス温度Ｔと既燃ガスの比熱比κｅとの関係を規定するテーブルＭａｐκｅ及び筒内ガス温度Ｔｃから求めることができる（κｅ＝Ｍａｐκｅ（Ｔｃ））。Ｒｅは既燃ガスの気体定数であり既知である。従って、（５７）式の右辺の各変数は総べて得ることができるので、同（５７）式に基づいて筒内ガスの比熱比κｃを求めることができる。
【０２３６】
次に、（５７）式の導出過程について説明する。ガスの内部エネルギーをｕ、エンタルピーをｈとするとともに、成分ｊのガスの単位質量あたりの内部エネルギー及びエンタルピーをそれぞれｕｊ及びｈｊとする。また、ガス全体に対する成分ｊのガスの質量割合をｇｊとすると、下記（５８）式及び下記（５９）式が成立する。
【０２３７】
【数５８】

【０２３８】
【数５９】

【０２３９】
従って、（５８）式及び（５９）式から、それぞれ（６０）式及び（６１）式が導かれる。
【０２４０】
【数６０】

【０２４１】
【数６１】

【０２４２】
これらの関係から、筒内ガスの定容比熱Ｃｖｃ及び定圧比熱Ｃｐｃについて、それぞれ下記（６２）式及び下記（６３）式が成立する。
【０２４３】
【数６２】

【０２４４】
【数６３】

【０２４５】
従って、筒内ガスの比熱比κｃについて下記（６４）式が成立する。
【０２４６】
【数６４】

【０２４７】
（６４）式に上記（５０）式及び上記（５１）式を適用すると下記の（６５）式が得られる。
【０２４８】
【数６５】

【０２４９】
いま、ガスの成分を、新気と既燃ガスの２成分と仮定する。この仮定は極めて現実的で実際的な仮定である。このとき、筒内ガスの比熱比κｃは上記（６５）式から下記（６６）式のように表される。
【０２５０】
【数６６】

【０２５１】
一方、新気の質量割合ｇｉ及び既燃ガスの質量割合ｇｅは、それぞれ下記（６７）式及び下記（６８）式により表される。
【０２５２】
【数６７】

【０２５３】
【数６８】

【０２５４】
（６７）式及び（６８）式を（６６）式に代入すると、以下の（６９）式のようにして上記（５７）式が得られる。
【０２５５】
【数６９】

【０２５６】
以上が、（５７）式の導出過程である。次に、この第１変形例の燃料噴射量制御装置の実際の作動について図１４に示したフローチャートを参照しながら説明する。この燃料噴射量制御装置のＣＰＵ８１は、第１実施形態のＣＰＵ８１が実行する図１０のステップ１０３０とステップ１０３２の間に、図１４に示したステップを追加したルーチンを実行する点において、同第１実施形態に係る燃料噴射量制御装置と相違している。従って、以下、係る相違点を主として説明する。
【０２５７】
ＣＰＵ８１は、図１０のステップ１０３０を実行すると図１４のステップ１４０２に進み、前回の特定の気筒の燃焼行程がフューエルカット状態であったか否かを判定する。いま、前回の特定の気筒の燃焼行程がフューエルカット状態であったとすると、ＣＰＵ８１はステップ１４０２にて「Ｙｅｓ」と判定してステップ１４０４に進み、上述したテーブルＭａｐκｉ及び吸気管空気温度Ｔｍから吸気通路内の新気の比熱比κｉ（ｉｍ）を求める。
【０２５８】
次いで、ＣＰＵ８１は、ステップ１４０６に進み、上述したテーブルＭａｐκｉ及び排気温度Ｔｅから排気通路内の排気ガスの比熱比κｅ（ｅｘ）を求める。ステップ１４０６にてテーブルＭａｐκｅではなくテーブルＭａｐκｉを使用するのは、前回の特定の気筒の燃焼行程がフューエルカット状態であったから、排気通路内には既燃ガスでなく新気が存在しているからである。
【０２５９】
次に、ＣＰＵ８１はステップ１４０８に進み、テーブルＭａｐκｉ及び筒内ガス温度Ｔｃからシリンダ２１内での新気の比熱比（筒内新気比熱比）κｉを求め、ステップ１４１０及びステップ１４１２にて筒内ガス比熱比κｃ及びシリンダ２１内での既燃ガスの比熱比（筒内既燃ガス比熱比）κｅをそれぞれ筒内新気比熱比κｉと等しい値に設定し、その後、図１０のステップ１０３２に進む。
【０２６０】
このように、前回の特定の気筒の燃焼行程がフューエルカット状態であった場合、シリンダ２１内には大気（実質的に既燃ガスを含まないので新気と言うことができる）のみが存在しているから、ＣＰＵ８１は筒内ガス比熱比κｃを筒内新気比熱比κｉと等しい値に設定する。また、前回の燃焼行程では実際に燃焼が生じていなかったので、筒内既燃ガス比熱比κｅも筒内新気比熱比κｉと等しい値に設定する。
【０２６１】
一方、前回の特定の気筒の燃焼行程がフューエルカット状態ではなかったとすると、ＣＰＵ８１はステップ１４０２にて「Ｎｏ」と判定してステップ１４１４に進み、上述したテーブルＭａｐκｉ及び筒内ガス温度Ｔｃから筒内新気比熱比κｉを求めるとともに、続くステップ１４１６にて上述したテーブルＭａｐκｅ及び筒内ガス温度Ｔｃから筒内既燃ガス比熱比κｅを求める。次いで、ＣＰＵ８１は、ステップ１４１８に進み、上記（５７）式にしたがって筒内ガスの比熱比κｃを算出する。
【０２６２】
その後、ＣＰＵ８１は、ステップ１４２０にて、上述したテーブルＭａｐκｉ及び吸気管空気温度Ｔｍから吸気通路内の新気の比熱比κｉ（ｉｍ）を求め、ステップ１４２２にて、上述したテーブルＭａｐκｅ及び排気温度Ｔｅから排気通路内の排気ガスの比熱比κｅ（ｅｘ）を求め、その後、図１０のステップ１０３２に進む。
【０２６３】
このように、ＣＰＵ８１は、前回の特定の気筒の燃焼行程がフューエルカット状態ではなかった場合、シリンダ２１内には新気と既燃ガスとが混在しているから、上述した（５７）式に基づいて筒内ガスの比熱比κｃを求める。
【０２６４】
そして、ＣＰＵ８１は、例えば、次回の本ルーチンの実行時に、ステップ１０２８（シリンダモデルＭ５を表す（３７）式）にて上記のように求められた比熱比κｃ，κｉ及びκｅを使用して筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃを求める。この結果、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃをより精度良く求めることができる。
【０２６５】
また、ＣＰＵ８１は、ステップ１０２４にて既燃ガスが吸気弁３２の周囲を通過していると判定したときであって、同既燃ガスがシリンダ２１内に吸入されていると判定したとき（Ｐｉ＞Ｐｃ）、（３１）式（（３３）式及び（３５）式）で使用する比熱比κ（ガスの状態パラメータ）の値を排気ガスの比熱比κｅ（ｅｘ）に設定して吸気弁通過ガス流量ｍｉを算出する。
【０２６６】
ＣＰＵ８１は、この場合を除き、吸気弁モデルＭ４を表す（３１）式（（３３）式及び（３５）式）での比熱比κに新気比熱比κｉ（ｉｍ）を代入して吸気弁通過ガス流量ｍｉを求めるとともに、（３２）式（（３４）式及び（３５）式）での比熱比κに排気ガス比熱比κｅ（ｅｘ）を代入して吸気弁通過ガス流量ｍｉを求める。
【０２６７】
同様に、ＣＰＵ８１は、排気弁モデルＭ６を表す（５２）式（（５４）式及び（３５）式）での比熱比κに排気ガスの比熱比κｅ（ｅｘ）を代入して排気弁通過ガス流量ｍｅを求めるとともに、（５３）式（（５５）式及び（３５）式）での比熱比κに筒内ガス比熱比κｃを代入して排気弁通過ガス流量ｍｅを求める。
【０２６８】
以上、説明したように、第１変形例の燃料噴射量制御装置は、吸気通路内の新気ガスの比熱比κｉ（ｉｍ）を取得する新気ガス比熱比取得手段（ステップ１４０４及びステップ１４２０）と、排気通路内の排気ガスの比熱比κｅ（ｅｘ）を取得する排気ガス比熱比取得手段（ステップ１４０６及び１４２２）と、筒内ガスの比熱比κｃを取得する筒内ガス比熱比取得手段（ステップ１４１０及びステップ１４１８）とを備えていて、これらの比熱比をガスの状態パラメータとして各種の計算に使用する。従って、より精度良く各値を求めることができる。
【０２６９】
次に、本発明による燃料噴射量制御装置の第２変形例について説明する。第１変形の燃料噴射量制御装置は、実際に計算により求まる新気の質量割合ｇｉ及び既燃ガスの質量割合ｇｅに基づいて筒内ガスの比熱比κｃを求めていた。これに対し、第２変形例に係る燃料噴射量制御装置は、筒内ガスの成分の質量割合ｇｉを運転状態別に（具体的には、フューエルカットであったか否かに応じて）予め記憶しておき、これを実際の運転状態に応じて読み出す。そして、燃料噴射量制御装置は、読み出した重量割合ｇｉ、予め記憶しておいた新気の比熱比κｉ及び既燃ガスの比熱比κｅを上記（６５）式に適用することで比熱比κｃを算出する。
【０２７０】
以下、第２変形例の燃料噴射量制御装置の実際の作動について図１５に示したフローチャートを参照しながら説明する。この燃料噴射量制御装置のＣＰＵ８１は、第１実施形態のＣＰＵ８１が実行する図１０のステップ１０３０とステップ１０３２の間に、図１５に示したステップを実行する点において、同第１実施形態に係る燃料噴射量制御装置と相違している。従って、以下、係る相違点を主として説明する。
【０２７１】
ＣＰＵ８１は、図１０のステップ１０３０を実行すると図１５のステップ１５０２に進み、前回の特定の気筒の燃焼行程がフューエルカット状態であったか否かを判定する。いま、前回の特定の気筒の燃焼行程がフューエルカット状態であったとすると、ＣＰＵ８１はステップ１５０２にて「Ｙｅｓ」と判定してステップ１５０４に進み、標準的な大気中の成分組成より各ガスｊの質量割合ｇｊをＲＯＭ８２から読み出す。
【０２７２】
次いで、ＣＰＵ８１は、ステップ１５０６にて各ガスｊのシリンダ２１内における比熱比κｊを、ガスｊの温度Ｔと比熱比κｊとの関係を規定したテーブルＭａｐκｊ及びその時点で求められている筒内ガス温度Ｔｃから求める（κｊ＝Ｍａｐκｊ（Ｔｃ））。そして、ＣＰＵ８１はステップ１５０８に進み、ＲＯＭ８２に記憶してある各ガスｊの気体定数Ｒｊ、上記ステップ１５０６にて求めた各ガスｊのシリンダ２１内における比熱比κｊ、上記ステップ１５０４にて読み出した各ガスｊの質量割合ｇｊ、及び上記（６５）式を用いて筒内ガスの比熱比κｃを算出し、図１０のステップ１０３２に進む。
【０２７３】
一方、前回の特定の気筒の燃焼行程がフューエルカット状態ではなかったとすると、ＣＰＵ８１はステップ１５０２にて「Ｎｏ」と判定してステップ１５１０に進み、既燃ガス中の各ガスｊの質量割合ｇｊをＲＯＭ８２から読み出す。この質量割合ｇｊは、標準的な燃料が完全燃焼した場合に生成される既燃ガス中の成分組成を予め調査することにより得た値である。その後、ＣＰＵ８１はステップ１５０６及びステップ１５０８を実行して筒内ガスの比熱比κｃを算出し、ステップ１０３２へと進む。
【０２７４】
そして、ＣＰＵ８１は、例えば、次回の本ルーチンの実行時に、ステップ１０２８にて上記のように求められた比熱比κｃを使用して筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃを求める。この結果、筒内ガス圧力Ｐｃ及び筒内ガス温度Ｔｃをより精度良く求めることができる。
【０２７５】
以上、説明したように、第２変形例の燃料噴射量制御装置は、筒内ガスの比熱比κｃを取得する筒内ガス比熱比取得手段（図１５を参照。）を備えていて、求めた比熱比κｃを筒内圧力Ｐｃ等の計算に使用する。従って、より精度良く各値を求めることができる。
【０２７６】
以上、説明したように、本発明による筒内吸入新気量推定装置の各実施形態及び各変形例によれば、筒内吸入新気量を精度良く求めることができるので、内燃機関１０を適切に制御することができる。
【０２７７】
なお、本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記シリンダモデルＭ５の（３７）式において伝達熱流量Ｑｗ’は無視されていたが、定義式である下記（７０）式及び下記（７１）式により求めてもよい。（７０）式及び（７１）式において、Ａは燃焼室を構成するシリンダ壁の面積、ｈｗはシリンダ壁面での熱伝達率、Ｔｃは筒内ガス温度、Ｔｗは冷却水温ＴＨＷに基づいて定められるシリンダ壁温度である。ｄはシリンダボア径、Ｋ１，Ｋ２は定数、Ｖｃｓｐｄはピストンの平均速度である。
【０２７８】
【数７０】
Ｑｗ’ ＝Ａ・ｈｗ・（Ｔｃ−Ｔｗ） …（７０）
【０２７９】
【数７１】
ｈｗ ＝Ｋ１・ｄ^−０．２・Ｐｃ^０．８・（Ｋ２・Ｖｃｓｐｄ）^０．８・Ｔｗ^{−０．５３} …（７１）
【図面の簡単な説明】
【図１】本発明の第１実施形態に係る内燃機関の筒内吸入新気量推定装置を含む燃料噴射量制御装置を火花点火式多気筒内燃機関に適用したシステムの概略構成図である。
【図２】図１に示した電気制御装置が筒内吸入新気量等を推定するために採用した各種モデルの接続関係を示した機能ブロック図である。
【図３】図１に示したＣＰＵが参照するアクセルペダル操作量と目標スロットルバルブ開度との関係を規定したテーブルを示す図である。
【図４】図１に示したＣＰＵが参照するスロットルバルブ開度と流量係数との関係を規定したテーブルを示す図である。
【図５】図１に示したＣＰＵが参照するスロットルバルブ開度と開口面積との関係を規定したテーブルを示す図である。
【図６】図１に示したＣＰＵが参照する、スロットルバルブ開度と、流量係数と開口面積の積値との関係を規定したテーブルを示す図である。
【図７】図１に示したＣＰＵが各値の計算時に使用する変数を説明するためシリンダ及びその近傍を概念的に示した図である。
【図８】図１に示したＣＰＵが参照する吸気弁リフト量と、流量係数と開口面積の積値との関係を規定したテーブルを示す図である。
【図９】図１に示したＣＰＵが参照する排気弁リフト量と、流量係数と開口面積の積値との関係を規定したテーブルを示す図である。
【図１０】図１に示したＣＰＵが実行するルーチンを示したフローチャートである。
【図１１】（Ａ）は時間経過に対する吸気弁及び排気弁のバルブリフト量を示したグラフであり、（Ｂ）は時間経過に対する吸気弁通過ガス流量の積算値の変化を示したグラフである。
【図１２】本発明の第２実施形態に係る筒内吸入新気量推定装置のＣＰＵが実行するルーチンを示したフローチャートである。
【図１３】本発明の第３実施形態に係る筒内吸入新気量推定装置のＣＰＵが実行するルーチンを示したフローチャートである。
【図１４】本発明の筒内吸入新気量推定装置の第１変形例に係るＣＰＵが実行するルーチンを示したフローチャートである。
【図１５】本発明の筒内吸入新気量推定装置の第２変形例に係るＣＰＵが実行するルーチンを示したフローチャートである。
【図１６】（Ａ）は時間経過に対する吸気弁及び排気弁のバルブリフト量を示したグラフであり、（Ｂ）は時間経過に対する吸気弁通過ガス流量の積算値の変化を示したグラフである。
【符号の説明】
１０…内燃機関、２１…シリンダ、２５…燃焼室、３２…吸気弁、３１…吸気ポート、３３…吸気弁制御装置、３４…排気ポート、３９…インジェクタ、４１…吸気管、５２…排気管、８０…電気制御装置、８１…ＣＰＵ。

Claims (7)

1. 内燃機関の吸気弁の周囲を通過して吸気通路からシリンダ内に吸入されているガスが新気であるか同シリンダ内から同吸気通路へ吹き返された既燃ガスであるかを判定する吸気弁通過ガス種判定手段と、
   前記シリンダ内に吸入されているガスが前記新気であると判定されているとき同新気の状態パラメータを使用して前記吸気弁の周囲を通過するガスの流量である吸気弁通過ガス流量を算出するとともに同シリンダ内に吸入されているガスが前記既燃ガスであると判定されているとき同既燃ガスの状態パラメータを使用して同吸気弁通過ガス流量を推定する吸気弁通過ガス流量算出手段と、
   前記算出された吸気弁通過ガス流量に基づいて前記シリンダ内に吸入される新気の量である筒内吸入新気量を推定する筒内吸入新気量推定手段と、
   を備えた内燃機関の筒内吸入新気量推定装置。
2. 請求項１に記載の内燃機関の筒内吸入新気量推定装置において、
   前記吸気弁通過ガス流量算出手段は、前記吸気通路から前記シリンダ内にガスが吸入されているときに前記吸気弁通過ガス流量を正及び負の何れか一方の符号付きガス流量として、前記シリンダ内から前記吸気通路にガスが吹き返されているときに前記吸気弁通過ガス流量を正及び負の何れか他方の符号付きガス流量として、所定時間の経過毎に算出するように構成され、
   前記吸気弁通過ガス種判定手段は、前記算出された吸気弁通過ガス流量の前記吸気弁の開弁時からの積算値が同積算値の初期値より大きいか否かに応じて前記シリンダ内に吸入されているガスが前記新気であるか前記既燃ガスであるかを判定するように構成された筒内吸入新気量推定装置。
3. 請求項１又は請求項２に記載の内燃機関の筒内吸入新気量推定装置であって、
   前記吸気通路内の新気の温度を取得する新気温取得手段と、
   前記内燃機関の排気通路内の排気ガスの温度を取得する排気ガス温取得手段と、
   を備えるとともに、
   前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の温度を使用するとともに前記既燃ガスの状態パラメータとして前記取得された排気ガスの温度を使用するように構成された筒内吸入新気量推定装置。
4. 請求項１又は請求項２に記載の内燃機関の筒内吸入新気量推定装置であって、
   前記吸気通路内の新気の温度を取得する新気温取得手段と、
   前記シリンダ内のガスの温度を取得する筒内ガス温取得手段と、
   を備えるとともに、
   前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の温度を使用するとともに前記既燃ガスの状態パラメータとして前記取得されたシリンダ内のガスの温度を使用するように構成された筒内吸入新気量推定装置。
5. 請求項１又は請求項２に記載の内燃機関の筒内吸入新気量推定装置であって、
   前記吸気通路内の新気の温度を取得する新気温取得手段と、
   前記吸気弁の周囲を通過して前記シリンダから前記吸気通路へ吹き返された既燃ガスのエネルギーを取得するとともに同取得した既燃ガスのエネルギーに基づいて同既燃ガスの温度を取得する吹き返しガス温取得手段と、
   を備えるとともに、
   前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の温度を使用するとともに前記既燃ガスの状態パラメータとして前記取得された既燃ガスの温度を使用するように構成された筒内吸入新気量推定装置。
6. 請求項１乃至請求項５の何れか一項に記載の内燃機関の筒内吸入新気量推定装置であって、
   前記吸気通路内の新気の比熱比を取得する新気比熱比取得手段と、
   前記排気通路内の排気ガスの比熱比を取得する排気ガス比熱比取得手段と、
   を備えるとともに、
   前記吸気弁通過ガス流量算出手段は、前記新気の状態パラメータとして前記取得された新気の比熱比を使用するとともに前記既燃ガスの状態パラメータとして前記取得された排気ガスの比熱比を使用するように構成された筒内吸入新気量推定装置。
7. 請求項１乃至請求項６の何れか一項に記載の内燃機関の筒内吸入新気量推定装置であって、
   前記吸気弁通過ガス流量算出手段は、前記吸気弁の周囲に形成されるガス通路を絞り部とみなしたときに成立する同絞り部を通過するガス流量を算出する式に基づいて前記吸気弁通過ガス流量を算出するように構成された筒内吸入新気量推定装置。

Images