(57)【要約】【課題】知られたモデルよりはるかに信頼することができ、現在使用されているエンジン中央制御ユニットにおいて実際に実行することができる微粒子フィルタに蓄積された微粒子の量を決定する方法を提供する。【解決手段】微粒子フィルタに蓄積された微粒子量（ｍS）を決定する方法であって、前記微粒子フィルタにおける微粒子の空間的な分布の変化および／または前記微粒子フィルタにおけるエンジンの動作状態および過去の微粒子蓄積の関数として微粒子の物理化学的特性の変化を決定する方法。前記決定工程は、マップの形態で、前記微粒子フィルタに蓄積された微粒子量（ｍS）と前記微粒子フィルタにおける圧力降下（ΔDPF）との関係を定義するパラメータ（β）の基準値（βMAP）の数を決定する工程であって、各規準値（βMAP）が、前記微粒子フィルタにおける特定の微粒子蓄積条件と対応する個々の一定のエンジン動作条件に関係する工程を含む。【選択図】図４
2008286208.tif 000002
【特許請求の範囲】
【請求項１】
微粒子フィルタ（９）に蓄積された微粒子量（ｍS）を決定する方法であって、
前記微粒子フィルタ（９）における微粒子の空間的な分布の変化および／または前記微粒子フィルタ（９）におけるエンジンの動作状態および過去の微粒子蓄積の関数として、微粒子の物理化学的特性の変化を決定する工程を含むことを特徴とする方法。
【請求項２】
前記微粒子フィルタ（９）における微粒子の空間的な分布の変化および／または微粒子の物理化学的特性の変化が、排気ガス温度および流量と前記微粒子フィルタ（９）における過去の微粒子蓄積との関数として決定される請求項１記載の方法。
【請求項３】
前記決定する工程が、
前記微粒子フィルタ（９）に蓄積された微粒子量（ｍS）と前記微粒子フィルタ（９）における圧力降下（ΔDPF）との関係を定義するパラメータ（β）の基準値（βMAP）の数を決定する工程であって、各規準値（βMAP）が、前記微粒子フィルタ（９）における特定の微粒子蓄積条件に対応する個々の一定のエンジン動作状態に関係する工程と、
同一の動作状態に対する前記パラメータ（β）の基準値（βMAP）および前記微粒子フィルタ（９）における過去の微粒子蓄積の関数として、現在のエンジン動作状態に対する動作値（βMOD）を決定する工程と、
前記現在のエンジン動作状態のために算出されたパラメータ（β）の動作値（βMOD）の関数として現在のエンジン動作状態で微粒子フィルタ（９）に蓄積された微粒子量を決定する工程
とを含む請求項１または２記載の方法。
【請求項４】
前記パラメータ（β）の基準値（βMAP）のそれぞれが、初期の較正工程において、所定長さの時間のあいだ、前記基準値（βMAP）に関連した安定したエンジン動作状態で、微粒子フィルタ（９）の微粒子を集積することによって決定される請求項３記載の方法。
【請求項５】
前記微粒子フィルタ（９）内の微粒子の状態インデックスＩPARTを決定する工程と、
前記微粒子フィルタ（９）内の微粒子の前記状態インデックスＩPARTに基づいて、前記微粒子フィルタ（９）の再生を実行する工程
とを含む請求項１、２、３または４記載の方法。
【請求項６】
たとえばＮＯｘベースの自然発生的な再生のような部分的な再生および停車効果のような延長された車両の廃止による、前記微粒子フィルタ（９）内の微粒子の部分的な分布の変化および／または前記微粒子の物理化学的特性の変化を、決定する工程を備える請求項１、２、３、４または５記載の方法。
【請求項７】
微粒子フィルタ（９）内に蓄積された微粒子の最終的な補正値ｍS_CORRを評価する工程であって、オフセットΔｍの計算により、たとえばＮＯｘベースの自然発生的な再生のような部分的な再生および停車効果のような延長された車両の廃止による、前記微粒子フィルタ（９）内の微粒子の部分的な分布の変化および／または前記微粒子の物理化学的特性の変化を、考慮することなく評価された微粒子フィルタに蓄積された微粒子の値（ｍS）に加える工程を含む請求項１、２、３、４、５または６記載の方法。
【請求項８】
キーオン後のある大きさの時間ののち、オフセットΔｍを以下の式
Δｍ＝ｍS_CORR_ECU−ｍS
（ここで、ｍS_CORR_ECUは、以前のキーオフにおいて、電子制御ユニットに記憶されたｍS_CORRの値であり、
ｍSは、停車効果のような延長された車両の廃止による、前記微粒子フィルタ（９）内の微粒子の部分的な分布の変化および／または前記微粒子の物理化学的特性の変化を考慮することなく評価された微粒子フィルタに蓄積された微粒子の値である）
で決定する請求項１、２、３、４、５、６または７記載の方法。
【請求項９】
たとえばＮＯｘベースの自然発生的な再生のような部分的な再生ののちにオフセットΔｍを以下の式
Δｍ＝ｍS_CORR−ｍS
（ここで、ｍS_CORRは、たとえばＮＯｘベースの自然発生的な再生のような部分的な再生の前に微粒子フィルタ内に蓄積された微粒子の値であり、毎回の工程において、前記部分的な再生状態のあいだ、エンジン動作状態の関数である、実験的に評価された瞬間的な再生された微粒子の値ｍNOxであり、
ｍSは、たとえばＮＯｘベースの自然発生的な再生のような部分的な再生による、前記微粒子フィルタ（９）内の微粒子の部分的な分布の変化および／または前記微粒子の物理化学的特性の変化を考慮することなく評価された微粒子フィルタに蓄積された微粒子の値である）
で決定する工程を備える請求項１、２、３、４、５、６、７または８記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は微粒子フィルタに蓄積された微粒子の量を決定する方法に関する。
【背景技術】
【０００２】
本発明は、自動車への応用、とりわけ内燃エンジン（とくにディーゼルエンジン）だけに使用されるわけではないが、これらに利するために使用することができるため、以下の説明は、これらを例にとってなされている。
【０００３】
本発明は、実際、燃料バーナーなど、微粒子フィルターが設けられたシステムによって発生された微粒子をフィルタにかけるなど、自動車への応用以外において使用することができる。
【０００４】
多くの国で知られているように、大気汚染を規制する規則は、内燃機関の排気ガスの成分についてますます厳しくなってきている。
【０００５】
ディーゼルエンジンの場合、一酸化炭素（ＣＯ）および炭化水素（ＨＣ）によってもたらされる問題は酸化窒素（ＮＯｘ）および排気ガス中の微粒子によってもたらされる問題ほど多くはない。
【０００６】
１９９０年代の終わりに市場に出された最近の直接噴射ディーゼルエンジンでさえ重大な微粒子の発生によって特徴づけられていることを種々の出典が示している。
【０００７】
大気中に発生された排気ガスの微粒子の含有量を最小にする種々の方法が提案されている。もちろん、排気管（ないしは、排気パイプ）にフィルタを設けることは、ディーゼルエンジンの微粒子排出に対する疑いのない最終的な解決策としてエンジンの技術において長いあいだ認められてきた。
【０００８】
微粒子トラップ（煤キャッチャーまたはスートトラップ）としても知られた微粒子フィルタは、通常、孔を有する壁に交互に塞がれた多数の通路を備えている。
【０００９】
より詳しくは、障害物（obstruction）がガスを通路の側壁を介して流し、その結果、微粒子を構成している未燃焼の粒子は、まず側壁の孔に保持される。孔が塞がれると、該未燃焼の粒子は、通路壁の内面上に有孔層を蓄積し形成する。
【００１０】
微粒子が通路壁の内面上に蓄積すると、フィルタを通して圧力が下がり、フィルタによって発生した背圧が増加する。
【００１１】
除去されない場合、微粒子の過度の蓄積は、
性能、駆動の快適さ、およびエンジンの燃費が損なわれ、ついにはエンジンが停止すること、および
微粒子の自発火および制御不能燃焼の場合にフィルタ自体が破壊すること
を引き起こす。とくに駆動条件において、微粒子の多大の蓄積は、突然の、制御不能の微粒子燃焼を引き起こし、ひいてはフィルタのセラミックマトリックスのオーバーヒートおよびフィルタ自体への損傷を引き起こすといった「臨界的な」再生（regeneration）現象を引き起こす。
【００１２】
捕捉された微粒子は、微粒子フィルタを「再生する」ことによって規則的に除去されなければならない。「再生する」とは、エンジンの技術において、蓄積された微粒子（主として炭素Ｃ）を燃焼することを意味し、該蓄積された微粒子は、排気ガス中の酸素と接触し、一酸化炭素および二酸化炭素になる。
【００１３】
しかしながら、この反応は、約６００℃を超える温度において（添加剤を使用することなく）自然に起こるが、この温度は通常のエンジン動作状態におけるフィルタ入口温度より高い。
【００１４】
ある条件（すなわち、フィルタにおける微粒子の一定の蓄積で）、フィルタ入口における排気ガス温度は微粒子の燃焼を開始するために人工的に上昇しなければならない。
【００１５】
多くの微粒子の燃焼を開始するためにフィルタ入口における排気ガス温度を人工的に上昇する方法は提案され、および／または実際に実施されている。
【００１６】
現在の車両の微粒子燃焼開始方法の主な欠点の１つは、マイルベース（たとえば、５００または１０００ｋｍ毎に）周期的に開始される微粒子フィルタの再生にある。
【００１７】
マイルベース微粒子フィルタの再生は、不必要なフィルタの再生を頻繁に開始したり、逆に実際に必要とされるときに、フィルタの再生を開始しないことによって、疑う余地なく、全体として不充分であることが最近分かった。
【００１８】
フィルタにおける微粒子蓄積を示すテストはマイルベースではないが、車両の種別プロファイル（たとえば、都市の走行用または高速道路の走行用）、運転スタイル（たとえば、競争用）など車両の他の要因によってかなり重大な影響を受ける。そのうえ、微粒子排出の増大をかなり引き起こすエンジンの欠点が考慮されていない、その場合、微粒子がフィルタに著しく蓄積し、エンジンおよびフィルタの信頼性を損なう。
【００１９】
一定の量を超えたときに、フィルタの再生を開始し得るように、フィルタに蓄積された微粒子の量を正確に決定するために、フィルタを通過する排気ガス流量および温度と圧力降下にもとづいて蓄積された微粒子の量が算出される、多くの物理的モデルが提案されている。提案されているすべてのモデルは、エンジンの種々動作状態および過去の微粒子の蓄積から独立して、実質的に、フィルタチャネル内の微粒子の分布および微粒子の物理化学的特性が一様かつ一定であるという仮定に基いている。
【００２０】
たとえば、非特許文献１においては、チャネル内部の（軸方向および半径方向）空間的に一様な分布を仮定して、フィルタの幾何学的特徴、フィルタ材料の物理的特徴、微粒子の特徴などの種々の要因を考慮する分析的な微粒子フィルタモデルが提案されており、当該要因はエンジンの種々の動作状態および過去の微粒子蓄積から独立して一定であると仮定されている。
【００２１】
明確性の観点から、つぎの式は、微粒子フィルタのモデルとして非特許文献１においてなされた考慮から本発明者によって開発されたものであり、フィルタの上流に触媒コンバータをもつことができる。
【００２２】
【数１】
2008286208.tif 000003
【００２３】
ここに、Ｖtrapはフィルタの容積、Ｈはセルの大きさ、Ｌは長さ、Ｎは開放セルの数、ｗは壁の厚さをそれぞれ示し、
ｋmは線形透過性を示し、ｋ’mは非線形透過性を示し、
Ｖcatは容積、Ｈcatはセルの大きさ、Ｌcatは長さ、ｗcatは壁の長さをそれぞれ示し、
Ｒはガス定数（８．３１４Ｊ／ｋｍｏｌ）、Ｆは正方形断面の導管内のガス摩擦係数（約２８．４５４）、ξは慣性係数（約３）をそれぞれ示し、
Ｐatmはフィルタの下流における絶対圧力（大雑把には、大気圧に等しいと考えられる）、Ｍwはガスの平均分子量、Ｔは温度、μoは粘性係数をそれぞれ示し、
ｍsは質量、ｋsは透過率、ρsは密度をそれぞれ示し、
ΔDPEはフィルタを横切る全圧力降下、Ｑmは排気ガスの質量流量をそれぞれ示す。
【００２４】
フィルタに蓄積された微粒子の質量ｍsを決定するために、エンジン制御ユニットにおいて式１を実行することは、実行にかかわる他の変数の関数として微粒子の質量ｍsを定義し計算することが現在使用されている自動車のエンジン制御ユニットの能力をはるかに超えた計算能力を要求する点でとくに複雑である。
【００２５】
たとえ、式１がエンジン制御ユニットで実行できたとしても、その結果は全体として満足できないものになろう。本発明者によってなされた机上および車載テストは、実際に、エンジンの種々の動作特性および過去の微粒子の蓄積から独立してフィルタチャネルにおいて一様かつ一定の微粒子分布を仮定することは、実際の動作状態でフィルタに蓄積された微粒子の量を正しく評価することを不可能にしていることを示した。これが、排気ガス流量および温度と、フィルタの圧力降下を測定することに基づく再生制御システムが自動車産業において採用されなかった理由である。
【００２６】
【非特許文献１】コンスタンドポウロス  ア．  ジ．、コストグロウ  エンメ．、スカペルダス  エ．、パパイオアンノウ  エ．、ツァルバレス  デイ．およびクラドポウロウ  エ．共著、「ディーゼル微粒子フィルタの基礎研究：一時的蓄積、再生および経年劣化」、ＳＡＥ２０００−０１−１０１６、２０００
【発明の開示】
【発明が解決しようとする課題】
【００２７】
前述の欠点を解消するように構成され、知られた微粒子量を含む場合に微粒子フィルタの再生を提供する、フィルタに蓄積された微粒子の量を決定する方法を提供することが本発明の目的である。
【課題を解決するための手段】
【００２８】
本発明によれば、請求項１に記載されたとおりの微粒子フィルタに蓄積された微粒子の量を決定する方法が提供される。
【発明の効果】
【００２９】
本発明によれば、微粒子フィルタ内の微粒子分布および微粒子の物理化学的特性がエンジンの動作状態および過去の蓄積に依存して変化することを仮定して、微粒子に蓄積された微粒子の量を計算することができるため、知られたモデルよりはるかに信頼することができ、自動車産業において現在使用されているエンジン中央制御ユニットにおいて実際に実行することができる。
【発明を実施するための最良の形態】
【００３０】
添付図面を参照しながら本発明の好ましい実施の形態を以下に説明するが、これらに限定されるものではない。
【００３１】
図１は、微粒子フィルタと、排気ガス粒量および温度と該微粒子フィルタの圧力降下を測定するための個々の装置とを備えた排気ガスシステムの概略説明図である。
【００３２】
図２は、異なるエンジン動作状態による微粒子フィルタチャネルにおける微粒子分布の変化を示す概略説明図である。
【００３３】
図３は、異なるエンジン動作状態により微粒子フィルタチャネルに蓄積された微粒子の物理化学的特性の変化を示す概略説明図である。
【００３４】
図４は、エンジン動作状態の可能な変化に依存する微粒子フィルタチャネルにおける微粒子の分布変化および物理化学的変化の概略説明図である。
【００３５】
図５は、部分再生（たとえば、ＮＯｘベースの同時再生）後の微粒子フィルタチャネルにおける微粒子の分布変化および物理化学的変化の概略説明図である。
【００３６】
図６は、長いあいだの自動車の不使用（「パーキング」効果）および続いておこる蓄積の後の微粒子フィルタチャネルにおける微粒子の物理化学的特性の変化を示す概略説明図である。
【００３７】
図７は、本発明の微粒子フィルタモデルを示すブロックダイアグラムである。
【００３８】
図８は、図７における「β定義ブロック」のより詳細なブロックダイアグラムである。
【００３９】
図９は、図８におけるステートマシーンのステートダイアグラムである。
【００４０】
図１０は図７における「物理的モデル補正」ブロックのより詳細なブロックダイアグラムである。
【００４１】
図１において参照符号１は、全体としてエンジン２（とりわけ、ディーゼルエンジン）の排気ガスシステムを示している。この例に限られることはないが、図１のディーゼルエンジン２は、過給器（ないしはターボチャージャ）が設けられており、給気導管５に沿って設けられた圧縮機４と該圧縮機４に接続されたタービン６によって特徴づけられたターボチャージャ３を備えている。
【００４２】
排気システム１は、ターボチャージャ３に接近し排気ガス導管７に沿って設けられた酸化触媒コンバータ（以下、「前触媒コンバータ(pre-cat)」という）８と、前触媒コンバータ８の下流に排気ガス導管７に沿って設けられた微粒子フィルタ９と、微粒子フィルタ９の上流に排気ガス導管７に沿って設けられた、さらなる酸化触媒コンバータ（以下、「フロント触媒コンバータ（front-cat）」という）１０とを備えている。
【００４３】
フロント触媒コンバータ１０および微粒子フィルタ９の機能は、触媒微粒子フィルタ（図示されていない）として知られた単一の構成要素によっても実行され得る。
【００４４】
微粒子フィルタ９は、フロント触媒コンバータ１０と別体であれ、一体であれ、微粒子燃焼温度を低減するための化学的な成分が設けられる。
【００４５】
また、排気システム１は、とりわけ各サイクルにおいてエンジンに噴射される全量を決定するための電子制御システム１１を備えている。
【００４６】
より詳しくは、電子制御システム１１は、
給気導管５に沿って設けられ、給気導管５に沿う空気の流れ（ないしは、エアフロー）を示す信号を発生するエアフロー測定装置１２と、前記フロント触媒コンバータ１０の入口に接近して接続された第１および第２入力部と、圧力降下信号を供給する出力部とを有する差動圧量センサ１３であって、本発明に記載された方法によって、微粒子フィルタ９に蓄積した微粒子の量と関係づけられる差動圧力センサ１３と、
微粒子フィルタ９の出口に設けられ、微粒子フィルタ９の出口における排気ガスの温度を示す第１温度信号を出す第１温度センサ１４と、微粒子フィルタ９の入口に設けられ、微粒子フィルタ９の入口における排気ガスの温度を示す第２温度信号を出す第２温度センサ１５と、大気圧を測定するためのセンサ１６と、前述のセンサ１３、１４および１５に接続され、後述する本発明の方法を実行する電子中央制御ユニット１７
とを備えている。
【００４７】
なお、第２温度センサ１５は前触媒コンバータ８の出口に接続してもよい。
【００４８】
本発明は、本発明者による微粒子フィルタモデルの可能性についての深いリサーチの結果によるものであり、知られたモデルよりはるかに信頼することができ、自動車産業において現在使用されているエンジン中央制御ユニットにおいて実際に実行することができる。
【００４９】
より詳しくは、本発明者のリサーチが、知られているモデルが基いている仮定、すなわち、微粒子フィルタチャネル内の微粒子分布および該微粒子の物理化学的特性が、エンジンの動作状態および過去の蓄積から独立して一定であるという仮定が間違っているとの仮定に基いている。
【００５０】
したがって、微粒子フィルタ内の微粒子分布および微粒子の物理化学的特性がエンジンの動作状態および過去の蓄積に依存して変化するものと仮定して、本発明者のつぎの式の定義へと導いた。この式は微粒子フィルタの圧力降下、排気ガス温度および流量、および４つのテストパラメータα、β、γ、δによって、フィルタに蓄積された微粒子の量にかかわるものである。
【００５１】
すなわち、
【００５２】
【数２】
2008286208.tif 000004
【００５３】
ここに、
ΔDPFは圧力降下を示し、Ｐatmは微粒子フィルタの下流における絶対圧力（大雑把には、大気圧に等しいと考えられる）を示し、Ｔは排気ガス温度を示し、Ｑmは排気ガス流量を示し（エンジン給気量および噴射された燃料の全量を加えることによって計算し得る）、
ｍSは微粒子フィルタに蓄積された微粒子量を示し、
α、β、γ、δは後述する４つのテスト動作パラメータを示している。
【００５４】
前記式２から、微粒子フィルタに蓄積された微粒子量が、つぎのとおり容易に計算できる。
すなわち、
【００５５】
【数３】
2008286208.tif 000005
【００５６】
より詳しくは、４つのテストパラメータα、β、γ、δについて、
α、γは、多孔性、孔の大きさなどのフロント触媒コンバータおよび微粒子フィルタの幾何学的特徴と、微粒子フィルタ材料の特性に依存し、
βは微粒子フィルタの幾何学的特徴、チャネル内の微粒子の軸方向および半径方向の分布および微粒子の物理化学的特性（たとえば、密度および透過性）に依存し、
δは排気ガス温度と粘度（典型的には、０．６５から０．７４）とのあいだの補正の指数項である。
【００５７】
微粒子フィルタのそれぞれのタイプについて、パラメータα、β、γ、δの値は、一定の動作速度におけるエンジンの特定の机上のテストによって決定され得る。
【００５８】
より詳しくは、パラメータα、γ、δの値は、きれいな微粒子フィルタ（すなわち、微粒子フィルタには微粒子がない）を用いて前記テストによって決定され、ｍSは零である。一方、パラメータβの値は所定の量の微粒子を含む（たとえば、７ｇ／ｄｍ3および１０ｇ／ｄｍ3）微粒子フィルタを用いて前記テストによって決定され得る。
【００５９】
より詳しくは、前記テストは、実質的に、前記微粒子ローディング条件における微粒子フィルタの圧力降下、排気ガス温度および排気ガス流量の時間パターンを決定することと、パラメータα、β、γ、δの値を算出することを含み、提案した式によって決定されたΔDPFの時間パターンは、テストによって決定されたパターンにできるだけ接近して対応する。
【００６０】
微粒子フィルタモデルは、エンジンの動作状態（たとえば、アイドリング速度、都市、特別の都市および高速道路走行に相当する条件、高トルクおよび高出力）の広い範囲で有効にならなければならない。
【００６１】
提案したモデルは、すでに述べたとおり、パラメータβの値において説明されている微粒子フィルタチャネルにおける軸方向および半径方向微粒子分布に有効である。これは、実際に、４つのテストパラメータのうちの１つだけであり、提案した式において微粒子の蓄積された量ｍSを乗じ、実際には微粒子フィルタの圧力降下ΔDPFと微粒子の蓄積された量ｍSとを関係づけている。
【００６２】
提案した式は、微粒子フィルタモデルを得るために提供しており、微粒子フィルタチャネル内の微粒子分布はエンジン動作状態と微粒子フィルタにおける過去の微粒子蓄積との関数として可変であると仮定される。
【００６３】
より詳しくは、動作パラメータβは、異なるエンジン動作状態（すなわち、とくに一定の微粒子蓄積条件）の関数としてマッピングされる。発明者によってなされたテストは、動作パラメータβと異なるエンジン動作状態とのあいだの関係についてつぎの仮定へと導く。
【００６４】
各微粒子フィルタチャネルの異なる微粒子分布および／または微粒子の異なる物理化学的特性（透過性および密度）は識別することができ、それぞれは特定のエンジン動作状態（すなわち、特定の蓄積条件）に関係し得るものであり、いわゆる「エンジンマップ」上で特定され、エンジン速度および正味平均有効圧力（ＢＭＥＰ）によって定義される。
【００６５】
図２は、各微粒子フィルタチャネルにおける微粒子の一定量の４つの可能な分布、相対的なエンジン動作状態、エンジン動作状態に対応するエンジンマップ領域、およびパラメータβの相対的な性質の値を例にとって示している。
【００６６】
図に示されるように、
エンジンが実質的にエンジンマップの左下領域で動作し、排気ガスが流量および音頭が低い都市の走行条件において、微粒子はチャネル内を一様に蓄積する傾向があり、
エンジンが実質的にエンジンマップの中間領域において動作し、排気ガス流量および温度が「中の低（medium-low）」である第１の特別都市走行条件（Ａ）において、微粒子はチャネルに沿って中ほどまで蓄積する傾向があり、「中の高（medium-high）」の値のパラメータβを与え（微粒子の一定量は微粒子フィルタの「中の低」の圧力降下に対応する）、
エンジンが、実質的にエンジンマップの右下領域において動作し、排ガス流量および温度が「中の高」である、第２の特別都市走行条件（Ｂ）において、微粒子はチャネルの端部において蓄積する傾向があり、「中の低」の値のパラメータβを与え（微粒子の一定量が微粒子フィルタの「中の低」の圧力降下に対応する）、
エンジンがエンジンマップの実質的に上部領域において動作し、排気ガスの流量および温度が「高い」トルクおよび出力走行条件において、微粒子がチャネルの実質的に端部において蓄積する傾向があり、「低い」値のパラメータβを与えている（微粒子の一定量が微粒子フィルタの低圧力降下に対応する）。
【００６７】
同様に、図３は、例示のために、一定量の微粒子の４つの起こりうる物理的特性値、相対的なエンジン動作状態、該エンジン動作状態に対応するエンジンマップ範囲、および相対的性質値のパラメータβを示している。
【００６８】
以上のように、
−都市走行状態では、エンジンは実質的にエンジンマップの左のもっとも下の領域内で動作し、排気ガスの流量および温度は低く、蓄積した微粒子は、低い密度値および透過性値により特徴づけられ、それにより、（一定量の微粒子が微粒子フィルタを通るときの高い圧力低下に対応する）高いパラメータ値βを与える。
−第１の特別都市走行状態（Ａ）では、エンジンは実質的にエンジンマップの中間の領域内で動作し、排気ガスの流量および温度は「中の低」であり、蓄積した微粒子は、「中の低」の密度値および透過性値により特徴づけられ、それにより、（一定量の微粒子が微粒子フィルタを通るときの「中の高」の圧力低下に対応する）「中の高」のパラメータ値βを与える。
−第２の特別都市走行状態（Ｂ）では、エンジンは実質的にエンジンマップの右のもっとも下の領域内で動作し、排気ガスの流量および温度は「中の高」であり、蓄積した微粒子は、「中の高」の密度値および透過性値により特徴づけられ、それにより、（一定量の微粒子が微粒子フィルタを通るときの「中の低」の圧力低下に対応する）「中の低」のパラメータ値βを与える。
−トルクおよびパワー駆動状態では、エンジンは実質的にエンジンマップのもっとも上の領域内で動作し、排気ガスの流量および温度は高く、蓄積した微粒子は、高い密度値および透過性値により特徴づけられ、それにより、（一定量の微粒子が微粒子フィルタを通るときの低い圧力低下に対応する）低いパラメータ値βを与える。
【００６９】
前記想定が与えられれば、テストデータ、関連する微粒子の分布、および物理化学的特性の変化に基づいて、以下の推察がなされ得る。
ａ）排気ガスの流量および温度が増加するにつれて、微粒子は微粒子フィルタのチャネルの終端へ向かって蓄積し、および／または高い密度値および透過性値によって特徴づけられるという傾向がより大きく、そのため、パラメータβ値はより低い。また逆もまた真である。
ｂ）微粒子の分布および物理化学的特性は、排気ガス流量および温度に対して不可逆的に敏感である。
ｂ１）微粒子が低い排気ガス流量および温度状態において蓄積されており、エンジンがフルパワーの動作（高い流量および温度）に切り換わるならば、すでに蓄積された微粒子は、エンジン動作状態に依存している時間には、チャネルの終端へ向かって移動し、および／または高い透過性値および密度値（微粒子がフィルタの壁に押圧される）によって特徴づけられる傾向にある。
ｂ２）微粒子の一部が高い排気ガス流量および温度状態において蓄積されており、かつ低い排気ガス流量および温度の状態に切り換わるならば、新しい状態で蓄積された微粒子の量に依存する微粒子の分布は、両方の状態の全分布に等しい傾向にあり、および／または新しい状態で蓄積された微粒子の量に依存する蓄積された微粒子は、２つの状態のあいだの中間の透過性値および密度値によって特徴づけられる傾向にある。
【００７０】
ｂ１）およびｂ２）における主張は、図４に概略的に示されており、低い排気ガス流量および温度の状態から高い排気ガス流量および温度の状態への通り道、およびその逆方向の通り道が矢印ｂ１およびｂ２によってそれぞれ示される、１つの例示であるエンジンマップを示している。
【００７１】
微粒子フィルタチャネルにおけるパラメータβ値と微粒子分布とのあいだの関係に関する前記想定は、エンジン制御ユニットレベルにおいて、たとえば以下のように実施され得る。
【００７２】
エンジン動作状態、すなわち微粒子フィルタ内の微粒子の安定した蓄積状態に依存する参照パラメータβ値のマトリックス（以下βMAPで示される）は、マップの形態で生成され、記憶される。
【００７３】
現在のエンジン動作状態を示すパラメータβの基準値βMAPは、蓄積された微粒子の量（以下βMOD_OLDで示される）を計算するための式（式２）で記載されたモデルで以前に用いられたパラメータβ値と比較される。
【００７４】
もし、たとえば都市走行状態から特別都市走行状態への切換えの結果として、βMAP＜βMOD_OLDならば、そのとき、エンジン動作状態の関数として変化する特性時間間隔ののちに、蓄積された微粒子の量を計算するための式（式２）で記載されたモデルで、基準値βMAPが用いられる。そのため、βMOD＝βMAPである（微粒子は微粒子フィルタのチャネルの終端へ向かって移動し、および／または圧縮される傾向にあり、すなわち高い透過性値および密度値によって特徴づけられる）。
【００７５】
もし、たとえば特別都市走行状態から都市走行状態への切換えの結果として、βMAP＞βMOD_OLDならば、そのとき、微粒子フィルタ内に蓄積された微粒子の量の最後の評価値（以下、ｍS_OLDで示される）が、記憶され、２つの駆動状態（ｍS−ｍS_OLD）のあいだの切換えののちに蓄積された微粒子の量が所定のしきい値を超えるとき、式（式２）で記載されたモデルで用いられるパラメータβの新しい動作値βMODが以下の式３にしたがって計算される。
【００７６】
【数４】
2008286208.tif 000006
【００７７】
以上のように、ｍS＝２・ｍS_OLDならば、すなわち、新しい状態で蓄積された微粒子の量が以前の状態の場合と同じならば、そのとき、
【００７８】
【数５】
2008286208.tif 000007
【００７９】
（２つの駆動状態の参照パラメータβ値の算術平均が用いられる）。
【００８０】
ｍS＞＞ｍS_OLD、すなわち、新しい状態で蓄積された微粒子の量が以前の状態の場合よりもより一層高いならば、そのとき、βMOD＝βMAPである（新しい駆動状態の参照パラメータβ値が用いられる）。
【００８１】
微粒子フィルタ内に蓄積される微粒子の量は、したがってエンジン動作状態およびまた微粒子フィルタ内の過去の微粒子蓄積の関数として計算される。
【００８２】
このとき、微粒子フィルタの再生が、単に微粒子フィルタ内に蓄積された微粒子の量に基づいて、または微粒子状態インデックスＩPART（すなわち、微粒子分布の均一性、または相対的物理化学的特性の状態の語で定義される）に基づいてより正確に、実行ないしは活性化（activate）される。微粒子状態インデックスＩPARTは、以下の式４にしたがって、パラメータβの動作値βMODの関数として計算される。
【００８３】
【数６】
2008286208.tif 000008
【００８４】
ここで、βMAXおよびβMINは、基準値βMAPマップに記憶された最大および最小のパラメータβ値である。
【００８５】
微粒子分布均一性インデックスＩPARTは、微粒子フィルタ内の微粒子の分布または物理化学的特性の状態が臨界になり、それにより、微粒子フィルタを通るときの圧力低下の効果を奏するまでの範囲の正確な決定を提供する。すなわち、ＩPART値が増加するにしたがって、微粒子フィルタの再生がより臨界的になる。
【００８６】
以下の２つの場合では、たとえば、
１）排気ガスの流量および温度が高い、パワーおよびトルク駆動状態：微粒子は主として微粒子フィルタの終端において集積し、および／または低い動作値のパラメータβを与える。
２）排気ガスの流量および温度が「中の低」である、特別都市走行状態：微粒子は主として微粒子フィルタのチャネルの中間で集積し、および／または「中の低」の動作値のパラメータβを与える。
【００８７】
所定量の蓄積された微粒子のために、これらの２つのエンジン動作状態の第１の場合は、第２の場合よりもさらに一層臨界的であり、ＩPART値に反映される。第１の場合には、実際に、
［外１］
2008286208.tif 000009
一方、第２の場合には
［外２］
2008286208.tif 000010
である。
【００８８】
合計するため、微粒子フィルタの各チャネルにおけるパラメータβ値と微粒子蓄積とのあいだ、または微粒子のパラメータβ値と物理化学的特性とのあいだ、または異なる範囲に対して両方同時に存在する関係が示される。
【００８９】
いいかえれば、エンジン動作状態の変化によって、微粒子フィルタのチャネル内の微粒子蓄積は、変化してもよいが、相対的物理化学的特性は一定のままであり、または実質的に変わらなくてもよく、それは変化する微粒子の物理化学的特性である。
【００９０】
前記関係するマッピング、およびパラメータβの動作値をエンジン動作状態の関数としての計算を行なうものは、両方の場合、またはあいだの場合のいずれかに当てはまる。
【００９１】
部分的な再生（たとえば、ＮＯｘベースの自然発生的な再生）および延長された車両の廃止（prolonged vehicle disuse）（「停車」効果）を考慮するため、特定のアルゴリズムが定義される。部分的な再生または延長された車両の廃止ののちに、微粒子の分布および／または物理化学的特性が、変わることができ、結果として、同じ排気ガス流量および温度における、フィルタに大量に蓄積された同じ微粒子の場合に、フィルタ通過の異なる（概して低い）圧力降下が生じる。それにもかかわらず、本出願人によって実施されたテストは、のちに蓄積された微粒子がフィルタ内にすでに存在する残留物に影響されないことが示された。いいかえれば、のちに蓄積された微粒子は、同じエンジン動作状態（同じ値βによって特徴づけられることを意味する）におけるクリーンフィルタから始まる同じ物理化学的特性および分布によって特徴づけられる。これらの現象は図５〜６によって概略的に示されている。
【００９２】
したがって、部分的な再生（ＮＯｘベースの自然発生的な再生）および延長された車両の廃止（「停車」効果）の効果を考慮することなく評価されたフィルタ内に蓄積された微粒子の質量ｍSを補正（概して増加）するために、式５にしたがって、オフセットΔｍが計算される。
ｍS_CORR＝ｍS＋Δｍ            （式５）
【００９３】
ここで、ｍS_CORRは、フィルタ内に蓄積された微粒子の実質量である。オフセットΔｍは以下の場合に決定される。
【００９４】
１）エンジンのスイッチオンに続くある待ち時間ののち、（もし低ければ）現在の値ｍSが以前のキーオフ（ｍS_CORR_ECU）において、電子制御ユニットに記憶された値ｍS_CORRに対して引かれる：
Δｍ＝ｍS_CORR_ECU−ｍS
それにより、延長された車両の廃止（「停車」効果）ののちに微粒子の物理化学的特性の変化を補償することが可能になる。
【００９５】
２）部分的な再生（たとえば、ＮＯｘベースの自然発生的な再生）ののち、（もし低ければ）現在の値ｍSが現在の値ｍS_CORRに対して引かれる：
Δｍ＝ｍS_CORR−ｍS
現在の値ｍS_CORRが、再生を始める前のｍS_CORRの値と再生された微粒子の実験的に評価された値ｍREG（時間とエンジン動作状態との関数）とのあいだの差として評価される。それにより、部分的な再生（たとえば、ＮＯｘベースの自然発生的な再生）ののちの微粒子分布の変化を補償することが可能になる。
【００９６】
オフセットΔｍが評価されると、現在のｍSの値に対して加えられ、そのため、ｍS_CORR＝ｍS＋Δｍになる。
【００９７】
最後に、主として制限されたセンサの精度（アイドリング状態における測定は、通常、信頼性なしとみなされる）の効果を制限するため、フィルタを通る排気ガス流量および圧力低下の値がある許容値よりも高い場合のみ、ｍS_CORRの値が信頼性ありとみなされる。フィルタを通る排気ガスの流量および圧力低下の値がある許容値より低ければ、シミュレーション状態が始められ（Ｓｉｍ＝オン）、最終の許容値ｍS_CORRは、毎回の工程で、微粒子の排出に評価された瞬間的な値ｍEXP（エンジン動作状態の関数）が増加される。
【００９８】
図７に示されるブロックダイヤグラムは、微粒子フィルタ内に蓄積された微粒子の量の評価およびフィルタの再生の実行に関連することに関して、すでに述べられたことについてまとめる。
【００９９】
さらに詳細にいえば、図７において、符号２０は、以前に示された式（式２）に記載された微粒子フィルタの物理モデルを実施するためのブロックを示し、微粒子フィルタ内に蓄積された微粒子の量を供給する。符号２１は、ブロック２０内で用いられ、現在のエンジン動作状態および微粒子フィルタ内の過去の微粒子蓄積に相対的なパラメータβの動作値βMODを供給するブロックを示す。符号２２は、微粒子状態インデックスＩPARTを計算するブロックを示す。
【０１００】
さらに詳細にいえば、ブロック２０は、微粒子フィルタを通るときの圧力低下ΔＰDPF、排気ガス温度Ｔ、排気ガス流量Ｑｍ、大気圧Ｐatm、およびブロック２１によって供給されるパラメータβの動作値βMODを受け取り、現在のエンジン動作状態における微粒子フィルタ内に蓄積された微粒子量ｍSを供給し、当該微粒子量ｍSは、叙上の式（式２）を用いて計算される。
【０１０１】
ブロック２１は、現在のエンジン動作状態を定義するエンジンスピードおよびエンジン負荷を受け取り、ブロック２０から微粒子フィルタ内に蓄積された微粒子量ｍSを受け取り、パラメータβの動作値βMODを、蓄積された微粒子量ｍSを計算するためのブロック２０および微粒子状態インデックスＩPARTを計算するためのブロック２１へ供給する。
【０１０２】
ブロック２２は、ブロック２１からのパラメータβの動作値βMODを受け取り、微粒子フィルタ内の微粒子の状態インデックスＩPARTを供給する。
【０１０３】
最後に、ブロック２５は、微粒子フィルタ内に蓄積された微粒子量ｍSを受け取り、部分的な再生（ＮＯｘベースの自然発生的な再生）および延長された車両の廃止（「停車」効果）の効果を考慮することなく評価され、フィルタ内に蓄積された微粒子の最終的に補正された質量ｍS_CORRを供給する。
【０１０４】
図８〜９には、微粒子フィルタ内に蓄積された微粒子量を評価するために用いられるパラメータβの動作値βMODを計算する図７のブロック３２内で行なわれる動作を示すブロックダイヤグラムおよび状態ダイヤグラムを示す。
【０１０５】
さらに詳細にいえば、図８において、符号２３は、パラメータβの基準値βMAPのマップをエンジン動作状態の関数として記憶するブロックを示し、該ブロック２３は、エンジンスピードおよびエンジン負荷を受け取り、現在のエンジン状態に関連するパラメータβの基準値βMAPを提供する。符号２４は、パラメータβの動作値βMODを計算するためのステートマシン（state machine）を示し、該ステートマシン２４は、ブロック２３から基準値βMAPおよび記憶された微粒子量ｍS_OLDを受け取り、を動作値βMODを供給する。
【０１０６】
ステートマシン２４は、図９における状態プログラムを参照しながら以下に記載された動作を実行する。
【０１０７】
図９に示されるように、エンジンのターンオンまたは点火は、βMOD＝βECU、βMOD_OLD＝βECU、およびｍS＝ｍS_ECUである状態Ａに対応する。ここで、βECUおよびｍS_ECUは、最後にターンオフ（エンジンが切られる）前のモデルに用いられ、エンジン制御ユニットに記憶されたパラメータβおよび蓄積された微粒子の質量ｍSの値である。
【０１０８】
もし、βMAP≦βMOD_OLD_LDであれば、そのとき、状態Ａは、ｍS_OLD＝ｍSを入力（エンター）することによって蓄積された微粒子の最後に評価された値ｍSが記憶され、かつβMOD＝βMAPおよびβMOD_OLD＝βMAPが入力された状態Ｂへ移行する。ところが、βMAP＞βMOD_OLDであれば、そのとき、状態Ａは、スタンドバイ状態Ｃへ移行する。状態Ｂにおいて、βMAP＞βMOD_OLDであれば、状態Ｂは状態Ｃへ移行するが、βMAP≦βMOD_OLDであれば、状態Ｂを維持する。
【０１０９】
状態Ｃでは、βMAP≦βMOD_OLDであれば、状態Ｃは状態Ｂへ移行するが、計算された、蓄積された微粒子ｍSは、記憶された、蓄積された微粒子ｍS_OLDを所定のしきい値だけ超えれば、すなわち、ｍS−ｍS_OLD≧（しきい値）であれば、状態Ｃは、パラメータβの値βMODが以前に示された式（式３）にしたがって計算され、かつβMOD_OLD＝βMODおよびｍS_OLD＝ｍSが入力された状態Ｄへ移行する。
【０１１０】
状態Ｄでは、βMAP＞βMOD_OLDであれば、状態Ｄは状態Ｃへ戻るが、βMAP≦βMOD_OLDであれば、状態Ｄは状態Ｂへ移行する。
【０１１１】
エンジンがターンオフまたは切られたとき、すべての状態Ａ、Ｂ、Ｃ、Ｄはいずれも状態Ｅへ移行し、状態Ｅでは、現在のβMODおよびｍSの値がエンジン制御ユニットにそれぞれβECUおよびｍS_ECUとして記憶され、そのため、それらは、エンジンをターンオンまたは点火した次のときに利用可能である。
【０１１２】
図１０は、フィルタに蓄積された微粒子の質量の最終的な補正値ｍS_CORRを計算するための図７のブロック２５で実行される動作を示す状態ダイアグラムを示している。
【０１１３】
図１０に示されるように、エンジンのターンオンまたは点火は、ｍS_CORR＝ｍS_CORR_ECUである状態Ａ状態Ａに対応する。ここで、ｍS_CORR_ECUは、最後にターンオフまたはエンジンが切られる前に、エンジン制御ユニットに記憶された、蓄積された微粒子の質量ｍS_CORRの値である。エンジンのスイッチオンに引き続くある待ち時間（状態Ｂ）ののち、すなわち、タイマが数値ｔに達したとき、現在のｍSの値が、以前のキーオフ（またはエンジンが切られた）のときに電子制御ユニットに記憶されたｍS_CORRの値ｍS_CORR_ECUに対して（もし低ければ）引かれる：Δｍ＝ｍS_CORR−ｍS（状態Ｃ）。
【０１１４】
オフセットΔｍは、通常の状態において、引き続く値ｍSに対して加えられる：ｍS_CORR＝ｍS＋Δｍ（状態Ｄ）。ステートマシンのこの部分は、２つの引き続くキーオフ（またはエンジンを切る）およびキーオン（またはエンジンを点火する）（すなわち、車両の停車）のあいだ、微粒子フィルタ内の微粒子の分布の起こりうる変化、および／または微粒子の物理化学的特性における変化を考慮に入れることができる。
【０１１５】
ｍSの値が信頼性ありとみなされる場合のみ、すなわち、変数Ｓｉｍ＝オフの場合に、オフセット補正（状態Ｄ）は実行される。それは、フィルタを通る排気ガスの流量および圧力低下の値がある許容値よりも高いことを意味する。Ｓｉｍ＝オンならば、シミュレーション状態（状態Ｆ）が始まり、ステートマシンの出力ｍS_CORRは、毎回の工程において、微粒子排出の実験的に評価された瞬間的な値ｍEXP（エンジン動作状態の関数）が増加される。
【０１１６】
ＮＯｘベースの自然発生の再生状態が認識され、すなわちＮＯｘ＿ｒｅｇ＝オン（ＮＯｘ  ｒｅｇはエンジン動作状態の関数である）であれば、状態Ｅが始まり、ステートマシンの出力ｍS_CORRは、毎回の工程において、微粒子再生の実験的に評価された瞬間的な値ｍNOx（エンジン動作状態の関数）が減らされる。そのとき、ｍS_CORRは、その値とｍS＋Δｍ（すなわち、状態ＤのｍS_CORRの値）との最大値として計算される。
【０１１７】
ＮＯｘベースの自然発生の再生、すなわち、ＮＯｘ＿ｒｅｇ＝オフののち、現在の値ｍSは、現在の値ｍS_CORRに対して（低ければ）引かれ、オフセット値を計算する：Δｍ＝ｍS_CORR−ｍS（状態Ｃ）。
【０１１８】
ステートマシンのこの部分は、フィルタのＮＯｘベースの自然発生の再生により、微粒子フィルタ内の微粒子の空間的な分布の起こりうる変化、および／または微粒子の物理化学的特性における変化を考慮に入れることができる。
【０１１９】
キーオフ（またはエンジンを切る）（状態Ｇ）のとき、ｍS_CORRの値は、エンジン制御ユニット内に記憶される（ｍS_CORR_ECU＝ｍS_CORR）。
【０１２０】
本発明の利点は、以下の記載から明らかになる。
【０１２１】
本発明は、とりわけ、自動車産業で一般的に採用されるエンジン制御ユニットを用いて、微粒子フィルタ内に蓄積された微粒子の量を計算することを提供し、すなわち、高出力でそれゆえ高いコストの計算手段の必要がない。
【０１２２】
その結果、本発明は、（蓄積された微粒子の実際の量または微粒子状態インデックスに基づいて）前記微粒子フィルタの再生の評価基準の採用を可能とする。この評価基準は、現在用いられている評価基準（車両マイル距離）よりも効果的であり、絶対に必要なときの微粒子フィルタの再生のみを提供し、そのため、微粒子フィルタの再生に一般的に関連する燃料消費および満載の性能という点における欠点を最小にする。
【０１２３】
出願人によって実行されたテストは、実際に、微粒子フィルタの１０秒の再生は２００％程度まで燃料消費を増加するが、本明細書で示される微粒子フィルタの再生は平均的な燃料消費の１〜２％のみの上昇する結果となることを示す。
【０１２４】
また、提案される微粒子フィルタのモデルは、単純な分析式に基づいており、一定スピードで動作するエンジンを用いた特定のセットのベンチテストによって、容易に正当であると確認され、温度および圧力低下センサは例外であるが、車両に普通設けられるセンサ以外の追加のセンサの必要なく、実施され得る。
【０１２５】
なお、添付されている請求項に定義されているような本発明の範囲から逸脱することなく、本明細書および図面に記載された事項を変更してもよいことは明らかである。
【図面の簡単な説明】
【０１２６】
【図１】微粒子フィルタと、排気ガス粒量および温度と該微粒子フィルタの圧力降下を測定するための個々の装置とを備えた排気ガスシステムの概略説明図である。
【図２】異なるエンジン動作状態による微粒子フィルタチャネルにおける微粒子分布の変化を示す概略説明図である。
【図３】異なるエンジン動作状態により微粒子フィルタチャネルに蓄積された微粒子の物理化学的特性の変化を示す概略説明図である。
【図４】エンジン動作状態の可能な変化に依存する微粒子フィルタチャネルにおける微粒子の分布変化および物理化学的変化の概略説明図である。
【図５】部分再生（たとえば、ＮＯｘベースの同時再生）後の微粒子フィルタチャネルにおける微粒子の分布変化および物理化学的変化の概略説明図である。
【図６】長いあいだの自動車の不使用（「パーキング」効果）および続いておこる蓄積の後の微粒子フィルタチャネルにおける微粒子の物理化学的特性の変化を示す概略説明図である。
【図７】本発明の微粒子フィルタモデルを示すブロックダイアグラムである。
【図８】図７における「β定義ブロック」のより詳細なブロックダイアグラムである。
【図９】図８におけるステートマシーンのステートダイアグラムである。
【図１０】図７における「物理的モデル補正」ブロックのより詳細なブロックダイアグラムである。
【図１】
2008286208.tif 000011
【図２】
2008286208.tif 000012
【図３】
2008286208.tif 000013
【図４】
2008286208.tif 000014
【図５】
2008286208.tif 000015
【図６】
2008286208.tif 000016
【図７】
2008286208.tif 000017
【図８】
2008286208.tif 000018
【図９】
2008286208.tif 000019
【図１０】
2008286208.tif 000020