JP3688533B2 - 排ガス浄化用触媒装置の劣化状態評価方法 - Google Patents
排ガス浄化用触媒装置の劣化状態評価方法 Download PDFInfo
- Publication number
- JP3688533B2 JP3688533B2 JP32242899A JP32242899A JP3688533B2 JP 3688533 B2 JP3688533 B2 JP 3688533B2 JP 32242899 A JP32242899 A JP 32242899A JP 32242899 A JP32242899 A JP 32242899A JP 3688533 B2 JP3688533 B2 JP 3688533B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- fuel ratio
- air
- catalyst device
- exhaust gas
- 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
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01N—GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR MACHINES OR ENGINES IN GENERAL; GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR INTERNAL COMBUSTION ENGINES
- F01N11/00—Monitoring or diagnostic devices for exhaust-gas treatment apparatus, e.g. for catalytic activity
-
- 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/1401—Introducing closed-loop corrections characterised by the control or regulation method
-
- 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/1439—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the position of the sensor
- F02D41/1441—Plural sensors
-
- 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/146—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 an NOx content or concentration
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01N—GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR MACHINES OR ENGINES IN GENERAL; GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR INTERNAL COMBUSTION ENGINES
- F01N2550/00—Monitoring or diagnosing the deterioration of exhaust systems
- F01N2550/02—Catalytic activity of catalytic converters
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01N—GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR MACHINES OR ENGINES IN GENERAL; GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR INTERNAL COMBUSTION ENGINES
- F01N2560/00—Exhaust systems with means for detecting or measuring exhaust gas components or characteristics
- F01N2560/02—Exhaust systems with means for detecting or measuring exhaust gas components or characteristics the means being an exhaust gas sensor
- F01N2560/023—Exhaust systems with means for detecting or measuring exhaust gas components or characteristics the means being an exhaust gas sensor for measuring or detecting HC
-
- 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/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1413—Controller structures or design
- F02D2041/143—Controller structures or design the control loop including a non-linear model or compensator
-
- 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/1454—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 an oxygen content or concentration or the air-fuel ratio
- F02D41/1456—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 an oxygen content or concentration or the air-fuel ratio with sensor output signal being linear or quasi-linear with the concentration of oxygen
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Description
【発明の属する技術分野】
本発明は、内燃機関の排気通路に設けた排ガス浄化用触媒装置の劣化状態を評価する方法に関する。
【0002】
【従来の技術】
内燃機関では、一般に、三元触媒等により構成される排ガス浄化用触媒装置が排気通路に設けられ、この触媒装置により、排ガス中に含まれるHC(ハイドロカーボン)やNOx(窒素酸化物)等のガス成分(以下、ここでは浄化成分という)を浄化するようにしている。
【0003】
尚、この種の触媒装置による各種類の浄化成分の浄化率は、基本的には触媒装置の上流側の排ガスの空燃比の空燃比の値がある値もしくはその近傍の値であるときにほぼ最大(極大)となり、そのような空燃比の値から実際の空燃比がリーン側あるいはリッチ側にずれると、上記浄化成分の浄化率は低下する傾向がある。ここで、前記触媒装置の上流側の排ガスの空燃比は、より詳しくは、該触媒装置に進入する排ガスの酸素濃度から把握される空燃比で、内燃機関で燃焼して排ガスとなった混合気の空燃比を表すものである(以下、ここでは、この空燃比を単に内燃機関の空燃比という)。
【0004】
この種の触媒装置は、その経時的な使用により徐々に劣化し、浄化性能が低下していく。このため、触媒装置による排ガスの適正な浄化状態を確保するためには、該触媒装置がある程度劣化したときに、該触媒装置を新品のものに交換する等の処置をとる必要がある。そして、このような処置をとるためには、触媒装置の劣化状態を適正に評価し得る技術が望まれる。
【0005】
一方、近年では、HCやNOx等、触媒装置により浄化される各種類のガス成分(以下、ここでは浄化成分という)の濃度を比較的精度よく検出することが可能な排ガスセンサ(具体的にはHCセンサやNOxセンサ)が開発されてきている。尚、この種の排ガスセンサは、検出する浄化成分の濃度の増加に伴い、ほぼリニアに該センサの出力レベルが大きくなるものが一般的であるが、これと逆に、浄化成分の濃度の増加に伴い、該センサの出力レベルが小さくなっていくものもある。以下の説明では、前者のタイプの排ガスセンサの出力特性をポジティブ特性、後者のタイプの排ガスセンサの出力特性をネガティブ特性ということがある。
【0006】
このような排ガスセンサを用いると、触媒装置が実際に浄化するガス成分、すなわち、該触媒装置の劣化状態が直接的に影響するガス成分の濃度を検出することができるため、触媒装置の劣化状態を把握するために有用に活用することが可能であると考えられる。
【0007】
例えば、このような排ガスセンサを触媒装置の上流側及び下流側にそれぞれ配置し、内燃機関の運転中のそれらの排ガスセンサの出力の比をとることで、HC、NOx等、該排ガスセンサが検出するガス成分の触媒装置による浄化率を直接的に把握することが考えられる。この場合、触媒装置により浄化されるガス成分の浄化率は、触媒装置の劣化の進行に伴って減少する。従って、そのガス成分の浄化率を上記のように把握すれば、触媒装置の劣化状態を適正に判別することが可能であると考えられる。
【0008】
しかしながら、この手法では、触媒装置が有する応答遅れ等の影響で、触媒装置によるガス成分(排ガスセンサが検出するガス成分)の浄化率を実際に精度よく把握することは難しい。また、触媒装置が有する応答遅れ等の影響に対する対策を施して、把握し得る浄化率の精度を高めようとすると、そのための処理が煩雑なものとなりやすい。
【0009】
このため、上記のような排ガスセンサ(HCセンサやNOxセンサ)を用いた比較的簡易な手法で、適正に触媒装置の劣化状態を判別し得る新規な技術が望まれていた。
【0010】
【発明が解決しようとする課題】
本発明はかかる背景に鑑みてなされたものであり、排ガス浄化用触媒装置により浄化するHCやNOx等のガス成分の濃度を検出する排ガスセンサを用いた簡易な手法によって、触媒装置の劣化状態を適正に評価することができる排ガス浄化用触媒装置の劣化状態評価方法を提供することを目的とする。
【0011】
【課題を解決するための手段】
本発明の排ガス浄化用触媒装置の劣化状態評価方法は、かかる目的を達成するために、内燃機関の排気通路に設けた排ガス浄化用の触媒装置の下流側に、前記触媒装置を通過した排ガスのうちの該触媒装置により浄化される特定成分の濃度を検出する排ガスセンサを配置し、内燃機関の運転中における前記排ガスセンサの出力のデータを用いて前記触媒装置の劣化状態を評価する方法であって、前記排ガスセンサが検出する前記特定成分の濃度の、前記触媒装置の上流側の排ガスの空燃比に対する相関特性を、該空燃比の変化に対して該特定成分の濃度が極小値を持つ相関特性として表現すべくあらかじめ定めた形式の非線形関数に含まれるパラメータの値を、前記触媒装置の上流側の排ガスの空燃比を表すデータと前記排ガスセンサの出力のデータを用いて同定する同定工程と、その同定したパラメータの値に基づき、前記触媒装置の劣化状態を評価する評価工程とを備えたことを特徴とする。
【0012】
尚、前記非線形関数は、詳しくは、前記触媒装置の上流側の排ガスの空燃比を表すもの(例えば該空燃比を検出するセンサの出力やその出力から把握される空燃比、あるいは該空燃比の目標値等)を独立変数、前記排ガスセンサの出力やその出力から把握される前記特定成分の濃度を従属変数とする関数である。そして、本発明に適した非線形関数の形式としては例えば二次関数や三次関数等の高次関数が挙げられる。また、前記パラメータは、詳しくは、前記非線形関数のグラフ形状を規定するパラメータである。例えば、該非線形関数の形式を二次関数や三次関数等の高次関数としたとき、その各次数の項に係る係数や定数項が前記パラメータである。
【0013】
また、触媒装置の上流側の排ガスの空燃比は詳しくは、前述の通り該排ガスの酸素濃度から把握される空燃比である。以下の本発明の説明では、この空燃比を単に内燃機関の空燃比ということがある。
【0014】
また、前記排ガスセンサが検出する前記特定成分は、具体的にはHC(ハイドロカーボン)やNOx(窒素酸化物)等である。
【0015】
かかる本発明によれば、前記同定工程において、前記内燃機関の空燃比(触媒装置の上流側の排ガスの空燃比)を表すデータと、排ガスセンサの出力のデータ(前記特定成分の濃度の検出値のデータ)を用いて前記非線形関数のパラメータの値を同定することで、前記排ガスセンサが検出する前記特定成分の濃度の、前記内燃機関の空燃比に対する相関特性を近似的に数式表現する非線形関数が確定される(非線形関数のグラフ形状が確定される)こととなる。
【0016】
この場合、前記排ガスセンサが検出する前記特定成分の濃度(触媒装置による浄化後の特定成分の濃度)は、一般に、内燃機関の空燃比の変化に対して極小値をもつ傾向がある。そして、そのような相関特性は例えば二次関数等、適切な形式の非線形関数により的確に表現することができる。
【0017】
上記のように同定工程で非線形関数のパラメータの値を同定したとき、その同定値により規定される非線形関数のグラフ形状等は、触媒装置の劣化の進行に伴い、ある特徴的な変化を呈する傾向がある。つまり、前記非線形関数のパラメータの同定値には、触媒装置の劣化状態が反映される。そこで、本発明では、前記評価工程は、前記同定工程で同定されたパラメータの値に基づき、前記触媒装置の劣化状態を評価する。
【0018】
このような本発明によれば、前記非線形関数のパラメータの同定値は統計的な処理によって求められることとなるので、該非線形関数のパラメータの同定値は触媒装置の応答遅れ特性等の影響を受けにくく、触媒装置の劣化状態が適正に反映される。また、該非線形関数のパラメータの値は、最小二乗法等、既知の手法によって同定することが可能である。従って、本発明によれば簡便な手法によって、適正に触媒装置の劣化状態を評価することが可能となる。
【0019】
かかる本発明において、前記同定工程で前記非線形関数のパラメータの値の同定のために用いる前記内燃機関の空燃比(触媒装置の上流側の排ガスの空燃比)を表すデータは、例えば該空燃比の値そのものを表すデータであってもよいが、該内燃機関の空燃比と所定の基準値との偏差を表すデータであることが好ましい。
【0020】
このように内燃機関の空燃比と所定の基準値との偏差を表すデータを用いて前記非線形関数のパラメータの値を同定することで、その同定値の精度を高めることができる。
【0021】
また、本発明では、例えば前記内燃機関の空燃比を所要の目標空燃比に制御するような場合には、前記非線形関数のパラメータの値の同定のために、前記内燃機関の空燃比として上記目標空燃比を使用することも可能であるが、好ましくは、前記内燃機関の空燃比を検出する空燃比センサを該触媒装置の上流側に設ける。そして、前記同定工程は、前記非線形関数のパラメータの値の同定のために、前記内燃機関の空燃比として前記空燃比センサにより検出される空燃比を用いる。
【0022】
このように空燃比センサにより検出される空燃比、すなわち、内燃機関の空燃比(触媒装置の上流側の排ガスの空燃比)の実際の値を用いて前記非線形関数のパラメータの値を同定することで、その同定値を精度よく求めることができる。
【0023】
そして、このように空燃比センサを備えた本発明では、前記同定工程で前記パラメータの値を同定してなる前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値を該非線形関数のパラメータの同定値を用いて求め、その求めた空燃比の値を前記触媒装置の上流側の排ガスの目標空燃比として得る目標空燃比算出工程と、該目標空燃比に前記空燃比センサが検出する空燃比を収束させるようにフィードバック制御により前記内燃機関で燃焼させる混合気の空燃比を操作する空燃比操作工程とを備える。
【0024】
すなわち、本発明では、前述のように同定されるパラメータの値により規定される非線形関数は、前記相関特性、すなわち、排ガスセンサが検出する前記特定成分の濃度の前記内燃機関の空燃比に対する相関特性を近似的に表現するものとなる。このため、上記のように内燃機関の目標空燃比(触媒装置の上流側の排ガスの目標空燃比)を求めたとき、その目標空燃比は、基本的には排ガスセンサが検出する特定成分の濃度、すなわち、触媒装置による浄化後の該特定成分の濃度が極小値となるような目標空燃比となる。従って、この目標空燃比に前記空燃比センサが検出する空燃比を収束させるようにフィードバック制御により内燃機関で燃焼させる混合気の空燃比を操作したとき、触媒装置による前記特定成分の浄化率が極大となるような空燃比状態に前記内燃機関の空燃比が制御されることとなる。これにより、触媒装置の劣化状態の評価を行ないながら、それと並行して触媒装置による前記特定成分の良好な浄化(特定成分の浄化率が最大となるような浄化)を行うことが可能となる。
【0025】
尚、内燃機関の空燃比のフィードバック制御は、適応制御器等の漸化式形式のフィードバック制御手段により行なうことが好適である。ここで、漸化式形式のフィードバック制御手段は、空燃比のフィードバック操作量(例えば燃料供給量の補正量)の現在以前の所定数の時系列データを含む所定の漸化式によって、新たなフィードバック操作量を求めるものである。
【0026】
また、本発明では、前記非線形関数のパラーメータの同定は、例えば、前記内燃機関の空燃比を表すデータ(前記空燃比センサの出力等)と排ガスセンサの出力のデータを収集・蓄積した上で行なうようにすることも可能であるが、前記同定工程は、逐次型同定アルゴリズムにより前記非線形関数のパラメータの値を逐次同定する。
【0027】
このように逐次型同定アルゴリズムを用いることで、そのアルゴリズムの実行に要するメモリ容量が少なくて済む。さらに、非線形関数のパラメータがリアルタイムで逐次更新されることから、瞬時瞬時の触媒装置や内燃機関の挙動状態に則して非線形関数のパラメータの同定値をリアルタイムで得ることが可能となる。その結果、該非線形関数のパラメータの同定値に基づく触媒装置の劣化状態の評価結果の信頼性を高めることができる。
【0028】
尚、前記逐次型同定アルゴリズムとしては、逐次型最小二乗法、逐次型重み付き最小二乗法、固定ゲイン法、漸減ゲイン法等のアルゴリズムが挙げられる。これらのアルゴリズムによれば、前記パラメータの現在の同定値を用いて前記非線形関数により求められる排ガスセンサの出力の値と、該排ガスセンサの出力の実際の値との間の誤差を最小化するように新たな前記パラメータの値が求められる(パラメータの値が更新される)。
【0029】
このように逐次型同定アルゴリズムを用いる本発明では、前記同定工程は、前記非線形関数のパラメータのうちの少なくとも一つのパラメータの値を所定の条件を満たす値に制限しつつ同定することが好ましい。
【0030】
すなわち、非線形関数のパラメータのうちには、その値が外乱等の影響で本来あるべき値の範囲からはずれた値に誤って同定されると、実際の前記相関特性のグラフ形状と、非線形関数のグラフ形状とが大きく異なるものとなることがある。そこで、本発明では、上記のように非線形関数のパラメータのうちの少なくとも一つを同定するに際しては、そのパラメータの値をあらかじめ定めた条件を満たす値に制限する。これにより、該パラメータの値の信頼性を確保し、該パラメータの値から求められる評価パラメータの信頼性を高めることができる。
【0031】
例えば非線形関数を二次関数とした場合、前記同定工程では、少なくとも該二次関数の最大次数の項に係る係数をその値を制限する前記パラメータとし、該二次関数の関数値により表される前記特定成分の濃度が極小値をもつように前記係数の値を制限しつつ同定することが好ましい。
【0032】
すなわち、二次関数は、その最大次数の項に係る係数の値が正であるか負であるかによって、その関数値が極小値をもつか極大値をもつかが異なる。従って、上記係数の同定値の正負の極性が不適正であると、前記非線形関数としての二次関数の関数値により表される前記特定成分の濃度が、本来もつべき極小値をもたずに、極大値をもつものとなってしまう。このため、二次関数の関数値により表される前記特定成分の濃度が極小値をもつように前記係数の同定値を制限する。具体的には、該係数の同定値を正又は負のいずれか一方の極性の値に制限する。これにより、前記パラメータの同定値、ひいては前記評価パラメータの信頼性を高めることが可能となる。
【0033】
尚、前記非線形関数の関数値により表される特定成分の濃度が極小値をもつとき、該非線形関数の関数値そのものは、極小値をもつ場合及び極大値をもつ場合のいずれの場合もありうる。すなわち、前記非線形関数の従属変数を前記排ガスセンサの出力から把握される特定成分の濃度とした場合、非線形関数の関数値により表される特定成分の濃度が極小値となる空燃比の値において、該非線形関数の関数値そのものも極小値となる。また、非線形関数の従属変数を、例えば前記排ガスセンサの出力とした場合にあっては、該排ガスセンサの出力の特定成分の濃度に対する特性が前述のポジティブ特性(濃度の増加に伴い出力が大きくなる特性)である場合には、非線形関数の関数値により表される特定成分の濃度が極小値となる空燃比の値において、非線形関数の関数値も極小値となる。一方、非線形関数の従属変数を、前記排ガスセンサの出力とした場合において、該排ガスセンサの出力の特定成分の濃度に対する特性が前述のネガティブ特性(濃度の増加に伴い出力が小さくなる特性)である場合には、非線形関数の関数値が表す特定成分の濃度は、該関数値が大きい程、小さくなる。このため、この場合には、該非線形関数の関数値が表す特定成分の濃度が極小値となる空燃比の値において、非線形関数の関数値そのものは極大値をとることとなる。
【0034】
従って、上記のように二次関数の最大次数の項の係数の値(同定値)を制限する場合において、非線形関数としての二次関数の従属変数を排ガスセンサの出力から把握される特定成分の濃度、あるいは前述のポジティブ特性の排ガスセンサの出力とした場合には、該二次関数の最大次数の項の係数の値は正の値に制限し、該二次関数の従属変数を前述のネガティブ特性の排ガスセンサの出力とした場合には、該係数の値は、負の値に制限することとなる。、
かかる本発明では、前記非線形関数のパラメータの同定値に基づく触媒装置の劣化状態の評価はより具体的には、次のように行われる。
【0035】
すなわち、前記評価工程は、前記触媒装置の劣化の進行に応じた変化を呈するものとして前記非線形関数のパラメータを用いてあらかじめ定義した評価パラメータの値を、前記同定工程により同定された前記非線形関数のパラメータの値から求める工程を備え、該評価パラメータの値に基づき、前記触媒装置の劣化状態を評価する。
【0036】
上記のように非線形関数のパラメータの同定値により定まる評価パラメータを導入したとき、該評価パラメータの値が触媒装置の劣化の進行に応じた変化(劣化の進行に伴い値が増加もしくは減少するなど)を呈する。このため、該評価パラメータの値を、あらかじめ定めた所定値と比較する等することによって、容易に触媒装置の劣化状態を評価することができる。
【0037】
この場合、より具体的には、前述のようにパラメータの値を同定する前記非線形関数は、触媒装置の劣化が進行するに伴い、該非線形関数のグラフの全体が空燃比のリーン側に変移していく傾向がある。
【0038】
従って、本発明における前記評価パラメータの一つの例としては、前記同定工程で前記パラメータの値を同定してなる前記非線形関数のグラフが前記触媒装置の劣化の進行に伴い空燃比のリーン側に変移していくことに呼応して変化する評価パラメータ(以下、ここでは第1評価パラメータという)が挙げられる。
【0039】
そして、このような第1評価パラメータとしては、より具体的には、例えば前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値が挙げられる。
【0040】
すなわち、前述のように、前記排ガスセンサが検出する特定成分の濃度が内燃機関の空燃比の変化に対して一般に極小値をもつことに対応して、前記非線形関数の関数値により表される特定成分の濃度も極小値をもつ。そして、前記のように触媒装置の劣化の進行に伴い、非線形関数のグラフが空燃比のリーン側に変移していくと、前記非線形関数の関数値により表される特定成分の濃度が極小値となるような空燃比の値もリーン側の値に変化していくこととなる。
【0041】
従って、このような空燃比の値が前記第1評価パラメータとして採用し得る。そして、このような第1評価パラメータに基づけば、触媒装置の劣化状態を適正に評価することができる。
【0042】
また、前述のようにパラメータの値を同定する前記非線形関数は、触媒装置の劣化が進行するに伴い、該非線形関数のグラフの全体が、該非線形関数の関数値により表される特定成分の濃度が増加する方向に変移する傾向がある。これは、触媒装置の劣化の進行に伴い、前記特定成分の浄化率が低下し、排ガスセンサが検出する特定成分の濃度が増加するからである。
【0043】
従って、本発明における前記評価パラメータの他の一つの例としては、前記同定工程で前記パラメータの値を同定してなる前記非線形関数のグラフが前記触媒装置の劣化の進行に伴い該非線形関数の関数値により表される前記特定成分の濃度が増加する方向に変移していくことに呼応して変化する評価パラメータ(以下、ここでは第2評価パラメータという)が挙げられる。
【0044】
そして、このような第2評価パラメータとしては、より具体的には、例えば前記非線形関数の関数値により表される前記特定成分の濃度の極小値に対応する該非線形関数の関数値の極値(極小値又極大値)が挙げられる。
【0045】
すなわち、触媒装置の劣化の進行に伴い、前記非線形関数の関数値により表される特定成分の濃度が増加するため、その濃度の極小値に対応する非線形関数の関数値の極値も同様の変化を生じることとなる。この場合、非線形関数の極値が極小値である場合には、該極値が触媒装置の劣化の進行に伴い、大きくなっていく。また、該極値が極大値である場合には、該極値は、触媒装置の劣化の進行に伴い小さくなっていく。
【0046】
従って、このような非線形関数の極値が前記第2評価パラメータとして採用し得る。そして、このような第2評価パラメータに基づけば、触媒装置の劣化状態を適正に評価することができる。
【0047】
また、前述のようにパラメータの値を同定する前記非線形関数にあっては、触媒装置の劣化が進行するに伴い、該非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値よりもリーン側又はリッチ側の空燃比域(基本的には両側の空燃比域)において該非線形関数のグラフの傾きが前記触媒装置の劣化の進行に伴い大きくなる。これは、触媒装置の劣化が進行すると、前記特定成分を良好に浄化し得る内燃機関の空燃比域が狭まるためである。
【0048】
従って、本発明における前記評価パラメータのさらに他の一つの例としては、前記同定工程で前記パラメータの値を同定してなる前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値よりもリーン側又はリッチ側の空燃比域において該非線形関数のグラフの傾きが前記触媒装置の劣化の進行に伴い大きくなることに呼応して変化する評価パラメータ(以下、ここでは第3評価パラメータという)が挙げられる。
【0049】
そして、このような第3評価パラメータとしては、より具体的には、例えば前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値からあらかじめ定めた所定量だけリーン側又はリッチ側に変化させた空燃比の値における前記非線形関数の関数値が挙げられる。
【0050】
すなわち、前記のように非線形関数のグラフの傾きが触媒装置の劣化の進行に伴い大きくなると、前記特定成分の濃度が極小値となる空燃比の値からあらかじめ定めた所定量だけリーン側又はリッチ側に変化させた空燃比の値における前記非線形関数の関数値は増加、あるいは減少していく。この場合、特定成分の濃度の極小値に対応する非線形関数の極値が極小値である場合には、上記の関数値は触媒装置の劣化の進行に伴い増加する。また、非線形関数の極値が極大値である場合には、上記の関数値は、触媒装置の劣化の進行に伴い減少する。
【0051】
従って、このような非線形関数の関数値が前記第3評価パラメータとして採用し得る。そして、このような第3評価パラメータに基づけば、触媒装置の劣化状態を適正に評価することができる。
【0052】
かかる本発明では、例えば上記の第1〜第3評価パラメータのいずれか一つのみに基づいて触媒装置の劣化状態を評価してもよいが、好ましくは、前記評価工程で求める評価パラメータを複数種類とする(例えば上記の第1〜第3の各評価パラメータを求める)。そして、該評価工程では、各種類の評価パラメータに基づく前記触媒装置の劣化状態の評価結果を総合して、該触媒装置の劣化状態の評価結果を確定することが好ましい。
【0053】
このようにすることで、触媒装置の劣化状態の評価結果の信頼性を高めることができる。
【0054】
【発明の実施の形態】
本発明の排ガス浄化用触媒装置の劣化状態評価方法の第1の実施形態を図1〜図15を参照して説明する。
【0055】
図1は本実施形態における装置の全体的システム構成を示すブロック図である。図中、1は例えば自動車あるいはハイブリッド車に車両の推進源として搭載された4気筒エンジン(内燃機関)である。このエンジン1の各気筒毎に燃料及び空気の混合気の燃焼により生成される排ガスは、エンジン1の近傍で共通の排気管2(排気通路)に集合され、この排気管2を介して大気側に放出される。
【0056】
排気管2には、例えば三元触媒により構成された排ガス浄化用の触媒装置3が介装されている。この触媒装置3は、排ガス中のHC(ハイドロカーボン)やNOx(窒素酸化物)を酸化/還元作用等により浄化する。
【0057】
また、排気管2には、触媒装置3の上流側(詳しくはエンジン1の各気筒毎の排ガスの集合箇所)で空燃比センサ4が装着され、さらに触媒装置3の下流側で排ガスセンサ5が装着されている。
【0058】
空燃比センサ4(以下、LAFセンサ4という)は、エンジン1で燃焼した混合気の空燃比、詳しくはエンジン1の各気筒毎の排ガスの集合箇所における排ガス中の酸素濃度により把握される空燃比(以下、単にエンジン1の空燃比という)の検出値を表す出力KACTを生成するものである。このLAFセンサ4は、例えば本願出願人が特開平4−369471号公報にて詳細に説明した広域空燃比センサにより構成されたものである。このLAFセンサ4の出力KACTは、排ガス中の酸素濃度に対応した空燃比の比較的広範囲にわたって、それに比例したレベルのものとなる(空燃比の変化に対して出力KACTがリニアに変化する)。
【0059】
排ガスセンサ5は、本実施形態では、触媒装置3を通過した排ガス中のNOx濃度(触媒装置3によるNOxの浄化後の濃度)の検出値を表す出力VSを生成するNOxセンサである。この場合、排ガスセンサ5は、その出力VSが、図2に実線で示すように、NOx濃度の増加に伴い、それにほぼ比例して該出力VSのレベルが増加していく特性、すなわち、前述のポジティブ特性を有するものである。
【0060】
尚、NOxセンサには、図2に仮想線で示すように、NOx濃度の増加に伴い、NOxセンサの出力がリニアに減少していくネガティブ特性を有するものもあるが、本実施形態では、実線で示すポジティブ特性を有するものである。
【0061】
本実施形態の装置は、基本的には、触媒装置3によるNOxの浄化率が該触媒装置3の劣化状態等によらずに最大(極大)となるように、換言すれば、排ガスセンサ5の出力VSが表すNOx濃度(触媒装置3による浄化後のNOx濃度)が最小(極小)となるようにエンジン1の空燃比を操作する制御を行う。また、この制御と並行して、触媒装置3の劣化状態の評価する処理を逐次行なう。
【0062】
このような制御及び処理を行なうために、本実施形態の装置は、主演算処理制御器6及び燃料供給制御器7を備えている。
【0063】
ここで、主演算処理制御器6は、前記LAFセンサ4の出力KACT及び排ガスセンサ5の出力VSのサンプリングデータを用いて後述の二次関数(非線形関数)のパラメータの値を逐次同定する同定手段6aとしての機能と、この同定手段6aによって同定されたパラメータの値を用いて排ガスセンサ5の出力VSが表すNOx濃度が最小(極小)となるようなエンジン1の目標空燃比KCMD(これはLAFセンサ4の出力KACTの目標値でもある)を逐次算出する目標空燃比算出手段6bとしての機能と、同定手段6aによって同定されたパラメータの値を用いて触媒装置3の劣化状態を評価する処理を行なう劣化状態評価手段6cとしての機能とを具備するものである。
【0064】
そして、主演算処理制御器6には、その劣化状態評価手段6cによる触媒装置3の劣化状態の評価結果に応じた報知を行なう劣化報知器8が接続されている。この劣化報知器8は、例えばランプやブザー、あるいは文字、図形等を表示可能な表示器等により構成されたものである。
【0065】
また、燃料供給制御器7は、主演算処理制御器6がその目標空燃比算出手段6bとして機能によって求めた目標空燃比KCM DにLAFセンサ4の出力KACT(エンジン1の空燃比の検出値)を収束させるようにエンジン1の燃料供給量(燃料噴射量)を調整することで、エンジン1の空燃比を操作する空燃比操作手段としての機能を有するものである。
【0066】
尚、燃料供給制御器7には、LAFセンサ4の出力KACTの他、排ガスセンサ5の出力VSや、エンジン1の回転数、吸気圧(吸気管内圧)、冷却水温等を検出するための図示しない各種センサの出力も与えられるようになっている。
【0067】
これらの制御器6,7は、マイクロコンピュータを用いて構成されたもので、それぞれの制御処理を所定の制御サイクルで実行する。この場合、本実施形態では、燃料供給制御器7が実行する処理(燃料噴射量の調整処理)の制御サイクルは、エンジン1のクランク角周期(所謂TDC)に同期した周期とされている。また、目標空燃比算出制御器6が実行する目標空燃比の算出処理の制御サイクルは、あらかじ定めた一定周期(前記クランク角周期よりも長い周期)としている。
【0068】
ここで、これらの制御器6,7をさらに説明する前に、触媒装置3によるNOxの浄化特性を説明しておく。
【0069】
触媒装置3は、エンジン1の空燃比の変化に対して基本的には図3に示すようなNOxの浄化特性(NOxの浄化率の特性)を有する。図3の実線は、新品の触媒装置3に関するグラフである。また、図3の破線及び一点鎖線は、ある程度劣化が進行した触媒装置3に関するグラフである。さらに詳しくは、破線のグラフは、劣化の進行度合いが比較的小さい触媒装置3に関するもの、一点鎖線のグラフ劣化の進行度合いが比較的大きい触媒装置3に関するものである。
【0070】
図3に実線のグラフで示すように、新品の触媒装置3によるNOxの浄化率は、エンジン1の空燃比がある値AF1よりもリッチ側の空燃比域にあるときには、ほぼ一定(略100%)の最大の浄化率となる。そして、エンジン1の空燃比が上記の値AF1よりもリーン側の空燃比域になると、NOxの浄化率が急激に低下する。
【0071】
また、触媒装置3の劣化がある程度進行すると、図3の破線や一点鎖線のグラフで示すように、NOxの浄化率は、ある値AF2(:破線),AF3(:一点鎖線)の空燃比において極大値をもつ(ピークをもつ)。そして、エンジン1の空燃比がNOxの浄化率の極大値に対応する値AF2,AF3から、リッチ側及びリーン側のいずれに変化しても、NOxの浄化率が減少していく。この場合、NOxの浄化率の極大値からの減少は、リーン側では比較的急激に生じ(グラフの傾きが大きい)、リッチ側では、比較的緩やかに生じる(グラフの傾きが小さい)。尚、触媒装置3が劣化した状態におけるNOxの浄化率の極大値は、触媒装置3の新品状態における浄化率の最大値よりも小さい。また、該極大値は、触媒装置3の劣化の進行に伴い小さくなる。
【0072】
このように触媒装置3によるNOxの浄化率は、触媒装置3の新品状態を除いて、エンジン1の空燃比の変化に対して基本的には極大値を持つような特性を呈する。
【0073】
尚、触媒装置3の新品状態を含めた該触媒装置3の各劣化状態において、NOxの浄化率が最大(極大)となるような空燃比の値(前記AF1,AF2,AF3等)は、該触媒装置3が浄化するNOxの他の成分、例えばHCの浄化率(これについては後述の他の実施形態で説明する)も概ね最大となるような空燃比である。そして、そのような空燃比の値(AF1,AF2,AF3等)は、基本的には触媒装置3の劣化が進行するに伴い、リーン側にシフトする。さらにそのような空燃比の値(AF1,AF2,AF3等)の付近におけるNOxの浄化率の変化は、基本的には、触媒装置3の劣化が進行するに伴い、急峻になる(浄化率のグラフの凸形状がより先鋭になる)。
【0074】
このような触媒装置3の浄化特性を考慮し、エンジン1の空燃比の変化に対する前記排ガスセンサ5(NOxセンサ)の出力VSの特性をさらに説明する。尚、以下の説明では、触媒装置3の浄化率や排ガスセンサ5の出力VS等の、エンジン1の空燃比に対する特性を単に対空燃比特性ということがある。
【0075】
触媒装置3が上述のようなNOxの浄化特性を有するため、該触媒装置3の下流側の排ガスセンサ5(NOxセンサ)が検出するNOx濃度、すなわち、触媒装置3による浄化後の排ガス中のNOx濃度は、図4に示す如く、エンジン1の空燃比の変化に対して、基本的には、前記図3のグラフの上下を反転させたような形態の特性を呈することとなる。つまり、該NOx濃度の対空燃比特性は、触媒装置3の新品状態(図4の実線のグラフ)を除いて、基本的には、極小値を持つ特性(図4の破線又は一点鎖線のような凹形状のグラフ)となる。
【0076】
そして、本実施形態では、NOx濃度に対する前記排ガスセンサ5の出力特性は前述したポジティブ特性であるため、排ガスセンサ5の出力VSの対空燃比特性も、図4に示したNOx濃度そのものの特性と同じになる。
【0077】
すなわち、その特性を改めて説明すれば、排ガスセンサ5の出力VSの対空燃比特性は、触媒装置3の新品状態を除いて(触媒装置3がある程度劣化した状態)、図4の破線あるいは一点鎖線のグラフで示す如く、触媒装置3によるNOxの浄化率が極大となるようなエンジン1の空燃比(AF2,AF3等)において排ガスセンサ5の出力VSが極小値を採る特性(凹形状のグラフ特性)となる。また、触媒装置3の新品状態では、図4の実線のグラフで示す如く、ある値AF1の空燃比よりもリッチ側の空燃比域、すなわち、NOxの浄化率がほぼ一定の最大値(略100%)に維持されるような空燃比域において、排ガスセンサ5の出力VSはほぼ一定の最小レベル(略0V)に維持される。そして、エンジン1の空燃比が上記値AF1よりもリーン側の空燃比になると、排ガスセンサ5の出力VSが急激に増加していく。
【0078】
また、排ガスセンサ5の出力VSが極小もしくは最小となるようなエンジン1の空燃比(AF1,AF2,AF3)は、触媒装置3の劣化の進行に伴い、リーン側に変移する。さらに、排ガスセンサ5の出力VSの極小値を含めて、エンジン1の空燃比の各値における排ガスセンサの出力VSは触媒装置3の劣化の進行に伴い、大きくなる。さらに、排ガスセンサ5の出力VSが極小もしくは最小となるようなエンジン1の空燃比(AF1,AF2,AF3)の値の前後の空燃比域における排ガスセンサ5の出力VSの対空燃比特性のグラフの傾きは、触媒装置3の劣化の進行に伴い大きくなる。
【0079】
以上説明した、触媒装置3の浄化特性や排ガスセンサ5(NOxセンサ)の出力VSの対空燃比特性を考慮しつつ、前記主演算処理制御器6及び燃料供給制御器7をさらに説明する。
【0080】
まず、前記主演算処理制御器6に関し、その処理の概要を説明すると、本実施形態では、目標空燃比KCMDの算出や触媒装置3の劣化状態の評価のために、前述のような排ガスセンサ5の出力VSの対空燃比特性を非線形関数としての二次関数により近似的に表現することとする(このことは排ガスセンサ5が検出するNOxの対空燃比特性を二次関数により表現することと実質的に等価である)。
【0081】
そして、主演算処理制御器6は、その同定手段6aとしての機能によって、上記二次関数のグラフ形状を規定するパラメータ、すなわち、該二次関数の各次数の項に係る係数や定数項の値を、該制御器6の制御サイクル毎に、エンジン1の運転時における前記LAFセンサ4の出力KACT(エンジン1の空燃比の検出値)及び排ガスセンサ5の出力VSのサンプリングデータを用いられる。
【0082】
さらに、主演算処理制御器6は、その目標空燃比算出手段6bとしての機能によって、上記パラメータの値を同定した二次関数を用いてエンジン1の目標空燃比KCMDを逐次求める。また、主演算処理制御器6は、その劣化状態評価手段6cとしての機能によって、上記二次関数のパラメータの同定値から、触媒装置3の劣化状態を評価するための複数種類(本実施形態では3種類)の劣化評価パラメータを求め、その劣化評価パラメータの値に基づいて触媒装置3の劣化状態の評価を行なう。
【0083】
さらに詳説すると、本実施形態で排ガスセンサ5の出力VSの対空燃比特性を近似的に表現する前記二次関数は、エンジン1の空燃比の検出値を表すLAFセンサ4の出力KACTを独立変数、排ガスセンサ5の出力VSを従属変数とする二次関数により与えられる。そして、本実施形態では、この二次関数を数式表現するために、独立変数としてのLAFセンサ4の出力KACT(空燃比の検出値)をそのまま用いる代わりに、該出力KACTとあらかじめ定めた所定の基準値FLAF/BASE(以下、空燃比基準値FLAF/BASEという)との偏差kact(=KACT−FLAF/BASE。以下、偏差空燃比kactという)を用い、次式(1)により、該二次関数を定義しておく。尚、上記空燃比基準値FLAF/BASEは、本実施形態では例えば理論空燃比の値としている。
【0084】
【数1】
【0085】
ここで、式(1)においては、この二次関数の従属変数の値である関数値(式(1)の右辺の演算結果の値)と、排ガスセンサ5の実際の出力VSの値とを区別するため、二次関数の従属変数としての排ガスセンサ5の出力の参照符号(変数記号)として、「VS」の代わりに「VSH1」を用いている。以下の説明では、この参照符号「VSH1」を付した排ガスセンサ5の出力を関数出力VSH1と称する。
【0086】
この二次関数の式(1)の二次の項「kact2」に係る係数a1と、一次の項「kact」に係る係数b1と、定数項c1とが、主演算処理制御器6の同定手段6aとしての機能によって、その値を同定するパラメータである。そして、主演算処理制御器6の同定手段6aが、上記パラメータa1,b1,c1の値を同定するアルゴリズムは、次のように構築された逐次型同定アルゴリズムである。
【0087】
すなわち、この逐次型同定アルゴリズムでは、主演算処理制御器6の制御サイクル毎に、LAFセンサ4の出力の現在値KACT(k)(kは制御サイクルの番数を示す。以下、同様)から前記空燃比基準値FLAF/BASEを減算してなる偏差空燃比kact(k)と、前記パラメータa1,b1,c1の同定値の現在値a1(k),b1(k),c1(k)(これは、基本的には前回の制御サイクルで求めたパラメータa1,b1,c1の同定値である)とを用いて次式(2)の演算を行うことで、各制御サイクルにおける排ガスセンサ5の関数出力VSH1(k)が求められる。
【0088】
【数2】
【0089】
尚、式(2)中の「Θ1」、「ξ1」は、同式(2)の但し書きで定義した通りのベクトルである。そして、式(2)やその但し書きで用いている「T」は転置を意味する(以下、同様)。
【0090】
さらに、次式(3)のように、上記関数出力VSH1(k)と、現在の制御サイクルにおける排ガスセンサ5の実際の出力VS(k)との間の偏差として与えらる同定誤差ID/E1(k)が求められる。
【0091】
【数3】
【0092】
主演算処理制御器6の同定手段6aが実行する逐次型同定アルゴリズムは、上記同定誤差I D/E1を最小にするように、前記パラメータa1,b1,c1の新たな同定値a1(k+1),b1(k+1),c1(k+1)、換言すれば、これらの同定値を成分とする新たな前記ベクトルΘ1(k+1)(以下、同定パラメータベクトルΘ1という)を求めるものである。そして、この新たな同定パラメータベクトルΘ1(k+1)の算出、すなわち、パラメータa1,b1,c1の同定値の更新は、次式(4)により行なわれる。
【0093】
【数4】
【0094】
すなわち、パラメータa1,b1,c1の現在の同定値a1(k),b1(k),c1(k)(前回の制御サイクルで求められた同定値)を、今回の制御サイクルにおける同定誤差ID/E1(k)に比例させた量だけ変化させることで、パラメータa1,b1,c1の新たな同定値a1(k+1),b1(k+1),c1(k+1)が求められる。
【0095】
ここで、式(4)中の「Kp1(k)」は、制御サイクル毎に次式(5)により決定される三次のベクトルで、各パラメータa1,b1,c1の同定値の同定誤差ID/E1に応じた変化度合い(ゲイン)を規定するものである。
【0096】
【数5】
【0097】
また、式(5)中の「P1(k)」は、制御サイクル毎に次式(6)の漸化式により更新される三次の正方行列である。
【0098】
【数6】
【0099】
尚、式(6)の行列P1(k)の初期値P1(0)は、その各対角成分を正の数とした対角行列である。また、式(6)中の「λ1」、「λ2」は0<λ1≦1及び0≦λ2<2の条件を満たすように設定される。
【0100】
この場合、上記λ1,λ2の設定の仕方によって、逐次型最小ニ乗法、逐次型重み付き最小二乗法、固定ゲイン法、漸減ゲイン法等、各種の具体的な同定アルゴリズムが定まる。いずれの手法を用いてもよいが、本実施形態では、例えば逐次型最小二乗法(この場合、λ1=λ2=1)の同定アルゴリズムを採用している。
【0101】
以上説明したアルゴリズムが、前記同定手段6aによって、二次関数のパラメータa1,b1,c1の値を同定するための逐次型同定アルゴリズムの内容である。
【0102】
尚、本実施形態では、パラメータa1,b1,c1の値の同定に際して、特に二次関数の二次の項の係数であるパラメータa1の同定値を制限する等の付加的な処理も行なうのであるが、これらについては後述する。
【0103】
次に、主演算処理制御器6の目標空燃比算出手段6bに関し、該目標空燃比算出手段6bは、上述のようにしてパラメータa1,b1,c1の値を同定した二次関数を用いて、次のように、エンジン1の目標空燃比KCMD(LAFセンサ4が検出する空燃比の目標値)を制御サイクル毎に求める。
【0104】
すなわち、上記のように二次関数のパラメータa1,b1,c1の値を同定したとき、その同定値により定まる二次関数のグラフ形状は、排ガスセンサ5の出力VSの対空燃比特性のグラフ形状を近似するものとなる。このとき、排ガスセンサ5の出力VSの対空燃比特性は、前述のように、基本的には極小値を有する特性(グラフ形状が凹形状となる特性)であるので、パラメータa1,b1,c1の同定値により定まる二次関数のグラフ形状も基本的には極小値をもつ凹形状となる。
【0105】
例えば、触媒装置3のある劣化状態において、排ガスセンサ5の出力VSと、LAFセンサ4による偏差空燃比kact(=KACT−FLAF/BASE)とのサンプリングデータは、図5に点描部分で示すような傾向で得られ、このとき、パラメータa1,b1,c1の同定値により定まる二次関数のグラフは、同図5に実線で示すような凹形状のものとなる。
【0106】
また、本実施形態では、基本的には、排ガスセンサ5の出力VSが極小となるようなエンジン1の空燃比(これは排ガスセンサ5の出力VSから把握されるNOx濃度自体が極小となるような空燃比である)を目標空燃比KCMDとして求める。
【0107】
そこで、主演算処理制御器6の目標空燃比算出手段6bは、基本的には、前記二次関数の関数値(前記関数出力VSH1)が極小値となるような偏差空燃比kactの値(図5の点Aにおける偏差空燃比kactの値)を、LAFセンサ4が検出する空燃比と前記空燃比基準値FLAF/BASEとの偏差の目標値、すなわち目標空燃比KCMDと空燃比基準値FLAF/BASEとの偏差(=KCMD−FLAF/BASE。以下、目標偏差空燃比k cmdという)として求める。
【0108】
この目標偏差空燃比kcmdは、該二次関数のパラメータa1,b1,c1のうちの、パラメータa1,b1の同定値(詳しくは最新の同定値a1(k+1),b1(k+1))を用いて次式(7)により求めることができる。
【0109】
【数7】
【0110】
そして、目標空燃比算出手段6bは、基本的には、この式(7)により得られる目標偏差空燃比kcmdに前記空燃比基準値FLAF/BASEを加算することで、制御サイクル毎の目標空燃比KCMDを決定する。
【0111】
ところで、前述したように、触媒装置3の新品状態では、排ガスセンサ5の出力VSの対空燃比特性は極小値をもたない。このような場合、排ガスセンサ5の出力VSと、LAFセンサ4による偏差空燃比kact(=KACT−FLAF/BASE)とのサンンプリングデータは例えば図6に点描部分で示すような傾向で得られ、このとき、パラメータa1,b1,c1の同定値により定まる二次関数のグラフは、同図6に実線で示すようなものとなる。
【0112】
このような場合、仮に、二次関数の関数値の極小値に対応する偏差空燃比kactの値を目標偏差空燃比kcmdとして求めると、この目標偏差空燃比kcmdに、空燃比基準値F LAF/BASEを加算してなる目標空燃比KCMDは、排ガスセンサ5の出力VSがほぼ一定に維持される空燃比域、すなわち、触媒装置3によるNOxの浄化率がほぼ一定の最大値(略100%)に維持される空燃比域内の値となる。
【0113】
このとき、このような値の目標空燃比KCMDは、NOxの浄化性能を確保する上では支障はないが、一般に触媒装置3による他のガス成分の浄化率が低下する。例えば、新品状態の触媒装置3によるHCの浄化率は、詳細は後述の他の実施形態で説明するが、前記図3に示した空燃比の値AF1にほぼ等しい値よりもリーン側の空燃比域ではほぼ一定の最大浄化率となる。そして、当該値よりもリッチ側の空燃比域では、浄化率が低下する。
【0114】
また、上記のように触媒装置3が新品状態であって、排ガスセンサ5の出力VSの対空燃比特性が極小値を持たないような場合に、パラメータa1,b1,c1の同定値により定まる二次関数においては、その極小値は、図6に示した如く、排ガスセンサ5の実際の出力VSが採り得ない負の値(NOx濃度が略「0」のときの排ガスセンサ5の出力VSよりも小さい値)となる。
【0115】
このようなことから、本実施形態では、主演算処理制御器6の目標空燃比算出手段6bは、パラメータa1,b1,c1の同定値により定まる二次関数の極小値が負の値となるような場合には、該二次関数の関数値(排ガスセンサ5の関数出力VSH1)が「0」となるような偏差空燃比kactの値(この値は二つある)のうち、リーン側の値(図6の点Bにおける偏差空燃比kactの値)を、目標偏差空燃比kcmdとして求める。
【0116】
この場合、二次関数の極小値が負の値となるのは、パラメータa1,b1,c1の同定値が、b12−4・a1・c1>0という不等式を満たす場合である。また、この場合の目標偏差空燃比kcmdは、各パラメータa1,b1,c1の同定値(詳しくは、最新の同定値a1(k+1),b1(k+1),c1(k+1))を用いて、次式(8)により求められる。
【0117】
【数8】
【0118】
そして、該二次関数の極小値が、「0」もしくは正の値となるような場合に、前述のように、該二次関数の極小値に対応する偏差空燃比kactの値を目標偏差空燃比kcmdとして求める。
【0119】
つまり、本実施形態では、目標空燃比算出手段6bは、制御サイクル毎に、その同定手段6aにより前述の如く求めた二次関数のパラメータa1,b1,c1の同定値a1(k+1),b1(k+1),c1(k+1)により定まる判別式(b12−4・a1・c1)の値に応じて、次式(9)により、目標偏差空燃比kcmd(k)を求める。
【0120】
【数9】
【0121】
そして、この式(9)に従って求めた目標偏差空燃比kcmd(k)に、次式(10)のうように空燃比基準値FLAF/BASEを加算することで、制御サイクル毎の目標空燃比KCDM(k)を決定する。
【0122】
【数10】
【0123】
以上説明した処理が、主演算処理制御器6の目標空燃比算出手段6の基本的処理である。
【0124】
次に、主演算処置制御器6の劣化状態評価手段6cは、二次関数のパラメータa1,b1,c1の同定値を用いて、次のように構築されたアルゴリズムによって触媒装置3の劣化状態を主演算処理制御器6の制御サイクル毎に評価する。
【0125】
排ガスセンサ5の出力VSあるいはその出力VSが表すNOx濃度の対空燃比特性は、触媒装置3の劣化の進行に伴い、前述のように変化するため(図4を参照)、その特性を近似する二次関数のグラフは、触媒装置3の劣化の進行に伴い、次のような特徴的な変化を呈する。図7、触媒装置3の劣化が比較的小さい場合、該劣化が比較的大きい場合にそれぞれ対応する二次関数のグラフp,qを用いて、触媒装置3の劣化の進行に伴う二次関数のグラフの変化の様子を概念的に示す。
【0126】
すなわち、図7参照して、触媒装置3の劣化が進行すると、前記パラメータa1,b1,c1の値を同定してなる二次関数のグラフは、矢印Xで表す如く、全体的に空燃比のリーン側に変移する。このため、例えば二次関数のグラフの極小点Cにおける偏差空燃比kactの値をAGDP1(図では、グラフp,qにそれぞれ対応させて添え字(p),(q)を付している)とおくと、この値AGDP1、あるいはこの値AGDP1に前記空燃比基準値FLAF/BASEを加算してなるエンジン1の空燃比の値は、触媒装置3の劣化の進行に伴い、リーン側に変移する(値AGDP1が小さくなる)。このことは、排ガスセンサ5が検出する実際のNOx濃度(触媒装置3による浄化後のNOx濃度)が極小となるようなエンジン1の空燃比が触媒装置3の劣化の進行に伴い、リーン側に変移することに対応するものである。
【0127】
尚、値AGDP1は、基本的には前述の如く、目標偏差空燃比kcmdとなる値である。
【0128】
また、二次関数のグラフは、触媒装置3の劣化が進行すると、矢印Yで表す如く、全体的に該二次関数の関数値(関数出力V SH1)が大きくなる方向(該関数値が表すNOx濃度が増加する方向)に変移する。このため、例えば二次関数のグラフの極小点Cにおける関数値、すなわち該関数値の極小値(関数出力VSH1の極小値)をAGDP2(図では、グラフp,qにそれぞれ対応させて添え字(p),(q)を付している)とおくと、この極小値AGDP2は、触媒装置3の劣化の進行に伴い、大きくなる。このことは排ガスセンサ5が検出する実際のNOx濃度の極小値が触媒装置3の劣化の進行に伴い大きくなることに対応するものである。
【0129】
また、二次関数のグラフの極小点Cの両側の傾き(詳しくは、極小点Cにおける偏差空燃比kactの値よりもリーン側及びリッチ側の両者の空燃比域におけるグラフの傾き)に着目すると、その傾きは、触媒装置3の劣化が進行すると、矢印Zで表す如く、大きくなる。換言すれば、二次関数のグラフの幅が触媒装置3の劣化の進行に伴い小さくなる。このため、例えば二次関数の極小点Cにおける偏差空燃比kactの値AGDP1から例えば、図示の如く所定の空燃比変化量αだけ、リッチ側に変化させた空燃比の値(AGDP1+α)における二次関数のグラフの傾き、すなわち図の点Dにおけるグラフの接線の傾きをAGDP3(図ではグラフp,qにそれぞれ対応させて添え字(p),(q)を付している)とおくと、該傾きAGDP3は、触媒装置3の劣化の進行に伴い、大きくなる。このことは、排ガスセンサ5が検出する実際のNOx濃度が極小となるエンジン1の空燃比の前後の空燃比域における該NOx濃度の対空燃比特性のグラフの傾きが、触媒装置3の劣化の進行に伴い、大きくなることに対応するものである。換言すれば、排ガスセンサ5が検出する実際のNOx濃度が極小値もしくはその近傍の値となるような空燃比域が狭まることに対応するものである。
【0130】
本実施形態における劣化状態評価手段6cのアルゴリズムは、上述にしたような触媒装置3の劣化の進行に伴う二次関数のグラフの変化に着目して構築されたものである。
【0131】
さらに詳細には、本実施形態では、二次関数のグラフの極小点Cにおける偏差空燃比kactの値AGDP1と、二次関数の関数値(関数出力VSH1)の極小値AGDP2と、二次関数の点Dにおけるグラフの傾きAGDP3とを、触媒装置3の劣化状態を評価するための三種類の劣化評価パラメータとして用いる。この場合、各劣化評価パラメータAGDP1,AGDP 2,AGDP3の値は、二次関数のパラメータa1,b1,c1の同定値を用いて、それぞれ次式(11)、(12)、(13)により求められる。
【0132】
【数11】
【0133】
【数12】
【0134】
【数13】
【0135】
このように、劣化評価パラメータAGDP1,AGDP2,AGDP3(以下、これらを順に第1、第2、第3の劣化評価パラメータという)を定めたとき、前述のように、第1の劣化評価パラメータAGDP1は、触媒装置3の劣化の進行に伴い、その値が小さくなる。また、第2及び第3の劣化評価パラメータAGDP2,AGDP3はいずれも、触媒装置3の劣化の進行に伴い、その値が大きくなる。そして、本実施形態では、触媒装置3の劣化状態の評価は、該触媒装置3が例えばその交換が必要か、もしくはその交換時期が近い程度に劣化した状態(以下、この劣化状態を劣化進行状態という)であるか否か(以下、劣化進行状態でない触媒装置3の状態を未劣化状態という)を判断する評価としている。
【0136】
そこで、本実施形態では、主演算処理制御器6の劣化状態評価手段6cは、主演算処理制御器6の制御サイクル毎に、二次関数のパラメータa1,b1,c1の同定値(詳しくは最新の同定値a1(k+1),b1(k+1) ,c1(k+1))から前記式(11)〜(13)によりそれぞれ求められる第1〜第3の劣化評価パラメータAGDP1,AGDP2,AGDP3の値を各劣化評価パラメータに対応させてあらかじめ定めた閾値LMT1,LMT2,LMT3と比較する。これにより、第1〜第3の各劣化評価パラメータAGDP1,AGDP2,AGDP3毎に、触媒装置3が劣化進行状態であるか未劣化状態であるかを暫定的に判断する。そして、それらの判断結果がいずれも「劣化進行状態」である場合に、触媒装置3の劣化状態を「劣化進行状態」であると決定する。
【0137】
これが、劣化状態評価手段6cによる触媒装置3の劣化状態の評価のアルゴリズムである。
【0138】
次に、前記燃料供給制御器7をさらに説明する。
【0139】
燃料供給制御器7は、図8に示すように、その機能的構成として、エンジン1の基本燃料噴射量Timを求める基本燃料噴射量算出部10と、基本燃料噴射量Timを補正するための第1補正係数KTOTAL及び第2補正係数KCMDMをそれぞれ求める第1補正係数算出部11及び第2補正係数算出部12とを具備する。
【0140】
前記基本燃料噴射量算出部10は、エンジン1の回転数NEと吸気圧PBとから、それらに応じたエンジン1の基準の燃料噴射量(燃料供給量)をあらかじめ設定されたマップを用いて求め、その基準の燃料噴射量をエンジン1の図示しないスロットル弁の有効開口面積に応じて補正することで基本燃料噴射量Timを算出するものである。
【0141】
また、第1補正係数算出部11が求める第1補正係数KTOTALは、エンジン1の排気還流率(エンジン1の吸入空気中に含まれる排ガスの割合)や、エンジン1の図示しないキャニスタのパージ時にエンジン1に供給される燃料のパージ量、エンジン1の冷却水温、吸気温等を考慮して前記基本燃料噴射量Timを補正するためのものである。
【0142】
また、第2補正係数算出部12が求める第2補正係数KCMDMは、前記目標空燃比算出制御器6が生成した目標空燃比KCMDに対応してエンジン1へ流入する燃料の冷却効果による吸入空気の充填効率を考慮して基本燃料噴射量Timを補正するためのものである。
【0143】
これらの第1補正係数KTOTAL及び第2補正係数KCMDMによる基本燃料噴射量Timの補正は、第1補正係数KTOTAL及び第2補正係数KCMDMを基本燃料噴射量Timに乗算することで行われ、この補正によりエンジン1の要求燃料噴射量Tcy lが得られる。
【0144】
尚、基本燃料噴射量Timや、第1補正係数KTOTAL、第2補正係数KCMDMのより具体的な算出手法は、特開平5−79374号公報等に本願出願人が開示しているので、ここでは詳細な説明を省略する。
【0145】
燃料供給制御器7は、上記の機能的構成の他、さらに、目標空燃比算出制御器6が逐次生成する目標空燃比KCMDにLAFセンサ4の出力KACT(エンジン1の空燃比の検出値)を収束させるように、フィードバック制御により3に対する燃料噴射量を調整するフィードバック制御部13を具備している。
【0146】
このフィードバック制御部13は、本実施形態では、エンジン1の各気筒の全体的な空燃比を制御する大局的フィードバック制御部14と、エンジン1の各気筒毎の空燃比を制御する局所的フィードバック制御部15とに分別される。
【0147】
前記大局的フィードバック制御部14は、LAFセンサ4の出力KACTを前記目標空燃比KCMDに収束させるように前記要求燃料噴射量Tcylを補正する(要求燃料噴射量Tcylに乗算する)フィードバック補正係数KFBを逐次求めるものである。
【0148】
この大局的フィードバック制御部14は、LAFセンサ4の出力KACTと目標空燃比KCMDとの偏差に応じて周知のPID制御を用いて前記フィードバック補正係数KFBとしてのフィードバック操作量KLAFを生成するPID制御器16と、LAFセンサ4の出力KACTと目標空燃比KCMDとからエンジン1の運転状態の変化や特性変化等を考慮して前記フィードバック補正係数KFBを規定するフィードバック操作量KSTRを適応的に求める適応制御器17(図ではSTRと称している)とをそれぞれ独立的に具備している。
【0149】
ここで、本実施形態では、前記PID制御器16が生成するフィードバック操作量KLAFは、LAFセンサ4の出力KACT(エンジン1の空燃比の検出値)が目標空燃比KCMDに一致している状態で「1」となり、該操作量KLAFをそのまま前記フィードバック補正係数KFBとして使用できるようになっている。一方、適応制御器17が生成するフィードバック操作量KSTRはLAFセンサ4の出力KACTが目標空燃比KCMDに一致する状態で「目標空燃比KCMD」となるものである。このため、該フィードバック操作量KSTRを除算処理部18で目標空燃比KCMDにより除算してなるフィードバック操作量kstr(=KSTR/KCMD)が前記フィードバック補正係数KFBとして使用できるようになっている。
【0150】
そして、大局的フィードバック制御部14は、PID制御器16により生成されるフィードバック操作量KLAFと、適応制御器17が生成するフィードバック操作量KSTRを目標空燃比KCMDにより除算してなるフィードバック操作量kstrとを切換部19で適宜、択一的に選択する。さらに、その選択したいずれか一方のフィードバック操作量KLAF又はkstrを前記フィードバック補正係数KFBとして使用し、該補正係数KFBを前記要求燃料噴射量Tcylに乗算することにより該要求燃料噴射量Tcylを補正する。尚、かかる大局的フィードバック制御部14(特に適応制御器17)については後にさらに詳細に説明する。
【0151】
前記局所的フィードバック制御部15は、LAFセンサ4の出力KACTからエンジン1の各気筒毎の実空燃比#nA/F(n=1,2,3,4)を推定するオブザーバ20と、このオブザーバ20により推定された各気筒毎の実空燃比#nA/Fから各気筒毎の空燃比のばらつきを解消するよう、PID制御を用いて各気筒毎の燃料噴射量のフィードバック補正係数#nKLAFをそれぞれ求める複数(気筒数個)のPID制御器21とを具備する。
【0152】
ここで、オブザーバ21は、それを簡単に説明すると、各気筒毎の実空燃比#nA/Fの推定を次のように行うものである。すなわち、エンジン1からLAFセンサ4の箇所(各気筒毎の排ガスの集合部)にかけての系を、エンジン1の各気筒毎の実空燃比#nA/FからLAFセンサ4が検出する空燃比を生成する系と考え、これを、LAFセンサ4の検出応答遅れ(例えば一次遅れ)や、LAFセンサ4が検出する空燃比に対するエンジン1の各気筒毎の空燃比の時間的寄与度を考慮してモデル化する。そして、そのモデルの基で、LAFセンサ4の出力KACTから、逆算的に各気筒毎の実空燃比#nA/Fを推定する。
【0153】
尚、この種のオブザーバ21は、本願出願人が例えば特開平7−83094号公報にて詳細に説明しているので、ここでは、さらなる説明を省略する。
【0154】
また、局所的フィードバック制御部15の各PID制御器21は、LAFセンサ4の出力KACTを、燃料供給制御器7における前回の制御サイクルで各PID制御器21により求められたフィードバック補正係数#nKLAFの全気筒についての平均値により除算してなる値を各気筒の空燃比の目標値とする。そして、その目標値とオブザーバ21により求められた各気筒毎の実空燃比#nA/Fの推定値との偏差が解消するように、今回の制御サイクルにおける各気筒毎のフィードバック補正係数#nKLAFを求める。
【0155】
さらに、局所的フィードバック制御部15は、前記要求燃料噴射量Tcylに大局的フィードバック制御部14のフィードバック補正係数KFBを乗算してなる値に、各気筒毎のフィードバック補正係数#nKLAFを乗算することで、各気筒の出力燃料噴射量#nTout(n=1,2,3,4)を求める。
【0156】
このようにして求められる各気筒の出力燃料噴射量#nToutは、燃料供給制御器7に備えた各気筒毎の付着補正部22により吸気管の壁面付着を考慮した補正が各気筒毎になされた後、各気筒に対する燃料噴射量の指令値として、エンジン1の図示しない燃料噴射装置に与えられる。そして、その付着補正がなされた出力燃料噴射量#nTout(燃料噴射量の指令値)に従って、気筒群4の各気筒への燃料噴射が行われるようになっている。
【0157】
尚、上記付着補正については、本願出願人が例えば特開平8−21273号公報にて詳細に説明しているので、ここではさらなる説明を省略する。
【0158】
前記大局的フィードバック制御部14、特に前記適応制御器17をさらに説明する。
【0159】
前記大局的フィードバック制御部14は、前述のようにLAFセンサ4の出力KACT(エンジン1の空燃比の検出値)を目標空燃比KCMDに収束させるようにフィードバック制御を行うものである。このとき、このようなフィードバック制御を周知のPID制御だけで行うようにすると、エンジン1の運転状態の変化や経年的特性変化等、動的な挙動変化に対して、安定した制御性を確保することが困難である。
【0160】
前記適応制御器17は、上記のようなエンジン1の動的な挙動変化を補償したフィードバック制御を可能とする漸化式形式の制御器であり、I.D.ランダウ等により提唱されているパラメータ調整則を用いて、図9に示すように、複数の適応パラメータを設定するパラメータ調整部23と、設定された適応パラメータを用いて前記フィードバック操作量KSTRを算出する操作量算出部24とにより構成されている。
【0161】
ここで、パラメータ調整部23について説明すると、ランダウ等の調整則では、離散系の制御対象の伝達関数B(Z-1)/A(Z-1)の分母分子の多項式を一般的に下記の式(14),(15)のようにおいたとき、パラメータ調整部23が設定する適応パラメータθハット(j)(jは制御サイクルの番数を示す)は、式(16)のようにベクトル(転置ベクトル)で表される。また、パラメータ調整部23への入力ζ(j)は、式(17)のように表される。この場合、本実施形態では、大局的フィードバック制御部14の制御対象であるエンジン1が一次系で3制御サイクル分の無駄時間dp(エンジン1の燃焼サイクルの3サイクル分の時間)を持つプラントと考え、式(14)〜式(17)でm=n=1,dp=3とし、設定する適応パラメータはs0,r1,r2,r3,b0の5個とした(図9参照)。尚、式(17)の上段式及び中段式におけるus,ysは、それぞれ、制御対象への入力(操作量)及び制御対象の出力(制御量)を一般的に表したものであるが、本実施形態では、上記入力をフィードバック操作量KSTR、制御対象(エンジン1)の出力を前記LAFセンサ4の出力KACT(空燃比の検出値)とし、パラメータ調整部23への入力ζ(j)を、式(17)の下段式により表す(図9参照)。
【0162】
【数14】
【0163】
【数15】
【0164】
【数16】
【0165】
【数17】
【0166】
ここで、前記式(16)に示される適応パラメータθハットは、適応制御器17のゲインを決定するスカラ量要素b0ハット-1(Z-1,j)、操作量を用いて表現される制御要素BRハット(Z-1,j)、及び制御量を用いて表現される制御要素S(Z-1,j)からなり、それぞれ、次式(18)〜(20)により表現される(図9の操作量算出部24のブロック図を参照)。
【0167】
【数18】
【0168】
【数19】
【0169】
【数20】
【0170】
パラメータ調整部23は、これらのスカラ量要素や制御要素の各係数を設定して、それを式(16)に示す適応パラメータθハットとして操作量算出部24に与えるもので、現在から過去に渡るフィードバック操作量KSTRの時系列データとLAFセンサ4の出力KACTとを用いて、該出力KACTが前記目標空燃比KCMDに一致するように、適応パラメータθハットを算出する。
【0171】
この場合、具体的には、適応パラメータθハットは、次式(21)により算出する。
【0172】
【数21】
【0173】
同式(21)において、Γ(j)は、適応パラメータθハットの設定速度を決定するゲイン行列(この行列の次数はm+n+dp)、eアスタリスク(j)は、適応パラメータθハットの推定誤差を示すもので、それぞれ次式(22),(23)のような漸化式で表される。
【0174】
【数22】
【0175】
【数23】
【0176】
ここで、式(23)中の「D(Z-1)」は、収束性を調整するための、漸近安定な多項式であり、本実施形態ではD(Z-1)=1としている。
【0177】
尚、式(22)のλ1(j),λ2(j)の選び方により、漸減ゲインアルゴリズム、可変ゲインアルゴリズム、固定トレースアルゴリズム、固定ゲインアルゴリズム等の種々の具体的なアルゴリズムが得られる。エンジン1の燃料噴射あるいは空燃比等の時変プラントでは、漸減ゲインアルゴリズム、可変ゲインアルゴリズム、固定ゲインアルゴリズム、および固定トレースアルゴリズムのいずれもが適している。
【0178】
前述のようにパラメータ調整部23により設定される適応パラメータθハット(s0,r1,r2,r3,b0)と、前記目標空燃比KCMDとを用いて、操作量算出部24は、次式(24)の漸化式により、フィードバック操作量KSTRを求める。図9の操作量算出部24は、同式(24)の演算をブロック図で表したものである。
【0179】
【数24】
【0180】
尚、式(24)により求められるフィードバック操作量KSTRは、LAFセンサ4の出力KACTが目標空燃比KCMDに一致する状態において、「目標空燃比KCMD」となる。このために、前述の如く、フィードバック操作量KSTRを除算処理部18によって目標空燃比KCMDで除算することで、前記フィードバック補正係数KFBとして使用できるフィードバック操作量kstrを求めるようにしている。
【0181】
このように構築された適応制御器17は、前述したことから明らかなように、制御対象であるエンジン1の動的な挙動変化を考慮した漸化式形式の制御器であり、換言すれば、エンジン1の動的な挙動変化を補償するために、漸化式形式で記述された制御器である。そして、より詳しくは、漸化式形式の適応パラメータ調整機構を備えた制御器と定義することができる。
【0182】
尚、この種の漸化式形式の制御器は、所謂、最適レギュレータを用いて構築する場合もあるが、この場合には、一般にはパラメータ調整機構は備えられておらず、エンジン1の動的な挙動変化を補償する上では、前述のように構成された適応制御器17が好適である。
【0183】
以上が適応制御器17の詳細である。
【0184】
尚、適応制御器17と共に、大局的フィードバック制御部14に具備したPID制御器16は、一般のPID制御と同様に、LAFセン4の出力KACTと、目標空燃比KCMDとの偏差から、比例項(P項)、積分項(I項)及び微分項(D項)を算出し、それらの各項の総和をフィードバック操作量KLAFとして算出する。この場合、本実施形態では、積分項(I項)の初期値を「1」とすることで、LAFセンサ4の出力KACTが目標空燃比KCMDに一致する状態において、フィードバック操作量KLAFが「1」になるようにし、該フィードバック操作量KLAFをそのまま燃料噴射量を補正するための前記フィードバック補正係数KFBとして使用することができるようしている。また、比例項、積分項及び微分項のゲインは、エンジン1の回転数と吸気圧とから、あらかじめ定められたマップを用いて決定される。
【0185】
また、大局的フィードバック制御部14の前記切換部19は、エンジン1の冷却水温の低温時や、高速回転運転時、吸気圧の低圧時等、エンジン1の燃焼が不安定なものとなりやすい場合、あるいは、目標空燃比KCMDの変化が大きい時や、空燃比のフィードバック制御の開始直後等、これに応じたLAFセンサ4の出力KACTが、そのLAFセンサ4の応答遅れ等によって、信頼性に欠ける場合、あるいは、エンジン1のアイドル運転時のようエンジン1の運転状態が極めて安定していて、適応制御器17による高ゲイン制御を必要としない場合には、PID制御器16により求められるフィードバック操作量KLAFを燃料噴射量を補正するためのフィードバック補正量数KFBとして出力する。そして、上記のような場合以外の状態で、適応制御器17により求められるフィードバック操作量KSTRを目標空燃比KCMDで除算してなるフィードバック操作量kstrを燃料噴射量を補正するためのフィードバック補正係数KFBとして出力する。これは、適応制御器17が、高ゲイン制御で、LAFセンサ4の出力KACTを急速に目標空燃比KCMDに収束させるように機能するため、上記のようにエンジン1の燃焼が不安定となったり、LAFセンサ4の出力KACTの信頼性に欠ける等の場合に、適応制御器17のフィードバック操作量KSTRを用いると、かえって空燃比の制御が不安定なものとなる虞れがあるからである。
【0186】
このような切換部19の作動は、例えば特開平8−105345号公報にて本願出願人が詳細に説明しているので、ここでは、さらなる説明を省略する。
【0187】
次に、本実施形態のシステムの全体の作動を詳細に説明する。
【0188】
まず、図10及び図11のフローチャートを参照して、前記燃料供給制御器7によるエンジン1の燃料噴射量の制御処理について説明する。燃料供給制御器7は、この処理をエンジン1のクランク角周期(TDC)と同期した制御サイクルで次のように行なう。
【0189】
燃料供給制御器7は、まず、エンジン1の回転数NE、吸気圧PB等を検出する図示しないセンサや、LAFセンサ4、排ガスセンサ5等、各種センサの出力を読み込む(STEPa)。
【0190】
次いで、基本燃料噴射量算出部10によって、前述の如くエンジン1の回転数NE及び吸気圧PBに対応する燃料噴射量をスロットル弁の有効開口面積に応じて補正してなる基本燃料噴射量Timが求められる(STEPb)。さらに、第1補正係数算出部11によって、エンジン1の冷却水温やキャニスタのパージ量等に応じた第1補正係数KTOTALが算出される(STEPc)。
【0191】
次いで、燃料供給制御器7は、前記主演算処理制御器6がその目標空燃比算出手段6bにより求める目標空燃比KCMDをエンジン1の空燃比を操作するために使用するか否か(ここでは空燃比操作のON/OFFという)の判別処理を行って、この空燃比操作のON/OFFを規定するフラグf/btc/onの値を設定する(STEPd)。このフラグf/ btc/onの値は、それが「0」であるとき、主演算処理制御器6が求める目標空燃比KCMDを使用しないこと(OFF)を意味し、「1」であるとき、主演算処理制御器6が求める目標空燃比KCMDを使用すること(ON)を意味する。
【0192】
上記の判別処理では、図11に示すように、前記排ガスセンサ5(NOxセンサ)が活性化しているか否かの判別(STEPd−1)、並びにLAFセンサ4が活性化しているか否かの判別(STEPd−2)が行なわれる。この判別は、それらのセンサの出力電圧等に基づいて行なわれる。
【0193】
このとき、排ガスセンサ5及びLAFセンサ4のいずれかが活性化していない場合には、それらセンサ5,4の出力データ(検出データ)を精度よく得ることができないため、フラグf/btc/onの値を「0」にセットする(STEPd−10)。また、エンジン1のリーン運転中(希薄燃焼運転)であるか否か、エンジン1の始動直後の触媒装置3の早期活性化を図るためにエンジン1の点火時期が遅角側に制御されているか否か、エンジン1のスロットル弁が略全開状態であるか否か、及びエンジン1のフュエルカット中(燃料供給の停止中)であるか否かの判別が行われる(STEPd−3〜d−6)。そして、これらのいずれかの条件が成立している場合には、主演算処理制御器6が求める目標空燃比KCMDを使用してエンジン1の空燃比を操作することは好ましくないか、もしくは操作することができないので、フラグf/btc/onの値を「0」にセットする(STEPd−10)。
【0194】
さらに、エンジン1の回転数NE及び吸気圧PBがそれぞれ所定範囲内(正常な範囲内)にあるか否かの判別が行なわれ(STEPd−7,d−8)、いずれかが所定範囲内に無い場合には、主演算処理制御器6が算出する目標空燃比KCMDを使用してエンジン1の空燃比を操作することは好ましくないので、フラグf/btc/onの値を「0」にセットする(STEPd−10)。
【0195】
そして、STEPd−1,d−2,d−7,d−8の条件が満たされ、且つ、STEPd−3〜d−6の条件が成立していない場合に(このような場合はエンジン1の通常的な運転状態である)、主演算処理制御器6が求める目標空燃比KCMDをエンジン1の空燃比の操作に使用すべく、フラグf/btc/onの値を「1」にセットする(STEPd−9)。
【0196】
図10の説明に戻って、上記のようにフラグf/btc/onの値を設定した後、燃料供給制御器7は、フラグf/btc/onの値を判断し(STEPe)、f/btc/on=1である場合には、主演算処理制御器6が求めた最新の目標空燃比KCMDを読み込む(STEPf)。また、f/btc/on=0である場合には、目標空燃比KCMDを所定値に設定する(STEPg)。この場合、目標空燃比KCMDとして設定する所定値は、例えばエンジン1の回転数NEや吸気圧PBからあらかじめ定めたマップ等を用いて決定される。
【0197】
次いで、燃料供給制御器7は、前記局所的フィードバック制御部15において、前述の如くオブザーバ20によりLAFセンサ4の出力KACTから推定した各気筒毎の実空燃比#nA/Fに基づき、PID制御器21により、各気筒毎の空燃比のばらつきを解消するためのフィードバック補正係数#nKLAFを算出する(STEPh)。さらに、大局的フィードバック制御部14により、フィードバック補正係数KFBを算出する(STEPi)。
【0198】
この場合、大局的フィードバック制御部14は、前述の如く、PID制御器16により求められるフィードバック操作量KLAFと、適応制御器17により求められるフィードバック操作量KSTRを目標空燃比KCMDで除算してなるフィードバック操作量kstrとから、切換部19によってエンジン1の運転状態等に応じていずれか一方のフィードバック操作量KLAF又はkstrを選択する(通常的には適応制御器17側のフィードバック操作量kstrが選択される)。そして、それを燃料噴射量を補正するためのフィードバック補正係数KFBとして決定する。
【0199】
尚、フィードバック補正係数KFBを、PID制御器16側のフィードバック操作量KLAFから適応制御器17側のフィードバック操作量kstrに切り換える際には、該補正係数KF Bの急変を回避するために、適応制御器17は、その切り換えの際の制御サイクルに限り、補正係数KFBを前回の制御サイクルにおける補正係数KFB(=KLAF)に保持するようにフィードバック操作量KSTRを求める。同様に、補正係数KFBを、適応制御器17側のフィードバック操作量kstrからPID制御器16側のフィードバック操作量KLAFに切り換える際には、PID制御器16は、自身が前回の制御サイクルで求めたフィードバック操作量KLAFが、前回の制御サイクルにおける補正係数KFB(=kstr)であったものとして、今回の制御サイクルにおけるフィードバック操作量KLAFを算出する。
【0200】
上記のようにしてフィードバック補正係数KFBを算出した後、燃料供給制御器7は、さらに、前記STEPf又はSTEPgで決定された目標空燃比KCMDに応じた前記第2補正係数KCMDMを第補正係数算出部12により算出する(STEPj)。
【0201】
次いで、燃料供給制御器7は、前述のように求められた基本燃料噴射量Timに、第1補正係数KTOTAL、第2補正係数KCMDM、フィードバック補正係数KFB、及び各気筒毎のフィードバック補正係数#nKLAFを乗算することで、各気筒毎の出力燃料噴射量#nToutを求める(STEPk)。そして、この各出力燃料噴射量#nToutが、付着補正部22によって、エンジン1の吸気管の壁面への燃料の付着を考慮した補正を施された後(STEPm)、各気筒毎の最終的な燃料噴射量の指令値として、エンジン1の図示しない燃料噴射装置に出力される(STEPn)。
【0202】
そして、エンジン1にあっては、各気筒毎の出力燃料噴射量#nToutに従って、各気筒への燃料噴射が行なわれる。
【0203】
以上のような各気筒毎の出力燃料噴射#nToutの算出及びそれに応じたエンジン1の燃料噴射がエンジン1のクランク角周期(TDC)に同期した制御サイクルで行なわれる。これにより、LAFセンサ4の出力KACT(空燃比の検出値)が、目標空燃比KCMDに収束するようにエンジン1の空燃比が操作される。この場合、特に、フィードバック補正係数KFBとして、適応制御器17側のフィードバック操作量kstrを使用している状態では、エンジン1の運転状態の変化や特性の変化等の挙動変化に対して、高い安定性を有してLAFセンサ4の出力KACTを迅速に目標空燃比KCMDに収束制御することができる。また、エンジン1が有する応答遅れの影響も適正に補償することができる。
【0204】
一方、前述のようなエンジン1の空燃比の操作(燃料噴射量の調整制御)と並行して、前記主演算処理制御器6は、一定周期の制御サイクルで図12のフローチャートに示すメインルーチン処理を実行する。
【0205】
すなわち、目標空燃比算出制御器6は、まず、自身の演算処理(前記パラメータa1,b1,c1の同定演算や目標空燃比KCMDの算出処理、触媒装置3の劣化状態の評価処理)を実行するか否かの判別処理を行なって、その実行の可否をそれぞれ値「1」、「0」で示すフラグf/btc/ calの値を設定する(STEP1)。また、このSTEP1では、パラメータa1,b1,c1の同定演算に関する後述の初期化を行なうか否かをそれぞれ値「1」、「0」で示すフラグf/id/resetの値も設定する。
【0206】
このSTEP1の処理は、図13のフローチャートに示すように行なわれる。
【0207】
すなわち、前記図10のSTEPdの場合と同様に、排ガスセンサ5及びLAFセンサ4が活性化しているか否かの判別が行なわれる(STEP1−1)。このとき、いずれかが活性化していない場合には、目標空燃比算出制御器6の演算処理に使用する排ガスセンサ5及びLAFセンサ4の検出データを精度よく取得することはできないため、フラグf/btc/calの値を「0」にセットする(STEP1−6)。
【0208】
さらにこのとき、パラメータa1,b1,c1の同定演算に関する後述の初期化を行なうために、フラグf/id/resetの値を「1」に設定する(STEP1−7)。
【0209】
また、エンジン1のリーン運転中(希薄燃焼運転)であるか否か(STEP1−3)、及びエンジン1の始動直後の触媒装置3の早期活性化を図るためにエンジン1の点火時期が遅角側に制御されているか否か(STEP1−4)の判別が行なわれる。これらのいずれかの条件が成立している場合には、主演算処理制御器6が算出する目標空燃比KCDMを前記燃料供給制御器7がエンジン1の燃料制御に使用することはないので、前記STEP1−6,1−7の処理が実行され、フラグf/btc/cal及びf/id/resetの値がそれぞれ「0」、「1」にセットされる。
【0210】
そして、STEP1−1,1−2の条件が満たされ、且つSTEP1−3,1−4の条件が成立していない場合には、フラグf/btc/calの値を「1」にセットする(STEP1−5)。
【0211】
尚、このようにフラグf/btc/calの値を設定することで、主演算処理制御器6が算出する目標空燃比KCMDを燃料供給制御器7が使用しない状況(図11を参照)であっても、例えばエンジン1のフュエルカット中やスロットル弁の全開時には、フラグf/btc/calの値が「1」に設定される。従って、エンジン1のフュエルカット中やスロットル弁の全開時には、主演算処理制御器6は目標空燃比KCMDの算出処理や触媒装置3の劣化状態の評価処理を行なうこととなる。これは、このようなエンジン1の運転状況は基本的には一時的なものであるからである。
【0212】
図12の説明に戻って、上記のような判別処理を行なった後、主演算処理制御器6は、次に、STEP1で設定したフラグf/btc/calの値を判断する(STEP2)。このとき、f/btc/cal=0である場合には、目標空燃比KCMDの算出や触媒装置3の劣化状態の評価等を行なうことなく、直ちに今回の制御サイクルの処理を終了する。
【0213】
一方、STEP2の判断で、f/btc/cal=1である場合には、LAFセンサ4の出力KACTの現在値KACT(k)のデータと排ガスセンサ5の出力VSの現在値VS(k)のデータとをそれぞれのセンサから取得する(STEP3)。
【0214】
さらに主演算処理制御器6は、STEP3で取得したLAFセンサ4の出力KACT(k)から、前記空燃比基準値F/LAF/BASEを減算することで、偏差出力kact(k)を求める(STEP4)。
【0215】
次いで、主演算処理制御器6は、その同定手段6aとしての機能によって、前記パラメータa1,b1,c1の同定演算に関する処理を実行する(STEP5)。
【0216】
この処理は、図14のフローチャートに示すように行なわれる。
【0217】
すなわち、主演算処理制御器6の同定手段6aは、まず、前記フラグf/id/resetの値を判断する(STEP5−1)。このとき、f/id/reset=である場合には、同定演算に関する初期化を行なう(STEP5−2)。この初期化では、前記パラメータa1,b1,c1の同定値の各値があらかじめ定めた初期値に設定される(前記同定パラメータベクトルΘ1が初期化される)。また前記式(6)の行列P1(対角行列)の各成分があらかじめ定めた初期値に設定される。さらに、フラグf/id/resetの値は「0」にリセットされる。
【0218】
次いで、主演算処理制御器6の同定手段6aは、パラメータa1 ,b1,c1の現在値a1(k),b1(k),c1(k)(これは基本的には、前回の制御サイクルで求められた同定値である)と、前記STEP4で求めたLAFセンサ4の偏差出力kactの現在値kact(k)とを用いて前記式(2)の演算を行うことで、排ガスセンサ5の関数出力VSH1(k)を求める(STEP5−3)。
【0219】
さらに同定手段6aは、パラメータa1,b1,c1の新たな同定値を求めるために使用する前記ベクトルKp1(k)と前記同定誤差ID/E1(k)とをそれぞれ算出する(STEP5−4,5−5)。
【0220】
この場合、同定誤差ID/E(k)は、前記STEP3で取得した排ガスセンサ5の現在の出力VS(k)とSTEP5−3で求めた関数出力VSH1(k)とから前記式(3)に従って算出される。また、ベクトルKp1(k)は、LAFセンサ4の現在の偏差出力kact(k)を用いて前記式(2)の但し書きに従って定義されるベクトルξ1(k)と、前記行列P1の現在値P1(k)(これは基本的には、前回の制御サイクルで最終的に決定される行列Pである)とから、前記式(5)の演算により算出される。
【0221】
次いで、同定手段6aは、エンジン1のフュエルカット中であるか否かを判断する(STEP5−6)。このとき、フュエルカット中である場合には、パラメータa1,b1,c1の値を同定する上で適正なLAFセンサ4の出力KACTや排ガスセンサ5の出力VSを得ることができないので、前記同定パラメータベクトルΘ1を現在値に保持する(パラメータa1,b1,c1の新たな同定値a1(k+1),b1(k+1),c1(k+1)を現在値a1(k),b1(k),c1(k)に維持する)と共に、前記行列P1を現在値に保持する(STEP5−7)。そして、STEP5の処理を終了し、図12の処理に復帰する。
【0222】
一方、STEP5−6でフュエルカット中でない場合には、前記STEP5−4,5−5で求めたベクトルKp1(k)及び同定誤差ID/E1(k)を用いて、前記式(4)に従って同定パラメータベクトルΘ1を更新することで、パラメータa1,b1,c1の新たな同定値a1(k+1),b1(k+1),c1(k+1)を求める(STEP5−8)。
【0223】
次いで、同定手段6aは、パラメータa1の新たな同定値a1(k+1)の値が正の値であるか否かを判断する(STEP5−9)。
【0224】
このとき、前記二次関数により表現する排ガスセンサ5の出力VSの対空燃比特性は、前述のように基本的には極小値をもつため、その二次関数の二次の項に係る係数であるパラメータa1の同定値a1(k+1)は、正の値でなければならない。また、外乱等の影響により、該同定値a1(k+1)が誤って負の値になると、この同定値a1(k+1)を二次の項の係数とした二次関数は極小値ではなく、極大値をもつものとなってしまう。そして、このような二次関数を用いても、触媒装置3によるNOxの浄化率を最大化(極大化)するための目標空燃比KCMDを適正に求めることはできない。
【0225】
このため、本実施形態では、STEP5−9の判断で、a1(k+1)≦0となっていた場合には、パラメータa1の最新の同定値a1(k+1)の値を、あらかじめ十分に小さい正の値(例えば0.01)に定めた所定値ε1に強制的に設定しなおす(STEP5−10)。これにより、パラメータa1の同定値は、二次関数の関数値である関数出力VSH1や該関数出力VSH1に対応するNOx濃度が極小値をもつような正の値に制限される。尚、STEP5−9の判断で、a1(k+1)>0であった場合には、パラメータa1の最新の同定値a1(k+1)はそのままに維持される。
【0226】
このようなパラメータa1の同定値a1(k+1)の制限処理を行なった後、同定手段6aは、次回の制御サイクルにおける演算処理のために、前記行列P1を前記式(6)に従って更新する(STEP5−11)。そして、STEP5の処理を終了し、図12の処理に復帰する。
【0227】
図12の説明に戻って、上述のように二次関数のパラメータa1,b1,c1の値の同定に関する演算処理を行なった後、主演算処理制御器6は、目標空燃比算出手段6bとしての機能によって、パラメータa1,b1,c1の最新の同定値a1(k+1),b1(k+ 1),c1(k+1)から、前記式(9)に従って、今回の制御サイクルにおける目標偏差空燃比kcmd(k)を算出する(STEP6)。
【0228】
すなわち、最新の同定値a1(k+1),b1(k+1),c1(k+1)により定まる判別式(b12−4・a1・c1)が「0」以下である場合(通常の場合)には、式(9)の上段側の式に従って目標偏差空燃比kcmd(k)を算出する。また、上記判別式(b12−4・a1・c1)の値が「0」よりも大きい正の値である場合(触媒装置3の新品状態に対応する場合)には、式(9)の下段側の式に従って目標偏差空燃比kcmd(k)を算出する。
【0229】
次いで、目標空燃比算出手段6bは、この目標偏差空燃比kcmd(k)に、前記空燃比基準値FLAF/BASEを加算することで、今回の制御サイクルにおける目標空燃比KCMD(k)を決定する(STEP7)。
【0230】
このようにして求められた目標空燃比KCMD(k)は図示しなメモリに時系列的に記憶保持される。そして、前記燃焼供給制御器7がその目標空燃比KCMD(k)を使用してエンジン1の空燃比を前述のように操作する場合にあっては、該燃料供給制御器7の制御サイクル毎に目標空燃比KCMD (k)の時系列データの中から、最新のものが選択される。
【0231】
尚、前述のように二次関数のパラメータa1の同定値a1(k+1)を強制的に正の値ε1に制限した場合(STEP5−10の処理を行った場合)において、そのパラメータa1の値とパラメータb1の同定値b1(k+1)とを用いて前記STEP6で算出される目標偏差空燃比kcmdは、本来あるべき目標偏差空燃比(触媒装置3の下流側のNOx濃度が極小となるような目標偏差空燃比)に対して誤差を生じている可能性がある。しかるに、前記同定手段8の逐次型同定アルゴリズムによれば、パラメータa1の同定値を正の値ε1に制限した後の制御サイクルにおいて、基本的には、速やかにパラメータa1の同定値が正の値となるようなパラメータa1,b1,c1の同定値の組を適正に得ることができる。また、二次関数のパラメータa1の同定値が負の値となるような状況は、基本的には、触媒装置3の新品状態で生じやすい。そして、このような触媒装置3の新品状態では、触媒装置3によるNOxの浄化率は、エンジン1の空燃比の広い範囲にわたって高いため、パラメータa1の同定値の強制的な制限の際に生じる目標偏差空燃比kcmdの誤差は、実用上、問題とならない。
【0232】
上記のように目標空燃比KCMDを求めた後、主演算処理制御器6は、その劣化状態評価手段6cとしての機能によって、触媒装置3の現在の劣化状態を評価する処理を実行する(STEP8)。
【0233】
この評価処理は、図15のフローチャートに示すように行なわれる。
【0234】
すなわち、主演算処理制御器6の劣化状態評価手段6cは、前記STEP5で求めた二次関数のパラメータの最新の同定値a1(k+1),b1(k+1),c1(k+1)から、前記式(11)、(12)、(13)によって、それぞれ前記第1〜第3の劣化評価パラメータAGDP1,AGDP2,AGDP3を算出する(STEP8−1)。
【0235】
次いで、劣化状態評価手段6cは、第1の劣化評価パラメータAGDP1の値を、この劣化評価パラメータに対応させてあらかじめ定めた閾値LMT1と比較する(STEP8−2)。このとき、前述の通り、第1の劣化評価パラメータAGDP1は、触媒装置3の劣化の進行に伴い、その値が小さくなる(図7参照)。従って、劣化状態評価手段6cは、AGDP1<LMT1である場合には、第1の劣化評価パラメータAGDP1に基づく触媒装置3の劣化状態の評価結果が「劣化進行状態」であるとして、フラグf/agd1の値を「1」にセットする(STEP8−3)。このフラグf/agd1は第1の劣化評価パラメータAGD1に基づく評価結果が「劣化進行状態」であるか「未劣化状態」であるかをそれぞれ値「1」、「0」で表すものであり、その初期値は「0」である。尚、STEP8−2でAGDP1≧LMT1である場合には、フラグf/agd1の値は現状に保持される。
【0236】
次いで、劣化状態評価手段6cは、第2の劣化評価パラメータAGDP2の値を、この劣化評価パラメータに対応させてあらかじめ定めた閾値L MT2と比較する(STEP8−4)。このとき、第2の劣化評価パラメータAGDP2は、触媒装置3の劣化の進行に伴い、その値が大きくなる(図7参照)。従って、劣化状態評価手段6cは、AGDP2>LMT2である場合には、第2の劣化評価パラメータAGDP2に基づく触媒装置3の劣化状態の評価結果が「劣化進行状態」であるとして、フラグf/agd2の値を「1」にセットする(STEP8−5)。このフラグf/agd2は第2の劣化評価パラメータAGD2に基づく評価結果が「劣化進行状態」であるか「未劣化状態」であるかをそれぞれ値「1」、「0」で表すものであり、その初期値は「0」である。尚、STEP8−4でAGDP2≦LMT2である場合には、フラグf/agd2の値は現状に保持される。
【0237】
さらに、劣化状態評価手段6cは、第3の劣化評価パラメータAGDP3の値を、この劣化評価パラメータに対応させてあらかじめ定めた閾値LMT3と比較する(STEP8−6)。このとき、第3の劣化評価パラメータAGDP3は、触媒装置3の劣化の進行に伴い、その値が大きくなる(図7参照)。従って、劣化状態評価手段6cは、AGDP3>LMT3である場合には、第3の劣化評価パラメータAGDP3に基づく触媒装置3の劣化状態の評価結果が「劣化進行状態」であるとして、フラグf/agd3の値を「1」にセットする(STEP8−7)。このフラグf/agd3は第3の劣化評価パラメータAGD3に基づく評価結果が「劣化進行状態」であるか「未劣化状態」であるかをそれぞれ値「1」、「0」で表すものであり、その初期値は「0」である。尚、STEP8−6でAGDP3≦LMT3である場合には、フラグf/agd3の値は現状に保持される。
【0238】
以上のようにて第1〜第3の各劣化評価パラメータAGD1,AGD 2,AGD3に基づく触媒装置3の劣化状態の評価を行なった後、劣化状態評価手段6cは、前記フラグf/agd1,f/agd2,f/agd3のいずれもが「1」であるか否かを判断する(STEP8−8)。このとき、フラグf/agd1=1且つf/agd2=1且つf/agd3=1である場合、すなわち、第1〜第3の各劣化評価パラメータAGD1,AGD2,AGD3に基づく評価結果(暫定的な評価結果)がいずれも「劣化進行状態」である場合には、劣化状態評価手段6cは、触媒装置3の現在の劣化状態が「劣化進行状態」であると確定して、前記劣化報知器8による報知を行なわしめる(STEP8−9)。すなわち、触媒装置3の交換を行なうべき時期である旨、もしくはその交換時期が近い旨の報知を劣化報知器8により行なう。そして、主演算処理制御器6は、図12の処理に復帰して今回の制御サイクルの処理を終了する。
【0239】
一方、STEP8−8で、フラグf/agd1,f/agd2,f/agd3のいずれかが「0」である場合には、主演算処理制御器6は、触媒装置3の現在の劣化状態が「未劣化状態」であると確定し、劣化報知器8による報知を行なうことなく、今回の制御サイクルの処理を終了する。
【0240】
以上説明した内容が本実施形態の装置の作動の詳細である。
【0241】
このような本実施形態の装置によれば、排ガスセンサ5(NOxセンサ)の出力VSの対空燃比特性(これは排ガスセンサ5が検出するNOx濃度の対空燃比特性でもある)を近似する二次関数のパラメータa1,b1,c1が排ガスセンサ5の実際の出力VSとLAFセンサ4の実際の出力KACT(エンジン1の実際の空燃比の検出値)とを用いて、逐次型同定アルゴリズム(本実施形態では逐次型最小二乗法のアルゴリズム)により主演算処理制御器6の制御サイクル毎に逐次同定される。
【0242】
このため、二次関数の各パラメータa1,b1,c1の値は、瞬時瞬時のエンジン1の運転状態や触媒装置3の挙動状態に則してリアルタイムで同定されることとなる。その結果、各パラメータa1,b1,c1の同定値により規定される二次関数は、上記の対空燃比特性を良好に近似し得るものとなる。
【0243】
さらに、このとき、二次関数の二次の項に係る係数であるパラメータa1の同定値は、排ガスセンサ5(NOxセンサ)の実際の出力VSの対空燃比特性が基本的に極小値をもつ特性であることを考慮し、二次関数が極小値をもつような正の値に制限される。このため、上記の対空燃比特性を近似するものとして、パラメータa1,b1,c1の同定値により規定される二次関数の信頼性を高めることができる。
【0244】
また、本実施形態では、パラメータa1,b1,c1の値を同定するために、LAFセンサ4の出力KACTをそのまま用いるのではなく、該出力KACTと所定の空燃比基準値FLAF/BASEとの偏差である前記偏差出力kactを用いている。このため、パラメータa1,b1,c1の同定値の精度を高めることができる。
【0245】
そして、本実施形態では、このようなパラメータa1,b1, c1の同定値により定まる劣化評価パラメータであって、触媒装置3の劣化の進行に伴いその値が増減する第1〜第3の劣化評価パラメータAGD1,AGD2,AGD3をそれぞれ、所定の閾値LMT1,LMT2,LMT3と比較することで、各劣化評価パラメータAGD1,AGD2,AGD3に基づく触媒装置3の劣化状態の評価(「劣化進行状態」であるか「未劣化状態」であるかの評価)を暫定的に行なう。さらにその暫定的な評価結果を総合することとで、触媒装置3の劣化状態の評価結果を確定する。
【0246】
これにより、触媒装置3の劣化状態の評価を、劣化状態評価手段6cの比較的簡単なアルゴリズムによって、高い信頼性で適正に行なうことができる。そして、その評価結果が「劣化進行状態」である場合に、劣化報知器8による報知を行なうことで、触媒装置3の交換を行なうべき時期である旨、あるいは触媒装置3の交換時期が近い旨を使用者に的確に認識させることができる。
【0247】
さらに、本実施形態では、目標空燃比算出手段6bの演算処理によって、基本的には、前述のような二次関数の関数値である排ガスセンサ5の関数出力VSH1が極小値となるような空燃比の値(これは該関数出力VSから把握されるNOx濃度が極小となるよう空燃比の値でもある)が目標空燃比KCMDとして設定されることとなる。
【0248】
このため、該目標空燃比KCMDは、排ガスセンサ5の実際の出力VSあるいはそれが表す実際のNOx濃度(触媒装置3による浄化後のNOx濃度)が極小となるようなエンジン1の空燃比の値にほぼ合致する。
【0249】
そして、本実施形態では、このような目標空燃比KCMDにLAFセンサ4の出力KACT、すなわち、エンジン1の実際の空燃比の検出値を収束させるようにフィードバック制御によってエンジン1の燃料噴射量(燃料供給量)を調整する。特に、このフィードバック制御は、前記適応制御器17を主体として行なわれる。これにより、エンジン1の空燃比を高い速応性で精度よく目標空燃比KCMDに制御することができる。
【0250】
この結果、触媒装置3によるNOxの最適な浄化、すなわち、NOxの浄化率が最大(極大)となるような浄化を行なうことができる。
【0251】
つまり、本実施形態では、触媒装置3によるNOxの最適な浄化を行ないながら、触媒装置3の劣化状態の評価を適正に行なうことができる。
【0252】
さらには、本実施形態では、触媒装置3の新品状態では、排ガスセンサ5の出力VSの対空燃比特性が極大値をもたないことを考慮する。そして、パラメータa1,b1,c1の同定値により規定される二次関数の関数値である排ガスセンサ5の関数出力VSH1が負の値になるような状況(前記の判別式(b12−4・a1・c1)の値が正の値となる状況)では、二次関数の関数値(関数出力VSH1)が「0」となるような空燃比の値のうちのリーン側の値が目標空燃比KCMDとして設定される。
【0253】
これにより、触媒装置3の新品状態でも、触媒装置3によるNOx以外のガス成分の浄化状態を損なったりすることなく、NOxを含めた種々のガス成分の触媒装置3による浄化を良好に行なうことができる。従って、触媒装置3の劣化状態によらずに、触媒装置3の良好な浄化性能を確保することができる。
【0254】
尚、本実施形態では、排ガスセンサ5は、その出力VSがNOx濃度に対して前述のポジティブ特性となるものを使用したが、図2に仮想線で示すようなネガティブ特性(NOxの増加に伴い、出力VSが小さくなる特性)のものを使用してもよい。この場合には、排ガスセンサが検出するNOx濃度(排ガスセンサ5の出力から把握されるNOx濃度)が、エンジン1の空燃比に対して前述のように極小値をもつことに対応して、排ガスセンサの出力自体は、エンジン1の空燃比に対して極大値をもつこととなる。そして、この場合であっても、本実施形態と同様に、二次関数のパラメータを用いて、触媒装置3の劣化状態の評価を行うことが可能である。
【0255】
この場合には、例えば、前記実施形態と同様に排ガスセンサの出力の対空燃比特性を二次関数により表現し、その二次関数のパラメータの値を同定する。このようにしたとき、該二次関数のグラフは、前述の実施形態と同様、触媒装置3の劣化の進行に伴い、空燃比のリーン側に変移し、また、その傾きが大きくなっていく。さらに、該二次関数のグラフは、触媒装置3の劣化の進行に伴い、その関数値が減少する方向(これは二次関数の関数値が表すNOx濃度が増加する方向である)に変移する。従って、例えば、その二次関数の関数値が極大値となる空燃比の値(偏差空燃比の値でもよい)と、該二次関数の関数値の極大値と、該二次関数の関数値が極大値となる空燃比の値からリッチあるいはリーン側に所定の空燃比変化量だけ変化させた空燃比の値における二次関数のグラフの傾き(より正確にはその絶対値)とを、それぞれ前記第1〜第3の劣化評価パラメータの代わりの評価パラメータとして用いて、触媒装置3の劣化状態の評価を行うようにすればよい。但し、この場合、二次関数の関数値の極大値は、前述の実施形態における第2の劣化評価パラメータ(二次関数の極小値)と逆に、触媒装置3の劣化の進行に伴い、値が小さくなる。このため、この二次関数の極大値に基づき触媒装置3の劣化状態を評価する場合、該極大値が所定の閾値よりも小さくなったときに、触媒装置3の劣化状態が「劣化進行状態」であると評価する必要がある。
【0256】
また、前記実施形態では、排ガスセンサ5の出力VSを二次関数の従属変数として用いているが、該排ガスセンサ5の出力VSから把握されるNOx濃度(図2の特性に従って出力VSから換算されるNOx濃度)を二次関数の従属変数として用いてもよい。この場合には、排ガスセンサ5の実際の出力VSのデータをNOx濃度に換算したデータを用いて二次関数のパラメータの値を同定する。そして、その二次関数のパラメータから前述の実施形態と全く同様に定めた劣化評価パラメータを用いて触媒装置3の劣化状態の評価を行うようにすればよい。尚、この場合、排ガスセンサ5はポジティブ特性及びネガティブ特性のいずれであっても触媒装置3の劣化状態を評価するためのアルゴリズムは同じになる。
【0257】
次に本発明の第2の実施形態を図16〜図19等を参照して説明する。尚、本実施形態のシステムは、前述の第1の実施形態のものと、排ガスセンサの種類(検出するガス成分の種類)と、主演算処理制御器の目標空燃比算出手段の処理の一部(詳しくは図12のSTEP6の処理)とが相違し、他のシステム構成及び作動は、第1の実施形態と同一である。従って、本実施形態の説明では、第1の実施形態と同一部分については、第1の実施形態と同一の参照符号を使用し、説明を省略する。
【0258】
本実施形態では、前記図1の触媒装置3の下流側に配置する排ガスセンサ5として、触媒装置3が浄化するガス成分のうちのHCの濃度を検出するHCセンサを用いている。この場合、排ガスセンサ5の出力VSのHC濃度に対する特性は、図16に実線で示すように、HC濃度の増加に伴い、ほぼリニアに出力VSのレベルが大きくなるポジティブ特性のものである。
【0259】
ここで、触媒装置3によるHCの浄化特性を説明しておく。
【0260】
触媒装置3は、エンジン1の空燃比の変化に対して基本的には図17に示すようなHCの浄化特性(HCの浄化率の特性)を有する。この場合、図17の実線、破線、及び一点鎖線のグラフは、それぞれ、前記図2の実線、破線、及び一点鎖線のグラフに対応する劣化状態の触媒装置3に関するものである。
【0261】
図17に実線のグラフで示すように、新品の触媒装置3によるHCの浄化率は、エンジン1の空燃比がある値AF4(この値は、図2に示した空燃比の値AF1とほぼ一致する)よりもリーン側の空燃比域にあるときには、ほぼ一定(略100%)の最大の浄化率となる。そして、エンジン1の空燃比が上記の値AF4よりもリッチ側の空燃比域になると、HCの浄化率が低下する。
【0262】
また、触媒装置3の劣化がある程度進行すると、図17の破線や一点鎖線のグラフで示すように、HCの浄化率は、ある値AF5(:破線),AF6(:一点鎖線)の空燃比において極大値を採る(ピークを持つ)。そして、エンジン1の空燃比がHCの浄化率の極大値に対応する値AF5,AF6から、リッチ側及びリーン側のいずれに変化しても、HCの浄化率が減少していく。尚、空燃比の上記値AF5,AF6は、それぞれ、図2に示した空燃比の値AF2,AF3と概ね一致する。
【0263】
この場合、HCの浄化率の極大値からの減少は、リッチ側よりもリーン側の方うが若干、緩やかに生じる(グラフの傾きがリッチ側よりもリーン側の方が若干小さい)。尚、触媒装置3が劣化した状態におけるHCの浄化率の極大値は、触媒装置3の新品状態における浄化率の最大値よりも小さい。また、該極大値は、触媒装置3の劣化の進行に伴い小さくなる。
【0264】
このように触媒装置3によるHCの浄化率は、NOxの浄化率と同様、触媒装置3の新品状態を除いて、エンジン1の空燃比の変化に対して基本的には極大値を持つような特性を呈する。
【0265】
尚、触媒装置3の新品状態を含めた該触媒装置3の各劣化状態において、HCの浄化率が最大(極大)となるような空燃比の値(前記AF4,AF5,AF6等)は、基本的には触媒装置3の劣化が進行するに伴い、リーン側にシフトする。さらにそのような空燃比の値(AF4,AF5,AF6等)の付近におけるHCの浄化率の変化は、基本的には、触媒装置3の劣化が進行するに伴い、急峻になる(浄化率のグラフの凸形状がより先鋭になる)。
【0266】
また、触媒装置3が上述のようなHCの浄化特性を有するため、該触媒装置3の下流側の排ガスセンサ5(HCセンサ)が検出するHC濃度、すなわち、触媒装置3による浄化後の排ガス中のHC濃度は、図18に示す如く、エンジン1の空燃比の変化に対して、基本的には、前記図17のグラフの上下を反転させたような形態の特性を呈することとなる。つまり、該HC濃度は、触媒装置3による浄化後のNOxと同様、触媒装置3の新品状態(図18の実線のグラフ)を除いて、基本的には、エンジン1の空燃比の変化に対して極小値を持つ特性(図18の破線又は一点鎖線のグラフ)となる。
【0267】
そして、本実施形態では、HC濃度に対する前記排ガスセンサ5の出力特性はポジティブ特性であるため、エンジン1の空燃比の変化に対する排ガスセンサ5の出力VSの特性も、図18に示したNOx濃度そのものの特性と同じになる。すなわち、その特性を改めて説明すれば、エンジン1の空燃比の変化に対する排ガスセンサ5(HCセンサ)の出力VSの特性は、触媒装置3の新品状態を除いて(ある程度劣化した状態では)、図18の破線あるいは一点鎖線のグラフで示す如く、触媒装置3によるHCの浄化率が極大となるようなエンジン1の空燃比(AF5,AF6等)において排ガスセンサ5の出力VSが極小値を採る特性(凹形状のグラフ特性)となる。また、触媒装置3の新品状態では、図18の実線のグラフで示す如く、ある値AF4の空燃比よりもリーン側の空燃比域、すなわち、HCの浄化率がほぼ一定の最大値(略100%)に維持されるような空燃比域において、排ガスセンサ5の出力VSはほぼ一定の最小レベル(略0V)に維持される。そして、エンジン1の空燃比が上記値AF4よりもリッチ側の空燃比になると、排ガスセンサ5の出力VSが増加していく。
【0268】
尚、排ガスセンサ5(HCセンサ)の出力VSが極小もしくは最小となるようなエンジン1の空燃比(AF4,AF5,AF6)は、排ガスセンサとしてNOxセンサを用いた場合と同様に、触媒装置3の劣化の進行に伴い、リーン側に変移する。さらに、排ガスセンサ5の出力VSの極小値を含めて、エンジン1の空燃比の各値における排ガスセンサの出力VSは、排ガスセンサとしてNOxセンサを用いた場合と同様に、触媒装置3の劣化の進行に伴い、大きくなる。さらに、排ガスセンサ5の出力VSが極小もしくは最小となるようなエンジン1の空燃比(AF1,AF2,AF3)の値の前後の空燃比域における排ガスセンサ5の出力VSの対空燃比特性のグラフの傾きは、排ガスセンサとしてNOxセンサを用いた場合と同様に、触媒装置3の劣化の進行に伴い大きくなる。つまり、排ガスセンサ5(HCセンサ)の出力VSの対空燃比特性のグラフの、触媒装置3の劣化の進行に伴う変化の形態は、排ガスセンサとしてNOxセンサを用いた場合と基本的には同じである。
【0269】
このような触媒装置3の浄化特性を考慮しつつ、本実施形態における主演算処理制御器6の処理、特に目標空燃比算出手段6bの処理についてさらに説明する。
【0270】
上述のように、排ガスセンサ5としてのHCセンサの出力VSの、エンジン1の空燃比に対する特性(対空燃比特性)は、基本的には、前記第1の実施形態で説明したNOxセンサの出力の対空燃比特性と同様に、極小値をもつ(グラフ形状が凹形状となる)。
【0271】
従って、本実施形態における排ガスセンサ5(HCセンサ)の出力VSの対空燃比特性(これは排ガスセンサ5の出力VSから把握されるHC濃度の対空燃比特性でもある)は、前記第1の実施形態と同様に前記式(1)の形の二次関数のパラメータa1,b1,c1の値を同定することで、該二次関数により近似的に表現することができる。そして、触媒装置3によるHCの浄化率を最大(極大)とするようなエンジン1の目標空燃比KCMDは、基本的には、上記二次関数の関数値(排ガスセンサ5の関数出力VSH1)が極小値となるような空燃比の値として与えられる。
【0272】
また、排ガスセンサ5(HCセンサ)の出力VSの対空燃比特性のグラフは、触媒装置3の劣化の進行に対して、第1の実施形態における排ガスセンサとしてのNOxセンサの出力の対空燃比特性と同様の変化を呈するので、上記二次関数も第1の実施形態の場合と同様の変化を呈することとなる。つまり、該二次関数のグラフは、触媒装置3の劣化の進行に伴い、空燃比のリーン側に変移し、また、該二次関数の関数値もしくはその関数値が表すHC濃度が増加する方向に変移する。さらに、該二次関数のグラフの極小点における空燃比の値の前後の空燃比域における該グラフの傾きは、触媒装置3の劣化の進行に伴い、大きくなる。
【0273】
そこで、本実施形態では、エンジン1の目標空燃比KCMDを求め、また触媒装置3の劣化状態を評価するために、排ガスセンサ5(HCセンサ)の出力VSの対空燃比特性を表現する非線形関数として、前記第1の実施形態と同様に、前記式(1)の形の二次関数を用いる。そして、その二次関数のパラメータa1,b1,c1の値を主演算処理制御器6の同定手段6aが第1の実施形態で説明した通りのアルゴリズム(逐次型同定アルゴリズム)によって、同定する(前記式(2)〜式(6)、並びに図14を参照)。
【0274】
そして、本実施形態においても、主演算処理制御器6の目標空燃比算出手段6bは、基本的には、パラメータa1,b1,c1の値を同定した二次関数の関数値が極小値となるような空燃比の値を目標空燃比KCMDとして求める。この場合、この目標空燃比KCMDを規定する目標偏差空燃比kcmd(=KCMD−FLAF/BASE)を、パラメータa1,b1の同定値を用いて前記式(7)により求めることができ、この目標偏差空燃比kcmdに空燃比基準値FLAF/BASEを加算すれば、目標空燃比KCMDが得られる。
【0275】
但し、本実施形態では、触媒装置3の新品状態では、前記第1の実施形態と異なる演算によって、目標空燃比KCMDを求めることとする。
【0276】
すなわち、触媒装置3の新品状態では、排ガスセンサ5(HCセンサ)の対空燃比特性は、図18に実線のグラフで示したようなものとなり、極小値をもたないと共に、リーン寄りの空燃比において出力VSがほぼ一定(略「0」)となる。そして、このような場合、例えば、図19に点描部分で示すような傾向で、排ガスセンサ5の出力V Sと、LAFセンサ4による偏差空燃比kact(=KACT−FLAF/BASE)とのサンンプリングデータが得られ、このとき、パラメータa1,b1,c1の同定値により定まる二次関数のグラフは、同図19に実線で示すようなものとなる。このとき、該二次関数の関数値(関数出力VSH1)の極小値は、図19に示した如く、排ガスセンサ5の実際の出力VSが採り得ない負の値(HC濃度が略「0」のときの排ガスセンサ5の出力VSよりも小さい値)となる。
【0277】
また、触媒装置3の新品状態で、HCを含めて種々のガス成分(NOx等)の浄化率が高いものとなる空燃比は、二次関数の関数値が極小値となるような空燃比の値ではなく、概ね図19の点Eにおける空燃比の値の近傍の空燃比である。
【0278】
このため、本実施形態では、目標空燃比算出手段6bは、パラメータa1,b1,c1の同定値(最新の同定値)により定まる二次関数の極小値が負の値となるような場合、すなわち、パラメータa1,b1,c1の同定値が、b12−4・a1・c1>0という不等式を満たす場合には、該二次関数の関数値(排ガスセンサ5の関数出力VSH1)が「0」となるような偏差空燃比kactの値(この値は二つある)のうち、リッチ側の値(図19の点Eにおける偏差空燃比kactの値)を、目標偏差空燃比kcmdとして求める。そして、その目標偏差空燃比kcmdに空燃比基準値FLAF/BASEを加算することで、目標空燃比KCMDを求める。
【0279】
つまり、本実施形態では、主演算処理制御器6の目標空燃比算出手段6bは、二次関数のパラメータa1,b1,c1の同定値(詳しくは最新の同定値a1(k+1),b1(k+1),c1(k+1))により定まる判別式(b12−4・a1・c1)の値に応じて、次式(25)により、目標偏差空燃比kcmd( k)を求める。
【0280】
【数25】
【0281】
そして、この式(25)に従って求めた目標偏差空燃比kcmd(k)に、前記式(10)のように空燃比基準値FLAF/BASEを加算することで、制御サイクル毎の目標空燃比KCDM(k)を決定する。
【0282】
この場合、前記式(25)では、b12−4・a1・c1≦0である場合は、前記第1の実施形態に関する式(9)と同じであるが、b12−4・a1・c1>0である場合(触媒装置3が新品状態である場合)は、式(9)と異なるものとなる。
【0283】
また、触媒装置3の劣化状態の評価に関しては、触媒装置3の劣化の進行に伴う二次関数のグラフの変化の形態は、前述のように第1の実施形態の場合と同様である。そこで、本実施形態においても、主演算処理制御器6の劣化状態評価手段6cは、二次関数のパラメータa1,b1,c1の同定値から前記式(11)〜(13)によりそれぞれ求められる第1〜第3の劣化評価パラメータAGDP1,AGDP2,AGDP3を用いる。そして、これらの第1〜第3の劣化評価パラメータAGDP1,AGDP2,AGDP3を前記図15のフローチャートに示した如く、あらかじめ定めた閾値LMT1,LMT2,LMT3と比較することで、触媒装置3の劣化状態が「劣化進行状態」であるか「未劣化状態」であるかを評価し、「劣化進行状態」である場合に、劣化報知器8による報知を行う。但し、この場合において、触媒装置3の劣化の進行に伴う二次関数の変化の度合いは、第1の実施形態と必ずしも同じではないので、第1〜第3の各劣化評価パラメータAGDP1,AGDP2,AGDP3と比較する閾値LMT1,LMT2,LMT3は、本実施形態に専用的なものである。
【0284】
上述のように、本実施形態における主演算処理制御器6の処理は、目標偏差空燃比kcmdの算出の演算処理、すなわち、前記図12のSTEP6における演算処理のみが第1の実施形態の実施形態と相違するものであり、その他の処理は、第1の実施形態と全く同一である。
【0285】
また、上述した以外の構成及び作動、例えば燃料供給制御器7の構成及び作動も第1の実施形態と全く同一である。
【0286】
かかる本実施形態の装置によれば、前記第1の実施形態で説明したものと同様の効果を奏することができる。すなわち、それを要約すれば、前記パラメータa1,b1,c1の同定によって、排ガスセンサ5(HCセンサ)の出力VSの対空燃比特性を良好に近似する二次関数を得ることができる。そして、この二次関数のパラメータa1,b1,c1から第1〜第3の三種類の劣化評価パラメータAGDP1,AGDP2,AGDP3を求め、これらの劣化評価パラメータAGDP1,AGDP2,AGDP3に基づいて、触媒装置3の劣化状態の評価を行うことによって、その評価を比較的簡単なアルゴリズムにより適正に行うことができる。
【0287】
また、上記二次関数を用いることで、排ガスセンサ5(HCセンサ)が検出するHC濃度、すなわち、触媒装置3による浄化後のHC濃度が極小となるようなエンジン1の目標空燃比KCMDを求めることができる。そして、この目標空燃比KCMDにエンジン1の空燃比を適応制御器17を用いたフィードバック制御によって、操作することで、触媒装置3によるHCの良好な浄化、すなわち、HCの浄化率が極大となるような浄化を行なうことができる。
【0288】
さらに、パラメータa1,b1,c1の同定値により規定される二次関数の関数値である排ガスセンサ5の関数出力VSH1が負の値になるような触媒装置3の新品状態では、二次関数の関数値(関数出力VSH1)が「0」となるような空燃比の値のうちのリーン側の値が目標空燃比KCMDとして設定されるので、触媒装置3の劣化状態によらずに、HCを含めた種々のガス成分の触媒装置3による浄化を良好に行なうことができる。
【0289】
尚、本実施形態にあっても、前記第1の実施形態に関連して説明した変形態様が可能である。すなわち、排ガスセンサ5は、図16に仮想線で示すように、HC濃度の増加に伴い出力VSが小さくなるようなネガティブ特性のものを用いることも可能である。さらに、二次関数の独立変数は、排ガスセンサ5の出力VSに代えて、該出力VSから図16の特性に従って把握されるHC濃度を用いてもよい。
【0290】
また、以上説明した第1及び第2の実施形態では、触媒装置3の劣化状態の評価のために、第1〜第3の三種類の劣化評価パラメータAGDP1,AGDP2,AGDP3を用いたが、それらのうちの一つもしくは二つの劣化評価パラメータを用いて触媒装置3の劣化状態の評価を行うようにしてもよい。
【0291】
また、前記各実施形態では、各劣化評価パラメータAGDP1,AGDP2,AGDP3に基づく暫定的な評価結果が、いずれも「劣化進行状態」である場合に限り、触媒装置3の劣化状態を「劣化進行状態」として確定するようにした。但し、例えば各劣化評価パラメータAG DP1,AGDP2,AGDP3に基づく暫定的な評価結果のうちのいずれか一つ、もしくは二つが「劣化進行状態」である場合に、触媒装置3の劣化状態を「劣化進行状態」として確定するようにしてもよい。
【0292】
また、前記各実施形態では、触媒装置3の劣化状態を「劣化進行状態」と「未劣化状態」との二つの劣化度合いに分けて評価するようにしたが、例えば各劣化評価パラメータAGDP1,AGDP2,AGDP3と比較する閾値をさらに多くすれば、触媒装置3の劣化状態を三段階以上の劣化度合いに分類して評価するようにすることも可能である。そして、この場合、触媒装置3の劣化度合いに応じて異なる報知を行うようにしてもよい。
【0293】
また、前記各実施形態では、触媒装置3の劣化の進行に伴う二次関数のグラフの傾きの変化を、前記第3の劣化評価パラメータAGDP3、すなわち、二次関数のグラフの極小点における偏差空燃比の値からリッチ側に所定の空燃比変化量αだけ変化させた偏差空燃比における二次関数のグラフの傾きによって把握するようにした。但し、二次関数のグラフの極小点における偏差空燃比の値からリーン側に所定の空燃比変化量αだけ変化させた偏差空燃比における二次関数のグラフの傾きによって、触媒装置3の劣化の進行に伴う二次関数のグラフの傾きの変化を把握するようにしてもよい。
【0294】
さらには、触媒装置3の劣化の進行に伴い、二次関数のグラフの傾きが大きくなると、該グラフの幅が狭くなることから、例えば、二次関数の関数値の極小値よりも所定量(一定量)だけ大きな関数値となる偏差空燃比の値(この値は二つある)を二次関数のパラメータa1,b1,c1の同定値を用いて求める。そして、それらの偏差空燃比の値の差(これは二次関数のグラフの幅に相当するものとなる)に基づいて、二次関数のグラフの傾きの変化を把握するようにすることも可能である。
【0295】
また、前記各実施形態では、排ガスセンサ5としてNOxセンサあるいはHCセンサを用いたが、触媒装置3により浄化する他のガス成分、例えば、CO(一酸化炭素)の濃度を検出するCOセンサ等を用いてもよい。このような排ガスセンサを用いた場合であっても、前記各実施形態と同様に、その排ガスセンサの出力の対空燃比特性を二次関数等の非線形関数により近似的に表現し、その非線形関数のパラメータを同定するようにすれば、触媒装置3の劣化状態を評価するようにすることも可能である。さらにはその非線形関数を用いて、排ガスセンサが検出するガス成分の濃度が極小となるようにエンジン1の空燃比を制御するようにすることも可能である。
【0296】
また、前記各実施形態では、排ガスセンサの出力あるいはその出力から把握されるガス成分の濃度の対空燃比特性を表現する非線形関数として、二次関数を用いた。但し、例えば三次関数や双極線関数等の他の形式の非線形関数を用いて、触媒装置3の劣化状態を評価し、さらには、排ガスセンサが検出するガス成分の濃度が極小となるようにエンジン1の空燃比を制御することも可能である。
【0297】
また、前述の各実施形態では、二次関数のパラメータa1,b1,c1の値を同定するために、LAFセンサ4の出力KACTのサンプリングデータを用いたが、LAFセンサ4の出力KACTは基本的には、目標空燃比KCMDに一致するように制御されるので、LAFセンサ4の出力KACTや偏差出力kactの代わりに、目標空燃比KCMDや目標偏差空燃比kcmdを用いて二次関数のパラメータa1,b1,c1の値を同定するようにすることも可能である。この場合には、例えば前記式(1)の「kact」の代わりに、目標空燃比算出制御器6が過去に求めた目標偏差空燃比k cmd(詳しくは、目標空燃比KCMDからLAFセンサ4が検出する空燃比を生成する系としてのエンジン1及び燃料供給制御器7からなる系が有する無駄時間前の制御サイクルで目標空燃比算出制御器6が求めた目標偏差空燃比kcmd)用いて排ガスセンサ5の関数出力VSH1を求める。そして、前述の各実施形態と同様の逐次型同定アルゴリズムによって、その関数出力VSH1が最小となるように二次関数のパラメータa1,b1,c1の値を同定すればよい。
【図面の簡単な説明】
【図1】本発明の第1の実施形態の装置の全体的システム構成を示すブロック図。
【図2】図1の装置で使用する排ガスセンサ(NOxセンサ)のNOx濃度に対する出力特性を示す線図。
【図3】図1の装置の触媒装置によるNOxの浄化特性を示す線図。
【図4】図1の装置で使用する排ガスセンサ(NOxセンサ)の出力あるいはその排ガスセンサが検出するNOx濃度の空燃比に対する特性を示す線図。
【図5】図1の装置の主演算処理制御器目標空燃比算出手段の処理を説明するための線図。
【図6】図1の装置の主演算処理制御器の目標空燃比算出手段の処理を説明するための線図。
【図7】図1の装置の主演算処理制御器の劣化状態評価手段の処理を説明するための線図。
【図8】図1の装置の燃料供給制御器の基本構成を示すブロック図。
【図9】図8の燃料供給制御器が備える適応制御器の基本構成を示すブロック図。
【図10】図1の装置の燃料供給制御器の処理を説明するためのフローチャート。
【図11】図10のフローチャートのサブルーチン処理を示すフローチャート。
【図12】図1の装置の主演算処理制御器の処理を説明するためのフローチャート。
【図13】図12のフローチャートのサブルーチン処理を示すフローチャート。
【図14】図12のフローチャートのサブルーチン処理を示すフローチャート。
【図15】図12のフローチャートのサブルーチン処理を示すフローチャート。
【図16】本発明の第2の実施形態で使用する排ガスセンサ(HCセンサ)のHC濃度に対する出力特性を示す線図。
【図17】第2の実施形態における触媒装置(図1の触媒装置)によるHCの浄化特性を示す線図。
【図18】第2の実施形態で使用する排ガスセンサ(HCセンサ)の出力あるいはその排ガスセンサが検出するHC濃度の空燃比に対する特性を示す線図。
【図19】第2の実施形態における目標空燃比の算出処理を説明するための線図。
【符号の説明】
1…エンジン、2…排気管(排気通路)、3…触媒装置、4…LAFセンサ(空燃比センサ)、5…排ガスセンサ。
Claims (15)
- 内燃機関の排気通路に設けた排ガス浄化用の触媒装置の下流側に、前記触媒装置を通過した排ガスのうちの該触媒装置により浄化される特定成分の濃度を検出する排ガスセンサを配置し、内燃機関の運転中における前記排ガスセンサの出力のデータを用いて前記触媒装置の劣化状態を評価する方法であって、
前記排ガスセンサが検出する前記特定成分の濃度の、前記触媒装置の上流側の排ガスの空燃比に対する相関特性を、該空燃比の変化に対して該特定成分の濃度が極小値を持つ相関特性として表現すべくあらかじめ定めた形式の非線形関数に含まれるパラメータの値を、前記触媒装置の上流側の排ガスの空燃比を表すデータと前記排ガスセンサの出力のデータを用いて同定する同定工程と、その同定したパラメータの値に基づき、前記触媒装置の劣化状態を評価する評価工程とを備えたことを特徴とする排ガス浄化用触媒装置の劣化状態評価方法。 - 前記同定工程で前記非線形関数のパラメータの値の同定のために用いる前記触媒装置の上流側の排ガスの空燃比を表すデータは、該排ガスの空燃比と所定の基準値との偏差を表すデータであることを特徴とする請求項1記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記触媒装置の上流側の排ガスの空燃比を検出する空燃比センサが該触媒装置の上流側に設けられ、前記同定工程は、前記非線形関数のパラメータの値の同定のために、前記触媒装置の上流側の排ガスの空燃比として前記空燃比センサにより検出される空燃比を用いることを特徴とする請求項1又は2記載の触媒装置の劣化状態評価方法。
- 前記同定工程で前記パラメータの値を同定してなる前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値を該非線形関数のパラメータの同定値を用いて求め、その求めた空燃比の値を前記触媒装置の上流側の排ガスの目標空燃比として得る目標空燃比算出工程と、該目標空燃比に前記空燃比センサが検出する空燃比を収束させるようにフィードバック制御により前記内燃機関で燃焼させる混合気の空燃比を操作する空燃比操作工程とを備えたことを特徴とする請求項3記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記同定工程は、逐次型同定アルゴリズムにより前記非線形関数のパラメータの値を逐次同定することを特徴とする請求項1〜4のいずれか1項に記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記同定工程は、前記非線形関数のパラメータのうちの少なくとも一つのパラメータの値を所定の条件を満たす値に制限しつつ同定することを特徴とする請求項5記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記非線形関数は二次関数であり、前記同定工程は、少なくとも該二次関数の最大次数の項に係る係数をその値を制限する前記パラメータとし、該二次関数の関数値により表される前記特定成分の濃度が極小値をもつように前記係数の値を制限しつつ同定することを特徴とする請求項6記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価工程は、前記触媒装置の劣化の進行に応じた変化を呈するものとして前記非線形関数のパラメータを用いてあらかじめ定義した評価パラメータの値を、前記同定工程により同定された前記非線形関数のパラメータの値から求める工程を備え、該評価パラメータの値に基づき、前記触媒装置の劣化状態を評価することを特徴とする請求項1〜7のいずれか1項に記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価パラメータは、前記同定工程で前記パラメータの値を同定してなる前記非線形関数のグラフが前記触媒装置の劣化の進行に伴い空燃比のリーン側に変移していくことに呼応して変化する評価パラメータであることを特徴とする請求項8記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価パラメータは、前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値であることを特徴とする請求項9記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価パラメータは、前記同定工程で前記パラメータの値を同定してなる前記非線形関数のグラフが前記触媒装置の劣化の進行に伴い該非線形関数の関数値により表される前記特定成分の濃度が増加する方向に変移していくことに呼応して変化する評価パラメータであることを特徴とする請求項8記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価パラメータは、前記非線形関数の関数値により表される前記特定成分の濃度の極小値に対応する該非線形関数の関数値の極値であることを特徴とする請求項11記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価パラメータは、前記同定工程で前記パラメータの値を同定してなる前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値よりもリーン側又はリッチ側の空燃比域において該非線形関数のグラフの傾きが前記触媒装置の劣化の進行に伴い大きくなることに呼応して変化する評価パラメータであることを特徴とする請求項8記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価パラメータは、前記非線形関数の関数値により表される前記特定成分の濃度が極小値となる空燃比の値からあらかじめ定めた所定量だけ変化させた空燃比の値における前記非線形関数の関数値であることを特徴とする請求項13記載の排ガス浄化用触媒装置の劣化状態評価方法。
- 前記評価工程で求める前記評価パラメータは複数種類あり、該評価工程は、各種類の評価パラメータに基づく前記触媒装置の劣化状態の評価結果を総合して、該触媒装置の劣化状態の評価結果を確定することを特徴とする請求項8記載の排ガス浄化用触媒装置の劣化状態評価方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP32242899A JP3688533B2 (ja) | 1999-11-12 | 1999-11-12 | 排ガス浄化用触媒装置の劣化状態評価方法 |
US09/708,707 US6427438B1 (en) | 1999-11-12 | 2000-11-09 | Method of evaluating deteriorated state of catalytic converter for purifying exhaust gas |
EP00310055A EP1099836B1 (en) | 1999-11-12 | 2000-11-13 | Method of evaluating deteriorated state of catalytic converter for purifying exhaust gas |
DE60027986T DE60027986T2 (de) | 1999-11-12 | 2000-11-13 | Verfahren zum Auswerten des degradierten Zustands eines Katalysators zur Reinigung von Abgasen |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP32242899A JP3688533B2 (ja) | 1999-11-12 | 1999-11-12 | 排ガス浄化用触媒装置の劣化状態評価方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2001140633A JP2001140633A (ja) | 2001-05-22 |
JP3688533B2 true JP3688533B2 (ja) | 2005-08-31 |
Family
ID=18143568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP32242899A Expired - Fee Related JP3688533B2 (ja) | 1999-11-12 | 1999-11-12 | 排ガス浄化用触媒装置の劣化状態評価方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US6427438B1 (ja) |
EP (1) | EP1099836B1 (ja) |
JP (1) | JP3688533B2 (ja) |
DE (1) | DE60027986T2 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004293338A (ja) * | 2003-03-25 | 2004-10-21 | Mitsubishi Fuso Truck & Bus Corp | NOx吸蔵量の推定方法 |
DE602005019954D1 (de) | 2004-05-27 | 2010-04-29 | Panasonic Corp | Schaltung zur regelung des sendeausgangs und dieselbe verwendende drahtlose vorrichtung |
US8215098B2 (en) * | 2005-05-02 | 2012-07-10 | Cummins Inc. | Method and apparatus for diagnosing exhaust gas aftertreatment component degradation |
US8386547B2 (en) * | 2008-10-31 | 2013-02-26 | Intel Corporation | Instruction and logic for performing range detection |
DE112012006716T5 (de) * | 2012-07-17 | 2015-09-10 | Delphi Technologies Holding S.A.R.L. | Steuervorrichtung für Verbrennungsmotor |
CN115853624A (zh) * | 2023-02-10 | 2023-03-28 | 潍柴动力股份有限公司 | 氮氧化物排放超标的原因检测方法、装置及可读存储介质 |
Family Cites Families (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4432243A (en) * | 1979-02-28 | 1984-02-21 | Ocean Research Equipment, Inc. | Flow calculator with velocity curve fitting circuit means |
US4394351A (en) * | 1981-09-08 | 1983-07-19 | General Motors Corporation | Dual-monolith catalytic converter with secondary air injection |
JP3065127B2 (ja) | 1991-06-14 | 2000-07-12 | 本田技研工業株式会社 | 酸素濃度検出装置 |
JP2678985B2 (ja) | 1991-09-18 | 1997-11-19 | 本田技研工業株式会社 | 内燃エンジンの空燃比制御装置 |
IT1257100B (it) * | 1992-09-14 | 1996-01-05 | Fiat Auto Spa | Sistema di monitoraggio dell'efficienza di un catalizzatore, particolarmente per autoveicoli. |
JP3179920B2 (ja) * | 1993-02-12 | 2001-06-25 | 株式会社日立製作所 | 内燃機関の触媒診断装置 |
US5539638A (en) * | 1993-08-05 | 1996-07-23 | Pavilion Technologies, Inc. | Virtual emissions monitor for automobile |
JP3169298B2 (ja) * | 1993-09-08 | 2001-05-21 | 株式会社日立製作所 | 内燃機関の故障診断装置 |
JP3162553B2 (ja) | 1993-09-13 | 2001-05-08 | 本田技研工業株式会社 | 内燃機関の空燃比フィードバック制御装置 |
JP3325392B2 (ja) | 1994-07-06 | 2002-09-17 | 本田技研工業株式会社 | 内燃エンジンの燃料噴射制御装置 |
US5452576A (en) * | 1994-08-09 | 1995-09-26 | Ford Motor Company | Air/fuel control with on-board emission measurement |
JP3269945B2 (ja) | 1994-08-12 | 2002-04-02 | 本田技研工業株式会社 | 内燃機関の燃料噴射制御装置 |
US5703777A (en) * | 1994-10-20 | 1997-12-30 | Anr Pipeline Company | Parametric emissions monitoring system having operating condition deviation feedback |
JP3422393B2 (ja) | 1995-02-24 | 2003-06-30 | 本田技研工業株式会社 | 内燃機関の空燃比制御装置 |
DE19606652B4 (de) * | 1996-02-23 | 2004-02-12 | Robert Bosch Gmbh | Verfahren der Einstellung des Kraftstoff-Luftverhältnisses für eine Brennkraftmaschine mit nachgeschaltetem Katalysator |
JP3674017B2 (ja) * | 1996-03-19 | 2005-07-20 | 株式会社デンソー | 排出ガス浄化用触媒劣化検出装置 |
US5852930A (en) | 1996-04-05 | 1998-12-29 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio control system for internal combustion engines |
JP3373724B2 (ja) | 1996-04-05 | 2003-02-04 | 本田技研工業株式会社 | 内燃機関の空燃比制御装置 |
JP3299109B2 (ja) | 1996-04-05 | 2002-07-08 | 本田技研工業株式会社 | スライディングモード制御方法 |
JP3261038B2 (ja) | 1996-04-05 | 2002-02-25 | 本田技研工業株式会社 | 内燃機関の空燃比制御装置 |
JP3300598B2 (ja) | 1996-04-05 | 2002-07-08 | 本田技研工業株式会社 | 内燃機関の空燃比制御装置 |
US6151888A (en) * | 1996-06-12 | 2000-11-28 | Robert Bosch Gmbh | Method of diagnosing a catalytic converter |
JPH1068346A (ja) * | 1996-06-21 | 1998-03-10 | Ngk Insulators Ltd | エンジン排ガス系の制御法 |
JP3331161B2 (ja) * | 1996-11-19 | 2002-10-07 | 本田技研工業株式会社 | 排気ガス浄化用触媒装置の劣化判別方法 |
US5732551A (en) * | 1996-12-11 | 1998-03-31 | Motorola Inc. | Method for monitoring the performance of a catalytic converter using mass transfer coefficients |
US5941928A (en) * | 1997-07-31 | 1999-08-24 | Motorola Inc. | System and method for measuring hydrocarbon conversion efficiency of a catalytic converter |
JP3331159B2 (ja) | 1997-09-16 | 2002-10-07 | 本田技研工業株式会社 | プラントの制御装置 |
JP3592519B2 (ja) | 1997-09-16 | 2004-11-24 | 本田技研工業株式会社 | 内燃機関の排気系の空燃比制御装置及びプラントの制御装置 |
JP3354088B2 (ja) | 1997-09-16 | 2002-12-09 | 本田技研工業株式会社 | 内燃機関の排気系の空燃比制御装置 |
ES2296392T3 (es) * | 1998-07-17 | 2008-04-16 | Honda Giken Kogyo Kabushiki Kaisha | Metodo de determinar un estado deteriorado de un convertidor catalitico para purificar los gases de escape. |
-
1999
- 1999-11-12 JP JP32242899A patent/JP3688533B2/ja not_active Expired - Fee Related
-
2000
- 2000-11-09 US US09/708,707 patent/US6427438B1/en not_active Expired - Fee Related
- 2000-11-13 EP EP00310055A patent/EP1099836B1/en not_active Expired - Lifetime
- 2000-11-13 DE DE60027986T patent/DE60027986T2/de not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
EP1099836B1 (en) | 2006-05-17 |
EP1099836A3 (en) | 2003-02-05 |
DE60027986D1 (de) | 2006-06-22 |
DE60027986T2 (de) | 2006-10-19 |
EP1099836A2 (en) | 2001-05-16 |
JP2001140633A (ja) | 2001-05-22 |
US6427438B1 (en) | 2002-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3967524B2 (ja) | 内燃機関の空燃比制御装置 | |
JP3992925B2 (ja) | 排ガス浄化用触媒装置の劣化判別方法 | |
JP3592519B2 (ja) | 内燃機関の排気系の空燃比制御装置及びプラントの制御装置 | |
US5715796A (en) | Air-fuel ratio control system having function of after-start lean-burn control for internal combustion engines | |
US6925372B2 (en) | Control apparatus, control method, and engine control unit | |
JP3621839B2 (ja) | プラントの制御装置 | |
JP4265704B2 (ja) | 内燃機関の空燃比制御装置及びプラントの制御装置 | |
JP4354068B2 (ja) | 内燃機関の排ガスの空燃比制御装置 | |
JP3331159B2 (ja) | プラントの制御装置 | |
JP3773684B2 (ja) | 内燃機関の空燃比制御装置 | |
JP3782269B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4312325B2 (ja) | 排ガス浄化用触媒装置の劣化状態評価方法 | |
JPH1193740A (ja) | 内燃機関の排気系の空燃比制御装置 | |
JP4490000B2 (ja) | 内燃機関の空燃比制御装置 | |
JP3655145B2 (ja) | 多気筒内燃機関の空燃比制御装置 | |
JP4437626B2 (ja) | 内燃機関の空燃比制御装置 | |
JP3688533B2 (ja) | 排ガス浄化用触媒装置の劣化状態評価方法 | |
JP2000234550A (ja) | 内燃機関の空燃比制御装置 | |
JP3696570B2 (ja) | プラントの制御装置 | |
JP3748418B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4074142B2 (ja) | 内燃機関の空燃比制御装置 | |
JPH08232726A (ja) | 内燃機関の制御装置 | |
JP3743591B2 (ja) | 内燃機関の空燃比制御装置 | |
JP2003020983A (ja) | 内燃機関の排気系の空燃比制御装置及びプラントの制御装置 | |
JP2006177369A (ja) | 制御装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20041125 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050201 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050316 |
|
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: 20050607 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050608 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080617 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090617 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |