特許 6835317 副室式エンジンの燃焼状態判別方法、副室式エンジン、及びエンジンシステム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6835317

(24)【登録日】2021年2月8日

(45)【発行日】2021年2月24日

(54)【発明の名称】副室式エンジンの燃焼状態判別方法、副室式エンジン、及びエンジンシステム

(51)【国際特許分類】

   F02D 45/00 20060101AFI20210215BHJP

【ＦＩ】

   F02D45/00 368A

   F02D45/00 368S

【請求項の数】7

【全頁数】19

(21)【出願番号】特願2017-233730(P2017-233730)

(22)【出願日】2017年12月5日

(65)【公開番号】特開2019-100282(P2019-100282A)

(43)【公開日】2019年6月24日

【審査請求日】2020年2月21日

(73)【特許権者】

【識別番号】000000284

【氏名又は名称】大阪瓦斯株式会社

(73)【特許権者】

【識別番号】504137912

【氏名又は名称】国立大学法人 東京大学

(74)【代理人】

【識別番号】110001818

【氏名又は名称】特許業務法人Ｒ＆Ｃ

(72)【発明者】

【氏名】▲高▼島 良胤

(72)【発明者】

【氏名】佐古 孝弘

(72)【発明者】

【氏名】山▲崎▼ 由大

(72)【発明者】

【氏名】高居 明弘

(72)【発明者】

【氏名】金子 成彦

【審査官】 沼生 泰伸

(56)【参考文献】

【文献】 特開２００７−１７０４０５（ＪＰ，Ａ）

【文献】 特開２００４−３５３５３１（ＪＰ，Ａ）

【文献】 特開２０１２−２２５２６８（ＪＰ，Ａ）

【文献】 米国特許出願公開第２０１５／０３１６０２０（ＵＳ，Ａ１）

【文献】 中国特許出願公開第１０５２７５６１５（ＣＮ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｆ０２Ｄ ４５／００

Ｇ０１Ｍ １５／１２

(57)【特許請求の範囲】

【請求項1】

シリンダとピストンとにより包囲される主燃焼室と、前記主燃焼室に連通孔を介して連通する副室とを備え、前記副室にて形成される火炎ジェットを前記連通孔を介して前記副室から前記主燃焼室へ噴射させて、前記主燃焼室の内部の混合気を燃焼させる副室式エンジンの燃焼状態判別方法であって、
前記主燃焼室の筒内圧力から導出される複数の第１指標に基づき、筒内圧力検出手段にて検出される１サイクル毎の前記主燃焼室の筒内圧力データから低周波成分及び高周波成分をカットしたフィルタ処理後データを分析するクラスタ分析により、１サイクル毎の前記主燃焼室の燃焼状態を、通常燃焼状態が分類される通常燃焼群と、ノッキングが発生するノッキング燃焼状態とノッキングとは別の異常燃焼状態とを含む非通常燃焼が分類される非通常燃焼群とに分類するクラスタ分析工程と、
前記非通常燃焼群から、前記フィルタ処理後データの閾値である教師データ閾値に基づいて、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを抽出する教師データ抽出工程と、
前記第１教師データと前記第２教師データとを用い、前記主燃焼室の筒内圧力から導出される複数の第２指標に基づいて前記フィルタ処理後データを分析する判別分析により、前記非通常燃焼群を、前記ノッキング燃焼状態が分類されるノッキング燃焼群及び前記異常燃焼状態が分類される異常燃焼群の何れかに分類する判別分析工程とを含む副室式エンジンの燃焼状態判別方法。

【請求項2】

前記第１指標は、前記フィルタ処理後データのＲＭＳ値と、前記フィルタ処理後データの最大値から最小値を減算したΔＰ値とを含むものである請求項１に記載の副室式エンジンの燃焼状態判別方法。

【請求項3】

前記教師データ抽出工程では、
前記フィルタ処理後データが初めて所定の圧力振動判別閾値以上となる圧力振動開始時期のデータを前記教師データ閾値として、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを絞り込む第１絞込工程と、
前記フィルタ処理後データの絶対値の最大値を前記教師データ閾値として、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを絞り込む第２絞込工程とを記載の順に実行する請求項１又は２に記載の副室式エンジンの燃焼状態判別方法。

【請求項4】

前記第２指標は、
前記フィルタ処理後データのＲＭＳ値を、前記フィルタ処理後データの最大値から最小値を減算したΔＰ値で除算したＲＭＳ／ΔＰ値と、
１サイクルにおいて、前記副室の内部の混合気に点火する点火手段の点火時期から判別時期までの圧力振幅積分値に対する、前記点火手段の点火時期から前記フィルタ処理後データが最大となる時期までの圧力振幅積分値の割合の平方根である二乗和平方根割合とを含む請求項１〜３の何れか一項に記載の副室式エンジンの燃焼状態判別方法。

【請求項5】

ネットワーク回線を介して複数の前記副室式エンジンの夫々から取得した１サイクル毎の前記主燃焼室の筒内圧力データを用いて、前記クラスタ分析工程と、前記教師データ抽出工程と、前記判別分析工程とを実行する請求項１〜４の何れか一項に記載の副室式エンジンの燃焼状態判別方法。

【請求項6】

シリンダとピストンとにより包囲される主燃焼室と、前記主燃焼室に連通孔を介して連通する副室とを備え、前記副室にて形成される火炎ジェットを前記連通孔を介して前記副室から前記主燃焼室へ噴射させて、前記主燃焼室の内部の混合気を燃焼させて、前記主燃焼室の内部の混合気を燃焼させたときの燃焼状態の判別を行う制御装置とを有する副室式エンジンであって、
前記制御装置は、前記主燃焼室の筒内圧力から導出される複数の第１指標に基づき、筒内圧力検出手段にて検出される１サイクル毎の前記主燃焼室の筒内圧力データから低周波成分及び高周波成分をカットしたフィルタ処理後データを分析するクラスタ分析により、１サイクル毎の前記主燃焼室の燃焼状態を、通常燃焼状態が分類される通常燃焼群と、ノッキングが発生するノッキング燃焼状態とノッキングとは別の異常燃焼状態とを含む非通常燃焼が分類される非通常燃焼群とに分類し、
前記非通常燃焼群から、前記フィルタ処理後データの閾値である教師データ閾値に基づいて、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを抽出し、
前記第１教師データと前記第２教師データとを用い、前記主燃焼室の筒内圧力から導出される複数の第２指標に基づいて前記フィルタ処理後データを分析する判別分析により、前記非通常燃焼群を、前記ノッキング燃焼状態が分類されるノッキング燃焼群及び前記異常燃焼状態が分類される異常燃焼群の何れかに分類する副室式エンジン。

【請求項7】

