(57)【要約】【課題】  燃焼室内にて燃料と空気を含む混合気の燃焼により発生するＮＯｘ量を、同混合気の周辺状態量を考慮して精度良く推定できるＮＯｘ発生量推定方法を提供すること。【解決手段】  この方法は、混合気が占める領域（燃焼領域）において燃焼により発生するＮＯｘ発生量に大きく影響を与える同混合気の周辺状態量として、ＮＯｘを発生させる材料となる吸気中のガスの濃度（吸気酸素濃度RO2c）、機関に与えられる負荷の程度を表す負荷指標値（燃料噴射量qfinc）、燃焼室内の燃料の霧化の程度を表す霧化指標値（燃料噴射圧力Pcrc）、及び最高火炎温度Tflameの４つを選択する。そして、この４つの周辺状態量と、同４つの周辺状態量と単位燃料量あたりの燃焼発生ＮＯｘ量（燃焼発生ＮＯｘ率RNOx_burn）との関係を規定する所定の実験式と、に基づいてRNOx_burnを求め、同RNOx_burnに前記燃料噴射量qfincを乗じた値をＮＯｘ発生量として推定する。【選択図】  図１０
2005139984.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の燃焼室内にて燃料と空気を含む混合気の燃焼により発生する燃焼発生ＮＯｘ量を、前記燃焼発生ＮＯｘ量に影響を与える、前記混合気の周辺状態量に基づいて推定する内燃機関のＮＯｘ発生量推定方法。
【請求項２】
請求項１に記載の内燃機関のＮＯｘ発生量推定方法において、
前記周辺状態量としての、前記機関に与えられる負荷の程度を表す負荷指標値に少なくとも基づいて前記燃焼発生ＮＯｘ量を推定する内燃機関のＮＯｘ発生量推定方法。
【請求項３】
請求項１又は請求項２に記載の内燃機関のＮＯｘ発生量推定方法において、
前記周辺状態量としての、前記燃焼室内の燃料の霧化の程度を表す霧化指標値に少なくとも基づいて前記燃焼発生ＮＯｘ量を推定する内燃機関のＮＯｘ発生量推定方法。
【請求項４】
内燃機関の燃焼室内にて燃料と空気を含む混合気の燃焼が発生する領域である燃焼領域を推定し、
請求項１乃至請求項３の何れか一項に記載のＮＯｘ発生量推定方法により前記燃焼領域内にて前記混合気の燃焼により発生する前記燃焼発生ＮＯｘ量を推定し、
前記燃焼室内における前記燃焼領域を除いた領域である非燃焼領域内におけるＮＯｘ量を推定するとともに、
前記燃焼発生ＮＯｘ量と、前記非燃焼領域内におけるＮＯｘ量とに基づいて、前記排気通路から外部へ排出される排ガス中のＮＯｘ量を推定する内燃機関のＮＯｘ排出量推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の燃焼室内にて燃料と空気を含む混合気の燃焼により発生する燃焼発生ＮＯｘ量を推定するＮＯｘ発生量推定方法に関する。
【背景技術】
【０００２】
火花点火式内燃機関、ディーゼル機関等の内燃機関においては、排気通路から外部へ排出される排ガス中のＮＯｘ量（以下、「ＮＯｘ排出量」と云うこともある。）を低減する必要があり、そのためには、燃焼室内にて（噴射により霧化した）燃料と空気を含む混合気の燃焼により発生する燃焼発生ＮＯｘ量を低減する必要がある。この燃焼発生ＮＯｘ量を低減するためには、例えば、ＥＧＲ装置により還流されるＥＧＲガス量を増大させること、或いは、燃料噴射時期を遅らせること等により最高火炎温度（最高燃焼温度）を低下させることが有効である。
【０００３】
しかしながら、一方では、燃焼発生ＮＯｘ量を低減するためにＥＧＲガス量を増大させていくと（ディーゼル機関において）微粒子状物質（パティキュレート・マター、ＰＭ）の発生量が増大していくというトレードオフがあり、また、燃焼発生ＮＯｘ量を低減するために燃料噴射時期を遅らせていくと燃費が悪化していくというトレードオフがあることが知られている。
【０００４】
従って、ＰＭ排出量の増大の抑制、燃費の悪化の抑制等を考慮した上で燃焼発生ＮＯｘ量を出来る限り低減するためには、同燃焼発生ＮＯｘ量を機関の運転状態に応じた所定の目標値に制御することが好ましい。他方、燃焼発生ＮＯｘ量を実測することは非常に困難である。よって、燃焼発生ＮＯｘ量を所定の目標値に精度良く制御するためには同燃焼発生ＮＯｘ量を精度良く推定する必要がある。
【０００５】
このため、下記特許文献１に記載の内燃機関の制御装置は、燃焼圧力、及び吸気酸素濃度を、筒内圧力センサ、及び吸気酸素濃度センサによりそれぞれ検出するとともに、これらに基づいて算出した燃焼温度、及び混合気濃度とに基づいて、代表的な公知の燃焼モデルの一つである拡大ゼルドビッチ（ＺＥＬＤＯＶＩＣＨ）機構を用いて上記燃焼発生ＮＯｘ量を推定する。そして、推定された燃焼発生ＮＯｘ量が所定の目標値になるようにＥＧＲガス量、或いは燃料噴射時期等を制御するようになっている。
【特許文献１】特開２００２−３７１８９３号公報
【０００６】
ところで、燃焼室内にて混合気の燃焼により発生する実際の燃焼発生ＮＯｘ量は、例えば、機関に与えられている負荷（駆動トルク）、燃焼前の混合気を構成している燃料の霧化の程度等、のような混合気の周辺状態量に大きく依存する。しかしながら、上記従来の装置においては、上記燃焼発生ＮＯｘ量を推定するにあたり、係る混合気の周辺状態量に対する考慮が全くなされていない。従って、上記燃焼発生ＮＯｘ量を精度良く推定することができず、この結果、上記（実際の）燃焼発生ＮＯｘ量を所定の目標値に精度良く制御することができないという問題があった。
【発明の開示】
【０００７】
本発明は、かかる課題に対処するためになされたものであって、その目的は、内燃機関の燃焼室内にて燃料と空気を含む混合気の燃焼により発生する燃焼発生ＮＯｘ量を、同混合気の周辺状態量を考慮して精度良く推定することができるＮＯｘ発生量推定方法を提供することにある。
【０００８】
本発明による内燃機関のＮＯｘ発生量推定方法は、上記燃焼発生ＮＯｘ量を、同燃焼発生ＮＯｘ量に影響を与える、前記混合気の周辺状態量に基づいて推定する方法である。これによれば、上記燃焼発生ＮＯｘ量に影響を与える混合気の周辺状態量が考慮されて同燃焼発生ＮＯｘ量が推定され得るから、同燃焼発生ＮＯｘ量が精度良く推定され得る。
【０００９】
この場合、前記周辺状態量としての、前記機関に与えられる負荷の程度を表す負荷指標値に少なくとも基づいて前記燃焼発生ＮＯｘ量を推定することが好適である。ここにおいて、前記負荷指標値は前記機関に与えられる負荷の程度を表す値であって、例えば、（一作動サイクルあたりの）燃料噴射量、機関の駆動トルク、燃焼室の内壁面の温度等である。
【００１０】
機関に与えられる負荷が大きくなるほど、同負荷に対抗するため燃焼室内においてより多くの爆発エネルギーを発生させる必要がある。この結果、燃焼室の内壁面温度が上昇し、同内壁面からの輻射熱により燃焼前の混合気が加熱され易くなる。従って、機関に与えられる負荷が大きくなるほど、上記最高火炎温度（最高燃焼温度）が高くなって燃焼発生ＮＯｘ量が増加する。
【００１１】
従って、上記のように、少なくとも前記負荷指標値に基づいて燃焼発生ＮＯｘ量を推定することにより、例えば、同負荷指標値が表す負荷が大きくなるほど燃焼発生ＮＯｘ量がより大きい値になるように推定され得、この結果、同燃焼発生ＮＯｘ量が精度良く推定され得る。
【００１２】
また、上記本発明によるＮＯｘ発生量推定方法においては、前記周辺状態量としての、前記燃焼室内の（噴射された）燃料の霧化の程度を表す霧化指標値に少なくとも基づいて前記燃焼発生ＮＯｘ量を推定することが好適である。ここにおいて、前記霧化指標値は燃焼室内の燃料の霧化の程度を表す値であって、例えば、燃料の噴射圧力、スワール比、燃焼が発生する領域における空気過剰率等である。
【００１３】
噴射された燃料の霧化の程度が大きくなるほど、同噴射された燃料の量に対する同燃料と混ざり合って混合気となる空気量の割合が多くなる。従って、噴射された燃料の霧化の程度が大きくなるほど、混合気が占める領域（即ち、燃焼が発生する領域（燃焼領域））における空気過剰率が大きくなって燃焼発生ＮＯｘ量が増加する。従って、上記のように、少なくとも前記霧化指標値に基づいて燃焼発生ＮＯｘ量を推定することにより、例えば、同霧化指標値が表す霧化の程度が大きくなるほど燃焼発生ＮＯｘ量がより大きい値になるように推定され得、この結果、同燃焼発生ＮＯｘ量が精度良く推定され得る。
【００１４】
また、上記本発明による内燃機関のＮＯｘ発生量推定方法を利用して同内燃機関の排気通路から外部へ排出される排ガス中のＮＯｘ量（以下、「ＮＯｘ排出量」と称呼する。）を推定する本発明による内燃機関のＮＯｘ排出量推定方法は、上述の燃焼が発生する領域である燃焼領域を推定し、上記本発明によるＮＯｘ発生量推定方法により前記燃焼領域内にて前記混合気の燃焼により発生する前記燃焼発生ＮＯｘ量を推定し、前記燃焼室内における前記燃焼領域を除いた領域である非燃焼領域内におけるＮＯｘ量を推定するとともに、前記燃焼発生ＮＯｘ量と、前記非燃焼領域内におけるＮＯｘ量とに基づいて、前記ＮＯｘ排出量を推定する方法である。
【００１５】
排気通路を流れる排ガスの一部を吸気通路へ還流せしめるＥＧＲ装置を備えた内燃機関においては、ＥＧＲ装置を介してＥＧＲガス中のＮＯｘが燃焼室内に還流してくる。加えて、上記燃焼発生ＮＯｘ量は、燃焼室の一部であって燃焼が発生する領域（即ち、上記燃焼領域）内において発生するＮＯｘの量である。従って、燃焼室内における燃焼領域を除いた領域（以下、「非燃焼領域」と称呼する。）においては燃焼後においても前記還流されてきたＮＯｘが残存しているから、排気通路から外部へ排出される排ガス中のＮＯｘ量を精度良く推定するためには、上記燃焼発生ＮＯｘ量のみならず、前記「非燃焼領域内に残存しているＮＯｘ量」をも考慮することが必要である。
【００１６】
以上の知見に基づき、上記のように、前記推定された燃焼領域内にて燃焼により発生する上記燃焼発生ＮＯｘ量のみならず、（燃焼後における）非燃焼領域内におけるＮＯｘ量、従って、上述した「非燃焼領域内に残存しているＮＯｘ量」を考慮してＮＯｘ排出量を推定すると、ＮＯｘ排出量を精度良く推定することができる。
【発明を実施するための最良の形態】
【００１７】
以下、本発明の実施形態に係る内燃機関のＮＯｘ発生量推定方法、及びＮＯｘ排出量推定方法を実施する内燃機関（ディーゼル機関）の制御装置について図面を参照しつつ説明する。
【００１８】
図１は、係る内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）１０に適用したシステム全体の概略構成を示している。このシステムは、燃料供給系統を含むエンジン本体２０、エンジン本体２０の各気筒の燃焼室（筒内）にガスを導入するための吸気系統３０、エンジン本体２０からの排ガスを放出するための排気系統４０、排気還流を行うためのＥＧＲ装置５０、及び電気制御装置６０を含んでいる。
【００１９】
エンジン本体２０の各気筒の上部には燃料噴射弁（噴射弁、インジェクタ）２１が配設されている。各燃料噴射弁２１は、図示しない燃料タンクと接続された燃料噴射用ポンプ２２に燃料配管２３を介して接続されている。燃料噴射用ポンプ２２は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（後述する（指令）基本燃料噴射圧力Pcrbaseに応じた指令信号）により燃料の実際の噴射圧力（吐出圧力）が同指令基本燃料噴射圧力Pcrbaseになるように同燃料を昇圧するようになっている。
【００２０】
これにより、燃料噴射弁２１には、燃料噴射用ポンプ２２から前記基本燃料噴射圧力Pcrbaseまで昇圧された燃料が供給されるようになっている。また、燃料噴射弁２１は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（後述する（指令）燃料噴射量qfinに応じた指令信号）により所定時間だけ開弁し、これにより各気筒の燃焼室内に前記指令基本燃料噴射圧力Pcrbaseにまで昇圧された燃料を前記指令燃料噴射量qfinだけ直接噴射するようになっている。
【００２１】
吸気系統３０は、エンジン本体２０の各気筒の燃焼室にそれぞれ接続された吸気マニホールド３１、吸気マニホールド３１の上流側集合部に接続され同吸気マニホールド３１とともに吸気通路を構成する吸気管３２、吸気管３２内に回動可能に保持されたスロットル弁３３、電気制御装置６０からの駆動信号に応答してスロットル弁３３を回転駆動するスロットル弁アクチュエータ３３ａ、スロットル弁３３の上流において吸気管３２に順に介装されたインタクーラー３４と過給機３５のコンプレッサ３５ａ、及び吸気管３２の先端部に配設されたエアクリーナ３６とを含んでいる。
【００２２】
排気系統４０は、エンジン本体２０の各気筒にそれぞれ接続された排気マニホールド４１、排気マニホールド４１の下流側集合部に接続された排気管４２、排気管４２に配設された過給機３５のタービン３５ｂ、及び排気管４２に介装されたディーゼルパティキュレートフィルタ（以下、「ＤＰＮＲ」と称呼する。）４３を含んでいる。排気マニホールド４１及び排気管４２は排気通路を構成している。
【００２３】
ＤＰＮＲ４３は、コージライト等の多孔質材料から形成されたフィルタ４３ａを備え、通過する排気ガス中のパティキュレートを細孔表面にて捕集するフィルタである。ＤＰＮＲ４３は、担体としてのアルミナに、カリウムＫ，ナトリウムＮａ，リチウムＬｉ，セシウムＣｓのようなアルカリ金属、バリウムＢａ，カルシウムＣａのようなアルカリ土類金属、及びランタンＬａ、イットリウムＹのような希土類金属から選ばれた少なくとも一つを白金とともに担持し、ＮＯｘを吸収した後に同吸収したＮＯｘを放出して還元する吸蔵還元型ＮＯｘ触媒としても機能するようになっている。
【００２４】
ＥＧＲ装置５０は、排気ガスを還流させる通路（ＥＧＲ通路）を構成する排気還流管５１と、排気還流管５１に介装されたＥＧＲ制御弁５２と、ＥＧＲクーラー５３とを備えている。排気還流管５１はタービン３５ｂの上流側排気通路（排気マニホールド４１）とスロットル弁３３の下流側吸気通路（吸気マニホールド３１）を連通している。ＥＧＲ制御弁５２は電気制御装置６０からの駆動信号に応答し、再循環される排気ガス量（排気還流量、ＥＧＲガス流量）を変更し得るようになっている。
【００２５】
電気制御装置６０は、互いにバスで接続されたＣＰＵ６１、ＣＰＵ６１が実行するプログラム、テーブル（ルックアップテーブル、マップ）、及び定数等を予め記憶したＲＯＭ６２、ＣＰＵ６１が必要に応じてデータを一時的に格納するＲＡＭ６３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ６４、並びにＡＤコンバータを含むインターフェース６５等からなるマイクロコンピュータである。
【００２６】
インターフェース６５は、空気流量（新気流量）計測手段であって吸気管３２に配置された熱線式エアフローメータ７１、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に設けられた吸気温センサ７２、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に配設された吸気管圧力センサ７３、クランクポジションセンサ７４、アクセル開度センサ７５、及び、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に設けられた吸気酸素濃度センサ７６と接続されていて、これらのセンサからの信号をＣＰＵ６１に供給するようになっている。また、インターフェース６５は、燃料噴射弁２１、燃料噴射用ポンプ２２、スロットル弁アクチュエータ３３ａ、及びＥＧＲ制御弁５２と接続されていて、ＣＰＵ６１の指示に応じてこれらに駆動信号を送出するようになっている。
【００２７】
熱線式エアフローメータ７１は、吸気通路内を通過する吸入空気（新気）の質量流量（単位時間あたりの吸入新気量）を計測し、同質量流量Ga（吸入新気流量Ga）を表す信号を発生するようになっている。吸気温センサ７２は、前述した吸気の温度を検出し、同吸気温度Tbを表す信号を発生するようになっている。吸気管圧力センサ７３は、吸気の圧力（即ち、吸気管圧力）を検出し、同吸気管圧力Pbを表す信号を発生するようになっている。
【００２８】
クランクポジションセンサ７４は、各気筒の絶対クランク角度を検出し、クランク角度CAを表すとともにエンジン１０の回転速度であるエンジン回転速度NEをも表す信号を発生するようになっている。アクセル開度センサ７５は、アクセルペダルＡＰの操作量を検出し、アクセル操作量Accpを表す信号を発生するようになっている。吸気酸素濃度センサ７６は、吸気中の酸素濃度（即ち、吸気酸素濃度）を検出し、同吸気酸素濃度RO2_inを表す信号を発生するようになっている。
【００２９】
（ＮＯｘ発生量推定方法の概要）
次に、上記のように構成された内燃機関の制御装置（以下、「本装置」と云うこともある。）による本発明の実施形態に係るＮＯｘ発生量推定方法の概要について説明する。図２は、機関１０の或る一つの気筒のシリンダ内(筒内)に吸気マニホールド３１からガス（即ち、吸気）が吸入され、筒内に吸入されたガス（筒内ガス）が排気マニホールド４１へ排出される様子を模式的に示した図である。
【００３０】
図２に示したように、吸気（従って、筒内ガス）には、吸気管３２の先端部からスロットル弁３３を介して吸入された新気と、排気還流管５１からＥＧＲ制御弁５２を介して吸入されたＮＯｘを含むＥＧＲガスが含まれる。吸入される新気量（新気質量）と吸入されるＥＧＲガス量（ＥＧＲガス質量）の和に対するＥＧＲガス質量の割合（即ち、ＥＧＲ率）は、運転状態に応じて電気制御装置６０（ＣＰＵ６１）により適宜制御されるスロットル弁３３の開度、及びＥＧＲ制御弁５２の開度に応じて変化する。
【００３１】
吸気（即ち、新気、及びＮＯｘを含むＥＧＲガスから構成されるガス）は、吸気行程において開弁している吸気弁Vinを介してピストンの下降に伴って筒内に吸入されて筒内ガスとなる。筒内ガスは、ピストンが下死点に達した時点（以下、「ATDC-180°」と称呼する。）で吸気弁Vinが閉弁することにより筒内に密閉され、その後の圧縮行程においてピストンの上昇に伴って圧縮される。そして、ピストンが上死点近傍に達っすると（具体的には、後述する最終燃料噴射時期finjfinが到来すると）、本装置は、前記指令燃料噴射量qfinに応じた所定時間だけ燃料噴射弁２１を開弁することで燃料を筒内に直接噴射する。この結果、噴射された燃料は、時間の経過に伴って同筒内ガスと混ざり合いながら混合気となって筒内において拡散していき、所定のタイミングにて自己着火が発生することに起因して燃焼する。
【００３２】
係る燃焼は、本実施例では、後述するように推定される燃焼領域（以下、「Ｂ領域」と云うこともある。図２を参照。）においてのみ発生し、燃焼室内におけるＢ領域を除いた領域である非燃焼領域（以下、「Ａ領域」と云うこともある。図２を参照。）では発生しないと仮定する。そして、燃焼後に燃焼室内に存在する筒内ガスは、排ガスとなって、排気行程において開弁している排気弁Voutを介してピストンの上昇に伴って排気マニホールド４１へ排出され、係る排ガスは、排気管４２を介して外部へと排出されていく。
【００３３】
以下、本装置が実施する具体的なＮＯｘ発生量推定方法について説明していく。このＮＯｘ発生量推定方法では、燃料が噴射される対象となる気筒（以下、「燃料噴射気筒」と称呼する。）について上記最終燃料噴射時期finjfinが到来する毎に、その直後の爆発行程において上記Ｂ領域において燃焼により発生するＢ領域燃焼発生ＮＯｘ量NOxBが推定されていく。
【００３４】
Ｂ領域燃焼発生ＮＯｘ量NOxBは、単位燃料量あたりの燃焼発生ＮＯｘ量（以下、「燃焼発生ＮＯｘ率RNOx_burn」と称呼する。）に今回の指令燃料噴射量qfinc(＝qfin)を乗じた値として、下記(1)式に従って求めることができる。
【００３５】
NOxB＝RNOx_burn・qfinc  ・・・(1)
【００３６】
ここで、上記(1)式における燃焼発生ＮＯｘ率RNOx_burnは下記(2)式にて推定される。
【００３７】
RNOx_burn＝eK0・(RO2c)K1・(qfinc)K2・(Pcrc)K3・e(K4/Tflame)  ・・・(2)
【００３８】
上記(2)式において、eは自然対数の底である。RO2cは、吸気弁Vin閉弁時（即ち、ATDC-180°）において吸気酸素濃度センサ７６により検出される吸気酸素濃度RO2_inである下死点時吸気酸素濃度である。qfincは上記(1)式と同様、今回の指令燃料噴射量(=qfin)である。Pcrcは今回の指令燃料噴射圧力(=Pcrbase)である。
【００３９】
また、上記(2)式において、Tflameは今回の爆発行程における最高火炎温度である。最高火炎温度Tflameは、混合気が燃焼を開始してから同燃焼を終了するまでの間における火炎温度のピーク値であって、エンジン回転速度NE、及び今回の指令燃料噴射量qfincを引数とする所定の関数に基づいて推定することができる。K0〜K4は代表的な公知の所定の重回帰分析に基づいて後述するように決定された適合定数である。
【００４０】
即ち、上記(2)式は、燃焼発生ＮＯｘ率RNOx_burnを求めるための実験式である。上記(2)式にて推定される燃焼発生ＮＯｘ率RNOx_burnは、上記下死点時吸気酸素濃度RO2c、上記今回の指令燃料噴射量qfinc、上記今回の指令燃料噴射圧力Pcrc、及び上記最高火炎温度Tflameの関数であって、より具体的には、下死点時吸気酸素濃度RO2cのべき乗、今回の指令燃料噴射量qfincのべき乗、今回の指令燃料噴射圧力Pcrcのべき乗、及び、最高火炎温度Tflameに応じて決定される値を指数とする指数関数、の積に基づいて計算される。
【００４１】
適合定数K0〜K4は、例えば、以下のような実験を行うことにより決定され得る。即ち、先ず、ＥＧＲ制御弁５２を閉状態に維持して機関１０を運転させることで、排気弁Voutから排出された排ガス（従って、同排ガス中のＮＯｘ）の総てが排気通路から外部へ排出されるようにしておく。これにより、排気通路から外部へ排出される排ガス中のＮＯｘ量（即ち、上記ＮＯｘ排出量）は上記Ｂ領域燃焼発生ＮＯｘ量NOxBと等しくなるから、ＮＯｘ排出量を所定のＮＯｘ濃度センサの出力に基づいて計測することでＢ領域燃焼発生ＮＯｘ量NOxB（従って、燃焼発生ＮＯｘ率RNOx_burn（＝NOxB/qfinc））を計測できるようになる。
【００４２】
次に、この状態にて、上記下死点時吸気酸素濃度RO2c、上記今回の指令燃料噴射量qfinc、上記今回の指令燃料噴射圧力Pcrc、及び上記最高火炎温度Tflame（即ち、エンジン回転速度NEと今回の指令燃料噴射量qfinc）の各値の組み合わせを種々の所定パターンになるように順次変更していき、そのパターン毎に、燃焼発生ＮＯｘ率RNOx_burnを順次計測していく。
【００４３】
そして、このような作業（実験）の結果得られた上記各値の組み合わせと、計測された燃焼発生ＮＯｘ率RNOx_burnの値との関係に関する多数のデータに基づいて、上記公知の所定の重回帰分析を行って上記適合定数K0〜K4を求めることができる。ここにおいて、少なくとも、適合定数K1〜K3は正の値に決定され、適合定数K4は負の値に決定される。
【００４４】
従って、上記(2)式から理解できるように、(2)式に従って推定・計算される燃焼発生ＮＯｘ率RNOx_burn（従って、Ｂ領域燃焼発生ＮＯｘ量NOxB）は、下死点時吸気酸素濃度RO2c、今回の指令燃料噴射量qfinc、今回の指令燃料噴射圧力Pcrc、及び最高火炎温度Tflameのうちのどの値が増加しても増大する。このことは、以下の実際の現象に沿うものである。
【００４５】
即ち、先ず、吸気酸素濃度RO2_inが増加するとＢ領域燃焼発生ＮＯｘ量NOxBが増大する。これは、酸素はＮＯｘを発生させる材料であって燃焼室内の酸素量が多くなると当然にＮＯｘが発生し易くなることに基づく現象である。
【００４６】
また、燃料噴射量qfinが増加するとＢ領域燃焼発生ＮＯｘ量NOxBが増大する。これは、燃料噴射量qfinの増加に応じて機関１０に与えられる負荷の大きさが増大することから、上述した燃焼室の内壁面温度の上昇に起因して、燃料噴射量qfinが増大するほど（即ち、機関に与えられる負荷が増加するほど）ＮＯｘが発生し易くなることに基づく現象である。
【００４７】
また、燃料噴射圧力Pcrが増加するとＢ領域燃焼発生ＮＯｘ量NOxBが増大する。これは、燃料噴射圧力Pcrの増加に応じて燃料の噴射速度が増大すること等により同燃料の霧化の程度が大きくなるから、上述した空気過剰率の増大に起因して、燃料噴射圧力Pcrが増大するほど（即ち、噴射された燃料の霧化の程度が大きくなるほど）ＮＯｘが発生し易くなることに基づく現象である。
【００４８】
また、最高火炎温度Tflameが増加するとＢ領域燃焼発生ＮＯｘ量NOxBが増大する。これは、ガス温度が高くなるほど窒素からＮＯｘが生成される化学反応が促進されることに基づく現象である。以上のことから、上記(2)式に従って燃焼発生ＮＯｘ率RNOx_burnを計算すれば、少なくとも上記４つの実際の現象に沿うように同燃焼発生ＮＯｘ率RNOx_burn（従って、上記(1)式に従うＢ領域燃焼発生ＮＯｘ量NOxB）を精度良く推定することができる。以上が、ＮＯｘ発生量推定方法の概要である。
【００４９】
（ＮＯｘ排出量推定方法の概要）
次に、図２に示したようなＥＧＲ装置５０を備えた内燃機関１０に適用される、本装置が実施する具体的なＮＯｘ排出量推定方法について説明していく。このＮＯｘ排出量推定方法では、燃料噴射気筒について上記最終燃料噴射時期finjfinが到来する毎に、その直後の排気行程において同燃料噴射気筒の排気弁Voutから排気通路を介して外部へ排出されることになる排ガス中のＮＯｘの質量（即ち、ＮＯｘ排出量。実ＮＯｘ排出量NOxact。）が推定されていく。
【００５０】
この方法では、実ＮＯｘ排出量NOxactを推定するにあたり、燃焼室内に吸入されたＮＯｘの全質量に対する上記燃焼前に上記Ｂ領域に存在するＮＯｘの質量の割合（以下、「ＮＯｘ量割合RatioNOx」と称呼する。）、上記燃焼後に上記Ａ領域内に残存するＮＯｘ質量、及び、上記燃焼後に上記Ｂ領域内に残存するＮＯｘ質量をそれぞれ推定する必要がある。従って、先ず、これらの求め方について図２を参照しながら説明する。
【００５１】
<ＮＯｘ量割合RatioNOxの求め方>
燃焼室内に吸入された吸気（筒内ガス）中の酸素分子、ＮＯｘを含む各ガスの成分は同燃焼室内全域に渡ってそれぞれ均一に分布するものと仮定する。また、この状態においてＢ領域内に存在する酸素は総て燃焼により消費されると仮定する。そうすると、「燃焼室内に吸入された酸素の全質量」に対する「燃焼により消費される酸素の質量」の割合（以下、「酸素量割合」と称呼する。）は、燃焼室の体積に対するＢ領域の体積の割合を表すとともに、燃焼室内に吸入されたＮＯｘの全質量に対する、燃焼前にＢ領域に存在するＮＯｘの質量の割合（従って、上記ＮＯｘ量割合RatioNOx）をも表すことになる。換言すれば、上記酸素量割合を利用してＢ領域が推定され得る。
【００５２】
従って、ＮＯｘ量割合RatioNOxを求めるためには上記酸素量割合を求めればよく、このためには、「燃焼室内に吸入された酸素の全質量」と「燃焼により消費される酸素の質量」とを求める必要がある。
【００５３】
ここで、「燃焼室内に吸入された酸素の全質量」は、吸入された燃焼室内の総ガス質量（以下、「筒内総ガス量Gcyl」と称呼する。）に、燃焼前の筒内ガスの酸素濃度を乗じることで求めることができる。筒内総ガス量Gcylは、ATDC-180°における気体の状態方程式に基づく下記(3)式に従って求めることができる。
【００５４】
Gcyl＝(Pa0・Va0)/(R・Ta0)  ・・・(3)
【００５５】
上記(3)式において、Pa0は、ATDC-180°における下死点時筒内ガス圧力である。ATDC-180°において筒内ガス圧力は吸気管圧力Pbと略等しいと考えられるから、下死点時筒内ガス圧力Pa0は、ATDC-180°において吸気管圧力センサ７３により検出される吸気管圧力Pbとして取得することができる。Va0は、ATDC-180°における下死点時燃焼室内容積である。燃焼室内容積Vaは機関１０の設計諸元に基づいてクランク角度CAの関数として表すことができるから、この関数に基づいて下死点時燃焼室内容積Va0も求めることができる。Ta0は、ATDC-180°における下死点時筒内ガス温度である。ATDC-180°において筒内ガス温度は吸気温度Tbと略等しいと考えられるから、下死点時筒内ガス温度Ta0は、ATDC-180°において吸気温センサ７２により検出される吸気温度Tbとして取得することができる。Rは筒内ガスのガス定数である。
【００５６】
また、燃焼前の筒内ガスの酸素濃度は吸気弁Vin閉弁時（即ち、ATDC-180°）における吸気酸素濃度RO2_inと略等しいと考えられるから、燃焼前の筒内ガスの酸素濃度は、上記(2)式にて使用される下死点時吸気酸素濃度RO2cとして取得することができる。以上より、「燃焼室内に吸入された酸素の全質量」は「Gcyl・RO2c」と表すことができる。
【００５７】
一方、「燃焼により消費される酸素の質量」は、噴射された燃料の総て（即ち、前記燃料噴射量qfinの燃料）が理論空燃比stoichをもって完全燃焼するとの仮定のもと、「K・qfinc」と表すことができる。ここで、Kは所定の係数であって、例えば、大気中に含まれる酸素の質量割合０．２３に理論空燃比stoich（例えば、１４．６）を乗じた値「０．２３・stoich」である。qfincは上記(1)式,(2)式と同様、今回の指令燃料噴射量(＝qfin)である。以上のことから、ＮＯｘ量割合RatioNOxは、下記(4)式に従って求めることができる。
【００５８】
RatioNOx＝(K・qfinc)/(Gcyl・RO2c)  ・・・(4)
【００５９】
<燃焼後にＡ領域内に残存するＮＯｘ質量の求め方>
上述のごとく、本実施例では燃焼室内における燃焼はＢ領域内でのみ発生し、Ａ領域では燃焼が発生しないと仮定するから、ＥＧＲ装置５０を介して燃焼室内に還流してきたＮＯｘ量のうち燃焼前にＡ領域内に存在するＮＯｘの質量（以下、「Ａ領域還流ＮＯｘ量NOxA」と称呼する。）は、燃焼後においてもそのままＡ領域内に保存（保持）されると考えることができる。換言すれば、Ａ領域還流ＮＯｘ量NOxAは、そのまま「燃焼後にＡ領域内に残存するＮＯｘ質量」となる。
【００６０】
ここで、燃焼室内に吸入されたＮＯｘの全質量に対する上記燃焼前に上記Ａ領域に存在するＮＯｘの質量の割合は、上記ＮＯｘ量割合RatioNOxを用いて「1-RatioNOx」と表すことができるから、Ａ領域還流ＮＯｘ量NOxAは、上記筒内総ガス量Gcylに、燃焼前の筒内ガスのＮＯｘ濃度と（1-RatioNOx）とを乗じることで求めることができる。燃焼前の筒内ガスのＮＯｘ濃度は吸気中のＮＯｘ濃度（吸気ＮＯｘ濃度RNOx_in）と略等しいと考えられるから、Ａ領域還流ＮＯｘ量NOxAは、下記(5)式にて表すことができる。
【００６１】
NOxA＝RNOx_in・(1-RatioNOx)・Gcyl  ・・・(5)
【００６２】
上記(5)式において、吸気ＮＯｘ濃度RNOx_inは、上記筒内総ガス量Gcylに対する、ＥＧＲ装置５０から還流されてきたＥＧＲガス中のＮＯｘ質量の質量割合である。ＥＧＲガス中のＮＯｘ濃度は前回（の燃料噴射時期に）演算された後述する排気ＮＯｘ濃度RNOx_exと等しいと仮定すると、吸気ＮＯｘ濃度RNOx_inは下記(6)式に従って求めることができる。
【００６３】
RNOx_in＝(RNOx_ex・Gegr)/Gcyl  ・・・(6)
【００６４】
上記(6)式において、Gegrは今回の吸気行程で吸気の一部としてＥＧＲ装置５０から燃焼室内に吸入されてきたＥＧＲガスの質量であって、下記(7)式に従って取得され得る。
【００６５】
Gegr＝Gcyl-Gm  ・・・(7)
【００６６】
上記(7)式において、Gmは今回の吸気行程で吸気管３２の先端部から吸気の一部として燃焼室内に吸入された新気の質量（吸入新気量）であって、エアフローメータ７１により計測される単位時間あたりの吸入新気量（吸入新気流量Ga）と、クランクポジションセンサ７４の出力に基づくエンジン回転速度NEと、吸入新気流量Ga及びエンジン回転速度NEを引数とする一吸気行程あたりの吸入新気量を求めるための関数ｆ(Ga,NE)とに基づいて計算される。吸入新気流量Ga及びエンジン回転速度NEとしては、ATDC-180°において各センサにより取得される下死点時吸入新気流量Ga0、及び下死点時エンジン回転速度NE0がそれぞれ使用される。以上のように、Ａ領域還流ＮＯｘ量NOxA、従って、「燃焼後にＡ領域内に残存するＮＯｘ質量」は上記(5)式に従って求めることができる。
【００６７】
<燃焼後にＢ領域内に残存するＮＯｘ質量の求め方>
燃焼後においてはＢ領域において、上記(1)式、(2)式により推定される上記Ｂ領域燃焼発生ＮＯｘ量NOxBが発生する。そして、この推定されたＢ領域燃焼発生ＮＯｘ量NOxBは、そのまま「燃焼後にＢ領域内に残存するＮＯｘ質量」に略一致し得ると考えることができる。従って、「燃焼後にＢ領域内に残存するＮＯｘ質量」は、上記(1)式、(2)式に従って推定されるＢ領域燃焼発生ＮＯｘ量NOxBとして求めることができる。
【００６８】
このようにして、「燃焼後にＡ領域内に残存するＮＯｘ質量」であるＡ領域還流ＮＯｘ量NOxA、及び、「燃焼後にＢ領域内に残存するＮＯｘ質量」であるＢ領域燃焼発生ＮＯｘ量NOxBが求められると、燃焼後に燃焼室内に残存するＮＯｘの全質量は「NOxA＋NOxB」として求めることができる。また、燃焼後に燃焼室内に残存するガスの全質量は、「Gcyl＋qfinc」と表すことができる。
【００６９】
従って、排気行程において燃焼室から排気弁Voutを介して排気通路（排気マニホールド４１）に排出される排ガス中のＮＯｘ濃度（排気ＮＯｘ濃度RNOx_ex）は、「燃焼後に燃焼室内に残存するガスの全質量」に対する、「燃焼後に燃焼室内に残存するＮＯｘの全質量」の質量割合であって、下記(8)式に従って求めることができる。
【００７０】
RNOx_ex＝(NOxA＋NOxB)/(Gcyl＋qfinc)  ・・・(8)
【００７１】
上記(8)式に従って求められる排気ＮＯｘ濃度RNOx_exの前回値は、排気通路（排気マニホールド４１）を流れる排ガス中のＮＯｘ濃度が排気還流管５１内を流れるＥＧＲガス中のＮＯｘ濃度に等しいとの仮定のもと、前述のごとく、上記吸気ＮＯｘ濃度RNOx_in（の
今回値）を求める際に上記(6)式にて使用される。
【００７２】
また、排気通路（排気マニホールド４１、及び排気管４２）を流れる排ガス中のＮＯｘ濃度は同排気通路全域に渡って一定であると仮定すると、上記排気ＮＯｘ濃度RNOx_exの値は、排気通路（具体的には、排気管４２の端部）から外部に排出される際の排ガス中のＮＯｘ濃度と等しくなる。
【００７３】
更には、通常の機関１０の運転状態（特に、定常運転状態）においては、排気通路（排気管４２）から外部へ排出される一作動サイクル（一排気行程）あたりの排ガスの質量は上記吸入新気量Gmと略等しくなる。以上のことから、排気通路を介して外部へ排出される排ガス中の一作動サイクルあたりのＮＯｘの質量（上記実ＮＯｘ排出量NOxact）は、下記(9)式に従って求めることができる。これにより、ＥＧＲガス量Gegrが増加すると吸入新気量Gmが減少することで実ＮＯｘ排出量NOxactが低下することになる。従って、ＥＧＲガス量Gegrが増加すると実ＮＯｘ排出量NOxactが低下するという現象が的確に表され得る。
【００７４】
NOxact＝RNOx_ex・Gm  ・・・(9)
【００７５】
以上、(1)〜(9)式を利用して、本装置は、燃料噴射気筒についての上記最終燃料噴射時期finjfinが到来する毎に、その直後の排気行程において同燃料噴射気筒の排気弁Voutから排出されることになるＮＯｘの質量、即ち実ＮＯｘ排出量NOxactを推定する。以上が、ＮＯｘ排出量推定方法の概要である。
【００７６】
（燃料噴射制御の概要）
上記ＮＯｘ排出量推定方法を実施する本装置は、一作動サイクルあたりの目標ＮＯｘ排出量NOxtを、上記燃料噴射量qfinとエンジン回転速度NEとに基づいて逐次計算する。そして、本装置は、前回推定された実ＮＯｘ排出量NOxactが目標ＮＯｘ排出量NOxtに一致するように、最終燃料噴射開始時期finjfin、及びＥＧＲ制御弁５２の開度をフィードバック制御する。
【００７７】
具体的には、前回推定された上記実ＮＯｘ排出量NOxactの値が上記目標ＮＯｘ排出量NOxtよりも大きいとき、今回の燃料噴射気筒についての最終燃料噴射開始時期finjfinを基本燃料噴射時期finjbaseよりも所定量だけ遅らせ、且つ、ＥＧＲ制御弁５２の開度を現時点での値から所定開度だけ大きくする。これにより、今回の燃料噴射気筒についての最高火炎温度が低くなるように制御され、この結果、今回の燃料噴射気筒から外部へ排出される実ＮＯｘ排出量NOxactが上記目標ＮＯｘ排出量NOxtに一致せしめられる。
【００７８】
一方、前回推定された上記実ＮＯｘ排出量NOxactの値が上記目標ＮＯｘ排出量NOxtよりも小さいとき、今回の燃料噴射気筒についての最終燃料噴射開始時期finjfinを基本燃料噴射時期finjbaseよりも所定量だけ早め、且つ、ＥＧＲ制御弁５２の開度を現時点での値から所定開度だけ小さくする。これにより、今回の燃料噴射気筒についての最高火炎温度が高くなるように制御され、この結果、今回の燃料噴射気筒から外部へ排出される実ＮＯｘ排出量NOxactが上記目標ＮＯｘ排出量NOxtに一致せしめられる。以上が、燃料噴射制御の概要である。
【００７９】
（燃焼発生ＮＯｘ率RNOx_burnの実際の演算方法）
上記(2)式に従って燃焼発生ＮＯｘ率RNOx_burnを演算するためには、「べき乗」の演算、及び「掛け算」の演算が必要となる。しかしながら、一般に、マイクロコンピュータを使用して「べき乗」の演算を実行すると計算負荷が増大し、同マイクロコンピュータを使用して「掛け算」の演算を実行すると計算精度が低下する傾向がある。よって、「べき乗」の演算、及び「掛け算」の演算を避けるため、本装置（ＣＰＵ６１）は、実際には、上記(2)式において両辺の自然対数をとることで得られる下記(10)式を利用して、テーブル検索と「足し算」の演算のみで燃焼発生ＮＯｘ率RNOx_burnを計算する。
【００８０】
log(RNOx_burn)＝K0＋K1・log(RO2c)＋K2・log(qfinc)＋K3・log(Pcrc)＋K4/Tflame  ・・・(10)
【００８１】
即ち、本装置は、上記(10)式の右辺の第２項〜第５項までの各項を求めるために予めＲＯＭ６２にそれぞれ記憶されているテーブルMaplog1(RO2c)、Maplog2(qfinc)、Maplog3(Pcrc)、及びMapinvpro(Tflame)に基づいて、テーブル検索値dataMap1(＝K1・log(RO2c))、dataMap2(＝K2・log(qfinc))、dataMap3(＝K3・log(Pcrc))、及びdataMap4(＝K4/Tflame)をそれぞれ決定し、「足し算」を伴う下記(11)式に従って「log(RNOx_burn)」を求める。
【００８２】
log(RNOx_burn)＝K0＋dataMap1＋dataMap2＋dataMap3＋dataMap4  ・・・(11)
【００８３】
そして、本装置は、上記（11）式に従って求められた「log(RNOx_burn)」から燃焼発生ＮＯｘ率RNOx_burnを求めるために予めＲＯＭ６２に記憶されているテーブルMapinvlog(log(RNOx_burn))に基づいて同燃料発生ＮＯｘ率RNOx_burnを求める。これにより、ＣＰＵ６１の計算負荷を軽減できるとともに計算精度の低下を防止できる。
【００８４】
（実際の作動）
次に、上記のように構成された内燃機関の制御装置の実際の作動について説明する。
<燃料噴射量等の制御>
ＣＰＵ６１は、図３にフローチャートにより示した燃料噴射量等の制御を行うためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ３００から処理を開始し、ステップ３０５に進んでアクセル開度Accp、エンジン回転速度NE、及び図４に示したテーブル（マップ）Mapqfinから（指令）燃料噴射量qfinを求める。テーブルMapqfinは、アクセル開度Accp及びエンジン回転速度NEと燃料噴射量qfinとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００８５】
次いで、ＣＰＵ６１はステップ３１０に進み、燃料噴射量qfin、エンジン回転速度NE、及び図５に示したテーブルMapfinjbaseから基本燃料噴射時期finjbaseを決定する。テーブルMapfinjbaseは、燃料噴射量qfin及びエンジン回転速度NEと基本燃料噴射時期finjbaseとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００８６】
その後、ＣＰＵ６１はステップ３１５に進んで、燃料噴射量qfin、エンジン回転速度NE、及び図６に示したテーブルMapPcrbaseから基本燃料噴射圧力Pcrbaseを決定する。テーブルMapPcrbaseは、燃料噴射量qfin及びエンジン回転速度NEと基本燃料噴射圧力Pcrbaseとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００８７】
次に、ＣＰＵ６１はステップ３２０に進み、燃料噴射量qfin、エンジン回転速度NE、及び図７に示したテーブルMapNOxtから目標ＮＯｘ排出量NOxtを決定する。テーブルMapNOxtは、燃料噴射量qfin及びエンジン回転速度NEと目標ＮＯｘ排出量NOxtとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００８８】
次いで、ＣＰＵ６１はステップ３２５に進んで、前記決定した目標ＮＯｘ排出量NOxtから後述するルーチンにより求められている最新の（具体的には、前回の燃料噴射時期に演算された）実ＮＯｘ排出量NOxactを減じた値をＮＯｘ排出量偏差ΔNOxとして格納する。
【００８９】
続いて、ＣＰＵ６１はステップ３３０に進んで、噴射時期補正値Δθを、前記ＮＯｘ排出量偏差ΔNOxと図８に示したテーブルMapΔθとから決定する。テーブルMapΔθは、ＮＯｘ排出量偏差ΔNOxと噴射時期補正値Δθとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００９０】
次いで、ＣＰＵ６１はステップ３３５に進み、基本噴射時期finjbaseを噴射時期補正値Δθで補正して最終燃料噴射時期finjfinを決定する。これにより、ＮＯｘ排出量偏差ΔNOxに応じて噴射時期が補正されることになる。この場合、図８から明らかなように、ＮＯｘ排出量偏差ΔNOxが正の大きい値になるほど噴射時期補正値Δθが正の大きな値となって最終燃料噴射時期finjfinが進角側となり、同ＮＯｘ排出量偏差ΔNOxが負の大きい値（絶対値が大きい値）になるほど噴射時期補正値Δθは負の大きな値となって最終燃料噴射時期finjfinが遅角側に移行される。
【００９１】
続いて、ＣＰＵ６１はステップ３４０に進み、燃料噴射気筒についての噴射開始時期（即ち、上記決定された最終燃料噴射時期finjfin）が到来したか否かを判定し、「Ｎｏ」と判定する場合はステップ３９５に直ちに進んで本ルーチンを一旦終了する。
【００９２】
一方、ステップ３４０の判定において「Ｙｅｓ」と判定する場合、ＣＰＵ６１はステップ３４５に進んで上記決定された（指令）燃料噴射量qfinの燃料を燃料噴射気筒についての燃料噴射弁２１から上記決定された基本燃料噴射圧力Pcrbaseをもって噴射するとともに、続くステップ３５０にて上記ＮＯｘ排出量偏差ΔNOxが正の値であるか否かを判定し、「Ｙｅｓ」と判定する場合、ステップ３５５に進んでＥＧＲ制御弁５２の開度を現時点での値よりも所定開度だけ小さくした後、ステップ３７０に進む。
【００９３】
ステップ３５０の判定において「Ｎｏ」と判定する場合、ＣＰＵ６１はステップ３６０に進み、上記ＮＯｘ排出量偏差ΔNOxが負の値であるか否かを判定する。ステップ３６０の判定において、ＣＰＵ６１は「Ｙｅｓ」と判定する場合、ＥＧＲ制御弁５２の開度を現時点での値よりも所定開度だけ大きくした後にステップ３７０に進む一方、「Ｎｏ」と判定する場合（即ち、ＮＯｘ排出量偏差ΔNOxの値が「０」のとき）、ＥＧＲ制御弁５２の開度を変更することなくステップ３７０に進む。
【００９４】
このようにして、ＮＯｘ排出量偏差ΔNOxに応じてＥＧＲ制御弁５２の開度が変更される。そして、ステップ３７０に進むと、ＣＰＵ６１は実際に噴射した上記燃料噴射量qfinの値を今回の燃料噴射量qfincとして格納し、続くステップ３７５にて実際の噴射圧力である上記基本燃料噴射圧力Pcrbaseの値を今回の燃料噴射圧力Pcrcとして格納した後、ステップ３９５に進んで本ルーチンを一旦終了する。以上により、燃料噴射量、燃料噴射時期、燃料噴射圧力、及びＥＧＲ制御弁５２の開度の制御が達成される。
【００９５】
<ＮＯｘ排出量の計算>
また、ＣＰＵ６１は、図９にフローチャートにより示した実ＮＯｘ排出量NOxactの計算を行うためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ９００から処理を開始し、ステップ９０５に進んで現時点がATDC-180°に一致しているか否かを判定する。
【００９６】
いま、現時点がATDC-180°より前であるものとして説明を続けると、ＣＰＵ６１はステップ９０５にて「Ｎｏ」と判定してステップ９３５に直ちに進み、燃料噴射気筒についての燃料噴射開始時期（即ち、上記最終燃料噴射時期finjfin）が到来したか否かを判定する。現時点はATDC-180°よりも前であるから、ＣＰＵ６１はステップ９３５にて「Ｎｏ」と判定し、ステップ９９５に直ちに進んで本ルーチンを一旦終了する。
【００９７】
以降、ＣＰＵ６１はATDC-180°が到来するまでの間、ステップ９００、９０５、９３５、９９５の処理を繰り返し実行する。そして、ATDC-180°が到来すると、ＣＰＵ６１はステップ９０５に進んだとき「Ｙｅｓ」と判定してステップ９１０に進むようになり、ステップ９１０にて、吸気温センサ７２、吸気管圧力センサ７３、エアフローメータ７１、及びクランクポジションセンサ７４によりそれぞれ検出される現時点（即ち、ATDC-180°）での吸気温度Tb、吸気管圧力Pb、吸入新気流量Ga、及びエンジン回転速度NEを、それぞれ、下死点時筒内ガス温度Ta0、下死点時筒内ガス圧力Pa0、下死点時吸入新気流量Ga0、及び下死点時エンジン回転速度NE0として格納する。
【００９８】
次いで、ＣＰＵ６１はステップ９１５に進み、吸気酸素濃度センサ７６により検出される現時点（即ち、ATDC-180°）での吸気酸素濃度RO2_inを下死点時吸気酸素濃度RO2cとして格納し、続くステップ９２０にて、上記(1)式に従って筒内総ガス量Gcylを求める。ここで、下死点時筒内ガス圧力Pa0、及び下死点時筒内ガス温度Ta0としては、ステップ９１０にて格納されている値が用いられる。
【００９９】
続いて、ＣＰＵ６１はステップ９２５に進んで、上記下死点時吸入新気流量Ga0と、上記下死点時エンジン回転速度NE0と、上記関数ｆとに基づいて吸入新気量Gmを求め、続くステップ９３０にて、ステップ９２０にて求めた筒内総ガス量Gcylと、前記吸入新気量Gmと、上記(5)式とに基づいてＥＧＲガス量Gegrを求める。そして、ＣＰＵ６１はステップ９３５に進んで「Ｎｏ」と判定してステップ９９５に直ちに進んで本ルーチンを一旦終了する。
【０１００】
以降、ＣＰＵ６１は燃料噴射時期（即ち、上記最終燃料噴射時期finjfin）が到来するまでの間、ステップ９００、９０５、９３５、９９５の処理を再び繰り返し実行する。そして、上記最終燃料噴射時期finjfinが到来すると、ＣＰＵ６１はステップ９３５にて「Ｙｅｓ」と判定してステップ９４０に進むようになり、ステップ９４０にて上記(6)式に従って吸気ＮＯｘ濃度RNOx_inを求める。ここで、ＥＧＲガス量Gegr、及び筒内総ガス量Gcylとして、ステップ９３０、及びステップ９２０にて求められた値がそれぞれ使用される。排気ＮＯｘ濃度RNOx_exとしては、前回の燃料噴射開始時期において後述するステップ９６５にて求められている値が使用される。
【０１０１】
続いて、ＣＰＵ６１はステップ９４５に進んで、上記(4)式、及び上記(5)式に相当するステップ９４５内に記載の式に基づいてＡ領域還流ＮＯｘ量NOxAを求める。ここで、今回の燃料噴射量qfincとしては、図３のステップ３７０にて格納されている最新の値が使用される。次に、ＣＰＵ６１はステップ９５０を経由して図１０のステップ１０００に進み、燃焼発生ＮＯｘ率RNOx_burnの計算を行うための処理を開始する。
【０１０２】
即ち、ＣＰＵ６１はステップ１０００からステップ１００５に進むと、先ず、現時点でのエンジン回転速度NE、前記今回の燃料噴射量qfincと、及び最高火炎温度Tflameを求めるためのステップ１００５内に記載されたテーブルと、に基づいて最高火炎温度Tflameを推定・決定する。
【０１０３】
次に、ＣＰＵ６１はステップ１０１０に進んで図９のステップ９１５にて求められている下死点時吸気酸素濃度RO2cの最新値と前記テーブルMaplog1とに基づいて前記テーブル検索値dataMap1(＝K1・log(RO2c))を求める。
【０１０４】
同様に、ＣＰＵ６１はステップ１０１５に進んで図３のステップ３７０にて格納されている今回の燃料噴射量qfincと前記テーブルMaplog2とに基づいて前記テーブル検索値dataMap2(＝K2・log(qfinc))を求め、続くステップ１０２０にて図３のステップ３７５にて格納されている今回の燃料噴射圧力Pcrcと前記テーブルMaplog3とに基づいて前記テーブル検索値dataMap3(＝K3・log(Pcrc))を求め、続くステップ１０２５にてステップ１００５にて求めた最新の最高火炎温度Tflameと前記テーブルMapinvproとに基づいて前記テーブル検索値dataMap4(＝K4/Tflame)を求める。
【０１０５】
次いで、ＣＰＵ６１はステップ１０３０に進んで、上記(11)式に従って「log(RNOx_burn)」を求め、続くステップ１０３５にて、同log(RNOx_burn)と前記テーブルMapinvlogとに基づいて燃焼発生ＮＯｘ率RNOx_burnを求めた後、ステップ１０９５を経由して、図９のステップ９５５に進む。
【０１０６】
ＣＰＵ６１はステップ９５５に進むと、上記(1)式に従ってＢ領域燃焼発生ＮＯｘ量NOxBを求める。次いで、ＣＰＵ６１はステップ９６０に進んで、上記(8)式に従って排気ＮＯｘ濃度RNOx_exを求めるとともに、続くステップ９６５にて上記(9)式に従って実ＮＯｘ排出量NOxactを求め、ステップ９９５に進んで本ルーチンを一旦終了する。以降、ＣＰＵ６１は次の燃料噴射気筒についてのATDC-180°が到来するまでの間、ステップ９００、９０５、９３５、９９５の処理を繰り返し実行する。
【０１０７】
以上のようにして、新たな実ＮＯｘ排出量NOxactが燃料噴射開始時期が到来する毎に求められていく。そして、この新たな実ＮＯｘ排出量NOxactは、前述のごとく、図３のステップ３２５にて使用され、この結果、次の燃料噴射気筒についての最終燃料噴射開始時期finjfin、及びＥＧＲ制御弁５２の開度が同新たな実ＮＯｘ排出量NOxactに基づいてフィードバック制御されていく。
【０１０８】
以上、説明したように、本発明の実施形態に係る内燃機関のＮＯｘ発生量推定方法によれば、燃焼領域（Ｂ領域）において燃焼により発生するＮＯｘ発生量（Ｂ領域燃焼発生ＮＯｘ量NOxB）に大きく影響を与える４つの周辺状態量、即ち、ＮＯｘを発生させる材料となる吸気中のガスの濃度（吸気酸素濃度RO2c）と、機関に与えられる負荷の程度を表す負荷指標値（燃料噴射量qfinc）と、燃焼室内の燃料の霧化の程度を表す霧化指標値（燃料噴射圧力Pcrc）と、及び最高火炎温度Tflameと、に基づいて同Ｂ領域燃焼発生ＮＯｘ量NOxBを推定する。従って、上記４つの周辺状態量とＢ領域燃焼発生ＮＯｘ量NOxBとの間における実際の関係に沿うように、Ｂ領域燃焼発生ＮＯｘ量NOxBを精度良く推定することができた。
【０１０９】
本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記実施形態においては、機関に与えられる負荷の程度を表す負荷指標値として燃料噴射量qfincを採用しているが、同負荷指標値として所謂機関の出力トルクを採用してもよい。また、負荷指標値として、燃焼室の内壁面の温度を採用してもよい。
【０１１０】
また、上記実施形態においては、燃焼室内の霧化の程度を表す霧化指標値として燃料噴射圧力Pcrcを採用しているが、同霧化指標値として所謂スワール比を採用してもよい。また、霧化指標値として燃焼領域（Ｂ領域）における空気過剰率を採用してもよい。
【０１１１】
また、上記実施形態においては、エンジン回転速度NEと燃料噴射量qfincとに基づいて最高火炎温度Tflameを推定しているが、エンジン回転速度NEと、機関の出力トルクとに基づいて最高火炎温度Tflameを推定してもよい。
【図面の簡単な説明】
【０１１２】
【図１】本発明の実施形態に係る内燃機関のＮＯｘ排出量推定方法を実施する内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）に適用したシステム全体の概略構成図である。
【図２】或る一つの気筒のシリンダ内(筒内)に吸気マニホールドからガスが吸入され、筒内に吸入された筒内ガスが排気マニホールドへ排出される様子を模式的に示した図である。
【図３】図１に示したＣＰＵが実行する燃料噴射量等を制御するためのルーチンを示したフローチャートである。
【図４】図１に示したＣＰＵが図３に示したルーチンを実行する際に参照する燃料噴射量を決定するためのテーブルである。
【図５】図１に示したＣＰＵが図３に示したルーチンを実行する際に参照する基本燃料噴射時期を決定するためのテーブルである。
【図６】図１に示したＣＰＵが図３に示したルーチンを実行する際に参照する基本燃料噴射圧力を決定するためのテーブルである。
【図７】図１に示したＣＰＵが図３に示したルーチンを実行する際に参照する目標ＮＯｘ排出量を決定するためのテーブルである。
【図８】図１に示したＣＰＵが図３に示したルーチンを実行する際に参照する噴射時期補正値を決定するためのテーブルである。
【図９】図１に示したＣＰＵが実行するＮＯｘ排出量（実ＮＯｘ排出量）の計算を実行するためのルーチンを示したフローチャートである。
【図１０】図１に示したＣＰＵが実行する燃焼発生ＮＯｘ率の計算を実行するためのルーチンを示したフローチャートである。
【符号の説明】
【０１１３】
２１…燃料噴射弁、２２…燃料噴射用ポンプ、３１…吸気マニホールド、３２…吸気管、４１…排気マニホールド、４２…排気管、５０…ＥＧＲ装置、５２…ＥＧＲ制御弁、６０…電気制御装置、６１…ＣＰＵ、７１…エアフローメータ、７２…吸気温センサ、７３…吸気管圧力センサ、７４…クランクポジションセンサ、７６…吸気酸素濃度センサ
【図１】
2005139984.tif 000003
【図２】
2005139984.tif 000004
【図３】
2005139984.tif 000005
【図４】
2005139984.tif 000006
【図５】
2005139984.tif 000007
【図６】
2005139984.tif 000008
【図７】
2005139984.tif 000009
【図８】
2005139984.tif 000010
【図９】
2005139984.tif 000011
【図１０】
2005139984.tif 000012