JP2008544155A - 粒子フィルタの再生制御方法 - Google Patents

Abstract

内燃エンジンの排気管に設けられる粒子フィルタの再生を制御する方法であって、本方法は、粒子フィルタの動作データアイテムを求めるステップであって、前記データアイテムが、時刻ｔでのフィルタの少なくとも一つの代表内部温度Ｔ（ｔ）である構成のステップと、フィルタの温度を、フィルタに堆積した煤の第１燃焼閾値温度よりも高い温度に上昇させることにより、フィルタを再生するステップと、時間と共に変化し、かつ時刻ｔにおいて再生プロセスによって生成されるエネルギーの量を表わすパラメータｄ_ｉｎｓｔ（ｔ）を計算するステップとを含む。本発明による再生制御方法は更に、関数ｄ_ｉｎｓｔ（ｔ）^＊Ｇｔｅｍｐ（ｔ）又はＧｔｅｍｐ（ｔ）によって少なくとも部分的に定義される重み付けパラメータＤ_ｉｎｓｔ ｗｅｉｇｈｔ（ｔ）を計算するステップを含み、Ｇｔｅｍｐ（ｔ）は前記代表温度Ｔ（ｔ）に従い時刻（ｔ）とともに変化する。

Description

技術分野
本発明は概して、粒子フィルタの再生を制御する方法に関する。
更に詳細には、本発明は、内燃エンジンの排気管に設けられる粒子フィルタの再生を制御する方法に関連するものであり、本方法は、
− 粒子フィルタの動作データを求めるステップであって、これらのデータが、時刻ｔでのフィルタ内の温度を表わす少なくとも一つの温度Ｔ（ｔ）である構成のステップと、
− フィルタを再生するステップであって、フィルタの再生において、前記フィルタの温度を、前記フィルタに堆積した煤が燃焼する第１閾値温度よりも高い温度に上昇させるステップと、
− 時間と共に変化することができ、かつ時刻ｔにおいて再生プロセスに伴うエネルギーの量を表わすパラメータｄ_ｉｎｓｔ（ｔ）を計算するステップとを含む。

背景技術
内燃エンジンによって外気に放出される粒子の量を減らすために、粒子フィルタを排気管に設けて燃焼ガスをエンジンから排出する方法が知られている。この粒子フィルタは、フィルタの濾過機能を良好に維持するために、再生する必要があり、このような再生では、フィルタ内の温度を上昇させて、フィルタに堆積した煤を燃焼させる。

粒子フィルタの再生によって、粒子フィルタにダメージが加わり、そしてフィルタ寿命が短くなる恐れがある。
このような現象が発生するので、粒子フィルタを装着した排気管を組み込んだ内燃エンジンの多くの製造業者は、フィルタ再生プロセスを制御する種々の解決方法を開発しているのである。

このような再生制御を可能にする、上述したタイプの再生制御方法が、例えば特許文献仏国特許発明第２８５５２１３号に記載されている。
この文献はフィルタの再生を制御する方法を提示しており、フィルタ再生の間、時間と共に変化することができ、かつ再生に伴うエネルギーの量を表わすパラメータｄ_ｉｎｓｔ（ｔ）が計算される。このパラメータｄ_ｉｎｓｔ（ｔ）は比で求まり、この比は、１秒当たりに燃焼する煤の質量としての煤の燃焼速度Ｖ_ｒ（ｔ）を、排気ガス質量としての排気ガス流量Ｑ_ｍ（ｔ）で割った計算値である。このような比で表わされるので、燃焼速度Ｖ_ｒ（ｔ）が大きくなると、煤の燃焼によって所定時刻ｔに生成されるエネルギーの量が大きくなる。同様に、排気ガス流量Ｑ_ｍ（ｔ）が大きくなると、煤の燃焼エネルギーをより迅速に除去することができる。

従って、この計算パラメータｄ_ｉｎｓｔ（ｔ）の変化を観察して、再生プロセスを中断してフィルタを劣化させる危険を回避する必要があるかどうかを判断する。
このような背景から、本発明の目的は再生を制御する方法を提案することにあり、この方法によって、フィルタの再生処理中に粒子フィルタにダメージを与える危険を小さくすることができる。