請求項６の副室式エンジンの複数で燃焼状態の判別を行うエンジンシステムであって、
前記制御装置は、ネットワーク回線を介して複数の前記副室式エンジンに接続され、前記副室式エンジンの夫々から１サイクル毎の前記主燃焼室の筒内圧力データを取得可能に構成されているエンジンシステム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、シリンダとピストンとにより包囲される主燃焼室と、前記主燃焼室に連通孔を介して連通する副室とを備え、前記副室にて形成される火炎ジェットを前記連通孔を介して前記副室から前記主燃焼室へ噴射させて、前記主燃焼室の内部の混合気を燃焼させる副室式エンジンの燃焼状態判別方法、及びそれを用いた副室式エンジン、エンジンシステムに関する。

【背景技術】

【0002】

コジェネレーションに用いられるガスエンジンの発電効率の向上技術の一つに、希薄燃焼方式がある。当該希薄燃焼方式とは、当量比が１以下の希薄状態下における燃焼を指し、低ＮＯｘ化と高効率化を同時に達成できるメリットがある。一方で、サイクル変動が大きく失火し易い等のデメリットも生じ、その解決のために副室が使用されている。副室とは、主燃焼室以外に、３〜５％程度の体積である予燃焼室を指すものである。
当該副室を備えた副室式エンジンでは、運転時に、副室へ量論混合比付近の点火し易い混合気を導入し、主燃焼室には当量比が低い混合気を導入した上で、副室内の混合気を火花点火させる。その後、副室内での燃焼により発生した既燃ガスが、急激に膨張し高温の火炎ジェットとして主燃焼室に流入することで、単室エンジンでは困難な主燃焼室内の希薄燃焼の安定化を実現している（特許文献１を参照）。
更に、エンジンでは、ノッキングが発生しているときの主燃焼室内圧力の振動成分を解析する場合、筒内圧力の経時データをフーリエ変換し、高周波成分と低周波成分とをカットするフィルタリング処理を実行した後に、筒内圧力のピーク値や、ＲＭＳ値が所定の閾値を超えているか否かにより、ノッキングが発生しているか否かを判別する技術が知られている（特許文献２を参照）。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２００１−２６３０６９号公報

【特許文献2】特開２００５−９０２５０号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

しかしながら、近年、副室式エンジンでは、希薄燃焼の安定化を図ることができるものであるが、副室式エンジンにおいて更なる高効率を求めて過給圧を高めると、主燃焼室内の筒内圧力値に、ノッキングとは異なる原因不明の振動成分が観測されることが判明している。
そして、当該原因不明の振動成分を伴う異常燃焼を、ノッキングを伴うノッキング燃焼と峻別して抽出することは、副室式エンジンの燃焼メカニズムを解明する上で重要であるが、上述したこれまで知られている周波数解析では難しく、新たな判別方法が望まれていた。

【0005】

本発明は、上述の課題に鑑みてなされたものであり、その目的は、ノッキング燃焼及び当該ノッキング燃焼とは異なる異常燃焼を峻別して、副室式エンジンにおける燃焼状態を判別できる副室式エンジンの燃焼状態判別方法、副室式エンジン、及びエンジンシステムを提供することにある。

【課題を解決するための手段】

【0006】

上記目的を達成するための副室式エンジンの燃焼状態判別方法は、
シリンダとピストンとにより包囲される主燃焼室と、前記主燃焼室に連通孔を介して連通する副室とを備え、前記副室にて形成される火炎ジェットを前記連通孔を介して前記副室から前記主燃焼室へ噴射させて、前記主燃焼室の内部の混合気を燃焼させる副室式エンジンの燃焼状態判別方法であって、その特徴構成は、
前記主燃焼室の筒内圧力から導出される複数の第１指標に基づき、筒内圧力検出手段にて検出される１サイクル毎の前記主燃焼室の筒内圧力データから低周波成分及び高周波成分をカットしたフィルタ処理後データを分析するクラスタ分析により、１サイクル毎の前記主燃焼室の燃焼状態を、通常燃焼状態が分類される通常燃焼群と、ノッキングが発生するノッキング燃焼状態とノッキングとは別の異常燃焼状態とを含む非通常燃焼が分類される非通常燃焼群とに分類するクラスタ分析工程と、
前記非通常燃焼群から、前記フィルタ処理後データの閾値である教師データ閾値に基づいて、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを抽出する教師データ抽出工程と、
前記第１教師データと前記第２教師データとを用い、前記主燃焼室の筒内圧力から導出される複数の第２指標に基づいて前記フィルタ処理後データを分析する判別分析により、前記非通常燃焼群を、前記ノッキング燃焼状態が分類されるノッキング燃焼群及び前記異常燃焼状態が分類される異常燃焼群の何れかに分類する判別分析工程とを含む点にある。

【0007】

上記目的を達成するための副室式エンジンは、
シリンダとピストンとにより包囲される主燃焼室と、前記主燃焼室に連通孔を介して連通する副室とを備え、前記副室にて形成される火炎ジェットを前記連通孔を介して前記副室から前記主燃焼室へ噴射させて、前記主燃焼室の内部の混合気を燃焼させて、前記主燃焼室の内部の混合気を燃焼させたときの燃焼状態の判別を行う制御装置とを有する副室式エンジンであって、その特徴構成は、
前記制御装置は、前記主燃焼室の筒内圧力から導出される複数の第１指標に基づき、筒内圧力検出手段にて検出される１サイクル毎の前記主燃焼室の筒内圧力データから低周波成分及び高周波成分をカットしたフィルタ処理後データを分析するクラスタ分析により、１サイクル毎の前記主燃焼室の燃焼状態を、通常燃焼状態が分類される通常燃焼群と、ノッキングが発生するノッキング燃焼状態とノッキングとは別の異常燃焼状態とを含む非通常燃焼が分類される非通常燃焼群とに分類し、
前記非通常燃焼群から、前記フィルタ処理後データの閾値である教師データ閾値に基づいて、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを抽出し、
前記第１教師データと前記第２教師データとを用い、前記主燃焼室の筒内圧力から導出される複数の第２指標に基づいて前記フィルタ処理後データを分析する判別分析により、前記非通常燃焼群を、前記ノッキング燃焼状態が分類されるノッキング燃焼群及び前記異常燃焼状態が分類される異常燃焼群の何れかに分類する点にある。

【0008】

上記特徴構成によれば、まず、筒内圧力の振動強度（筒内圧力の経時変化の振幅）により、通常燃焼群と非通常燃焼群とに容易に分類することができるクラスタ分析を採用している。クラスタ分析は、１サイクル毎の主燃焼室の筒内圧力データから低周波成分及び高周波成分をカットしたフィルタ処理後データを分析することにより、１サイクル毎の主燃焼室の燃焼状態を、通常燃焼状態が分類される通常燃焼群とそれ以外の燃焼群とに分類する。ここで、それ以外の燃焼群は、ノッキングが発生するノッキング燃焼状態と、ノッキングとは別の異常燃焼状態とを含む非通常燃焼が分類される非通常燃焼群である。
通常燃焼群と非通常燃焼群との分類は、筒内圧力の振動強度により、クラスタ分析のうち、特に、非階層的分類法の１つでありｎ個の対象を予め定められたｋ個（本発明の場合は、ｋ＝２）のクラスタに分ける分類手法であるｋ−Ｍｅａｎｓクラスタ分析を用いることにより、良好に実現できる。
さて、非通常燃焼群を成すノッキング燃焼状態と異常燃焼状態とは、特に、筒内圧力の振動強度（筒内圧力の経時変化の振幅）が共に大きく、両者を分類するためにはより精密な分析による分類が必要である。本発明にあっては、発明者が定めた教師データ閾値に基づいて抽出する教師データに基づいて、より精密な分析による分類が可能である判別分析を行うことで、非通常燃焼群が分類されるフィルタ処理後データを、ノッキング燃焼群と異常燃焼群とに分類する。
当該分類方法を用いることにより、まずもって、非通常燃焼群に分類されるフィルタ処理後データから、予め設定される教師データ閾値により、ノッキング燃焼状態としての第１教師データと、異常燃焼状態としての第２教師データとを、抽出する。
当該教師データ抽出工程において、予め設定される教師データ閾値により、非通常燃焼群としてのフィルタ処理後データが、第１教師データと第２教師データとに分類されるが、教師データ閾値を比較的厳しく設定するため、第１教師データと第２教師データの何れにも分類されない未分類のデータが多く残ることになる。
そこで、多変量解析の一手法として知られる判別分析により、当該第１教師データと第２教師データとを用いて、非通常燃焼群としてのフィルタ処理後データを未分類のデータを出すことなく、適切に分類することができる。尚、当該判別分析では、分類することにより、相関比が最大となる線形判別式を導出し、当該線形判別式により判別する。
以上より、ノッキング燃焼及び当該ノッキング燃焼とは異なる異常燃焼を峻別して、副室式エンジンにおける燃焼状態を判別できる副室式エンジンの燃焼状態判別方法、及び副室式エンジンを実現できる。

【0009】

副室式エンジンの燃焼状態判別方法の更なる特徴構成は、
前記第１指標は、前記フィルタ処理後データのＲＭＳ値と、前記フィルタ処理後データの最大値から最小値を減算したΔＰ値とを含むものである点にある。

【0010】

クラスタ分析を用いる場合、各クラスタの重心位置から所属群が決定されるため、散布図における各軸（複数の第１指標）の設定が重要になる。即ち、筒内圧力の振動の大小により各クラスタを分類するためには、筒内圧力の振動強度を示すパラメータを、各軸（複数の第１指標）として設定する必要がある。また、散布図上のユークリッド距離を用いるため、軸として用いる指標は、質的変数ではなく量的変数である必要がある。
以上の点を鑑みて、発明者らは、第１指標として、フィルタ処理後データのＲＭＳ値と、フィルタ処理後データの最大値から最小値を減算したΔＰ値とを設定している。

【0011】

副室式エンジンの燃焼状態判別方法の更なる特徴構成は、
前記教師データ抽出工程では、
前記フィルタ処理後データが初めて所定の圧力振動判別閾値以上となる圧力振動開始時期のデータを前記教師データ閾値として、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを絞り込む第１絞込工程と、
前記フィルタ処理後データの絶対値の最大値を前記教師データ閾値として、前記ノッキング燃焼状態としての第１教師データと、前記異常燃焼状態としての第２教師データとを絞り込む第２絞込工程とを記載の順に実行する点にある。

【0012】

ノッキング燃焼状態と異常燃焼状態とを判別する有効な指標としては、例えば、ＲＭＳ／ΔＰ値等が挙げられるが、当該値では閾値による両者の判別が難しく、教師データを作成する際の指標としては適していない。
そこで、上記特徴構成では、圧力振動開始時期のデータを教師データ閾値とする絞込みと、フィルタ処理後データの絶対値の最大値を教師データ閾値とする絞り込みとを記載の順に実行することで、教師データを良好に抽出できる。

【0013】

副室式エンジンの燃焼状態判別方法の更なる特徴構成は、
前記第２指標は、
前記フィルタ処理後データのＲＭＳ値を、前記フィルタ処理後データの最大値から最小値を減算したΔＰ値で除算したＲＭＳ／ΔＰ値と、
１サイクルにおいて、前記副室の内部の混合気に点火する点火手段の点火時期から判別時期までの圧力振幅積分値に対する、前記点火手段の点火時期から前記主燃焼室の筒内圧力データが最大となる時期までの圧力振幅積分値の割合の平方根である二乗和平方根割合とを含む点にある。

【0014】

さて、上記特徴構成の如く、ＲＭＳ／ΔＰ値と、二乗和平方根割合とを、判別分析の第２指標として採用した理由を以下に示す。
本発明の発明者らは、判別分析時に用いる指標を決定するため、異常燃焼状態のサイクル波形とノッキング燃焼状態のサイクル波形の特徴を把握するべく、複数のＲＭＳ値毎において、両者のサイクル波形を比較した。図６、７、８において、ＲＭＳ値が１６９ｋＰａ，１０７ｋＰａ，７５ｋＰａ付近となる波形を示す。因みに、図６、７、８で、上段は、筒内圧力データの経時変化を示すグラフ図、中段は、異常燃焼状態のフィルタ処理後データを示すグラフ図、下段は、ノッキング燃焼状態のフィルタ処理後データを示すグラフ図である。

【0015】

因みに、当該データを取得したときの副室式エンジンの仕様及び運転条件は下記に示すとおりである。
副室式エンジン１００は、ボア径１６５ｍｍ、ストローク２１５ｍｍ、排気量４．６Ｌの、都市ガス１３Ａを燃料とする単気筒副室式エンジンを使用した。また、圧縮比１４．５、回転数１２００ｒｐｍ、連通孔数６個で一定とした。
操作パラメータは、主燃焼室６に対する副室８の体積比率（２．３５、２．９、３．４７）、トーチ強度Ｉ（０．８、１．３１、１．８８）、連通孔狭角φ（１２０、１４０、１６０ｄｅｇ）、主燃焼室６と副室８との合計空気過剰率λ（１．８〜２．５）、ＩＭＥＰ（１．４３、２．０、２．４ＭＰａ）、副室点火時期（−１９．２〜−２．０）の６つとした。なお、トーチ強度Ｉは副室体積Ｖ_ｐｒｅ、連通孔数Ｎ、連通孔径Φ、主燃焼室ボア径Ｒを用いて、以下の〔数１〕で表される。また、計測は測定刻み幅（以下、ステップと呼ぶ場合がある）０．２ｄｅｇ、測定範囲−３６０〜３６０ｄｅｇＡＴＤＣ、１条件あたりの計測サイクル数２００サイクルとし、計３３８条件、６７、６００サイクル分を解析対象とした。

【0016】

【数1】

【0017】

図６、７、８の比較から、ＲＭＳ値の大小に関わらず、異常燃焼状態及びノッキング燃焼状態の筒内圧力データ及びフィルタ処理後データを示す波形には、以下の差異がある。
１点目は、筒内圧力の圧力振動の開始時期である。異常燃焼状態では、副室点火時期（θ＝０ｄｅｇ）の直後に圧力振動が開始するのに対して、ノッキング燃焼状態では、筒内圧力のピーク付近で圧力振動が開始している。
２点目は、圧力振動の振幅である。フィルタ処理後データの比較から、異常燃焼状態よりもノッキング燃焼状態の方が、圧力最大振幅ΔＰは大きい傾向にあることが確認される。
３点目は、圧力振幅の推移である。フィルタ処理後データから、異常燃焼状態では、筒内圧力のピークまでは圧力振動の振幅が略一定で推移し、筒内圧力のピークを過ぎて徐々に振幅が減衰する傾向である。これに対し、ノッキング燃焼状態では、筒内圧力の振動開始直後に振幅が最大となり、その後は単調に減少し、異常燃焼状態のように振幅が一定となる期間は観測されない。

【0018】

これらの比較結果から、異常燃焼状態とノッキング燃焼状態とを判別するには、以下の要件を満たす指標が有効である。
第１要件は、フィルタ処理後データの圧力振幅の推移から、両者を区別できることである。
即ち、異常燃焼状態では振動開始時期が早く、圧力振幅は振動開始から筒内圧力のピークにかけて略一定の値を保ちながら推移した後、筒内圧力のピークを過ぎると減少する。一方、ノッキング燃焼状態では、振動開始時期が筒内圧力のピーク付近であり異常燃焼状態に比べて遅く、圧力振幅は振動開始直後に最大に達した後、単調減少する。これらの各燃焼の圧力振幅の推移の特徴を表すことができる指標は、両者を判別できる指標として有効である。
上記第１要件を満たす指標として、フィルタ処理後データの圧力が所定の閾値に初めて達した時のクランク角度としての圧力振動開始時期が考えられる。しかしながら、筒内圧力の経時変化は、サイクル変動を伴うため、異常燃焼状態とノッキング燃焼状態との両燃焼を区別するための振動開始時期を定義する適切な閾値を設定することは難しい。一方、両燃焼の振動開始時期が異なることは、１サイクルの圧力振動成分の積分値に対する筒内圧力がピークに達するまでの圧力振動成分の積分値の割合の大小に対応するため、当該割合を指標として用いることでより適切に、両燃焼を峻別することができる。
そこで、当該指標として、１サイクルにおいて、副室の内部の混合気に点火する点火手段の点火時期から判別時期までの圧力振幅積分値に対する、点火手段の点火時期から主燃焼室の筒内圧力データが最大となる時期までの圧力振幅積分値の割合の平方根である二乗和平方根割合を採用した。

【0019】

第２要件は、フィルタ処理後のＲＭＳ値とΔＰ値の大きさの比の点から、両グループを区別できることである。
上述したように、ノッキング燃焼状態の方が異常燃焼状態よりもΔＰ値が大きくなる傾向がみられる。また、振動開始時期は、異常燃焼状態では早く、ノッキング燃焼状態では遅くなることからも、ＲＭＳ値が一定の条件では、異常燃焼状態のΔＰ値は、ノッキング燃焼状態のΔＰ値に比べ相対的に小さくなることが予想される。このことから、ＲＭＳ値とΔＰ値の大きさの比は両者を判別するための指標として有効である。
このことから、本発明にあっては、フィルタ処理後データのＲＭＳ値を、フィルタ処理後データの最大値から最小値を減算したΔＰ値で除算したＲＭＳ／ΔＰ値を、指標として採用した。

【0020】

副室式エンジンの燃焼状態判別方法の更なる特徴構成は、
ネットワーク回線を介して複数の前記副室式エンジンの夫々から取得した１サイクル毎の前記主燃焼室の筒内圧力データを用いて、前記クラスタ分析工程と、前記教師データ抽出工程と、前記判別分析工程とを実行する点にある。

【0021】

以上の特徴構成を有する燃焼状態判別方法は、単一の副室式エンジンから得られる筒内圧力のデータのみならず、複数の副室式エンジンから得られる筒内圧力のデータに基づいて分析を実行することができるから、例えば、当該分析結果に基づいて、より精度の高いエンジン制御を実行できる。

【0022】

上記目的を達成するためのエンジンシステムは、
これまで説明してきた副室式エンジンの複数で燃焼状態の判別を行うエンジンシステムであって、その特徴構成は、
前記制御装置は、ネットワーク回線を介して複数の前記副室式エンジンに接続され、前記副室式エンジンの夫々から１サイクル毎の前記主燃焼室の筒内圧力データを取得可能に構成されている点にある。

【0023】

以上のように構成することで、制御装置では、単一の副室式エンジンから得られる筒内圧力のデータのみならず、複数の副室式エンジンから得られる筒内圧力のデータに基づいて分析を実行することができるから、例えば、当該分析結果に基づいて、より精度の高いエンジン制御を実行できる。

【図面の簡単な説明】

【0024】

【図1】実施形態に係る副室式エンジンの概略構成図

【図2】燃焼状態判別の制御フロー図

【図3】ΔＰ値の定義を示すグラフ図

【図4】ＲＭＳ値の定義を示すグラフ図

【図5】教師データを抽出するためのフロー図

【図6】ＲＭＳ値が高い場合のグラフ図であり、上段は、筒内圧力データの経時変化を示すグラフ図、中段は、異常燃焼状態のフィルタ処理後データを示すグラフ図、下段は、ノッキング燃焼状態のフィルタ処理後データを示すグラフ図

【図7】ＲＭＳ値が中程度の場合のグラフ図であり、上段は、筒内圧力データの経時変化を示すグラフ図、中段は、異常燃焼状態のフィルタ処理後データを示すグラフ図、下段は、ノッキング燃焼状態のフィルタ処理後データを示すグラフ図

【図8】ＲＭＳ値が低い場合のグラフ図であり、上段は、筒内圧力データの経時変化を示すグラフ図、中段は、異常燃焼状態のフィルタ処理後データを示すグラフ図、下段は、ノッキング燃焼状態のフィルタ処理後データを示すグラフ図

【図9】図５の＃０２のステップでＹＥＳと判別されたデータを重畳させたグラフ図

【図10】図５の＃０３のステップでＹＥＳと判別されたデータを重畳させたグラフ図

