(57)【要約】グラフィックシグナチャーに基づいて内燃機関の試験用に有用な情報を抽出する方法。機関の運転に関する少なくとも１つの情報項目を含んでいる搬送信号が取得される。次に、各機関サイクルについて、この信号がグラフィックシグナチャーに変換され、信号の特徴の集合が、解析が容易なグラフの形態となる。次に、このシグナチャーの少なくとも１つの属性が求められ、これから、属性と情報項目との所定の関係を使用して関連情報が推定される。最後に、情報が機関の試験に使用される。内燃機関の試験に適用される。
2011501038.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の作動に関連した少なくとも１つの情報を含む少なくとも１つの信号を得るための少なくとも１つの検出器を備えた内燃機関の制御方法であって、該少なくとも１つの信号は離散的な計測値のセットからなる、内燃機関の制御方法において、
信号をグラフィックシグナチャーに変換可能な関数であって、該関数は、スライドする時間窓に含まれている信号計測値を、多次元空間からこれより小次元の空間に投影することによって、信号をグラフィックシグナチャーに変換可能である関数を選択するステップと、
各機関サイクルについて、
前記関数によって前記信号をグラフィックシグナチャーに変換するステップと、
前記情報を前記シグナチャーから抽出するステップと、
前記情報を使用して前記機関を制御するステップと、
を有することを特徴とする方法。
【請求項２】
前記情報は、
前記シグナチャーの少なくとも１つの属性を選択するステップと、
前記属性と前記情報との間の関係を求めるステップと、
各機関サイクルについて、
前記シグナチャーの前記属性の値を計算するステップと、
前記値と前記関係とによって前記情報を抽出するステップと、
を実行することによって抽出される、請求項１に記載の方法。
【請求項３】
前記これより小次元の空間は２次元である、請求項１または２に記載の方法。
【請求項４】
前記グラフィックシグナチャーは、任意の時刻ｔにおいて、
時刻ｔにおける計測値ｙ（ｔ）と時刻ｔ以前のＮ個の計測値とからなるベクトルＹm（ｔ）を構成するステップと、
前記ベクトルＹm（ｔ）を２次元平面内の点を表す対（ｙ１，ｙ２）に変換するステップと、
を実行することによって得られる、請求項１から３のいずれか１項に記載の方法。
【請求項５】
前記属性と前記情報との間の前記関係は、
様々な情報の値に対してグラフィックシグナチャーを構成するステップと、
各前記シグナチャーについて前記属性を計算するステップと、
各前記シグナチャーについて前記属性と前記情報の組を比較することによって前記関係を推定するステップと、
を機関テストベンチで実行することによって定められる、請求項２から４のいずれか１項に記載の方法。
【請求項６】
前記信号は、前記機関に備えられたコモンレール内の圧力計測値である、請求項１から５のいずれか１項に記載の方法。
【請求項７】
前記情報は噴射の検出に相当する、請求項６に記載の方法。
【請求項８】
前記属性は前記シグナチャーの表面積である、請求項７に記載の方法。
【請求項９】
前記噴射は、前記シグナチャーの前記表面積を所定のしきい値と比較することによって検出される、請求項８に記載の方法。
【請求項１０】
前記信号は瞬間機関回転数の計測値である、請求項１から５のいずれか１項に記載の方法。
【請求項１１】
前記情報は機関トルクである、請求項１０に記載の方法。
【請求項１２】
前記属性は前記グラフィックシグナチャーの横方向の直径及び縦方向の直径に依存している、請求項１０または１１に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関制御の分野に関する。本発明は特に、機関に配置されている検出器から取得された信号を解析し、機関制御に関連する情報を抽出することができる方法に関する。
【背景技術】
【０００２】
内燃機関はますます複雑化し、それらの制御はますます洗練されてきている。先進的な機関制御と診断システムは、機関内で発生する様々な事象（噴射、燃焼等）に関する詳細な情報を必要とする。これらの事象を監視し、それらに依存した機関パラメータを知ることによって、機関性能と排出物低減とを相当に改善することができる。
【０００３】
機関制御に関連する情報を取得する様々な方法が知られている。
【０００４】
たとえば、以下のような、シリンダ内圧力の再現方法が知られている（非特許文献１，２）。
【０００５】
燃焼パラメータの推定方法も知られている（特許文献１，非特許文献３）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】J．，Chauvin、Y．，Bentolila、及びO．，Grondin：（2006）、「Methode d'estimation de parametres de combustion a partir de signaux vibratoires（振動の予兆に基づく燃焼パラメータの評価方法）」、仏国特許出願第２００６／０２１１１号、Institut Francais du petrole
【非特許文献】
【０００７】
【非特許文献１】J．Antoni、J．Daniere、F．Guillet、「Effective vibration analysis of IC Engines using Cyclostationarity．Part I：A Methodology for Condition Monitoring（サイクル定常性を利用した内燃機関の効果的な振動解析  第１部：状態監視の方法）」、Journal of Sound and Vibration． Vol．257、No．5、2002年11月、815〜837ページ
【非特許文献２】J．Antoni、J．Daniere、F．Guillet、R.B. Randall、「Effective vibration analysis of IC Engines using Cyclostationarity．Part II：New Results on the Reconstruction of the Cylinder Pressure（サイクル定常性を利用した内燃機関の効果的な振動解析  第２部：シリンダ圧力の再現についての新しい結果）」、Journal of Sound and Vibration． Vol．257、No．5、2002年11月、839〜856ページ
【非特許文献３】「Rapid−Prototyping Multi−Sensors Processing Platform for Real Time Engine Control and Diagnosis（実時間機関制御と診断のためのマルチセンサ処理プラットフォームの高速なプロトタイピング）」、Olivier Grondin、Laurent Duval、Fabrice Guillemin、Stephan Ker、Gilles Corde、Christian Vigild、Fifth IFAC Symposium on Advances in Automotive Control Seascape Resort、Aptos、California、USA、2007年8月
【発明の概要】
【発明が解決しようとする課題】
【０００８】
これらの方法全てに固有の問題は、多量の情報を含んでいる信号の複雑さである。この情報のいくつかは機関の作動に直接的に関連し、いくつかは間接的に関連し、そしていくつかは攪乱、雑音または背景信号からなっている。
【０００９】
機関制御に関連した情報の取得は、大量の他の「寄生」情報から有用な情報、すなわち機関の作動に関連した情報を抽出するための、複雑な信号の解析を伴う。
【００１０】
機関制御には実時間の適用が必要なことを考慮すると、この解析技法は単純で、高速で、正確でなければならない。
【００１１】
本発明の目的は、内燃機関の制御のための有用な情報を抽出する代替の方法である。この方法は、様々な機関検出器から得られた高周波信号から生成されたグラフィックシグナチャーに基づいている。
【課題を解決するための手段】
【００１２】
本発明は、内燃機関の作動に関連した少なくとも１つの情報を含む少なくとも１つの信号を得るための少なくとも１つの検出器を備えた内燃機関の制御方法であって、該少なくとも１つの信号は離散的な計測値のセットからなる、内燃機関の制御方法に関する。この方法は、信号をグラフィックシグナチャーに変換可能な関数を選択するステップと、各機関サイクルについて、前記関数によって前記信号をグラフィックシグナチャーに変換するステップと、前記情報を前記シグナチャーから抽出するステップと、前記情報を使用して前記機関を制御するステップと、を有している。
【００１３】
情報は、前記シグナチャーの少なくとも１つの属性を選択するステップと、前記属性と前記情報との間の関係を求めるステップと、各機関サイクルについて、前記シグナチャーの前記属性の値を計算するステップと、前記値と前記関係とによって前記情報を抽出するステップと、を実行することによって抽出することができる。
【００１４】
信号が離散的な計測値のセットである場合、グラフィックシグナチャーは、スライドする時間窓に含まれている信号計測値を、多次元空間からこれより小次元、例えば２次元のの空間に投影可能である関数によって得ることができる。そのため、任意の時刻ｔにおいて、時刻ｔにおける計測値ｙ（ｔ）と時刻ｔ以前のＮ個の計測値とからなるベクトルＹm（ｔ）を構成するステップと、前記ベクトルＹm（ｔ）を２次元平面内の点を表す対（ｙ１，ｙ２）に変換するステップと、を実行することができる。
【００１５】
前記属性と前記情報との間の前記関係は、様々な情報の値に対してグラフィックシグナチャーを構成するステップと、各前記シグナチャーについて前記属性を計算するステップと、各前記シグナチャーについて前記属性と前記情報の組を比較することによって前記関係を推定するステップと、を機関テストベンチで実行することによって定めることができる。
【００１６】
一実施態様によれば、信号は機関に備えられたコモンレール内の圧力計測値である。抽出される情報は噴射の検出であってよい。この場合、シグナチャーの表面積を属性として使用することができる。次に、噴射が行われたかどうかが、シグナチャーの表面積を所定のしきい値と比較することによって検出される。
【００１７】
他の実施態様によれば、信号は瞬間機関回転数の計測値である。抽出される情報は機関トルクの推定値であってよい。この場合、グラフィックシグナチャーの横方向の直径及び縦方向の直径に基づく属性を使用することができる。
【００１８】
本発明の方法のその他の特徴及び利点とは、添付の図面を参照して、非限定的な例により説明する実施形態の以降の説明を読むことで明らかになろう。
【図面の簡単な説明】
【００１９】
【図１】グラフィックシグナチャーを使用して、高周波の計測値から機関情報を抽出する方法のステップを示す図である。
【図２】本発明のグラフィックシグナチャーの構成方法を示す図である。
【図３】グラフィックシグナチャーの構成中における中間点の定義の例である。
【図４】２回噴射の場合のレール内の実際の圧力計測値から得られたシグナチャーの例を示す図である。
【図５】様々なＭＩＰ値（Ｎ＝１５００ｒｐｍ）について、瞬間機関回転数の実際の計測値から得られたシグナチャーの例を示す図である。
【図６】ＭＩＰとグラフィックシグナチャーとの間の相関の例を示す図である。
【図７】図５のシグナチャーから抽出された属性を使用した、ＭＩＰのオンライン推定を示す図である。
【発明を実施するための形態】
【００２０】
図１は、検出器から得られた計測値から機関情報を抽出する方法を示している。本方法は４つのステップを有している。
１−検出器からの信号の取得（ＡＣＱ）
２−グラフィックシグナチャーの生成（ＳＩＧＮ）
３−グラフィックシグナチャーの解析による、相関のある属性の決定（ＡＮＡ）
４−情報の抽出（ＩＮＦ）。
【００２１】
１−検出器からの信号の取得（ＡＣＱ）
内燃機関（ＭＯＴ）には様々な検出器が備えられている。これらの検出器としては、例えば、シリンダ内に配置されている圧力検出器、機関シャフトに配置されている瞬間機関回転数検出器、コモンレールを備えた機関用のレール圧力検出器が考えられる。
【００２２】
好ましくは、高周波信号が（６クランク角度毎、１クランク角度毎などで）計測される。
【００２３】
これらの信号は、「寄生的」と呼ばれる他の情報の中に、機関作動に関する情報を含んでいる。その情報は、噴射の発生等の事象に関する情報である場合がある。その情報は、トルクなどの機関パラメータを知るための情報であることもある。よって、目的はこの情報を抽出することである。そのため、グラフィックシグナチャーが使用される。
【００２４】
２−グラフィックシグナチャーの構成（ＳＩＧＮ）
シグナチャーと呼ばれるものは、ある１つのものを他の１つのものに割り当てることができるようにする、認識可能な特徴のセットである。本発明の範囲において、それはあるものを内燃機関の作動に関連した特定の事象に割り当てることができるようにする、認識可能な特徴のセットである。
【００２５】
本発明によれば、グラフィックシグナチャーが構成される。グラフィックシグナチャーとは、認識可能な特徴がグラフの形態で表されるシグナチャーである。
【００２６】
このグラフは、幾何学形状を離散的に構成する一式の点からなっている。一例を図４に示す。同図では、図示のために連続的な線が点と点の間に付け加えられている。
【００２７】
これらのシグナチャーは、信号をグラフィックシグナチャーに変換可能な関数によって求められる。これを行うためのいくつかの技術は、専門家には公知である。
【００２８】
好ましい実施形態によれば、シグナチャーは、オンラインで（実時間で）得られスライドする時間窓内に含まれる計測値を、多次元空間からこれより小次元の空間（一例では２次元平面）に投影することを可能にする関数によって得られる。この次元数の減少によって、信号をより容易に解析することができる。
【００２９】
図２は、２次元のグラフィックシグナチャーを生成する方法を示している。シグナチャーは、１つの機関サイクルと少なくとも１つの信号（ｙ）とに関連付けられている。このようなグラフィックシグナチャーを生成するために、各サイクルで以下のステップが実行される。
【００３０】
−各時刻ｔで、現在の計測値ｙ（ｔ）とＮ個の過去の計測値とからなるベクトルＹm（ｔ）が生成される。
【００３１】
【数１】
2011501038.tif 000003
【００３２】
ベクトルと呼ばれるものは、ｎ次元空間において所与の順番に配置されているｎ個のスカラ量によって特徴付けられる量である。従って、この場合（Ｎ＋１）個の組み(uplet)となる。
【００３３】
整数Ｎは「シグナチャー次元」と呼ばれる。スライドする時間窓（ＦＴＧ）は時間間隔（ｔ−Ｎ．τ；ｔ）によって定められる。パラメータτは、信号ｙの、２つの計測値の間の時間間隔を表している。
【００３４】
−次に、ベクトルＹm（ｔ）は、シグナチャー平面と呼ばれる２次元平面（Ｙ1；Ｙ2）上の、１つの点を表す対（ｙ１，ｙ２）に変換される。この変換は以下の写像Ｐ(application)によって実行される。
【００３５】
【数２】
2011501038.tif 000004
【００３６】
−こうして、図２に示すように、時間の経過とともにベクトルＹm（ｔ），Ｙm（＋τ），．．．，Ｙm（ｔ＋４τ），．．．が順次得られる。写像Ｐは、それらの各々をシグナチャー平面内の点に関連付ける。こうして、２次元のグラフィックシグナチャーが得られる。
【００３７】
写像Ｐの定義
任意の時刻ｔにおいて、スライドする時間窓の中に含まれる計測値は、以下の計測値ベクトルの生成に使用されることを想起されたい。
【００３８】
【数３】
2011501038.tif 000005
【００３９】
シグナチャー平面内の点を各ベクトルＹm（ｔ）に関連付けるために、以下のステップが実行される。
【００４０】
ａ）まず過去の計測値ベクトルＹ（ｔ）が正規化され、以下の正規化ベクトルが得られる。
【００４１】
【数４】
2011501038.tif 000006
【００４２】
ここでεは一定の正規化定数である。
【００４３】
ｂ）次に、超立方体の各点［−１、＋１］Nを関連付ける写像が以下のように定義される。
【００４４】
【数５】
2011501038.tif 000007
【００４５】
ここで「image」は、この複素数のイメージポイント(point image)、つまり２次元平面においてそれに対応する点を指す。
【００４６】
図３は、正規化されたベクトルが以下のように与えられる特定の場合における、点Ψi（Ｙ）（注：原文はＹの上に横線がある。）の位置の例を示している。
【００４７】
【数６】
2011501038.tif 000008
【００４８】
ｃ）中間点Ψi（Ｙ）（注：原文はＹの上に横線がある。）は、次に、シグナチャー平面における以下の２つの点を計算するために使用される。
【００４９】
【数７】
2011501038.tif 000009
【００５０】
これらは、それまでに計算された中間点（それぞれ重み付けされていたり、されていなかったりする。）の重心である。そして、写像Ｐは以下のように与えられる。
【００５１】
【数８】
2011501038.tif 000010
【００５２】
ベクトル（Ｙ，ｙ）（注：原文はＹの上に横線がある。）は、最後のＮ個の要素に正規化手順を適用することによって、計測値ベクトルＹmから得られることに留意されたい。
【００５３】
このようにして、各計測値ベクトルＹm（ｔ）（時刻ｔにおいて、スライドする窓内に記録されている過去の計測値を使用して生成される。）をシグナチャー平面の点に関連付けることができる写像が完全に定義される。
【００５４】
結果の例
図４と５は、本発明の方法によって得られたグラフィックシグナチャーを示している。
【００５５】
図４は、２回噴射の場合にレール内の実際の圧力計測値から得られたシグナチャーを示している。
【００５６】
図５は、様々なＭＩＰ値（Ｎ＝１５００ｒｐｍ）について実際の瞬間機関回転数の計測値から得られたシグナチャーを示している。
【００５７】
３−グラフィックシグナチャーの解析による相関のある属性の決定（ＡＮＡ）
グラフィックシグナチャーは、検出器から送られてくる各信号について、サイクル毎にこのようにして生成される。従って、グラフィックシグナチャーは検出器によって与えられる情報を表現することになる。このグラフィックシグナチャーは、いくつかの点からなる形状（１つの点が１つの信号計測値を表し、実際には、各点はＮ＋１個の計測値の組に相当している。）を表している。この形状は、多くの属性を有している。
【００５８】
次に、機関作動（事象、機関パラメータ）に関する情報と相関のある属性が、グラフィックシグナチャーを特徴付けるそれらの属性の中から探索される。そのため、様々な試験が実施される機関テストベンチが使用され、
−グラフィックシグナチャーが情報の様々な値について生成され、
−情報の値によって変化する属性が求められる。
【００５９】
グラフィックシグナチャーから得られた形状が正確な幾何学形状を有している場合、たとえば円形の場合には、属性を直接計算することができる。つまり直径、表面積、周長、またはいくつかの属性の組み合わせを計算することができる。これらの属性は、シグナチャーを構成している点だけを使用して計算され、点を結ぶ曲線からは計算されないことが好ましい。
【００６０】
以降に例を示す。
【００６１】
４−機関制御に有用な情報の抽出
これらの相関のある属性によって、噴射などの事象を検出したり、ＭＩＰなどのパラメータを推定したりすることができる。
【００６２】
このステップは、機関制御のための有用な情報を求めるために、グラフィックシグナチャーの属性を活用することからなる。
【００６３】
情報は直接的に抽出することが可能である。すなわち、属性の値によってはその値から直接的に抽出することができ、または、属性の値を所定のしきい値と比較することによって抽出することができる。あるいは、属性の値から情報を推定できるようにあらかじめ関係を定めておくことによって、情報を間接的に抽出することができる。
【００６４】
１つの属性とある１つの情報との間の関係を求めるために、以下のステップを機関テストベンチ上で実行することができる。
−様々な情報の値に対してグラフィックシグナチャーを構成するステップ
−各シグナチャーについて属性を計算するステップ
−各シグナチャーについて属性と情報の組を比較することによって関係を推定するステップ。
【００６５】
提案している方法は、従来の機関制御機構内で容易に実施でき、実時間で実行することができる。
【００６６】
情報抽出方法の適用
Ｉ−レール内の圧力計測値からのパイロット噴射の検出
コモンレール噴射装置は、様々な噴射の仕方（マルチ噴射(multi-injection)）によって必要な量の燃料を生成できるようにする高圧噴射装置である。短いパイロット噴射が主噴射に先行する。このパイロット噴射は、特に冷間始動条件下で、燃焼ノイズを減少させるために使用される。パイロット噴射は持続時間が短いため、常に行われるわけではない。ある条件下では、噴射ノズルは制御されるが、燃料は噴射されない。噴射が行われないと、機関制御が影響を受ける。
【００６７】
そのため、パイロット噴射を検出することによって、機関の性能を改善することができる。
【００６８】
この状況において、レール内の圧力計測値から生成されたグラフィックシグナチャーが、パイロット噴射の有無を検出するために用いられる。
【００６９】
図４は、本発明によってレール内の圧力計測値から得られたグラフィックシグナチャーの例を示している。比較的円形に近い形状を各々が構成している多数の点の集合が見られる。機関テストベンチ上で主噴射（ＰＲＩ）は最大の円に相当することがわかる。パイロット噴射（ＰＩＬ）は、円が点線の円（ＳＥＵ）よりも大きい時に生じている。この円（ＳＥＵ）がしきい値である。しきい値は機関テストベンチ上で定められ、そしてグラフィックシグナチャーとの比較によって、パイロット噴射を検出するようにされる。
【００７０】
具体的には、機関の作動中に、シグナチャーがサイクル毎に計算され、（解析的に表現可能な）しきい値と比較される。たとえば、グラフィックシグナチャーの表面積がしきい値の表面積よりも大きいとき、パイロット噴射が発生したとみなすことができる。
【００７１】
ＩＩ−瞬間機関回転数の計測値からの、機関トルクの推定
各シリンダ内で燃焼によって発生するトルクを知ることは、機関制御の重要な要素である。しかし、コストと実現性のために、現在の車輛にはそれを直接計測するのに適した検出器は備えられていない。他方、機関トルクの情報をもたらす瞬間機関回転数は計測することが可能である。
【００７２】
測定可能な瞬間機関回転数の計測値を使用して、これに対応するトルクが推定される。瞬間機関回転数から機関トルクを推定する方法は、文献よりいくつか知られている。たとえば、周波数領域での逆畳み込み(deconvolution)手法を使用することができる。
【００７３】
本発明によれば、グラフィックシグナチャーに基づいて、各シリンダが発生するトルクについての定量的な情報を得ることができる。シグナチャーは瞬間機関回転数の計測値から生成され、トルクの情報を含んでいるため、トルクに対して敏感である。
【００７４】
従って、得られたシグナチャーとトルクの値との間の相関を求め、トルクの推定に有用な属性をシグナチャーから抽出することができる。
【００７５】
図５は、指示平均圧力（ＭＩＰ）の様々な値について、実際の瞬間機関回転数の計測値から得られたシグナチャーを示している。このシグナチャーによって、瞬間機関回転数から、ＭＩＰについての定量的な情報を、そしてその結果、各シリンダが発生するトルクについての定量的な情報を得ることができる。実際、トルクとＭＩＰとは、以下の関係で結ばれている。
【００７６】
【数９】
2011501038.tif 000011
【００７７】
ここで、ＭＩＰ  ＰＭＩはバールで表され、トルクＣはＮ・ｍで表され、シリンダ容積Ｖcylはリットルで表される。
【００７８】
様々なシグナチャーの点によって構成されている楕円の大きさがＭＩＰと相関があることが分かっている。実際、ＭＩＰが大きくなると、楕円は大きくなる。
【００７９】
次に、以下のシグナチャーの属性ＡＴＲ、すなわち横方向の直径＋縦方向の直径が定められる。この属性は、シグナチャーの横方向の両端２点の横座標の差と、シグナチャーの縦方向の両端２点の縦座標の差と、の和を求めることによって計算される。
【００８０】
横方向の直径＝ｍａｘ（ｘｉ）−ｍｉｎ（ｘｉ）
縦方向の直径＝ｍａｘ（ｙｉ）−ｍｉｎ（ｙｉ）
ＡＴＲ＝ｍａｘ（ｘｉ）−ｍｉｎ（ｘｉ）＋ｍａｘ（ｙｉ）−ｍｉｎ（ｙｉ）。
【００８１】
図６は、ＭＩＰとグラフィックシグナチャーの属性ＡＴＲとの間の相関を示している。連続している線（ＲＥＬ）は、ＭＩＰと図５のシグナチャーから抽出された属性ＡＴＲとの間の関係の推定値である。
【００８２】
従って、走行中の車輛の機関を制御する状況においては、グラフィックシグナチャーが実際の瞬間機関回転数の計測値から計算され、次に属性ＡＴＲが計算され、関係ＲＥＬがＭＩＰを推定するために適用される。
【００８３】
グラフィックシグナチャーに基づく方法は、最適化またはフィルタ処理のための複雑なアルゴリズムと比べてシグナチャーの計算コストが低い（単純な算術演算である）限り、容易に実時間で適用することができる。
【００８４】
図７は、図５のシグナチャーから実時間で抽出された属性を使用して得られたオンラインＭＩＰ推定の結果を示している。薄色の曲線はＭＩＰの実際の値を示しており、濃色の曲線は推定値を示している。
【図１】
2011501038.tif 000012
【図２】
2011501038.tif 000013
【図３】
2011501038.tif 000014
【図４】
2011501038.tif 000015
【図５】
2011501038.tif 000016
【図６】
2011501038.tif 000017
【図７】
2011501038.tif 000018
【国際調査報告】
2011501038.tif 000019
2011501038.tif 000020
2011501038.tif 000021
2011501038.tif 000022
2011501038.tif 000023
2011501038.tif 000024