課題を解決するための手段
この目的を達成するために、上の序文に示される一般的な定義も満たす本発明の再生制御方法は基本的に、前記方法が、関数ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）によって少なくとも部分的に定義される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）を計算するステップを含み、Ｇｔｅｍｐ（ｔ）は時刻（ｔ）とともに変化することができ、かつ前記代表温度Ｔ（ｔ）の関数として表わされることを特徴とする。

重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）を、パラメータｄ_ｉｎｓｔ（ｔ）に、代表温度Ｔ（ｔ）とともに変化することができる関数である関数Ｇｔｅｍｐ（ｔ）を乗算して計算することにより、再生の変化をより正確に見積もって、例えばフィルタにダメージを与える全ての危険を防止することができる。これは、ｄ_ｉｎｓｔ（ｔ）を温度Ｔ（ｔ）の関数としての重み付けする処理では、フィルタからガス流によって実際に除去することができる熱量をより正確に取り込むことができるからである。

例えば構成は、パラメータｄ_ｉｎｓｔ（ｔ）を次の公式によって得られるように決定することができる。

上の式では、Ｖｒ（ｔ）は時刻ｔにおいて粒子フィルタに堆積した煤の燃焼速度であり、そしてＱ_ｍ（ｔ）は時刻ｔにおいて推定される、または測定される排気ガス質量流量である。このパラメータｄ_ｉｎｓｔ（ｔ）は特に、関数Ｇｔｅｍｐ（ｔ）と組み合わせて表現して、再生の変化を観察し、そして当該変化の挙動を予想する。

構成は、ｄ_ｉｎｓｔ（ｔ）が次の公式によって確実に得られるように決定することもできる。

上の式では、
− ｍ_ｓｏｏｔ（ｔ）は時刻ｔにおいてフィルタに堆積した煤の質量の推定値であり、
− βは煤の燃焼に関する反応の順序であり、
− ｋ_０は一定の頻度因子であり、
− Ｅ_ａは、系の触媒活性によって変わる活性化エネルギーであり、この活性化エネルギーは、フィルタ内の煤燃焼反応の温度に対する感度を表わし、
− Ｒは完全ガス定数であり、
− Ｐ_０２（ｔ）は時刻ｔにおいて、フィルタに含まれるガス中の酸素の量を圧力の単位で表わした酸素分圧、すなわちＰ_Ｏ２（ｔ）＝Ｐ_Ｔ×Ｘ_Ｏ２であり、Ｐ_Ｔは合計圧力であり、そしてＸ_Ｏ２は酸素のモル分率であり、そして
− αは酸素に関する反応の順序である。

構成は、Ｇｔｅｍｐ（ｔ）が次の等式、
− Ｇｔｅｍｐ＝Ｔ_ｓ／Ｔ_ｅ（この場合、Ｔ_ｓ及びＴ_ｅはそれぞれ、時刻ｔにおける粒子フィルタのガス排出口温度、及びガス流入口温度である）
− Ｇｔｅｍｐ＝１／Ｔ_ｅ
− Ｇｔｅｍｐ＝１／（Ｔ_ｓ−Ｔ_ｅ）
− Ｇｔｅｍｐ＝（１／Ｔ_ｅ ^２）
− Ｇｔｅｍｐ＝１／Ｔ_ｉ（この場合、Ｔ_ｉは時刻ｔでのフィルタの内部温度である、及び
− Ｇｔｅｍｐ＝１／Ｔ_ｉ ^２
から成るグループに属する複数の公式の内のいずれか一つの公式によって得られるように決定することもできる。

Ｇｔｅｍｐに関するこれらの全ての式によって、Ｇｔｅｍｐを、所定時刻ｔでのフィルタの内部の温度の変化を表わす温度の関数として見積もることができる。
従って、代表温度はＴ_ｅまたはＴ_ｉ、あるいはＴ_ｓとすることができ、Ｔ_ｅは、粒子フィルタの上流の、すなわちフィルタの流入口でのガスの温度であり、Ｔ_ｉは、フィルタ内のガスの温度であり、Ｔ_ｓは、フィルタの排出口でのガスの温度である。

次に、Ｇｔｅｍｐ＝Ｔ_ｓ／Ｔ_ｅの場合、時刻ｔにおいて再生反応に伴うエネルギーの量を定義する温度比を考慮に入れる。この比を使用してフィルタの両端の温度差Ｔ_ｓ−Ｔ_ｅを定量化して、流入口と排出口との間の温度差に起因してフィルタにダメージを与える危険を評価することもできる。
Ｇｔｅｍｐ＝Ｔ_ｅの場合、フィルタの流入口でのガスの温度を考慮に入れるので、この温度Ｔ_ｅとフィルタの最大許容温度との差を求めることができる。温度Ｔ_ｅが高くなると、除去できる再生に伴う熱量が少なくなる。