【図11】二乗和平方根割合の定義を説明するためのグラフ図

【図12】異常燃焼状態としての教師データと、ノッキング燃焼状態としての教師データと、教師データ以外のデータとを重畳表示したグラフ図

【図13】通常燃焼群以外のデータに関し、判別分析による線形判別関数により、異常燃焼群とノッキング燃焼群とに判別したグラフ図

【発明を実施するための形態】

【0025】

本発明の実施形態に係る副室式エンジンの燃焼状態判別方法、副室式エンジン、エンジンシステムは、ノッキング燃焼及び当該ノッキング燃焼とは異なる異常燃焼を峻別して、副室式エンジンにおける燃焼状態を判別できるものである。
以下、図面に基づいて、その構成及び制御について説明する。

【0026】

図１に示すように、副室式エンジン１００は、シリンダ１７とピストン５とにより包囲される主燃焼室６と当該主燃焼室６に連通孔７を介して連通する副室８とを燃焼室４として備え、副室８にて形成される火炎ジェットを連通孔７を介して副室８から主燃焼室６へ噴出させて、主燃焼室６の内部の混合気Ｍを燃焼させるものであり、燃料ガスＧとして、気体燃料である都市ガス（１３Ａ）を利用する。
副室式エンジン１００は、主燃焼室６の吸気路９に燃料ガスＧを供給するための燃料ガス流量調整弁４１及び燃料ガスＧと空気とを混合するミキサ４０と、吸気路９に設けられた吸気弁１０と、主燃焼室６の排気路１１に設けられた排気弁１２と、副室８に燃料ガスＧを供給する副室燃料ガス通路１３と、副室８に形成された混合気Ｍを火花点火する点火プラグ１４と、燃料ガスＧの供給圧力を調整可能な調圧弁３３と、副室燃料ガス通路１３において調圧弁３３の下流に設けられ、燃料ガスＧの噴射量及び噴射タイミングを調整可能な燃料ガス噴射調整装置３４と、当該燃料ガス噴射調整装置３４の下流に設けられる逆止弁３５とを備えている。そして、吸気路９には吸気圧センサＳ２が設けられ、副室８には副室圧センサＳ１が設けられ、主燃焼室６には筒内圧センサＳ６が設けられている。

【0027】

副室式エンジン１００では、ピストン５がシリンダ１７内で往復運動し、それと共に吸気弁１０及び排気弁１２を開閉動作し、主燃焼室６において吸気、圧縮、膨張（燃焼）、排気の各行程が行われる。ピストン５の往復運動が、連結棒１８によって出力軸１９の回転運動として出力される。更に、副室式エンジン１００の出力軸１９に連結駆動されて発電を行うと共に、商用電力系統（図示せず）に連系可能に構成された同期発電機Ｅ（発電機に相当）が備えられている。

【0028】

副室式エンジン１００は、吸気行程において吸気弁１０を開弁状態として、燃料ガス流量調整弁４１およびミキサ４０によって生成された空気と燃料ガスＧとの混合気Ｍが吸気路９から主燃焼室６に供給されると共に、副室８に燃料ガスＧを供給し、圧縮行程にて連通孔７を通して主燃焼室６から副室８に混合気Ｍを流入させて、副室８内に直接供給された燃料ガスＧと混合気Ｍとの副室混合気を形成し、副室８の点火プラグ１４での火花点火によって副室８内で燃焼させた副室混合気を、連通孔７を介して主燃焼室６に火炎ジェットとして噴射するように構成されている。これにより、副室式エンジン１００は、主燃焼室６に供給される空気と燃料ガスＧとの混合気Ｍが、燃料希薄の場合においても良好に燃焼させることを可能にしている。

【0029】

シリンダヘッド２０に形成された副室８を形成する円柱状の凹部の上方開口部には、その上方開口部に嵌合する形態で有底筒状の口金２１が設けられている。口金２１には、副室８に連通する副室燃料ガス通路１３が形成されていると共に、口金２１の底部には、副室燃料ガス通路１３の通路端である燃料供給口１３ａが形成されている。
また、口金２１には、ピストン５の摺動方向に沿って配設される点火プラグ１４が、その点火点を副室８に位置する状態で、嵌合されている。

【0030】

副室式エンジン１００では、吸気工程のタイミングで燃料ガス噴射調整装置３４が動作することで、副室８へ燃料ガスＧが供給される。

【0031】

副室燃料ガス通路１３において燃料ガス噴射調整装置３４よりも上流側に設けられる調圧弁３３が備えられ、当該調圧弁３３により、副室燃料ガス通路１３における燃料ガスＧの圧力が調整される。従って、調圧弁３３により調整される副室燃料ガス通路１３の圧力の調整により、副室燃料ガス通路１３から副室８への燃料ガスＧの供給量が所望の供給量となる。

【0032】

副室式エンジン１００には、主燃焼室６に供給する燃料ガスＧの流量を調整する燃料ガス流量調整弁４１や、副室８へ供給する燃料ガスＧの流量を調整する調圧弁３３及び燃料ガス噴射調整装置３４や、副室８における点火プラグ１４の作動等を制御して副室式エンジン１００を制御するエンジン・コントロール・ユニットとしての制御装置２７が設けられている。当該制御装置２７は、ＣＰＵやＧＰＵ等から成る演算装置及びＲＡＭ等から成る記憶部を含むハードウェア群と、ソフトウェア群とが協働する形で、各種制御及び通信を実行可能に構成されている。
当該制御装置２７には、出力軸１９の回転角を検出するクランク角センサＳ３、エンジンの回転数を検出する回転数センサＳ５、酸素センサＳ４、副室圧センサＳ１、吸気圧センサＳ２、主燃焼室６の内部の筒内圧力を検出する筒内圧センサＳ６（筒内圧力検出手段の一例）等の検出信号が入力されるように構成されている。なお、クランク角センサＳ３の検出情報に基づいて、副室燃料ガス通路１３の燃料ガス噴射調整装置３４を噴射状態に切り換える期間を所望の期間とすると共に、点火プラグ１４を所望のタイミングで作動させて火花点火するようにしている。

【0033】

