JP4123811B2 - Particulate emission estimation device for diesel engine - Google Patents

Particulate emission estimation device for diesel engine Download PDF

Info

Publication number
JP4123811B2
JP4123811B2 JP2002111925A JP2002111925A JP4123811B2 JP 4123811 B2 JP4123811 B2 JP 4123811B2 JP 2002111925 A JP2002111925 A JP 2002111925A JP 2002111925 A JP2002111925 A JP 2002111925A JP 4123811 B2 JP4123811 B2 JP 4123811B2
Authority
JP
Japan
Prior art keywords
particulate
oxygen concentration
exhaust
diesel engine
amount
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2002111925A
Other languages
Japanese (ja)
Other versions
JP2003307110A (en
Inventor
欣悟 陶山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toyota Motor Corp filed Critical Toyota Motor Corp
Priority to JP2002111925A priority Critical patent/JP4123811B2/en
Publication of JP2003307110A publication Critical patent/JP2003307110A/en
Application granted granted Critical
Publication of JP4123811B2 publication Critical patent/JP4123811B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D41/00—Electrical control of supply of combustible mixture or its constituents
    • F02D41/02—Circuit arrangements for generating control signals
    • F02D41/14—Introducing closed-loop corrections
    • F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
    • F02D41/1444—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases
    • F02D41/1466—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being a soot concentration or content
    • F02D41/1467—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being a soot concentration or content with determination means using an estimation
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D2200/00—Input parameters for engine control
    • F02D2200/02—Input parameters for engine control the parameters being related to the engine
    • F02D2200/08—Exhaust gas treatment apparatus parameters
    • F02D2200/0812—Particle filter loading

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Exhaust Gas After Treatment (AREA)

Description

【０００１】
【発明の属する技術分野】
本発明は、ディーゼル機関から排出されるパティキュレートの量（パティキュレート排出量）を推定するディーゼル機関のパティキュレート排出量推定装置に関する。
【０００２】
【従来の技術】
従来より、ディーゼル機関の排気通路にパティキュレートフィルタを配設した排気浄化装置が知られている。パティキュレートフィルタは、排気ガス中の煤（黒鉛、Soot）、未燃燃料や未燃潤滑油を主成分とする有機材可溶分の微粒子状物質（ＳＯＦ）、及び燃料中の硫黄分が反応して生成される有機材不溶分の微粒子状物質（サルフェート）等からなる微粒子状物質（これらの微粒子状物質は「パティキュレート・マター」と称呼されるが、以下、単に「パティキュレート」と称呼する。）を捕集するフィルタである。
【０００３】
ところで、パティキュレートフィルタは、捕集したパティキュレートにより目詰まりを起こすことがあり、目詰まりによって機関の排気抵抗を増大させる。このため、パティキュレートフィルタのパティキュレート捕集量を推定し、同パティキュレート捕集量が所定量以上とならないように同パティキュレートをフィルタ内で燃焼させる等のフィルタ再生制御を行うことが必要である。
【０００４】
このようなパティキュレート捕集量を算出（推定）する方法は、例えば、特開２０００−１７０５２１に開示されている。開示された方法においては、先ず、ディーゼル機関から排出されるパティキュレート排出量が推定される。
【０００５】
具体的に説明すると、種々の燃料噴射量及び機関の回転速度に対するパティキュレート排出量を予め測定し、その結果をテーブル１として記憶しておく。また、排気還流量（ＥＧＲ量）によりパティキュレート排出量が変化することから、排気還流量を決定するパラメータである機関回転速度、燃料噴射量、及び機関の冷却水温毎にＥＧＲ補正係数を求めておき、これをテーブル２として記憶しておく。
【０００６】
そして、実際の運転にあっては、燃料噴射量、機関回転速度、及び前記テーブル１からパティキュレート基本排出量を求めるとともに、燃料噴射量、機関回転数、冷却水温、及び前記テーブル２からＥＧＲ補正係数を求め、同ＥＧＲ補正係数でパティキュレート基本排出量を補正することによりパティキュレート排出量を推定する。
【０００７】
また、上記公報に開示された方法では、機関回転速度等に基づいてパティキュレートフィルタの捕集効率を求め、前記推定したパティキュレート排出量と同捕集効率とからパティキュレート捕集量を推定する。
【０００８】
【発明が解決しようとする課題】
しかしながら、燃料噴射量及び機関回転数が同一であっても、例えば、燃焼温度や燃焼圧力が異なればパティキュレート排出量は異なるので、燃料噴射量と機関回転数のみを引数としてパティキュレート排出基本量を求めたのでは、ＥＧＲ補正係数が正しい場合であっても、パティキュレート排出量を精度良く推定することができないという問題がある。
【０００９】
また、上記開示された従来の方法においては、排気還流量が燃料噴射量、機関回転速度、及び機関の冷却水温に応じて一義的に決るとの前提下でパティキュレート排出量を推定していて、機関が過渡運転状態となったときの排気還流量の変動（例えば、過給機の過給遅れに起因する排気還流量の上昇遅れ）を考慮していないので、ＥＧＲ補正係数を正しく求めることが困難であり、その結果、パティキュレート排出量を精度良く推定できないという問題がある。
【００１０】
従って、本発明の目的の一つは、ディーゼル機関から排出されるパティキュレート排出量を、簡単に且つ精度良く推定し得るパティキュレート排出量推定装置を提供することにある。
【００１１】
【本発明の概要】
上記目的を達成するための本発明に係るパティキュレート排出量推定装置は、ディーゼル機関の燃料噴射量及び機関回転速度を一定に維持した定常運転状態にて排気酸素濃度を変更するとともに、同変更された異なる排気酸素濃度のそれぞれに対してパティキュレート排出量を測定し、同変更した各排気酸素濃度及び同各排気酸素濃度に対して測定された同各パティキュレート排出量からなるデータに基いて、同一定に維持した燃料噴射量及び同一定に維持した機関回転速度で定まる運転状態に対するパティキュレートの生成定数及びパティキュレート酸化定数を決定する作業を、それぞれが所定の燃料噴射量及び所定の機関回転速度で定まる複数の運転状態について実行することにより形成された、燃料噴射量及び機関回転速度の組み合わせとパティキュレート生成定数との関係を規定した生成定数テーブルと、燃料噴射量及び機関回転速度の組み合わせとパティキュレート酸化定数との関係を規定した酸化定数テーブルとを記憶した記憶手段と、実際の燃料噴射量を表す値を取得する燃料噴射量取得手段と、実際の機関回転速度を検出する機関回転速度検出手段と、実際の排気酸素濃度を表す値を取得する排気酸素濃度取得手段と、前記取得された燃料噴射量を表す値と、前記検出された機関回転速度と、前記生成定数テーブル及び前記酸化定数テーブルと、に基いてパティキュレート排出量を推定するために使用する生成定数と酸化定数とを決定するとともに、同決定された生成定数と、同決定された酸化定数と、前記取得された排気酸素濃度を表す値とに基いてパティキュレート排出量を推定するパティキュレート排出量推定手段と、を備える。
【００１２】
更に、前記パティキュレート排出量推定手段は、単位時間あたりのパティキュレート排出量をＰＭ、前記決定されたパティキュレート生成定数をＡsft、前記決定されたパティキュレート酸化定数をＡsot、前記取得された燃料噴射量を表す値に基く単位時間あたりの燃料噴射量をＱ、前記取得された排気酸素濃度を表す値をＸ02、及び定数をＣとするとき、ＰＭ＝Ａsft・Ｑ／（Ａsot・ＸO2）＋Ｃ／exp（Ａsot・ＸO2）なる式に基いて前記パティキュレート排出量を推定するように構成されている。
【００１３】
本発明者は、種々の検討の結果、ディーゼル機関が定常運転状態にあるか過渡運転状態にあるかに関わらず、単位時間あたりのパティキュレート排出量ＰＭが、ＰＭ＝Ａsft・Ｑ／（Ａsot・ＸO2）＋Ｃ／exp（Ａsot・ＸO2）なる式に実質的に従って変化するとの知見を得た。ここで、Ａsft及びＡsotはそれぞれ１ストロークあたりの燃料噴射量とエンジン回転数とにより一義的に定まるパティキュレート生成定数及びパティキュレート酸化定数であり、Ｑは１ストロークあたりの燃料噴射量と機関回転速度とにより定まる量である単位時間あたりの燃料噴射量、Ｘ02は排気酸素濃度、及びＣは定数である。
【００１４】
このように、パティキュレート排出量ＰＭは、パティキュレート生成定数Ａsft、及びパティキュレート酸化定数Ａsotを特定できれば、排気酸素濃度ＸO2の関数として簡易的に求めることができるのである。
【００１５】
そこで、先ず、本発明のパティキュレート排出量推定装置は、燃料噴射量及び機関回転速度の組み合わせとパティキュレート生成定数との関係を規定した生成定数テーブルと、燃料噴射量及び機関回転速度の組み合わせとパティキュレート酸化定数との関係を規定した酸化定数テーブルとを作成し、記憶手段に記憶させておく。
【００１６】
これらのテーブルは、ディーゼル機関の燃料噴射量及び機関回転速度を一定に維持した定常運転状態（この場合、機関回転速度が一定であるから、燃料噴射量はストロークあたりの燃料噴射量Ｑfinであっても、単位時間あたりの燃料噴射量Ｑであってもよい。）において排気酸素濃度を変更するとともに、同変更された異なる排気酸素濃度のそれぞれに対してパティキュレート排出量を測定し、同変更した各排気酸素濃度及び同各排気酸素濃度に対して測定された同各パティキュレート排出量からなるデータ点を上述した式が通るようにパティキュレート生成定数Ａsft及びパティキュレート酸化定数Ａsotを決定（同定又はフィッティング）し、かかる作業を、種々の燃料噴射量及び機関回転速度に対して繰り返し行うことにより、作成される。
【００１７】
そして、実際の燃料噴射量を表す値を取得するとともに実際の機関回転速度を検出し、取得した燃料噴射量を表す値と検出した機関回転速度と前記記憶している生成定数テーブル及び前記記憶している酸化定数テーブルとに基いて、実際のパティキュレート排出量を推定するために使用するパティキュレート生成定数とパティキュレート酸化定数とを決定し、これら決定した定数と単位時間当たりの燃料噴射量と前記取得された排気酸素濃度とを上記数式に代入してパティキュレート排出量を推定する。即ち、少なくとも、上記決定したパティキュレート生成定数と、上記決定したパティキュレート酸化定数と、上記取得された排気酸素濃度とに基いてパティキュレート排出量を推定する。
【００１８】
このように、上記本発明の特徴を備えたパティキュレート排出量推定装置はパティキュレート排出量を推定する際の変数として排気酸素濃度を用いていて、同排気酸素濃度は燃焼温度や燃焼圧力等に比べて比較的測定精度を確保し易い変数であり、しかも、上記パティキュレート生成定数Ａsftと上記パティキュレート酸化定数Ａsotは、燃料噴射量と機関回転速度とを一定に維持した場合において実際に測定したパティキュレート排出量により決定されるので、燃料噴射量と機関回転速度とで定まる各運転状態毎での燃焼温度、燃焼圧力、燃焼割合、燃焼形式、噴射弁の形状、及び機関の排気量等のパティキュレート排出量に影響を及ぼす総べてのパラメータが反映された値となっている。この結果、本発明のパティキュレート排出量推定装置は、パティキュレート排出量を精度良く推定することができる。
【００１９】
また、本発明のパティキュレート排出量推定装置は、排気還流量（ＥＧＲ量）の変動分によっても変動する実際の排気酸素濃度を表す値に基いてパティキュレート排出量を推定するので、機関の運転状態が過渡運転状態となって排気還流量が予定通りに得られないような場合であっても、パティキュレート排出量を精度良く求めることができる。
【００２０】
このようなパティキュレート排出量推定装置においては、前記パティキュレート排出量推定手段が、前記取得された排気酸素濃度を表す値が所定濃度以上の値となるときにのみ前記パティキュレート排出量の推定を行うように構成されることが好適である。
【００２１】
上記の式は、ディーゼル機関に吸入されるガスの空燃比が例えば約２０よりも大きな空燃比である場合にのみ有効であって、同空燃比Ａ／Ｆを２０以下として燃焼温度を低温とすることでパティキュレート排出量を抑制する燃焼（所謂、低温燃焼）には対応していない。従って、上記構成のように、取得された排気酸素濃度を表す値が所定濃度以上の値となるときにのみ前記パティキュレート排出量の推定を行うように構成することで、パティキュレート排出量の推定精度の低下を回避することができる。
【００２２】
ところで、一般に、機関の排気酸素濃度（パティキュレート排出量推定に用いるべき排気酸素濃度）に対して、排気酸素濃度検出手段により検出される排気酸素濃度は時間的に遅れて変化する。また、前記ディーゼル機関が排気還流管を介して排気還流を行うとともに同排気還流の量を制御する排気還流手段（ＥＧＲ装置）を備えている場合、運転状態の変化に基いて排気酸素濃度が変化した排気は、排気還流遅れ時間だけ遅れて同機関に吸入されることになる。従って、運転状態の変化に伴う排気酸素濃度変化が検出されてから排気還流遅れ時間だけ経過した時点においては、機関の排気酸素濃度（パティキュレート排出量推定に用いるべき排気酸素濃度）は更に変化を開始しているはずである。
【００２３】
従って、このようなパティキュレート排出量推定装置であって、前記ディーゼル機関が排気還流管を介して排気還流を行うとともに同排気還流の量を制御する排気還流手段を備えている場合、前記排気酸素濃度取得手段は、同ディーゼル機関の排気通路内の酸素濃度を検出する排気酸素濃度検出手段を含むとともに、同ディーゼル機関が備える排気還流手段によって同機関の排気が還流されて同機関に吸入されるまでに要する排気還流遅れ時間に応じ、同検出された排気通路内の酸素濃度を補正して前記排気酸素濃度を表す値を取得するように構成されることが好適である。
【００２４】
これによれば、運転状態の変化に伴う排気酸素濃度の変化が排気還流により更なる排気酸素濃度の変化を発生させる場合、同更なる変化を排気酸素濃度検出手段が検出する前の適切な時点で、パティキュレート排出量推定に用いる排気酸素濃度に反映することができ、その結果、同パティキュレート排出量推定用排気酸素濃度が適切な値となるので、パティキュレート排出量の推定精度を向上することができる。
【００２５】
この場合、前記排気酸素濃度取得手段は、前記還流される排気と前記ディーゼル機関に新たに吸入される新気とが混合される直前の前記排気還流管内の酸素濃度を検出する排気還流遅れ検出用酸素濃度検出手段を含むとともに、前記検出される排気通路内の酸素濃度と同検出される排気還流管内の酸素濃度とに基いて前記排気還流遅れ時間を推定するように構成されることが好適である。
【００２６】
これによれば、排気通路内の検出される酸素濃度と、還流される排気とディーゼル機関に新たに吸入される新気とが混合される直前の排気還流管内の検出される酸素濃度と、に基いて排気還流遅れ時間が推定されるので、同排気還流遅れ時間を実測値に基いて精度良く決定することができ、その結果、パティキュレート排出量推定用排気酸素濃度を適切な値とすることができるので、パティキュレート排出量の推定精度を向上することができる。
【００２７】
また、前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼温度に関連する値を取得するとともに、同取得された燃焼温度に関連する値に基いて前記パティキュレート排出量を推定するように構成されることが好適である。
【００２８】
燃焼温度はパティキュレート排出量に比較的強い影響を及ぼすから、同燃焼温度をパティキュレート排出量の推定に用いることにより、一層精度良くパティキュレート排出量を推定することができる。
【００２９】
また、前記パティキュレート排出量推定手段は、前記ディーゼル機関の吸気通路内の酸素濃度を検出する酸素濃度検出手段を含み、同検出される吸気通路内の酸素濃度を前記燃焼温度に関連する値として扱うように構成されることが好適である。吸気通路内の酸素濃度は燃焼温度と密接な関係にあるから、同吸気通路内の酸素濃度を燃焼温度に関連する値として扱うことにより、パティキュレート排出量を精度良く推定することが可能となる。
【００３０】
また、前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼室壁面の温度を検出する燃焼室壁面温度検出手段を含み、同検出される燃焼室壁面温度を前記燃焼温度に関連する値として扱うように構成されることが好適である。燃焼室壁面温度は、例えば、燃焼室壁面温度に熱伝達効率を乗じることにより燃焼平均温度が求まるように、燃焼温度と密接な関係にあるから、同燃焼室壁面温度を燃焼温度に関連する値として扱うことにより、パティキュレート排出量を精度良く推定することが可能となる。
【００３１】
また、前記パティキュレート排出量推定手段は、前記ディーゼル機関の吸気温度を検出する吸気温度検出手段を含み、同検出される吸気温度を前記燃焼温度に関連する値として扱うように構成されることが好適である。吸気温度は燃焼温度と密接な関係にあるから、同吸気温度を燃焼温度に関連する値として扱うことにより、パティキュレート排出量を精度良く推定することが可能となる。
【００３２】
また、前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼室内の圧力を検出する燃焼室内圧力検出手段を含み、同検出される燃焼室内圧力を前記燃焼温度に関連する値として扱うように構成されることが好適である。例えば、圧縮上死点での燃焼室内圧力は燃焼温度と密接な関係にあるから、同燃焼室内圧力を燃焼温度に関連する値として扱うことにより、パティキュレート排出量を精度良く推定することが可能となる。
【００３３】
更に、このようなパティキュレート排出量推定装置であって、前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼中における燃焼室内圧力に関連する値を取得するとともに、同取得された燃焼室内圧力に関連する値に基いて前記パティキュレート排出量を推定するように構成されることが好適である。
【００３４】
燃焼中における燃焼室内圧力、特に燃焼圧と呼ばれる燃焼室内圧力の極大値は、パティキュレート排出量に比較的強い影響を及ぼすから、同燃焼中における燃焼圧力に関する値を取得してパティキュレート排出量の推定に用いることにより、一層精度良くパティキュレート排出量を推定することができる。
【００３５】
更に、このようなパティキュレート排出量推定装置であって、前記ディーゼル機関の燃焼室内の煤発生量を測定する煤発生量測定手段と、前記測定された燃焼室内の煤発生量に基いて基準パティキュレート排出量を算定する基準パティキュレート排出量算出手段とを備えるとともに、前記パティキュレート排出量推定手段は、前記基準パティキュレート排出量が算出された運転状態において推定された前記パティキュレート排出量と同算出された基準パティキュレート排出量との差に応じて、他の運転状態において推定されるパティキュレート排出量を修正するように構成されることが好適である。
【００３６】
これによれば、煤の発生量が実際に測定される。パティキュレートは、煤のみでなく上述したＳＯＦ及びサルフェートを含んでいるが、ＳＯＦの量とサルフェートの量は、煤の量に対して少ないこと、及び煤の量と相関関係があることから、例えば、エンジンが定常運転状態にあるときなどの特定の運転状態においては、測定された煤の発生量に基いてパティキュレート排出量の基準量、即ち、基準パティキュレート排出量が比較的精度良く算出され得る。
【００３７】
そして、前記パティキュレート排出量推定手段により、同じ運転状態における前記推定されたパティキュレート排出量と前記算出された基準パティキュレート排出量との差に応じて、他の運転状態で推定されるパティキュレート排出量が修正される。換言すると、パティキュレート排出量推定に用いる上記式で表されたモデルの誤差が実測値に基いた値で修正される。従って、パティキュレート排出量を精度良く推定することが可能となる。
【００３８】
【発明の実施の形態】
以下、本発明によるディーゼル機関のパティキュレート排出量推定装置の各実施形態について図面を参照しつつ説明する。
【００３９】
（１．第１実施形態）
図１は、本発明の第１実施形態に係るパティキュレート排出量推定装置を含む内燃機関の排気浄化装置を４気筒内燃機関（ディーゼル機関）１０に適用したシステムの概略構成を示している。このシステムは、燃料供給系統を含むエンジン本体２０、エンジン本体２０の各気筒の燃焼室に大気（新気）を導入するための吸気系統３０、エンジン本体２０からの排気ガスを放出するための排気系統４０、排気還流を行うためのＥＧＲ装置５０、及び電気制御装置６０を含んでいる。
【００４０】
エンジン本体２０の各気筒の上部には燃料噴射弁２１が配設されている。各燃料噴射弁２１は、図示しない燃料タンクと接続された燃料噴射用ポンプ２２に燃料配管２３を介して接続されている。これにより、燃料噴射弁２１には、燃料噴射用ポンプ２２から噴射圧力まで昇圧された燃料が供給されるようになっている。また、燃料噴射弁２１は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（指令信号）に応じて所定時間だけ開弁し、これにより各気筒の燃焼室内に前記昇圧された燃料を噴射するようになっている。
【００４１】
吸気系統３０は、エンジン本体２０の各気筒の燃焼室にそれぞれ接続された吸気マニホールド３１、吸気マニホールド３１の上流側集合部に接続され同吸気マニホールド３１とともに吸気通路を構成する吸気管３２、吸気管３２内に回動可能に保持されたスロットル弁３３、電気制御装置６０からの駆動信号に応答してスロットル弁３３を回転駆動するスロットル弁アクチュエータ３３ａ、スロットル弁３３の上流で吸気管３２に介装されたインタクーラー３４と過給機３５のコンプレッサ３５ａ、吸気管３２の先端部に配設されたエアクリーナ３６とを含んでいる。
【００４２】
排気系統４０は、エンジン本体２０の各気筒にそれぞれ接続された排気マニホールド４１、排気マニホールド４１の下流側集合部に接続された排気管４２、排気管４２に配設された過給機３５のタービン３５ｂ、排気管４２に介装されるとともにＮＯｘ触媒及び捕集したパティキュレートを酸化する触媒を担持したディーゼルパティキュレートフィルタ（以下、「ＤＰＮＲ」と称呼する。）４３、ＤＰＮＲ４３の上流の排気管４２に配設され添加燃料（軽油）をＤＰＮＲ４３に供給する添加燃料供給ノズル４４、添加燃料供給ノズル４４と燃料噴射用ポンプ２２とを接続する添加燃料供給管４５、添加燃料供給管４５に介装され電気制御装置６０からの駆動信号に応答して同添加燃料供給管４５の燃料通路を開閉する電磁開閉弁４４ａ、ＤＰＮＲ４３の下流の排気管内に回動可能に保持された排気絞り弁４６、及び電気制御装置６０からの駆動信号に応答して排気絞り弁４６を回転駆動する排気絞り弁アクチュエータ４６ａを含んでいる。排気マニホールド４１及び排気管４２は排気通路を構成している。
【００４３】
ＤＰＮＲ４３は、コージライト（セラミックの一種）のような多数の細孔を有する材料（即ち、多孔質材料）から形成されたフィルタ４３ａを備え、通過する排気ガス中のパティキュレートを細孔表面にて捕集する、ハニカム構造のウオールフロー型フィルタである。
【００４４】
また、ＤＰＮＲ４３のフィルタの表面には、触媒の担体としてのアルミナ（Ａｌ2Ｏ3）がコートされていて、同アルミナに捕集したパティキュレートを酸化させる（及び／又は、酸化作用を有する）例えばＰｔ等の貴金属からなる触媒が担持されている。
【００４５】
ところで、このＤＰＮＲ４３内で酸化されることにより同ＤＰＮＲ４３の外部に放出されるパティキュレートの量よりも同ＤＰＮＲ４３に捕集されるパティキュレートの量が多くなると、同ＤＰＮＲ４３内にパティキュレートが堆積する。その結果、ＤＰＮＲ４３のパティキュレート捕集能力が低下するとともに、機関１０の排気抵抗が増大する。
【００４６】
このため、本排気浄化装置は、前述した添加燃料供給ノズル４４からＤＰＮＲ４３に添加燃料を供給し、同添加燃料がＤＰＮＲ４３の内部で燃焼することに伴って発生する熱により触媒床温を上昇せしめ、パティキュレートの酸化速度を上昇させることで堆積したＤＰＮＲ４３をより多く酸化して除去する。即ち、本排気浄化装置は、添加燃料の供給によりフィルタの再生を定期的に行うことで、ＤＰＮＲ４３のパティキュレート捕集能力を維持するとともに、排圧の上昇を防止するようになっている。なお、パティキュレートが酸化されるためには酸素が必要であるから、添加燃料供給ノズル４４からの添加燃料の供給量は、ＤＰＮＲ４３の全体の空燃比がリッチとならない程度（即ち、ＤＰＮＲ４３内にパティキュレートを酸化させるための酸素が十分に存在するリーン空燃比が得られる程度）の量に設定される。
【００４７】
更に、ＤＰＮＲ４３は、前記担体としてのアルミナに、カリウムＫ，ナトリウムＮａ，リチウムＬｉ，セシウムＣｓのようなアルカリ金属、バリウムＢａ，カルシウムＣａのようなアルカリ土類金属、及びランタンＬａ、イットリウムＹのような希土類金属から選ばれた少なくとも一つを、前記担持した白金とともに担持している。
【００４８】
これにより、ＤＰＮＲ４３は、同ＤＰＮＲ４３に流入するガスの空燃比が理論空燃比よりもリーンのときにＮＯｘを吸収し、同流入するガスの空燃比が理論空燃比、又は理論空燃比よりもリッチのときに吸収したＮＯｘをＮＯ2又はＮＯとして放出するようになっている。また、この放出されたＮＯ2又はＮＯはＤＰＮＲ４３に流入するリッチ空燃比のガス中の未燃ＨＣや未燃ＣＯと反応してＮ2に還元せしめられて外部に放出される。このように、ＤＰＮＲ４３は、ＮＯｘを吸収した後に同吸収したＮＯｘを放出して還元する吸蔵還元型ＮＯｘ触媒としても機能するようになっている。
【００４９】
ＥＧＲ装置５０は、排気ガスを還流させる通路（再循環通路）を構成する排気還流管５１と、排気還流管５１に介装されたＥＧＲ弁５２と、ＥＧＲクーラー５３とを備えている。排気還流管５１はタービン３５ｂの上流側排気通路（排気マニホールド４１）とスロットル弁３３の下流側吸気通路（吸気マニホールド３１）を連通している。ＥＧＲ弁５２は電気制御装置６０からの駆動信号に応答し、再循環される排気ガス量（排気還流量、ＥＧＲ量）が、機関回転速度NEと一吸気行程あたりの筒内燃料噴射量Ｑfin（又は、アクセル操作量Accpに基づく要求トルクTQ）とに応じて決定される目標ＥＧＲ量と等しくなるように、排気還流管５１により構成されるＥＧＲ通路の開口断面積を変更するようになっている。なお、ＥＧＲ弁により、ＥＧＲ率を制御するように構成してもよい。
【００５０】
電気制御装置６０は、互いにバスで接続されたＣＰＵ６１、ＣＰＵ６１が実行するプログラム、テーブル（ルックアップテーブル、マップ）、定数等を予め記憶したＲＯＭ６２、ＣＰＵ６１が必要に応じてデータを一時的に格納するＲＡＭ６３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ６４、及びＡＤコンバータを含むインターフェース６５等からなるマイクロコンピュータである。
【００５１】
インターフェース６５は、新気量計測手段であって吸気管３２に配置された熱線式エアフローメータ７１、スロットル弁３３の下流の吸気管３２に設けられた吸気温センサ７２、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路（吸気マニホールド３１）に配置された吸気酸素濃度センサ７３、第３気筒の燃焼室上部に配置された燃焼圧センサ７４、エンジン本体２０のシリンダボアライナー部に配設された壁面温度センサ７５、第１気筒の燃焼室に対して配置された光学式煤センサ７６、及びタービン３５ｂの下流側でＤＰＮＲ４３の上流側の排気管４２内に配置された排気酸素濃度センサ７７、機関回転速度センサ７８、排気還流遅れ検出用酸素濃度センサ７９、及びアクセル開度センサ８１と接続されていて、これらのセンサからの信号をＣＰＵ６１に供給するようになっている。
【００５２】
また、インターフェース６５は、燃料噴射弁２１、スロットル弁アクチュエータ３３ａ、添加燃料供給管４５の燃料通路を開閉する電磁開閉弁４４ａ、排気絞り弁アクチュエータ４６ａ、及びＥＧＲ弁５２と接続されていて、ＣＰＵ６１の指示に応じてこれらに駆動信号を送出するようになっている。
【００５３】
熱線式エアフローメータ７１は、吸気通路内を通過する吸入空気の質量流量（単位時間当りの吸入空気量、単位時間あたりの新気量）を計測し、同質量流量を表す信号Ｇaを発生するようになっている。吸気温センサ７２は、エンジン１０の気筒（即ち、燃焼室）に吸入される新気の温度（即ち、吸気温度）を検出し、同吸気温度を表す信号Ｔｂを発生するようになっている。吸気酸素濃度センサ７３は、エンジン１０の気筒（即ち、燃焼室）に吸入されるガス中の酸素濃度（即ち、吸気酸素濃度）を検出し、同吸気酸素濃度を表す信号ＸO2INを発生するようになっている。
【００５４】
燃焼圧センサ７４は、第３気筒の燃焼室内の圧力を検出し、同燃焼室内圧力Ｐａを表す信号を発生するようになっている。壁面温度センサ７５は、燃焼室壁面温度を検出し、同燃焼室壁面温度を表す信号Ｔhekiを表す信号を発生するようになっている。
【００５５】
光学式煤センサ７６は、第１気筒のシリンダ壁面であって、同第１気筒のピストンが上死点にあるときに、同第１気筒のピストンの最上位置に配された第１ピストンリングの下方であって、同第１気筒のピストンが下死点に向けて僅かに移動したときに同第１ピストンリングの上方に位置する部分に透光性の窓を備えるとともに、同窓の反射率を測定することにより同第１気筒の燃焼室内に発生した煤量Ｓsootを測定するようになっている。
【００５６】
排気酸素濃度センサ７７は、エンジン１０の気筒（即ち、燃焼室）から排出された排気ガス中の酸素濃度（即ち、排気酸素濃度）を検出し、同排気酸素濃度を表す信号ＸO2Rを発生するようになっている。機関回転速度センサ７８は、エンジン１０の回転速度を検出し、機関回転速度NEを表す信号を発生するとともに、各気筒の絶対クランク角度を検出し得るようになっている。排気還流遅れ検出用酸素濃度センサ７９は、エンジン１０に新たに吸入される新気と還流された排気とが混合される直前の排気還流管５１内の酸素濃度を検出する排気還流遅れ検出用酸素濃度検出手段を構成するもので、排気還流管５１のＥＧＲ弁５２の近傍に設けられていて、ＥＧＲ弁５２に導入される直前の排気管流管５１内の酸素濃度（還流排気中酸素濃度）を検出し、同還流排気中酸素濃度ＸO2INRを表す信号を発生するようになっている。アクセル開度センサ８１は、アクセルペダルＡＰの操作量を検出し、アクセル操作量を表す信号Accpを発生するようになっている。
【００５７】
＜第１実施形態のパティキュレート排出量推定方法＞
次に、上記のように構成された排気浄化装置に含まれるパティキュレート排出量推定装置によるパティキュレート排出量の推定方法について説明する。このパティキュレート排出量推定装置は、下記数１（数１で記述されたパティキュレート排出量推定モデル）に基いて単位時間あたりにエンジン本体２０から排出されるパティキュレート排出量ＰＭを推定する。なお、数１において、Ａsftはパティキュレート生成定数（以下、単に「生成定数」とも呼ぶ。）、Ａsotはパティキュレート酸化定数（以下、単に「酸化定数」とも呼ぶ。）、Ｑは単位時間あたりの燃料噴射量、Ｘ02は排気酸素濃度（酸素O2のモル分率）、及びＣは予め定められた定数である。
【００５８】
【数１】
ＰＭ＝Ａsft・Ｑ／（Ａsot・ＸO2）＋Ｃ／exp（Ａsot・ＸO2）
【００５９】
ここで、上記数１の根拠について簡単に説明する。下記数２，数３は、煤の生成及び酸化についてのモデルとしてＳＡＥペーパー９３０６１２（タイトル「Approach to Low Nox and Smoke Emission Engines by Using Phenomenological Simulation」）にて紹介された実験式である広安モデル（Sootモデル）を記述した式である。数２，数３において、Ｍsfは煤の生成量、Ａsfは定数、Ｍfvは１ストロークあたりの燃料噴射量Ｍfの燃料のうち気化率（燃焼割合）RHRで気化して燃焼した燃料量、Ｐは燃焼圧力、Ｒは気体定数、Ｔは燃焼温度、Ｍscは煤の酸化量、Ａsoは定数、及びＭsは煤の排出量である。
【００６０】
【数２】
dＭsf／dt＝Ａsf・Ｍfv・Ｐ0.5・exp（−12500／(Ｒ・Ｔ)）
【００６１】
【数３】
dＭsc／dt＝Ａso・Ｍs・ＸO2・Ｐ1.8・exp（−14000／(Ｒ・Ｔ)）
【００６２】
煤の生成量の変化量dＭs／dtは、煤の生成量（dＭsf／dt）から煤の酸化量（dＭsc／dt）を減算した値であるから、下記数４が成立する。
【００６３】
【数４】
dＭs／dt＝dＭsf／dt−dＭsc／dt
【００６４】
ここで、燃焼温度Ｔ、燃焼圧力Ｐ、及び燃焼割合RHRを、機関回転速度NEと単位時間あたりの燃料噴射量Ｑが変化した場合であっても、それぞれ一定であると仮定し、Ｍfv＝Ａrhr・Ｍf、Ｐ0.5＝Ａpf、exp（−12500／(Ｒ・Ｔ)）＝Ａtf、Ｐ1.8＝Ａpo、及びexp（−14000／(Ｒ・Ｔ)）＝Ａtoとおくと、上記数２〜上記数４より下記数５が得られる。なお、値Ａrhr、Ａpf、Ａtf、Ａpo、Ａtoは、定数である。
【００６５】
【数５】
dＭs／dt＝Ａsf・Ａrhr・Ｍf・Ａpf・Ａtf−Ａso・Ｍs・ＸO2・Ａpo・Ａto
【００６６】
そこで、パティキュレートの生成定数Ａsft、及びパティキュレートの酸化定数Ａsotをそれぞれ下記数６，数７のように定義すると、上記数５は下記数８のように書き換えられる。
【００６７】
【数６】
Ａsft＝Ａsf・Ａrhr・Ａpf・Ａtf
【００６８】
【数７】
Ａsot＝Ａso・Ａpo・Ａto
【００６９】
【数８】
dＭs／dt＝Ａsft・Ｍf−Ａsot・Ｍs・ＸO2
【００７０】
数８で表されたＭsについての微分方程式を解くと、下記数９が得られる。
【００７１】
【数９】
Ｍs＝Ａsft・Ｍf／（Ａsot・ＸO2）＋Ｃ１／exp（Ａsot・ＸO2・t）
【００７２】
いま、エンジン１０の運転状態が、燃料噴射量Ｑfin、及び機関回転速度NEが変化しない定常運転状態であると仮定すると、上記数９の右辺第２項は時間に依存しないと考えることができる。また、上記数２，数３により表される広安モデルは煤についてのモデルであるところ、パティキュレートにはＳＯＦやサルフェートが含まれているから、この煤以外のパティキュレートの排出量を生成定数Ａsft及び酸化定数Ａsotの値に反映させる。そして、これらの考えと、上記数９の単位を変換することにより、上記数１と同一である下記数１０が得られる。この数１０が、本パティキュレート排出量推定装置が用いるモデルを記述した式である。
【００７３】
【数１０】
ＰＭ＝Ａsft・Ｑ／（Ａsot・ＸO2）＋Ｃ／exp（Ａsot・ＸO2）
【００７４】
次に、上記数１０（数１）を使用してパティキュレート排出量ＰＭを求める方法について説明する。上記数１０において、Ｑは単位時間あたりの燃料噴射量（単位は、g/h）であり、１ストロークあたりの燃料噴射量Ｑfin（単位は、mm3/ストローク）と機関回転速度NEとに基いて求めることができる。例えば、４気筒ディーゼル機関であれば、クランク軸の一回転（エンジン一回転）あたり２回の爆発行程が存在するから、Ｑ＝Ｋ・Ｑfin・２・NE（Ｋは定数）である。
【００７５】
そこで、燃料噴射量Ｑfinと機関回転速度NEとを一定の状態（Ｑfin＝Ｑfin０，NE＝NE０）に維持し、例えば、ＥＧＲ弁５２の開度を変更することにより排気酸素濃度ＸO2を変更し、そのときの排気酸素濃度ＸO2とパティキュレート排出量ＰＭとを実測する。そして図２に示したように、異なる排気酸素濃度ＸO2に対するパティキュレート排出量ＰＭの実測データ点を少なくとも３点Ｐ１〜Ｐ３ほど得て、上記数１０で表される曲線が同３点Ｐ１〜Ｐ３を通るように、生成定数Ａsftと酸化定数Ａsotとを、例えば、最小二乗法により決定する。
【００７６】
以降、同様に、燃料噴射量Ｑfinと機関回転速度NEとを異なる値に維持しながら排気酸素濃度ＸO2を変更してパティキュレート排出量ＰＭを測定し、そのデータに基いて同燃料噴射量Ｑfinと同機関回転速度NEとにおける生成定数Ａsftと酸化定数Ａsotとを決定し、図３のフローチャートにおけるステップ３０５、及びステップ３１０中にそれぞれ示したように、燃料噴射量Ｑfin（実質的に燃料噴射量Ｑfinに応じた値であればよい）及び機関回転速度NEに対するパティキュレート生成定数Ａsftの関係を規定した生成定数テーブルと、燃料噴射量Ｑfin（実質的に燃料噴射量Ｑfinに応じた値であればよい）及び機関回転速度NEとパティキュレート酸化定数Ａsotの関係を規定した酸化定数テーブルを作成し、ＲＯＭ６２内に格納しておく。
【００７７】
そして、実際のエンジン運転においては、その時点における燃料噴射量Ｑfin及び機関回転速度NEと上記ＲＯＭ６２内に格納した二つのテーブルとから生成定数Ａsftと酸化定数Ａsotとを決定し、決定した生成定数Ａsftと酸化定数Ａsotとを上記数１０（数１）に適用するとともに、実際の単位時間当たりの燃料噴射量Ｑを数１に代入し、排気酸素濃度センサ７７が検出する実際の排気酸素濃度ＸO2Rを排気酸素濃度ＸO2として同数１０に代入しパティキュレート排出量ＰＭを推定する。以上が、本実施形態のパティキュレート排出量推定方法である。
【００７８】
なお、上記数１０によれば、パティキュレート排出量ＰＭは排気酸素濃度ＸO2に対して単調減少する。しかしながら、この傾向はエンジン１０に吸入されるガスの空燃比が例えば約２０よりも大きな空燃比である場合にのみ有効であって、同空燃比Ａ／Ｆを２０以下として燃焼温度を低温とすることでパティキュレート排出量を抑制する燃焼（所謂、低温燃焼）には対応していない。従って、上記パティキュレート排出量推定方法は、エンジン１０に吸入されるガスの空燃比が低温燃焼をもたらす所定の空燃比の最大値（約２０）よりも大きい場合にのみ適用することとする。
【００７９】
＜第１実施形態の作動＞
次に、本パティキュレート排出量推定装置を含む排気浄化装置の作動について説明すると、電気制御装置６０のＣＰＵ６１はアクセル開度センサ８１により検出されるアクセル操作量Accpと機関回転速度センサ７８により検出される機関回転速度NEとＲＯＭ６２内に格納されているテーブルとに基いて燃料噴射量Ｑfinを決定し、同決定された燃料噴射量Ｑfinの燃料が爆発行程を迎えた気筒に対して噴射されるように、同気筒に対する燃料噴射弁２１に駆動信号を発生する。
【００８０】
同様に、ＣＰＵ６１は前記検出されるアクセル操作量Accpと前記検出される機関回転速度NEとＲＯＭ６２内に格納されているテーブルとに基いて目標ＥＧＲ量（又は、目標ＥＧＲ率）を決定し、排気還流管５１を介して還流される排気量が同決定された目標ＥＧＲ量となるように、ＥＧＲ弁５２に駆動信号を発生する。
【００８１】
また、ＣＰＵ６１は、図３にフローチャートにより示したパティキュレート排出量ＰＭ推定ルーチンを所定時間TSAMPLEの経過毎に実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ３００から処理を開始してステップ３０５に進み、その時点の機関回転速度NE及び燃料噴射量Ｑfinと上述した生成定数テーブルとに基づいて生成定数Ａsftを決定し、次のステップ３１０にて、同機関回転速度NE及び同燃料噴射量Ｑfinと上述した酸化定数テーブルとに基づいて酸化定数Ａsotを決定する。
【００８２】
次に、ＣＰＵ６１は、ステップ３１５に進んで排気酸素濃度センサ７７が検出している排気酸素濃度ＸO2Rを排気酸素濃度ＸO2として取得し、ステップ３２０に進んで、ステップ３０５にて決定した生成定数Ａsft、ステップ３１０にて決定した酸化定数Ａsot、ステップ３１５にて取得した排気酸素濃度ＸO2、及び上記数１（数１０）に基いてパティキュレート排出量ＰＭを推定する。なお、ステップ３２０に記載したように、上記数１における単位時間あたりの燃料噴射量Ｑは下記数１１により求められる。
【００８３】
【数１１】
Ｑ＝Ｋ・Ｑfin・２・NE
【００８４】
次いで、ＣＰＵ６１はステップ３２５に進み、ステップ３２０にて推定したパティキュレート排出量ＰＭに、本ルーチンの実行時間間隔TSAMPLE（秒）を３６００（秒）で除した値を乗じて、その結果を現時点のパティキュレート排出量合計値SUMPMに加え、新たなパティキュレート排出量合計値SUMPMを得る。その後、ＣＰＵ６１は、ステップ３９５にて本ルーチンを一旦終了する。以降、ＣＰＵ６１は本ルーチンを所定時間TSAMPLEの経過毎に繰り返し実行する。この結果、パティキュレート排出量ＰＭとパティキュレート排出量合計値SUMPMが更新されて行く。
【００８５】
また、ＣＰＵ６１は、別途図示しないＤＰＮＲ４３を再生させるフィルタ再生制御ルーチンを実行し、上記パティキュレート排出量合計値SUMPMが所定値以上となると、添加燃料供給ノズル４４からＤＰＮＲ４３に添加燃料を供給し、同添加燃料がＤＰＮＲ４３の内部で燃焼することに伴って発生する熱により触媒床温を上昇せしめ、パティキュレートの酸化速度を上昇させることで堆積したＤＰＮＲ４３をより多く酸化する。
【００８６】
以上、説明したように、第１実施形態に係るパティキュレート推定装置によれば、排気酸素濃度ＸO2と燃料噴射量Ｑのみを変数とした関数（上記数１のモデル）に基いてパティキュレート排出量ＰＭが推定される。即ち、数１にて示されるモデルは変数が少なく簡素なモデルであって、その計算も単純であるから、通常の制御用コンピュータによってもパティキュレート排出量ＰＭをリアルタイムに推定することができる。しかも、そのモデルの変数である排気酸素濃度は、燃焼温度や燃焼圧力等に比べて比較的測定精度を確保し易い変数であるから、パティキュレート排出量ＰＭの推定を精度良く行うことができる。
【００８７】
また、上記数１のモデルにおけるパティキュレート生成定数Ａsftとパティキュレート酸化定数Ａsotは、燃料噴射量Ｑfinと機関回転速度NEとを一定に維持した場合において、実際に測定したパティキュレート排出量ＰＭにより決定される。従って、生成定数Ａsftと酸化定数Ａsotは、燃料噴射量Ｑfinと機関回転速度NEとで定まる各運転状態毎での燃焼温度、燃焼圧力、燃焼割合、燃焼形式、燃料噴射弁２１の形状、及びエンジン１０の排気量等のパティキュレート排出量に影響を及ぼす総べてのパラメータが反映された値となる。この結果、適合に要する工数を少なくしつつ、精度良くパティキュレート排出量を推定することができる。
【００８８】
また、排気酸素濃度ＸO2の関数としてパティキュレート排出量ＰＭを推定するから、例えば、過給機付きエンジンにおいて顕著に表れる過渡運転（加速運転）状態での新気の燃焼室流入遅れ、或いは過渡運転状態への移行に伴うＥＧＲ量変化遅れによるパティキュレート排出量への影響に対しても、これを排気酸素濃度で捉えることができるので、過渡運転状態におけるパティキュレート排出量を精度良く推定することが可能となる。
【００８９】
更に、ディーゼル機関１０の適合過程において、燃料噴射量、機関回転速度、及び冷却水温等に対する目標ＥＧＲ量の設定を変更した場合であっても、ＥＧＲ量の変化は排気酸素濃度に現われるから、上記パティキュレート生成定数Ａsft、及びパティキュレート酸化定数Ａsotを再び適合し直す必要がない。従って、第１実施形態によれば、このような場合にも適合工数を低減することができる。
【００９０】
（２．第２実施形態）
次に、本発明の第２実施形態に係るパティキュレート排出量推定装置を含む内燃機関の排気浄化装置について説明する。この排気浄化装置は、ＥＧＲ装置５０による排気還流の遅れを考慮して現時点の排気酸素濃度（数１に使用する排気酸素濃度）ＸO2の値の精度を向上することで、パティキュレート排出量ＰＭの推定精度を向上する点のみにおいて、第１実施形態の排気浄化装置と異なっている。従って、以下、かかる相違点を中心に説明する。
【００９１】
＜第２実施形態のパティキュレート排出量推定方法＞
排気酸素濃度ＸO2は、エンジン１０の運転状態に応じて変化する。従って、エンジンの運転状態がある運転状態から別の運転状態へと急変すると、排気酸素濃度ＸO2も急変する。一方、排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rは、燃焼室から排出された直後の排気酸素濃度ＸO2に対して時間的に遅れた挙動を示す。
【００９２】
他方、排気還流管５１を介して還流される排気が再び燃焼室に流入するまでには所定の時間を要する。換言すると、エンジン１０の運転状態の急変に伴う排気酸素濃度ＸO2の急変が、排気還流によって新たな排気酸素濃度ＸO2の変化をもたらすまでには所定の時間（以下、この時間を「排気還流遅れ時間」と称呼する。）を要する。
【００９３】
そこで、本実施形態は、かかる排気還流遅れ時間Ｔｄを計測・学習しておき、エンジンの運転状態が変化して排気酸素濃度ＸO2が急変した場合（例えば、これは、排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rが急変することで検出できる。）には、同運転状態の変化から排気還流遅れ時間Ｔｄの経過後に、同運転状態の変化に伴って酸素濃度が変化した排気が燃焼室に吸気され、その時点で排気酸素濃度が更に急変するとの考えに基いて同排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rを補正することにより、上記数１にて使用する排気酸素濃度ＸO2の精度を高める。
【００９４】
ここで、排気還流遅れ時間Ｔｄの計測方法について図４（Ａ）を参照しながら述べる。図４（Ａ）において、実線は排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rを示し、破線は排気還流遅れ検出用酸素濃度センサ７９が検出する還流排気中酸素濃度ＸO2INRを示している。
【００９５】
エンジン１０の運転状態が時刻ｔ１にて急変すると、燃焼室から排出された直後の排気酸素濃度ＸO2は直ちに変化を開始し、排気酸素濃度センサ７７の検出する排気酸素濃度ＸO2は時刻ｔ１よりも遅れた時刻ｔ２にて変化を開始する。また、排気還流遅れ検出用酸素濃度センサ７９が検出する還流排気中酸素濃度ＸO2INRは時刻ｔ２よりも遅い時刻ｔ３から変化を開始する。従って、時刻ｔ２と時刻ｔ３を知ることにより、排気還流遅れ時間Ｔｄ（＝ｔ３−ｔ２）を求めることができる。
【００９６】
本実施形態は、このようにして排気還流遅れ時間Ｔｄを求め、更に、図４（Ｂ）に示したように、エンジン１０の運転状態の急変に伴い排気酸素濃度センサ７７によって検出される排気酸素濃度ＸO2Rが時刻ｔ４にて変化を開始した場合、同時刻ｔ４から排気還流遅れ時間Ｔｄの経過後である時刻ｔ５には時刻ｔ４において検出された排気酸素濃度ＸO2の変化が現われているはずであると考え、同時刻ｔ５以降において、実線にて示した排気酸素濃度センサ７７の検出する排気酸素濃度ＸO2を一点鎖線で示したように補正する。
【００９７】
この補正は、初期値を排気酸素濃度補正初期値ｋ１とし、その後、時間経過とともに「０」に向けて減衰する補正量を求め、同補正量を排気酸素濃度センサ７７の検出する排気酸素濃度ＸO2Rに加えることにより行う。なお、排気酸素濃度補正初期値ｋ１は、定常運転状態から過渡運転状態に変化した時点の燃料噴射量Ｑfinと機関回転速度NEとに基いて定められ、正の所定値又は負の所定値をとる。
【００９８】
＜第２実施形態の作動＞
次に、第２実施形態に係るパティキュレート排出量推定装置を含む排気浄化装置の作動について説明する。本装置のＣＰＵ６１は図４〜図６に示した各ルーチンを追加的に実行する点を除き、同第１実施形態のＣＰＵ６１と同一の処理を行うようになっている。従って、以下においては、かかる相違点を中心として説明する。
【００９９】
先ず、排気還流遅れ時間Ｔｄを求めるための作動から説明すると、ＣＰＵ６１は図５に示した排気還流遅れ時間学習ルーチンを所定時間の経過毎に実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ５００から処理を開始してステップ５０５に進み、学習完了フラグＧの値が「０」であるか否かを判定する。この学習完了フラグＧの値は、図示しないイグニッションスイッチが「オフ」状態から「オン」状態へと変更されたときに図示しないイニシャルルーチンにて「０」に設定されるようになっている。
【０１００】
従って、いま、イグニッションスイッチが「オフ」状態から「オン」状態へと変更された直後であるとして説明を続けると、学習完了フラグＧの値は「０」となっている。このため、ＣＰＵ６１はステップ５０５に「Ｙｅｓ」と判定してステップ５１０に進み、排気還流遅れ時間計測中フラグＦの値が「１」であるか否かを判定する。このフラグＦの値は、上述したイニシャルルーチンにて「０」に設定されるようになっている。
【０１０１】
このため、ＣＰＵ６１はステップ５１０にて「Ｙｅｓ」と判定してステップ５１５に進み、排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rと、所定時間前の同排気酸素濃度ＸO2ROLDとの差の絶対値（即ち、排気酸素濃度ＸO2Rの所定時間内の変化量）が所定の閾値ＴＨ１より大きいか否かを判定する。この場合、エンジン１０の運転状態が定常状態にあれば、排気酸素濃度ＸO2Rの所定時間内の変化量は小さいから、ＣＰＵ６１はステップ５１５にて「Ｎｏ」と判定してステップ５９５に進み、本ルーチンを一旦終了する。
【０１０２】
以降、ＣＰＵ６１は本ルーチンを所定時間の経過毎に繰り返す。従って、エンジン１０の運転状態が変化して排気酸素濃度ＸO2Rの所定時間内の変化量が閾値ＴＨ１より大きくなると、ＣＰＵ６１はステップ５１５に進んだとき「Ｙｅｓ」と判定してステップ５２０に進み、排気還流遅れ時間計測中フラグＦの値を「１」に設定するとともに、カウンタＣ１の値を「０」に設定し、本ルーチンを一旦終了する。
【０１０３】
そして、所定の時間が経過してＣＰＵ６１が本ルーチンを再び実行すると、同ＣＰＵ６１はステップ５０５にて「Ｙｅｓ」と判定するとともに、排気還流遅れ時間計測中フラグＦの値が「１」に変更されているから、ステップ５１０にて「Ｎｏ」と判定してステップ５３０に進み、同ステップ５３０にて排気還流遅れ検出用酸素濃度センサ７９が検出する還流排気中酸素濃度ＸO2INRと、所定時間前の同還流排気中酸素濃度ＸO2INROLDとの差の絶対値（即ち、還流排気中酸素濃度ＸO2INRの所定時間内の変化量）が所定の閾値ＴＨ２より大きいか否かを判定する。
【０１０４】
現時点においては、排気酸素濃度ＸO2の所定時間内の変化量が閾値ＴＨ１より大きくなった直後であるから、その排気酸素濃度ＸO2Rが変化した排気が排気還流遅れ検出用酸素濃度センサ７９の配設位置にまで到達していない。従って、還流排気中酸素濃度ＸO2INRの所定時間内の変化量は閾値ＴＨ２よりも小さいので、ＣＰＵ６１はステップ５３０にて「Ｎｏ」と判定してステップ５３５に進み、カウンタＣ１の値を「１」だけ増大した後、ステップ５９５にて本ルーチンを一旦終了する。
【０１０５】
以降、ＣＰＵ６１は、還流排気中酸素濃度ＸO2INRの所定時間内の変化量が閾値ＴＨ２より大きくなるまで、ステップ５００〜ステップ５１０、ステップ５３０、及びステップ５３５を繰り返し実行するので、カウンタＣ１の値は所定時間の経過毎に「１」だけ増大して行く。そして、所定の時間が経過して還流排気中酸素濃度ＸO2INRの所定時間内の変化量が閾値ＴＨ２より大きくなると、ＣＰＵ６１はステップ５３０に進んだとき、同ステップ５３０にて「Ｙｅｓ」と判定してステップ５４０に進み、排気還流遅れ時間Ｔｄとして上記カウンタＣ１の値を格納し、続くステップ５４５にて学習完了フラグＧの値を「１」に設定した後、ステップ５９５にて本ルーチンを一旦終了する。
【０１０６】
このように、ＣＰＵ６１は、排気酸素濃度ＸO2Rの所定時間内の変化量が閾値ＴＨ１より大きくなってから、還流排気中酸素濃度ＸO2INRの所定時間内の変化量が閾値ＴＨ２より大きくなるまでの時間をカウンタＣ１により計測し、同計測した時間を排気還流遅れ時間Ｔｄとして学習する。
【０１０７】
また、ＣＰＵ６１は、図６に示した排気酸素濃度補正初期値ｋ１決定ルーチンを所定時間TSAMPLEの経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ６００から処理を開始してステップ６０５に進み、燃料噴射量Ｑfinと所定時間前の燃料噴射量Ｑfinoldとの差の絶対値が所定の閾値ＴＨ３より大きくなったか否かを判定する。このステップは、エンジン１０の運転状態が過渡状態となったか否かを判定するために設けられたステップである。
【０１０８】
いま、エンジン１０の運転状態が過渡運転状態でなく定常運転状態であるとして説明を続けると、この場合燃料噴射量Ｑfinと所定時間前の燃料噴射量Ｑfinoldとの差の絶対値は所定の閾値ＴＨ３より小さいので、ＣＰＵ６１はステップ６０５にて「Ｎｏ」と判定してステップ６９５に進み、本ルーチンを一旦終了する。このように、エンジン１０が定常運転状態にあるとき、排気酸素濃度補正初期値ｋ１は更新されない。
【０１０９】
その後、エンジン１０の運転状態が過渡運転状態となると、燃料噴射量Ｑfinと所定時間前の燃料噴射量Ｑfinoldとの差の絶対値は所定の閾値ＴＨ３より大きくなる。従って、ＣＰＵ６１はステップ６０５に進んだとき、同ステップ６０５にて「Ｙｅｓ」と判定してステップ６１０に進み、過渡運転状態フラグＨ１の値を「１」に設定する。この過渡運転状態フラグＨ１は、その値が「１」のときエンジン１０の運転状態が過渡運転状態となったことを示し、図示しない前述したイニシャルルーチンにより「０」に設定されるようになっている。
【０１１０】
次いで、ＣＰＵ６１はステップ６１５に進み、その時点の燃料噴射量Ｑfinと機関回転速度NEと予めＲＯＭ６２に記憶されているテーブルとに基いて、同燃料噴射量Ｑfinと同機関回転速度NEとが維持された場合の排気酸素濃度を予測排気酸素濃度ＳＸO2Rとして予測し、続くステップ６２０にて排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rが前記予測排気酸素濃度ＳＸO2Rより大きいか否かを判定する。
【０１１１】
そして、ＣＰＵ６１は、予測排気酸素濃度ＳＸO2Rが現時点の排気酸素濃度センサ７７の検出する排気酸素濃度ＸO2Rよりも大きいとき、ステップ６２０にて「Ｙｅｓ」と判定してステップ６２５に進み、排気酸素濃度補正初期値ｋ１を所定値Ａ１の絶対値に設定し、ステップ６９５にて本ルーチンを一旦終了する。これに対し、予測排気酸素濃度ＳＸO2Rが現時点の排気酸素濃度センサ７７の検出する排気酸素濃度ＸO2Rよりも小さいとき、ステップ６２０にて「Ｎｏ」と判定してステップ６３０に進み、排気酸素濃度補正初期値ｋ１を所定値Ａ１の絶対値の符号を負とした値に設定し、ステップ６９５にて本ルーチンを一旦終了する。
【０１１２】
このように、ＣＰＵ６１は、エンジン１０の運転状態が変化した後の排気酸素濃度が如何なる濃度になるかを予測排気酸素濃度ＳＸO2Rとして予測し、現時点の排気酸素濃度ＸO2R（即ち、エンジン１０の運転状態が変化する前の運転状態により得られた排気酸素濃度ＸO2R）との大小比較を行う（ステップ６２０）ことにより排気酸素濃度補正初期値ｋ１を正の値とすべきか負の値にすべきかを決定し、同大小比較に応じた符号を有する同排気酸素濃度補正初期値ｋ１を決定する。即ち、エンジン１０の運転状態の変化により排気酸素濃度ＸO2Rが増大すると予測される場合は排気酸素濃度補正初期値ｋ１を正の値（|Ａ|）とし、排気酸素濃度ＸO2Rが減少するであろうと予測される場合は、排気酸素濃度補正初期値ｋ１を負の値（−|Ａ|）とする。
【０１１３】
一方、ＣＰＵ６１は、図３に示したルーチンを所定時間の経過毎に実行しているから、所定のタイミングとなるとステップ３００〜３１０を実行して生成定数Ａsftと酸化定数Ａsotを決定し、ステップ３１５に進んで排気酸素濃度ＸO2を取得するため、酸素濃度ＸO2に排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rを設定してから、図７に示した酸素濃度ＸO2算出ルーチンの処理をステップ７００より開始する。
【０１１４】
次いで、ＣＰＵ６１はステップ７０５に進み、補正必要フラグＨ２の値が「１」であるか否かを判定する。この補正必要フラグＨ２は、その値が「１」のとき排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rが補正される必要があることを示すとともに、その値が「０」のとき排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rが補正される必要がないことを示す。なお、補正必要フラグＨ２の値は、図示しない前述したイニシャルルーチンにより「０」に設定されるようになっている。
【０１１５】
いま、エンジン１０の運転状態が定常運転状態であるとして説明を続けると、このような場合、排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rは略一定の値を維持するから、パティキュレート排出量推定のための排気酸素濃度ＸO2として使用することができる。従って、補正必要フラグＨ２の値は「０」に維持されているので（後述するステップ７５５を参照）、ＣＰＵ６１はステップ７０５にて「Ｎｏ」と判定してステップ７１０に進み、過渡運転状態フラグＨ１の値が「１」であるか否かを判定する。
【０１１６】
図６に示したルーチンのステップ６０５及びステップ６１０にて説明したように、過渡運転状態フラグＨ１の値は、エンジン１０が定常運転状態にある場合は「０」に維持されている。従って、ＣＰＵ６１はステップ７１０にて「Ｎｏ」と判定してステップ７９５に進み、本ルーチンを一旦終了する。この結果、排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rがパティキュレート排出量推定用の酸素濃度ＸO2となる。
【０１１７】
次に、エンジン１０の運転状態が定常運転状態から過渡運転状態へと変化した場合について説明する。この場合、燃料噴射量Ｑfinが急変するから、ＣＰＵ６１は図６のステップ６０５を実行するとき、同ステップ６０５にて「Ｙｅｓ」と判定し、ステップ６１０に進んで過渡運転状態フラグＨ１の値を「１」に設定する。ただし、この段階では、補正必要フラグＨ２の値は「０」に維持されている。
【０１１８】
この結果、ＣＰＵ６１が図３のルーチンのステップ３１５の実行に際し図７のルーチンのステップ７０５に至ると、同ステップ７０５にて「Ｎｏ」と判定してステップ７１０に進み、同ステップ７１０にて「Ｙｅｓ」と判定してステップ７１５に進んで排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rと所定時間前の排気酸素濃度ＸO2ROLDとの差の絶対値が所定の閾値ＴＨ４より大きくなったか否かを判定する。換言すると、ステップ７１５において、ＣＰＵ６１は、エンジン１０の運転状態の変化に伴う酸素濃度の急変が排気酸素濃度センサ７７により検出されたか否かを判定する。
【０１１９】
現時点は、エンジン１０の運転状態が過渡運転状態に変化した直後であるから、同エンジン１０の燃焼室から排出された直後の排気中の排気酸素濃度は変化しても、同変化は排気酸素濃度センサ７７により検出されない。従って、ＣＰＵ６１はステップ７１５にて「Ｎｏ」と判定して本ルーチンを一旦終了する。
【０１２０】
その後、所定の時間が経過すると、運転状態の変化に伴う排気酸素濃度の変化が排気酸素濃度センサ７７により検出される。このとき、ＣＰＵ６１がステップ７１５に至ると、同ＣＰＵ６１は同ステップ７１５にて「Ｙｅｓ」と判定してステップ７２０に進み、カウンタＣ２の値を「０」に設定するとともに、続くステップ７２５にて補正必要フラグＨ２の値を「１」に設定し、その後ステップ７９５にて本ルーチンを一旦終了する。
【０１２１】
その後、ＣＰＵ６１が図３に示したパティキュレート排出量推定ルーチンを実行すると、ステップ３１５の処理に伴うステップ７０５の処理において「Ｙｅｓ」と判定してステップ７３０に進み、同ステップ７３０にてカウンタＣ２の値を「１」だけ増大し、続くステップ７３５にてカウンタＣ２の値が前述した排気還流遅れ時間Ｔｄ以上となったか否かを判定する。この時点では、ＣＰＵ６１はステップ７３５にて「Ｎｏ」と判定してステップ７９５に進み本ルーチンを一旦終了する。この結果、ステップ３１５にて取得される排気酸素濃度ＸO2は依然として排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rに維持される。
【０１２２】
その後、ＣＰＵ６１はステップ７３０を繰り返し実行するようになるから、カウンタＣ２の値は次第に増大し、排気還流遅れ時間Ｔｄが経過するとカウンタＣ２の値が排気還流遅れ時間Ｔｄと一致する。従って、ＣＰＵ６１はステップ７３５にて「Ｙｅｓ」と判定してステップ７４０に進み、下記数１２に従って排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rを補正する。ここで、ｋ１は上述した図６に示したルーチンにより決定される排気酸素濃度補正初期値ｋ１である。
【０１２３】
【数１２】
ＸO2＝ＸO2R＋ｋ１／exp（Ｃ２−Ｔ）
【０１２４】
次いで、ＣＰＵ６１はステップ７４５にて、上記数１２の第２項である排気酸素濃度補正値ｋ１／exp（Ｃ２−Ｔ）の絶対値が所定の閾値ＴＨ５より小さくなったか否かを判定する。閾値ＴＨ５は、補正初期値ｋ１よりも小さい値に選択されている。現時点においてはＣ２＝Ｔであるから、排気酸素濃度補正値ｋ１／exp（Ｃ２−Ｔ）は初期値ｋ１と等しく、閾値ＴＨ５より大きいので、ＣＰＵ６１はステップ７４５にて「Ｎｏ」と判定してステップ７９５に進み、本ルーチンを一旦終了する。この結果、ＣＰＵ６１が図３に示したルーチンのステップ３２０にて使用する酸素濃度ＸO2は、排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rに初期値ｋ１を加えた値となる。
【０１２５】
以降、このような処理が継続されてカウンタＣ２の値がステップ７３０にて増大し続けるので、排気酸素濃度補正値ｋ１／exp（Ｃ２−Ｔ）の絶対値は次第に減少するとともに、ステップ７４０にて求められる酸素濃度ＸO2がパティキュレート排出量ＰＭの推定に用いられる。そして、排気酸素濃度補正値ｋ１／exp（Ｃ２−Ｔ）の絶対値が所定の閾値ＴＨ５より小さくなると、ＣＰＵ６１はステップ７４５にて「Ｙｅｓ」と判定してステップ７５０に進み、過渡運転状態フラグＨ１の値を「０」に設定するとともに、続くステップ７５５にて補正必要フラグＨ２の値を「０」に設定して、ステップ７９５にて本ルーチンを一旦終了する。
【０１２６】
この結果、ＣＰＵ６１が再び図７に示したルーチンを実行すると、同ＣＰＵ６１はステップ７０５にて「Ｎｏ」、ステップ７１０にて「Ｎｏ」と判定するようになるので、排気酸素濃度センサ７７が検出する排気酸素濃度ＸO2Rが補正されることなくパティキュレート排出量ＰＭを推定するための排気酸素濃度ＸO2として使用されるようになる。
【０１２７】
以上、説明したように、第２実施形態のパティキュレート排出量推定装置は、パティキュレート排出量ＰＭの推定に用いる排気酸素濃度ＸO2を、燃焼室から排出された直後の排気酸素濃度により近い値となるように、排気還流遅れ時間を考慮して求める。従って、パティキュレート排出量ＰＭの推定精度がより一層向上する。
【０１２８】
（３．第３実施形態）
次に、本発明の第３実施形態に係るパティキュレート排出量推定装置を含む内燃機関の排気浄化装置について説明する。この排気浄化装置は、吸気酸素濃度を考慮してパティキュレート排出量ＰＭの推定を行うことで、同パティキュレート排出量ＰＭの推定精度を向上する点、具体的には、ＣＰＵ６１が図３に示したルーチンに代えて図８に示したパティキュレート排出量推定ルーチンを所定時間TSAMPLEの経過毎に実行する点においてのみ、第１実施形態の排気浄化装置と異なっている。従って、以下、かかる相違点を中心に説明する。なお、図８において、図３と同一のステップには同一の符号を付し、その詳細説明を省略する。
【０１２９】
ＣＰＵ６１は、所定のタイミングになると、ステップ８００から処理を開始し、ステップ３０５からステップ３１５を実行することで、生成定数Ａsft、及び酸化定数Ａsotを決定するとともに、排気酸素濃度ＸO2を取得する。なお、排気酸素濃度ＸO2は、上記第２実施形態のように、排気還流遅れ時間Ｔｄを考慮して取得されるように構成してもよい。
【０１３０】
次いで、ＣＰＵ６１はステップ８０５に進み、吸気酸素濃度センサ７３が検出する吸気酸素濃度ＸO2INを読み込み、続くステップ８１０にて同吸気酸素濃度ＸO2INに応じてＲＯＭ６２内に格納されている吸気酸素濃度補正係数テーブルに基いて吸気酸素濃度補正係数ｆXO2を算出する。吸気酸素濃度補正係数テーブル内の吸気酸素濃度補正係数ｆXO2は実験的に適合されている値である。そして、ＣＰＵ６１は、下記数１３に基づいてパティキュレート排出量ＰＭを推定する。
【０１３１】
【数１３】
ＰＭ=Ａsft・(Ｋ・Ｑfin・2・NE)・ｆXO2/(Ａsot・ＸO2)+Ｃ/exp(Ａsot・ＸO2）
【０１３２】
次いで、ＣＰＵ６１はステップ３２５にてパティキュレート排出量合計値SUMPMを求め、ステップ８９５にて本ルーチンを一旦終了する。以降、ＣＰＵ６１は本ルーチンを所定時間TSAMPLEの経過毎に繰り返し実行する。この結果、パティキュレート排出量ＰＭとパティキュレート排出量合計値SUMPMが更新されて行く。
【０１３３】
このように、第３実施形態によれば、吸気酸素濃度ＸO2INにも基いてパティキュレート排出量ＰＭが推定される。吸気酸素濃度ＸO2INは燃焼温度に大きな影響を及ぼすから、パティキュレート排出量ＰＭも吸気酸素濃度ＸO2INによる影響を受ける。従って、この第３実施形態のように、吸気酸素濃度補正係数ｆXO2を導入してパティキュレート排出量ＰＭを推定するようにすれば、同パティキュレート排出量ＰＭの推定精度を向上することができる。なお、吸気酸素濃度ＸO2INに応じて第２吸気酸素濃度補正係数ｆ２XO2を求め、数１３の右辺第２項にも乗じるようにしてもよい。
【０１３４】
（４．第４実施形態）
次に、本発明の第４実施形態に係るパティキュレート排出量推定装置を含む内燃機関の排気浄化装置について説明する。この排気浄化装置は、燃焼室壁面温度に基いて推定された燃焼温度を考慮してパティキュレート排出量ＰＭの推定を行うことで、同パティキュレート排出量ＰＭの推定精度を向上する点、具体的には、ＣＰＵ６１が図３に示したルーチンに代えて図９に示したパティキュレート排出量推定ルーチンを所定時間TSAMPLEの経過毎に実行する点においてのみ、第１実施形態の排気浄化装置と異なっている。従って、以下、かかる相違点を中心に説明する。なお、図９において、図３と同一のステップには同一の符号を付し、その詳細説明を省略する。
【０１３５】
ＣＰＵ６１は、所定のタイミングになると、ステップ９００から処理を開始し、ステップ３０５からステップ３１５を実行することで、生成定数Ａsft、及び酸化定数Ａsotを決定するとともに、排気酸素濃度ＸO2を取得する。なお、排気酸素濃度ＸO2は、上記第２実施形態のように、排気還流遅れ時間Ｔｄを考慮して取得されるように構成してもよい。
【０１３６】
次いで、ＣＰＵ６１はステップ９０５に進み、壁面温度センサ７５が検出する燃焼室壁面温度Ｔhekiを読み込み、続くステップ９１０にて燃焼室壁面温度Ｔhekiに熱伝達効率ηを乗じて燃焼平均温度Ｔnenを求める。そして、ＣＰＵ６１はステップ９１５にて燃焼平均温度Ｔnenに基く補正係数ｇ（Ｔnen）を求めるとともに、これをステップ３０５にて求めた生成定数Ａsftに乗じることにより補正後の生成定数Ａsftを求め、同様に、燃焼平均温度Ｔnenに基く補正係数ｈ（Ｔnen）を求めるとともに、これをステップ３１０にて求めた酸化定数Ａsotに乗じることにより補正後の酸化定数Ａsotを求め、続くステップ３２０にてこれら補正後の定数を用いてパティキュレート排出量ＰＭを推定し、ステップ３２５にてパティキュレート排出量合計値SUMPMを求めた後、ステップ９９５にて本ルーチンを一旦終了する。なお、補正係数ｇ（Ｔnen）、及びｈ（Ｔnen）を求めるための関数ｇ，ｈは予め実験により定められていて、テーブルとしてＲＯＭ６２内に格納されている。
【０１３７】
以降、ＣＰＵ６１は図９に示したルーチンを所定時間TSAMPLEの経過毎に繰り返し実行する。この結果、パティキュレート排出量ＰＭとパティキュレート排出量合計値SUMPMが更新されて行く。
【０１３８】
以上、説明したように、第４実施形態によれば、燃焼平均温度Ｔnenに基いて生成定数Ａsft、及び酸化定数Ａsotが修正され、この修正された生成定数Ａsft、及び酸化定数Ａsotに基いてパティキュレート排出量ＰＭが推定される。燃焼平均温度Ｔnenは、パティキュレート排出量ＰＭに影響を与える。従って、この第４実施形態のように、燃焼平均温度Ｔnenに応じてパティキュレート排出量ＰＭを推定するように構成すれば、同パティキュレート排出量ＰＭの推定精度を向上することができる。
【０１３９】
なお、第４実施形態においては、燃焼平均温度Ｔnenに基いて生成定数Ａsft、及び酸化定数Ａsotを修正することによりパティキュレート排出量ＰＭの推定精度を向上していたが、ステップ９１５を省略するとともに、ステップ３２０にて求めたパティキュレート排出量ＰＭを燃焼平均温度Ｔnenに応じて修正するように構成してもよい。
【０１４０】
（５．第５実施形態）
次に、本発明の第５実施形態に係るパティキュレート排出量推定装置を含む内燃機関の排気浄化装置について説明する。この排気浄化装置は、吸気温度を考慮してパティキュレート排出量ＰＭの推定を行うことで、同パティキュレート排出量ＰＭの推定精度を向上する点、具体的には、ＣＰＵ６１が図３に示したルーチンに代えて図１０に示したパティキュレート排出量推定ルーチンを所定時間TSAMPLEの経過毎に実行する点においてのみ、第１実施形態の排気浄化装置と異なっている。従って、以下、かかる相違点を中心に説明する。なお、図１０において、図３と同一のステップには同一の符号を付し、その詳細説明を省略する。
【０１４１】
ＣＰＵ６１は、所定のタイミングになると、ステップ１０００から処理を開始し、ステップ３０５からステップ３１５を実行することで、生成定数Ａsft、及び酸化定数Ａsotを決定するとともに、排気酸素濃度ＸO2を取得する。なお、排気酸素濃度ＸO2は、上記第２実施形態のように、排気還流遅れ時間Ｔｄを考慮して取得されるように構成してもよい。
【０１４２】
次いで、ＣＰＵ６１はステップ１００５に進み、吸気温センサ７２が検出する吸気温度Ｔｂを読み込み、続くステップ１０１０にて読み込んだ吸気温度Ｔｂに応じて吸気温度補正係数ｆTbを関数ｉ（Ｔｂ）に基いて求める。そして、ＣＰＵ６１はステップ１０１５にて、下記数１３に従ってパティキュレート排出量ＰＭを推定する。なお、関数ｉ（Ｔｂ）は予め実験により定められていて、テーブルとしてＲＯＭ６２内に格納されている。
【０１４３】
【数１４】
ＰＭ=ｆTb・(Ａsft・(Ｋ・Ｑfin・２・NE)/(Ａsot・ＸO2)＋Ｃ/exp(Ａsot・ＸO2）)
【０１４４】
その後、ＣＰＵ６１は、ステップ３２５にてパティキュレート排出量合計値SUMPMを求めた後、ステップ１０９５にて本ルーチンを一旦終了する。以降、ＣＰＵ６１は本ルーチンを所定時間TSAMPLEの経過毎に繰り返し実行する。この結果、パティキュレート排出量ＰＭとパティキュレート排出量合計値SUMPMが更新されて行く。
【０１４５】
以上、説明したように、第５実施形態によれば、吸気温度Ｔｂにも基いてパティキュレート排出量ＰＭが推定される。吸気温度Ｔｂは燃焼温度に大きな影響を及ぼすから、パティキュレート排出量ＰＭも吸気温度Ｔｂによる影響を受ける。従って、この第５実施形態のように、吸気温度補正係数ｆTbを導入してパティキュレート排出量ＰＭを推定するようにすれば、同パティキュレート排出量ＰＭの推定精度を向上することができる。
【０１４６】
（６．第６実施形態）
次に、本発明の第６実施形態に係るパティキュレート排出量推定装置を含む内燃機関の排気浄化装置について説明する。この排気浄化装置は、圧縮端温度（圧縮上死点近傍における燃焼室内ガス温度）と燃焼圧力（最高燃焼圧力）とを考慮してパティキュレート排出量ＰＭの推定を行うことで、同パティキュレート排出量ＰＭの推定精度を向上する点、具体的には、ＣＰＵ６１が図１１に示した圧縮端温度推定及び燃焼圧力検出ルーチンを実行する点、及び、図３に示したルーチンに代えて図１２に示したパティキュレート排出量推定ルーチンを所定時間TSAMPLEの経過毎に実行する点においてのみ、第１実施形態の排気浄化装置と異なっている。従って、以下、かかる相違点を中心に説明する。なお、図１２において、図３と同一のステップには同一の符号を付し、その詳細説明を省略する。
【０１４７】
先ず、圧縮端温度Ｔtdc、及び燃焼圧力Ｐｎを検出するための作動について説明すると、ＣＰＵ６１は、図１１に示したルーチンを所定クランク角度（例えば、クランク角度１°）の経過毎に実行するようになっている。従って、エンジン１０のクランク角度が所定のクランク角度になると、ＣＰＵ６１はステップ１１００から処理を開始し、ステップ１１０５に進んで現時点のクランク角が第３気筒の圧縮上死点か否かを判定する。
【０１４８】
そして、現時点のクランク角が第３気筒の圧縮上死点であれば、ＣＰＵ６１はステップ１１１０に進んでその時点にて燃焼圧センサ７４が検出している燃焼室内圧力Ｐａと関数ｑとに基いて推定圧縮端温度Ｔtdcを求め（圧縮端温度を推定し）、ステップ１１１５に進む。一方、現時点のクランク角が第３気筒の圧縮上死点でなければ、ＣＰＵ６１は直接ステップ１１１５に進む。
【０１４９】
次に、ＣＰＵ６１は、ステップ１１１５にて第３気筒が圧縮上死点となった後に、燃焼圧センサ７４が検出している燃焼室内圧力Ｐａが極大値をとったか否かを判定する。そして、燃焼室内圧力Ｐａが極大値をとっていれば、その極大値である燃焼室内圧力Ｐａを燃焼圧力Ｐｎ（燃焼中の燃焼室内圧力に関連した値）として設定し、ステップ１１９５に進んで本ルーチンを一旦終了する。一方、燃焼室内圧力Ｐａが極大値をとっていなければ、そのままステップ１１９５に進んで本ルーチンを一旦終了する。以上により、圧縮端温度Ｔtdcと燃焼圧力Ｐｎが求められる。
【０１５０】
また、ＣＰＵ６１は、所定のタイミングになると、図１２のステップ１２００から処理を開始し、ステップ３０５からステップ３１５を実行することで、生成定数Ａsft、及び酸化定数Ａsotを決定するとともに、排気酸素濃度ＸO2を取得する。なお、排気酸素濃度ＸO2は、上記第２実施形態のように、排気還流遅れ時間Ｔｄを考慮して取得されるように構成してもよい。
【０１５１】
次いで、ＣＰＵ６１は、ステップ１２０５及びステップ１２１０にて、前述した図１１に示したルーチンにより求められている圧縮端温度Ｔtdc、及び燃焼圧力Ｐｎをそれぞれ読み込み、続くステップ１２１５にて圧縮端温度Ｔtdc及び燃焼圧力Ｐｎにより定まる補正係数ｍ（Ｔtdc，Ｐｎ）をＲＯＭ６２内に予め格納されているテーブルｍと同圧縮端温度Ｔtdcと同燃焼圧力Ｐｎとから求め、この補正係数ｍ（Ｔtdc，Ｐｎ）をステップ３０５にて求めた生成定数Ａsftに乗じることにより補正後の生成定数Ａsftを求めるとともに、同様に、圧縮端温度Ｔtdc及び燃焼圧力Ｐｎにより定まる補正係数ｎ（Ｔtdc，Ｐｎ）をＲＯＭ６２内に予め格納されているテーブルｎと同圧縮端温度Ｔtdcと同燃焼圧力Ｐｎとから求め、この補正係数ｎ（Ｔtdc，Ｐｎ）をステップ３１０にて求めた酸化定数Ａsotに乗じることにより補正後の酸化定数Ａsotを求める。
【０１５２】
次いで、ＣＰＵ６１はステップ３２０にてパティキュレート排出量ＰＭを推定し、ステップ３２５にてパティキュレート排出量合計値SUMPMを求めた後、ステップ１２９５にて本ルーチンを一旦終了する。
【０１５３】
以降、ＣＰＵ６１は図１１に示したルーチンを所定クランク角度の経過毎に繰り返し実行するとともに、図１２に示したルーチンを所定時間TSAMPLEの経過毎に繰り返し実行する。この結果、パティキュレート排出量ＰＭとパティキュレート排出量合計値SUMPMが更新されて行く。
【０１５４】
以上、説明したように、第６実施形態によれば、燃焼圧センサ７４が検出する燃焼室内圧力Ｐａに基いて圧縮端温度Ｔtdcと燃焼圧力Ｐｎとが求められ、この圧縮端温度Ｔtdcと燃焼圧力Ｐｎとに基いて生成定数Ａsft、及び酸化定数Ａsotが修正され、この修正された生成定数Ａsft、及び酸化定数Ａsotに基いてパティキュレート排出量ＰＭが推定される。圧縮端温度Ｔtdcは燃焼温度に強く影響するから、パティキュレート排出量ＰＭも同圧縮端温度Ｔtdcによる影響を受ける。また、燃焼圧力Ｐｎはパティキュレート排出量ＰＭに強く影響する。従って、この第６実施形態のように、圧縮端温度Ｔtdcと燃焼圧力Ｐｎとに応じてパティキュレート排出量ＰＭを推定するように構成すれば、同パティキュレート排出量ＰＭの推定精度を向上することができる。
【０１５５】
なお、第６実施形態においては、圧縮端温度Ｔtdcと燃焼圧力Ｐｎとに基いて生成定数Ａsft、及び酸化定数Ａsotを修正することによりパティキュレート排出量ＰＭの推定精度を向上していたが、ステップ１２１５を省略し、ステップ３２０にて求めたパティキュレート排出量ＰＭを圧縮端温度Ｔtdcと燃焼圧力Ｐｎとに応じて修正するように構成してもよい。
【０１５６】
また、第６実施形態においては、所定の微小クランク角の経過毎に燃焼室内の圧力Ｐをサンプリングしていたが、これを所定時間毎にサンプリングするように構成してもよい。この場合、圧縮端温度Ｔtdcと燃焼圧力Ｐｎとを精度良く求められなければ、ステップ１２１５にて燃焼行程中における燃焼室内圧力Ｐａの時間平均値Ｐaveに応じて生成定数Ａsft、及び酸化定数Ａsotを修正するように構成したり、又は、ステップ１２１５を省略してステップ３２０にて求めたパティキュレート排出量ＰＭを同平均圧力Ｐaveに応じて修正するように構成してもよい。
【０１５７】
（７．第７実施形態）
次に、本発明の第７実施形態に係るパティキュレート排出量推定装置を含む内燃機関の排気浄化装置について説明する。この排気浄化装置は、例えば、エンジン１０の運転状態が定常運転状態にあって、光学式煤センサ７６が比較的精度良く煤量を測定しうる場合に、同光学式煤センサ７６の出力である煤量Ｓsootに基いて実際のパティキュレート排出量である基準パティキュレート排出量ＳＰＭを求め、この基準パティキュレート排出量ＳＰＭと上記数１（数１０）によって推定したパティキュレート基本排出量ＰＭ０との差を学習値DLPMとして学習し、他の運転状態において同学習値DLPMで上記数１（数１０）によって推定したパティキュレート基本排出量ＰＭ０を修正してパティキュレート排出量ＰＭを求める点においてのみ、第１実施形態の排気浄化装置と異なっている。従って、以下、かかる相違点を中心に説明する。
【０１５８】
この排気浄化装置のＣＰＵ６１は、図１３に示した学習値更新ルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ１３００から処理を開始してステップ１３０５に進み、燃料噴射量Ｑfinと所定時間前の燃料噴射量Ｑfinoldとの差の絶対値が所定の閾値ＴＨ６より小さい状態が所定の時間以上継続したか否かを判定する。即ち、ステップ１３０５では、エンジン１０の運転状態が定常運転状態にあって光学式煤センサ７６が比較的精度良く煤量Ｓsootを測定しうる状態となっているか否かが判定される。
【０１５９】
いま、燃料噴射量Ｑfinと所定時間前の燃料噴射量Ｑfinとの差の絶対値が所定の閾値ＴＨ６より小さい状態が所定の時間以上継続していないと仮定して説明を続けると、ＣＰＵ６１はステップ１３０５にて「Ｎｏ」と判定してステップ１３９５に進み、本ルーチンを一旦終了する。従って、この場合、学習値DLPMは更新されない。
【０１６０】
これに対し、燃料噴射量Ｑfinと所定時間前の燃料噴射量Ｑfinoldとの差の絶対値が所定の閾値ＴＨ６より小さい状態が所定の時間以上継続していると、ＣＰＵ６１はステップ１３０５に進んだとき、同ステップ１３０５にて「Ｙｅｓ」と判定してステップ１３１０に進み、光学式煤センサ７６が測定している煤量Ｓsootを読み込む。次いで、ＣＰＵ６１はステップ１３１５に進み、基準パティキュレート排出量ＳＰＭを、前記読み込んだ煤量Ｓsootと関数ｒとに基いて求める。なお、関数ｒは予め実験により求められていて、ＲＯＭ６２内にテーブルとして格納されている。
【０１６１】
次に、ＣＰＵ６１は、ステップ１３２０に進み、前記基準パティキュレート排出量ＳＰＭと、後述する図１４に示したステップ１４０５にて上記数１（数１０）に従って推定されているパティキュレート基本排出量ＰＭ０との差を学習値DLPMとして学習し、ステップ１３９５に進んで本ルーチンを一旦終了する。
【０１６２】
一方、ＣＰＵ６１は、所定時間TSAMPLEの経過毎に図１４に示したＰＭ推定ルーチンを繰り返し実行するようになっている。なお、図１４において、図３と同一のステップには同一の符号を付し、その詳細説明を省略する。
【０１６３】
従って、ＣＰＵ６１は所定のタイミングになると、図１４のステップ１４００から処理を開始し、ステップ３０５からステップ３１５を実行することで、生成定数Ａsft、及び酸化定数Ａsotを決定するとともに、排気酸素濃度ＸO2を取得する。なお、排気酸素濃度ＸO2は、上記第２実施形態のように、排気還流遅れ時間Ｔｄを考慮して取得されるように構成してもよい。
【０１６４】
次いで、ＣＰＵ６１は、ステップ１４０５にてパティキュレート基本排出量ＰＭ０を上記数１（数１０）に従って推定し、続くステップ１４１０にてパティキュレート基本排出量ＰＭ０に前記学習値DLPMを加えることで同パティキュレート基本排出量ＰＭ０を補正する。そして、ステップ３２５にてパティキュレート排出量合計値SUMPMを求めた後、ステップ１４９５にて本ルーチンを一旦終了する。
【０１６５】
以降、ＣＰＵ６１は図１３及び図１４に示したルーチンを繰り返し実行するので、学習値DLPMが更新されるとともに、同学習値DLPMによりパティキュレート基本排出量ＰＭ０が修正されることにより、パティキュレート排出量ＰＭが精度良く求められて行く。
【０１６６】
以上、説明したように、第７実施形態によれば、光学式煤センサ７６の出力である煤量Ｓsoot、即ち、実測煤量に基いて実際のパティキュレート排出量である基準パティキュレート排出量ＳＰＭを推定し、この基準パティキュレート排出量ＳＰＭと上記数１（数１０）のモデルにより推定されたパティキュレート基本排出量ＰＭ０との差を学習値DLPMとして学習し、この学習値DLPMによってパティキュレート基本排出量ＰＭ０が修正される。
【０１６７】
従って、数１（数１０）のパティキュレート排出量推定モデルが有する誤差が実測値に基いた値で修正されるので、パティキュレート排出量を精度良く推定することが可能となる。
【０１６８】
以上、説明したように、本発明の各実施形態によれば、パティキュレート排出量ＰＭがリアルタイムに、且つ、定常運転状態及び過渡運転状態において、精度良く推定され得る。なお、本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上述したフィルタ再生制御においては、添加燃料の供給のみでなく、燃料噴射弁２１からのメイン噴射の後に再度噴射を行うポスト噴射の噴射量を変更したり、排気絞り弁４６の開度を変更して過給機３５の過給圧を変更してもよい。
【図面の簡単な説明】
【図１】 本発明の第１実施形態に係るパティキュレート排出量推定装置を備えた排気浄化装置を４気筒内燃機関（ディーゼル機関）に適用したシステムの概略構成図である。
【図２】 パティキュレート排出量推定モデルのパティキュレート生成定数及びパティキュレート酸化定数を決定する方法を説明するためのグラフである。
【図３】 図１に示したＣＰＵが実行するパティキュレート排出量を推定するために実行するルーチンを示したフローチャートである。
【図４】 図４（Ａ）は本発明の第２実施形態に係るパティキュレート排出量推定装置による排気還流遅れ時間の計測方法について説明するために、排気酸素濃度センサが検出する排気酸素濃度を実線で、排気還流遅れ検出用酸素濃度センサが検出する還流排気中酸素濃度を破線で示したタイムチャートであり、図４（Ｂ）は排気酸素濃度センサが検出する排気酸素濃度を実線で、学習した排気還流遅れ時間を用いて補正した排気酸素濃度を一点鎖線で示したタイムチャートである。
【図５】 第２実施形態のパティキュレート排出量推定装置のＣＰＵが実行する排気還流遅れ時間学習ルーチンを示したフローチャートである。
【図６】 第２実施形態のパティキュレート排出量推定装置のＣＰＵが実行する排気酸素濃度補正初期値決定ルーチンを示したフローチャートである。
【図７】 第２実施形態のパティキュレート排出量推定装置のＣＰＵが実行する排気酸素濃度算出ルーチンを示したフローチャートである。
【図８】 第３実施形態のパティキュレート排出量推定装置のＣＰＵが実行するパティキュレート排出量推定ルーチンを示したフローチャートである。
【図９】 第４実施形態のパティキュレート排出量推定装置のＣＰＵが実行するパティキュレート排出量推定ルーチンを示したフローチャートである。
【図１０】 第５実施形態のパティキュレート排出量推定装置のＣＰＵが実行するパティキュレート排出量推定ルーチンを示したフローチャートである。
【図１１】 第６実施形態のパティキュレート排出量推定装置のＣＰＵが実行する圧縮端温度推定及び燃焼圧力検出ルーチンを示したフローチャートである。
【図１２】 第６実施形態のパティキュレート排出量推定装置のＣＰＵが実行するパティキュレート排出量推定ルーチンを示したフローチャートである。
【図１３】 第７実施形態のパティキュレート排出量推定装置のＣＰＵが実行する学習値更新ルーチンを示したフローチャートである。
【図１４】 第７実施形態のパティキュレート排出量推定装置のＣＰＵが実行するパティキュレート排出量推定ルーチンを示したフローチャートである。
【符号の説明】
２１…燃料噴射弁、３１…吸気マニホールド、３２…吸気管、３３…スロットル弁、３５…過給機、４１…排気マニホールド、４２…排気管、４３…ＤＰＮＲ、４４…添加燃料供給ノズル、５０…ＥＧＲ装置、５１…排気管流管、７２…吸気温センサ、７３…吸気酸素濃度センサ、７４…燃焼圧センサ、７５…壁面温度センサ、７６…光学式煤センサ、７７…排気酸素濃度センサ、７８…機関回転速度センサ、７９…排気還流遅れ検出用酸素濃度センサ、Ａsft…パティキュレート生成定数、Ａsot…パティキュレート酸化定数。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a particulate matter emission amount estimation device for a diesel engine that estimates the amount of particulate matter (particulate emission amount) discharged from the diesel engine.
[0002]
[Prior art]
2. Description of the Related Art Conventionally, an exhaust emission control device in which a particulate filter is disposed in an exhaust passage of a diesel engine is known. The particulate filter reacts with soot (graphite, soot) in the exhaust gas, particulate matter (SOF) soluble in organic materials mainly composed of unburned fuel and unburned lubricating oil, and sulfur in the fuel. Particulate matter consisting of particulate matter (sulfate) insoluble in organic materials produced by the above (these particulate matter is referred to as “particulate matter”, but hereinafter simply referred to as “particulate”) It is a filter that collects.
[0003]
By the way, the particulate filter may cause clogging due to the collected particulate, and the exhaust resistance of the engine is increased by the clogging. For this reason, it is necessary to estimate the particulate collection amount of the particulate filter and perform filter regeneration control such as burning the particulate in the filter so that the particulate collection amount does not exceed a predetermined amount. is there.
[0004]
A method for calculating (estimating) such a particulate collection amount is disclosed in, for example, Japanese Patent Laid-Open No. 2000-170521. In the disclosed method, first, the particulate discharge amount discharged from the diesel engine is estimated.
[0005]
More specifically, particulate discharge amounts with respect to various fuel injection amounts and engine speeds are measured in advance, and the results are stored as table 1. Further, since the particulate discharge amount varies depending on the exhaust gas recirculation amount (EGR amount), an EGR correction coefficient is obtained for each engine speed, fuel injection amount, and engine cooling water temperature, which are parameters for determining the exhaust gas recirculation amount. This is stored as Table 2.
[0006]
In actual operation, the fuel injection amount, the engine speed, and the basic particulate discharge amount are obtained from the table 1, and the fuel injection amount, the engine speed, the cooling water temperature, and the EGR correction from the table 2. A coefficient is obtained, and the particulate discharge is estimated by correcting the basic particulate discharge with the EGR correction coefficient.
[0007]
Further, in the method disclosed in the above publication, the collection efficiency of the particulate filter is obtained based on the engine speed and the like, and the particulate collection amount is estimated from the estimated particulate discharge amount and the same collection efficiency. .
[0008]
[Problems to be solved by the invention]
However, even if the fuel injection amount and the engine speed are the same, for example, if the combustion temperature and the combustion pressure are different, the particulate discharge amount will be different, so the particulate discharge basic amount with only the fuel injection amount and the engine speed as arguments. Therefore, there is a problem that the particulate discharge amount cannot be accurately estimated even when the EGR correction coefficient is correct.
[0009]
In the conventional method disclosed above, the particulate discharge amount is estimated on the assumption that the exhaust gas recirculation amount is uniquely determined according to the fuel injection amount, the engine rotational speed, and the engine coolant temperature. Since the engine does not take into account fluctuations in the exhaust gas recirculation amount when the engine is in a transient operation state (for example, the increase in exhaust gas recirculation amount due to the supercharging delay of the supercharger), the EGR correction coefficient should be obtained correctly. As a result, there is a problem that the particulate discharge amount cannot be accurately estimated.
[0010]
Accordingly, one of the objects of the present invention is to provide a particulate emission estimation device that can easily and accurately estimate the particulate emission discharged from a diesel engine.
[0011]
[Outline of the present invention]
  The present invention for achieving the above objectThe particulate emission estimation device according toThe exhaust oxygen concentration was changed in a steady operation state in which the fuel injection amount and engine speed of the diesel engine were kept constant, and the particulate discharge amount was measured for each of the changed different exhaust oxygen concentrations. Based on the data consisting of the changed exhaust oxygen concentrations and the same particulate emissions measured for each exhaust oxygen concentration, the fuel injection amount maintained at the same constant and the engine speed maintained at the same constant are determined. Fuel injection amount formed by performing the operation of determining the particulate generation constant and the particulate oxidation constant for the operation state for a plurality of operation states each determined by a predetermined fuel injection amount and a predetermined engine speed. And a generation constant table that defines the relationship between the combination of the engine speed and the particulate generation constant, and fuel A storage means for storing an oxidation constant table that defines the relationship between the combination of the injection quantity and the engine rotational speed and the particulate oxidation constant; a fuel injection quantity acquisition means for acquiring a value representing the actual fuel injection quantity; Engine rotational speed detecting means for detecting the engine rotational speed, exhaust oxygen concentration acquiring means for acquiring a value representing the actual exhaust oxygen concentration, a value representing the acquired fuel injection amount, and the detected engine rotational speed And a production constant and an oxidation constant used for estimating the particulate emission amount based on the production constant table and the oxidation constant table, and the production constant determined and the same Particulate emission estimating means for estimating the particulate emission based on the oxidation constant and the value representing the acquired exhaust oxygen concentration.Yeah.
[0012]
  MoreThe particulate emission estimation means is PM for particulate emission per unit time, Asft for the determined particulate generation constant, Asot for the determined particulate oxidation constant, and the acquired fuel injection amount.Value representingPM = Asft · Q / (Asot · XO2) + C / exp (Asot · Assume) where Q is the fuel injection amount per unit time based on the above, X02 is the value representing the acquired exhaust oxygen concentration, and C is the constant. XO2) is configured to estimate the particulate emissions based on the formulaIt is.
[0013]
As a result of various studies, the present inventor has determined that the particulate matter emission amount PM per unit time is PM = Asft · Q / (Asot · regardless of whether the diesel engine is in a steady operation state or a transient operation state. XO2) + C / exp (Asot · XO2) was found to change substantially according to the equation. Here, Asft and Asot are respectively a particulate generation constant and a particulate oxidation constant that are uniquely determined by the fuel injection amount per one stroke and the engine speed, and Q is the fuel injection amount per one stroke and the engine rotational speed. The fuel injection amount per unit time, which is an amount determined by the above, X02 is the exhaust oxygen concentration, and C is a constant.
[0014]
As described above, the particulate matter emission amount PM can be easily obtained as a function of the exhaust oxygen concentration XO2 if the particulate generation constant Asft and the particulate oxidation constant Asot can be specified.
[0015]
Therefore, first, the particulate emission estimation apparatus of the present invention includes a generation constant table that defines a relationship between a combination of the fuel injection amount and the engine rotational speed and a particulate generation constant, and a combination of the fuel injection amount and the engine rotational speed. An oxidation constant table defining the relationship with the particulate oxidation constant is created and stored in the storage means.
[0016]
These tables show the steady operation state in which the fuel injection amount of the diesel engine and the engine rotational speed are kept constant (in this case, since the engine rotational speed is constant, the fuel injection amount is the fuel injection amount Qfin per stroke). May be the fuel injection amount Q per unit time.) And the particulate emission amount was measured for each of the changed different exhaust oxygen concentrations. The particulate generation constant Asft and the particulate oxidation constant Asot are determined (identified or identified) so that the above-described equation passes through the data points consisting of the respective exhaust oxygen concentrations and the respective particulate emissions measured for the respective exhaust oxygen concentrations. Fitting), and such work is repeatedly performed for various fuel injection amounts and engine speeds.
[0017]
Then, the value representing the actual fuel injection amount is acquired and the actual engine rotation speed is detected, the value representing the acquired fuel injection amount, the detected engine rotation speed, the stored generation constant table, and the storage The particulate generation constant and the particulate oxidation constant used for estimating the actual particulate emission amount are determined based on the oxidation constant table, and the determined constant and the fuel injection amount per unit time are determined. The obtained exhaust oxygen concentration is substituted into the above formula to estimate the particulate emission amount. That is, the particulate emission amount is estimated based on at least the determined particulate generation constant, the determined particulate oxidation constant, and the acquired exhaust oxygen concentration.
[0018]
Thus, the particulate emission estimation device having the above-described features of the present invention uses the exhaust oxygen concentration as a variable when estimating the particulate emission, and the exhaust oxygen concentration depends on the combustion temperature, the combustion pressure, etc. In comparison, it is a variable that is relatively easy to ensure measurement accuracy, and the particulate generation constant Asft and the particulate oxidation constant Asot were actually measured when the fuel injection amount and the engine speed were kept constant. Since it is determined by the particulate emission amount, the combustion temperature, combustion pressure, combustion ratio, combustion type, shape of the injection valve, engine exhaust amount, etc. for each operating state determined by the fuel injection amount and the engine speed The values reflect all parameters that affect particulate emissions. As a result, the particulate discharge amount estimation apparatus of the present invention can accurately estimate the particulate discharge amount.
[0019]
In addition, the particulate emission estimation device of the present invention estimates the particulate emission based on the value representing the actual exhaust oxygen concentration that varies depending on the variation of the exhaust gas recirculation amount (EGR amount). Even when the state becomes a transient operation state and the exhaust gas recirculation amount cannot be obtained as scheduled, the particulate discharge amount can be accurately obtained.
[0020]
In such a particulate emission estimation apparatus, the particulate emission estimation means estimates the particulate emission only when the acquired value representing the exhaust oxygen concentration is a value equal to or higher than a predetermined concentration. Suitably configured to do.
[0021]
The above formula is effective only when the air-fuel ratio of the gas sucked into the diesel engine is an air-fuel ratio greater than about 20, for example, and the combustion temperature is lowered by setting the air-fuel ratio A / F to 20 or less. Therefore, it does not correspond to combustion (so-called low temperature combustion) for suppressing particulate emission. Accordingly, as in the above configuration, the particulate discharge amount is estimated only when the acquired value representing the exhaust oxygen concentration is a value equal to or higher than a predetermined concentration, thereby estimating the particulate discharge amount. A decrease in accuracy can be avoided.
[0022]
By the way, in general, the exhaust oxygen concentration detected by the exhaust oxygen concentration detection means changes with a time delay with respect to the exhaust oxygen concentration of the engine (exhaust oxygen concentration to be used for estimating the particulate emission amount). Further, when the diesel engine is provided with exhaust gas recirculation means (EGR device) for performing exhaust gas recirculation through the exhaust gas recirculation pipe and controlling the amount of the exhaust gas recirculation, the exhaust oxygen concentration changes based on the change in the operating state. The exhausted air is sucked into the engine with a delay of the exhaust gas recirculation delay time. Therefore, at the time when the exhaust gas recirculation delay time has elapsed since the detection of the exhaust oxygen concentration change accompanying the change in the operating state, the engine exhaust oxygen concentration (exhaust oxygen concentration to be used for estimating the particulate emission amount) further changes. Should have started.
[0023]
Accordingly, in such a particulate emission estimation device, when the diesel engine includes exhaust gas recirculation means for performing exhaust gas recirculation through the exhaust gas recirculation pipe and controlling the exhaust gas recirculation amount, The concentration acquisition means includes exhaust oxygen concentration detection means for detecting the oxygen concentration in the exhaust passage of the diesel engine, and exhaust gas from the engine is recirculated by the exhaust gas recirculation means provided in the diesel engine and sucked into the engine. It is preferable that the detected oxygen concentration in the exhaust passage is corrected according to the exhaust gas recirculation delay time required to obtain a value representing the exhaust oxygen concentration.
[0024]
According to this, when the change in the exhaust oxygen concentration accompanying the change in the operating state causes a further change in the exhaust oxygen concentration due to the exhaust gas recirculation, an appropriate time point before the exhaust oxygen concentration detecting means detects the further change Therefore, it can be reflected in the exhaust oxygen concentration used for estimating the particulate emission amount. As a result, the exhaust oxygen concentration for estimating the particulate emission amount becomes an appropriate value, so that the estimation accuracy of the particulate emission amount is improved. be able to.
[0025]
In this case, the exhaust oxygen concentration acquisition means is for exhaust gas recirculation delay detection for detecting the oxygen concentration in the exhaust recirculation pipe immediately before the recirculated exhaust gas and fresh air newly sucked into the diesel engine are mixed. It is preferable that the exhaust gas recirculation delay time is estimated based on the detected oxygen concentration in the exhaust passage and the detected oxygen concentration in the exhaust gas recirculation pipe. is there.
[0026]
According to this, the detected oxygen concentration in the exhaust passage and the detected oxygen concentration in the exhaust recirculation pipe immediately before the exhaust gas to be recirculated and the fresh air newly sucked into the diesel engine are mixed. Since the exhaust gas recirculation delay time is estimated based on this, the exhaust gas recirculation delay time can be accurately determined based on the actual measurement value, and as a result, the exhaust gas oxygen concentration for estimating the particulate emission amount should be an appropriate value. Therefore, the estimation accuracy of the particulate discharge amount can be improved.
[0027]
The particulate emission estimation means is configured to acquire a value related to the combustion temperature of the diesel engine and to estimate the particulate emission based on a value related to the acquired combustion temperature. It is preferred that
[0028]
Since the combustion temperature has a relatively strong influence on the particulate discharge amount, the particulate discharge amount can be estimated with higher accuracy by using the combustion temperature for the estimation of the particulate discharge amount.
[0029]
The particulate emission estimation means includes oxygen concentration detection means for detecting the oxygen concentration in the intake passage of the diesel engine, and the detected oxygen concentration in the intake passage is set as a value related to the combustion temperature. Suitably configured to handle. Since the oxygen concentration in the intake passage is closely related to the combustion temperature, it is possible to accurately estimate the particulate emission amount by treating the oxygen concentration in the intake passage as a value related to the combustion temperature. .
[0030]
The particulate emission estimation means includes combustion chamber wall surface temperature detection means for detecting the temperature of the combustion chamber wall surface of the diesel engine, and treats the detected combustion chamber wall surface temperature as a value related to the combustion temperature. It is preferable to be configured as described above. The combustion chamber wall temperature is closely related to the combustion temperature so that, for example, the combustion average temperature can be obtained by multiplying the combustion chamber wall temperature by the heat transfer efficiency, so the combustion chamber wall temperature is a value related to the combustion temperature. It is possible to accurately estimate the particulate discharge amount.
[0031]
Further, the particulate emission amount estimation means includes intake air temperature detection means for detecting an intake air temperature of the diesel engine, and is configured to handle the detected intake air temperature as a value related to the combustion temperature. Is preferred. Since the intake air temperature is closely related to the combustion temperature, it is possible to accurately estimate the particulate emission amount by treating the intake air temperature as a value related to the combustion temperature.
[0032]
The particulate emission estimation means includes combustion chamber pressure detection means for detecting a pressure in the combustion chamber of the diesel engine, and is configured to handle the detected combustion chamber pressure as a value related to the combustion temperature. It is preferred that For example, the pressure in the combustion chamber at the compression top dead center is closely related to the combustion temperature. By treating the pressure in the combustion chamber as a value related to the combustion temperature, it is possible to accurately estimate the particulate emissions. It becomes.
[0033]
Further, in such a particulate emission estimation device, the particulate emission estimation means acquires a value related to the pressure in the combustion chamber during combustion of the diesel engine, and the acquired pressure in the combustion chamber Preferably, the particulate emission is estimated based on a value related to.
[0034]
The pressure in the combustion chamber during combustion, particularly the maximum value of the pressure in the combustion chamber called combustion pressure, has a relatively strong effect on particulate emissions. By using it for estimation, it is possible to estimate the particulate discharge amount with higher accuracy.
[0035]
Furthermore, in such a particulate emission estimation device, a soot generation measuring means for measuring the soot generation amount in the combustion chamber of the diesel engine, and a reference particulate based on the measured soot generation amount in the combustion chamber. A reference particulate discharge amount calculating means for calculating a curate discharge amount, and the particulate discharge amount estimating means is the same as the particulate discharge amount estimated in the operation state in which the reference particulate discharge amount is calculated. It is preferable to be configured to correct the particulate discharge amount estimated in another operating state in accordance with the difference from the calculated reference particulate discharge amount.
[0036]
According to this, the amount of soot generated is actually measured. The particulates contain not only soot but also the above-mentioned SOF and sulfate. Since the amount of SOF and the amount of sulfate are small relative to the amount of soot and correlate with the amount of soot, for example, In a specific operating state such as when the engine is in a steady operating state, the particulate discharge reference amount, that is, the reference particulate discharge amount is calculated with relatively high accuracy based on the measured amount of soot generated. obtain.
[0037]
Then, the particulate discharge amount estimation means estimates the particulate matter in another operation state according to the difference between the estimated particulate discharge amount and the calculated reference particulate discharge amount in the same operation state. Emissions are corrected. In other words, the error of the model expressed by the above formula used for estimating the particulate discharge amount is corrected with a value based on the actual measurement value. Therefore, it is possible to accurately estimate the particulate discharge amount.
[0038]
DETAILED DESCRIPTION OF THE INVENTION
Embodiments of a diesel engine particulate emission estimation device according to the present invention will be described below with reference to the drawings.
[0039]
(1. First embodiment)
FIG. 1 shows a schematic configuration of a system in which an exhaust purification device for an internal combustion engine including a particulate emission estimation device according to a first embodiment of the present invention is applied to a four-cylinder internal combustion engine (diesel engine) 10. This system includes an engine main body 20 including a fuel supply system, an intake system 30 for introducing the atmosphere (fresh air) into the combustion chamber of each cylinder of the engine main body 20, and an exhaust for releasing exhaust gas from the engine main body 20. A system 40, an EGR device 50 for performing exhaust gas recirculation, and an electric control device 60 are included.
[0040]
A fuel injection valve 21 is disposed above each cylinder of the engine body 20. Each fuel injection valve 21 is connected to a fuel injection pump 22 connected to a fuel tank (not shown) via a fuel pipe 23. Thus, the fuel that has been boosted to the injection pressure is supplied from the fuel injection pump 22 to the fuel injection valve 21. Further, the fuel injection valve 21 is electrically connected to the electric control device 60, and is opened for a predetermined time according to a drive signal (command signal) from the electric control device 60, whereby the combustion of each cylinder is performed. The pressurized fuel is injected into the room.
[0041]
The intake system 30 includes an intake manifold 31 connected to a combustion chamber of each cylinder of the engine body 20, an intake pipe 32 connected to an upstream side assembly of the intake manifold 31 and constituting an intake passage together with the intake manifold 31, an intake pipe The throttle valve 33 rotatably held in the throttle 32, the throttle valve actuator 33a that rotationally drives the throttle valve 33 in response to a drive signal from the electric control device 60, and the intake pipe 32 upstream of the throttle valve 33 The intercooler 34, the compressor 35a of the supercharger 35, and the air cleaner 36 disposed at the tip of the intake pipe 32 are included.
[0042]
The exhaust system 40 includes an exhaust manifold 41 connected to each cylinder of the engine body 20, an exhaust pipe 42 connected to a downstream gathering portion of the exhaust manifold 41, and a turbine of the supercharger 35 disposed in the exhaust pipe 42. 35 b, a diesel particulate filter (hereinafter referred to as “DPNR”) 43 that is interposed in the exhaust pipe 42 and carries a catalyst that oxidizes the NOx catalyst and the collected particulates, and the exhaust pipe 42 upstream of the DPNR 43. The additional fuel supply nozzle 44 that supplies the added fuel (light oil) to the DPNR 43, the added fuel supply pipe 45 that connects the added fuel supply nozzle 44 and the fuel injection pump 22, and the added fuel supply pipe 45. An electromagnetic on-off valve 44a for opening and closing the fuel passage of the added fuel supply pipe 45 in response to a drive signal from the electric control device 60, DP R43 downstream exhaust pipe rotatably held in the exhaust throttle valve 46, and a throttle valve 46 exhaust in response to a drive signal from the electronic control device 60 includes an exhaust throttle valve actuator 46a for rotating. The exhaust manifold 41 and the exhaust pipe 42 constitute an exhaust passage.
[0043]
The DPNR 43 includes a filter 43a formed of a material having a large number of pores (ie, a porous material) such as cordierite (a kind of ceramic), and the particulates in the exhaust gas passing through the filter are formed on the surface of the pores. It is a wall-flow type filter having a honeycomb structure to be collected.
[0044]
Further, the surface of the filter of DPNR43 has alumina (Al2OThree) And a catalyst made of a noble metal such as Pt, which oxidizes (and / or has an oxidizing action) the particulates collected on the alumina, is supported.
[0045]
By the way, when the amount of particulates collected in the DPNR 43 becomes larger than the amount of particulates released to the outside of the DPNR 43 by being oxidized in the DPNR 43, the particulates are deposited in the DPNR 43. As a result, the particulate collection ability of the DPNR 43 is reduced, and the exhaust resistance of the engine 10 is increased.
[0046]
For this reason, the present exhaust purification device supplies the added fuel from the aforementioned added fuel supply nozzle 44 to the DPNR 43, and raises the catalyst bed temperature by the heat generated when the added fuel burns inside the DPNR 43, By increasing the oxidation rate of the particulates, the deposited DPNR 43 is more oxidized and removed. In other words, the exhaust purification device periodically regenerates the filter by supplying the added fuel, thereby maintaining the particulate collection ability of the DPNR 43 and preventing the exhaust pressure from rising. Since oxygen is required for the particulates to be oxidized, the supply amount of the added fuel from the added fuel supply nozzle 44 is such that the entire air-fuel ratio of the DPNR 43 does not become rich (that is, the particulates in the DPNR 43). The amount of the air / fuel ratio is sufficient to obtain a lean air-fuel ratio in which oxygen for oxidizing the curate is sufficiently present.
[0047]
Further, DPNR 43 is composed of alumina as the carrier, alkali metal such as potassium K, sodium Na, lithium Li, and cesium Cs, alkaline earth metal such as barium Ba and calcium Ca, and lanthanum La and yttrium Y. At least one selected from rare earth metals is supported together with the supported platinum.
[0048]
Thus, the DPNR 43 absorbs NOx when the air-fuel ratio of the gas flowing into the DPNR 43 is leaner than the stoichiometric air-fuel ratio, and the air-fuel ratio of the gas flowing into the DPNR 43 is richer than the stoichiometric air-fuel ratio or the stoichiometric air-fuel ratio. Sometimes NOx absorbed is NO2Or it releases as NO. Also, this released NO2Alternatively, NO reacts with unburned HC or unburned CO in the rich air-fuel ratio gas flowing into DPNR 43, and N2It is reduced to be released to the outside. Thus, the DPNR 43 functions also as an NOx storage reduction catalyst that absorbs NOx and then releases and reduces the absorbed NOx.
[0049]
The EGR device 50 includes an exhaust gas recirculation pipe 51 constituting a passage (recirculation passage) for recirculating exhaust gas, an EGR valve 52 interposed in the exhaust gas recirculation pipe 51, and an EGR cooler 53. The exhaust gas recirculation pipe 51 communicates the upstream exhaust passage (exhaust manifold 41) of the turbine 35b and the downstream intake passage (intake manifold 31) of the throttle valve 33. The EGR valve 52 responds to the drive signal from the electric control device 60, and the amount of exhaust gas to be recirculated (exhaust gas recirculation amount, EGR amount) is determined by the engine speed NE and the in-cylinder fuel injection amount Qfin per intake stroke. Alternatively, the opening cross-sectional area of the EGR passage constituted by the exhaust gas recirculation pipe 51 is changed so as to be equal to the target EGR amount determined according to the required torque TQ based on the accelerator operation amount Accp. . In addition, you may comprise so that an EGR rate may be controlled by an EGR valve.
[0050]
The electric control device 60 is connected to each other by a bus 61, a ROM 62 that stores programs executed by the CPU 61, tables (lookup tables, maps), constants, and the like, and the CPU 61 temporarily stores data as necessary. The microcomputer includes a RAM 63, a backup RAM 64 that stores data while the power is turned on and holds the stored data while the power is shut off, and an interface 65 including an AD converter.
[0051]
The interface 65 is a new air amount measuring means, and is a hot-wire air flow meter 71 disposed in the intake pipe 32, an intake air temperature sensor 72 provided in the intake pipe 32 downstream of the throttle valve 33, and downstream of the throttle valve 33. Thus, the intake oxygen concentration sensor 73 disposed in the intake passage (intake manifold 31) downstream of the portion where the exhaust gas recirculation pipe 51 is connected, the combustion pressure sensor 74 disposed in the upper combustion chamber of the third cylinder, and the engine body 20 Wall temperature sensor 75 disposed in the cylinder bore liner portion, optical soot sensor 76 disposed in the combustion chamber of the first cylinder, and disposed in the exhaust pipe 42 upstream of the DPNR 43 downstream of the turbine 35b. The exhaust oxygen concentration sensor 77, the engine rotation speed sensor 78, the exhaust gas recirculation delay detection oxygen concentration sensor 79, and the accelerator opening sensor 81 are in contact with each other. It has been, and supplies the signals from these sensors to the CPU 61.
[0052]
The interface 65 is connected to the fuel injection valve 21, the throttle valve actuator 33 a, the electromagnetic on-off valve 44 a that opens and closes the fuel passage of the added fuel supply pipe 45, the exhaust throttle valve actuator 46 a, and the EGR valve 52. Drive signals are sent to these in response to instructions.
[0053]
The hot-wire air flow meter 71 measures the mass flow rate of intake air (intake air amount per unit time, fresh air amount per unit time) passing through the intake passage, and generates a signal Ga representing the mass flow rate. It has become. The intake air temperature sensor 72 detects the temperature of fresh air (ie, intake air temperature) drawn into the cylinder (ie, combustion chamber) of the engine 10 and generates a signal Tb representing the intake air temperature. The intake oxygen concentration sensor 73 detects the oxygen concentration (ie, intake oxygen concentration) in the gas sucked into the cylinder (ie, combustion chamber) of the engine 10 and generates a signal XO2IN representing the intake oxygen concentration. It has become.
[0054]
The combustion pressure sensor 74 detects the pressure in the combustion chamber of the third cylinder, and generates a signal representing the pressure in the combustion chamber. The wall surface temperature sensor 75 detects the combustion chamber wall surface temperature and generates a signal representing a signal Theki representing the combustion chamber wall surface temperature.
[0055]
The optical soot sensor 76 is a cylinder wall surface of the first cylinder, and when the piston of the first cylinder is at the top dead center, the first piston ring disposed at the uppermost position of the piston of the first cylinder. When the piston of the first cylinder is slightly moved toward the bottom dead center, a translucent window is provided in a portion located above the first piston ring, and the reflectance of the window is increased. By measuring, the soot amount Ssoot generated in the combustion chamber of the first cylinder is measured.
[0056]
The exhaust oxygen concentration sensor 77 detects the oxygen concentration (that is, the exhaust oxygen concentration) in the exhaust gas discharged from the cylinder (that is, the combustion chamber) of the engine 10 and generates a signal XO2R representing the exhaust oxygen concentration. It has become. The engine rotation speed sensor 78 detects the rotation speed of the engine 10, generates a signal representing the engine rotation speed NE, and can detect the absolute crank angle of each cylinder. The exhaust gas recirculation delay detection oxygen concentration sensor 79 detects the oxygen concentration in the exhaust gas recirculation pipe 51 immediately before the fresh air newly sucked into the engine 10 and the recirculated exhaust gas are mixed. Concentration detection means, which is provided in the vicinity of the EGR valve 52 of the exhaust gas recirculation pipe 51, and the oxygen concentration in the exhaust pipe flow pipe 51 just before being introduced into the EGR valve 52 (oxygen concentration in the recirculated exhaust gas) And a signal representing the oxygen concentration XO2INR in the recirculated exhaust gas is generated. The accelerator opening sensor 81 detects the operation amount of the accelerator pedal AP, and generates a signal Accp representing the accelerator operation amount.
[0057]
<Particulate Emission Estimation Method of First Embodiment>
Next, a method for estimating the particulate discharge amount by the particulate discharge amount estimation device included in the exhaust gas purification apparatus configured as described above will be described. This particulate matter emission amount estimation apparatus estimates the particulate matter emission amount PM discharged from the engine main body 20 per unit time based on the following equation 1 (the particulate emission amount estimation model described by the equation 1). In Equation 1, Asft is a particulate generation constant (hereinafter also simply referred to as “generation constant”), Asot is a particulate oxidation constant (hereinafter also simply referred to as “oxidation constant”), and Q is a unit time. Fuel injection amount, X02 is the exhaust oxygen concentration (oxygen O2) And C are predetermined constants.
[0058]
[Expression 1]
PM = Asft.Q / (Asot.XO2) + C / exp (Asot.XO2)
[0059]
Here, the basis of Equation 1 will be briefly described. Equations 2 and 3 below are Guang'an models (Soot) which are empirical formulas introduced in SAE paper 930612 (title “Approach to Low Nox and Smoke Emission Engines by Using Phenomenological Simulation”) as models for soot formation and oxidation. Model). In equations (2) and (3), Msf is the amount of soot generated, Asf is a constant, Mfv is the amount of fuel that has been vaporized at the vaporization rate (combustion rate) RHR of the fuel injection amount Mf per stroke, and P is Combustion pressure, R is a gas constant, T is a combustion temperature, Msc is the amount of oxidation of soot, Aso is a constant, and Ms is the amount of soot discharged.
[0060]
[Expression 2]
dMsf / dt = Asf ・ Mfv ・ P0.5・ Exp (-12500 / (R ・ T))
[0061]
[Equation 3]
dMsc / dt = Aso, Ms, XO2, P1.8・ Exp (−14000 / (R ・ T))
[0062]
The change amount dMs / dt of the soot production amount is a value obtained by subtracting the soot oxidation amount (dMsc / dt) from the soot production amount (dMsf / dt).
[0063]
[Expression 4]
dMs / dt = dMsf / dt-dMsc / dt
[0064]
Here, it is assumed that the combustion temperature T, the combustion pressure P, and the combustion ratio RHR are constant even when the engine rotational speed NE and the fuel injection amount Q per unit time are changed, and Mfv = Arhr.・ Mf, P0.5= Apf, exp (-12500 / (RT)) = Atf, P1.8= Apo and exp (−14000 / (R · T)) = Ato, the following formula 5 is obtained from the above formula 2 to the above formula 4. The values Arhr, Apf, Atf, Apo, and Ato are constants.
[0065]
[Equation 5]
dMs / dt = Asf, Arhr, Mf, Apf, Atf-Aso, Ms, XO2, Apo, Ato
[0066]
Therefore, when the particulate generation constant Asft and the particulate oxidation constant Asot are defined as shown in the following equations 6 and 7, the above equation 5 is rewritten as the following equation 8.
[0067]
[Formula 6]
Asft = Asf ・ Arhr ・ Apf ・ Atf
[0068]
[Expression 7]
Asot = Aso ・ Apo ・ Ato
[0069]
[Equation 8]
dMs / dt = Asft ・ Mf−Asot ・ Ms ・ XO2
[0070]
Solving the differential equation for Ms expressed by Equation 8, the following Equation 9 is obtained.
[0071]
[Equation 9]
Ms = Asft · Mf / (Asot · XO2) + C1 / exp (Asot · XO2 · t)
[0072]
Assuming that the operating state of the engine 10 is a steady operating state in which the fuel injection amount Qfin and the engine rotational speed NE do not change, it can be considered that the second term on the right side of Equation 9 does not depend on time. In addition, the Guangan model represented by the above formulas 2 and 3 is a model for soot, and since particulates include SOF and sulfate, the emission amount of particulates other than this soot is generated constant Asft. And reflected in the value of the oxidation constant Asot. Then, by converting these ideas and the unit of the formula 9, the following formula 10 which is the same as the formula 1 is obtained. This formula 10 is an equation describing a model used by the particulate emission estimation device.
[0073]
[Expression 10]
PM = Asft.Q / (Asot.XO2) + C / exp (Asot.XO2)
[0074]
Next, a method of obtaining the particulate discharge amount PM using the above equation (10) will be described. In the above formula 10, Q is the fuel injection amount per unit time (unit: g / h), and fuel injection amount Qfin per stroke (unit: mm)Three/ Stroke) and the engine speed NE. For example, in the case of a four-cylinder diesel engine, there are two explosion strokes per revolution of the crankshaft (engine revolution), so Q = K · Qfin · 2 · NE (K is a constant).
[0075]
Therefore, the fuel injection amount Qfin and the engine rotational speed NE are maintained in a constant state (Qfin = Qfin0, NE = NE0), for example, the exhaust oxygen concentration XO2 is changed by changing the opening of the EGR valve 52, The exhaust oxygen concentration XO2 and the particulate emission amount PM at that time are actually measured. Then, as shown in FIG. 2, at least three points P1 to P3 are measured data points of the particulate emission amount PM for different exhaust oxygen concentrations XO2, and the curve represented by the above equation 10 is the same three points P1 to P3. The generation constant Asft and the oxidation constant Asot are determined by, for example, the least square method.
[0076]
Thereafter, similarly, while maintaining the fuel injection amount Qfin and the engine speed NE at different values, the exhaust gas concentration PM02 is changed to measure the particulate emission amount PM, and the fuel injection amount Qfin is determined based on the data. The production constant Asft and the oxidation constant Asot at the engine speed NE are determined, and as shown in steps 305 and 310 in the flowchart of FIG. 3, respectively, the fuel injection amount Qfin (substantially the fuel injection amount Qfin). And a generation constant table that defines the relationship between the particulate generation constant Asft and the engine speed NE, and the fuel injection amount Qfin (substantially the value corresponding to the fuel injection amount Qfin). ) And an oxidation constant table that defines the relationship between the engine rotational speed NE and the particulate oxidation constant Asot is created and stored in the ROM 62.
[0077]
In actual engine operation, the generation constant Asft and the oxidation constant Asot are determined from the fuel injection amount Qfin and the engine speed NE at that time and the two tables stored in the ROM 62, and the determined generation constant Asft is determined. And the oxidation constant Asot are applied to Equation 10 (Equation 1), the actual fuel injection amount Q per unit time is substituted into Equation 1, and the actual exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 is obtained. By substituting the exhaust oxygen concentration XO2 into the same number 10, the particulate emission amount PM is estimated. The above is the particulate discharge amount estimation method of the present embodiment.
[0078]
According to the above formula 10, the particulate matter discharge amount PM monotonously decreases with respect to the exhaust oxygen concentration XO2. However, this tendency is effective only when the air-fuel ratio of the gas sucked into the engine 10 is an air-fuel ratio larger than about 20, for example, the air-fuel ratio A / F is set to 20 or less, and the combustion temperature is lowered. Therefore, it does not correspond to combustion (so-called low temperature combustion) for suppressing particulate emission. Therefore, the particulate emission estimation method is applied only when the air-fuel ratio of the gas sucked into the engine 10 is larger than the predetermined maximum value (about 20) of the predetermined air-fuel ratio that causes low-temperature combustion.
[0079]
<Operation of First Embodiment>
Next, the operation of the exhaust emission control device including the particulate emission amount estimating device will be described. The CPU 61 of the electric control device 60 is detected by the accelerator operation amount Accp detected by the accelerator opening sensor 81 and the engine rotational speed sensor 78. The fuel injection amount Qfin is determined based on the engine rotational speed NE and the table stored in the ROM 62, and the fuel of the determined fuel injection amount Qfin is injected into the cylinder that has reached the explosion stroke. In addition, a drive signal is generated in the fuel injection valve 21 for the cylinder.
[0080]
Similarly, the CPU 61 determines a target EGR amount (or a target EGR rate) based on the detected accelerator operation amount Accp, the detected engine rotational speed NE, and a table stored in the ROM 62, and the exhaust gas is exhausted. A drive signal is generated in the EGR valve 52 so that the exhaust amount recirculated through the reflux pipe 51 becomes the determined target EGR amount.
[0081]
Further, the CPU 61 executes the particulate discharge amount PM estimation routine shown by the flowchart in FIG. 3 every elapse of a predetermined time TSAMPLE. Accordingly, when the predetermined timing is reached, the CPU 61 starts processing from step 300 and proceeds to step 305 to determine the generation constant Asft based on the engine rotational speed NE and the fuel injection amount Qfin at that time and the above-described generation constant table. In the next step 310, the oxidation constant Asot is determined based on the engine rotational speed NE and the fuel injection amount Qfin and the above-described oxidation constant table.
[0082]
Next, the CPU 61 proceeds to step 315 to acquire the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 as the exhaust oxygen concentration XO2, and proceeds to step 320 to generate the production constant Asft, The particulate emission amount PM is estimated based on the oxidation constant Asot determined in step 310, the exhaust oxygen concentration XO2 acquired in step 315, and the above equation (10). As described in step 320, the fuel injection amount Q per unit time in the above equation 1 is obtained by the following equation 11.
[0083]
[Expression 11]
Q = K ・ Qfin ・ 2 ・ NE
[0084]
Next, the CPU 61 proceeds to step 325 and multiplies the particulate discharge amount PM estimated in step 320 by the value obtained by dividing the execution time interval TSAMPLE (seconds) of this routine by 3600 (seconds), and the result is obtained as of the present time. In addition to the total particulate emission amount SUMPM, a new total particulate emission amount SUMPM is obtained. Thereafter, the CPU 61 once ends this routine at step 395. Thereafter, the CPU 61 repeatedly executes this routine every elapse of a predetermined time TSAMPLE. As a result, the particulate discharge amount PM and the particulate discharge amount total value SUMPM are updated.
[0085]
Further, the CPU 61 executes a filter regeneration control routine for separately regenerating the DPNR 43 (not shown), and when the total particulate discharge amount value SUMPM exceeds a predetermined value, the additional fuel is supplied from the additional fuel supply nozzle 44 to the DPNR 43. The catalyst bed temperature is raised by the heat generated as the added fuel burns inside the DPNR 43, and the deposited DPNR 43 is more oxidized by increasing the oxidation rate of the particulates.
[0086]
As described above, according to the particulate estimator according to the first embodiment, the particulate emission amount is based on the function (the above-mentioned equation 1) using only the exhaust oxygen concentration XO2 and the fuel injection amount Q as variables. PM is estimated. That is, the model represented by Equation 1 is a simple model with few variables, and its calculation is also simple. Therefore, the particulate emission PM can be estimated in real time even by a normal control computer. In addition, since the exhaust oxygen concentration, which is a variable of the model, is a variable that is relatively easy to ensure measurement accuracy compared to the combustion temperature, combustion pressure, and the like, the particulate emission amount PM can be estimated with high accuracy.
[0087]
Further, the particulate generation constant Asft and the particulate oxidation constant Asot in the above equation 1 are determined by the actually measured particulate emission PM when the fuel injection amount Qfin and the engine speed NE are kept constant. Is done. Therefore, the generation constant Asft and the oxidation constant Asot are determined by the combustion temperature, the combustion pressure, the combustion rate, the combustion type, the shape of the fuel injection valve 21 and the engine in each operating state determined by the fuel injection amount Qfin and the engine speed NE. The value reflects all the parameters that affect the particulate emission amount, such as 10 exhaust amount. As a result, it is possible to accurately estimate the particulate discharge amount while reducing the man-hours required for adaptation.
[0088]
Further, since the particulate emission amount PM is estimated as a function of the exhaust oxygen concentration XO2, for example, a delay in the flow of fresh air into the combustion chamber in a transient operation (acceleration operation) state that appears remarkably in an engine with a supercharger, or a transient operation The influence on the particulate emission due to the delay in changing the EGR amount due to the transition to the state can also be grasped by the exhaust oxygen concentration, so that the particulate emission in the transient operation state can be accurately estimated. It becomes possible.
[0089]
Further, even when the setting of the target EGR amount with respect to the fuel injection amount, the engine rotational speed, the cooling water temperature, and the like is changed in the adaptation process of the diesel engine 10, the change in the EGR amount appears in the exhaust oxygen concentration. There is no need to re-adjust the particulate generation constant Asft and the particulate oxidation constant Asot. Therefore, according to the first embodiment, the matching man-hours can be reduced even in such a case.
[0090]
(2. Second Embodiment)
Next, an exhaust emission control device for an internal combustion engine that includes a particulate emission estimation device according to a second embodiment of the present invention will be described. This exhaust purification device improves the accuracy of the value of the current exhaust oxygen concentration (exhaust oxygen concentration used in Equation 1) XO2 in consideration of the delay in exhaust gas recirculation by the EGR device 50, thereby reducing the particulate matter emission amount PM. It differs from the exhaust emission control device of the first embodiment only in improving the estimation accuracy. Therefore, the difference will be mainly described below.
[0091]
<Particulate Emission Estimation Method of Second Embodiment>
The exhaust oxygen concentration XO2 changes according to the operating state of the engine 10. Therefore, when the engine operating state suddenly changes from one operating state to another operating state, the exhaust oxygen concentration XO2 also changes suddenly. On the other hand, the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 shows a behavior that is delayed with respect to the exhaust oxygen concentration XO2 immediately after being discharged from the combustion chamber.
[0092]
On the other hand, a predetermined time is required until the exhaust gas recirculated through the exhaust gas recirculation pipe 51 flows into the combustion chamber again. In other words, a predetermined time (hereinafter referred to as “exhaust gas recirculation delay time”) until a sudden change in the exhaust gas oxygen concentration XO2 due to a sudden change in the operating state of the engine 10 causes a new change in the exhaust gas oxygen concentration XO2 due to exhaust gas recirculation. ").
[0093]
Therefore, in the present embodiment, the exhaust gas recirculation delay time Td is measured and learned, and when the engine operating state changes and the exhaust oxygen concentration XO2 changes suddenly (for example, this is detected by the exhaust oxygen concentration sensor 77). This can be detected by a sudden change in the exhaust oxygen concentration XO2R.) After the elapse of the exhaust gas recirculation delay time Td from the change in the operating state, the exhaust gas whose oxygen concentration has changed with the change in the operating state is transferred to the combustion chamber. By correcting the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 based on the idea that the exhaust oxygen concentration changes abruptly at that time, the accuracy of the exhaust oxygen concentration XO2 used in the above equation 1 is corrected. To increase.
[0094]
Here, a method of measuring the exhaust gas recirculation delay time Td will be described with reference to FIG. 4A, the solid line indicates the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77, and the broken line indicates the recirculated exhaust gas oxygen concentration XO2INR detected by the exhaust gas recirculation delay detection oxygen concentration sensor 79.
[0095]
When the operating state of the engine 10 changes suddenly at time t1, the exhaust oxygen concentration XO2 immediately after being discharged from the combustion chamber immediately starts to change, and the exhaust oxygen concentration XO2 detected by the exhaust oxygen concentration sensor 77 is later than time t1. The change starts at time t2. Further, the oxygen concentration XO2INR in the recirculated exhaust gas detected by the exhaust gas recirculation delay detecting oxygen concentration sensor 79 starts to change at a time t3 later than the time t2. Therefore, the exhaust gas recirculation delay time Td (= t3−t2) can be obtained by knowing the time t2 and the time t3.
[0096]
In the present embodiment, the exhaust gas recirculation delay time Td is obtained in this way, and further, as shown in FIG. 4B, the exhaust oxygen concentration detected by the exhaust oxygen concentration sensor 77 when the operating state of the engine 10 changes suddenly. When the concentration XO2R starts to change at time t4, the change in the exhaust oxygen concentration XO2 detected at time t4 should appear at time t5 after elapse of the exhaust gas recirculation delay time Td from time t4. Therefore, after the time t5, the exhaust oxygen concentration XO2 detected by the exhaust oxygen concentration sensor 77 indicated by the solid line is corrected as indicated by the alternate long and short dash line.
[0097]
In this correction, the initial value is set to the exhaust oxygen concentration correction initial value k1, and thereafter, a correction amount that attenuates toward “0” with time is obtained, and this correction amount is detected by the exhaust oxygen concentration sensor 77. By adding to The exhaust oxygen concentration correction initial value k1 is determined based on the fuel injection amount Qfin and the engine rotational speed NE at the time of change from the steady operation state to the transient operation state, and takes a positive predetermined value or a negative predetermined value. .
[0098]
<Operation of Second Embodiment>
Next, the operation of the exhaust emission control device including the particulate emission estimation device according to the second embodiment will be described. The CPU 61 of this apparatus performs the same processing as the CPU 61 of the first embodiment except that the routines shown in FIGS. 4 to 6 are additionally executed. Therefore, in the following, the description will focus on such differences.
[0099]
First, the operation for obtaining the exhaust gas recirculation delay time Td will be described. The CPU 61 executes the exhaust gas recirculation delay time learning routine shown in FIG. 5 every elapse of a predetermined time. Accordingly, when the predetermined timing is reached, the CPU 61 starts processing from step 500 and proceeds to step 505 to determine whether or not the value of the learning completion flag G is “0”. The value of the learning completion flag G is set to “0” in an initial routine (not shown) when an ignition switch (not shown) is changed from the “off” state to the “on” state.
[0100]
Therefore, if the description is continued assuming that the ignition switch has just been changed from the “off” state to the “on” state, the value of the learning completion flag G is “0”. Therefore, the CPU 61 determines “Yes” in step 505 and proceeds to step 510 to determine whether or not the value of the exhaust gas recirculation delay time measuring flag F is “1”. The value of the flag F is set to “0” in the above-described initial routine.
[0101]
Therefore, the CPU 61 makes a “Yes” determination at step 510 to proceed to step 515, where the absolute value of the difference between the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 and the exhaust oxygen concentration XO2ROLD a predetermined time ago. It is determined whether (ie, the amount of change in the exhaust oxygen concentration XO2R within a predetermined time) is greater than a predetermined threshold TH1. In this case, if the operating state of the engine 10 is in a steady state, the amount of change in the exhaust oxygen concentration XO2R within a predetermined time is small, so the CPU 61 makes a “No” determination at step 515 to proceed to step 595, and this routine Is temporarily terminated.
[0102]
Thereafter, the CPU 61 repeats this routine every elapse of a predetermined time. Accordingly, when the operating state of the engine 10 changes and the amount of change in the exhaust oxygen concentration XO2R within the predetermined time becomes larger than the threshold value TH1, the CPU 61 determines “Yes” when it proceeds to step 515, proceeds to step 520, and exhausts The value of the reflux delay time measuring flag F is set to “1”, the value of the counter C1 is set to “0”, and this routine is temporarily ended.
[0103]
When the CPU 61 executes this routine again after a predetermined time has elapsed, the CPU 61 makes a “Yes” determination at step 505 and changes the value of the exhaust gas recirculation delay time measuring flag F to “1”. Therefore, it is determined as “No” in Step 510 and the process proceeds to Step 530. In Step 530, the oxygen concentration XO2INR in the recirculated exhaust gas detected by the exhaust gas recirculation delay detecting oxygen concentration sensor 79 is the same as the predetermined time before. It is determined whether or not the absolute value of the difference from the recirculated exhaust oxygen concentration XO2INROLD (that is, the amount of change in the recirculated exhaust oxygen concentration XO2INR within a predetermined time) is greater than a predetermined threshold TH2.
[0104]
At the present time, since the amount of change in the exhaust oxygen concentration XO2 within a predetermined time is immediately after the threshold value TH1, the exhaust gas whose exhaust oxygen concentration XO2R has changed is disposed at the position where the oxygen concentration sensor 79 for detecting exhaust gas recirculation delay is provided. Not reached up to. Therefore, since the amount of change in the oxygen concentration XO2INR in the recirculated exhaust gas within the predetermined time is smaller than the threshold value TH2, the CPU 61 makes a “No” determination at step 530 to proceed to step 535 and set the value of the counter C1 to “1”. After the increase, in step 595, this routine is once ended.
[0105]
Thereafter, the CPU 61 repeatedly executes Step 500 to Step 510, Step 530, and Step 535 until the amount of change in the oxygen concentration XO2INR in the recirculated exhaust gas within the predetermined time becomes larger than the threshold value TH2, so the value of the counter C1 is predetermined. It increases by “1” as time passes. When the predetermined amount of time has elapsed and the amount of change in the oxygen concentration XO2INR in the recirculated exhaust gas within the predetermined time becomes greater than the threshold value TH2, when the CPU 61 proceeds to step 530, it determines “Yes” at step 530. In step 540, the value of the counter C1 is stored as the exhaust gas recirculation delay time Td. In the subsequent step 545, the value of the learning completion flag G is set to “1”. Then, in step 595, this routine is temporarily terminated. .
[0106]
In this way, the CPU 61 determines the time from when the change amount of the exhaust oxygen concentration XO2R within the predetermined time becomes larger than the threshold value TH1 until the change amount of the reflux exhaust gas oxygen concentration XO2INR within the predetermined time becomes larger than the threshold value TH2. The time measured by the counter C1 is learned as the exhaust gas recirculation delay time Td.
[0107]
Further, the CPU 61 is configured to repeatedly execute the exhaust oxygen concentration correction initial value k1 determination routine shown in FIG. 6 every elapse of a predetermined time TSAMPLE. Accordingly, when the predetermined timing comes, the CPU 61 starts processing from step 600 and proceeds to step 605, where the absolute value of the difference between the fuel injection amount Qfin and the fuel injection amount Qfinold before the predetermined time becomes larger than the predetermined threshold value TH3. It is determined whether or not. This step is a step provided to determine whether or not the operating state of the engine 10 has become a transient state.
[0108]
Now, assuming that the operation state of the engine 10 is not a transient operation state but a steady operation state, in this case, the absolute value of the difference between the fuel injection amount Qfin and the fuel injection amount Qfinold before a predetermined time is a predetermined threshold value TH3. Therefore, the CPU 61 makes a “No” determination at step 605 to proceed to step 695 to end the present routine tentatively. Thus, when the engine 10 is in a steady operation state, the exhaust oxygen concentration correction initial value k1 is not updated.
[0109]
Thereafter, when the operation state of the engine 10 becomes a transient operation state, the absolute value of the difference between the fuel injection amount Qfin and the fuel injection amount Qfinold before a predetermined time becomes larger than a predetermined threshold value TH3. Accordingly, when the CPU 61 proceeds to step 605, it determines “Yes” in step 605, proceeds to step 610, and sets the value of the transient operation state flag H1 to “1”. The transient operation state flag H1 indicates that the operation state of the engine 10 is in a transient operation state when the value is “1”, and is set to “0” by the above-described initial routine (not shown). Yes.
[0110]
Next, the CPU 61 proceeds to step 615, and the fuel injection amount Qfin and the engine rotational speed NE are maintained based on the fuel injection amount Qfin and the engine rotational speed NE at that time and the table stored in the ROM 62 in advance. In this case, the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 is determined to be greater than the predicted exhaust oxygen concentration SXO2R.
[0111]
When the predicted exhaust oxygen concentration SXO2R is larger than the exhaust oxygen concentration XO2R detected by the current exhaust oxygen concentration sensor 77, the CPU 61 determines “Yes” in step 620 and proceeds to step 625 to correct the exhaust oxygen concentration. The initial value k1 is set to the absolute value of the predetermined value A1, and this routine is once terminated in step 695. On the other hand, when the predicted exhaust gas oxygen concentration SXO2R is smaller than the exhaust gas oxygen concentration XO2R detected by the current exhaust gas oxygen concentration sensor 77, it is determined as “No” in step 620, and the process proceeds to step 630, where the exhaust gas oxygen concentration correction initial stage is set. The value k1 is set to a value in which the sign of the absolute value of the predetermined value A1 is negative, and in step 695 this routine is temporarily ended.
[0112]
In this way, the CPU 61 predicts what the exhaust oxygen concentration after the operating state of the engine 10 has changed will be the predicted exhaust oxygen concentration SXO2R, and the current exhaust oxygen concentration XO2R (that is, the operating state of the engine 10). (Step 620) to determine whether the exhaust oxygen concentration correction initial value k1 should be a positive value or a negative value. Then, the exhaust oxygen concentration correction initial value k1 having a sign corresponding to the magnitude comparison is determined. That is, if it is predicted that the exhaust oxygen concentration XO2R will increase due to a change in the operating state of the engine 10, the exhaust oxygen concentration correction initial value k1 will be a positive value (| A |) and the exhaust oxygen concentration XO2R will decrease. When predicted, the exhaust gas oxygen concentration correction initial value k1 is set to a negative value (-| A |).
[0113]
On the other hand, since the CPU 61 executes the routine shown in FIG. 3 every elapse of a predetermined time, when the predetermined timing is reached, the CPU 61 executes steps 300 to 310 to determine the generation constant Asft and the oxidation constant Asot, and then step 315. In order to obtain the exhaust oxygen concentration XO2, the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 is set as the oxygen concentration XO2, and the processing of the oxygen concentration XO2 calculation routine shown in FIG. Start.
[0114]
Next, the CPU 61 proceeds to step 705 to determine whether or not the value of the correction necessary flag H2 is “1”. The correction necessary flag H2 indicates that the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 when the value is “1” needs to be corrected, and when the value is “0”, the exhaust oxygen concentration. This indicates that the exhaust oxygen concentration XO2R detected by the sensor 77 does not need to be corrected. The value of the correction necessary flag H2 is set to “0” by the above-described initial routine (not shown).
[0115]
Now, assuming that the operation state of the engine 10 is a steady operation state, in such a case, the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 maintains a substantially constant value. It can be used as the exhaust oxygen concentration XO2 for estimation. Therefore, since the value of the correction necessary flag H2 is maintained at “0” (see step 755 described later), the CPU 61 makes a “No” determination at step 705 to proceed to step 710 to enter the transient operation state flag H1. It is determined whether or not the value of “1” is “1”.
[0116]
As described in step 605 and step 610 of the routine shown in FIG. 6, the value of the transient operation state flag H1 is maintained at “0” when the engine 10 is in the steady operation state. Therefore, the CPU 61 makes a “No” determination at step 710 to proceed to step 795 to end the present routine tentatively. As a result, the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 becomes the oxygen concentration XO2 for estimating the particulate discharge amount.
[0117]
Next, a case where the operation state of the engine 10 changes from the steady operation state to the transient operation state will be described. In this case, since the fuel injection amount Qfin changes abruptly, when executing step 605 in FIG. 6, the CPU 61 determines “Yes” in step 605 and proceeds to step 610 to set the value of the transient operation state flag H1 to “ Set to “1”. However, at this stage, the value of the correction necessary flag H2 is maintained at “0”.
[0118]
As a result, when the CPU 61 reaches step 705 of the routine of FIG. 7 when executing step 315 of the routine of FIG. 3, it determines “No” at step 705 and proceeds to step 710, and “Yes” at step 710. And the routine proceeds to step 715, where it is determined whether or not the absolute value of the difference between the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 and the exhaust oxygen concentration XO2ROLD before a predetermined time is greater than a predetermined threshold TH4. To do. In other words, in step 715, the CPU 61 determines whether or not the exhaust oxygen concentration sensor 77 has detected a sudden change in oxygen concentration accompanying a change in the operating state of the engine 10.
[0119]
Since the present time is immediately after the operating state of the engine 10 is changed to the transient operating state, even if the exhaust oxygen concentration in the exhaust immediately after being discharged from the combustion chamber of the engine 10 is changed, the change is the exhaust oxygen concentration. It is not detected by the sensor 77. Accordingly, the CPU 61 makes a “No” determination at step 715 to end the present routine tentatively.
[0120]
Thereafter, when a predetermined time elapses, the exhaust oxygen concentration sensor 77 detects a change in the exhaust oxygen concentration accompanying a change in the operating state. At this time, when the CPU 61 reaches step 715, the CPU 61 determines “Yes” at step 715 and proceeds to step 720 to set the value of the counter C2 to “0” and to correct at step 725. The value of the necessary flag H2 is set to “1”, and then this routine is once ended in step 795.
[0121]
After that, when the CPU 61 executes the particulate discharge amount estimation routine shown in FIG. 3, it determines “Yes” in the process of step 705 accompanying the process of step 315 and proceeds to step 730. The value is increased by “1”, and in the next step 735, it is determined whether or not the value of the counter C2 is equal to or longer than the exhaust gas recirculation delay time Td. At this time, the CPU 61 makes a “No” determination at step 735 to proceed to step 795 to end the present routine tentatively. As a result, the exhaust oxygen concentration XO2 acquired in step 315 is still maintained at the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77.
[0122]
Thereafter, since the CPU 61 repeatedly executes step 730, the value of the counter C2 gradually increases, and when the exhaust gas recirculation delay time Td elapses, the value of the counter C2 coincides with the exhaust gas recirculation delay time Td. Accordingly, the CPU 61 determines “Yes” at step 735 and proceeds to step 740 to correct the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77 according to the following formula 12. Here, k1 is the exhaust gas oxygen concentration correction initial value k1 determined by the routine shown in FIG.
[0123]
[Expression 12]
XO2 = XO2R + k1 / exp (C2-T)
[0124]
Next, in step 745, the CPU 61 determines whether or not the absolute value of the exhaust oxygen concentration correction value k1 / exp (C2-T), which is the second term of the formula 12, is smaller than a predetermined threshold value TH5. The threshold value TH5 is selected to be smaller than the corrected initial value k1. Since C2 = T at the present time, the exhaust oxygen concentration correction value k1 / exp (C2-T) is equal to the initial value k1 and is larger than the threshold value TH5. Therefore, the CPU 61 makes a “No” determination at step 745 to execute the step. Proceed to 795 to end the present routine tentatively. As a result, the oxygen concentration XO2 used in step 320 of the routine shown in FIG. 3 by the CPU 61 becomes a value obtained by adding the initial value k1 to the exhaust oxygen concentration XO2R detected by the exhaust oxygen concentration sensor 77.
[0125]
Thereafter, since such processing is continued and the value of the counter C2 continues to increase in step 730, the absolute value of the exhaust oxygen concentration correction value k1 / exp (C2-T) gradually decreases and in step 740. The required oxygen concentration XO2 is used for estimating the particulate matter discharge amount PM. When the absolute value of the exhaust oxygen concentration correction value k1 / exp (C2-T) becomes smaller than the predetermined threshold value TH5, the CPU 61 makes a “Yes” determination at step 745 to proceed to step 750, where the transient operation state flag H1. In step 755, the value of the correction necessity flag H2 is set to “0”, and in step 795, this routine is temporarily ended.
[0126]
As a result, when the CPU 61 executes the routine shown in FIG. 7 again, the CPU 61 determines “No” in step 705 and “No” in step 710, so that the exhaust oxygen concentration sensor 77 detects. The exhaust oxygen concentration XO2R is used as the exhaust oxygen concentration XO2 for estimating the particulate emission amount PM without being corrected.
[0127]
As described above, the particulate matter emission amount estimation apparatus of the second embodiment has an exhaust oxygen concentration XO2 used for estimating the particulate matter discharge amount PM that is closer to the exhaust oxygen concentration immediately after being discharged from the combustion chamber. Thus, the exhaust gas recirculation delay time is taken into consideration. Accordingly, the estimation accuracy of the particulate discharge amount PM is further improved.
[0128]
(3. Third embodiment)
Next, an exhaust emission control device for an internal combustion engine including a particulate emission estimation device according to a third embodiment of the present invention will be described. In this exhaust purification device, by estimating the particulate discharge amount PM in consideration of the intake oxygen concentration, the estimation accuracy of the particulate discharge amount PM is improved. Specifically, the CPU 61 is shown in FIG. 8 is different from the exhaust purification apparatus of the first embodiment only in that the particulate discharge amount estimation routine shown in FIG. 8 is executed every elapse of a predetermined time TSAMPLE. Therefore, the difference will be mainly described below. In FIG. 8, the same steps as those in FIG. 3 are denoted by the same reference numerals, and detailed description thereof is omitted.
[0129]
At a predetermined timing, the CPU 61 starts the process from step 800 and executes steps 305 to 315 to determine the generation constant Asft and the oxidation constant Asot and obtain the exhaust oxygen concentration XO2. The exhaust oxygen concentration XO2 may be acquired in consideration of the exhaust gas recirculation delay time Td as in the second embodiment.
[0130]
Next, the CPU 61 proceeds to step 805 to read the intake oxygen concentration XO2IN detected by the intake oxygen concentration sensor 73, and in step 810, the intake oxygen concentration correction coefficient table stored in the ROM 62 according to the intake oxygen concentration XO2IN. Based on this, the intake oxygen concentration correction coefficient fXO2 is calculated. The intake oxygen concentration correction coefficient fXO2 in the intake oxygen concentration correction coefficient table is an experimentally adapted value. Then, the CPU 61 estimates the particulate discharge amount PM based on the following equation (13).
[0131]
[Formula 13]
PM = Asft ・ (K ・ Qfin ・ 2 ・ NE) ・ fXO2 / (Asot ・ XO2) + C / exp (Asot ・ XO2)
[0132]
Next, the CPU 61 obtains the particulate discharge total value SUMPM at step 325 and once ends this routine at step 895. Thereafter, the CPU 61 repeatedly executes this routine every elapse of a predetermined time TSAMPLE. As a result, the particulate discharge amount PM and the particulate discharge amount total value SUMPM are updated.
[0133]
Thus, according to the third embodiment, the particulate matter discharge amount PM is estimated based on the intake oxygen concentration XO2IN. Since the intake oxygen concentration XO2IN has a great influence on the combustion temperature, the particulate emission PM is also affected by the intake oxygen concentration XO2IN. Therefore, if the intake oxygen concentration correction coefficient fXO2 is introduced and the particulate discharge amount PM is estimated as in the third embodiment, the estimation accuracy of the particulate discharge amount PM can be improved. The second intake oxygen concentration correction coefficient f2XO2 may be obtained according to the intake oxygen concentration XO2IN and multiplied by the second term on the right side of Equation 13.
[0134]
(4. Fourth embodiment)
Next, an exhaust emission control device for an internal combustion engine including a particulate emission estimation device according to a fourth embodiment of the present invention will be described. This exhaust emission control device improves the estimation accuracy of the particulate discharge amount PM by estimating the particulate discharge amount PM in consideration of the combustion temperature estimated based on the combustion chamber wall surface temperature. Only the point that the CPU 61 executes the particulate discharge amount estimation routine shown in FIG. 9 every elapse of a predetermined time TSAMPLE instead of the routine shown in FIG. 3 is different from the exhaust purification apparatus of the first embodiment. Yes. Therefore, the difference will be mainly described below. In FIG. 9, the same steps as those in FIG. 3 are denoted by the same reference numerals, and detailed description thereof is omitted.
[0135]
At a predetermined timing, the CPU 61 starts processing from step 900 and executes steps 305 to 315 to determine the generation constant Asft and the oxidation constant Asot and obtain the exhaust oxygen concentration XO2. The exhaust oxygen concentration XO2 may be acquired in consideration of the exhaust gas recirculation delay time Td as in the second embodiment.
[0136]
Next, the CPU 61 proceeds to step 905, reads the combustion chamber wall surface temperature Theki detected by the wall surface temperature sensor 75, and obtains the combustion average temperature Tnen by multiplying the combustion chamber wall surface temperature Theki by the heat transfer efficiency η in the subsequent step 910. In step 915, the CPU 61 obtains a correction coefficient g (Tnen) based on the combustion average temperature Tnen, and obtains a corrected generation constant Asft by multiplying it by the generation constant Asft obtained in step 305. Similarly, Then, a correction coefficient h (Tnen) based on the combustion average temperature Tnen is obtained, and this is multiplied by the oxidation constant Asot obtained in step 310 to obtain a corrected oxidation constant Asot. The particulate discharge amount PM is estimated using a constant, and the particulate discharge amount total value SUMPM is obtained in step 325. Then, in step 995, this routine is temporarily ended. The functions g and h for obtaining the correction coefficients g (Tnen) and h (Tnen) are determined in advance by experiments and stored in the ROM 62 as a table.
[0137]
Thereafter, the CPU 61 repeatedly executes the routine shown in FIG. 9 every elapse of a predetermined time TSAMPLE. As a result, the particulate discharge amount PM and the particulate discharge amount total value SUMPM are updated.
[0138]
As described above, according to the fourth embodiment, the generation constant Asft and the oxidation constant Asot are corrected based on the combustion average temperature Tnen, and the patty based on the corrected generation constant Asft and the oxidation constant Asot. The curate discharge amount PM is estimated. The combustion average temperature Tnen affects the particulate matter emission amount PM. Therefore, if the particulate discharge amount PM is estimated according to the combustion average temperature Tnen as in the fourth embodiment, the estimation accuracy of the particulate discharge amount PM can be improved.
[0139]
In the fourth embodiment, the accuracy of estimating the particulate matter emission amount PM is improved by correcting the generation constant Asft and the oxidation constant Asot based on the combustion average temperature Tnen. However, the step 915 is omitted. The particulate matter emission amount PM obtained in step 320 may be modified according to the combustion average temperature Tnen.
[0140]
(5. Fifth embodiment)
Next, an exhaust emission control device for an internal combustion engine including a particulate emission estimation device according to a fifth embodiment of the present invention will be described. This exhaust purification device improves the estimation accuracy of the particulate discharge amount PM by estimating the particulate discharge amount PM in consideration of the intake air temperature. Specifically, the CPU 61 has shown in FIG. It differs from the exhaust emission control device of the first embodiment only in that the particulate emission amount estimation routine shown in FIG. 10 is executed every elapse of a predetermined time TSAMPLE instead of the routine. Therefore, the difference will be mainly described below. 10, the same steps as those in FIG. 3 are denoted by the same reference numerals, and detailed description thereof is omitted.
[0141]
At a predetermined timing, the CPU 61 starts processing from step 1000 and executes steps 305 to 315 to determine the generation constant Asft and the oxidation constant Asot and obtain the exhaust oxygen concentration XO2. The exhaust oxygen concentration XO2 may be acquired in consideration of the exhaust gas recirculation delay time Td as in the second embodiment.
[0142]
Next, the CPU 61 proceeds to step 1005, reads the intake air temperature Tb detected by the intake air temperature sensor 72, and obtains the intake air temperature correction coefficient fTb based on the function i (Tb) according to the intake air temperature Tb read in the subsequent step 1010. . In step 1015, the CPU 61 estimates the particulate discharge amount PM according to the following equation (13). Note that the function i (Tb) is determined in advance by experiment and is stored in the ROM 62 as a table.
[0143]
[Expression 14]
PM = fTb ・ (Asft ・ (K ・ Qfin ・ 2 ・ NE) / (Asot ・ XO2) + C / exp (Asot ・ XO2))
[0144]
Thereafter, the CPU 61 obtains the particulate discharge total value SUMPM in step 325, and then temporarily ends this routine in step 1095. Thereafter, the CPU 61 repeatedly executes this routine every elapse of a predetermined time TSAMPLE. As a result, the particulate discharge amount PM and the particulate discharge amount total value SUMPM are updated.
[0145]
As described above, according to the fifth embodiment, the particulate discharge amount PM is estimated based on the intake air temperature Tb. Since the intake air temperature Tb greatly affects the combustion temperature, the particulate matter discharge amount PM is also affected by the intake air temperature Tb. Therefore, if the particulate discharge amount PM is estimated by introducing the intake air temperature correction coefficient fTb as in the fifth embodiment, the estimation accuracy of the particulate discharge amount PM can be improved.
[0146]
(6. Sixth embodiment)
Next, an exhaust emission control device for an internal combustion engine including a particulate emission estimation device according to a sixth embodiment of the present invention will be described. This exhaust purification device estimates the particulate discharge amount PM in consideration of the compression end temperature (combustion chamber gas temperature in the vicinity of compression top dead center) and the combustion pressure (maximum combustion pressure), thereby releasing the particulate discharge. The point where the estimation accuracy of the quantity PM is improved, specifically, the point where the CPU 61 executes the compression end temperature estimation and combustion pressure detection routine shown in FIG. 11, and the routine shown in FIG. 12 instead of the routine shown in FIG. It differs from the exhaust emission control device of the first embodiment only in that the particulate emission estimation routine shown is executed every elapse of a predetermined time TSAMPLE. Therefore, the difference will be mainly described below. In FIG. 12, the same steps as those in FIG. 3 are denoted by the same reference numerals, and detailed description thereof is omitted.
[0147]
First, the operation for detecting the compression end temperature Ttdc and the combustion pressure Pn will be described. The CPU 61 executes the routine shown in FIG. 11 every elapse of a predetermined crank angle (for example, a crank angle of 1 °). It has become. Therefore, when the crank angle of the engine 10 reaches a predetermined crank angle, the CPU 61 starts processing from step 1100 and proceeds to step 1105 to determine whether or not the current crank angle is the compression top dead center of the third cylinder.
[0148]
If the current crank angle is the compression top dead center of the third cylinder, the CPU 61 proceeds to step 1110 and based on the combustion chamber pressure Pa and the function q detected by the combustion pressure sensor 74 at that time. The estimated compression end temperature Ttdc is obtained (the compression end temperature is estimated), and the process proceeds to Step 1115. On the other hand, if the current crank angle is not the compression top dead center of the third cylinder, the CPU 61 proceeds directly to step 1115.
[0149]
Next, the CPU 61 determines whether or not the combustion chamber pressure Pa detected by the combustion pressure sensor 74 has reached the maximum value after the third cylinder has reached compression top dead center in step 1115. If the combustion chamber pressure Pa takes a maximum value, the combustion chamber pressure Pa, which is the maximum value, is set as the combustion pressure Pn (a value related to the pressure in the combustion chamber during combustion), and the process proceeds to step 1195. The routine is temporarily terminated. On the other hand, if the combustion chamber pressure Pa does not take the maximum value, the routine proceeds to step 1195 as it is, and this routine is once ended. Thus, the compression end temperature Ttdc and the combustion pressure Pn are obtained.
[0150]
At a predetermined timing, the CPU 61 starts processing from step 1200 in FIG. 12 and executes step 305 to step 315 to determine the generation constant Asft and the oxidation constant Asot, and the exhaust oxygen concentration XO2. To get. The exhaust oxygen concentration XO2 may be acquired in consideration of the exhaust gas recirculation delay time Td as in the second embodiment.
[0151]
Next, the CPU 61 reads the compression end temperature Ttdc and the combustion pressure Pn obtained by the routine shown in FIG. 11 in steps 1205 and 1210, respectively. In the subsequent step 1215, the CPU 61 reads the compression end temperature Ttdc and the combustion pressure. A correction coefficient m (Ttdc, Pn) determined by the pressure Pn is obtained from the table m stored in advance in the ROM 62, the compression end temperature Ttdc, and the combustion pressure Pn, and the correction coefficient m (Ttdc, Pn) is obtained in step 305. The generation constant Asft after correction is obtained by multiplying the generation constant Asft obtained in step 5, and similarly, a correction coefficient n (Ttdc, Pn) determined by the compression end temperature Ttdc and the combustion pressure Pn is stored in the ROM 62 in advance. The table n, the compression end temperature Ttdc, and the combustion pressure Pn are obtained, and the correction coefficient n (Ttdc, Pn) is calculated. Request oxidation constants Asot corrected by multiplying the oxidation constants Asot determined by flop 310.
[0152]
Next, the CPU 61 estimates the particulate discharge amount PM in step 320, obtains the particulate discharge amount total value SUMPM in step 325, and then temporarily ends this routine in step 1295.
[0153]
Thereafter, the CPU 61 repeatedly executes the routine shown in FIG. 11 every elapse of a predetermined crank angle, and repeatedly executes the routine shown in FIG. 12 every elapse of a predetermined time TSAMPLE. As a result, the particulate discharge amount PM and the particulate discharge amount total value SUMPM are updated.
[0154]
As described above, according to the sixth embodiment, the compression end temperature Ttdc and the combustion pressure Pn are obtained based on the combustion chamber pressure Pa detected by the combustion pressure sensor 74, and the compression end temperature Ttdc and the combustion pressure are obtained. The production constant Asft and the oxidation constant Asot are corrected based on Pn, and the particulate discharge amount PM is estimated based on the corrected production constant Asft and oxidation constant Asot. Since the compression end temperature Ttdc strongly affects the combustion temperature, the particulate matter discharge amount PM is also affected by the compression end temperature Ttdc. Further, the combustion pressure Pn strongly affects the particulate discharge amount PM. Accordingly, if the particulate discharge amount PM is estimated according to the compression end temperature Ttdc and the combustion pressure Pn as in the sixth embodiment, the estimation accuracy of the particulate discharge amount PM is improved. Can do.
[0155]
In the sixth embodiment, the accuracy of estimating the particulate matter discharge amount PM is improved by correcting the generation constant Asft and the oxidation constant Asot based on the compression end temperature Ttdc and the combustion pressure Pn. 1215 may be omitted, and the particulate discharge amount PM obtained in step 320 may be corrected according to the compression end temperature Ttdc and the combustion pressure Pn.
[0156]
In the sixth embodiment, the pressure P in the combustion chamber is sampled every elapse of a predetermined minute crank angle. However, the pressure P may be sampled every predetermined time. In this case, if the compression end temperature Ttdc and the combustion pressure Pn cannot be obtained with high accuracy, the generation constant Asft and the oxidation constant Asot are corrected in step 1215 according to the time average value Pave of the combustion chamber pressure Pa during the combustion stroke. Alternatively, step 1215 may be omitted and the particulate discharge amount PM obtained in step 320 may be corrected according to the average pressure Pave.
[0157]
(7. Seventh embodiment)
Next, an exhaust emission control device for an internal combustion engine including a particulate emission estimation device according to a seventh embodiment of the present invention will be described. This exhaust purification device is an output of the optical soot sensor 76 when, for example, the operating state of the engine 10 is in a steady operating state and the optical soot sensor 76 can measure the soot amount with relatively high accuracy. Based on the soot amount Ssoot, a reference particulate discharge amount SPM which is an actual particulate discharge amount is obtained, and the difference between the reference particulate discharge amount SPM and the basic particulate discharge amount PM0 estimated by the above equation 1 (Equation 10). Is obtained as the learning value DLPM, and the particulate basic emission amount PM0 estimated by the above equation 1 (Equation 10) is corrected with the same learning value DLPM in other operating states to obtain the particulate emission amount PM. This is different from the exhaust gas purification apparatus of one embodiment. Therefore, the difference will be mainly described below.
[0158]
The CPU 61 of this exhaust purification apparatus repeatedly executes the learning value update routine shown in FIG. 13 every elapse of a predetermined time. Accordingly, when the predetermined timing comes, the CPU 61 starts processing from step 1300 and proceeds to step 1305, where the absolute value of the difference between the fuel injection amount Qfin and the fuel injection amount Qfinold before a predetermined time is smaller than the predetermined threshold value TH6. It is determined whether or not has continued for a predetermined time or more. That is, in step 1305, it is determined whether or not the operating state of the engine 10 is in a steady operating state and the optical soot sensor 76 can measure the soot amount Ssoot with relatively high accuracy.
[0159]
Now, assuming that the state in which the absolute value of the difference between the fuel injection amount Qfin and the fuel injection amount Qfin before the predetermined time is smaller than the predetermined threshold value TH6 has not continued for a predetermined time or longer, the CPU 61 proceeds to a step. In step 1305, it is determined as “No”, the process proceeds to step 1395, and this routine is temporarily ended. Therefore, in this case, the learning value DLPM is not updated.
[0160]
On the other hand, if the state where the absolute value of the difference between the fuel injection amount Qfin and the fuel injection amount Qfinold before the predetermined time is smaller than the predetermined threshold value TH6 continues for a predetermined time or longer, the CPU 61 proceeds to step 1305 In step 1305, “Yes” is determined, and the process proceeds to step 1310 to read the soot amount Ssoot measured by the optical soot sensor 76. Next, the CPU 61 proceeds to step 1315 and obtains the reference particulate discharge amount SPM based on the read soot amount Ssoot and the function r. Note that the function r is obtained in advance by experiments and stored in the ROM 62 as a table.
[0161]
Next, the CPU 61 proceeds to step 1320, where the reference particulate discharge amount SPM and the particulate basic discharge amount PM0 estimated according to the above equation 1 (equation 10) in step 1405 shown in FIG. Is learned as a learning value DLPM, and the routine proceeds to step 1395 to end the present routine tentatively.
[0162]
On the other hand, the CPU 61 repeatedly executes the PM estimation routine shown in FIG. 14 every elapse of a predetermined time TSAMPLE. In FIG. 14, the same steps as those in FIG. 3 are denoted by the same reference numerals, and detailed description thereof is omitted.
[0163]
Therefore, the CPU 61 starts the processing from step 1400 in FIG. 14 at a predetermined timing, and executes the processing from step 305 to step 315 to determine the generation constant Asft and the oxidation constant Asot and to set the exhaust oxygen concentration XO2. get. The exhaust oxygen concentration XO2 may be acquired in consideration of the exhaust gas recirculation delay time Td as in the second embodiment.
[0164]
Next, in step 1405, the CPU 61 estimates the particulate basic discharge amount PM0 according to the above equation 1 (equation 10), and in the subsequent step 1410, adds the learned value DLPM to the particulate basic discharge amount PM0 to thereby increase the particulate discharge rate. The basic emission amount PM0 is corrected. Then, after obtaining the particulate discharge total value SUMPM in step 325, this routine is once ended in step 1495.
[0165]
Thereafter, since the CPU 61 repeatedly executes the routines shown in FIGS. 13 and 14, the learning value DLPM is updated, and the particulate basic discharge amount PM0 is corrected by the learned value DLPM. PM is required with high accuracy.
[0166]
As described above, according to the seventh embodiment, the soot amount Ssoot that is the output of the optical soot sensor 76, that is, the reference particulate discharge amount SPM that is the actual particulate discharge amount based on the actually measured soot amount. And the difference between the reference particulate emission amount SPM and the particulate basic emission amount PM0 estimated by the model of the above equation 1 (Equation 10) is learned as a learning value DLPM. The discharge amount PM0 is corrected.
[0167]
Therefore, since the error of the particulate discharge amount estimation model of Formula 1 (Formula 10) is corrected with a value based on the actual measurement value, the particulate discharge amount can be accurately estimated.
[0168]
As described above, according to each embodiment of the present invention, the particulate emission amount PM can be accurately estimated in real time and in a steady operation state and a transient operation state. In addition, this invention is not limited to the said embodiment, A various modification can be employ | adopted within the scope of the present invention. For example, in the above-described filter regeneration control, not only the supply of added fuel but also the amount of post-injection that is re-injected after main injection from the fuel injection valve 21 is changed, or the opening of the exhaust throttle valve 46 is changed. The supercharging pressure of the supercharger 35 may be changed by changing.
[Brief description of the drawings]
FIG. 1 is a schematic configuration diagram of a system in which an exhaust emission control device provided with a particulate emission estimation device according to a first embodiment of the present invention is applied to a four-cylinder internal combustion engine (diesel engine).
FIG. 2 is a graph for explaining a method of determining a particulate generation constant and a particulate oxidation constant of a particulate emission estimation model.
3 is a flowchart showing a routine that is executed by the CPU shown in FIG. 1 to estimate the particulate discharge amount. FIG.
FIG. 4 (A) shows the exhaust oxygen concentration detected by the exhaust oxygen concentration sensor in order to explain the method for measuring the exhaust gas recirculation delay time by the particulate emission estimation device according to the second embodiment of the present invention. FIG. 4B is a time chart in which the oxygen concentration in the recirculated exhaust gas detected by the exhaust gas recirculation delay detection oxygen concentration sensor is indicated by a broken line with a solid line, and FIG. 5 is a time chart showing the exhaust oxygen concentration corrected using the exhaust gas recirculation delay time indicated by a one-dot chain line.
FIG. 5 is a flowchart showing an exhaust gas recirculation delay time learning routine executed by the CPU of the particulate discharge amount estimation apparatus of the second embodiment.
FIG. 6 is a flowchart showing an exhaust oxygen concentration correction initial value determination routine executed by a CPU of the particulate emission estimation device of the second embodiment.
FIG. 7 is a flowchart showing an exhaust oxygen concentration calculation routine executed by the CPU of the particulate emission estimation device of the second embodiment.
FIG. 8 is a flowchart showing a particulate discharge amount estimation routine executed by a CPU of the particulate discharge amount estimation apparatus of the third embodiment.
FIG. 9 is a flowchart showing a particulate discharge amount estimation routine executed by a CPU of the particulate discharge amount estimation apparatus of the fourth embodiment.
FIG. 10 is a flowchart showing a particulate discharge amount estimation routine executed by a CPU of the particulate discharge amount estimation apparatus of the fifth embodiment.
FIG. 11 is a flowchart showing a compression end temperature estimation and combustion pressure detection routine executed by the CPU of the particulate discharge amount estimation apparatus of the sixth embodiment.
FIG. 12 is a flowchart showing a particulate discharge amount estimation routine executed by a CPU of the particulate discharge amount estimation apparatus of the sixth embodiment.
FIG. 13 is a flowchart showing a learning value update routine executed by the CPU of the particulate discharge amount estimation apparatus of the seventh embodiment.
FIG. 14 is a flowchart showing a particulate discharge amount estimation routine executed by a CPU of the particulate discharge amount estimation apparatus of the seventh embodiment.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 21 ... Fuel injection valve, 31 ... Intake manifold, 32 ... Intake pipe, 33 ... Throttle valve, 35 ... Supercharger, 41 ... Exhaust manifold, 42 ... Exhaust pipe, 43 ... DPNR, 44 ... Addition fuel supply nozzle, 50 ... EGR device 51 ... exhaust pipe flow pipe 72 ... intake air temperature sensor 73 ... intake oxygen concentration sensor 74 ... combustion pressure sensor 75 ... wall surface temperature sensor 76 ... optical soot sensor 77 ... exhaust oxygen concentration sensor 78 ... engine rotation speed sensor, 79 ... oxygen concentration sensor for detecting exhaust gas recirculation delay, Asft ... particulate generation constant, Asot ... particulate oxidation constant.

Claims (11)

ディーゼル機関の燃料噴射量及び機関回転速度を一定に維持した定常運転状態にて排気酸素濃度を変更するとともに、同変更された異なる排気酸素濃度のそれぞれに対してパティキュレート排出量を測定し、同変更した各排気酸素濃度及び同各排気酸素濃度に対して測定された同各パティキュレート排出量からなるデータに基いて、同一定に維持した燃料噴射量及び同一定に維持した機関回転速度で定まる運転状態に対するパティキュレート生成定数及びパティキュレート酸化定数を決定する作業を、それぞれが所定の燃料噴射量及び所定の機関回転速度で定まる複数の運転状態について実行することにより形成された、燃料噴射量及び機関回転速度の組み合わせとパティキュレート生成定数との関係を規定した生成定数テーブルと、燃料噴射量及び機関回転速度の組み合わせとパティキュレート酸化定数との関係を規定した酸化定数テーブルと、を記憶した記憶手段と、
実際の燃料噴射量を表す値を取得する燃料噴射量取得手段と、
実際の機関回転速度を検出する機関回転速度検出手段と、
実際の排気酸素濃度を表す値を取得する排気酸素濃度取得手段と、
前記取得された燃料噴射量を表す値と、前記検出された機関回転速度と、前記生成定数テーブル及び前記酸化定数テーブルと、に基いてパティキュレート排出量を推定するために使用するパティキュレート生成定数とパティキュレート酸化定数とを決定するとともに、同決定されたパティキュレート生成定数と、同決定されたパティキュレート酸化定数と、前記取得された排気酸素濃度を表す値とに基いてパティキュレート排出量を推定するパティキュレート排出量推定手段と、
を備えてなるディーゼル機関のパティキュレート排出量推定装置において、
前記パティキュレート排出量推定手段は、単位時間あたりのパティキュレート排出量をＰＭ、前記決定されたパティキュレート生成定数をＡsft、前記決定されたパティキュレート酸化定数をＡsot、前記取得された燃料噴射量を表す値に基く単位時間あたりの燃料噴射量をＱ、前記取得された排気酸素濃度を表す値をＸ02、及び定数をＣとするとき、ＰＭ＝Ａsft・Ｑ／（Ａsot・ＸO2）＋Ｃ／exp（Ａsot・ＸO2）なる式に基いて前記パティキュレート排出量を推定するように構成されたディーゼル機関のパティキュレート排出量推定装置。The exhaust oxygen concentration was changed in a steady operation state in which the fuel injection amount and engine speed of the diesel engine were kept constant, and the particulate emissions were measured for each of the changed different exhaust oxygen concentrations. Based on the data consisting of the changed exhaust oxygen concentrations and the same particulate emissions measured for each exhaust oxygen concentration, the fuel injection amount maintained at the same constant and the engine speed maintained at the same constant are determined. A fuel injection amount formed by executing the operation of determining the particulate generation constant and the particulate oxidation constant for the operating state for a plurality of operating states, each of which is determined by a predetermined fuel injection amount and a predetermined engine speed, and A generation constant table that defines the relationship between the combination of engine speed and particulate generation constant, and fuel injection Storage means for storing the oxidizing constant table that defines the relationship between the combination of the amount and the engine rotational speed and the particulate oxidation constants, and
Fuel injection amount acquisition means for acquiring a value representing the actual fuel injection amount;
An engine speed detecting means for detecting an actual engine speed;
Exhaust oxygen concentration acquisition means for acquiring a value representing the actual exhaust oxygen concentration;
The particulate generation constant used for estimating the particulate discharge amount based on the value representing the acquired fuel injection amount, the detected engine speed, the generation constant table and the oxidation constant table. And the particulate oxidation constant, and the particulate emission amount based on the determined particulate generation constant, the determined particulate oxidation constant, and the value representing the obtained exhaust oxygen concentration. A particulate emission estimating means for estimating;
In particulate emissions estimation equipment of a diesel engine comprising comprise,
The particulate emission estimation means is PM for particulate emission per unit time, Asft for the determined particulate generation constant, Asot for the determined particulate oxidation constant, and the acquired fuel injection amount . PM = Asft · Q / (Asot · XO 2) + C / exp () where Q is the fuel injection amount per unit time based on the value expressed , X 02 is the value representing the acquired exhaust oxygen concentration, and C is the constant. A particulate emission estimation device for a diesel engine configured to estimate the particulate emission based on the formula Asot · XO2). 請求項１に記載のディーゼル機関のパティキュレート排出量推定装置において、
前記パティキュレート排出量推定手段は、前記取得された排気酸素濃度を表す値が所定濃度以上の値となるときにのみ前記パティキュレート排出量の推定を行うように構成してなるディーゼル機関のパティキュレート排出量推定装置。 In the diesel engine particulate emission estimation apparatus according to claim 1 ,
The particulate emission estimating means is configured to estimate the particulate emission only when the acquired value representing the exhaust oxygen concentration becomes a value equal to or higher than a predetermined concentration. Emission estimation device. 請求項１又は請求項２に記載のディーゼル機関のパティキュレート排出量推定装置であって、
前記排気酸素濃度取得手段は、
前記ディーゼル機関の排気通路内の酸素濃度を検出する排気酸素濃度検出手段を含むとともに、同ディーゼル機関が備える排気還流管を介して排気還流を行う排気還流手段によって同機関の排気が還流されて同機関に吸入されるまでに要する排気還流遅れ時間に応じ、同検出された排気通路内の酸素濃度を補正して前記排気酸素濃度を表す値を取得するように構成されてなるディーゼル機関のパティキュレート排出量推定装置。A diesel engine particulate emission estimation apparatus according to claim 1 or claim 2 ,
The exhaust oxygen concentration acquisition means includes
Exhaust oxygen concentration detection means for detecting the oxygen concentration in the exhaust passage of the diesel engine is included, and the exhaust gas of the engine is recirculated by an exhaust gas recirculation means that performs exhaust gas recirculation through an exhaust gas recirculation pipe provided in the diesel engine. Particulates of a diesel engine configured to obtain a value representing the exhaust oxygen concentration by correcting the detected oxygen concentration in the exhaust passage according to the exhaust gas recirculation delay time required until being sucked into the engine Emission estimation device. 請求項３に記載のディーゼル機関のパティキュレート排出量推定装置であって、
前記排気酸素濃度取得手段は、
前記還流される排気と前記ディーゼル機関に新たに吸入される新気とが混合される直前の前記排気還流管内の酸素濃度を検出する排気還流遅れ検出用酸素濃度検出手段を含むとともに、前記検出される排気通路内の酸素濃度と同検出される排気還流管内の酸素濃度とに基いて前記排気還流遅れ時間を推定するように構成されてなるディーゼル機関のパティキュレート排出量推定装置。A diesel engine particulate emission estimation device according to claim 3 ,
The exhaust oxygen concentration acquisition means includes
An oxygen concentration detection means for detecting an exhaust gas recirculation delay for detecting an oxygen concentration in the exhaust gas recirculation pipe immediately before the exhaust gas to be recirculated and fresh air newly sucked into the diesel engine are mixed. A diesel engine particulate emission estimation device configured to estimate the exhaust gas recirculation delay time based on the oxygen concentration in the exhaust passage and the detected oxygen concentration in the exhaust gas recirculation pipe. 請求項１又は請求項２に記載のディーゼル機関のパティキュレート排出量推定装置であって、
前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼温度に関連する値を取得するとともに、同取得された燃焼温度に関連する値に基いて前記パティキュレート排出量を推定するように構成されてなるディーゼル機関のパティキュレート排出量推定装置。A diesel engine particulate emission estimation apparatus according to claim 1 or claim 2 ,
The particulate emission estimation means is configured to acquire a value related to the combustion temperature of the diesel engine and to estimate the particulate emission based on a value related to the acquired combustion temperature. Particulate emission estimation device for diesel engines. 請求項５に記載のディーゼル機関のパティキュレート排出量推定装置において、
前記パティキュレート排出量推定手段は、前記ディーゼル機関の吸気通路内の酸素濃度を検出する酸素濃度検出手段を含み、同検出される吸気通路内の酸素濃度を前記燃焼温度に関連する値として扱うように構成されてなるディーゼル機関のパティキュレート排出量推定装置。 In the diesel engine particulate emission estimation device according to claim 5 ,
The particulate emission estimation means includes oxygen concentration detection means for detecting an oxygen concentration in the intake passage of the diesel engine, and treats the detected oxygen concentration in the intake passage as a value related to the combustion temperature. A diesel engine particulate emission estimation device. 請求項５に記載のディーゼル機関のパティキュレート排出量推定装置において、
前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼室壁面の温度を検出する燃焼室壁面温度検出手段を含み、同検出される燃焼室壁面温度を前記燃焼温度に関連する値として扱うように構成されてなるディーゼル機関のパティキュレート排出量推定装置。 In the diesel engine particulate emission estimation device according to claim 5 ,
The particulate emission amount estimation means includes combustion chamber wall surface temperature detection means for detecting the temperature of the combustion chamber wall surface of the diesel engine, and treats the detected combustion chamber wall surface temperature as a value related to the combustion temperature. A diesel engine particulate emission estimation device. 請求項５に記載のディーゼル機関のパティキュレート排出量推定装置において、
前記パティキュレート排出量推定手段は、前記ディーゼル機関の吸気温度を検出する吸気温度検出手段を含み、同検出される吸気温度を前記燃焼温度に関連する値として扱うように構成されてなるディーゼル機関のパティキュレート排出量推定装置。 In the diesel engine particulate emission estimation device according to claim 5 ,
The particulate emission estimation means includes intake air temperature detection means for detecting an intake air temperature of the diesel engine, and is configured to handle the detected intake air temperature as a value related to the combustion temperature. Particulate emission estimation device. 請求項５に記載のディーゼル機関のパティキュレート排出量推定装置において、
前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼室内の圧力を検出する燃焼室内圧力検出手段を含み、同検出される燃焼室内圧力を前記燃焼温度に関連する値として扱うように構成されてなるディーゼル機関のパティキュレート排出量推定装置。 In the diesel engine particulate emission estimation device according to claim 5 ,
The particulate emission estimation means includes combustion chamber pressure detection means for detecting the pressure in the combustion chamber of the diesel engine, and is configured to handle the detected pressure in the combustion chamber as a value related to the combustion temperature. The diesel engine particulate emission estimation device. 請求項１又は請求項２に記載のディーゼル機関のパティキュレート排出量推定装置であって、
前記パティキュレート排出量推定手段は、前記ディーゼル機関の燃焼中における燃焼室内圧力に関連する値を取得するとともに、同取得された燃焼室内圧力に関連する値に基いて前記パティキュレート排出量を推定するように構成されてなるディーゼル機関のパティキュレート排出量推定装置。A diesel engine particulate emission estimation apparatus according to claim 1 or claim 2 ,
The particulate emission estimating means acquires a value related to the pressure in the combustion chamber during combustion of the diesel engine and estimates the particulate discharge based on the acquired value related to the pressure in the combustion chamber. A diesel engine particulate emission estimation device constructed as described above. 請求項１又は請求項２に記載のディーゼル機関のパティキュレート排出量推定装置であって、
前記ディーゼル機関の燃焼室内の煤発生量を測定する煤発生量測定手段と、
前記測定された燃焼室内の煤発生量に基いて基準パティキュレート排出量を算定する基準パティキュレート排出量算出手段とを備えるとともに、
前記パティキュレート排出量推定手段は、前記基準パティキュレート排出量が算出された運転状態において推定された前記パティキュレート排出量と同算出された基準パティキュレート排出量との差に応じて、他の運転状態において推定されるパティキュレート排出量を修正するように構成されてなるディーゼル機関のパティキュレート排出量推定装置。A diesel engine particulate emission estimation apparatus according to claim 1 or claim 2 ,
Soot generation amount measuring means for measuring the soot generation amount in the combustion chamber of the diesel engine;
A reference particulate emission amount calculating means for calculating a reference particulate emission amount based on the soot generation amount measured in the combustion chamber,
The particulate discharge amount estimation means is configured to perform another operation according to a difference between the calculated particulate discharge amount and the calculated reference particulate discharge amount in the operation state in which the reference particulate discharge amount is calculated. A particulate emission estimation device for a diesel engine configured to correct a particulate emission estimated in a state.
JP2002111925A 2002-04-15 2002-04-15 Particulate emission estimation device for diesel engine Expired - Fee Related JP4123811B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002111925A JP4123811B2 (en) 2002-04-15 2002-04-15 Particulate emission estimation device for diesel engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002111925A JP4123811B2 (en) 2002-04-15 2002-04-15 Particulate emission estimation device for diesel engine