次に、Ｔ_ｉが代表関数Ｇｔｅｍｐの分母として使用される場合、一方においてエンジンによって生成される熱、及び他方において再生処理中の煤の燃焼に伴う熱を考慮に入れる。
構成は、Ｄ_ｉｎｓｔｗｔ（ｔ）＝ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）が確実に成り立つように決定することもできる。この場合、ｄ_ｉｎｓｔ（ｔ）及びＧｔｅｍｐ（ｔ）のみを考慮に入れて再生の変化を観察する。

構成は、

が確実に成り立つように決定することもでき、この式では、ｔ１及びｔ２は固定期間である時間区間の限界値であり、時間と共に遷移する。
この関数によって、ｔ１とｔ２との間の時間区間に渡って時間ｔで積分した、ｄ_ｉｎｓｔ（ｔ）及びＧｔｅｍｐ（ｔ）の瞬時積の和を計算することにより、ｄ_ｉｎｓｔ（ｔ）及びＧｔｅｍｐ（ｔ）の経時変化が考慮され、この場合、ｔ１及びｔ２は時間と共に遷移する。

構成は、Ｄ_ｉｎｓｔｗｔ（ｔ）がｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）を時間ｔで微分した値に確実に等しくなるように決定することもできる。
この場合、ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）によって時間軸で定義される曲線の勾配が計算される。この勾配が急なほど、フィルタにダメージを与える危険が大きくなるが、これは、フィルタが非常に短い時間区間に渡って大きな熱変動及び温度変動に曝されることを意味するからである。

構成は、計算される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）が：
− 第１閾値を下回る場合に、再生が制御されることを示す信号が送出される処理、及び／又は、
− 前記第１閾値を上回り、かつ第２閾値を下回る場合に、フィルタを劣化させる危険があることを示す第１信号が送出される処理、及び／又は、
− 前記第２閾値を上回る場合に、再生が制御されないことを示す信号が送出される処理が確実に行なわれるように決定することもできる。
上に説明したように、所定の閾値を設定することにより、更にこれらの値を計算基準値Ｄ_ｉｎｓｔｗｔ（ｔ）と比較することにより、フィルタにダメージを与える危険が大きい、または小さいことを模式的に示す信号を送出することが可能になる。

構成は、確実に本方法が再生を中止するステップを含み、このステップにおいて、計算される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）が所定の中止閾値を上回る場合に、再生を中断するように決定することもできる。
この所定の中止閾値は、第１閾値と第２閾値との間の値であることが好ましく、かつ第１閾値と等しくすることができる。
この機能によって、再生を制御することができない状態になることによる粒子フィルタを劣化させる危険を小さくすることができ、再生を制御することができない状態は、過熱によって、そして／またはフィルタ内部の圧力／温度の急激かつ制御不能な上昇によって顕在化する。

構成は、計算される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）が時間ｔの経過と共に採る値が、確実にメモリに保存されるように決定することもできる。
これらの値をメモリに保存することにより、粒子フィルタの動作を表わす履歴を保持することが可能になる。このような履歴によって、粒子フィルタの磨耗の程度に関する指標が得られ、かつフィルタが誤動作する前にフィルタを使用するために残される平均時間に関する指標も得られる。

本発明の他の特徴及び利点は、以下に提示される記述から、完全に非制限的な例を通して、添付の図を参照することにより一層明確になると考えられる。

図１を参照すると、排気管２はエンジン（図示せず）に接続され、エンジンから燃焼ガス５が排気管に排出される。
フィルタ１には徐々に燃焼ガスに含まれる粒子３が堆積する。
コンピュータ４をセンサに接続して、粒子フィルタ１に関するデータを収集する。この場合、測定データはフィルタ上流のガス流入口温度Ｔ_e、フィルタ内のガスの内部温度Ｔ_ｉ，及びフィルタのガス排出口温度Ｔ_ｓである。これらの温度パラメータは、コンピュータ４に接続される３つの個別のセンサによって収集される。