特に、当該実施形態に係る制御装置２７は、燃焼室４の燃焼状態をより適切に判別するべく、以下のように構成されている。
即ち、制御装置２７は、主燃焼室６の筒内圧力から導出される複数の第１指標に基づき、筒内圧センサＳ６にて検出される１サイクル毎の主燃焼室６の筒内圧力データＰ（ｉ）から低周波成分及び高周波成分をカットしたフィルタ処理後データを分析するクラスタ分析により、１サイクル毎の主燃焼室６の燃焼状態を、通常燃焼状態が分類される通常燃焼群と、ノッキングが発生するノッキング燃焼状態とノッキングとは別の異常燃焼状態とを含む非通常燃焼が分類される非通常燃焼群とに分類するよう構成されている。
更に、制御装置２７は、非通常燃焼群から、フィルタ処理後データＰ’（ｉ）の閾値である教師データ閾値に基づいて、ノッキング燃焼状態としての第１教師データと、異常燃焼状態としての第２教師データとを抽出するように構成されると共に、第１教師データと第２教師データとを用い、主燃焼室６の筒内圧力から導出される複数の第２指標に基づいてフィルタ処理後データを分析する判別分析により、非通常燃焼群を、ノッキング燃焼状態が分類されるノッキング燃焼群及び異常燃焼状態が分類される異常燃焼群の何れかに分類するように構成されている。

【0034】

以下、具体的な筒内圧力データＰ（ｉ）又はフィルタ処理後データＰ’（ｉ）に基づいて説明する。尚、以下に、筒内圧力データＰ（ｉ）の取得に使用した副室式エンジン１００の仕様及び運転条件を示す。以下、断りのない限り、図面及び明細書に示す筒内圧力データＰ（ｉ）及びフィルタ処理後データＰ’（ｉ）は、以下の仕様及び運転条件にて取得したものとする。
副室式エンジン１００は、ボア径１６５ｍｍ、ストローク２１５ｍｍ、排気量４．６Ｌの、都市ガス１３Ａを燃料とする単気筒副室式エンジンを使用した。また、圧縮比１４．５、回転数１２００ｒｐｍ、連通孔数６個で一定とした。
操作パラメータは、主燃焼室６に対する副室８の体積比率（２．３５、２．９、３．４７）、トーチ強度Ｉ（０．８、１．３１、１．８８）、連通孔狭角φ（１２０、１４０、１６０ｄｅｇ）、主燃焼室６と副室８との合計空気過剰率λ（１．８〜２．５）、ＩＭＥＰ（１．４３、２．０、２．４ＭＰａ）、副室点火時期（−１９．２〜−２．０）の６つとした。なお、トーチ強度Ｉは副室体積Ｖ_ｐｒｅ、連通孔数Ｎ、連通孔径Φ、主燃焼室ボア径Ｒを用いて、上述の〔数１〕で表される。また、計測は測定刻み幅（以下、ステップと呼ぶ場合がある）０．２ｄｅｇ、測定範囲−３６０〜３６０ｄｅｇＡＴＤＣ、１条件あたりの計測サイクル数２００サイクルとし、計３３８条件、６７、６００サイクル分を解析対象とした。

【0035】

図２の制御フローに示すように、制御装置２７は、上述したクラスタ分析及び判別分析を実行するべく、複数の筒内圧力データＰ（ｉ）を取得し、取得したデータを記憶部（図示せず）に記憶する（＃１０１）。

【0036】

その後、制御装置２７は、筒内圧力データＰ（ｉ）をバンドパスフィルタにより、２０００Ｈｚ未満の低周波数及び１５０００Ｈｚより大きい高周波数をカットしその平均値を０ｋＰａに調整したフィルタ処理後データＰ’（ｉ）を生成する（＃１０２）。

【0037】

更に、制御装置２７は、Ｋ−ｍｅａｎｓクラスタ分析の指標（第１指標の例）として、副室点火直後から９０ｄｅｇまでのフィルタ処理後データＰ’（ｉ）の最大値から最小値を減算したΔＰ値（図３に図示）と、副室点火直後から９０ｄｅｇまでのフィルタ処理後データＰ’（ｉ）のＲＭＳ値を算出する。
ここで、ＲＭＳ値の定義について説明を加えると、ＲＭＳ値は、図４に示すように、フィルタ処理後データＰ’（ｉ）において、副室点火時期（θ＝０ｄｅｇＡＴＤＣ）から９０ｄｅｇ、計４５０ステップの夫々における値を、以下の〔数２〕に示すように、二乗平均平方根を計算したものである。ただし、ｉ_０は副室点火時のクランク角を、Δ９０はクランク角９０ｄｅｇ分のステップ数をそれぞれ示す。
当該ＲＭＳ値及び上述のΔＰ値はいずれも、筒内圧力データＰ（ｉ）（及びフィルタ処理後データＰ’（ｉ））の圧力振動の大きさ（振動強度）を示すものである。

【0038】

【数2】

【0039】

制御装置２７は、上記ΔＰ値とＲＭＳ値とを指標として、１サイクル毎のフィルタ処理後データＰ’（ｉ）を、予め与えられたｋ個（当該実施形態では、２個）のクラスタに分割（分類）する非階層的分類方法であるｋ−Ｍｅａｎｓクラスタ分析工程を実行する（＃１０３）。
当該ｋ−Ｍｅａｎｓクラスタ分析では、ΔＰ値とＲＭＳ値とを軸とする散布図上で、各個体とクラスタ中心との距離から暫定境界を定め、フィルタ処理後データＰ’（ｉ）のクラスタ間移動がなくなるまで収束計算をすることで最適な境界を出力する分析手法であり、分析の結果、ΔＰ値及びＲＭＳ値が小さいフィルタ処理後データＰ’（ｉ）が通常燃焼群に分類され、ΔＰ値及びＲＭＳ値が大きいフィルタ処理後データＰ’（ｉ）が非通常燃焼群（ノッキング燃焼状態のフィルタ処理後データＰ’（ｉ）と、異常燃焼状態のフィルタ処理後データＰ’（ｉ）とから成る群）として分類された。

【0040】

さて、非通常燃焼群としてのノッキング燃焼状態のフィルタ処理後データＰ’（ｉ）と異常燃焼状態のフィルタ処理後データＰ’（ｉ）とは、特に、筒内圧力の振動強度（筒内圧力の経時変化の振幅）が共に大きい。このため、両者を分類するためにはより精密な分析による分類が必要である。本発明にあっては、発明者が定めた教師データ閾値に基づいて抽出する教師データに基づいて、より精密な分析による分類が可能である判別分析を行うことで、非通常燃焼群が分類されるフィルタ処理後データＰ’（ｉ）を、ノッキング燃焼群と異常燃焼群とに分類する。

【0041】

そこで、制御装置２７は、非通常燃焼群のフィルタ処理後データＰ’（ｉ）から、教師データを抽出する処理を実行する（＃１０４）。
当該教師データ抽出工程では、制御装置２７は、図５の制御フローに示すように、フィルタ処理後データＰ’（ｉ）が初めて所定の圧力振動判別閾値（当該実施形態では、１２５ｋＰａ）以上となる圧力振動開始時期のデータを教師データ閾値として、ノッキング燃焼状態としての第１教師データと、異常燃焼状態としての第２教師データとを絞り込む第１絞込工程（＃２０２、２０３）と、フィルタ処理後データＰ’（ｉ）の絶対値の最大値を教師データ閾値として、ノッキング燃焼状態としての第１教師データと、異常燃焼状態としての第２教師データとを絞り込む第２絞込工程（＃２０４、２０５）とを記載の順に実行する。

