(57)【要約】【課題】ディーゼルパティキュレートフィルタ（ＤＰＦ）の再生時期を簡単かつ高精度に判じることができるようにすることである。【解決手段】ＤＰＦ３２の差圧ΔＰを検出する差圧センサ５４と、排気通路３の流速ｖを検出するエアフローメータ５２とを設け、ＥＣＵ５１で、排気微粒子堆積量ＭＬ　を、発明者らの知見に基づいてなされた次式ＭＬ　＝［ΔＰ−（Ａμｖ＋Ｃρｖ２　］／（Ｂμｖ＋Ｄρｖ２　）により算出し、算出堆積量ＭＬ　を基準の堆積量と比較する判定式に基づいて再生するか否かを判じるようにする。マップのような大規模化する手段によることなく高精度に堆積量を求められる。【選択図】　　　　図１
2004019523.tif 000002
【特許請求の範囲】
【請求項１】
排気通路の途中に、排気微粒子を捕集するパティキュレートフィルタを有する内燃機関の排気ガス浄化装置において、
前記パティキュレートフィルタの入口側と出口側との間の差圧を検出する差圧検出手段と、
前記パティキュレートフィルタを流通する排気ガスの流速を検出する流速検出手段と、
前記パティキュレートフィルタに堆積した排気微粒子の堆積量の多少を検出差圧と検出流速とに基づいて判定する判定式にしたがって、前記パティキュレートフィルタを再生するか否かを決定する再生決定手段とを具備し、
前記判定式は、前記差圧をΔＰ、前記流速をｖ、前記パティキュレートフィルタを流通する排気ガスの粘度をμ、前記排気ガスの密度をρとし、定数Ｍ，Ｎを用いて、ΔＰ≧Ｍμｖ＋Ｎρｖ２　と等価であり、これを満たしたときに前記パティキュレートフィルタの再生が許容されることを特徴とする内燃機関の排気ガス浄化装置。
【請求項２】
請求項１記載の内燃機関の排気ガス浄化装置において、前記再生決定手段は、堆積量を、これをＭＬ　とし、定数Ａ，Ｂ，Ｃ，Ｄを用いて、
ＭＬ　＝［ΔＰ−（Ａμｖ＋Ｃρｖ２　）］／（Ｂμｖ＋Ｄρｖ２　）により算出するように設定し、かつ、前記判定式を、算出された堆積量を予め設定した堆積量と大小を比較する判定式とした内燃機関の排気ガス浄化装置。
【請求項３】
請求項２記載の内燃機関の排気ガス浄化装置において、前記再生決定手段は、前記堆積量の算出式が、定数Ａ，Ｂ，Ｃ，Ｄがそれぞれ異なる値をとる２種類の算出式であって、一方の算出式が予め設定した基準値以下の範囲の堆積量に対して適合し、他方の算出式が前記基準値以上の範囲の堆積量に対して適合する算出式が記憶されており、
算出式により算出された堆積量がその算出式が適合する堆積量の範囲を外れると、堆積量を残りの算出式により算出し直すように設定した内燃機関の排気ガス浄化装置。
【請求項４】
請求項２または３いずれか記載の内燃機関の排気ガス浄化装置において、前記パティキュレートフィルタが、これに堆積した排気微粒子を酸化燃焼する酸化触媒を有するパティキュレートフィルタであり、
前記パティキュレートフィルタの再生中に、算出堆積量の低下状態を検出する堆積量低下状態検出手段と、
算出堆積量の低下が遅いほど、前記酸化触媒が劣化しているものと判定する触媒劣化判定手段とを具備せしめた内燃機関の排気ガス浄化装置。
【請求項５】
請求項１ないし４いずれか記載の内燃機関の排気ガス浄化装置において、パティキュレートフィルタを流通する排気ガスの温度を検出する温度検出手段を具備せしめ、
前記再生決定手段を、予め記憶した排気ガスの粘度と温度との対応関係に基づいて、検出温度から前記粘度μを求めるように設定した内燃機関の排気ガス浄化装置。
【請求項６】
請求項１ないし５いずれか記載の内燃機関の排気ガス浄化装置において、前記パティキュレートフィルタを流通する排気ガスの温度を検出する温度検出手段を具備せしめ、
前記再生決定手段を、予め記憶した排気ガスの密度と温度との対応関係に基づいて、検出温度から前記密度ρを求めるように設定した内燃機関の排気ガス浄化装置。
【請求項７】
請求項１ないし６いずれか記載の内燃機関の排気ガス浄化装置において、前記流速検出手段には、
内燃機関の吸入空気量を検出する吸入空気量検出手段と、吸入空気量に、噴射燃料の燃焼で増大する体積分を加算して、これを、前記パティキュレートフィルタを流通する排気ガスの体積流量とする補正手段とを具備せしめ、体積流量から流速に換算するようにした内燃機関の排気ガス浄化装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、内燃機関の排気ガス浄化装置に関し、特にパティキュレートフィルタを再生する技術に関する。
【０００２】
【従来の技術】
近年、自動車等に搭載される内燃機関では、排気エミッションの向上が要求されており、特に軽油を燃料とする圧縮着火式のディーゼルエンジンでは、ＣＯ、ＨＣ、ＮＯｘ　に加え、排気ガス中に含まれる煤やＳＯＦ等の排気微粒子を除去することが必要になる。このため、排気通路にパティキュレートフィルタを配置し、ここで、排気ガス中の排気微粒子を捕集している。
【０００３】
パティキュレートフィルタは、流入した排気ガスに多孔質の隔壁を透過させ、その際に、隔壁の表面や細孔で排気ガス中の排気微粒子を捕集する。捕集されて堆積する量が過剰に増えると、パティキュレートフィルタにおける流通抵抗の増大で内燃機関の背圧が上昇し、出力の低下等をもたらす。このため、パティキュレートフィルタに捕集された排気微粒子をパティキュレートフィルタから適宜、除去してパティキュレートフィルタを再生し、パティキュレートフィルタの排気微粒子捕集能力を回復させる必要がある。
【０００４】
パティキュレートフィルタの再生を内燃機関の運転中に可能としたものとして、パティキュレートフィルタに白金等の酸化触媒を設けて、酸化触媒の酸化作用を利用したものがある。このものでは、例えば排気行程において燃料を噴射するポスト噴射により燃料をパティキュレートフィルタに供給し、その燃焼熱を利用して、噴射燃料に比して酸化しにくい堆積排気微粒子を酸化、除去する。
【０００５】
パティキュレートフィルタの再生は頻繁に行うと燃費が悪化し、一方、次に再生するまでの間が空きすぎると、排気微粒子の堆積量が過剰で、再生処理において排気微粒子が急激に燃焼し、パティキュレートフィルタが異常な高温となり、破損するおそれがある。このため、排気微粒子の堆積量の多少を判断し、再生時期を決定するのが望ましい。特開平７−３３２０６５号には、パティキュレートフィルタへの排気微粒子の堆積量の増大による前記通気抵抗の増大で、パティキュレートフィルタの入口と出口との間の差圧が増大することを利用して、この差圧を検出し、検出差圧が所定値を越えると再生すべき時期と判じるものが開示されている。
【０００６】
【発明が解決しようとする課題】
しかしながら、前記特開平７−３３２０６５号の技術では、内燃機関の運転状態によって、同じ大きさの差圧が検出されても排気微粒子の堆積量が異なり、必ずしも十分な精度で排気微粒子の堆積量の多少を判断することができない。運転状態別の詳細なマップを記憶しておくことが考えられるが、大きな記憶容量が必要になり、簡単ではない。
【０００７】
本発明は前記実情に鑑みなされたもので、簡単な構成で、再生時期を適正に決定し得る内燃機関の排気ガス浄化装置を提供することを目的とする。
【０００８】
【課題を解決するための手段】
発明者らは、パティキュレートフィルタに関し、排気微粒子の堆積と、これが排気ガスの流通におよぼす影響について鋭意実験研究を重ねた結果、パティキュレートフィルタの入口側と出口側との間の差圧、すなわちパティキュレートフィルタにおける圧力損失をΔＰ、パティキュレートフィルタを流通する排気ガスの流速をｖ、排気ガスの粘度をμ、排気ガスの密度をρとすると、モデル式　ΔＰ＝Ｍμｖ＋Ｎρｖ２　の関係を満たし、係数Ｍ，Ｎが排気微粒子の堆積量が多いほど、大きな値をとることがわかった。詳しくは、Ｍ，Ｎは、実質的に、堆積量ＭＬ　の一次関数形をとり、堆積量がある値を越えると、堆積量ＭＬ　に対するＭ，Ｎの変化率が小さくなる。この、ある値を挟む前後で異なる傾向を示すのは、パティキュレートフィルタにおける圧力損失が、最初は、パティキュレートフィルタの細孔のうちの、排気微粒子が詰まっている細孔の割合で増大するのに対して、細孔が排気微粒子により略詰まってしまうと、排気微粒子の堆積層の厚さに応じて増大するためと認められる。本発明はかかる知見に基づきなされたものである。
【０００９】
請求項１記載の発明では、排気通路の途中に、排気微粒子を捕集するパティキュレートフィルタを有する内燃機関の排気ガス浄化装置において、
前記パティキュレートフィルタの入口側と出口側との間の差圧を検出する差圧検出手段と、
前記パティキュレートフィルタを流通する排気ガスの流速を検出する流速検出手段と、
前記パティキュレートフィルタに堆積した排気微粒子の堆積量の多少を検出差圧と検出流速とに基づいて判定する判定式にしたがって、前記パティキュレートフィルタを再生するか否かを決定する再生決定手段とを具備する構成とする。
かつ、前記判定式は、前記差圧をΔＰ、流速をｖ、前記パティキュレートフィルタを流通する排気ガスの粘度をμ、前記排気ガスの密度をρとし、定数Ｍ，Ｎを用いて、ΔＰ≧Ｍμｖ＋Ｎρｖ２　と等価であり、これを満たしたときに前記パティキュレートフィルタの再生が許容されるものである。
【００１０】
前記のごとく、パティキュレートフィルタにおける排気ガスの流通に関し、ΔＰ＝Ｍμｖ＋Ｎρｖ２　の関係が成り立つので、Ｍ，Ｎを基準の堆積量に対して求めておけば、前記判定式から排気微粒子の堆積量の多少について正確に判断することができる。しかも、式を記憶しておけばよいので、記憶するデータのサイズがマップのように大規模化しない。
【００１１】
請求項２記載の発明では、請求項１の発明の構成において、前記再生決定手段は、堆積量を、これをＭＬ　とし、定数Ａ，Ｂ，Ｃ，Ｄを用いて、
ＭＬ　＝［ΔＰ−（Ａμｖ＋Ｃρｖ２　）］／（Ｂμｖ＋Ｄρｖ２　）により算出するように設定し、かつ、前記判定式を、算出された堆積量を予め設定した堆積量と大小を比較する判定式とする。
【００１２】
前記のごとくＭ，Ｎが堆積量ＭＬ　の一次関数形をとるので、前記モデル式　ΔＰ＝Ｍμｖ＋Ｎρｖ２　において、前記ＭＬ　の算出式が成り立つ。排気微粒子の堆積量が客観的に知られる。
【００１３】
請求項３記載の発明では、請求項２の発明の構成において、前記再生決定手段は、前記堆積量の算出式が、定数Ａ，Ｂ，Ｃ，Ｄがそれぞれ異なる値をとる２種類の算出式であって、一方の算出式が予め設定した基準値以下の範囲の堆積量に対して適合し、他方の算出式が前記基準値以上の範囲の堆積量に対して適合する算出式が記憶されており、
算出式により算出された堆積量がその算出式が適合する堆積量の範囲を外れると、堆積量を残りの算出式により算出し直すように設定する。
【００１４】
パティキュレートフィルタにおける圧力損失に対する支配要因が変わっても、より適切な判定式が選択されるので、より正確に排気微粒子の堆積量が知られる。
【００１５】
請求項４記載の発明では、請求項２または３の発明の構成において、前記パティキュレートフィルタが、これに堆積した排気微粒子を酸化燃焼する酸化触媒を有するパティキュレートフィルタであり、
前記パティキュレートフィルタの再生中に、算出堆積量の低下状態を検出する堆積量低下状態検出手段と、
算出堆積量の低下が遅いほど、前記酸化触媒が劣化しているものと判定する触媒劣化判定手段とを具備せしめる。
【００１６】
酸化触媒が劣化しその酸化能力が低下すると、パティキュレートフィルタの再生中において、堆積した排気微粒子の燃焼速度が鈍るため、堆積量の低下速度が鈍る。したがって、算出堆積量の低下状態に基づいて、パティキュレートフィルタの故障の有無や劣化を判断することができる。
【００１７】
請求項５記載の発明では、請求項１ないし４の発明の構成において、パティキュレートフィルタを流通する排気ガスの温度を検出する温度検出手段を具備せしめ、
前記再生決定手段を、予め記憶した排気ガスの粘度と温度との対応関係に基づいて、検出温度から前記粘度μを求めるように設定する。
【００１８】
排気ガスの粘度μがより正確に知られるので、より適切な判定式を与えることができる。
【００１９】
請求項６記載の発明では、請求項１ないし５の発明の構成において、前記パティキュレートフィルタを流通する排気ガスの温度を検出する温度検出手段を具備せしめ、
前記再生決定手段を、予め記憶した排気ガスの密度と温度との対応関係に基づいて、検出温度から前記密度ρを求めるように設定する。
【００２０】
排気ガスの密度ρがより正確に知られるので、より適切な判定式を与えることができる。
【００２１】
請求項７記載の発明では、請求項１ないし６の発明の構成において、前記流速検出手段には、
内燃機関の吸入空気量を検出する吸入空気量検出手段と、吸入空気量に、噴射燃料の燃焼で増大する体積分を加算して、これを、前記パティキュレートフィルタを流通する排気ガスの体積流量とする補正手段とを具備せしめ、体積流量から流速に換算するようにする。
【００２２】
噴射燃料が気筒内で燃焼すると、燃焼に供された燃料の量に応じて、排気ガスの体積流量が、吸入空気量よりも増大する。したがって、排気微粒子を含む排気ガスを直接、測定することなく、排気ガスの体積流量を正確に求めることができる。パティキュレートフィルタを流通する排気ガスの流速は、パティキュレートフィルタの形状などで規定される係数で比例関係にあるから、体積流量から換算することができる。
【００２３】
気筒からの、排気微粒子を含む排気ガスの流量を直接、検出するとすれば、検出用のセンサが排気微粒子による汚れで検出特性への影響があるが、比較的清浄な吸入空気の流量を検出すればよいから、検出特性への影響がなく、排気ガスの流速について信頼性の高い検出が可能となる。
【００２４】
また、通常の内燃機関には吸入空気量検出手段であるエアフローメータが付設されており、ここで検出された吸入空気量をそのまま用いることができるから、構成が簡単である。
【００２５】
【発明の実施の形態】
図１に本発明を適用した第１実施形態になるディーゼルエンジンの構成を示す。ディーゼルエンジンは、４気筒を備えたエンジン本体１に、吸気通路２の最下流部である吸気マニホールド２１と、排気通路３の最上流部である排気マニホールド３１とが接続され、排気通路３は、排気マニホールド３１の集合部にパティキュレートフィルタ３２が連なっている。パティキュレートフィルタ３２は、コーディエライトや炭化珪素等の多孔質セラミック製のハニカム構造体の流路を目封じしてフィルタ本体４を形成したもので、入口３２ａから流入したエンジン本体１の各気筒からの排気ガスが、多孔質の隔壁を透り、出口３２ｂから下流へと流れていく。このとき、パティキュレートフィルタ３２には、排気ガスに含まれる排気微粒子が捕集され、走行距離に応じて堆積していく。また、パティキュレートフィルタ３２のフィルタ本体４の表面には白金やパラジウム等の貴金属を主成分とする酸化触媒が担持されており、所定の温度条件下で排気微粒子を酸化、燃焼し、除去する。
【００２６】
エンジン本体１のインジェクタ等、エンジン各部を制御するＥＣＵ５１が設けられている。
【００２７】
ＥＣＵ５１には、運転状態を示す種々の信号が入力している。この中には、パティキュレートフィルタ３２に堆積する排気微粒子の堆積量を知るための信号も含まれており、そのためのセンサが設けられている。すなわち、排気通路３には管壁を貫通して温度検出手段である温度センサ５３ａ，５３ｂが設けてあり、排気温度を検出するようになっている。温度センサ５３ａ，５３ｂはパティキュレートフィルタ３２の直上流と直下流とのそれぞれに設けられている。温度センサ５３ａの検出温度は、パティキュレートフィルタ３２の入口３２ａにおける、流通する排気ガスの温度であり、以下、ＤＰＦ入口温度という。温度センサ５３ｂの検出温度は、パティキュレートフィルタ３２の出口３２ｂにおける、流通する排気ガスの温度であり、以下、ＤＰＦ出口温度という。
【００２８】
また、排気通路３には、パティキュレートフィルタ３２の直上流側で分岐する第１の分岐通路３３ａと、パティキュレートフィルタ３２の直下流側で分岐する第２の分岐通路３３ｂとが接続され、両分岐通路３３ａ，３３ｂに介設された差圧検出手段である差圧センサ５４が、パティキュレートフィルタ入口３２ａとパティキュレートフィルタ出口３２ｂとの差圧を検出するようになっている。この差圧はパティキュレートフィルタ３２における圧力損失を示している。
【００２９】
また、吸気通路２には吸入空気量検出手段であるエアフローメータ５２が設けられ、吸入空気量を検出するようになっている。
【００３０】
その他、ＥＣＵ５１に、アクセル開度、冷却水温等の運転状態を示すパラメータが入力しているのは勿論である。
【００３１】
ＥＣＵ５１はマイクロコンピュータを中心に構成された一般的な構成のもので、そのＲＯＭには、内燃機関各部を制御するための制御プログラムの他、パティキュレートフィルタ３２の排気微粒子の堆積量を算出するプログラムや、この算出プログラムで用いられる判定式を特定する情報が格納されている。判定式については後述する。算出された堆積量に基づいてパティキュレートフィルタ３２を再生するか否かが判断される。
【００３２】
図２に、ＥＣＵ５１で実行されるパティキュレートフィルタ３２の再生に関する制御内容を示す。先ずステップＳ１０１で、ＤＰＦ温度、吸入空気量、圧力損失を取り込む。ＤＰＦ温度は、ＤＰＦ入口温度とＤＰＦ出口温度より演算処理して求める。ＤＰＦ入口温度は変動が大きいので一次遅れフィルタによる処理を経たものを用いるのがよい。吸入空気量はここでは質量流量である。
【００３３】
ステップＳ１０２は、ＥＣＵ５１の、補正手段としての処理で、ＥＣＵ５１はエアーフローメータ５２とともに流速検出手段を構成する。ステップＳ１０２では、先ず、排気ガスの体積流量を算出する。算出は式（１）にしたがって算出する。なお、エアーフローメータ５２は、吸気マニホールド２１の上流位置に配されている。
体積流量（ｍ３　／ｓｅｃ）
＝［吸入空気量（ｇ／ｓｅｃ）／２８．８（ｇ／ｍｏｌ）］
×２２．４×１０−３（ｍ３　／ｍｏｌ）
×［ＤＰＦ温度（Ｋ）／２７３（Ｋ）］
×［大気圧（ｋＰａ）／（大気圧（ｋＰａ）＋差圧（ｋＰａ））］
＋燃料噴射量（ｃｃ／ｓｅｃ）／２０７．３（ｇ／ｍｏｌ）］
×０．８４（ｇ／ｃｃ）×６．７５（ｍｏｌ）
×２２．４×１０−３（ｍ３　／ｍｏｌ）
×［大気圧（ｋＰａ）／（大気圧（ｋＰａ）＋差圧（ｋＰａ））］・・・（１）
【００３４】
式（１）の第１項は吸入吸気量を質量流量から体積流量に変換したものである。第２項は、噴射燃料の燃焼で、吸入空気量に対して増大する排気ガスの増量分である。第２項中、０．８４（ｇ／ｃｃ）は軽油の代表的な液密度である。６．７５（ｍｏｌ）は燃料噴射量１（ｍｏｌ）に対する増量分のモル数である。
【００３５】
増量分（６．７５（ｍｏｌ））については以下のように得ている。軽油の組成は代表的には、Ｃ１５Ｈ２７．３（分子量２０７．３）と表され、燃焼時の反応式は次の通りである。
【００３６】
Ｃ１５Ｈ２７．３＋２１．７５Ｏ２　→１５ＣＯ２　＋１３．５Ｈ２　Ｏ
【００３７】
したがって、燃料噴射量１ｍｏｌに対し、６．７５（＝（１５＋１３．５）−２１．７５）ｍｏｌのモル数の増加がある。
【００３８】
また、燃料噴射はＥＣＵ５１で決定される所定の噴射時期にのみ噴射され、間欠的な噴射となる。式（１）中の燃料噴射量は、非噴射期間も合わせた平均的な燃料噴射量である。
【００３９】
排気ガスの体積流量は、パティキュレートフィルタ３２の有効通路面積で除して、流速に換算される。
【００４０】
ステップＳ１０３では、排気ガスの粘度μをＤＰＦ温度に基づいて算出する。これは、所定の演算式若しくはマップに基づいて行う。図３に排気ガスの粘度と温度の対応関係を示す。
【００４１】
ステップＳ１０３ではまた、排気ガスの密度ρをＤＰＦ温度に基づいて算出する。これは、所定の演算式若しくはマップに基づいて行う。図４に排気ガスの密度と温度の対応関係を示す。
【００４２】
ステップＳ１０４では、排気微粒子の堆積量（以下、適宜、ＰＭ堆積量という）を算出する。ＰＭ堆積量の算出式として、次の式（２−１）、式（２−２）をＲＯＭに格納している。式中、ＭＬ　は堆積量、ΔＰは圧力損失、ｖは流速である。式（２−１）のＡ１　，Ｂ１　，Ｃ１　，Ｄ１　、式（２−２）のＡ２　，Ｂ２　，Ｃ２　，Ｄ２　は定数である。
【００４３】
ＭＬ　＝［ΔＰ−（Ａ１　μｖ＋Ｃ１　ρｖ２　）］／（Ｂ１　μｖ＋Ｄ１　ρｖ２　）
・・・（２−１）
ＭＬ　＝［ΔＰ−（Ａ２　μｖ＋Ｃ２　ρｖ２　）］／（Ｂ２　μｖ＋Ｄ２　ρｖ２　）
・・・（２−２）
【００４４】
式（２−１）、式（２−２）について説明する。
【００４５】
発明者らは、パティキュレートフィルタについて排気微粒子の堆積と、これが排気ガスの流通におよぼす影響について鋭意実験研究を重ねた結果、次の知見を得た。パティキュレートフィルタ３２における圧力損失ΔＰと流速ｖとの間には、モデル式として式（３）の関係があり、係数Ｍ，ＮがＰＭ堆積量が多いほど、大きな値をとることがわかった。
ΔＰ＝Ｍμｖ＋Ｎρｖ２　・・・・（３）
【００４６】
また、圧力損失ΔＰはＰＭ堆積量ＭＬ　に対して一次関数形で変化することが分かった。さらに、この変化は、図５に示すように、ＰＭ堆積量ＭＬ　がある大きさ（以下、適宜、遷移点堆積量という）になる点（以下、適宜、遷移点という）で傾きが不連続に変化する。
【００４７】
これは、以下のように、パティキュレートフィルタ３２における圧力損失についての理論的なモデルからも導き出される結論であることを発明者らは確認した。図６はパティキュレートフィルタ３２内を表したもので、入口３２ａから流入した排気ガスが出口３２ｂから流出するまでの間に、管摩擦による圧力損失、通路断面積が急拡大したり急縮小することによる圧力損失、フィルタ本体４の隔壁を透過する際の圧力損失がある。以下に各圧力損失ΔＰｉ　について説明する。
【００４８】
（通路の摩擦による圧力損失）
この圧力損失ΔＰｉ　について式（４）と書ける。式中、λは損失の程度を示す係数であり、ｌは通路の長さであり、ｄは通路の径である。ｇは重力定数である。
ΔＰｉ　／（ρｇ）＝λ（ｌ／ｄ）（ｖ２　／２ｇ）・・・（４）
【００４９】
ここで、Ｒｅをレイノルズ数として、層流であればλ＝６４／Ｒｅであり、Ｒｅ＝ｖｄ／ｒ、ｒ＝μ／ρであるから、式（４）は式（５）と書ける。式中、αは係数である。
ΔＰｉ　＝αμｖ・・・（５）
【００５０】
（通路断面積の急縮小による圧力損失）
この圧力損失ΔＰｉ　について式（６）と書ける。式中、ζは損失の程度を示す係数である。
ΔＰｉ　／（ρｇ）＝ζ（ｖ２　／２ｇ）・・・（６）
【００５１】
したがって、式（６）は式（７）と書ける。式中、βは係数である。
ΔＰｉ　＝βρｖ２　・・・（７）
【００５２】
（通路断面積の急縮小による圧力損失）
この圧力損失ΔＰｉ　について式（８）と書ける。式中、ξは損失の程度を示す係数である。
ΔＰｉ　／（ρｇ）＝ξ（ｖ２　／２ｇ）・・・（８）
【００５３】
したがって、式（８）は式（９）と書ける。式中、γは係数である。
ΔＰｉ　＝γρｖ２　・・・（９）
【００５４】
（フィルタ本体４の隔壁を透過する際の圧力損失）
この圧力損失ΔＰｉ　について式（１０）と書ける。これはエルガンの式として知られている。式中、ｋ１　、ｋ２　は係数、εは気孔率、Ｓは多孔質体の表面積／体積、Ｌ＊　は透過層の厚さである。
ΔＰｉ　ｇ＝ｋ１　［（１−ε）２　／ε３　］Ｓ２　Ｌ＊　μｖ
＋ｋ２　［（１−ε）／ε３　］ＳＬ＊　ρｖ２　・・・（１０）
【００５５】
式（１０）は、フィルタ本体４の隔壁について（差圧ΔＰｉＷ）と、ＰＭ堆積層について（差圧ΔＰｉＳ）とのそれぞれについて成立するので、式（１１−１）、（１１−２）と表せる。式中、ＬはＰＭ堆積層の厚さである。
ΔＰｉＷ＝δＷ　μｖ＋ηＷ　ρｖ２　・・・（１１−１）
ΔＰｉＳ＝（δＳ　μｖ＋ηＳ　ρｖ２　）Ｌ・・・（１１−２）
【００５６】
差圧センサ５４で検出される圧力損失は、前記各圧力損失ΔＰｉ　、ΔＰｉＷ　、ΔＰｉＳを加算したものと考えられ、前記式（３）が成り立つ。ＰＭ堆積量が多ければＰＭ堆積層も厚くなるから、ＰＭ堆積量が多いほどＭ，Ｎは大きくなる。なお、圧力損失ΔＰｉ　、ΔＰｉＷ　、ΔＰｉＳの表式において、パティキュレートフィルタ３２よりも上流、下流での流速の項を無視しているのは、パティキュレートフィルタ３２内部での流速に比してごく小さいからである。
【００５７】
また、Ｌ＝ａＭＬ　とすれば式（１２−１）、式（１２−２）となる。
Ｍ＝α＋δＷ　＋δＳ　Ｌ＝α＋δＷ　＋δＳ　ａＭＬ　・・・・（１２−１）
Ｎ＝β＋δＳ　＋ηＳ　Ｌ＝β＋δＳ　＋ηＳ　ａＭＬ　・・・・（１２−２）
【００５８】
式（１２−１）、式（１２−２）を式（１３−１）、式（１３−２）と書き換えると、式（３）は式（１４）となる。Ａ，Ｂ，Ｃ，Ｄは定数である。
Ｍ＝Ａ＋ＢＭＬ　・・・・（１３−１）
Ｎ＝Ｃ＋ＤＭＬ　・・・・（１３−２）
ΔＰ＝（Ａμｖ＋Ｃρｖ２　）＋（Ｂμｖ＋Ｄρｖ２　）ＭＬ　・・・・（１４）
【００５９】
式（１４）を変形すると式（１５）となる。
ＭＬ　＝［ΔＰ−（Ａμｖ＋Ｃρｖ２　）］／（Ｂμｖ＋Ｄρｖ２　）
・・・・（１５）
【００６０】
ここで、図５の、遷移点堆積量のときを挟んで圧力損失ΔＰの変化率が異なることについて説明する。図７（Ａ）、図７（Ｂ）、図７（Ｃ）は、パティキュレートフィルタ本体４の隔壁で排気微粒子の堆積が進行していく様子を示しており、この順に、ＰＭ堆積量が多くなる。
【００６１】
図７（Ａ）は排気微粒子が堆積していない新品若しくはパティキュレートフィルタ３２を完全に再生した直後の状態であり、パティキュレートフィルタ本体４の隔壁を排気微粒子が透過する際における圧力損失は、パティキュレートフィルタ３２の形状諸元で規定される。
【００６２】
この状態から図７（Ｂ）に示すように、排気微粒子が、隔壁表面に堆積していく他、細孔に詰まっていき、圧力損失が増大するが、図中、矢印で示すように、排気ガスは細孔に向かうように流れていくので、細孔が詰まることが圧力損失の支配要因となる。
【００６３】
細孔の多くが詰まり全面にＰＭ堆積層が形成されると、今度は、図７（Ｃ）に示すように、ＰＭ堆積層の厚さが増大していくことになる。ここでは、隔壁全面を覆うＰＭ堆積層が厚くなることが圧力損失の支配要因となる。
【００６４】
このように、細孔の多くが詰まり、全面にＰＭ堆積層が形成される遷移点の前の領域（第１の領域）と後の領域（第２の領域）とで圧力損失の支配要因が異なる。排気微粒子が詰まっていない状態では良好に流通が自在であった細孔が、排気微粒子が細孔で捕集されて細孔に詰まると、急激に圧力損失が増大するので、細孔の多くが詰まるまでは、前記図５に示すように、ＰＭ堆積量に対する圧力損失の変化率は比較的、大きい（直線▲１▼）。一方、細孔の多くが詰まってしまった以降では支配要因がＰＭ堆積層の厚さの増大に変わるから、ＰＭ堆積量に対する圧力損失の変化率は緩やかなものに変わることになる（直線▲２▼）。
【００６５】
式（１４）は、ＰＭ堆積量に対して圧力損失が直線的に変化していく様子をよく記述するものであるが、遷移点を挟む支配要因の切り換わりを考慮すれば、Ａ，Ｂ，Ｃ，Ｄは２種類用意し、これが異なる式を使い分けることで、広い範囲でＰＭ堆積量を高精度に求めることができることになる。これが式（２−１）、式（２−２）である。式（２−１）のＡ１　，Ｂ１　，Ｃ１　，Ｄ１　は第１の領域でＡ，Ｂ，Ｃ，Ｄを適合したものであり、式（２−２）のＡ２　，Ｂ２　，Ｃ２　，Ｄ２　は第２の領域でＡ，Ｂ，Ｃ，Ｄを適合したものである。
【００６６】
本制御フローでは、この式（２−１）により算出した堆積量が基準値を越えるか否かを判断した後に、越えた場合には次ステップとして式（２−２）を用いて堆積量を算出する。しかし、本制御フローに限定されず、式（２−２）による判断に次いで式（２−１）による算出が行われてもよい。
【００６７】
また、基準値である、遷移点となる堆積量のデータは、予め実験的に求められて記憶されている。
【００６８】
なお、図８のように、酸化触媒が担持されていないフィルタ本体４０の上流にフロースルーモノリス担体４１が配置された触媒前置き構造のものでも、フロースルーモノリス担体４１において管摩擦があり、フロースルーモノリス担体４１の入口では急縮小が、出口では急拡大が生じることから、パティキュレートフィルタの構成によらず、圧力損失ΔＰは式（２−１）、式（２−２）、式（１４）と同様に表される。
【００６９】
ステップＳ１０４では、第１の領域に適合した式（２−１）によりＰＭ堆積量ＭＬ　を算出する。
【００７０】
ステップＳ１０５では、算出されたＰＭ堆積量ＭＬ　を遷移点のＰＭ堆積量と比較し、ＰＭ堆積量ＭＬ　が遷移点堆積量よりも小さいか否かを判定する。
【００７１】
肯定判断されると、ステップＳ１０６で、判定式　ＰＭ堆積量ＭＬ　≧基準堆積量ＭＬｔｈ　が成り立つか否かを判定する。基準堆積量ＭＬｔｈ　は、機関背圧や出力の低下がさほどでなくパティキュレートフィルタ３２の再生をしないことが許容される堆積量の上限値を考慮して設定する。
【００７２】
ステップＳ１０６が否定判断されると、ステップＳ１０１に戻り、ステップＳ１０１以降の手順が実行される。
【００７３】
そして、パティキュレートフィルタ３２において排気微粒子の堆積が進むと、ステップＳ１０５が否定判断されることになる。否定判断されるとステップＳ１０８に進み、遷移点以後の第２領域に適合した式（２−２）によりＰＭ堆積量ＭＬ　を算出する。これにより、より正確なＰＭ堆積量が知られる。
【００７４】
そして、ステップＳ１０６が肯定判断されるとステップＳ１０７に進み、パティキュレートフィルタ３２の再生を実施する。これは例えばポスト噴射等が用いられる。
【００７５】
このように、本排気ガス浄化装置では、ＰＭ堆積量を正確に知ることで適正な時期にパティキュレートフィルタ３２の再生を実施することができるので、再生時期が早すぎて燃費が悪化したり、再生時期が遅すぎて内燃機関の出力の低下やパティキュレートフィルタ３２での異常な昇温を招いたりすることを防止することができる。
【００７６】
しかも、式（２−１）、式（２−２）を記憶していればよいから、マップのようにデータが大規模化しない。
【００７７】
次に、ＥＣＵ５１で実行される別の制御フローを図９に示す。これは、パティキュレートフィルタ３２の再生の実施中に実行される制御フローであり、再生の開始とともに実行される。ステップＳ２０１〜Ｓ２０５は、堆積量低下状態検出手段としての処理である。
【００７８】
先ず、ステップＳ２０１でタイマのカウントをスタートさせるとともに、ステップＳ２０２で式（２−２）によりＰＭ堆積量を算出し、これを変化前値ＭＬｉとする。
【００７９】
ステップＳ２０３ではタイマのカウント値ｔを予め設定した基準時間ｔ０　と比較し、カウント値ｔが基準時間ｔ０　に達したか否かを判定し、肯定判断、すなわち基準時間が経過すると、ステップＳ２０４に進む。
【００８０】
ステップＳ２０４では、式（２−２）によりＰＭ堆積量を算出し、これを変化後値ＭＬｆとする。
【００８１】
ステップＳ２０５では、ＰＭ堆積量の差分ΔＭＬ　を式（１６）により算出する。
ΔＭＬ　＝ＭＬｉ−ＭＬｆ・・・（１６）
【００８２】
ステップＳ２０６は触媒劣化判定手段としての処理であり、差分ΔＭＬ　をしきい値ΔＭＬ０と比較し、差分ΔＭＬ　がしきい値ΔＭＬ０よりも小さいか否かを判定し、肯定判断されると、ステップＳ２０７で警告灯５５を点灯せしめ、本フローを終了する。ステップＳ２０６が否定判断されると、ステップＳ２０７をスキップして、本フローを終了する。
【００８３】
ＰＭ堆積量の差分ΔＭＬ　は、一定時間ｔ０　内のＰＭ堆積量の低減量であり、再生の進行速度を示している。これはパティキュレートフィルタ本体４に担持した触媒が劣化していると遅くなる。例えば経時劣化により、図１０のように走行距離が進むと漸次、遅くなる。
【００８４】
しきい値ΔＭＬ０を、触媒が正常とみなせる差分ΔＭＬ　の下限値（図中の触媒性能規定値に対応する）に設定することで、触媒の異常の発生や、交換すべき時期を適正に運転者に報知することができる。
【００８５】
なお、本制御フローはパティキュレートフィルタ３２の再生開始とともに実行しており、また、再生はＰＭ堆積量が上限堆積量ＭＬｔｈ　に達すると開始されるから、変化前値ＭＬｉは略上限堆積量ＭＬｔｈ　である。したがって、常に略同じ条件でＰＭ堆積量の低下状態を検出することができるが、ＰＭ堆積量が予め設定した所定値になった時点でカウンタをスタートするようにしてもよい。この場合、カウンタをスタートさせる時のＰＭ堆積量を、上限堆積量ＭＬｔｈ　よりもある程度小さな値に設定し、ＰＭ堆積量の燃焼が十分に進行している状態で計測するようにしてもよい。この場合、カウンタをスタートさせる時のＰＭ堆積量が第１領域であってもよい。遷移点よりも前か後かでＰＭ堆積量の算出に用いる式を式（２−１）と式（２−２）とで切り換えるのは勿論である。
【００８６】
また、基準時間ｔ０　後のＰＭ堆積量の変化をみるのではなく、再生中に、予め設定したＰＭ堆積量になった時点から、別の予め設定したＰＭ堆積量になる時点までの経過時間をカウントし、この経過時間が基準の時間よりも長くなったら酸化触媒の劣化と判断するのもよい。
【００８７】
なお、本実施形態では、ＰＭ堆積量を、圧力損失ΔＰ、体積流量Ｖ、ＤＰＦ入口温度およびＤＰＦ出口温度から算出した粘度μおよび密度ρに基づいて算出し、ＰＭ堆積量の大きさで再生が必要か否かを判定しているが、式（３）を利用するものであればよい。例えば、Ｍ，Ｎを再生が必要な基準のＰＭ堆積量に対応する値として、ΔＰ（ｖ，μ，ρ）＝Ｍμｖ＋Ｎρｖ２　を算出するようにし、判定式　ΔＰ≧ΔＰ（ｖ，μ，ρ）が成り立つか否かを判定する。再生が必要でなければ、ΔＰ＜ΔＰ（ｖ，μ，ρ）で判定式は不成立であり、再生が必要であればΔＰ≧ΔＰ（ｖ，μ，ρ）で判定式は成立となる。
【００８８】
本実施形態ではμ，ρを２つの温度センサにより検出したＤＰＦ入口温度およびＤＰＦ出口温度に基づいて算出しているが、いずれか一つの温度センサのみの温度に基づいて算出してもよい。この場合、検出温度に所定のオフセット温度を加減したものを、ＤＰＦ温度としてもよい。あるいは、要求される仕様によっては、粘度μおよび密度ρのいずれかひとつまたは両方を固定値とするのもよい。いずれか一つを固定値とする場合には、ＰＭ堆積量への依存性の小さいものの方を固定値とするのは勿論である。
【００８９】
また、式（１）では燃料噴射による体積流量の増分を考慮しているが、要求される仕様によっては式（１）の第１項を排気ガスの体積流量としてもよい。また、排気ガスの体積流量を、エアフローメータ５２により検出された吸入空気量を用いて算出するのではなく、排気通路３に流速を検出するセンサを設けて、直接に流速を検出するのもよい。
【図面の簡単な説明】
【図１】本発明の排気ガス浄化装置を適用した内燃機関の構成図である。
【図２】前記内燃機関の各部を制御するＥＣＵで実行される制御内容を示す第１のフローチャートである。
【図３】内燃機関が排出する排気ガスの粘度と温度との関係の一例を示すグラフである。
【図４】内燃機関が排出する排気ガスの密度と温度との関係の一例を示すグラフである。
【図５】前記パティキュレートフィルタにおける排気微粒子の堆積量と圧力損失との関係を示すグラフである。
【図６】前記排気ガス浄化装置のパティキュレートフィルタ内部における排気ガスの流通を示す図である。
【図７】（Ａ）、（Ｂ）、（Ｃ）は、それぞれ前記排気ガス浄化装置のパティキュレートフィルタに排気微粒子が堆積していく状態を示す堆積量の異なる図である。
【図８】別の構造を有するパティキュレートフィルタ内部における排気ガスの流通を示す図である。
【図９】前記内燃機関の各部を制御するＥＣＵで実行される制御内容を示す第２のフローチャートである。
【図１０】走行距離と、パティキュレートフィルタの再生時の排気微粒子浄化速度との関係を示すグラフである。
【符号の説明】
１　エンジン本体
２　吸気通路
２１　吸気マニホールド
３　排気通路
３１　排気マニホールド
３２　パティキュレートフィルタ
３２ａ　入口
３２ｂ　出口
４　本体
５１　ＥＣＵ（再生決定手段、補正手段、堆積量低下状態検出手段、触媒劣化判定手段）
５２　エアフローメータ（吸入空気量検出手段）
５３　温度センサ（温度検出手段）
５４　差圧センサ（差圧検出手段）
５５　警告灯
【図１】
2004019523.tif 000003
【図２】
2004019523.tif 000004
【図３】
2004019523.tif 000005
【図４】
2004019523.tif 000006
【図５】
2004019523.tif 000007
【図６】
2004019523.tif 000008
【図７】
2004019523.tif 000009
【図８】
2004019523.tif 000010
【図９】
2004019523.tif 000011
【図１０】
2004019523.tif 000012