Publications (2)

Publication Number Publication Date
JP2003307110A JP2003307110A (en) 2003-10-31
JP4123811B2 true JP4123811B2 (en) 2008-07-23

Family

ID=29394588

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002111925A Expired - Fee Related JP4123811B2 (en) 2002-04-15 2002-04-15 Particulate emission estimation device for diesel engine

Country Status (1)

Country Link
JP (1) JP4123811B2 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4603951B2 (en) 2005-08-08 2010-12-22 トヨタ自動車株式会社 Soot generation amount estimation device for internal combustion engine
DE102006000845A1 (en) * 2006-01-05 2007-07-12 Volkswagen Ag Method for determining the soot load of a particulate filter
JP4910844B2 (en) * 2007-04-05 2012-04-04 トヨタ自動車株式会社 Exhaust gas purification device for internal combustion engine
JP5057453B2 (en) * 2007-09-21 2012-10-24 独立行政法人産業技術総合研究所 Method for measuring transient particle mass emission concentration from engine
JP5126554B2 (en) * 2008-05-16 2013-01-23 トヨタ自動車株式会社 Soot emission estimation device for internal combustion engine
ES2771358T3 (en) 2009-06-03 2020-07-06 Toyota Motor Co Ltd Device to estimate the amount of carbon discharged from an internal combustion engine
US8332124B2 (en) * 2009-11-24 2012-12-11 Corning Incorporated Mass based methods and systems for estimating soot load
JP5429408B2 (en) * 2011-01-20 2014-02-26 トヨタ自動車株式会社 Control device for compression ignition type internal combustion engine and method for determining smoke generation state in compression ignition type internal combustion engine
CN115405404A (en) * 2022-09-29 2022-11-29 潍柴动力股份有限公司 Construction method and device of Soot original emission model and engine aftertreatment system