【0042】

第１絞込工程について説明を追加すると、図６（ＲＭＳ値が高い場合：１６９〜１７０ｋＰａ程度）、図７（ＲＭＳ値が中程度の場合：１０６〜１０７ｋＰａ程度）、図８（ＲＭＳ値が低い場合：７４〜８０ｋＰａ程度）に示すように、どのＲＭＳ値の場合であっても、ノッキング燃焼状態では、圧力振動が筒内圧力のピーク付近（θ＝２０〜３０ｄｅｇ）以降で発生しているのに対し、異常燃焼状態では、圧力振動が副室点火時期（θ＝０ｄｅｇ）の直後から発生していることから、圧力振動開始時期で、両者を切り分けることができる。
このことから、当該実施形態にあっては、制御装置２７は、第１絞込工程において、振動開始時期が７ｄｅｇ以下の場合、異常燃焼状態としての教師データ候補として抽出し（＃２０２）、振動開始時期が１０ｄｅｇより大きい場合、ノッキング燃焼状態としての教師データ候補として抽出する（＃２０３）。

【0043】

次に、第１絞込工程にて絞り込まれた異常燃焼状態としての教師データ候補のフィルタ処理後データＰ’（ｉ）を重ね合わせた図９から判明するように、筒内圧力のピーク付近（ｉ＝１００〜１２０）で強いピークが含まれている。これは、図６中段、図７中段、図８中段の異常燃焼状態のフィルタ処理後データＰ'（ｉ）には存在しないピークであり、ノッキング燃焼状態のフィルタ処理後データＰ’（ｉ）の特徴であると考えられるから、第１絞込工程にて絞り込まれた異常燃焼状態としての教師データ候補には、ノッキング燃焼状態のデータが混在していると言える。
一方で、第１絞込工程にて絞り込まれたノッキング燃焼状態としての教師データ候補のフィルタ処理後データＰ’（ｉ）を重ね合わせた図１０から判明するように、副室点火時期（ｉ＝０）直後から継続する燃焼振動が含まれている。これは、図６下段、図７下段、図８下段のノッキング燃焼状態のフィルタ処理後データＰ’（ｉ）には存在しない燃焼振動であり、異常燃焼状態のフィルタ処理後データＰ’（ｉ）の特徴であると考えられるから、第１絞込工程にて絞り込まれたノッキング燃焼状態としての教師データ候補には、異常燃焼状態のデータが混在していると言える。

【0044】

そこで、当該実施形態の制御装置２７では、より適切な教師データを抽出するべく、第１絞込工程の後の第２絞込工程として、＃２０２で抽出された異常燃焼状態としての教師データ候補のうち、フィルタ処理後データＰ’（ｉ）の絶対値の最大値が７００ｋＰａ未満のものを、異常燃焼状態の教師データとして抽出し（＃２０４）、＃２０３で抽出されたノッキング燃焼状態としての教師データ候補のうち、フィルタ処理後データＰ’（ｉ）の絶対値の最大値が３５０ｋＰａを超えるものを、ノッキング燃焼状態の教師データとして抽出する（＃２０５）。

【0045】

以上の如く、第１絞込工程及び第２絞込工程を記載の順に実行することにより、異常燃焼状態としての教師データと、ノッキング燃焼状態としての教師データとを良好に抽出することができる。
しかしながら、上述のような閾値による判別では、非通常燃焼群としてのデータのうち、異常燃焼状態としての教師データにも、ノッキング燃焼状態としての教師データにも分類されない未分類のデータが多数発生することになる。
そこで、当該実施形態にあっては、多変量解析の一手法として知られる判別分析により、当該第１教師データと第２教師データとを用いて、非通常燃焼群としてのフィルタ処理後データＰ’（ｉ）を未分類のデータを出すことなく、適切に分類する。尚、当該判別分析では、線形判別式を用いて、相関比が最大となる線形判別式を導出し、当該線形判別式により判別する。当該判別分析については、例えば、「多変量解析入門（多変量解析入門、小西 貞則 (著)、 岩波書店）に示されるように公知の分析方法であるので、ここでは、その詳細な説明は割愛する。

【0046】

尚、判別分析の指標（第２指標の例）は、その指標を採用した場合の分布が正規分布とみなされる分布となる指標を好適に採用でき、当該実施形態にあっては、ＲＭＳ／ΔＰ値と、二乗和平方根割合とを用いる。
ここで、二乗和平方根割合とは、図１１下段に示すように、１サイクルにおいて、前記副室の内部の混合気に点火する点火手段の点火時期から判別時期までの圧力振幅積分値に対する、前記点火手段の点火時期から前記主燃焼室６の筒内圧力データＰ（ｉ）が最大となる時期までの圧力振幅積分値の割合の平方根であり、以下の〔数３〕で示される値である。

【0047】

【数3】

【0048】

上述の如く、ＲＭＳ／ΔＰ値と、二乗和平方根割合とを、判別分析の第２指標として採用した理由を以下に示す。
本発明の発明者らは、判別分析時に用いる指標を決定するため、異常燃焼状態のサイクル波形とノッキング燃焼状態のサイクル波形の特徴を把握するべく、複数のＲＭＳ値毎において、両者のサイクル波形を比較した。図６、７、８において、ＲＭＳ値が１６９ｋＰａ付近、１０７ｋＰａ付近、７５ｋＰａ付近となる波形を夫々示す。因みに、図６、７、８で、上段は、筒内圧力データの経時変化を示すグラフ図、中段は、異常燃焼状態のフィルタ処理後データを示すグラフ図、下段は、ノッキング燃焼状態のフィルタ処理後データを示すグラフ図である。

【0049】

図６、７、８の比較から、ＲＭＳ値の大小に関わらず、異常燃焼状態及びノッキング燃焼状態の筒内圧力データＰ（ｉ）及びフィルタ処理後データＰ’（ｉ）を示す波形には、以下の差異がある。
１点目は、筒内圧力の圧力振動の開始時期である。異常燃焼状態では、副室点火時期（θ＝０ｄｅｇ）の直後に圧力振動が開始するのに対して、ノッキング燃焼状態では、筒内圧力のピーク付近で圧力振動が開始している。
２点目は、圧力振動の振幅である。フィルタ処理後のデータの比較から、異常燃焼状態よりもノッキング燃焼状態の方が、圧力最大振幅ΔＰは大きい傾向にあることが確認される。
３点目は、圧力振幅の推移である。フィルタ処理後のデータから、異常燃焼状態では、筒内圧力のピークまでは圧力振動の振幅が略一定で推移し、筒内圧力のピークを過ぎて徐々に振幅が減衰する傾向である。これに対し、ノッキング燃焼状態では、筒内圧力の振動開始直後に振幅が最大となり、その後は単調に減少し、異常燃焼状態のように振幅が一定となる期間は観測されない。