他のパラメータも測定され、または推定され、そしてコンピュータ４に送信される。このような測定パラメータまたは推定パラメータは、エンジンから排出される排気ガスの質量流量、酸素分圧、及び粒子フィルタ内の煤の質量である。
これらのパラメータを使用して、コンピュータは、再生の変化を表わすパラメータｄ_ｉｎｓｔ（ｔ）を生成する。

このパラメータｄ_ｉｎｓｔ（ｔ）は、煤の燃焼速度Ｖ_ｒ（ｔ）、すなわち時刻ｔにおいて燃焼する煤の質量（例えば、グラム／秒、またはモル／秒）を、エンジンから排出され、そしてフィルタを通過する燃焼ガスの質量流量Ｑ_ｍ（ｔ）で割って得られる値を表わす、各時刻ｔにおける比である。
燃焼ガスの流れる量が大きくなると、フィルタから除去可能な熱量が大きくなり、従って、パラメータｄ_ｉｎｓｔ（ｔ）は質量流量Ｑ_ｍ（ｔ）によって最小化される。

一般的に、Ｑ_ｍ（ｔ）は時刻ｔにエンジンに吸入される空気の量の測定値、及びエンジンに噴射される燃料の量から推定される。Ｑ_ｍ（ｔ）の単位は、単位時間当たりに燃焼する煤の質量、すなわちｋｇ／ｓである。
従って、この比ｄ_ｉｎｓｔ（ｔ）は、時刻ｔにおいてフィルタから除去可能な熱量の大きさを示唆することができる。

本発明によれば、この比は、関数ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）によって定義される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）を計算することにより重み付けされ、この場合、Ｇｔｅｍｐ（ｔ）は代表温度Ｔ（ｔ）に従って時刻ｔとともに変化することができ、この代表温度は、例えばフィルタのガス流入口温度Ｔ_ｅ、及び／又はフィルタ内のガス温度Ｔ_ｉ、及び／又はフィルタのガス排出口温度Ｔ_ｓである。
詳細には、時刻ｔにおけるｄ_ｉｎｓｔ（ｔ）を、温度パラメータの関数として重み付けすることにより、フィルタの再生処理時の変化をより正確に表示することができる。

この重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）を計算するために、次の公式を使用する。

上の式では、Ｃｐ_ｅｘはケルビン／モルで表わされる排気ガスの比熱であり、そしてΔＨは、煤１モル当たりのジュール／ケルビンで表わされる煤の単位質量当たりの燃焼エンタルピーである。次の公式も使用する。

上の式では次式のような関係があり、Ｋは定数である。

上の式では、フィルタ内の煤の質量ｍ_ｓｏｏｔ（ｔ）はｇ（グラム）で表わされ、そしてコンピュータ４によって推定される値である。この質量は次のようにして推定することができる。
− エンジンが消費する燃料の量を分析し、そして／または
− フィルタの両端部での背圧を分析し、この場合、背圧は煤によるフィルタの目詰まりを表わし、従って背圧によってフィルタ内の煤の質量を間接的に推定する。

βは煤に関する反応の順序であり、研究室での実験により得られる。
ｋ_０は一定の頻度因子または前指数因子である。
Ｅ_ａは活性化エネルギーであり、系の触媒活性の関数である。このエネルギーは、炭素＋酸素→二酸化炭素で表わされる反応の温度に対する感度を表わす。Ｅ_ａは煤１モル当たりのジュールで表わされる。
Ｒは完全ガス定数であり、約８．３１４である。
Ｔ（ｔ）は時刻ｔでのフィルタ内の温度をケルビンで表わした値である。この温度Ｔ（ｔ）は内部温度であり、Ｔ_ｉとも表記される。
Ｐ_０２（ｔ）は時刻ｔでの酸素分圧であり、フィルタの流入口での燃焼ガス５に含まれる酸素の量を表わす。この量はセンサによって測定、またはエンジンパラメータに基づいて推定される。
αは酸素に関する反応の順序である。

