"外部から取り込んだ空気を気筒内に導入する吸気通路と、前記吸気通路内の空気を圧縮するコンプレッサを有する過給機と、前記コンプレッサの下流側の前記吸気通路内に配置されたスロットル弁と、前記スロットル弁の下流にて前記吸気通路と前記気筒との接続部を連通状態または遮断状態にする吸気弁とを備えた内燃機関に適用され、前記吸気通路内を流れる空気の挙動を表す物理モデルに基づいて前記気筒内に導入されている空気の量である筒内空気量を推定する内燃機関の空気量推定装置であって、前記コンプレッサから前記スロットル弁までの前記吸気通路であるスロットル弁上流部内の物理モデルであるスロットル弁上流部モデルと、前記スロットル弁から前記吸気弁までの前記吸気通路であるスロットル弁下流部内の物理モデルであるスロットル弁下流部モデルと、に基づいて前記スロットル弁上流部内の空気の圧力であるスロットル弁上流圧力および前記スロットル弁下流部内の空気の圧力であるスロットル弁下流圧力をそれぞれ推定する第1圧力推定手段と、前記コンプレッサから前記吸気弁までの前記吸気通路である結合部内の物理モデルである結合部モデルに基づいて前記結合部内の空気の圧力である結合部内圧力を推定する第2圧力推定手段と、前記第2圧力推定手段により推定された前記結合部内圧力に基づいて、該結合部内圧力より小さい値となるように前記スロットル弁下流圧力を推定する第3圧力推定手段と、所定の選択条件を満足するか否かを判定する選択条件判定手段と、前記選択条件を満足しないと判定されるときには前記第1圧力推定手段により推定された前記スロットル弁下流圧力に基づいて前記筒内空気量を推定し、前記選択条件を満足すると判定されるときには前記第3圧力推定手段により推定された前記スロットル弁下流圧力に基づいて前記筒内空気量を推定する筒内空気量推定手段と、を備えることを特徴とする内燃機関の空気量推定装置。前記第3圧力推定手段は、前記スロットル弁による圧力損失を算出する圧力損失算出手段と、前記結合部内圧力と前記圧力損失とに基づいて補正値を算出する補正値算出手段と、前記結合部内圧力から前記補正値を差し引くことにより前記スロットル弁下流圧力を算出するスロットル弁下流圧力算出手段と、を含むことを特徴とする請求項1記載の内燃機関の空気量推定装置。上記公報に開示された従来の空気量推定装置では、スロットル弁開度が閾値スロットル弁開度より小さいときは、インタークーラ部スロットル弁上流部内の空気に関する保存則に基づいて構築されたインタークーラモデルスロットル弁上流部モデルと、吸気管部スロットル弁下流部内の空気に関する保存則に基づいて構築された吸気管モデルスロットル弁下流部モデルと、に基づいて吸気管部内圧力スロットル弁下流圧力を推定する。一方、スロットル弁開度が閾値スロットル弁開度より大きいときは、過給機のコンプレッサから吸気弁までの吸気通路である結合部内の空気に関する保存則に基づいて構築された結合部モデルに基づいて吸気管部内圧力を推定する。これにより、スロットル弁開度が比較的大きいためにインタークーラ部内圧力や吸気管部内圧力の変化によりスロットル通過空気流量が短い時間内に大きく変化しやすい状態において、スロットル通過空気流量を所定時間の間一定であると仮定することを必要としない結合部モデルにより吸気管部内圧力を推定することができるので、計算負荷を増大させることなく、吸気管部内圧力を高い精度にて推定することができる。この結果、筒内空気量を高い精度にて推定することができる。しかしながら、本発明者らの研究によれば、上記公報に開示された従来の空気量推定装置を用いて内燃機関の空燃比を制御した場合、スロットル弁上流部モデルおよびスロットル弁下流部モデルを用いて筒内空気量を推定する状態と、結合部モデルを用いて筒内空気量を推定する状態との切り替え時に、空燃比の乱れを生ずる場合がある。この点で、上記従来の技術には、改善の余地が残されている。本発明は、上述のような課題を解決するためになされたもので、過給機を備える内燃機関において計算負荷を増大させることなく筒内空気量を高い精度にて推定することのできる内燃機関の空気量推定装置を提供することを目的とする。第1の発明は、上記の目的を達成するため、外部から取り込んだ空気を気筒内に導入する吸気通路と、前記吸気通路内の空気を圧縮するコンプレッサを有する過給機と、前記コンプレッサの下流側の前記吸気通路内に配置されたスロットル弁と、前記スロットル弁の下流にて前記吸気通路と前記気筒との接続部を連通状態または遮断状態にする吸気弁とを備えた内燃機関に適用され、前記吸気通路内を流れる空気の挙動を表す物理モデルに基づいて前記気筒内に導入されている空気の量である筒内空気量を推定する内燃機関の空気量推定装置であって、前記コンプレッサから前記スロットル弁までの前記吸気通路であるスロットル弁上流部内の物理モデルであるスロットル弁上流部モデルと、前記スロットル弁から前記吸気弁までの前記吸気通路であるスロットル弁下流部内の物理モデルであるスロットル弁下流部モデルと、に基づいて前記スロットル弁上流部内の空気の圧力であるスロットル弁上流圧力および前記スロットル弁下流部内の空気の圧力であるスロットル弁下流圧力をそれぞれ推定する第1圧力推定手段と、前記コンプレッサから前記吸気弁までの前記吸気通路である結合部内の物理モデルである結合部モデルに基づいて前記結合部内の空気の圧力である結合部内圧力を推定する第2圧力推定手段と、前記第2圧力推定手段により推定された前記結合部内圧力に基づいて、該結合部内圧力より小さい値となるように前記スロットル弁下流圧力を推定する第3圧力推定手段と、所定の選択条件を満足するか否かを判定する選択条件判定手段と、前記選択条件を満足しないと判定されるときには前記第1圧力推定手段により推定された前記スロットル弁下流圧力に基づいて前記筒内空気量を推定し、前記選択条件を満足すると判定されるときには前記第3圧力推定手段により推定された前記スロットル弁下流圧力に基づいて前記筒内空気量を推定する筒内空気量推定手段と、を備えることを特徴とする。また、第2の発明は、第1の発明において、前記第3圧力推定手段は、前記スロットル弁による圧力損失を算出する圧力損失算出手段と、前記結合部内圧力と前記圧力損失とに基づいて補正値を算出する補正値算出手段と、前記結合部内圧力から前記補正値を差し引くことにより前記スロットル弁下流圧力を算出するスロットル弁下流圧力算出手段と、を含むことを特徴とする。第1の発明によれば、スロットル弁上流部モデルおよびスロットル弁下流部モデルを用いて筒内空気量を推定する個別モデル状態と、結合部モデルを用いて筒内空気量を推定する結合部モデル状態との切り替え時に、スロットル弁下流圧力の推定値が不連続的に変化することを防止することができる。その結果、個別モデル状態と結合部モデル状態との切り替え時においても、筒内空気量をより高い精度で推定することができ、空燃比の乱れを確実に抑制することができる。第2の発明によれば、結合部内圧力に対して精度良く補正を施して、正確なスロットル弁下流圧力を算出することができる。実施の形態1.以下、図面を参照して本発明の内燃機関の空気量推定装置の実施の形態について説明するが、本実施形態の内燃機関の空気量推定装置は、特開2006-152899号公報の実施形態に記載された空気量推定装置と共通する部分が多いため、以下では、同公報記載の空気量推定装置と共通する事項については説明を簡略化または省略する。本実施形態の空気量推定装置は、以下の説明で言及していない点については、特開2006-152899号公報記載の空気量推定装置と同様である。図1は、本発明の実施の形態1のシステム構成を説明するための図である。図1に示すように、本実施形態のシステムは、内燃機関10と、この内燃機関10を制御するためのECUElectronicControlUnit70とを備えている。内燃機関10の各気筒には、吸気弁32と、排気弁35と、点火プラグ37と、燃料インジェクタ39とが設けられている。吸気系統40は、インテークマニホールド41と、インテークマニホールド41に連通したサージタンク42と、エアフィルタ44と、過給機91のコンプレッサ91aと、インタークーラ45と、スロットル弁46とを備えている。なお、コンプレッサ91aの出口下流からスロットル弁46までの吸気通路は、インタークーラ45とともに、スロットル弁上流部としてのインタークーラ部を構成している。更に、スロットル弁46から吸気弁32までの吸気通路は、スロットル弁下流部としての吸気管部を構成している。また、コンプレッサ91aの出口下流から吸気弁32までの吸気通路インタークーラ部および吸気管部は、結合部を構成している。排気系統50は、エキゾーストマニホールドを含む排気管51と、排気管51内に配設された過給機91のタービン91bと、タービン91bの下流の排気管51に配設された三元触媒装置52とを備えている。このような配置により、過給機91のタービン91bは排気ガスのエネルギーにより回転し、コンプレッサ91aを駆動する。これにより、コンプレッサ91aが吸気通路内の空気を圧縮する。また、図示を省略するが、本実施形態のシステムは、特開2006-152899号公報の図2に記載されたシステムと同様の各種のセンサを備えている。上記構成を有する本実施形態の空気量推定装置は、図2Aに示したように、所定の選択条件が満足されないときには、インタークーラ部内の空気の圧力Picおよび温度Ticと、吸気管部内の空気の圧力Pmおよび温度Tmと、を推定するための物理モデルとして、インタークーラ部内の空気に関する保存則に基づいて構築された物理モデルインタークーラモデルM5と、吸気管部内の空気に関する保存則に基づいて構築された物理モデル吸気管モデルM6と、を採用する。上記選択条件には、例えば、スロットル弁開度が所定の閾値スロットル弁開度より大きい、という条件が含まれる。一方、選択条件が満足されるときは、インタークーラ部内の空気の圧力や吸気管部内の空気の圧力の変化によりスロットル弁46の周囲を通過する空気の流量スロットル通過空気流量が短い時間内に大きく変化しやすい。そこで、本実施形態の空気量推定装置は、図2Bに示したように、スロットル弁開度が上記閾値スロットル弁開度より大きいとき、インタークーラ部内の空気の圧力Picおよび温度Ticと、吸気管部内の空気の圧力Pmおよび温度Tmと、を推定するための物理モデルとして、結合部内の空気に関する保存則に基づいて構築された物理モデルインタークーラ吸気管結合モデルIC吸気管結合モデルM8を採用する。より具体的に述べると、本実施形態の空気量推定装置は、選択条件が満足されないとき、図3に示した電子制御スロットル弁モデルM1、スロットルモデルM2、吸気弁モデルM3、コンプレッサモデルM4、インタークーラモデルスロットル弁上流部モデルM5、吸気管モデルスロットル弁下流部モデルM6、吸気弁モデルM7および電子制御スロットル弁ロジックA1を用いて筒内空気量を推定する。この場合には、スロットルモデルM2、インタークーラモデルM5および吸気管モデルM6により、インタークーラ部内圧力すなわち、過給機91からスロットル弁46までの吸気通路内の空気の圧力であるスロットル弁上流圧力Pic、インタークーラ部内温度すなわち、コンプレッサ91aからスロットル弁46までの吸気通路内の空気の温度であるスロットル弁上流温度Tic、吸気管部内圧力すなわち、スロットル弁46から吸気弁32までの吸気通路内の空気の圧力であるスロットル弁下流圧力Pmおよび吸気管部内温度すなわち、スロットル弁46から吸気弁32までの吸気通路内の空気の温度であるスロットル弁下流温度Tmが推定される。このような推定の処理が実行されることは、第1圧力推定手段の機能が達成されることに対応している。一方、本実施形態の空気量推定装置は、選択条件が満足されるとき、図4に示した電子制御スロットル弁モデルM1、吸気弁モデルM3、コンプレッサモデルM4、吸気弁モデルM7、IC吸気管結合モデル結合部モデルM8および電子制御スロットル弁ロジックA1を用いて筒内空気量を推定する。この場合には、IC吸気管結合モデルM8により、結合部内の空気の圧力である結合部内圧力Picmおよび結合部内の空気の温度である結合部内温度Ticmが推定される。このような推定の処理が実行されることは、第2圧力推定手段の機能が達成されることに対応している。以下の説明では、インタークーラモデルM5および吸気管モデルM6を用いて筒内空気量が推定される状態すなわち、選択条件が満足されないときを個別モデル使用状態と称し、IC吸気管結合モデルM8を用いて筒内空気量が推定される状態すなわち、選択条件が満足されるときを結合部モデル使用状態と称する。特開2006-152899号公報記載の空気量推定装置では、結合部モデル使用状態においては、インタークーラ部内圧力Picおよび吸気管部内圧力Pmが上記結合部内圧力Picmに設定されるとともに、インタークーラ部内温度Ticおよび吸気管部内温度Tmが上記結合部内温度Ticmに設定される。すなわち、インタークーラ部内圧力Picおよび吸気管部内圧力Pmが共に結合部内圧力Picmと等しいものとして推定され、また、インタークーラ部内温度Ticおよび吸気管部内温度Tmが共に結合部内温度Ticmに等しいものとして推定されている。個別モデル使用状態、結合部モデル使用状態の何れにおいても、推定された吸気管部内圧力Pmを用いて、吸気弁モデルM7に基づき、筒内流入空気流量mcが算出される。そして、この筒内流入空気流量mcに吸気弁開弁時間吸気弁32が開弁してから閉弁するまでの時間Tintを乗じて予測筒内空気量KLfwdが算出される。このような算出処理が実行されることは、筒内空気量推定手段の機能が達成されることに対応している。この予測筒内空気量KLfwdに基づいて燃料インジェクタ39からの燃料噴射量を算出することにより、空燃比が制御される。特開2006-152899号公報記載の空気量推定装置では、個別モデル使用状態から結合部モデル使用状態に切り替わるときに、吸気管部内圧力Pmの推定値が不連続的に変化することにより、予測筒内空気量KLfwdが急変して、空燃比が乱れる場合がある。本発明者らの最近の研究によれば、この原因は、次のようなものであると考えられる。上述したように、特開2006-152899号公報記載の空気量推定装置では、結合部モデル使用状態において、インタークーラ部内圧力Picおよび吸気管部内圧力Pmが共に結合部内圧力Picmと等しいものとして推定している。しかしながら、結合部モデル使用状態においても、実際には、吸気管部内圧力Pmがインタークーラ部内圧力Picより僅かに小さい場合がある。このような場合、結合部内圧力Picmは、実際の吸気管部内圧力Pmと実際のインタークーラ部内圧力Picとの間の値となる。すなわち、実際の吸気管部内圧力Pmは、結合部内圧力Picmより小さくなる。このため、吸気管部内圧力Pmが結合部内圧力Picmに等しいとする推定に基づいて筒内空気量を推定すると、実際の筒内空気量に対して誤差を生じ、その結果、空燃比が乱れると考えられる。上記の事項に鑑みて、本実施形態の空気量推定装置では、結合部モデル使用状態において、実際の吸気管部内圧力Pmが結合部内圧力Picmより僅かに小さいことを考慮し、結合部内圧力Picmに基づいて、実際の値により近い吸気管部内圧力Pmを推定することとした。そして、その実際の値により近い吸気管部内圧力Pm<結合部内圧力Picmに基づいて、筒内空気量を推定することとした。以下の説明では、実際の値により近いことを正確なと称する。以下、図5を参照して、結合部内圧力Picmに基づいて正確な吸気管部内圧力Pmを推定する方法について説明する。吸気弁32の周囲を通過して気筒内に流入する空気の流量である筒内流入空気流量mcは、吸気管部内圧力Pmに比例する。この比例関係は、吸気弁モデルM3により求めることができ、次式で表される。mc=...上記式において、cおよびdの値は、ECU70に予め記憶された吸気弁モデルM3が備えるテーブルに基づき、現在のエンジン回転速度NE等に応じて求められる。吸気温度Taは、図示しない温度センサにより検出される。吸気管部内温度Tmは、結合部内温度Ticmにほぼ等しいので、結合部内温度Ticmの値を用いることができる。このようにして求めることのできるc、d、TaおよびTmの値から算出される値であるaおよびbを用いると、上記式は、次式のように書き直すことができる。mc=aPmb...以下では、数式の記載を簡単化するための便宜上、上記式に代えて上記式を用いる。図5中の右上がりの直線は、グラフの横軸を吸気管部内圧力Pmとし、縦軸を筒内流入空気流量mcとして、上記式の関係を表したものである。上記式の関係を以下筒内流入空気流量特性と称する。一方、図5中の曲線は、インタークーラ部内圧力Picをグラフの横軸とし、コンプレッサ91aから流出する空気の流量であるコンプレッサ流出空気流量mcmをグラフの縦軸をとして、両者の関係を表したものである。この関係は、コンプレッサモデルM4により、求めることができる。すなわち、コンプレッサモデルM4では、ECU70に予め記憶されたテーブルに基づき、現在の吸気圧力Paおよびコンプレッサ回転速度Ncmに応じて、この関係を求めることができる。この関係、すなわち図5中の曲線で示す関係を以下コンプレッサ流量特性と称する。以下の説明では、求めるべき正確な吸気管部内圧力PmをPmacとし、正確なインタークーラ部内圧力PicをPicacとする。正確な吸気管部内圧力Pmacは結合部内圧力Picmより小さく、正確なインタークーラ部内圧力Picacは結合部内圧力Picmより大きい。正確なインタークーラ部内圧力Picacと正確な吸気管部内圧力Pmacとの差は、スロットル弁46による圧力損失圧力低下量ΔPlossに相当する。圧力損失ΔPlossの算出方法については後述する。図5に示すように、結合部内圧力Picmと正確な吸気管部内圧力Pmacとの差をαとし、正確なインタークーラ部内圧力Picacと結合部内圧力Picmとの差をβとすると、次式が成り立つ。αβ=ΔPloss...図5では、説明の便宜上、αおよびβを大きく誇張して描いているが、結合部モデル使用状態における実際のインタークーラ部内圧力Picと吸気管部内圧力Pmとは近いので、実際のαおよびβは小さい値である。図5に示すように、各記号を次のように定義する。結合部内圧力Picmをインタークーラ部内圧力Picとしてコンプレッサ流量特性に当てはめた場合のコンプレッサ流出空気流量mcmと、結合部内圧力Picmを吸気管部内圧力Pmとして筒内流入空気流量特性に当てはめた場合の筒内流入空気流量mcとの差をΔmとする。正確なインタークーラ部内圧力Picacをインタークーラ部内圧力Picとしてコンプレッサ流量特性に当てはめた場合のコンプレッサ流出空気流量mcm、すなわち正確なコンプレッサ流出空気流量mcmをmcmacとする。正確な吸気管部内圧力Pmacを吸気管部内圧力Pmとして筒内流入空気流量特性に当てはめた場合の筒内流入空気流量mc、すなわち正確な筒内流入空気流量mcをmcacとする。正確なコンプレッサ流出空気流量mcmacと、正確な筒内流入空気流量mcacとの差をΔm’とする。インタークーラ部内圧力Picが結合部内圧力Picmに等しい点でのコンプレッサ流量特性の傾きを-a’とするただし、a’>0。このように定義すると、図5に示す関係より、次式が成り立つ。Δm’=Δmaα-a’β...また、上述の定義から、次式が成り立つ。Δm’=mcmac-mcac...mcac=aPmacb=ab...インタークーラ部内圧力Picと吸気管部内圧力Pmとが近い場合には、インタークーラ部内圧力Picの変化量と吸気管部内圧力Pmの変化量とはほぼ等しくなる。このため、スロットル弁46の周囲を通過する空気の流量であるスロットル通過空気流量をmt、インタークーラ部の容積をVic、吸気管部の容積をVmとすると、次式が成り立つ。Vic=Vm...上記式、式および式から、次式が導かれる。mt=Δm’ab...スロットルモデルM2は、スロットル通過空気流量mtを下記式により推定することができる。下記式において、Ctは流量係数、Atはスロットル弁開度θtに応じて変化するスロットル開口断面積、Rは気体定数である。Φは所定の関数である。mt=CtAtPic√Φ...CtAtの値は、スロットル弁46の開度θtに基づいて決定できることが経験的に知られている。スロットルモデルM2は、スロットル弁開度θtとCtAtの値との関係を規定するテーブルをECU70に記憶させており、スロットル弁開度θtに基づいてCtAtの値を算出することができる。更に、スロットルモデルM2は、関数Φの値を算出するためのテーブルをECU70に記憶させており、PmPicの値に基づいてΦの値を算出することができる。上記式において、Picの値としてPicac=Picmβを用い、PmPicの値としてPmacPicacを用いることにより、次式が得られる。mt=CtAt√Φ...βの値が小さいため、PmacPicacの値は、次式のように近似できる。PmacPicac=≒∴PmacPicac≒Picm...また、結合部モデル使用状態における実際のインタークーラ部内温度Ticは、結合部内温度Ticmにほぼ等しいので、結合部内温度Ticmを近似的に用いることができる。これらの近似によれば、上記式は、次式に代えることができる。mt=CtAt√ΦPicm...上記式、式、式および式を連立してmt、βおよびΔm’を消去することにより、αの値を求めることができる。このαを結合部内圧力Picmから差し引くことにより、正確な吸気管部内圧力Pmacを算出することができる。本実施形態では、結合部モデル使用状態のとき、このようにして求められた正確な吸気管部内圧力Pmacを吸気管部内圧力Pmの推定値として用いる。これにより、個別モデル使用状態との切り替わり時に、吸気管部内圧力Pmの推定値が不連続的に変化することを防止することができる。その結果、筒内空気量をより高い精度で推定することができ、空燃比の乱れを確実に抑制することができる。図6は、上述した原理に基づき、結合部モデル使用状態において正確な筒内流入空気流量mcacを推定するためにECU70が実行するルーチンのフローチャートである。図6に示すルーチンは、結合部モデル使用状態において、第2圧力推定手段のIC吸気管結合モデルM8によって結合部内圧力Picmおよび結合部内温度Ticmが推定された後、引き続いて実行される。図6に示すルーチンによれば、まず、スロットル弁46による圧力損失ΔPlossが算出されるステップ100。このステップ100での圧力損失ΔPlossの算出原理は、特開2005-69020号公報で説明された原理を過給機付き内燃機関に適用したものであるので、ここでは説明を簡略化する。図7は、圧力損失ΔPlossの算出方法を説明するための図である。図7に示すグラフには、筒内流入空気流量特性およびスロットル通過空気流量特性が示されている。筒内流入空気流量特性は、上記式で表される。スロットル通過空気流量特性は、上記式の関係を表したものである。吸気管部内圧力Pmがインタークーラ部内圧力Picに等しい場合には、スロットル弁46の上流と下流との圧力差がゼロであるので、スロットル通過空気流量mtはゼロとなる。このため、図7に示すように、スロットル通過空気流量特性は、吸気管部内圧力Pmがインタークーラ部内圧力Picに等しいときにゼロとなる。内燃機関10の定常運転時には、筒内流入空気流量mcとスロットル通過空気流量mtとは一致する。したがって、図7に示すように、定常運転時における吸気管部内圧力Pmは、筒内流入空気流量特性とスロットル通過空気流量特性との交点として求めることができる。そして、圧力損失ΔPlossは、インタークーラ部内圧力Picから定常運転時における吸気管部内圧力Pmを差し引いた値に相当する。上記ステップ100では、結合部内圧力Picmおよび結合部内温度Ticmの値をインタークーラ部内圧力Picおよびインタークーラ部内温度Ticとしてそれぞれ用いて、上述した方法を適用することにより、圧力損失ΔPlossが算出される。次いで、αの値を算出するための準備として、前述したa、b、a’およびΔmの値を算出する処理が実行されるステップ102。前述したように、aおよびbの値は吸気弁モデルM3に基づいて求められ、a’の値はコンプレッサモデルM4に基づいて求められる。また、Δmの値は、コンプレッサ流量特性に結合部内圧力Picmを当てはめた値と、筒内流入空気流量特性に結合部内圧力Picmを当てはめた値との差として算出される。続いて、αの値を算出する処理が実行されるステップ104。このステップ104では、結合部内圧力Picm、結合部内温度Ticm、圧力損失ΔPloss、a、b、a’およびΔmの値を上記式、式、式および式に代入することにより、αの値が算出される。次いで、結合部内圧力Picmの値から、上記ステップ104で算出されたαの値を差し引くことにより、正確な吸気管部内圧力Pmacが算出されるステップ106。そして、この正確な吸気管部内圧力Pmacを筒内流入空気流量特性に当てはめることにより、正確な筒内流入空気流量mcacが算出されるステップ108。すなわち、次式により正確な筒内流入空気流量mcacが算出される。mcac=aPmacb...結合部モデル使用状態においては、上述した図6のルーチンの処理によって正確な筒内流入空気流量mcacを算出し、この正確な筒内流入空気流量mcacに基づいて予測筒内空気量KLfwdを算出し、この予測筒内空気量KLfwdに基づいて燃料インジェクタ39からの燃料噴射量を算出する。以上説明したように、本実施形態の空気量推定装置は、選択条件が満足されないときスロットル弁開度が閾値スロットル弁開度より小さいときは、インタークーラ部スロットル弁上流部内の空気に関する保存則に基づいて構築されたインタークーラモデルスロットル弁上流部モデルM5と、吸気管部スロットル弁下流部内の空気に関する保存則に基づいて構築された吸気管モデルスロットル弁下流部モデルM6と、に基づいて吸気管部内圧力スロットル弁下流圧力を推定する。すなわち、個別モデル使用状態とする。一方、選択条件が満足されるときスロットル弁開度が閾値スロットル弁開度より大きいときは、コンプレッサ91aから吸気弁32までの吸気通路である結合部内の空気に関する保存則に基づいて構築されたIC吸気管結合モデル結合部モデルM8に基づいて結合部内圧力Picmを推定する。すなわち、結合部モデル使用状態とする。そして、その結合部内圧力Picmに基づいて、吸気管部内圧力を推定する。更に、個別モデル使用状態、結合部モデル使用状態の何れの場合も、推定された吸気管部内圧力に基づいて筒内空気量を推定する。これにより、スロットル弁開度が比較的大きいためにインタークーラ部内圧力や吸気管部内圧力の変化によりスロットル通過空気流量が短い時間内に大きく変化しやすい状態において、スロットル通過空気流量を所定時間の間一定であると仮定することを必要としないIC吸気管結合モデルM8を用いて吸気管部内圧力を推定することができるので、計算負荷を増大させることなく、吸気管部内圧力を高い精度にて推定することができる。この結果、筒内空気量を高い精度にて推定することができる。また、本実施形態の空気量推定装置の実施形態は、結合部モデル使用状態においては、実際の吸気管部内圧力Pmが結合部内圧力Picmより僅かに小さいことを考慮し、結合部内圧力Picmに基づいて、正確な吸気管部内圧力Pmacを算出し、その正確な吸気管部内圧力Pmacを吸気管部内圧力Pmの推定値として用いて、筒内空気量予測筒内空気量KLfwdを推定する。これにより、個別モデル使用状態と結合部モデル使用状態とが切り替わるときに、吸気管部内圧力Pmの推定値が不連続的に変化することを確実に防止することができるので、予測筒内空気量KLfwdが急変することがない。このため、個別モデル使用状態と結合部モデル使用状態との切り替わり時においても、空燃比が乱れが生ずることを確実に防止することができる。以上説明した実施の形態では、個別モデル使用状態と結合部モデル使用状態との切り替えを判断するための選択条件が、スロットル弁開度が所定の閾値スロットル弁開度より大きい、との条件を含む所定の条件であるものとして説明したが、本発明における選択条件はこれに限定されるものではない。例えば、インタークーラモデルM5および吸気管モデルM6を用いて推定される筒内流入空気流量と、IC吸気管結合モデルM8を用いて推定される筒内流入空気流量とを同時に計算し、その両者の差が所定値以下になることを選択条件としてもよい。この場合、スロットル弁開度の増大に伴って両者の差が所定値以下になった場合には個別モデル使用状態から結合部モデル使用状態に切り替え、スロットル弁開度の縮小に伴って両者の差が所定値以下になった場合には結合部モデル使用状態から個別モデル使用状態に切り替えればよい。また、上述した実施の形態1においては、ECU70が、上記ステップ100~106の処理を実行することにより前記第1の発明における第3圧力推定手段が、上記ステップ100の処理を実行することにより前記第2の発明における圧力損失算出手段が、上記ステップ102,104の処理を実行することにより前記第2の発明における補正値算出手段が、上記ステップ106の処理を実行することにより前記第2の発明におけるスロットル弁下流圧力算出手段が、それぞれ実現されている。"