【0050】

これらの比較結果から、異常燃焼状態とノッキング燃焼状態とを判別するには、以下の要件を満たす指標が有効である。
第１要件は、フィルタ処理後の圧力の振幅の推移から、両者を区別できることである。
即ち、異常燃焼状態では振動開始時期が早く、圧力振幅は振動開始から筒内圧力のピークにかけて略一定の値を保ちながら推移した後、筒内圧力のピークを過ぎると減少する。一方、ノッキング燃焼状態では、振動開始時期が筒内圧力のピーク付近で遅く、圧力振幅は振動開始直後に最大に達した後、単調減少する。これらの各燃焼の圧力振幅の推移の特徴を表すことができる指標は、両者を判別できる指標として有効である。
上記第１要件を満たす指標として、フィルタ処理後の圧力が所定の閾値に初めて達した時のクランク角度としての圧力振動開始時期が考えられる。しかしながら、筒内圧力の経時変化は、サイクル変動を伴うため、両燃焼を区別するための振動開始時期を定義する適切な閾値を設定することは難しい。一方、両燃焼の振動開始時期が異なることは、１サイクルの圧力振動成分の積分値に対する筒内圧力がピークに達するまでの圧力振動成分の積分値の割合の大小に対応するため、当該割合を指標として用いることでより適切に、両燃焼を峻別することができる。
そこで、当該指標として、１サイクルにおいて、副室の内部の混合気に点火する点火手段の点火時期から判別時期までの圧力振幅積分値に対する、点火手段の点火時期から主燃焼室６の筒内圧力データＰ（ｉ）が最大となる時期までの圧力振幅積分値の割合である二乗和平方根割合を採用した。

【0051】

第２要件は、フィルタ処理後のＲＭＳ値とΔＰ値の大きさの比の点から、両グループを区別できることである。
上述したように、ノッキング燃焼状態の方が異常燃焼状態よりもΔＰ値が大きくなる傾向がみられた。また、振動開始時期は、異常燃焼状態では早く、ノッキング燃焼状態では遅くなることからも、ＲＭＳ値が一定の条件では、異常燃焼状態のΔＰ値は相対的に小さく、ノッキング燃焼状態のΔＰ値は大きくなることが予想される。このことから、ＲＭＳ値とΔＰ値の大きさの比は両者を判別するための指標として有効である。
このことから、本発明にあっては、フィルタ処理後データＰ’（ｉ）のＲＭＳ値を、フィルタ処理後データＰ’（ｉ）の最大値から最小値を減算したΔＰ値で除算したＲＭＳ／ΔＰ値を、指標として採用した。

【0052】

次に、これまで説明してきた手法により、非通常燃焼群を判別した結果について説明する。
図１２は、判別分析前の非通常燃焼群に属するデータとして、異常燃焼状態としての教師データ、ノッキング燃焼状態としての教師データ、及び教師データ以外のデータ（未分類データ）を、ＲＭＳ／ΔＰ値及び二乗和平方根割合を軸とした散布図上に図示したものであり、図１３は、判別分析後のデータを示すグラフ図である。
因みに、判別分析後の線形判別関数は、以下の〔数４〕で示される関数として導出された。

【0053】

【数4】

【0054】

図１２、１３の比較により、分析前に散布図上に広範囲に分布していた未分類データは、〔数４〕で仕切られた異常燃焼群とノッキング燃焼群との何れかに、適切に判別されていることが判る。

【0055】

〔別実施形態〕
（１）上記実施形態において、クラスタ分析工程は、非階層的分類法の一つで、ｎ個の対象を予め決定されたｋ個のクラスタに分ける分類手法であるｋ−ｍｅａｎｓクラスタ分析を採用する例を示した。
しかしながら、当該クラスタ分析工程では、小さなクラスタを関連付けて大きなクラスタを形成し、その過程を視覚的に把握する階層的分類方法を採用しても構わない。

【0056】

（２）上記実施形態では、副室式エンジン１００は、制御装置２７が、副室式エンジンに対し一体で取り付けられているものを例として説明した。しかしながら当該制御装置２７は、副室式エンジンに対してネットワーク回線等を介して接続されるものであっても構わない。即ち、当該制御装置２７は、監視センター等に設けられ、複数の副室式エンジン１００から1サイクル毎の筒内圧力データＰ（ｉ）を取得可能に構成され、取得した筒内圧力データＰ（ｉ）を一括分析して、当該分析結果に基づいて、夫々の副室式エンジンの制御支援を実行する構成を採用しても構わない。

【0057】

（３）上記実施形態に示す副室式エンジン１００では、副室８の内部に点火プラグ１４の点火点を配置させる構成例を示した。しかしながら、副室８の内部に点火点が設けられない構成も、本発明の権利範囲に含むものである。

【0058】

（４）上記実施形態では、副燃焼室８の混合気に点火する点火プラグ１４を備える構成例を示した。しかしながら、当該点火プラグ１４を設けない副室式エンジンも、本発明の権利範囲に含むものである。

【0059】

（５）上記実施形態において、第１指標として、フィルタ処理後データのＲＭＳ値及びΔＰ値を例示したが、これらの値に限定されない。
第１指標としては、筒内圧力の振動強度を示すパラメータが好適に用いられ、散布図上のユークリッド距離を用いるため、質的変数ではなく量的変数であれば、他の指標を用いても構わない。
また、第２指標としては、主燃焼室６の筒内圧力から導出され、且つ正規分布とみなされる分布となる指標を好適に用いることができる。

【0060】

尚、上記実施形態（別実施形態を含む、以下同じ）で開示される構成は、矛盾が生じない限り、他の実施形態で開示される構成と組み合わせて適用することが可能であり、また、本明細書において開示された実施形態は例示であって、本発明の実施形態はこれに限定されず、本発明の目的を逸脱しない範囲内で適宜改変することが可能である。

【産業上の利用可能性】

【0061】

本発明の副室式エンジンの燃焼状態判別方法、副室式エンジン、及びエンジンシステムは、ノッキング燃焼及び当該ノッキング燃焼とは異なる異常燃焼を峻別して、副室式エンジンにおける燃焼状態を判別できる副室式エンジンの燃焼状態判別方法、副室式エンジン、及びエンジンシステムとして、有効に利用可能である。

【符号の説明】

【0062】

４ ：燃焼室
５ ：ピストン
６ ：主燃焼室
７ ：連通孔
８ ：副室
１４ ：点火プラグ
２７ ：制御装置
１００ ：副室式エンジン
Ｓ６ ：筒内圧センサ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】