(57)【特許請求の範囲】
【請求項１】
内燃機関の吸気ポート内に燃料を噴射するポートインジェクタと、
前記内燃機関の筒内に燃料を噴射する筒内インジェクタと、
前記筒内から前記吸気ポートへの吸気の吹き返し中に含まれる燃料の量を所定のパラメータに基づいて推定する吹き返し燃料量推定手段と、
前記吹き返し燃料量推定手段により推定された量に基づいて、吹き返しによる誤差を補正するための吹き返し補正分を算入して、前記ポートインジェクタから噴射するポート噴射量と前記筒内インジェクタから噴射する筒内噴射量とを算出する噴射量算出手段と、
を備え、
前記所定のパラメータは、前記ポートインジェクタと前記筒内インジェクタとの噴射割合を含み、
前記吹き返し燃料量推定手段は、前記筒内インジェクタの噴射割合が大きい場合には小さい場合に比べて吹き返し燃料量を少なく推定することを特徴とする内燃機関の制御装置。
【請求項２】
内燃機関の吸気ポート内に燃料を噴射するポートインジェクタと、
前記内燃機関の筒内に燃料を噴射する筒内インジェクタと、
前記筒内から前記吸気ポートへの吸気の吹き返し中に含まれる燃料の量を所定のパラメータに基づいて推定する吹き返し燃料量推定手段と、
前記吹き返し燃料量推定手段により推定された量に基づいて、吹き返しによる誤差を補正するための吹き返し補正分を算入して、前記ポートインジェクタから噴射するポート噴射量と前記筒内インジェクタから噴射する筒内噴射量とを算出する噴射量算出手段と、
を備え、
前記所定のパラメータは、前記筒内インジェクタの噴射開始時期を含み、
前記吹き返し燃料量推定手段は、前記筒内インジェクタの噴射開始時期が遅い場合には早い場合に比べて吹き返し燃料量を少なく推定することを特徴とする内燃機関の制御装置。
【請求項３】
前記噴射量算出手段は、前記吹き返し補正分を前記ポート噴射量と前記筒内噴射量との双方に振り分けて算入する手段を含むことを特徴とする請求項１または２記載の内燃機関の制御装置。
【請求項４】
前記噴射量算出手段は、前記吹き返し補正分を前記筒内噴射量のみに算入する手段を含むことを特徴とする請求項１または２記載の内燃機関の制御装置。
【請求項５】
前記吹き返し補正分が前記筒内噴射量のみに算入された場合に、前記筒内インジェクタの噴射開始時期を吸気弁が閉じた後の時期にする噴射時期制御手段を更に備えることを特徴とする請求項４記載の内燃機関の制御装置。
【請求項６】
前記噴射量算出手段は、前記吹き返し補正分を前記筒内噴射量のみに算入する際に、前サイクルの吹き返し補正分に対する今サイクルの吹き返し補正分の変化が判定値より大きい場合、前記吹き返し補正分の一部のみを前記筒内噴射量に算入する一部算入手段を含むことを特徴とする請求項４または５記載の内燃機関の制御装置。
【請求項７】
前記吹き返し補正分が正の値であり、かつ、前記一部算入手段により前記吹き返し補正分の一部しか前記筒内噴射量に算入されなかった場合に、前記吹き返し補正分の残部に相当する量の燃料を排気行程において前記筒内インジェクタから噴射させる排気行程噴射手段を更に備えることを特徴とする請求項６記載の内燃機関の制御装置。
【請求項８】
前記吹き返し補正分が負の値であり、かつ、前記一部算入手段により前記吹き返し補正分の一部しか前記筒内噴射量に算入されなかった場合に、前記吹き返し補正分の残部を次サイクル以降の噴射量に算入することによって次サイクル以降の噴射量を減量させる噴射量減量手段を更に備えることを特徴とする請求項６または７記載の内燃機関の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、内燃機関の制御装置に係り、特に、ポートインジェクタと筒内インジェクタとの双方を備えた内燃機関の制御装置に関する。
【背景技術】
【０００２】
吸気ポート内に燃料を噴射するポートインジェクタと、筒内に燃料を噴射する筒内インジェクタとの双方を備えた内燃機関が知られている。特開平５−２３１２２１号公報には、このような内燃機関において、空燃比を精度良く制御するべく、ポートインジェクタから噴射された燃料のうち吸気ポート内壁に付着する分や、吸気ポート内壁に付着した燃料が筒内へ流入する分を考慮して、筒内インジェクタの噴射量を補正するシステムが開示されている。
【０００３】
【特許文献１】特開平５−２３１２２１号公報
【特許文献２】特許第２７５４７４４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、上記従来の空燃比制御技術では、実際上は空燃比を必ずしも十分な精度で制御することができず、その結果、排気エミッションやドライバビリティに悪影響を及ぼす場合があった。その原因は、ポートインジェクタと筒内インジェクタとの双方を備えた内燃機関に特有の現象に対する配慮が不十分であるためと考えられる。
【０００５】
本発明は、上述のような課題を解決するためになされたもので、ポートインジェクタと筒内インジェクタとの双方を備えた内燃機関における空燃比を精度良く制御することのできる内燃機関の制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
第１の発明は、上記の目的を達成するため、内燃機関の制御装置であって、
内燃機関の吸気ポート内に燃料を噴射するポートインジェクタと、
前記内燃機関の筒内に燃料を噴射する筒内インジェクタと、
前記筒内から前記吸気ポートへの吸気の吹き返し中に含まれる燃料の量を所定のパラメータに基づいて推定する吹き返し燃料量推定手段と、
前記吹き返し燃料量推定手段により推定された量に基づいて、吹き返しによる誤差を補正するための吹き返し補正分を算入して、前記ポートインジェクタから噴射するポート噴射量と前記筒内インジェクタから噴射する筒内噴射量とを算出する噴射量算出手段と、
を備え、
前記所定のパラメータは、前記ポートインジェクタと前記筒内インジェクタとの噴射割合を含み、
前記吹き返し燃料量推定手段は、前記筒内インジェクタの噴射割合が大きい場合には小さい場合に比べて吹き返し燃料量を少なく推定することを特徴とする。
【０００７】
また、第２の発明は、上記の目的を達成するため、内燃機関の制御装置であって、
内燃機関の吸気ポート内に燃料を噴射するポートインジェクタと、
前記内燃機関の筒内に燃料を噴射する筒内インジェクタと、
前記筒内から前記吸気ポートへの吸気の吹き返し中に含まれる燃料の量を所定のパラメータに基づいて推定する吹き返し燃料量推定手段と、
前記吹き返し燃料量推定手段により推定された量に基づいて、吹き返しによる誤差を補正するための吹き返し補正分を算入して、前記ポートインジェクタから噴射するポート噴射量と前記筒内インジェクタから噴射する筒内噴射量とを算出する噴射量算出手段と、
を備え、
前記所定のパラメータは、前記筒内インジェクタの噴射開始時期を含み、
前記吹き返し燃料量推定手段は、前記筒内インジェクタの噴射開始時期が遅い場合には早い場合に比べて吹き返し燃料量を少なく推定することを特徴とする。
【０００８】
また、第３の発明は、第１または第２の発明において、
前記噴射量算出手段は、前記吹き返し補正分を前記ポート噴射量と前記筒内噴射量との双方に振り分けて算入する手段を含むことを特徴とする。
【０００９】
また、第４の発明は、第１または第２の発明において、
前記噴射量算出手段は、前記吹き返し補正分を前記筒内噴射量のみに算入する手段を含むことを特徴とする。
【００１０】
また、第５の発明は、第４の発明において、
前記吹き返し補正分が前記筒内噴射量のみに算入された場合に、前記筒内インジェクタの噴射開始時期を吸気弁が閉じた後の時期にする噴射時期制御手段を更に備えることを特徴とする。
【００１１】
また、第６の発明は、第４または第５の発明において、
前記噴射量算出手段は、前記吹き返し補正分を前記筒内噴射量のみに算入する際に、前サイクルの吹き返し補正分に対する今サイクルの吹き返し補正分の変化が判定値より大きい場合、前記吹き返し補正分の一部のみを前記筒内噴射量に算入する一部算入手段を含むことを特徴とする。
【００１２】
また、第７の発明は、第６の発明において、
前記吹き返し補正分が正の値であり、かつ、前記一部算入手段により前記吹き返し補正分の一部しか前記筒内噴射量に算入されなかった場合に、前記吹き返し補正分の残部に相当する量の燃料を排気行程において前記筒内インジェクタから噴射させる排気行程噴射手段を更に備えることを特徴とする。
【００１３】
また、第８の発明は、第６または第７の発明において、
前記吹き返し補正分が負の値であり、かつ、前記一部算入手段により前記吹き返し補正分の一部しか前記筒内噴射量に算入されなかった場合に、前記吹き返し補正分の残部を次サイクル以降の噴射量に算入することによって次サイクル以降の噴射量を減量させる噴射量減量手段を更に備えることを特徴とする。
【発明の効果】
【００１４】
第１の発明によれば、ポートインジェクタと筒内インジェクタとの双方を備えた内燃機関における燃料噴射量を算出する場合において、吸気の吹き返しの影響による誤差を補正することができる。また、その吹き返し補正分の基礎とする吹き返し燃料量を、ポートインジェクタと筒内インジェクタとの噴射割合の影響を反映させた上で、推定することができる。このため、様々な運転状態の下で、吹き返し燃料量を、ひいては吹き返し補正分を、精度良く求めることができる。このようなことから、本発明によれば、上記内燃機関における空燃比を精度良く制御することができる。
【００１５】
第２の発明によれば、ポートインジェクタと筒内インジェクタとの双方を備えた内燃機関における燃料噴射量を算出する場合において、吸気の吹き返しの影響による誤差を補正することができる。また、その吹き返し補正分の基礎とする吹き返し燃料量を、筒内インジェクタの噴射開始時期の影響を反映させた上で、推定することができる。このため、様々な運転状態の下で、吹き返し燃料量を、ひいては吹き返し補正分を、精度良く求めることができる。このようなことから、本発明によれば、上記内燃機関における空燃比を精度良く制御することができる。
【００１６】
第３の発明によれば、吹き返し補正分がポート噴射量と筒内噴射量との双方に算入されるので、ポート噴射量と筒内噴射量との比率が、吹き返し補正分の算入前と変わらないようにすることができる。このため、本発明によれば、ポートインジェクタと筒内インジェクタとの噴射割合を目標値に容易かつ精度良く制御することができる。
【００１７】
第４の発明によれば、吹き返し燃料量の算出時期を遅らせ、実際の吹き返し量が確定する吸気弁閉弁時に近づけることが可能となる。このため、本発明によれば、吹き返し燃料量を、ひいては吹き返し補正分を、更に精度良く求めることができる。その結果、空燃比を更に精度良く制御することができる。
【００１８】
第５の発明によれば、吹き返し燃料量の算出時期を、実際の吹き返し量が確定する吸気弁閉弁時の近傍まで遅らせることが可能となる。このため、本発明によれば、吹き返し燃料量を、ひいては吹き返し補正分を、特に高い精度で求めることができる。その結果、空燃比を特に高い精度で制御することができる。
【００１９】
第６の発明によれば、吹き返し補正分の算入に起因してトルクショックが生じることを確実に回避することができる。
【００２０】
第７の発明によれば、トルクショック回避のために吹き返し補正分が一部しか燃料噴射量に算入されなかった場合において、その残部が算入されなかったことによる排気空燃比のずれを回復させることができる。このため、本発明によれば、触媒酸素吸蔵量のバランスを良好に保つことが容易となる。
【００２１】
第８の発明によれば、トルクショック回避のために吹き返し補正分が一部しか燃料噴射量に算入されなかった場合において、その残部が算入されなかったことによる排気空燃比のずれを回復させることができる。このため、本発明によれば、触媒酸素吸蔵量のバランスを良好に保つことが容易となる。
【発明を実施するための最良の形態】
【００２２】
実施の形態１．
［システム構成の説明］
図１は、本発明の実施の形態１のシステム構成を説明するための図である。図１に示すように、本実施形態のシステムは、内燃機関６を備えている。内燃機関６は、複数の気筒を有する多気筒式の機関であり、図１は、そのうちの一気筒の断面を示している。内燃機関６の各気筒には、それぞれ、ピストン８と、燃焼室１０と、吸気弁１２と、排気弁１４と、点火プラグ１６と、筒内に連通する吸気ポート１８および排気ポート２０とが設けられている。吸気弁１２は、燃焼室１０と吸気ポート１８とを導通状態または遮断状態とするように開閉する。排気弁１４は、燃焼室１０と排気ポート２０とを導通状態または遮断状態とするように開閉する。
【００２３】
また、内燃機関６の各気筒に対しては、それぞれ、吸気ポート１８内に燃料を噴射するポートインジェクタ２２と、筒内（燃焼室１０内）に燃料を直接噴射する筒内インジェクタ２４とが設けられている。ポートインジェクタ２２および筒内インジェクタ２４には、それぞれ、図示しないポンプにより加圧された燃料が送り込まれている。内燃機関６では、各気筒に対して、ポートインジェクタ２２と筒内インジェクタ２４との双方から燃料を供給することができる。
【００２４】
吸気ポート１８は、吸気通路３０に連通している。吸気通路３０の上流端にはエアクリーナ３２が設けられ、空気はエアクリーナ３２を介して吸気通路３０内に取り込まれる。エアクリーナ３２の下流には、エアフロメータ３３が配置されている。エアフロメータ３３は、吸気通路３０内を流れる吸入空気量GAを検出するセンサである。吸気通路３０の下流部は気筒毎（吸気ポート１８毎）に分岐しており、その分岐部にはサージタンク３４が設けられている。吸気通路３０のサージタンク３４の上流にはスロットルバルブ３６が配置されている。スロットルバルブ３６には、その開度を検出するためのスロットルポジションセンサ３７が付設されている。スロットルバルブ３６より下流には、吸気管圧力PMを検出する吸気圧センサ３８が設けられている。
【００２５】
排気ポート２０には、燃焼室１０内での燃焼により生成された燃焼ガスを排出ガスとして排出するための排気通路４０が接続されている。排気通路４０には、排出ガスを浄化するための触媒４２が設けられている。排気通路４０における触媒４２の上流には、排出ガスの空燃比を検出するための空燃比センサ４４が配置されている。
【００２６】
内燃機関６のクランク軸４５の近傍には、クランク角センサ４６が設置されている。クランク角センサ４６は、クランク軸４５が所定回転角だけ回転する毎に、Hi出力とLo出力を反転させるセンサである。クランク角センサ４６の出力によれば、クランク軸４５の回転位置や回転速度、更には、機関回転数NEなどを検知することができる。また、内燃機関６は、冷却水温TWを検出する水温センサ４８を更に備えている。
【００２７】
内燃機関６の吸気弁１２は、可変動弁機構５０により駆動される。可変動弁機構５０は、クランク軸４５の回転と同期して吸気弁１２を開閉させると共に、そのバルブタイミングVTおよびバルブリフト量VLを変更することができる。可変動弁機構５０の近傍には、可変動弁機構５０の状態を検知するセンサ５２が設けられている。センサ５２の出力によれば、吸気弁１２のバルブタイミングVTおよびバルブリフト量VLの現実値を検出することができる。
【００２８】
本実施形態のシステムは、ECU(Electronic Control Unit)６０を備えている。ECU６０には、上述したエアフロメータ３３、空燃比センサ４４等の各種センサと、上述したポートインジェクタ２２、筒内インジェクタ２４、可変動弁機構５０等の各種アクチュエータが接続されている。ECU６０は、それらの各センサの出力に基づいて、各アクチュエータを適当に駆動することにより、内燃機関６の運転状態を制御することができる。
【００２９】
［燃料挙動モデル］
図２および図３は、本実施形態で用いる燃料挙動モデルを説明するための図である。この燃料挙動モデルは、ポートインジェクタ２２および筒内インジェクタ２４から噴射された後の燃料の挙動を表すモデルである。
【００３０】
ポートインジェクタ２２から噴射された燃料は、その一部が吸気ポート１８の内壁や吸気弁１２の外側に付着し、その残部が筒内に吸入されると考えられる。本モデルでは、図２に示すように、ポートインジェクタ２２から噴射される燃料の量をポート噴射量fipとし、吸気ポート１８の内壁に付着している燃料（以下、「ポート付着燃料」と称する）の総量をポート付着量fwpとし、吸気弁１２に付着している燃料（以下、「吸気弁付着燃料」と称する）の総量を吸気弁付着量fwvとする。なお、図２においては、ポート付着量fwp、吸気弁付着量fwvがそれぞれ一箇所にまとまって描かれているが、実際にはポート付着燃料や吸気弁付着燃料が一箇所にまとまって付着するとは限らない。
【００３１】
一方、筒内インジェクタ２４から噴射された燃料は、その一部がシリンダ内壁や燃焼室１０の内壁に付着し、その残部が気化して燃焼に寄与すると考えられる。本モデルでは、図２に示すように、筒内インジェクタ２４から噴射される燃料の量を筒内噴射量fidとし、シリンダ内壁や燃焼室１０の内壁に付着している燃料（以下、「筒内付着燃料」と称する）の総量を筒内付着量fwcとする。なお、図２においては、筒内付着量fwcが一箇所にまとまって描かれているが、実際には筒内付着燃料が一箇所にまとまって付着するとは限らない。
【００３２】
図３中のfcpは、吸気ポート１８から筒内へ吸入される燃料の量（以下、「筒内燃料吸入量」と称する）を表している。以下、図３を参照しつつ、本モデルにおける筒内燃料吸入量fcpの算出の仕方について説明する。
【００３３】
本モデルでは、ポートインジェクタ２２から噴射された燃料のうち、吸気ポート１８の内壁に付着する割合を「ポート付着率Rp」と定義し、吸気弁１２に付着する割合を「吸気弁付着率Rv」と定義する。この定義によれば、ポート噴射量fipのうち、筒内に吸入されることなくポート付着燃料に加わる量は、「Rp・fip」で表される。また、ポート噴射量fipのうち、筒内に吸入されることなく吸気弁付着燃料に加わる量は、「Rv・fip」で表される。一方、ポート噴射量fipのうち、噴射された後にそのまま筒内に吸入される量は「(1-Rp-Rv)・fip」で表される。
【００３４】
筒内燃料吸入量fcpには、上記の演算式(1-Rp-Rv)・fipで表される量のほかに、ポート付着燃料の気化分と、吸気弁付着燃料の気化分とが加算される。本モデルでは、今サイクル開始時（吸気弁１２の開弁前）にあったポート付着燃料のうち、今サイクルにおいて筒内に吸入されずに付着したまま残る割合を「ポート残留率Pp」と定義する。この定義によれば、「(1-Pp)・fwp」で表される量の燃料が今サイクルにおいてポート付着燃料の気化分として筒内に吸入されることになる。一方、ポート付着量fwpのうち、「Pp・fwp」で表される量は、今サイクル終了後にそのまま残る。同様にして、本モデルでは、今サイクル開始時にあった吸気弁付着燃料のうち、今サイクルにおいて筒内に吸入されずに付着したまま残る割合を「吸気弁残留率Pv」と定義する。この定義によれば、「(1-Pv)・fwv」で表される量の燃料が今サイクルにおいて吸気弁付着燃料の気化分として筒内に吸入されることになる。一方、吸気弁付着量fwvのうち、「Pv・fwv」で表される量は、今サイクル終了後にそのまま残る。
【００３５】
上述したポート付着率Rp、吸気弁付着率Rv、ポート残留率Pp、吸気弁残留率Pvや、後述する筒内付着率Rc、筒内残留率Pcは、それぞれ、内燃機関６の運転状態に応じて変化する特性を有している。また、本モデルによれば、上述したポート付着量fwp、吸気弁付着量fwv、筒内付着量fwcは、それぞれ、新たな付着による増加と気化による減少との収支により、内燃機関６が１サイクル動作する毎に変化していく。そこで、以下では、第kサイクルにおけるそれらの値を表す際には、それらの符号の後に(k)を付すこととする。なお、ポート付着量fwp、吸気弁付着量fwvおよび筒内付着量fwcの初期値fwp(0)、fwc(0)およびfwc(0)は、いずれもゼロである。
【００３６】
以上より、第kサイクルにおける筒内燃料吸入量fcp(k)は、次式（１）のように表すことができる。また、第kサイクル終了後、つまり第k+1サイクル開始時における、ポート付着量fwp(k+1)および吸気弁付着量fwv(k+1)は、それぞれ、次式（２）および（３）のように表すことができる。
【００３７】
【数１】
0004470832.tif 000002
【００３８】
一方、図３中のfcdは、筒内で気化して燃焼に寄与する燃料の量（以下、「筒内燃料気化量」と称する）を表している。以下、図３を参照しつつ、本モデルにおける筒内燃料気化量fcdの算出の仕方について説明する。
【００３９】
本モデルでは、筒内インジェクタ２４から噴射された燃料のうち、シリンダ内壁や燃焼室１０の内壁に付着する割合を「筒内付着率Rc」と定義する。この定義によれば、筒内噴射量fidのうち、気化することなく筒内付着燃料に加わる量は、「Rc・fid」で表される。一方、筒内噴射量fidのうち、噴射後そのまま気化する量は「(1-Rc)・fid」で表される。
【００４０】
筒内燃料気化量fcdには、上記の演算式(1-Rc)・fidで表される量のほかに、筒内付着燃料の気化分が加算される。本モデルでは、今サイクル開始時にあった筒内付着燃料のうち、今サイクルにおいて気化することなく付着したまま残る割合を「筒内残留率Pc」と定義する。この定義によれば、「(1-Pc)・fwc」で表される量の燃料が今サイクルにおける筒内付着燃料の気化分に相当する。一方、筒内付着量fwcのうち、「Pc・fwc」で表される量は、今サイクル終了後にそのまま残る。
【００４１】
以上より、第kサイクルにおける筒内燃料気化量fcd(k)は、次式（４）のように表すことができる。また、第kサイクル終了後、つまり第k+1サイクル開始時における、筒内付着量fwc(k+1)は、次式（５）のように表すことができる。
【００４２】
【数２】
0004470832.tif 000003
【００４３】
以下、第kサイクルにおいて、筒内で燃焼に寄与し得る燃料の量（以下、「燃焼寄与燃料量」と称する）をfc(k)と表す。第kサイクルにおける燃焼寄与燃料量fc(k)は、以上説明した燃料挙動モデルによれば、筒内燃料吸入量fcp(k)と、筒内燃料気化量fcd(k)との和と考えられるので、次式（６）のように表すことができる。
【００４４】
【数３】
0004470832.tif 000004
【００４５】
［吹き返し燃料量fbの推定］
上述した燃料挙動モデルによれば、噴射された燃料の壁面への付着分や、壁面に付着している燃料の気化分を考慮することができる。このため、所望の空燃比を実現する上で、燃料噴射量を精度良く算出することができる。本実施形態では、更なる精度向上を図るため、このような燃料挙動モデルに、吸気の吹き返し分の補正を加えることとする。
【００４６】
ここで、吸気の吹き返しとは、筒内に一旦吸入された新気が、吸気弁１２が閉じる前に、吸気ポート１８へ逆流することである。一般に、内燃機関においては、吸気慣性を有効に利用して吸入効率を高める等の理由から、吸気弁１２は、下死点を過ぎた後に閉じることとされる。このため、吸気下死点を過ぎてピストン８が上昇し始め、かつ吸気弁１２がまだ閉じていない状態では、筒内から吸気ポート１８への新気の逆流が起こり得る。
【００４７】
このような吸気の吹き返しに伴い、その中に含まれる燃料も筒内から吸気ポート１８へ戻されることとなる。以下では、吸気の吹き返し中に含まれる燃料を「吹き返し燃料」と称する。また、吹き返し燃料の量を「吹き返し燃料量」と称し、符号fbで表すこととする。本実施形態では、吹き返し燃料量fbを推定し、その推定値に基づいて、吸気の吹き返しの影響による誤差を補正する。本実施形態のような、ポートインジェクタ２２および筒内インジェクタ２４の双方を備えた内燃機関６においては、吹き返し燃料には、ポートインジェクタ２２から噴射された燃料と、筒内インジェクタ２４から噴射された燃料とのいずれもが含まれ得る。
【００４８】
吹き返し燃料量fbは、内燃機関６の運転状態に応じて変化する。より具体的には、吹き返し燃料量fbは、内燃機関６の機関回転数NEや、負荷率KL（あるいは吸気管圧力PM）、吸気弁１２のバルブタイミングVTおよびバルブリフト量VL、更には、ポート噴射量fipと筒内噴射量fidとを合計した全噴射量fiといったパラメータに応じて変化する。例えば、全噴射量fiが多いほど、吹き返し燃料量fbは多くなる。また、吸気弁１２のバルブリフト量VLが大きいほど、吹き返し燃料量fbは多くなる。したがって、吹き返し燃料量fbを様々な運転状態において常に精度良く推定するためには、それらのパラメータの値に応じた吹き返し燃料量fbを算出する必要がある。
【００４９】
そこで、本実施形態では、機関回転数NE、負荷率KL、バルブタイミングVT、バルブリフト量VLおよび全噴射量fiの各々に対する吹き返し燃料量fbの依存性を予め把握しておき、それらの依存性を表す情報をECU６０に記憶させておく。そして、ECU６０は、それらの依存性を反映させて、吹き返し燃料量fbを算出する。これにより、様々な運転状態下において、吹き返し燃料量fbを精度良く推定することが可能となる。
【００５０】
ところで、本実施形態のシステムでは、筒内インジェクタ２４の燃料噴射開始時期（以下、「筒内噴射開始時期」と称する）TDは、固定されてはおらず、内燃機関６の運転状態に応じて最適となるように調整される。図４は、吹き返し燃料量fbと筒内噴射開始時期TDとの関係を、他のパラメータを一定にした条件下で表したグラフである。図４に示されるように、本実施形態の内燃機関６においては、吹き返し燃料量fbは、筒内噴射開始時期TDにも依存する。その訳は、次の通りである。
【００５１】
上述したように、吹き返し燃料には、ポートインジェクタ２２から噴射された分と、筒内インジェクタ２４から噴射された分とが含まれ得る。一方、筒内インジェクタ２４の噴射期間が吸気弁閉弁時をまたぐ場合、吸気弁閉弁時以後に噴射される燃料については、吹き返しによって吸気ポート１８に戻ることはない。このため、図４に示されるように、筒内インジェクタ２４の噴射期間のうち吸気弁開弁期間と重なる期間が短くなると、これに比例するようにして、筒内噴射量fidのうち吸気ポート１８へ吹き返す分が少なくなる。そして、筒内噴射開始時期TDが吸気弁閉弁時以後になると、筒内噴射量fidのうち吸気ポート１８へ吹き返す分はゼロになり、筒内噴射量fidの全部が筒内に残る。このため、図４中、筒内噴射開始時期TDが吸気弁閉弁時以後になる範囲においては、吹き返し燃料量fbは、ポートインジェクタ２２から噴射された分だけになるので、一定の値をとる。
【００５２】
このようにして、吹き返し燃料量fbは、他のパラメータを一定にした場合、筒内噴射開始時期TDが遅いほど少なくなる。本実施形態では、吹き返し燃料量fbの推定に際し、このような傾向を反映させるべく、吹き返し燃料量fbの推定の基礎とするパラメータに、筒内噴射開始時期TDを加えることとした。これにより、本実施形態では、吹き返し燃料量fbを更に精度良く推定することができる。
【００５３】
また、本実施形態では、吹き返し燃料量fbは、ポートインジェクタ２２と筒内インジェクタ２４との噴射割合にも依存する。本実施形態では、この噴射割合を表すのに、噴き分け率γを用いることとする。噴き分け率γは、ポート噴射量fipと筒内噴射量fidとの合計を１としたときのポート噴射量fipの比率である。すなわち、ポート噴射量fip：筒内噴射量fid＝γ:(1-γ)となる。本実施形態のシステムでは、この噴き分け率γは、固定されてはおらず、内燃機関６の運転状態に応じて最適となるように調整される。
【００５４】
図５は、吹き返し燃料量fbと噴き分け率γとの関係を、他のパラメータを一定にした条件下で表したグラフである。上述したように、筒内インジェクタ２４から噴射される燃料のうち、吸気弁閉弁時以後に噴射される燃料については、吹き返しによって吸気ポート１８に戻ることはない。つまり、筒内インジェクタ２４から噴射される燃料は、ポートインジェクタ２２から噴射される燃料に比べて、吹き返しにくいと言える。このため、図５に示されるように、筒内噴射の割合が高いほど、換言すれば噴き分け率γが０に近いほど、吹き返し燃料量fbは少なくなる。
【００５５】
そこで、本実施形態では、吹き返し燃料量fbの推定に際し、このような傾向を反映させるべく、吹き返し燃料量fbの推定の基礎とするパラメータに、噴き分け率γを加えることとした。これにより、本実施形態では、吹き返し燃料量fbを更に精度良く推定することができる。
【００５６】
［ポート噴射量fipおよび筒内噴射量fidの算出］
次に、所望の空燃比を実現するために必要なポート噴射量fipおよび筒内噴射量fidを、上述したモデルに基づいて算出する方法について説明する。
【００５７】
本実施形態では、ポート噴射量fipおよび筒内噴射量fidを更に良い精度で算出するために、筒内残留ガスの影響を考慮する。内燃機関においては、一般に、既燃ガスを筒内から完全に排出することは困難であるため、既燃ガスの一部は次回サイクルまで筒内に残留する。以下では、第ｋサイクル終了後、つまり第k+1サイクル開始時に筒内に残留するガス中に含まれる燃料（既燃、未燃を問わない）の量を「筒内残留燃料量fr(k)」と表す。
【００５８】
以下、第ｋサイクルにおける吹き返し燃料量をfb(k)と表す。また、第kサイクルの動作によって排気ポート２０へ送り出される排出ガス中に含まれる燃料（既燃、未燃を問わない）の量を「排出燃料量fex(k)」と表す。このfex(k)は、上述した燃焼寄与燃料量fc(k)に、吹き返しの影響に対する補正と、筒内残留ガスの影響に対する補正とを加えることにより、正確に算出することができる。すなわち、第ｋサイクルにおける排出燃料量fex(k)は、次式（７）のように表すことができる。
【００５９】
【数４】
0004470832.tif 000005
【００６０】
上記（７）式の右辺における-fb(k)なる項は、今サイクルにおいて吹き返した分が減算されることを意味する。また、+fb(k-1)なる項は、前サイクルにおいて吹き返した分が今サイクルにおいて吸入されることによって加算されることを意味する。また、-fr(k)なる項は、今サイクルにおいて筒内に残留して排出されなかった分が減算されることを意味する。また、+fr(k-1)なる項は、前サイクルから筒内に残留していた分が今サイクルにおいて排出されることによって加算されることを意味する。
また、上記（７）式の右辺の２〜３行目は、右辺１行目のfc(k)に上記（６）式を代入したものである。
【００６１】
以下、第ｋサイクルにおける噴き分け率をγ(k)と表す。これにより、噴き分け率γの定義から、次式（８）が成り立つ。
また、第ｋサイクルにおける目標空燃比をα(k)とし、第ｋサイクルにおいて筒内に吸入された空気量（以下、「筒内吸入空気量」と称する）をｍ(k)とすると、目標空燃比α(k)を実現するために要求される条件は、次式（９）で表すことができる。
【００６２】
【数５】
0004470832.tif 000006
【００６３】
上記（７）式に上記（８）式および（９）式を代入して、fip(k)とfid(k)についてそれぞれ解くことにより、次式（１０）および（１１）が得られる。
【００６４】
【数６】
0004470832.tif 000007
【００６５】
本実施形態では、上記（１０）式および（１１）式に基づいて、ポート噴射量fipおよび筒内噴射量fidを算出する。本実施形態では、上記（１０）式および（１１）式中のfb(k)、fb(k-1)を含む項で表される量が、吹き返しによる誤差を補正するための吹き返し補正分に相当する。本実施形態では、この吹き返し補正分を算入してポート噴射量fipおよび筒内噴射量fidを算出することにより、吹き返しの影響による空燃比の誤差（ずれ）を正確に補正することができるので、目標空燃比α(k)を精度良く実現することができる。
【００６６】
［実施の形態１における具体的処理］
図６は、上述した機能を実現するために本実施形態においてECU６０が実行するルーチンのフローチャートを示す。なお、このルーチンは、内燃機関６の気筒毎に別々に実行される。また、このルーチンは、対象とする気筒が１サイクル動作する毎に実行される。
【００６７】
図６に示すルーチンによれば、まず、現在の内燃機関６の状態を表す各種のパラメータが取得される（ステップ１００）。具体的には、ここでは、機関回転数NE、負荷率KL、バルブタイミングVT、バルブリフト量VL、全噴射量fi、筒内吸入空気量ｍ、吸気管圧力PM、冷却水温TWなど、内燃機関６の運転状態を表す物理量が、計測により、あるいは公知の演算により求められる。ここで、負荷率KLは、スロットルバルブ３６を全開にしたときに得られる吸入空気量と、実吸入空気量GAとの比率であり、エアフロメータ３３の出力等に基づいて算出することができる。また、筒内吸入空気量ｍは、エアフロメータ３３の出力、あるいは、吸気圧センサ３８の出力に基づいて算出することができる。
【００６８】
なお、ECU６０は、別のルーチンにおいて、内燃機関６の運転状態等に基づいて、目標空燃比αを逐次設定している。上記ステップ１００の処理において、全噴射量fiは、この目標空燃比αと、筒内吸入空気量ｍとに基づいて算出される。このようにして算出された全噴射量fiは、実際に噴射される量ではなく、図６に示すルーチン中の後述するステップ１０８の処理において、吹き返し燃料量fbを算出する基礎として近似的に用いられるものである。
【００６９】
図６に示すルーチンでは、次に、ポート付着率Rp、吸気弁付着率Rv、筒内付着率Rc、ポート残留率Pp、吸気弁残留率Pvおよび筒内残留率Pcをそれぞれ算出する処理が実行される（ステップ１０２）。ECU６０は、それらの値を、内燃機関６の状態を表す各種のパラメータ（例えば、機関回転数NE、負荷率KL、バルブタイミングVT、冷却水温TW等）との関係で定めたマップを記憶している。ここでは、そのマップを参照して、現在の運転状態に応じた値が算出される。
【００７０】
次に、目標とする筒内噴射開始時期TDを算出する処理が行われる（ステップ１０４）。ECU６０は、最適な筒内噴射開始時期TDを、内燃機関６の状態を表す各種のパラメータとの関係において定めたマップを記憶している。ここでは、そのマップを参照することにより、現在の運転状態に適した筒内噴射開始時期TDが目標として設定される。
【００７１】
次に、目標とする噴き分け率γを算出する処理が行われる（ステップ１０６）。ECU６０は、最適な噴き分け率γを、内燃機関６の状態を表す各種のパラメータとの関係において定めたマップを記憶している。ここでは、そのマップを参照することにより、現在の運転状態に適した噴き分け率γが目標として設定される。
【００７２】
次に、吹き返し燃料量fbを算出する処理が次のようにして行われる（ステップ１０８）。ECU６０は、所定の基準吹き返し燃料量を記憶している。更に、ECU６０は、吹き返し燃料量fbに影響を及ぼす複数のパラメータである機関回転数NE、負荷率KL、バルブタイミングVT、バルブリフト量VL、全噴射量fi、筒内噴射開始時期TDおよび噴き分け率γの各々について、それらの影響を反映させるための補正係数を定めたマップを記憶している。そして、ECU６０は、それらのマップを参照して各パラメータ毎の補正係数を取得し、それらの補正係数を上記基準吹き返し燃料量にすべて掛け合わせることにより、吹き返し燃料量fbを算出する。
【００７３】
上記ステップ１０８の処理において、筒内噴射開始時期TDについての補正係数を定めたマップとしては、図４に示すような傾向を吹き返し燃料量fbに反映させることが可能なマップが用いられる。これにより、筒内噴射開始時期TDの違いが及ぼす影響を正確に織り込むことができるので、現実の値により合致した吹き返し燃料量fbを求めることができる。
【００７４】
また、上記ステップ１０８の処理において、噴き分け率γについての補正係数を定めたマップとしては、図５に示すような傾向を吹き返し燃料量fbに反映させることが可能なマップが用いられる。これにより、噴き分け率γの違いが及ぼす影響を正確に織り込むことができるので、現実の値により合致した吹き返し燃料量fbを求めることができる。
【００７５】
図６に示すルーチンでは、次に、筒内残留燃料量frを算出する処理が行われる（ステップ１１０）。ECU６０は、筒内残留燃料量frの値を、内燃機関６の状態を表す各種のパラメータとの関係で定めたマップを記憶している。ここでは、そのマップを参照して、現在の運転状態に応じた値が算出される。
【００７６】
次に、ポート噴射量fipおよび筒内噴射量fidを次のようにして算出する処理が行われる（ステップ１１２）。ECU６０は、上記（１０）式および（１１）式を記憶している。ECU６０は、それらの式に、今回の処理サイクルの上記ステップ１０２〜１１０で算出された値と、前回の処理サイクル時の上記ステップ１１０で算出された値とを代入する。また、ECU６０は、後述するステップ１１８で算出される吹き返し燃料量fbの前回の処理サイクル時の値を、上記（１０）式および（１１）式中の「fb(k-1)」として代入する。更に、ECU６０は、後述するステップ１１４で算出されるポート付着量fwp、吸気弁付着量fwvおよび筒内付着量fwcの前回の処理サイクル時の値を上記（１０）式および（１１）式に代入する。このような処理により、ポート噴射量fipおよび筒内噴射量fidを算出することができる。
【００７７】
なお、ECU６０は、別のルーチンにおいて、現在の運転状態に適したポート噴射時期を設定している。そのポート噴射時期に合わせて、上記ステップ１１２で算出されたポート噴射量fipの分の燃料がポートインジェクタ２２から噴射される。また、上記ステップ１０４で算出された筒内噴射開始時期TDに合わせて、上記ステップ１１２で算出された筒内噴射量fidの分の燃料が筒内インジェクタ２４から噴射される。
【００７８】
図６に示すルーチンでは、次に、ポート付着量fwp、吸気弁付着量fwvおよび筒内付着量fwcを算出する処理が行われる（ステップ１１４）。ECU６０は、上記（２）式、（３）式および（５）式を記憶している。ECU６０は、それらの式に、今回の処理サイクルの上記ステップ１０２および１０８で算出された値と共に、前回の処理サイクル時に求めたポート付着量fwp、吸気弁付着量fwvおよび筒内付着量fwcを代入することにより、ポート付着量fwp、吸気弁付着量fwvおよび筒内付着量fwcを最新値に更新することができる。このようにして更新されたポート付着量fwp、吸気弁付着量fwvおよび筒内付着量fwcの値は、次回の処理サイクルでの吹き返し燃料量fbの算出（ステップ１０８）に用いられる。
【００７９】
図６に示すルーチンでは、次に、対象とする気筒の吸気弁１２が閉じたか否かを判断する（ステップ１１６）。この判断は、クランク角センサ４６によって検出されるクランク軸４５の回転位置や、センサ５２によって検出されるバルブタイミングVTなどに基づいて行われる。
【００８０】
上記ステップ１１６において吸気弁１２が閉じたと判断された後、吹き返し燃料量fbを最新の情報に基づいて再度算出する処理が行われる（ステップ１１８）。ここで再計算された吹き返し燃料量fbは、既述した通り、次回の処理サイクルでのポート噴射量fipおよび筒内噴射量fidの算出時（ステップ１１２）に、fb(k-1)として用いられる。
【００８１】
上記ステップ１１８において吹き返し燃料量fbを再計算する理由は、次の通りである。吸気の吹き返しは吸気弁１２が開弁しているときに生じるので、吸気弁１２が閉じた瞬間に実際の吹き返し量が確定する。このため、吸気弁閉弁時に算出した吹き返し燃料量fbは、現実の値に更に近いものとなる。よって、吸気弁閉弁時に再計算した吹き返し燃料量fbを、次回の処理サイクルでfb(k-1)として使用することにより、ポート噴射量fipおよび筒内噴射量fidを更に精度良く算出することができる。
【００８２】
以上説明した通り、図６に示すルーチンによれば、吹き返し燃料量fbに、筒内噴射開始時期TDの影響および噴き分け率γの影響を反映させることにより、吹き返し燃料量fbを精度良く推定することができる。そして、その吹き返し燃料量fbに基づく吹き返し補正分をポート噴射量fipおよび筒内噴射量fidの計算に算入することにより、様々な運転状態の下で、実際の空燃比を目標空燃比αに精度良く制御することができる。
【００８３】
また、本実施形態では、上記（１０）式および（１１）式を見て分かる通り、吹き返し補正分がポート噴射量fipと筒内噴射量fidとの双方に振り分けて算入される。その結果、目標とされた噴き分け率γを満足するようにポート噴射量fipおよび筒内噴射量fidを算出することができる。このため、実際の噴き分け率γを容易かつ精度良く制御することができる。
【００８４】
ところで、上述した実施の形態１では、吸気弁１２が閉じるのを待って（上記ステップ１１６）、吹き返し燃料量fbを再計算する（上記ステップ１１８）ようにしているが、本発明では、これらの処理は行わなくても良い。すなわち、上記ステップ１０８で算出された吹き返し燃料量fbを、次回の処理サイクルのステップ１１２において「fb(k-1)」として用いるようにしてもよい。
【００８５】
また、上述した実施の形態１では、吹き返し燃料量fbを推定する基礎とするパラメータに、筒内噴射開始時期TDと噴き分け率γとの双方を含ませることとしているが、本発明では、筒内噴射開始時期TDと噴き分け率γとの何れか一方のみを含ませることとしてもよい。その場合においても、吹き返し燃料量fbを十分な精度で推定することができる。
【００８６】
また、吹き返し燃料量fbを推定する基礎とするパラメータのうち筒内噴射開始時期TDおよび噴き分け率γ以外のものについては、上述した実施の形態１で述べたものが必ずしもすべて含まれていなくても良い。例えば、吸気弁１２の可変動弁機構５０を備えないシステムの場合には、バルブタイミングVTおよびバルブリフト量VLは、吹き返し燃料量fbを推定する基礎とするパラメータに含める必要はない。逆に、吹き返し燃料量fbを推定する基礎とするパラメータとして、上述した実施の形態１で述べていないものが更に含まれていても良い。
【００８７】
また、上述した実施の形態１においては、吹き返し燃料量fbの算出を、複数のマップを参照して得られる複数の補正係数を基準となる量に乗算する方法によって行っているが、吹き返し燃料量fbの算出方法は、このような方法に限定されるものではない。例えば、吸気ポート１８内の流量を算出することが可能な公知のシリンダモデルを使用する方法を採ってもよい。その場合、シリンダモデルにより算出される吸気ポート１８内の流量を、吹き返しが生じている期間中に渡って積算することにより、総吹き返し量を算出することができる。その総吹き返し量に基づいて、吹き返し燃料量fbを算出することができる。
【００８８】
また、上述した実施の形態１においては、機関回転数NE、負荷率KL、バルブタイミングVT、バルブリフト量VL、全噴射量fi、筒内噴射開始時期TDおよび噴き分け率γが前記第１および第２の発明における「所定のパラメータ」に相当している。また、ECU６０が、上記ステップ１０８の処理を実行することにより前記第１および第２の発明における「吹き返し燃料量推定手段」が、上記ステップ１１２の処理を実行することにより、前記第１および第２の発明における「噴射量算出手段」が、それぞれ実現されている。
【００８９】
また、上述した実施の形態１においては、ECU６０が、上記ステップ１１２の処理を実行することにより前記第３の発明における「振り分けて算入する手段」が実現されている。
【００９０】
実施の形態２．
次に、図７を参照して、本発明の実施の形態２について説明するが、上述した実施の形態１との相違点を中心に説明し、同様の事項については、その説明を省略または簡略する。本実施形態のシステムは、図１に示すハードウェア構成を用いて、ECU６０に、後述する図７に示すルーチンを実行させることにより実現することができる。
【００９１】
［実施の形態２の特徴］
既述した通り、実施の形態１においては、吹き返し補正分を、ポート噴射量fipと筒内噴射量fidとの双方に振り分けて算入する。これに対し、本実施形態においては、吹き返し補正分を筒内噴射量fidのみに一括して算入する。以下、本実施形態におけるポート噴射量fipおよび筒内噴射量fidの算出方法について、具体的に説明する。
【００９２】
［ポート噴射量fipの算出］
本実施形態では、既述した筒内燃料吸入量fcp(k)と、目標空燃比α(k)と、筒内吸入空気量ｍ(k)と、噴き分け率γ(k)との間に次式の関係が成り立つように、ポート噴射量fipを設定する。
fcp(k)＝γ(k)・ｍ(k)／α(k)              ・・・（１２）
上記（１）式および（１２）式から、ポート噴射量fipは、次式のように算出することができる。
【００９３】
【数７】
0004470832.tif 000008
【００９４】
［筒内噴射量fidの算出］
排出燃料量fex(k)は、既述した通り、吹き返しの影響と、筒内残留ガスの影響とを考慮することにより、次式のように表すことができる。
fex(k)＝fcp(k)+fcd(k)+fb(k-1)+fr(k-1)-fb(k)-fr(k)      ・・・（１４）
上記（４）式、（９）式、（１２）式および（１４）式から、筒内噴射量fidは、次式のように算出することができる。
【００９５】
【数８】
0004470832.tif 000009
【００９６】
［吹き返し補正分fbc］
上記（１５）式の右辺の２段目は、右辺１段目のfb(k)、fb(k-1)を含む項をまとめてfbcとして置き換えたものである。本実施形態では、このfbcが、吹き返しによる誤差を補正するための吹き返し補正分に相当する。すなわち、吹き返し補正分fbcは、次式で表される。
fbc(k)＝{fb(k-1)-fb(k)}／{1-Rc(k)}          ・・・（１６）
本実施形態では、ECU６０は、上記（１３）式、（１５）式および（１６）式を記憶している。
【００９７】
本実施形態では、上記吹き返し補正分fbcを筒内噴射量fidのみに算入することにより、次のような利点がある。吹き返し補正分fbcの基礎となるのは、吹き返し燃料量fbである。そして、実際の吹き返し量が確定するのは、既述した通り、吸気弁閉弁時である。このため、最も現実に合致した吹き返し燃料量fbを算出することができる時期は、吸気弁閉弁時の近傍である。一方、吸気弁１２が閉じた後では、ポートインジェクタ２２から噴射した燃料を筒内に供給することはできないので、ポート噴射時期は吸気弁閉弁時より前とされる。このため、吹き返し補正分fbcをポート噴射量fipに算入する場合には、吹き返し燃料量fbは、吸気弁閉弁時よりも前に算出しておかなければならない。
【００９８】
これに対し、筒内噴射の場合には、吸気弁１２が閉じた後も、つまり圧縮行程においても、筒内に燃料を供給することができる。したがって、吹き返し補正分fbcを筒内噴射量fidのみに算入する場合には、吸気弁閉弁時まで待ってから吹き返し燃料量fbを算出した場合でも、吸気弁閉弁後の圧縮行程時に、吹き返し補正分fbcを加味した筒内噴射を行うことが可能である。すなわち、本実施形態では、吸気弁閉弁時まで待ってから最新の情報に基づいて吹き返し燃料量fbを求めることができるので、吹き返し燃料量fbを最良の精度で算出することができる。このため、吹き返しの影響による空燃比のずれを更に精度良く補正することができる。
【００９９】
本実施形態では、上記の理由から、筒内噴射開始時期TDを吸気弁閉弁後の圧縮行程中とする。図４を参照して既述した通り、筒内噴射開始時期TDが吸気弁閉弁時以後になる範囲においては、吹き返し燃料量fbは、筒内噴射開始時期TDによらず一定の値をとる。従って、本実施形態では、吹き返し燃料量fbの推定の基礎とするパラメータに、筒内噴射開始時期TDを加える必要がない。このため、吹き返し燃料量fbを筒内噴射開始時期TDによらずに推定することとした。
【０１００】
ところで、吹き返し補正分fbcは、上記（１６）式を見て分かる通り、前サイクルに生じた吹き返し燃料量fb(k-1)と今サイクルに生じた吹き返し燃料量fb(k)との収支に基づくものである。このため、吹き返し補正分fbcは、正の値にも負の値にもなり得る。つまり、吹き返し補正分fbcは、増量補正となる場合と、減量補正となる場合との両方がある。
【０１０１】
定常運転時には、一般に、この吹き返し補正分fbcのサイクル毎の変動は小さいものとなる。これに対し、過渡運転時には、吹き返し補正分fbcのサイクル毎の変動が大きくなり易い。吹き返し補正分fbcが前サイクルでの値に対して大きく変化した場合には、今サイクルの燃料噴射量が前サイクルに比べて急増または急減し易いので、今サイクルと前サイクルとの間のトルク段差が生じ易い。このトルク段差が大きくなり過ぎると、トルクショックが生じ易くなる。
【０１０２】
そこで、本実施形態では、吹き返し補正に起因してトルクショックが生じることを確実に回避するため、吹き返し補正分fbcが前サイクルでの値に対して大きく変化した場合には、吹き返し補正分fbcの一部のみを筒内噴射量fidに算入することとした。これにより、今サイクルと前サイクルとの間のトルク段差に起因するトルクショックが生じるのを確実に回避することができる。
【０１０３】
ところで、触媒４２に良好な排気浄化作用を発揮させるためには、触媒４２の酸素吸蔵量（以下、「触媒酸素吸蔵量」と称する）を適度な範囲に保つ必要がある。このため、本実施形態のシステムにでは、触媒酸素吸蔵量を適度な範囲に制御するべく、空燃比を調整することが行われる。このようなシステムにおいて、上述したトルクショック防止のために吹き返し補正分fbcの一部しか筒内噴射量fidに算入しなかった場合には、目標空燃比αからのずれが大きくなり易いので、触媒酸素吸蔵量を制御し難くなり易い。
【０１０４】
例えば、吹き返し補正分fbcが正の値であり、かつ、吹き返し補正分fbcの一部しか筒内噴射量fidに算入されなかった場合には、燃料の総噴射量が目標空燃比αを達成するには足りないこととなるので、触媒酸素吸蔵量が過多になり易い。そこで、本実施形態では、このような場合であって、実際に触媒酸素吸蔵量に過多になっている場合には、吹き返し補正分fbcのうち筒内噴射量fidに算入されなかった分に相当する量の燃料を排気行程において筒内インジェクタ２４から噴射させることとした。この排気行程において噴射された燃料は、燃焼に寄与しないので、トルクを増大させることがない。よって、トルクの急増によるトルクショックを発生させることなく、触媒酸素吸蔵量のバランスを保つことができる。
【０１０５】
一方、吹き返し補正分fbcが負の値であり、かつ、吹き返し補正分fbcの一部しか筒内噴射量fidに算入されなかった場合には、燃料の総噴射量が目標空燃比αを達成する量より多くなることとなるので、触媒酸素吸蔵量が不足し易い。そこで、本実施形態では、このような場合であって、実際に触媒酸素吸蔵量に不足している場合には、吹き返し補正分fbcのうち筒内噴射量fidに算入されなかった分を次サイクルのポート噴射量fipに算入することにより、次サイクルのポート噴射量fipを減量させることとした。すなわち、吹き返し補正分fbcを今サイクルと次サイクルとに分散させることで、燃料噴射量の減少によるトルク段差を分散させて小さくすることができる。このため、トルクショックを発生させることなく、触媒酸素吸蔵量のバランスを保つことができる。
【０１０６】
［実施の形態２における具体的処理］
図７は、上記の機能を実現するために本実施形態においてECU６０が実行するルーチンのフローチャートである。なお、このルーチンは、内燃機関６の気筒毎に別々に実行される。また、このルーチンは、対象とする気筒が１サイクル動作する毎に実行される。
【０１０７】
図７に示すルーチンによれば、まず、フラグFに１がセットされているか否かが判別される（ステップ１３０）。このフラグFの意味するところについては後述する。フラグFの初期値は０とされている。このため、本ルーチンが初めて実行された際には、上記ステップ１３０においてF＝１の不成立が認められる。この場合には、次に、ポート噴射量fipが算出される（ステップ１３２）。このポート噴射量fipの算出は、次のようにして行われる。
【０１０８】
ECU６０は、既述した図６に示すルーチンと同様にして、付着量fwp、fwv、fwcをサイクル毎に最新値に更新している。また、上記ステップ１３２においては、図６に示すルーチン中のステップ１０２、１０６および１１０と同様の処理が行われ、付着率Rp、Rv、Rc、残留率Pp、Pv、Pc、噴き分け率γ、筒内残留燃料量frとが算出される。ECU６０は、これらのパラメータを上記（１３）式に代入することにより、ポート噴射量fipを算出する。
【０１０９】
図７に示すルーチンでは、次に、算出されたポート噴射量fipの分の燃料がポートインジェクタ２２から噴射される（ステップ１３４）。このポート噴射の時期は、吸気弁１２が開弁する前か、または吸気弁１２の開弁期間中とされる。
【０１１０】
次に、吹き返し補正分fbcの変化率Ｓが算出される（ステップ１３６）。この変化率Ｓは、前サイクルの吹き返し補正分fbc(k-1)に対する今サイクルの吹き返し補正分fbc(k)の変化が、トルクショックの原因になるおそれがあるほどに大きいか否かを判断する基礎とされるものである。具体的には、変化率Ｓは、次式に示す通り、今サイクルの吹き返し補正分fbc(k)と前サイクルの吹き返し補正分との差を全噴射量fi(k)で除した値として定義される。
Ｓ＝{fbc(k)-fbc(k-1)}／fi(k)        ・・・（１７）
【０１１１】
上記ステップ１３６の変化率Ｓの算出処理では、まず、今サイクルの吹き返し燃料量fb(k)が算出される。この吹き返し燃料量fb(k)の算出処理は、筒内噴射開始時期TDを算出の基礎としないこと以外は、既述した図６に示すルーチンのステップ１０８と同様にして行われる。ここでは、吹き返し燃料量fb(k)は、推定精度が最も高くなるタイミングである吸気弁閉弁時の近傍で、機関回転数NE、負荷率KL、バルブタイミングVT、バルブリフト量VL、全噴射量fiおよび噴き分け率γの最新の値を取得した上で、それらの値に基づいて算出される。これにより、本実施形態では、吹き返し燃料量fbを特に高い精度で算出することができる。
【０１１２】
上記ステップ１３６の処理においては、次に、ECU６０は、算出された今サイクルのfb(k)と、前回の処理サイクルで算出されたfb(k-1)とを上記（１６）式に代入することにより、今サイクルの吹き返し補正分fbc(k)を算出する。そして、ECU６０は、算出された今サイクルの吹き返し補正分fbc(k)と、前回の処理サイクルで算出されたfbc(k-1)とを上記（１７）式に代入することにより、変化率Ｓを算出する。
【０１１３】
図７に示すルーチンでは、次に、吹き返し補正分fbcの変化率Ｓの絶対値が所定の判定値より大きいか否かが判別される（ステップ１３８）。その結果、変化率Ｓの絶対値が判定値以下であることが認められた場合には、吹き返し補正分fbc(k)の全量を筒内噴射量fidに算入しても、トルクショックが発生するおそれはないと判断できる。よって、この場合には、次に、上記（１５）式に基づいて、筒内噴射量fidが算出され（ステップ１４０）、更に、その算出された量の燃料が筒内インジェクタ２４から実際に噴射される（ステップ１４２）。この筒内噴射時期は、吸気弁閉弁後の圧縮行程中とされる。
【０１１４】
一方、上記ステップ１３８において、変化率Ｓの絶対値が判定値より大きいことが認められた場合には、吹き返し補正分fbc(k)の全量を筒内噴射量fidに算入すると、トルクショックが発生するおそれがあると判断できる。この場合には、トルクショックを確実に回避するべく、吹き返し補正分fbc(k)の一部のみを算入して筒内噴射量fidが算出される（ステップ１４４）。具体的には、上記（１５）式中のfbc(k)を、０＜β＜１を満足する所定の係数βを乗じたβ・fbc(k)に変更した式に基づいて、筒内噴射量fidが算出される。その後、その算出された量の燃料が筒内インジェクタ２４から実際に噴射される（ステップ１４６）。この筒内噴射時期は、吸気弁閉弁後の圧縮行程中とされる。このような処理によれば、トルク段差が少なくなり、トルクショックの発生が確実に防止される。
【０１１５】
本実施形態では、上述したように、吹き返し燃料量fbを特に精度良く算出することができるので、これに基づいて得られる吹き返し補正分fbc(k)の精度も極めて高いものとなる。このため、吹き返しの影響による空燃比のずれを特に高い精度で補正することができる。
【０１１６】
なお、上述した吹き返し燃料量fbの算出処理や、筒内噴射量fidの算出処理においては、付着率Rp、Rv、Rc、残留率Pp、Pv、Pc、噴き分け率γ、筒内残留燃料量frについては上記ステップ１３２で算出された値を用いても良いが、更なる精度向上を図るべく、最新の情報に基づいて再計算した値を用いることとしてもよい。
【０１１７】
上記ステップ１４６の処理により、圧縮工程中の筒内噴射が、吹き返し補正分fbc(k)の一部（β・fbc(k)）のみが算入された筒内噴射量fidで行われた場合には、次に、吹き返し補正分fbc(k)が正であるか負であるかが判別される（ステップ１４８）。
【０１１８】
上記ステップ１４８において、吹き返し補正分fbc(k)が正の値であることが認められた場合には、次に、触媒酸素吸蔵量が過多であるか否かが判別される（ステップ１５０）。ECU６０は、別のルーチンにおいて、公知の手法により触媒酸素吸蔵量を推定しており、その推定値に基づいて上記ステップ１５０の判別を行う。その結果、触媒酸素吸蔵量が過多であることが認められた場合には、上記ステップ１４４において筒内噴射量fidに算入されなかった、吹き返し補正分の残部(1-β)・fbcに相当する量の燃料を、排気行程に合わせて、筒内インジェクタ２４から噴射する（ステップ１５２）。これにより、触媒酸素吸蔵量を減少させることができ、触媒酸素吸蔵量を好ましい範囲へ近づけることができる。
【０１１９】
なお、上記ステップ１５２における排気行程での筒内噴射は、排気下死点以後に限らず、排気弁開弁後なら排気下死点以前から開始してもよい。
また、ステップ１５０において触媒酸素吸蔵量が過多であるとは認められなかった場合には、排気行程での燃料噴射を行う必要はないと考えられるので、上記ステップ１５２の処理は飛ばされる。
【０１２０】
一方、上記ステップ１４８において、吹き返し補正分fbc(k)が負の値であることが認められた場合には、次に、触媒酸素吸蔵量が不足しているか否かが上記触媒酸素吸蔵量の推定値に基づいて判別される（ステップ１５４）。その結果、触媒酸素吸蔵量が不足していることが認められた場合には、上記ステップ１４４において筒内噴射量fidに算入されなかった、吹き返し補正分の残部(1-β)・fbcに相当する量を、次回の処理サイクルにおいてポート噴射量fipに算入することにより、内燃機関６の次サイクルのポート噴射量fipを減量させることとされる。そこで、この場合には、次回の処理サイクルにおいて吹き返し補正分の残部(1-β)・fbcをポート噴射量fipに算入すべきことを表すべく、フラグFに１がセットされ（ステップ１５６）、その後、今回の処理サイクルを終了する。
【０１２１】
なお、ステップ１５４において触媒酸素吸蔵量の不足が認められなかった場合には、内燃機関６の次サイクルのポート噴射量fipを減量させる必要はないと考えられるので、次回の処理サイクルにおいてポート噴射量fipを通常通り算出すべきことを表すべく、フラグFに０がセットされ（ステップ１５８）、その後、今回の処理サイクルを終了する。
【０１２２】
また、吹き返し補正分fbc(k)が正の値であった場合や、更には、吹き返し補正分fbc(k)の全部が筒内噴射量fidに算入されて筒内噴射噴射が行われた場合にも、次回の処理サイクルにおいてポート噴射量fipを通常通り算出すべきことを表すべく、フラグFに０がセットされ（ステップ１５８）、その後、今回の処理サイクルを終了する。
【０１２３】
さて、上記ステップ１５６においてフラグFに１がセットされた場合には、次回の処理サイクルの実行時、上記ステップ１３０においてF＝１の成立が認められる。この場合には、前サイクルの吹き返し補正分の残部(1-β)・fbcを算入してポート噴射量fipが算出される（ステップ１６０）。このステップ１６０の処理は、(1-β)・fbcを算入すること以外は、上記ステップ１３２と同様にして行われる。
【０１２４】
上記ステップ１６０の処理によれば、吹き返し補正分fbcが負の場合において、前サイクルで減量し足りなかった、吹き返し補正分の残部(1-β)・fbcを、今サイクルにおいてポート噴射量fipから差し引くことができる。これにより、触媒酸素吸蔵量を増加させることができ、触媒酸素吸蔵量を好ましい範囲へ近づけることができる。
【０１２５】
ところで、上述した実施の形態２では、吹き返し補正分fbcが負の場合であって、吹き返し補正分の残部(1-β)・fbcを次サイクルの燃料噴射量に算入する場合、これをポート噴射量fipのみに算入しているが、筒内噴射量fidに算入するようにしてもよい。
【０１２６】
また、上述した実施の形態２では、吹き返し補正分fbcが負の場合の吹き返し補正分の残部(1-β)・fbcを、次サイクルの燃料噴射量にのみ算入しているが、次サイクル以降の複数のサイクルに更に分散させて算入してもよい。これにより、トルク段差を更に分散し縮小することができる。
【０１２７】
また、上述した実施の形態２においては、機関回転数NE、負荷率KL、バルブタイミングVT、バルブリフト量VL、全噴射量fiおよび噴き分け率γが前記第１の発明における「所定のパラメータ」に相当している。また、ECU６０が、上記ステップ１３６の処理を実行することにより前記第１の発明における「吹き返し燃料量推定手段」が、上記ステップ１３２および１４０の処理を実行することにより、前記第１の発明における「噴射量算出手段」が、それぞれ実現されている。
【０１２８】
また、上述した実施の形態２においては、ECU６０が、上記ステップ１４０の処理を実行することにより前記第４の発明における「筒内噴射量のみに算入する手段」が、上記ステップ１４２および１４６の処理を実行することにより前記第５の発明における「噴射時期制御手段」が、上記ステップ１４４の処理を実行することにより前記第６の発明における「一部算入手段」が、上記ステップ１５２の処理を実行することにより前記第７の発明における「排気行程噴射手段」が、上記ステップ１６０の処理を実行することにより前記第８の発明における「噴射量減量手段」が、それぞれ実現されている。
【図面の簡単な説明】
【０１２９】
【図１】本発明の実施の形態１のシステム構成を説明するための図である。
【図２】本発明の実施の形態１で用いる燃料挙動モデルを説明するための図である。
【図３】本発明の実施の形態１で用いる燃料挙動モデルを説明するための図である。
【図４】吹き返し燃料量fbと筒内噴射開始時期TDとの関係を、他のパラメータを一定にした条件下で表したグラフである。
【図５】吹き返し燃料量fbと噴き分け率γとの関係を、他のパラメータを一定にした条件下で表したグラフである。
【図６】本発明の実施の形態１において実行されるルーチンのフローチャートである。
【図７】本発明の実施の形態２において実行されるルーチンのフローチャートである。
【符号の説明】
【０１３０】
６    内燃機関
１０  燃焼室
１２  吸気弁
１４  排気弁
１６  点火プラグ
１８  吸気ポート
２０  排気ポート
２２  ポートインジェクタ
２４  筒内インジェクタ
３０  吸気通路
３３  エアフロメータ
３６  スロットルバルブ
４０  排気通路
４２  触媒
４４  空燃比センサ
４６  クランク角センサ
６０  ＥＣＵ
fip  ポート噴射量
fid  筒内噴射量
fwp  ポート付着量
fwv  吸気弁付着量
fwc  筒内付着量
fcp  筒内燃料吸入量
fcd  筒内燃料気化量
Pp  ポート残留率
Rp  ポート付着率
Pv  吸気弁残留率
Rv  吸気弁付着率
Pc  筒内残留率
Rc  筒内付着率
fb  吹き返し燃料量
fr  筒内残留燃料量
fex  排出燃料量
fbc  吹き返し補正分
ｍ  筒内吸入空気量
α  目標空燃比
γ  噴き分け率
【図１】
0004470832.tif 000010
【図２】
0004470832.tif 000011
【図３】
0004470832.tif 000012
【図４】
0004470832.tif 000013
【図５】
0004470832.tif 000014
【図６】
0004470832.tif 000015
【図７】
0004470832.tif 000016