"エンジンの排気ガスから粒子状物質を除去する粒子状物質除去フィルタの入口側と出口側との差圧および排気ガス流量をもとに前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定した差圧堆積量を算出する差圧堆積量算出部と、前記粒子状物質除去フィルタに入力される前記排気ガス内の粒子状物質量から前記粒子状物質除去フィルタで燃焼する粒子状物質量を減算するモデルを用いて前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定したモデル堆積量を算出するモデル堆積量算出部と、前記差圧堆積量の時間変化率を算出する差圧堆積量変化率算出部と、第1係数と第2係数との和が一定値となる条件で、前記差圧堆積量に第1係数を乗算した値と前記モデル堆積量に第2係数を乗算した値との加算値を粒子状物質堆積量として算出する際、前記排気ガス流量が所定値を超える場合、前記第1係数の値を前記第2係数の値に比して大きく設定し、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、該設定された前記第1係数を小さくして前記粒子状物質堆積量を算出する粒子状物質堆積量算出部と、を備えたことを特徴とする粒子状物質堆積量推定装置。前記粒子状物質堆積量算出部は、前記第1係数と前記第2係数との和を1に設定し、前記排気ガス流量が所定値を超える場合、前記第1係数の値を1に設定し、前記排気ガス流量が所定値を超えない場合、前記第1係数の値を0に設定することを特徴とする請求項1に記載の粒子状物質堆積量推定装置。前記粒子状物質堆積量算出部は、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、前記差圧堆積量の時間変化率の増加に応じて、設定された前記第1係数を減少させることを特徴とする請求項1または2に記載の粒子状物質堆積量推定装置。前記粒子状物質堆積量算出部は、算出した前記粒子状物質堆積量が閾値を超える場合、前記粒子状物質除去フィルタに対する手動強制再生を行わせる指示を出力することを特徴とする請求項1~3のいずれか一つに記載の粒子状物質堆積量推定装置。エンジンの排気ガスから粒子状物質を除去する粒子状物質除去フィルタと、前記エンジンを制御するエンジン制御部と、前記粒子状物質除去フィルタの再生を制御する再生制御部と、前記粒子状物質除去フィルタの入口側と出口側との差圧を検出する差圧センサと、前記粒子状物質除去フィルタに対する手動強制再生を指示する再生指示部と、を備え、前記再生制御部は、前記差圧および前記エンジン制御部から入力される排気ガス流量をもとに前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定した差圧堆積量を算出する差圧堆積量算出部と、前記粒子状物質除去フィルタに入力される前記排気ガス内の粒子状物質量から前記粒子状物質除去フィルタで燃焼する粒子状物質量を減算するモデルを用いて前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定したモデル堆積量を算出するモデル堆積量算出部と、前記差圧堆積量の時間変化率を算出する差圧変化率算出部と、第1係数と第2係数との和が一定値となる条件で、前記差圧堆積量に第1係数を乗算した値と前記モデル堆積量に第2係数を乗算した値との加算値を粒子状物質堆積量として算出する際、前記排気ガス流量が所定値を超える場合、前記第1係数の値を前記第2係数の値に比して大きく設定し、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、該設定された前記第1係数を小さくして前記粒子状物質堆積量を算出し、この算出した前記粒子状物質堆積量が閾値を超える場合、前記手動強制再生を行うべき旨の指示を前記再生指示部に出力するする粒子状物質堆積量算出部と、を備えたことを特徴とする排気ガス浄化システム。前記再生制御部は、前記再生指示部から手動強制再生実行の指示を受け付けた場合、前記エンジン制御部を介して手動強制再生を実行することを特徴とする請求項5に記載の排気ガス浄化システム。エンジンの排気ガスから粒子状物質を除去する粒子状物質除去フィルタの入口側と出口側との差圧および排気ガス流量をもとに前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定した差圧堆積量を算出する差圧堆積量ステップと、前記粒子状物質除去フィルタに入力される前記排気ガス内の粒子状物質量から前記粒子状物質除去フィルタで燃焼する粒子状物質量を減算するモデルを用いて前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定したモデル堆積量を算出するモデル堆積量算出ステップと、前記差圧堆積量の時間変化率を算出する差圧変化率算出ステップと、第1係数と第2係数との和が一定値となる条件で、前記差圧堆積量に第1係数を乗算した値と前記モデル堆積量に第2係数を乗算した値との加算値を粒子状物質堆積量として算出する際、前記排気ガス流量が所定値を超える場合、前記第1係数の値を前記第2係数の値に比して大きく設定し、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、該設定された前記第1係数を小さくして前記粒子状物質堆積量を算出する粒子状物質堆積量算出ステップと、を含むことを特徴とする粒子状物質堆積量推定方法。ところで、図13上段に示すように、DPFを構成するセル内に捕集されたPM101はセル内で一様に堆積されるが、PM101が多くなると、図13下段に示すように、セル内で堆積されたPM101が所々剥がれ、剥がれたPM102によってセル閉塞が生じる場合があることが確認されている。ここで、上述した特許文献1のように、差圧を用いてPM堆積量を推定する場合では、セル閉塞が生じたときに問題が生じる。すなわち、セル閉塞が生じたときに、差圧を用いてPM堆積量を推定する場合では、PMの推定量が実際のPM堆積量より大きくなってしまうため、実際には、手動強制再生が必要なほどにはPMが堆積していないにもかかわらず、手動強制再生を促す警告を発報してしまうという問題点があった。この発明は、上記に鑑みてなされたものであって、PM堆積量を精度高く推定することができる粒子状物質堆積量推定装置、これを用いた排気ガス浄化システム、および粒子状物質堆積量推定方法を提供することを目的とする。上述した課題を解決し、目的を達成するために、この発明にかかる粒子状物質堆積量推定装置は、エンジンの排気ガスから粒子状物質を除去する粒子状物質除去フィルタの入口側と出口側との差圧および排気ガス流量をもとに前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定した差圧堆積量を算出する差圧堆積量算出部と、前記粒子状物質除去フィルタに入力される前記排気ガス内の粒子状物質量から前記粒子状物質除去フィルタで燃焼する粒子状物質量を減算するモデルを用いて前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定したモデル堆積量を算出するモデル堆積量算出部と、前記差圧堆積量の時間変化率を算出する差圧変化率算出部と、第1係数と第2係数との和が一定値となる条件で、前記差圧堆積量に第1係数を乗算した値と前記モデル堆積量に第2係数を乗算した値との加算値を粒子状物質堆積量として算出する際、前記排気ガス流量が所定値を超える場合、前記第1係数の値を前記第2係数の値に比して大きく設定し、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、該設定された前記第1係数を小さくして前記粒子状物質堆積量を算出する粒子状物質堆積量算出部と、を備えたことを特徴とする。また、この発明にかかる粒子状物質堆積量推定装置は、上記の発明において、前記粒子状物質堆積量算出部は、前記第1係数と前記第2係数との和を1に設定し、前記排気ガス流量が所定値を超える場合、前記第1係数の値を1に設定し、前記排気ガス流量が所定値を超えない場合、前記第1係数の値を0に設定することを特徴とする。また、この発明にかかる粒子状物質堆積量推定装置は、上記の発明において、前記粒子状物質堆積量算出部は、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、前記差圧堆積量の時間変化率の増加に応じて、設定された前記第1係数を減少させることを特徴とする。また、この発明にかかる粒子状物質堆積量推定装置は、上記の発明において、前記粒子状物質堆積量算出部は、算出した前記粒子状物質堆積量が閾値を超える場合、前記粒子状物質除去フィルタに対する手動強制再生を行わせる指示を出力することを特徴とする。また、この発明にかかる排気ガス浄化システムは、エンジンの排気ガスから粒子状物質を除去する粒子状物質除去フィルタと、前記エンジンを制御するエンジン制御部と、前記粒子状物質除去フィルタの再生を制御する再生制御部と、前記粒子状物質除去フィルタの入口側と出口側との差圧を検出する差圧センサと、前記粒子状物質除去フィルタに対する手動強制再生を指示する再生指示部と、を備え、前記再生制御部は、前記差圧および前記エンジン制御部から入力される排気ガス流量をもとに前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定した差圧堆積量を算出する差圧堆積量算出部と、前記粒子状物質除去フィルタに入力される前記排気ガス内の粒子状物質量から前記粒子状物質除去フィルタで燃焼する粒子状物質量を減算するモデルを用いて前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定したモデル堆積量を算出するモデル堆積量算出部と、前記差圧堆積量の時間変化率を算出する差圧変化率算出部と、第1係数と第2係数との和が一定値となる条件で、前記差圧堆積量に第1係数を乗算した値と前記モデル堆積量に第2係数を乗算した値との加算値を粒子状物質堆積量として算出する際、前記排気ガス流量が所定値を超える場合、前記第1係数の値を前記第2係数の値に比して大きく設定し、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、該設定された前記第1係数を小さくして前記粒子状物質堆積量を算出し、この算出した前記粒子状物質堆積量が閾値を超える場合、前記手動強制再生を行うべき旨の指示を前記再生指示部に出力するする粒子状物質堆積量算出部と、を備えたことを特徴とする。また、この発明にかかる排気ガス浄化システムは、上記の発明において、前記再生制御部は、前記再生指示部から手動強制再生実行の指示を受け付けた場合、前記エンジン制御部を介して手動強制再生を実行することを特徴とする。また、この発明にかかる粒子状物質堆積量推定方法は、エンジンの排気ガスから粒子状物質を除去する粒子状物質除去フィルタの入口側と出口側との差圧および排気ガス流量をもとに前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定した差圧堆積量を算出する差圧堆積量ステップと、前記粒子状物質除去フィルタに入力される前記排気ガス内の粒子状物質量から前記粒子状物質除去フィルタで燃焼する粒子状物質量を減算するモデルを用いて前記粒子状物質除去フィルタに堆積する粒子状物質堆積量を推定したモデル堆積量を算出するモデル堆積量算出ステップと、前記差圧堆積量の時間変化率を算出する差圧変化率算出ステップと、第1係数と第2係数との和が一定値となる条件で、前記差圧堆積量に第1係数を乗算した値と前記モデル堆積量に第2係数を乗算した値との加算値を粒子状物質堆積量として算出する際、前記排気ガス流量が所定値を超える場合、前記第1係数の値を前記第2係数の値に比して大きく設定し、前記排気ガス流量が前記所定値を超え、かつ前記差圧堆積量の時間変化率が所定閾値以上である場合、該設定された前記第1係数を小さくして前記粒子状物質堆積量を算出する粒子状物質堆積量算出ステップと、を含むことを特徴とする。この発明によれば、第1係数と第2係数との和が一定値となる条件で、差圧堆積量に第1係数を乗算した値とモデル堆積量に第2係数を乗算した値との加算値を粒子状物質堆積量として算出する際、排気ガス流量が所定値を超える場合、前記第1係数の値を前記第2係数の値に比して大きく設定し、前記排気ガス流量が前記所定値を超え、かつ差圧堆積量の時間変化率が所定閾値以上である場合、該設定された前記第1係数を小さくして前記粒子状物質堆積量を算出するようにしているので、精度の高い粒子状物質堆積量を推定することができる。これにより、DPF内のセル閉塞が生じた場合のように、実際の粒子状物質堆積量とは無関係に差圧が急上昇し、推定した粒子状物質堆積量が非常に大きくなることによって、手動強制再生を促す誤った警報を出してしまうことを防止することができる。以下、添付図面を参照してこの発明を実施するための形態について説明する。図1は、この発明の実施の形態にかかる排気ガス浄化システムを含むディーゼルエンジン100の概要構成を示す模式図である。図1において、ディーゼルエンジン100は、ディーゼルエンジン100は、内部に複数の燃焼室が形成されたエンジン本体1と、吸入される空気をフィルタを用いて濾過し、埃などの異物が燃焼室に混入するのを防止するエアクリーナ2と、エンジン本体1内部の各燃焼室に給気を供給する給気管路3と、エンジン本体1内部の各燃焼室から排出された排気ガスを排出する排気管路4と、冷却機構5と、排気タービン過給機6と、排気浄化装置7と、排気再循環システム8と、エンジンコントローラ30と、再生制御部40と、再生指示部50とを有する。を備える。エンジン本体1と給気管路3との間には、給気管路3からの給気がエンジン本体1内部の各燃焼室に分配されるように、給気マニホールド3Aが取り付けられている。エンジン本体1と排気管路4との間には、エンジン本体1内部の各燃焼室から排出された排気ガスがまとめて排気管路4に流入するように、排気マニホールド4Aが取り付けられている。給気管路3には、排気タービン過給機6によって圧縮された空気を冷却するためのアフタークーラ11が設けられている。冷却機構5は、エンジン本体2内に収められた図示しないクランクシャフト等により駆動されるポンプ12を備える。ポンプ12によって圧送された冷却水は、エンジン本体1,排気タービン過給機6,図示しないオイルクーラ等の冷却が必要な部位を冷却した後、冷却機構5に設けられたラジエータ13によって空冷されるようになっている。アフタークーラ11とラジエータ13とは、エンジン本体1に設けられ、且つ、図示しないクランクシャフト等によって回転駆動されるファン14によって、その冷却作用が促進されるようになっている。排気タービン過給機6は、排気管路4の途中に設けられたタービン21と、給気管路3の途中に設けられ、タービン21に連結されて駆動されるコンプレッサ22と、タービン21に供給される排気ガスの流速を制御するために可変ターボノズル23とを備える。排気タービン過給機6は、可変ターボノズル23の開度を制御することにより、タービン21の回転数を制御する。タービン21の回転によってコンプレッサ22が駆動し、エンジン本体1への給気過給が行われる。なお、可変ターボノズル23は、全閉時には、パイパス路24を介して排気浄化装置7側に排気するようにしている。すなわち、可変ターボノズル23が開時には、排気ガスをタービン翼車21aに供給して仕事をさせ、可変ターボノズル23が全閉時には、排気ガスを、パイパス路24を介して排気浄化装置7側に出力してタービン翼車21aへの仕事を小さくして排気温度を高めるようにしている。排気浄化装置2は、タービン21の下流側に設けられ、排気ガスに含まれるPMを除去するものであり、DOCディーゼル酸化触媒71、DPF72、差圧センサ73、温度センサ74を有する。DOC71およびDPF72は、円筒状の排気管路内部に設けられ、排気管路の上流側にDOC71が設けられ、排気管路の下流側にDPF72が設けられる。また、タービン21と排気浄化装置7との間には、ドージング燃料供給装置70から供給されるドージング燃料を噴射するドージングノズル70aが配置される。このドージング燃料の噴射は、強制再生が指示された場合に行われる。DOC71は、Ptなどによって実現され、排気ガスに含まれるCO一酸化炭素、HC炭化水素、PMに含まれるSOF有機可溶成分を酸化して除去する。さらに、DOC21は、排気ガスに含まれるNO一酸化窒素を酸化してNO2二酸化窒素に変化させ、さらにはドージングノズル70aから噴射されたドージング燃料を酸化することによって排気ガス温度を上昇させる。DPF72は、PMを捕集する。DPF72は、炭化珪素等を基材として実現される。排気ガスに含まれるPMは、DPF72に形成された微細な穴を通過する際に捕集される。そして、DPF72は、図14に示すように、排気ガスの流れ方向に沿った微細流路を有するセルが円筒状の排気管路内に密集配置される。そして、上流側端部を目封じさせたセルと、下流側端部を目封じさせたセルとを交互に配置したウォールフロー型DPFである。捕集されたPMは、排気ガスが酸化反応を進行させることができる温度であることを条件として、排気ガスに含まれる酸素およびDOC71で生成されたNO2によって酸化燃焼されることになる。差圧センサ73は、DPF72の上流側に配置されてDPF72の上流側の圧力を検出する圧力センサ73aと、DPF72の下流側に配置されてDPF72の下流側の圧力を検出する圧力センサ73bと、圧力センサ73aが検出した圧力から圧力センサ73bが検出した圧力を減算した差圧を再生制御部4に出力する差圧検出部73cとを有する。温度センサ74は、DPF72の上流側に配置され、DPF72の入口の排気温度を検出し、DPF温度として再生制御部4に出力する。排気再循環システム8は、排気マニホールド4Aと給気管路3とを連通する排気再循環通路31を備える。排気再循環通路31は、排気マニホールド4Aから排気ガスの一部を抽出して給気管路3に再循環させる。排気再循環通路31には、排気再循環通路31を開閉するEGRバルブ32と、排気マニホールド4Aからの排気ガスを冷却するEGRクーラ33とが設けられている。排気再循環システム8は、排気再循環通路31を介して排気ガスの一部を給気マニホールド3Aに還流させることによって、給気中の酸素濃度を低下させ、エンジン本体1の燃焼温度を下げる。これにより、排気ガス中に含まれる窒素酸化物の量を低減させることができる。ここで、ディーゼルエンジン100は、制御系として、エンジン回転速度センサ3aと、給気圧センサ3cと、排気圧センサ3dと、タービン回転速度センサ3e,流量センサ3fを備える。エンジン回転速度センサ3aは、エンジン本体1の図示しないクランクシャフトの回転速度を検出し、図示しないクランクシャフトの回転速度を示す信号をエンジンコントローラ30に入力する。給気圧センサ4cは、コンプレッサ22の出口通路と給気マニホールド3Aとの間の給気圧力を検出し、エンジンコントローラ30に入力する。排気圧センサ4dは、排気マニホールド4Aとタービン21の入口通路との間の排気圧力を検出し、排気圧力を示す信号をエンジンコントローラ30に入力する。タービン回転速度センサ4eは、タービン21の回転速度を検出し、タービン21の回転速度を示す信号をエンジンコントローラ30に入力する。エンジンコントローラ30は、図示しないアクセルペダルなどのオペレータの要求に応じた入力手段の入力値に応じて、燃料噴射量、燃料噴射タイミング、EGRバルブ32、可変ターボノズル23を調整してエンジン回転数やトルク制御を行うとともに、再生制御部40からの指示によって、燃料噴射量、燃料噴射タイミング、EGRバルブ32、可変ターボノズル23を調整して排気温度を上昇させ、その後ドージングノズル70aからドージング燃料を噴射することによって、強制再生制御を行う。特に、エンジンコントローラ30は、強制再生制御を行う場合、燃料噴射量を抑え、EGRバルブ32、可変ターボノズル23をそれぞれ閉にすることによって排気温度を上昇させる。再生制御部40は、差圧堆積量算出部41、モデル堆積量算出部42、差圧堆積量時間変化率算出部43、およびPM堆積量算出部44を有する。再生制御部40は、エンジンコントローラ30から得られる排気ガス流量情報、PM発生量情報、PM燃焼量情報、差圧センサ73から得られる差圧、および温度センサ74から得られるDPF温度をもとにDPF72に堆積するPM堆積量を推定し、PM堆積量が所定の閾値を越える場合に、再生指示部50に通知する。また、再生制御部40は、再生指示部50からの手動再生指示があった場合に、エンジンコントローラ30に手動強制再生を行わせる。また、手動強制再生が終了した場合、その旨を再生指示部50に通知する。再生指示部50は、運転席近傍に設けられた入出力操作が可能な表示パネルに設けられ、再生制御部40から通知があった場合、手動再生指示を促す警告表示を行うとともに、手動再生指示の操作があった場合に、手動再生指示を再生制御部40に出力する。また、再生指示部50は、手動強制再生指示の終了通知があった場合、その旨を表示出力する。ここで、再生制御部40による再生制御処理について説明する。図2は、再生制御部40の詳細構成を示すブロック図である。図2において、差圧堆積量算出部41には、差圧センサ73からのDPF72の差圧と、温度センサ74からのDPF温度と、エンジンコントローラ30からの排気ガス流量情報とが入力される。排気ガス流量情報は、DPF72に入力される排気ガス体積流量である。エンジンコントローラ30は、流量センサ3fが検出した給気質量流量と燃料噴射量センサ3bが検出した燃料噴射量から推定される排気質量流量とを加算した質量流量を差圧堆積量算出部41に出力する。そして、差圧堆積量算出部41は、エンジンコントローラ30から入力された排気ガス流量を、温度センサ74から入力されたDPF温度をもとに排気体積流量に変換する。差圧堆積量算出部41は、この変換した排気ガス流量と差圧とPM堆積量との関係を示す3次元マップを有し、変換された排気ガス流量と差圧とを入力値としてPM堆積量を出力する。図3は、PM堆積量をパラメータとした排気ガス流量と差圧との関係を示す図である。排気ガス流量が一定の場合、差圧の増大に伴ってPM堆積量は増大する。また、差圧が一定の場合、排気ガス流量の増大に伴ってPM堆積量は減少する。たとえば、図3に示すように、排気ガス流量がQV1で差圧がDP1~DP4の場合、それぞれPM堆積量は、P1~P4となる。すなわち、PM堆積量は、P1よりP4のほうが大きい。そして、差圧堆積量算出部41は、推定したPM堆積量を差圧堆積量PMaとしてPM堆積量算出部44および差圧堆積量時間変化率算出部43に出力する。一方、図4に示すように、モデル堆積量算出部42には、エンジンコントローラ30からのPM発生量情報およびPM燃焼量情報と、温度センサ74からのDPF温度とが入力される。エンジンコントローラ30は、燃料噴射量や空燃比などをもとに推定した排気ガス内のPM発生量を示すPM発生量情報と、DOC21から出力されると推定したNO2発生量およびDOC71から出力されるO2発生量を示すPM燃焼量情報とをモデル堆積量算出部42に出力する。モデル堆積量算出部42は、エンジンから発生するPM発生量からDPF72で燃焼するPM燃焼量を減算することによってPM堆積量を推定する燃焼モデルMを用いている。なお、燃焼モデルMは、PM堆積速度を算出している。NO2によるPM燃焼量算出部42aは、エンジンコントローラ30から入力されたNO2発生量とDPF温度とをもとに、あらかじめ実測等で求めておいたDPF72に堆積したPMがNO2によって酸化燃焼する反応速度をマップ等で記憶しておいたものを用いて、NO2によるPM燃焼速度を算出する。O2によるPM燃焼量算出部42bは、エンジンコントローラ30から入力されたO2発生量とDPF温度とをもとに、あらかじめ実測等で求めておいたDPF72に堆積したPMがO2によって酸化燃焼する反応速度をマップ等で記憶しておいたものを用いて、O2によるPM燃焼速度を算出する。ここで、エンジンコントローラ3から入力されるPM発生量は、PM発生速度であり、演算器42cは、PM発生速度からNO2によるPM燃焼速度およびO2によるPM燃焼速度を減算したPM堆積速度を積算部42dに出力する。積算部42dは、入力されたPM堆積速度をもとにPM堆積量を積算し、この積算したPM堆積量をモデル堆積量PMbとしてPM積算量算出部44に出力する。図2に戻り、PM堆積量算出部44は、差圧堆積量算出部41から入力された差圧堆積量PMaとモデル堆積量算出部42から入力されたモデル堆積量PMbとに対して重み係数をかけてPM堆積量PMsを算出する。すなわち、PM堆積量PMsは、係数α0≦α≦1として次式1で表せる。PMs=αPMa1-αPMb...1すなわち、差圧堆積量PMaにかかる重み係数である係数αと、モデル堆積量PMbにかかる重み係数である係数1-αとは、それぞれ加算して1となるようにしている。換言すれば、係数αによって、差圧堆積量PMaとモデル堆積量PMbとの配分比を決定している。係数αが1であれば、PM堆積量Psは、差圧堆積量PMaの値となり、係数αが0であれば、PM堆積量Psは、モデル堆積量PMbの値となる。また、PM堆積量算出部44は、係数αを決定する係数決定部45を有する。係数決定部45には、エンジンコントローラ30からの排気ガス流量情報と、温度センサ24からのDPF温度とが入力され、DPF温度で変換補正したDPF72での排気ガス流量を算出し、この排気ガス流量の値が所定値QVthを越えた場合には、係数αを1に設定し、排気ガス流量の値が所定値QVth未満である場合には、係数αを0に設定する。一方、差圧堆積量時間変化率算出部43は、差圧堆積量算出部41から入力される差圧堆積量の時間変化率である差圧堆積量時間変化率ΔPMaを算出し、係数決定部45内の係数修正部46に出力する。係数決定部46は、差圧堆積量時間変化率ΔPMaが閾値ΔPMathを越えた場合に、係数αの値を、係数αの値よりも小さい係数α’に修正する。図5および図6に示すフローチャートを参照して、PM堆積量算出部44の処理手順について説明する。まず、図5において、PM堆積量算出部44は、まず、差圧堆積量算出部41から差圧堆積量PMaを取得するステップS101。さらに、PM堆積量算出部44は、モデル堆積量算出部42からモデル堆積量PMbを取得するステップS102。その後、係数決定部45によって係数αを決定する係数決定処理を行わせるステップS103。その後、PM堆積量算出部44は、決定された係数αを用い、式1に従って、PM堆積量PMsを算出するステップS104。その後、PM堆積量算出部44は、算出されたPM堆積量PMsが閾値PMthを越えているか否かを判断するステップS105。PM堆積量PMsが閾値PMthを越えていない場合ステップS105,Noには、ステップS101に移行して上述した処理を繰り返す。一方、PM堆積量PMsが閾値PMthを越えている場合ステップS105,Yesには、再生指示部5に対して手動強制再生をすべき旨の指示を行いステップS106、その後、ステップS101に移行して上述した処理を繰り返す。ステップS103に示す係数決定処理は、図6に示したフローチャートに従って行われる。すなわち、係数決定部45は、まず、排気ガス流量が所定値QVthを越えているか否かを判断するステップS201。排気ガス流量が所定値QVthを越えていない場合ステップS201,Noには、係数αを0に設定してステップS203、ステップS103にリターンする。一方、排気ガス流量が所定値QVthを越えている場合ステップS201,Yesには、さらに差圧堆積量時間変化率ΔPMaが閾値ΔPMathを越えているか否かを判断するステップS202。差圧堆積量時間変化率ΔPMaが閾値ΔPMathを越えていない場合ステップS202,Noには、係数αを1に設定してステップS204、ステップS103にリターンする。一方、差圧堆積量時間変化率ΔPMaが閾値ΔPMathを越えている場合ステップS202,Yesには、係数αを1よりも小さい係数α’に修正されステップS205、ステップS103にリターンする。この結果、係数αは、図7に示すように決定される。すなわち、排気ガス流量QVが所定値QVth未満である場合、係数αは0となる。所定値QVthを越えている場合であって差圧堆積量時間変化率ΔPMaが閾値ΔPMath未満である場合、係数αは1となる。所定値QVthを越えている場合であって差圧堆積量時間変化率ΔPMaが閾値ΔPMathを越えている場合、係数αは1よりも小さい係数α’に修正される。ここで、排気ガス流量QVが所定値QVth未満図8に示す領域R1の場合に係数αを0とするのは、領域R1では、差圧堆積量PMaの推定精度が低くなるからである。このため、領域R1では、係数αを0としてPM堆積量PMsをモデル堆積量PMbで推定するようにしている。また、排気ガス流量QVが所定値QVthを越えた図8に示す領域R2の場合には係数αを1としてPM堆積量PMsを差圧堆積量PMaで推定するようにしている。一方、排気ガス流量QVが所定値QVthを越えた図8に示す領域R2の場合に係数αを1としてPM堆積量PMsをすべて差圧堆積量PMaで推定すると、図13に示したようにDPF72のセル内でPMが剥がれて、セルが閉塞される場合があり、この場合、実際のPM堆積量とは無関係に差圧が急上昇してPM堆積量の推定精度が極端に悪くなるとともに、無用な手動強制再生を指示することになる。このため、差圧堆積量時間変化率ΔPMaが閾値ΔPMathを越えた場合、係数αを1よりも小さい係数α’にして、セル閉塞による差圧急上昇に伴う差圧堆積量PMaの増大分の影響を小さくなるようにしている。たとえば、図9に示す時点t1でDPF72のセル閉塞が発生した場合、図9bに示すように、セル閉塞に伴って、破線部のように差圧DPは急上昇する。そのとき、図9cに示すように、時点t1で差圧堆積量時間変化率ΔPMaは、破線部のように閾値ΔPMathを越える。その場合、図9aに示すように、算出されたPM堆積量PMsは、手動強制再生指示を行うための閾値PMthを越えてしまう。この場合、図9aの実線で示す実際のPM堆積量と破線で示す算出されたPM堆積量との差である斜線分のPM堆積量増加分ΔSは、大きな誤差となる。このため、係数αを1よりも小さい係数α’にして、PM堆積量増加分ΔSを減少させ、モデル堆積量PMbの推定割合を増やして補うようにしている。したがって、この実施の形態では、排気ガス流量QVが所定値QVthを越えない領域R1では、差圧堆積量PMaに比して誤差の小さいモデル堆積量PMbを用いてPM堆積量PMsを推定し、排気ガス流量QVが所定値QVthを越える領域R2では、精度の高い差圧堆積量PMaを用いてPM堆積量PMsを推定するようにし、さらに、排気ガス流量QVが所定値QVthを越える領域R2では、係数αを小さくしてDPF22のセル閉塞に伴う差圧急上昇による誤差成分が小さくなるようにしているので、精度の高く無駄のない手動強制再生の指示を確実に行うことができる。なお、上述した実施の形態では、差圧堆積量時間変化率ΔPMaの変動にかかわらず、係数α’を予め設定された一定の値としていた図10のL1が、これに限らず、係数α’を差圧堆積量時間変化率ΔPMaの関数値としてもよい。すなわち、α‘=fΔPMaとし、図10のL2に示すように、差圧堆積量時間変化率ΔPMaの増加に伴って係数α‘の値が減少するようにすればよい。これによって、可能な限り差圧堆積値PMaを用いて推定精度を上げるとともに、無駄な手動強制再生の指示発生を抑えることができる。なお、差圧堆積量時間変化率ΔPMaに替えて、差圧時間変化率の関数としてもよい。また、上述した実施の形態では、図7に示すように、排気ガス流量QVの所定値QVthで急激に係数が変化するようにしていたが、これに限らず、図11に示すように、所定値QVthを越えてから、滑らかに係数αが変化するようにモジュレーションをかけるようにしてもよい。さらに、上述した実施の形態では、差圧堆積量時間変化率ΔPMaの急上昇の発生がない場合、領域R1で係数αを0に設定し、領域R2で係数αを1に設定することを前提としていたが、これに限らず、図12に示すように、領域R1の係数αは0でなく、領域R1では主にモデル堆積量PMbが用いられ、また、領域R2の係数αは1でなく、領域R2では主に差圧堆積量PMaが用いられるようにしてもよい。"