(57)【要約】【課題】特定のための段階を必要とせず、かつエンジンのすべての動作点についてロバスト性のより高い空燃比推定モデルを実現する、排気プロセスのより詳細なモデル化を可能にする。【解決手段】本発明は、空燃比を測定する検出器を含む排気回路を有する内燃エンジンの各シリンダ内の空燃比を推定する方法に関する。非線形適応フィルタに基づく推定器が、複数のシリンダからの気体の排出と、それらの気体の検出器までの排気回路内の移動とを表している物理モデルに結合されている。推定器は、排気ガスの総質量および新気の質量などの前記モデルの少なくとも１つの変数から測定される空燃比の推定値にも結合されている。【選択図】図３
2006336644.tif 000002
【特許請求の範囲】
【請求項１】
マニフォールドに連結されている少なくとも複数のシリンダと、前記マニフォールドの下流で空燃比（λ）を測定する検出器とを含む排気回路を有する内燃エンジンの各シリンダ内の空燃比を推定する方法において、
複数の前記シリンダの各々からの気体の排出と、前記排気回路内における前記検出器までの気体の移動とをリアルタイムで表す物理モデルを定めるステップと、
前記物理モデルの少なくとも１つの変数から、前記検出器によって測定される前記空燃比（λ）の推定値を定めるステップと、
前記物理モデルを、前記測定される空燃比（λ）の測定値の前記推定値が考慮される適応型の非線形推定器に結合するステップと、
前記適応型の非線形推定器によって前記各シリンダ内の前記空燃比の値をリアルタイムで推定するステップと、
を有することを特徴とする、内燃エンジンの各シリンダ内の空燃比を推定する方法。
【請求項２】
気体の通過時間と前記検出器の応答時間とによる遅延時間を、決められたシリンダ内で試験的な外乱を発生させることと、その影響を前記検出器の位置で測定することとによって求める、請求項１に記載の方法。
【請求項３】
前記物理モデルは、前記排気マニフォールド内の気体の総質量（ＭT）と、前記排気マニフォールド内の新気の質量（Ｍair）と、前記各シリンダ内の空燃比（λi）との、少なくとも３種類の変数を有している、請求項１または２に記載の方法。
【請求項４】
前記物理モデルは、前記排気マニフォールド内の気体の総質量（ＭT）と、前記シリンダから流出する質量流量（ｄi）との、少なくとも２種類の出力データを有している、請求項１から３のいずれか１項に記載の方法。
【請求項５】
前記測定される空燃比（λ）を、前記排気マニフォールド内の気体の総質量（ＭT）と前記排気マニフォールド内の新気の質量（Ｍair）との関数として推定する、請求項１から４のいずれか１項に記載の方法。
【請求項６】
前記各シリンダ内の前記空燃比の前記推定値は、前記排気マニフォールド内の気体の総質量（ＭT）の推定値、前記排気マニフォールド内の新気の質量（Ｍair）の推定値、およい前記各シリンダ内の前記空燃比（λi）の推定値の、リアルタイムの修正値を含んでいる、請求項１から５のいずれか１項に記載の方法。
【請求項７】
すべてのシリンダ内の空燃比を調節するために各シリンダ内に噴射される燃料質量を変更するエンジン制御への、請求項１から６のいずれか１項に記載の方法の適用。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、噴射内燃エンジンの各シリンダの空燃比をマニフォールドの下流の空燃比の測定値から推定する方法および適応非線形フィルタに関する。
【背景技術】
【０００２】
空気の質量に対する燃料の質量によって特徴づけられる空燃比についての理解は、空燃比が１に近いときに混合気が良好な燃焼状態になるガソリンエンジンを備えているか、希薄混合状態（１未満の空燃比）で動作する限り空燃比を知ることについての関心が異なるディーゼルエンジンを備えているかによらず、すべての乗物に対して重要である。特に、ＮＯｘトラップを使用している触媒は、時間の経過と共に効率が低下する。最適な効率を回復するには、希薄混合気での正常運転状態に戻る前に、空燃比を数秒間１付近に維持しなければならない。そのため、ＤｅＮＯｘ触媒による汚染物除去は、シリンダ毎の空燃比の正確な制御が必要である。
【０００３】
そのため、ＮＯｘトラップの上流の（過給エンジンの）タービンの出口に配置されているプローブは、排気過程の結果としての平均空燃比の測定値を出力する。高度にフィルタリングされる、ノイズの影響を受けているこの測定値は、空燃比が１の期間中にシリンダ内に噴射される燃料の質量の制御に使用され、そのため、各シリンダには同じ質量の燃料が噴射される。
【０００４】
シリンダ内への燃料の質量をより正確に、特に個別に制御するためには、シリンダ毎に空燃比の再現が必要である。空燃比プローブを乗物の各シリンダの出口に取り付けることは、空燃比プローブの原価を考えれば不可能であるので、１つのプローブから出力される測定値を基に動作する推定器は、各シリンダ内の空燃比を別個に知ることができるので有利である。
【０００５】
したがって、エンジン制御では、すべてのシリンダの空燃比が均衡するように、再現された複数の空燃比から各シリンダ内に噴射する燃料の質量を調節することができる。
【０００６】
特許文献１は、モデルの定義と、その観測と、カルマンフィルタによるフィルタリングとについて記述している。このモデルは、マニフォールド内の混合気の物理的記述を含んでおらず、脈動が非常に大きい流量の現象を考慮していない。
【０００７】
シリンダ内の空燃比の推定は、行列の複数の係数によってのみ調節されており、これらの係数は最適化アルゴリズムによってオフラインで特定することができる。さらに、行列のさまざまな調節、そのための行列のパラメータの特定は、各動作点（エンジン回転数／負荷）に対応している。したがって、この推定器は（５つの空燃比プローブを備えた）大規模な取得試験手段が必要で、エンジンの状態が変化する場合のロバスト性を有していない。
【特許文献１】仏国特許出願公開第２，８３４，３１４号明細書
【発明の開示】
【発明が解決しようとする課題】
【０００８】
本発明の目的は、特定のための段階を必要とせず、かつエンジンのすべての動作点についてロバスト性のより高い空燃比推定モデルを実現する、排気プロセスのより詳細なモデル化を可能にすることである。本発明は、さらに、クランク軸の６°の回転毎の測定を可能にし、したがって、測定ノイズに影響されることなく空燃比の測定情報を高頻度に得ることを可能にする。
【課題を解決するための手段】
【０００９】
したがって、本発明は、マニフォールドに連結されている少なくとも複数のシリンダと、前記マニフォールドの下流で空燃比（λ）を測定する検出器とを含む排気回路を有する内燃エンジンの各シリンダ内の空燃比を推定する方法において、
複数の前記シリンダの各々からの気体の排出と、前記排気回路内における前記検出器までの気体の移動とをリアルタイムで表す物理モデルを定めるステップと、
前記物理モデルの少なくとも１つの変数から、前記検出器によって測定される前記空燃比（λ）の推定値を定めるステップと、
前記物理モデルを、前記測定される空燃比（λ）の測定値の前記推定値が考慮される適応型の非線形推定器に結合するステップと、
前記適応型の非線形推定器によって前記各シリンダ内の前記空燃比の値をリアルタイムで推定するステップと、
を有することを特徴とする、内燃エンジンの各シリンダ内の空燃比を推定する方法に関する。
【００１０】
また、気体の通過時間と前記検出器の応答時間とによる遅延時間を、決められたシリンダ内で試験的な外乱を発生させることと、その影響を前記検出器の位置で測定することとによって求めてもよい。
【００１１】
一態様によれば、前記物理モデルは、前記排気マニフォールド内の気体の総質量（ＭT）と、前記排気マニフォールド内の新気の質量（Ｍair）と、前記各シリンダ内の空燃比（λi）との、少なくとも３種類の変数を有していてもよい。この態様はまた、前記排気マニフォールド内の気体の総質量（ＭT）と、前記シリンダから流出する質量流量（ｄi）との、少なくとも２種類の出力データを有していてもよい。
【００１２】
前記測定される空燃比（λ）を、前記排気マニフォールド内の気体の総質量（ＭT）と前記排気マニフォールド内の新気の質量（Ｍair）との関数として推定してもよい。
【００１３】
さらに、前記各シリンダ内の前記空燃比の前記推定値は、前記排気マニフォールド内の気体の総質量（ＭT）の推定値のリアルタイムの修正値と、前記排気マニフォールド内の新気の質量（Ｍair）の推定値のリアルタイムの修正値と、前記各シリンダ内の前記空燃比（λi）の推定値の修正値とを含んでいてもよい。
【００１４】
最後に、本方法を、すべてのシリンダ内の空燃比を調節するために各シリンダ内に噴射される燃料質量を変更するエンジン制御に適用してもよい。
【発明を実施するための最良の形態】
【００１５】
すべてのシリンダの空燃比の推定値の平均に関して各シリンダ内での個別の空燃比を推定する利点は、非常に多い。
−推定をタービンの出口の１つの空燃比プローブを基に行う場合のコスト価格
−排出物質の削減
−運転のしやすさ（伝達されるトルクの調節）の改善
−燃費の低減
−燃料噴射装置の診断（噴射ノズルのドリフトの検出や噴射装置の故障の検出）
排気プロセスの説明
排気プロセスは、気体が排気弁から外気（排気サイレンサの出口）まで移動する経路を有している。本実施形態のエンジンの例は、２２００ｃｍ3の４シリンダエンジンである。このエンジンには、可変容量ターボ過給器が備えられている。図１の模式図は、排気プロセスの具体的な構成部材を示している。
−λ1からλ4は、４つのシリンダの各々の内部の空燃比である。
−ＳＲは空燃比プローブである。
−ＣＥは排気マニフォールドに対応している。
−Ｔはターボ過給器のタービンに対応している。
−ＤＳ１からＤＳ４は、シリンダの出口での流量を表している。
【００１６】
空燃比プローブ（ＳＲ）は、タービン（Ｔ）の直後に配置されている。シリンダ内での燃焼後の気体は、以下の動作を行う。
−排気弁を通過する。排気弁はカム軸によって制御され、そのリフトローは釣り鐘形状である。流量は、弁が開いているときの高い値から、シリンダとマニフォールドの圧が等しくなったときの低い値まで変化し、ピストンが排気ガスを排出するように再び上に摺動し始めた時に再び増加する。
−マニフォールドをシリンダヘッドの出口に連結している短いパイプを通過する。
−４つのシリンダからの流れ（ＤＳ１からＤＳ４）が出会う排気マニフォールド（ＣＥ）内での混合相。ここで、各シリンダからの排気の一吹きが、マニフォールドの型（対称または非対称）と、流れの重なる割合を定めるＥＥＯ（Early Exhaust Opening）と、ＬＥＣ（Late Exhaust Closing）とによって混合される。
−吸気部の上流に配置されている過給器に必要なトルクを供給するタービンを通過する。流量に対するタービンの動作の影響はよくわかっていないが、さまざまなシリンダからの一吹きの排気をより多く混合することになると考えることができる。
−ＵＥＧＯ型のプローブによる測定。
【００１７】
排気ガスの組成は、燃焼室に供給される燃料および空気の量と、燃料の組成と、燃焼の進行とによって決まる。
【００１８】
実際には、空燃比プローブは多孔質の物質から作られている拡散バリアによって排気管に連結されている拡散室内のＯ2の濃度を測定する。この構成では、特に空燃比プローブの付近の温度および／または圧力の変動のため、プローブ用に選択された位置によって空燃比に差が生じることがある。
【００１９】
しかしながら、必要なのはシリンダ間の空燃比の差異の検出であり、平均値は推定器によって通常は維持されているため、圧力または温度に依存しているこの空燃比の変動現象は無視されている。
【００２０】
本発明の推定器によって用いられるリアルタイム物理モデルでは、測定される空燃比（λ）はプローブの周辺の空気の質量（すなわち空気の流れ）と総質量（すなわち総流量）とに関連している。このモデルは、空気、燃料および既燃ガスの３つの気体によるアプローチに基づいている。したがって、燃焼後に残ったすべての気体は、希薄混合気の場合、空気と既燃ガスとの混合物であると考えられる。過濃混合気については、燃料が過剰であり、燃焼後には未燃焼の燃料と既燃ガスとが残るのに対して、空気はすべて消滅している。実際には、燃焼は１００％完全ではないが、この推定器については完全と見なされている。
【００２１】
前述の３種類の質量への空燃比の関連付けを定式化する。希薄混合気の場合、空気が過剰で、燃焼後は燃料は残っていない。燃焼前に、以下の質量がシリンダ内に存在していると仮定する。
【００２２】
空気の質量Ｍair、燃料の質量Ｍcarb、および既燃ガスの質量ＭgazB
Ｍair＝ｘ；Ｍcarb＝ｙ；ＭgazB＝０
理論空燃比を達成するには燃料の１４．７倍の空気が必要なことがわかっていると、燃焼前後の各種の質量を示した以下の表を作成することができる。
【００２３】
【表１】
2006336644.tif 000003
【００２４】
以下の式の計算によって得られる、比Ｍcarb／Ｍairを表している空燃比λは、混合気が希薄混合気であるときにだけ有効である。
【００２５】
【数１】
2006336644.tif 000004
【００２６】
ここで、ＰＣＯは、混合が理論空燃比の時の比Ｍcarb／Ｍairに対応している。ＰＣＯは、燃料の発熱量である。
【００２７】
過濃混合気の場合、式は以下のようになる。
【００２８】
【数２】
2006336644.tif 000005
【００２９】
しかしながら、吸気部に既燃ガスが存在すると、３つの気体の濃度が排気部で変わるので、これらの式は混合気がＥＧＲを含んでいない場合に有効である。
【００３０】
本実施形態では、式（７）の空燃比を組み入れるために希薄混合気用の空燃比の式だけが推定器で使用され、非常にわずかな空気（＜３％）だけが無視される。しかし、本発明は、本実施形態にだけ限定されているわけではなく、実際、式は空燃比が１の付近では連続しており、式の逆数を取ったものは過濃混合気に対して何ら問題がない。
【００３１】
排気管内で気体が混合する様子をより良く理解するために、ディーゼルエンジンのモデルをＩＭＡＧＩＮＥ社（フランス）のＡＭＥＳｉｍソフトウェアと共に使用した。このモデルは、逆数を取ることはできないが、本発明のモデルの確認するための基準として使用されることになる。
【００３２】
ＡＭＥＳｉｍは、０Ｄモデル化のソフトウェアで、特に熱および液圧の現象に適している。ＡＭＥＳｉｍは、特に、容積部分、管、または絞りをモデル化することができる。
【００３３】
排気モデルは以下を有している。
−容積部分と管とによって表されている複数の排気管
−熱交換を備えた排気マニフォールド
−タービンとバイパス弁
−タービンの合流部の容積部分と弁の流量
−タービンと測定プローブとの間の管
−排気経路の容積部分と管
管、絞り、および容積部分をモデル化する基本ブロックは、ＡＭＥＳｉｍの取扱説明書、「Ｔｈｅｒｍａｌ Ｐｎｅｕｍａｔｉｃ Ｌｉｂｒａｒｙ」に記載されている。絞りを通過する流量、質量およびエネルギーの保存を計算するために、標準の式が使用される。さらに、モデルは気体の慣性も考慮しており、これは気体の組成の動力学の検討には重要である。
【００３４】
このモデルは０Ｄモデルのため、時間の次元は考慮されず、物理的アプローチによって遅れ時間をモデル化することはできない。入力変数を変化させると、出力は即座に変化する。したがって、移動時間は無視される。この制限は、リアルタイムでの取得に取り組む場合に、重要である。
【００３５】
本発明によれば、１つのリアルタイム物理モデルが、大局的なシステム、つまり排気ガスがシリンダからマニフォールドを通ってタービンの下流の排気部まで移動する経路全体をモデル化するために定められる。
【００３６】
Ｉ−リアルタイム物理モデルの定義
本実施形態では、温度の変動は１つのエンジンサイクルを通して低く、流量の変動への温度の影響は限られていると考えられる。圧力の変動は、流量に直接関連しているため、実際にはプロセスにおいて不可欠である。したがって、シリンダ、マニフォールド、およびタービンの各要素に対して一定の温度が設定されている。そのため、熱交換もモデル化されていない。この簡略化のための仮定による大きな影響はない。
【００３７】
第１のアプローチでは、新気と既燃ガスとの２つの気体を考慮する。従来の式は、容積部分内の気体の総質量の変動と、新気の質量の変動とを示している。そのため、既燃ガスをそれらから求めることができる。この手順は希薄混合状態の場合に有効であるが、過濃混合気の場合にも同様の式を燃料と既燃ガスについて示すことができる。
【００３８】
Ａ）排気マニフォールドの物理モデル
排気マニフォールドは、質量が保存される容積部分に従ってモデル化されている。温度は、実質的に一定であると仮定して、エンジン回転数と負荷との関数としてグラフから求められている。
【００３９】
本発明によれば、測定される空燃比をプローブの周辺の空気の質量と総質量とに関連付けることにしている。したがって、マニフォールド内の総質量の保存は、マニフォールド内の排気ガスの質量はマニフォールドに流入する排気ガスの質量（シリンダの出口の流量）からマニフォールドから流出する質量を引いたものに等しいという事実を表している。タービン内の流れの組成は、マニフォールドの出口と同じと仮定している。したがって、マニフォールドから流出する質量は、タービンを通過する流れに等しい。したがって、総質量については以下の式が得られる。
【００４０】
【数３】
2006336644.tif 000006
【００４１】
ここで、
Ｎe：エンジン回転数
α：クランク軸の角度
ＭT：排気マニフォールド内の総質量
ｄi：シリンダｉから流出する質量流量
ｄT：タービンを通過する全流量
である。
【００４２】
同様に、空気の質量の保存については、以下の式が得られる。
【００４３】
【数４】
2006336644.tif 000007
【００４４】
ここで、
Ｎe：エンジン回転数
α：クランク軸の角度
Ｍair：排気マニフォールド内の新気の質量
λi：各シリンダ内の空燃比
ｄi：シリンダｉから流出する質量流量
ｄair：タービンを通過する空気の流量
である。
【００４５】
シリンダの出口の流量とタービンを通過する流量とを求めることができる物理モデルについて説明する。
【００４６】
シリンダの出口での流量を求めることができるモデル：気体の排出
シリンダの出口での気体の流量は、排気弁の出口での流量を示す物理モデルによってモデル化することができる。３つの変数が、弁を通過する気体のこの排出モデルに使用される。
−クランク軸の角度（α）
−シリンダを通じて吸い込まれる流れｄasp（上流のエンジン制御によって推定される変数）
−１サイクルにわたってプローブによって測定される空燃比の平均値
流出する流れの平均値は、吸入される流れと噴射されたガソリンの流量から知ることができる。流出する流れの瞬間値は、吸入される流れに依存したテンプレートに基づいて得られる。このテンプレートは、エンジン回転数、クランク軸の角度、シリンダによって吸入される流れ、および１サイクルにわたってプローブによって測定される空燃比の平均値から、クランク軸の角度の関数として、シリンダについての平均流量を推定することを可能にする経験則に基づいた物理モデル（曲線）である。この物理法則の唯一の制約は、流出する流れの平均値（曲線の領域）を考慮し、以下の２つの現象を明らかにする曲線を提供することである。
−クランク角度の関数として流量のピークによって表されているシリンダ／排気の圧力の均衡
−振幅がより小さい流量の第２のピークで表される、排気弁の流れの断面に依存する流量
このテンプレート
【数５】
2006336644.tif 000008
【００４７】
は、つまり排気弁の出口での質量流量の経過（曲線）ｄi、つまりすべてのシリンダに共通の流量の推定値
【数６】
2006336644.tif 000009
【００４８】
を出力する。これは、エンジンのテストベンチでの測定値と関連させて求められる。それから、吸入される流れとプローブで測定される空燃比の平均値
【数７】
2006336644.tif 000010
【００４９】
とによって、テンプレートを相似変形し、クランク軸角度の関数として各シリンダに対して位相を移動させ、各シリンダの出口での気体の流量を推定する。
【００５０】
ｄi（α）＝ｄ（αi＋α）・α0
ここで、
−ｄi（α）はシリンダｉの出口で流出する気体の流量
−ｄ（α）はテンプレート、つまりシリンダの出口での流量の推定値
【００５１】
【数８】
2006336644.tif 000011
【００５２】
−αiはシリンダｉの位相ずれの角度
である。
【００５３】
テンプレート曲線の位相ずれは、図１（ＤＳ１からＤＳ４）において概略的に示されている。
【００５４】
タービンを通過する流量を求めることができる複数の物理モデルを以下に説明する。
【００５５】
タービンを通過する流量を求めることができるモデル：タービンモデル
タービンは、流量制限部を通過する流れによってモデル化されている。タービン内の流量は、タービンの回転数と、タービンの上流と下流の圧力の比との関数としてのマッピング（チャート）によって一般的に求められる。
【００５６】
タービンを通過する流量ｄTは、排気マニフォールド内の排気ガスの総質量（ＭT）と、排気マニフォールド内の温度と、ターボ過給器の回転数と、ターボ過給器の形状寸法との関数である。したがって、このモデルの入力データは、
−排気ガスの総質量（ＭT）
−空気の質量（Ｍair）
−エンジン回転数（Ｎe）
−（ターボ過給器の）タービンの回転数
である。
【００５７】
この流量は、総質量ＭTの凹関数から推定することができる。この関数をｐで表す。したがって、タービン内の流量は、次のように示される。
【００５８】
ｄT（ＭT）＝ＭT・ｐ（ＭT）
関数ｐは、タービンの回転数の関数として表され、かつ大気状態での排気マニフォールド内の質量（Ｍ0）に対する排気マニフォールド内の排気ガスの総質量（ＭT）の比として表される、ルートの形態の関数である。したがって、マッピングによって、ＭT／Ｍ0の比と、（ターボ過給器）のタービンの回転数との関数として、ｐ（ＭT）が得られる。このマッピングで使用される式は次のとおりである。
【００５９】
【数９】
2006336644.tif 000012
【００６０】
ここで、
−ｆは多項式関数
−ｇは定数
である。
【００６１】
関数ｆのパラメータは、タービンのマッピングとの関連付けによって最適化されている。
【００６２】
さらに、空気の組成は排気マニフォールド内と同じであると仮定している。したがって、タービンを通過する空気の流れは、次の式で表される。
【００６３】
【数１０】
2006336644.tif 000013
【００６４】
したがって、物理的な気体の排出とタービンのモデルとによって、式（１）と式（２）とは次のように示すことができる。
【００６５】
【数１１】
2006336644.tif 000014
【００６６】
これらの式（３）のシステムは、排気マニフォールドの物理モデルである。
【００６７】
このモデルの入力データは、
Ｎe：エンジンの回転数
α：クランク軸の角度
ｄi：シリンダｉから流出する質量流量
ｄT：タービンを通過する全流量
ｄair：タービンを通過する空気の流量
λi：各シリンダ内の空燃比
であり、システムの未知数は、
ＭT：排気マニフォールド内の総質量
Ｍair：排気マニフォールド内の新気の質量
である。
【００６８】
第１の式は１つの未知数ＭTを含んでいる。第２の式は２つの未知数Ｍairおよびλiを含んでいる。これは、以下に説明するさらなる仮定につながる。
【００６９】
Ｂ）シリンダの出口の空燃比についての仮定
排気マニフォールドのリアルタイム物理モデル（ＲＴＭ）を完成するために、シリンダの出口での空燃比は、動作点で一定であると仮定し、その結果以下の式が得られる。
【００７０】
【数１２】
2006336644.tif 000015
【００７１】
実際に、リアルタイムでの計算では、λiが定数であると推定されている。
【００７２】
Ｃ）リアルタイム物理モデルの表現
最後に、リアルタイム物理モデルＲＴＭは行列形式で以下のように表すことができる。
【００７３】
【数１３】
2006336644.tif 000016
【００７４】
物理モデルの未知数は、最終的に、ＭT，Ｍairおよびλiである。
【００７５】
物理モデルの出力データはＭTおよびｄiである。
【００７６】
Ｄ）排気の遅延時間
管内および様々な容積部分内における気体の移動による遅延時間と、測定プローブの「アイドル時間」とは、前述の物理モデル（式（５）のシステム）では考慮されていない。しかしながら、管内の移動が無視されるので、物理モデルはこれらの遅延時間に関して線形に構成されている。そのため、これらの遅延時間は、排気プロセス全体について１つの遅延時間にまとめることが可能であり、以下に説明するように遅延時間の影響を後で考慮することができるため、物理モデルの逆数を取ることができる。
【００７７】
ＩＩ−空燃比の推定器
上記の物理モデル（５）は、タービンの下流の空燃比（マニフォールド内の空燃比と同一であると考えられる）は、排気マニフォールドの入口での気体流の組成の関数として表されることを示している。
【００７８】
測定されるデータは、
−プローブで測定される空燃比：λ
である。
【００７９】
システムのそれ以外の既知のデータは、
−エンジン回転数：Ｎe
−クランク軸の角度：α
−（過給器の）タービンの回転数
−シリンダに吸入される流れ
である。
【００８０】
システムのモデル化されたデータは、
−シリンダｉからの質量流量：ｄi
−タービンを通過する全流量：ｄT
−タービンを通過する空気の流量：ｄair
−排気マニフォールド内の総質量：ＭT
である。
【００８１】
したがって、未知数は、
−４つのシリンダの各々における空燃比：λi
−排気マニフォールド内の新気の質量：Ｍair
である。
【００８２】
物理モデル（５）は非線形であり、そのようなシステムはリアルタイムに解くことはできない。そのため、システムの未知数を直接計算しようとするよりも、推定器を使用することが必要である。本発明による推定器の選択は、システムの構造がシリンダ内の空燃比λiの関数として線形である（空気の質量の変動がλiの関数として線形である）ことに基づいている。この条件において特に適切な技法は、適応フィルタを使用することである。したがって、物理モデルＲＴＭから未知数を推定するために、本発明の方法は、適応フィルタに基づいて推定器を構築することを提案している。この推定器によって、最終的に、空燃比の推定値をタービンの下流にある検出器によって測定される空燃比からシリンダ毎に求めることができる。
【００８３】
大まかに言えば、適応フィルタはノイズを含むデータに適用して特定の時刻ｔにおける有用な情報を取得するシステムであり、これらのシステムは以下の構成で使用される。
−フィルタリング、つまり、時刻ｔでの時刻ｔを含む時刻ｔまでに測定されるノイズを含むデータからの有用な情報の抽出
−時刻ｔの後に測定されるデータも用いる平滑化
−時刻ｔでの有用な情報を求めるための時刻ｔ−τまでのデータのみを使用する推定
この実施形態の場合、目的は、測定値と考えるｙ1＝ＭTおよびｙ2＝λの２つの要素から以下のデータを再現することである。
ＭT：排気マニフォールド内の気体の総質量
Ｍair：排気マニフォールド内の新気の質量：
λi：各シリンダ内の空燃比
推定器の、測定される入力パラメータまたはモデル化された入力パラメータは、次のようになる。
【００８４】
【数１４】
2006336644.tif 000017
【００８５】
λはプローブで測定されており、ＭTはリアルタイムモデルＲＴＭ（５）に基づいて推定されている。
【００８６】
次のように設定すると、
【００８７】
【数１５】
2006336644.tif 000018
【００８８】
推定器は以下のように示すことができる。
【００８９】
【数１６】
2006336644.tif 000019
【００９０】
言い換えれば、これは、推定される各パラメータに修正項，，が適用されるリアルタイムモデル（５）である。
【００９１】
推定器の原理は、物理モデル（５）と、その結果得られる空燃比λiとを実際の値に収束させることである。実際、物理モデル（５）は、ＭTとＭairとを出力し、さらに入力パラメータＹが存在している。したがって、推定器はリアルタイムモデルＲＴＭの複数の出力値を複数の入力値と比較して、適切な修正を行う。たとえば、空燃比λiは、ＭTとλとの誤差の関数として修正されなければならず、入力値ＭTおよびλと、対応する推定値
【数１７】
2006336644.tif 000020
および
【数１８】
2006336644.tif 000021
との間の誤差がマイナスの場合、推定値を増加させなければならず、また、その誤差がプラスの場合は減少させなければならない。したがって、次のように示される。
【００９２】
【数１９】
2006336644.tif 000022
【００９３】
同じ原理をＭTおよびＭairの修正に適用することによって、以下の式が得られる。
【００９４】
【数２０】
2006336644.tif 000023
【００９５】
ここで、Ｌ1，Ｌ2，Ｌλは、３つの未知数の解への収束速度を調節することを可能にする修正パラメータである。これらは、厳密に正の実数のパラメータである。これらのパラメータは、収束の速度と測定ノイズに対する低い感度との間の良好な妥協点が得られるように、マニュアルで調節することができる。
【００９６】
最後に、上述の空燃比の式を使用することによって、ＭTおよびＭairから測定される空燃比を推定することができる。
【００９７】
【数２１】
2006336644.tif 000024
【００９８】
Ｍair＝（１−λ）ＭTと推定する。
【００９９】
ｙ1＝ＭTおよびｙ2＝λという表記を使用することによって、適応フィルタと、ＭTおよびＭairに基づくの空燃比の推定値とを基にした物理モデルＲＴＭの推定器は次のように書くことができる。
【０１００】
【数２２】
2006336644.tif 000025
【０１０１】
ここで、δT＝１である。
【０１０２】
このように構成された推定器によって、リアルタイムモデルによって得られた最初の値ＭTとプローブによって測定される空燃比とから、ＭT，Ｍairおよびλをリアルタイムで修正することが可能である。
【０１０３】
システム（８）はリアルタイムで数値的に解くことができ、計算には当業者に公知の前進オイラー法による離散化を使用する。
【０１０４】
シミュレーションの結果：推定器の試験（８）
既知の個別の空燃比から、基準モデルＡＭＥＳｉｍによってプローブでの空燃比（λ）を推定した。この空燃比の値は、推定器の入力で使用される。プローブの動力学は考慮していない。噴射を不均衡にして、タービンの下流で測定される空燃比（λ）からシリンダ毎の空燃比の推定値（λi）を観測する。
【０１０５】
シミュレーションについては、８０μｓ長い噴射を１つのシリンダで行うことによって、４つのシリンダを連続的に不均衡にすることができ、したがってシリンダ１と４とが同様に不均衡になった。図２Ａおよび図２Ｂは、ＡＭＥＳｉｍによって得られた時間（Ｔ）の関数としての空燃比
【数２３】
2006336644.tif 000026
を下側に、時間（Ｔ）の関数としての推定器
【数２４】
2006336644.tif 000027
の推定結果を上側に示している。４つの曲線は４つのシリンダの各々に対応している。適応フィルタに基づく推定器の性能は非常に良好である。しかし、このモデルでは考慮していない気体の慣性によるわずかな位相差が認められる。そのため、排気の遅れ時間の推定器によって、このモデルおよび推定器を完成することとする。
【０１０６】
排気遅延時間の推定器
上述のように構築された推定器では、シリンダの排気部とプローブで取得された信号との間の遅延時間を考慮する推定方法を行うことができない。実際には、遅延時間には、管内および容積部分を通過する移動時間、測定プローブのアイドル時間などのいくつかの原因がある。
【０１０７】
推定器の入力において遅延時間Ｄをモデルからの変数に適用することによって、推定器を空燃比の測定値と同期させることができる。遅延時間を有する推定器の構成を図３に示しており、ここで、
−Ｎeおよびαは、式（５）で表されるリアルタイムモデルＲＴＭの入力データ
−ＭＭＢＯは、開ループ質量モデル（モデルＲＴＭ）
−Ｄは、モデルＲＴＭ（ＭＭＢＯ）の複数の出力変数に適用される遅延時間であり、この遅延時間は式（９）から得られる
−ＳＲは、式（７）によって推定器で使用される空燃比をタービンの下流で測定するプローブ
−ＥＲＦＡは、適応フィルタに基づいた、式（８）で示される空燃比の推定器
−λiは、推定器ＥＲＦＡによって推定されるシリンダｉの空燃比
である。
【０１０８】
遅延時間は、エンジンの回転数、負荷、排気マニフォールドの圧力などの動作条件によって決まる。遅延時間はモデル化が困難なため、推定器と測定値との間の遅延時間を他の計器を使用せずにリアルタイムに計算するための同定方法が開発された。この原理は、シリンダ１の噴射点の近傍でわずかな増量を行い、各シリンダについて空燃比の変動を計算することにある。そのため、シリンダ２，３および４の変動にペナルティを課すように、同定基準Ｊkを設定する。
【０１０９】
【数２５】
2006336644.tif 000028
ここで、
【０１１０】
【数２６】
2006336644.tif 000029
である。
【０１１１】
ペナルティは、βによって課される。シリンダ２について推定された空燃比の値に正の変動がある場合、推定器と測定値との間の遅延は正である。シリンダ３に変動がある場合、遅延は負であり、ペナルティは負である。シリンダ４の変動は、正または負の遅延時間の結果と考えられる。モデルＲＴＭの複数の出力変数に適用される遅延時間Ｄは、追加の遅延であり、Ｊkを最小化することによる最小２乗法によって計算される。
【０１１２】
規準Ｊkは、ＰＩ（比例積分）制御器によって、推定器の遅延時間でゼロに調節される。制御器が安定すると、推定された空燃比の変動はシリンダ１で最大で、シリンダ４で最小になる。そして、推定器は測定値と位相が一致する。
【０１１３】
結果
図４Ａおよび図４Ｂは、上述の推定器による、中負荷および１５００ｒｐｍでの空燃比の推定値をシリンダ毎に示している。これらの図では、時間（Ｔ）の関数としての基準空燃比
【数２７】
2006336644.tif 000030
を上側に、時間（Ｔ）の関数としての推定器による結果
【数２８】
2006336644.tif 000031
を下側に示している。４つの曲線は４つのシリンダの各々に対応している。
【０１１４】
本発明は、プローブによって測定される空燃比（λ）とマニフォールド内の総質量（ＭT）の情報とから、４つのシリンダの出口での空燃比（λi）の推定を可能にする推定器の構成を有する推定方法に関する。このようにして得られた推定器は能率的であり、特に、動作点が変化しても更なる調節は必要ない。同定段階は必要なく、１回のノイズの測定とモデルの調節だけを行えばよい。
【０１１５】
本発明による推定のロバスト性を動作状態によらずに高めるために、遅れ時間制御器が推定器と共に使用され、シリンダでの噴射時間を増加させた後に遅れ時間を再調節することを可能にしている。これによって、たとえば空燃比が１の段階の前に、推定器を最適に校正することができる。
【０１１６】
本発明によれば、クランク軸の６°の回転ごとに測定を実行し、それによって、測定ノイズに影響されずに空燃比の測定値の情報を高頻度に得ることが可能である。さらに、この高頻度の表現によって、システムの脈動効果を考慮することができる。モデル化されたシステムは、周期的であり、動力学的により良好な、つまり排気脈動が推定される推定器を得ることができる。
【０１１７】
さらに、本発明によれば、計算時間を従来の方法に対して８０分の１に短縮することができる。
【図面の簡単な説明】
【０１１８】
【図１】排気プロセスの具体的な構成部材を模式的に示す図である。
【図２Ａ】４つのシリンダの１つについて、基準空燃比
【数２９】
2006336644.tif 000032
を時間（Ｔ）の関数として示した図と、本発明の推定器の結果
【数３０】
2006336644.tif 000033
を時間（Ｔ）の関数として示した図である。
【図２Ｂ】４つのシリンダの１つについて、基準空燃比
【数３１】
2006336644.tif 000034
を時間（Ｔ）の関数として示した図と、本発明の推定器の結果
【数３２】
2006336644.tif 000035
を時間（Ｔ）の関数として示した図である。
【図３】推定器の構造を示す図である。
【図４Ａ】４つのシリンダの１つについて、基準空燃比
【数３３】
2006336644.tif 000036
を時間（Ｔ）の関数として示した図と、本発明による遅延を考慮した推定器の結果
【数３４】
2006336644.tif 000037
を時間（Ｔ）の関数として示した図である。
【図４Ｂ】４つのシリンダの１つについて、基準空燃比
【数３５】
2006336644.tif 000038
を時間（Ｔ）の関数として示した図と、本発明による遅延を考慮した推定器の結果
【数３６】
2006336644.tif 000039
を時間（Ｔ）の関数として示した図である。
【符号の説明】
【０１１９】
ＣＥ    排気マニフォールド
Ｔ    ターボ過給器のタービン
ＳＲ    プローブ
ＤＳ１〜ＤＳ４    シリンダの出口での流量
λ1〜λ4    シリンダ内の空燃比
Ｎｅ    エンジンの回転数
α    クランク軸の角度
ＭＮＢＯ    開ループ質量モデル
Ｄ    遅れ時間
ＥＲＦＡ    空燃比の推定器
【図１】
2006336644.tif 000040
【図２Ａ】
2006336644.tif 000041
【図２Ｂ】
2006336644.tif 000042
【図３】
2006336644.tif 000043
【図４Ａ】
2006336644.tif 000044
【図４Ｂ】
2006336644.tif 000045