Ｄ_ｉｎｓｔｗｔ（ｔ）を計算するために或る関数を使用するが、この関数では、Ｄ_ｉｎｓｔｗｔ（ｔ）はｄ_ｉｎｓｔ（ｔ）及びＧｔｅｍｐ（ｔ）によって表わされる。
この計算を行なうために、Ｇｔｅｍｐ（ｔ）は次式で表わされるＧｔｅｍｐから選択される。
− Ｇｔｅｍｐ＝１／Ｔ_ｅまたは、
− Ｇｔｅｍｐ＝１／（Ｔ_ｓ−Ｔ_ｅ）または、
− Ｇｔｅｍｐ＝１／Ｔ_ｅ ^２）または、
− Ｇｔｅｍｐ＝１／Ｔ_ｉ，この場合、Ｔ_ｉは時刻ｔでのフィルタの内部温度であり、または、
− Ｇｔｅｍｐ＝１／Ｔ_ｉ ^２

Ｄ_ｉｎｓｔｗｔ（ｔ）は次のようにして計算することができる。
− Ｄ_ｉｎｓｔｗｔ（ｔ）＝ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）
− または、

− または、ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）を時間ｔで微分する。
好ましくは、

であり、この式では、Ｔ_ｉはフィルタ内のガスの内部温度である。
後出の重み付けパラメータはフィルタの動作状態が極端な場合を表わす。

図２は、再生が制御される場合のパラメータの変動を示している。この場合、ｄ_ｉｎｓｔ（ｔ）は、再生が開始されるまで０に等しい。一旦、再生が開始されると、ｄ_ｉｎｓｔ（ｔ）は、最大値０．００１６に達するまで増え続ける。ｄ_ｉｎｓｔ（ｔ）のこの増加は、約１５０秒の時間区間に渡って発生し、そしてｄ_ｉｎｓｔ（ｔ）のこの最大値は、再生が制御されることを示している。同時に、フィルタのガス排出口温度は、約１００℃から当該温度が６６０℃に達するまで上昇する。この同じ時間区間に渡って、フィルタの流入口と排出口との間のミリバールで表わされる圧力差が小さくなり、フィルタを通過するガスの体積流量が増加する。時間が経過するにつれて、フィルタが徐々に目詰まりしていることが示される。

ｄ_ｉｎｓｔ（ｔ）の値が０．０２未満である限り、すなわちこの場合の図２に示す事例では、再生は制御され、そして粒子フィルタにダメージが加わる危険はない。
ｄ_ｉｎｓｔ（ｔ）の値が０．０２を超え、かつ０．０３未満である場合には、すなわち図３に示す事例では、フィルタ再生を制御することができない危険があり、熱量が高くなり過ぎることが分かる（７００℃に近い）。このような状況では、エンジン調整による介入を行なって、進行中の再生操作を停止する、または制限する必要がある場合が多い。
ｄ_ｉｎｓｔ（ｔ）の値が０．０３を超える場合には、すなわち図４に示す事例では、再生をもはや制御することができず、粒子フィルタにダメージが加わる危険が非常に大きくなることが分かる。

これは、パラメータｄ_ｉｎｓｔ（ｔ）が大きくなると、更に激しい反応がフィルタ内で発生するからである。従って、図４に示す事例では、１００秒未満の間に、ｄ_ｉｎｓｔ（ｔ）が０から０．０３４に変化する。同じ期間に渡って、フィルタから出て行くガスの温度Ｔ_ｓは、フィルタにダメージが加わる危険な温度にまで上昇し、７００℃の限界を超えて７８０℃にまで達する。
従って、本発明によれば、再生反応を解析することが更に容易になるようにｄ_ｉｎｓｔ（ｔ）に重み付けし、これにより図４に示す状況を回避することができる。

本発明の方法を実施するために適する装置を示す。 排気管の排気パラメータの変化を時間の関数として表わすグラフを示し、再生が正常に行なわれる様子を示す。 排気管の排気パラメータの変化を時間の関数として表わすグラフを示し、再生が正常状態から移行した様子を示す。 排気管の排気パラメータの変化を時間の関数として表わすグラフを示し、再生を制御することができず、粒子フィルタにダメージを与える危険が伴う様子を示す。

Claims (10)

1. 内燃エンジンの排気管に設けられる粒子フィルタの再生を制御する方法であって、前記方法は、
   − 粒子フィルタの動作データを求めるステップであって、これらのデータが時刻ｔでのフィルタ内の温度を表わす少なくとも一つの温度Ｔ（ｔ）である構成のステップと、
   − フィルタを再生するステップであって、フィルタの再生において、前記フィルタの温度を、前記フィルタに堆積した煤が燃焼する第１閾値温度よりも高い温度に上昇させるステップと、
   − 時間と共に変化することができ、かつ時刻ｔにおいて再生プロセスによって生成されるエネルギーの量を表わすパラメータｄ_ｉｎｓｔ（ｔ）を計算するステップとを含み、
   前記方法は、関数ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）によって少なくとも部分的に定義される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）を計算するステップを含み、Ｇｔｅｍｐ（ｔ）は時刻（ｔ）とともに変化することができ、かつ前記代表温度Ｔ（ｔ）の関数として表わされることを特徴とする、方法。
2. パラメータｄ_ｉｎｓｔ（ｔ）は次の公式：
   

   

   により得られることを特徴とし、上の式では、Ｖ_ｒ（ｔ）は時刻ｔにおいて粒子フィルタに堆積した煤の燃焼速度であり、Ｑ_ｍ（ｔ）は時刻ｔにおいて推定、または測定される排気ガス質量流量である、請求項１記載の方法。
3. ｄ_ｉｎｓｔ（ｔ）は次の公式：
   

   

   により得られることを特徴とし、上の式では、
   − ｍ_ｓｏｏｔ（ｔ）は時刻ｔにおいてフィルタに堆積した煤の質量の推定値であり、
   − βは煤の燃焼に関する反応の順序であり、
   − ｋ_０は一定の頻度因子であり、
   − Ｅ_ａは、系の触媒活性によって変わる活性化エネルギーであり、この活性化エネルギーは、フィルタ内の煤燃焼反応の温度に対する感度を表わし、
   − Ｒは完全ガス定数であり、
   − Ｐ_０２（ｔ）は時刻ｔでの酸素分圧であり、時刻ｔにおいてフィルタに含まれるガスの中の酸素の質量割合として表わされ、そして
   − αは酸素に関する反応の順序である、
   請求項１及び２のいずれか一項に記載の方法。
4. Ｇｔｅｍｐ（ｔ）は、次の等式：
   − Ｇｔｅｍｐ＝Ｔ_ｓ／Ｔ_ｅ（この場合、Ｔ_ｓ及びＴ_ｅはそれぞれ、時刻ｔにおける粒子フィルタのガス排出口温度、及びガス流入口温度である）
   − Ｇｔｅｍｐ＝１／Ｔ_ｅ
   − Ｇｔｅｍｐ＝１／（Ｔ_ｓ−Ｔ_ｅ）
   − Ｇｔｅｍｐ＝（１／Ｔ_ｅ ^２）
   − Ｇｔｅｍｐ＝１／Ｔ_ｉ，（この場合、Ｔ_ｉは時刻ｔでのフィルタの内部温度である）、及び
   − Ｇｔｅｍｐ＝１／Ｔ_ｉ ^２、
   から成るグループに属する複数の公式の内のいずれか一つの公式によって得られることを特徴とする、請求項１〜３のいずれか一項に記載の方法。
5. Ｄ_ｉｎｓｔｗｔ（ｔ）＝ｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）が成り立つことを特徴とする、請求項１〜４のいずれか一項に記載の方法。
6. 次の式：
   

   

   が成り立つことを特徴とし、この式では、ｔ１及びｔ２は固定期間である時間区間の限界値であり、時間と共に遷移することを特徴とする、請求項１〜４のいずれか一項に記載の方法。
7. Ｄ_ｉｎｓｔｗｔ（ｔ）はｄ_ｉｎｓｔ（ｔ）×Ｇｔｅｍｐ（ｔ）を時間ｔで微分した値に等しいことを特徴とする、請求項１〜４のいずれか一項に記載の方法。
8. 計算される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）が、
   − 第１閾値を下回る場合に、再生が制御されることを示す信号が送出され、そして／または、
   − 前記第１閾値を上回り、かつ第２閾値を下回る場合に、フィルタを劣化させる危険があることを示す第１信号が送出され、そして／または
   − 前記第２閾値を上回る場合に、再生が制御されないことを示す信号が送出されることを特徴とする、請求項１〜７のいずれか一項に記載の方法。
9. 前記方法が、再生を中止するステップを含み、このステップでは、計算される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）が所定の中止閾値を上回る場合に、再生を中断することを特徴とする、請求項１〜８のいずれか一項に記載の方法。
10. 計算される重み付けパラメータＤ_ｉｎｓｔｗｔ（ｔ）が時間ｔの経過と共に採る値がメモリに保存されることを特徴とする、請求項１〜９のいずれか一項に記載の方法。

Images