Also Published As

Publication number Publication date
JP2003307110A (en) 2003-10-31

Similar Documents

Publication Publication Date Title
JP4606965B2 (en) Exhaust gas purification device for internal combustion engine
KR100611550B1 (en) Egr control device and method for internal combustion engine
JP4042399B2 (en) Exhaust purification device
US8037675B2 (en) Exhaust gas purification system for internal combustion engine and method for exhaust gas purification
US7779680B2 (en) Estimation of engine-out NOx for real time input to exhaust aftertreatment controller
JP3925485B2 (en) NOx emission estimation method for internal combustion engine
US20080295491A1 (en) Exhaust gas purification device for internal combustion engine
JP4440823B2 (en) Exhaust gas purification device for internal combustion engine
CN110735697A (en) method, system and engine for controlling DPF regeneration
JP4044908B2 (en) Exhaust gas purification device for internal combustion engine
JP5846286B2 (en) Exhaust gas purification device for internal combustion engine
JP4123811B2 (en) Particulate emission estimation device for diesel engine
JP5776619B2 (en) Exhaust purification equipment
US7827783B2 (en) Exhaust gas purifying apparatus for internal combustion engine
JP4759496B2 (en) Exhaust gas purification device for internal combustion engine
JP4061995B2 (en) Exhaust gas purification device for internal combustion engine
JP4012043B2 (en) Particulate filter regeneration method
EP1536120B1 (en) Exhaust gas control apparatus for internal combustion engine and control method thereof
JP2008128217A (en) Exhaust emission control device of internal combustion engine
JP4692376B2 (en) Exhaust gas purification device for internal combustion engine
JP5699957B2 (en) Exhaust gas purification device for internal combustion engine
JP5370252B2 (en) Exhaust gas purification device for internal combustion engine
JP2013160106A (en) Exhaust emission control device for internal combustion engine
JP2011220260A (en) Engine control device
JP2010090875A (en) Exhaust gas control device for internal combustion engine

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20041215

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20071225

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080108

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080307

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20080415

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080428

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110516

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110516

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120516

Year of fee payment: 4

LAPS Cancellation because of no payment of annual fees