シリンダ毎のノック強度と相関のある物理量を検出するノックセンサと、前記ノックセンサの検出信号に対して周波数解析を行い、前記物理量に関するシリンダ毎のパワースペクトルを算出するパワースペクトル算出手段と、ノッキング周波数を含む所定の周波数帯域における前記パワースペクトルの部分オーバーオール値を算出するPOA算出手段と、少なくとも前記部分オーバーオール値に基づいて内燃機関を制御するコントローラとを備えることを特徴とする内燃機関の制御装置。前記コントローラは、いずれかのシリンダに関する前記部分オーバーオール値が第1閾値を超えた場合に内燃機関を停止することを特徴とする請求項1に記載の内燃機関の制御装置。前記コントローラは、いずれかのシリンダに関する前記部分オーバーオール値が前記第1閾値よりも小さい第2閾値を超えた場合に、そのシリンダの負荷を下げることを特徴とする請求項2に記載の内燃機関の制御装置。前記部分オーバーオール値が所定の閾値を超える頻度であるシビアリティをシリンダ毎に算出するシビアリティ算出手段をさらに備え、前記コントローラは、前記部分オーバーオール値および前記シビアリティに基づいて内燃機関を制御することを特徴とする請求項1乃至3のいずれか一項に記載の内燃機関の制御装置。前記コントローラは、いずれかのシリンダに関する前記シビアリティが第3閾値を超えた場合に、そのシリンダの負荷を下げることを特徴とする請求項4に記載の内燃機関の制御装置。前記コントローラは、いずれかのシリンダに関する前記シビアリティが前記第3閾値よりも小さい第4閾値を超えた場合に、そのシリンダの着火タイミングを遅らせることを特徴とする請求項5に記載の内燃機関の制御装置。前記部分オーバーオール値が所定の閾値を超える頻度であるシビアリティをシリンダ毎に算出するシビアリティ算出手段をさらに備え、前記コントローラは、前記シビアリティを前記部分オーバーオール値で重み付けたノック評価値に基づいて内燃機関を制御することを特徴とする請求項1に記載の内燃機関の制御装置。ノック強度と相関のある前記物理量はシリンダの筒内圧であり、前記ノックセンサは、各シリンダの筒内圧を検出する筒内圧センサであることを特徴とする請求項1乃至7のいずれか一項に記載の内燃機関の制御装置。ノック強度と相関のある前記物理量はシリンダの振動レベルであり、前記ノックセンサは、各シリンダの振動レベルを検出する加速度センサであることを特徴とする請求項1乃至7のいずれか一項に記載の内燃機関の制御装置。前記内燃機関がガスエンジンである請求項1乃至9のいずれか一項に記載の内燃機関の制御装置。シリンダ毎のノック強度と相関のある物理量を検出するノックセンサと、前記ノックセンサの検出信号に対して周波数解析を行い、前記物理量に関するシリンダ毎のパワースペクトルを算出するパワースペクトル算出手段と、ノッキング周波数を含む所定の周波数帯域における前記パワースペクトルの部分オーバーオール値を算出するPOA算出手段と、少なくとも前記部分オーバーオール値に基づいて異常燃焼の有無をシリンダ毎に判定する異常燃焼判定手段とを備えることを特徴とする内燃機関の燃焼診断装置。しかしながら、シビアリティは、緩やかな燃焼状態の変化を捉えるのに適しているが、突発的に生じる燃焼状態の大きな変化を捉えることができない。特に、ガスエンジンは、燃料の性状の影響等を受けて突発的に大きなノッキングを起こす場合があるから、シビアリティを用いた燃焼診断では十分と言えない。したがって、特許文献1のように評価ノッキング頻度シビアリティを用いる場合、ガスエンジンの燃焼状態に応じた細やかな制御を行うことができない。本発明は、上述の事情に鑑みてなされたものであり、突発的に生じる燃焼状態の大きな変化を把握しうる内燃機関の燃焼診断装置と、燃焼状態に応じた細やかな制御を行いうる内燃機関の制御装置を提供することを目的とする。本発明に係る内燃機関の制御装置は、シリンダ毎のノック強度と相関のある物理量を検出するノックセンサと、前記ノックセンサの検出信号に対して周波数解析を行い、前記物理量に関するシリンダ毎のパワースペクトルを算出するパワースペクトル算出手段と、ノッキング周波数を含む所定の周波数帯域における前記パワースペクトルの部分オーバーオール値を算出するPOA算出手段と、少なくとも前記部分オーバーオール値に基づいて内燃機関を制御するコントローラとを備えることを特徴とする。なお、内燃機関で起こるノッキングの周波数は、シリンダのボア径と音速によって推定できることが知られており、本明細書ではこれをノッキング周波数という。具体的には、ノッキング周波数=ボア径/音速の関係が成立する。この内燃機関の制御装置では、ノックセンサで検出されたノック強度と相関のある物理量からパワースペクトル算出手段によってパワースペクトルを求め、このパワースペクトルの部分オーバーオール値をPOA算出手段によって算出する。このようにして得られる部分オーバーオール値には、ノック強度の大きさが反映される。したがって、部分オーバーオール値を用いることによって、突発的に生じる燃焼状態の大きな変化に応じて細やかな制御が可能となり、安定した運転を行うことができる。また、狭い運転領域での安定運用が可能であるから、高出力高効率な運転を維持できる。上記内燃機関の制御装置において、前記コントローラは、いずれかのシリンダに関する前記部分オーバーオール値が第1閾値を超えた場合に内燃機関を停止してもよい。このように、部分オーバーオール値が第1閾値を超えたシリンダが一つでもあれば内燃機関を停止することで、突発的に生じる大きなノッキングに応じて即座に内燃機関を停止して、安全を確保できる。この場合、前記コントローラは、いずれかのシリンダに関する前記部分オーバーオール値が前記第1閾値よりも小さい第2閾値を超えた場合に、そのシリンダの負荷を下げてもよい。このように、部分オーバーオール値が第1閾値よりも小さい第2閾値を超えたシリンダすなわち、内燃機関全体を停止するまでもないが、ノッキングにつながりうる程度の比較的大きな振動が突発的に発生したシリンダについて負荷を下げることで、該シリンダのノッキングの発生を抑制できる。これにより、より一層の安定運転を実現できる。上記内燃機関の制御装置は、前記部分オーバーオール値が所定の閾値を超える頻度であるシビアリティをシリンダ毎に算出するシビアリティ算出手段をさらに備え、前記コントローラは、前記部分オーバーオール値および前記シビアリティに基づいて内燃機関を制御してもよい。このように部分オーバーオール値およびシビアリティに基づいて内燃機関を制御することで、突発的に生じる燃焼状態の大きな変化および緩やかな燃焼状態の変化に応じてより細やかな制御が可能となる。この場合、前記コントローラは、いずれかのシリンダに関する前記シビアリティが第3閾値を超えた場合に、そのシリンダの負荷を下げてもよい。このように、シビアリティが第3閾値を超えたシリンダすなわち、所定の強度を超える振動が頻発しているシリンダについて負荷を下げることで、該シリンダのノッキングの発生を抑制できる。これにより、より一層の安定運転を実現できる。上記内燃機関の制御装置において、前記コントローラは、いずれかのシリンダに関する前記シビアリティが前記第3閾値よりも小さい第4閾値を超えた場合に、そのシリンダの着火タイミングを遅らせてもよい。このように、シビアリティが第3閾値よりも小さい第4閾値を超えたシリンダすなわち、所定の強度を超える振動がやや高い頻度で発生しているシリンダについて着火タイミングを遅らせることで、負荷を下げることなく、該シリンダのノッキングの発生を抑制しつつ様子を見ることができる。上記内燃機関の制御装置は、前記部分オーバーオール値が所定の閾値を超える頻度であるシビアリティをシリンダ毎に算出するシビアリティ算出手段をさらに備え、前記コントローラは、前記シビアリティを前記部分オーバーオール値で重み付けたノック評価値に基づいて内燃機関を制御してもよい。このようにシビアリティを部分オーバーオール値で重み付けたノック評価値に基づいて内燃機関を制御することで、ノック評価値という一つのパラメータのみによって、突発的に生じる燃焼状態の大きな変化と緩やかな燃焼状態の変化との両方を把握して、燃焼状態に応じた細やかな制御を行うことができる。上記内燃機関の制御装置において、ノック強度と相関のある前記物理量はシリンダの筒内圧であり、前記ノックセンサは、各シリンダの筒内圧を検出する筒内圧センサであってもよい。筒内圧センサは、各シリンダの筒内圧を直接検出してシリンダ内の燃焼状態を正確に把握して、より一層細やかな制御を行いたい場合に好適である。あるいは、ノック強度と相関のある前記物理量はシリンダの振動レベルであり、前記ノックセンサは、各シリンダの振動レベルを検出する加速度センサであってもよい。加速度センサは、筒内圧センサに比べて安価であるから、上記内燃機関の制御装置を低コストで実現したい場合に好適である。上記内燃機関の制御装置はガスエンジンに適用してもよい。ガスエンジンは、燃料の性状の影響等を受けて突発的に大きなノッキングを起こす場合がある。この点、上記内燃機関の制御装置では、ノック強度の大きさが反映された部分オーバーオール値に基づいて内燃機関の制御を行うから、突発的に生じる燃焼状態の大きな変化に適切に対処して、高効率であり且つ安定した運転を実現できる。本発明に係る内燃機関の燃焼診断装置は、シリンダ毎のノック強度と相関のある物理量を検出するノックセンサと、前記ノックセンサの検出信号に対して周波数解析を行い、前記物理量に関するシリンダ毎のパワースペクトルを算出するパワースペクトル算出手段と、ノッキング周波数を含む所定の周波数帯域における前記パワースペクトルの部分オーバーオール値を算出するPOA算出手段と、少なくとも前記部分オーバーオール値に基づいて異常燃焼の有無をシリンダ毎に判定する異常燃焼判定手段とを備えることを特徴とする。この内燃機関の燃焼診断装置では、ノックセンサで検出されたノック強度と相関のある物理量からパワースペクトル算出手段によってパワースペクトルを求め、このパワースペクトルの部分オーバーオール値をPOA算出手段によって算出する。このようにして得られる部分オーバーオール値には、ノック強度の大きさが反映される。したがって、部分オーバーオール値に基づく異常燃焼の判定によって、突発的に生じる燃焼状態の大きな変化を把握することができる。本発明では、ノックセンサで検出されたノック強度と相関のある物理量からパワースペクトル算出手段によってパワースペクトルを求め、このパワースペクトルの部分オーバーオール値をPOA算出手段によって算出する。このようにして得られる部分オーバーオール値には、ノック強度の大きさが反映される。したがって、部分オーバーオール値を用いることによって、突発的に生じる燃焼状態の大きな変化に応じて細やかな制御が可能となり、高効率であり且つ安定した運転を行うことができる。また、部分オーバーオール値に基づく異常燃焼の判定によって、突発的に生じる燃焼状態の大きな変化を把握することができる。以下、添付図面に従って本発明の実施形態について説明する。ただし、この実施形態に記載されている構成部品の寸法、材質、形状、その相対的配置等は、特定的な記載がない限り本発明の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。[第1実施形態]図1は、第1実施形態に係るガスエンジンの構成例を示す図である。図2は、筒内圧センサで計測した筒内圧とクランク角度との関係の一例を示すグラフである。図3は、筒内圧センサの検出結果から得られるパワースペクトルの一例を示すグラフである。以下、内燃機関の一例として副室式ガスエンジンを採り上げるが、本発明が他の形式のガスエンジンやディーゼルエンジンを含むあらゆる内燃機関に適用しうることは言うまでもない。なお、副室式ガスエンジンは、シリンダ内に設けられた主室とは別に設けた副室内で混合ガス燃料を着火し、副室から主室に向けて火炎を噴出することで、主室内の希薄混合ガス燃料を燃焼させる方式のガスエンジンをいう。図1に示すように、ガスエンジン1は、シリンダ2と、シリンダ2内に収納され往復運動するピストン3と、給気ポート4及び排気ポート5と、給気ポート4を開閉する給気弁と、排気ポート5を開閉する排気弁と、シリンダ2内の燃焼状態を判定する燃焼診断装置30と、ガスエンジン1の各部を制御するコントローラ40とを有する。シリンダ2とピストン3上面とで囲まれた空間は主室6である。また、シリンダ2上方には、副室8が設けられており、主室6と副室8とは、複数個例えば4~10個の副室噴孔によって連通している。主室6の上部には、主室6内の圧力筒内圧を計測する筒内圧センサノックセンサに相当26が設けられている。筒内圧センサ26により検出された筒内圧Pは、後述する燃焼診断装置30のパワースペクトル算出部32に送られる。またピストン3により回転駆動されるクランク軸10には、クランク軸10のクランク角度θを検出するクランク軸角度検出器24が設けられている。クランク軸角度検出器24により検出されたクランク角度θは、後述する燃焼診断装置30のパワースペクトル算出部32に送られる。副室8の上部には、副室8内の混合気を火花点火するための着火装置としての点火プラグ9が設けられている。また副室8内には、副室ガス供給通路点火用補助ガス通路11を介して燃料ガスが供給されるようになっている。給気ポート4には、主室ガス供給弁22が介装された燃料ガス供給管が接続されており、不図示の過給機から供給された空気と主室ガス供給弁22を通過した燃料ガスとが給気ポート4内で混合され、希薄混合気となった気体が主室6内に給気されるようになっている。なお、図1に示すように、吸気ポート4に接続された燃料ガス供給管に主室ガス圧力調整弁23を設け、吸気ポート4内に供給される燃料ガスの流量及び/又は圧力を調整してもよい。一方、副室ガス供給通路11には、副室ガス供給弁20が介装された燃料ガス供給管が接続されており、副室ガス供給弁20を通過した燃料ガスが副室ガス供給通路11を介して副室8内に供給されるようになっている。なお、図1に示すように、副室ガス供給通路11に接続された燃料ガス供給管に副室ガス圧力調整弁21を設け、副室8に供給される燃料ガスの流量及び/又は圧力を調節してもよい。燃焼診断装置30は、パワースペクトル算出部32、POA算出部34および異常燃焼判定部36により構成される。燃焼診断装置30の各部の働きは、次のとおりである。パワースペクトル算出部32は、筒内圧センサ26で計測された筒内圧Pと、クランク軸角度検出器24で検出されたクランク角度θを受け取って、筒内圧Pとクランク角度θとの関係から周波数解析により筒内圧パワースペクトルを取得する。具体的には、図2に示す筒内圧Pとクランク角度θとの関係を示す筒内圧波形100を求め、この筒内圧波形100を任意のタイムウインドウΔθで周波数解析FFT分析し、筒内圧パワースペクトル102図3参照を算出する。タイムウインドウΔθは、上死点付近の任意の範囲に設定可能であり、例えば、上死点から40度のクランク角度θの範囲であってもよい。POA算出部34は、パワースペクトル算出部32により得られた筒内圧パワースペクトル102からノック強度の大きさを表す部分オーバーオール値以下、筒内圧POAと称する。を求める。すなわち、POA算出部34は、ノッキング周波数を含む任意の周波数帯域Δνにおける筒内圧パワースペクトル102の筒内圧POAを算出する。このときに用いる周波数帯域Δνは、例えば、ノッキング周波数±500Hzであってもよい。ここで、ノッキング周波数はシリンダのボア径と音速によって推定でき、ノッキング周波数=ボア径/音速で求まる。また、筒内圧POAは、周波数帯域Δνにおける筒内圧パワースペクトル102の二乗平均値であり、ノック強度を反映したものである。異常燃焼判定部36は、POA算出部34によって求めた筒内圧POAに基づいてシリンダ2の異常燃焼の有無を判定する。例えば、筒内圧POAが第1閾値>第2閾値を超えている場合には、ガスエンジン1の運転停止が必要な程度の大きな強度の振動が発生していると判断する。筒内圧POAが第2閾値を超えている場合には、ガスエンジン1の運転を停止する必要はないものの、何らかの対策が必要な程度の振動が発生していると判断する。なお、異常燃焼判定部36による異常燃焼の有無はサイクル毎に判定する。また、シリンダ2が複数ある場合には、各シリンダ2毎に異常燃焼の有無を判定する。この異常燃焼判定部36の判定結果に基づいて、コントローラ40はガスエンジン1の各部を制御し、必要に応じて、ガスエンジン1の運転を停止したり、シリンダ2の負荷を下げたり、点火プラグ9の点火タイミングを遅らせたりする。上記構成のガスエンジン1では、筒内圧センサ26で検出された筒内圧Pからパワースペクトル算出部32によって筒内圧パワースペクトル102を求め、この筒内圧パワースペクトル102の筒内圧POAをPOA算出部34によって算出する。このようにして得られる筒内圧POAには、ノック強度の大きさが反映される。したがって、筒内圧POAに基づく異常燃焼の判定によって、突発的に生じる燃焼状態の大きな変化を把握することができる。また、筒内圧POAを用いることによって、突発的に生じる燃焼状態の大きな変化に応じて細やかな制御が可能となり、安定した運転を行うことができる。また、狭い運転領域での安定運用が可能であるから、高出力高効率な運転を維持できる。なお、本実施形態では、シリンダ2の筒内圧Pを検出する筒内圧センサ26を用いる例を示したが、シリンダ2毎のノック強度と相関のある物理量を検出する任意のノックセンサを用いてもよい。例えば、筒内圧センサ26に替えてシリンダ2の振動レベルを検出する加速度センサを用いて、次のように、加速度センサの検出信号から求めた振動レベルスペクトルの部分オーバーオール値以下、振動レベルPOAと称する。に基づいて異常燃焼を判定してもよい。図4は、加速度センサで検出された振動レベルとクランク角度θとの関係の一例を示すグラフである。図5は、加速度センサの検出結果から得られるパワースペクトルの一例を示すグラフである。パワースペクトル算出部32において、図4に示す振動レベル波形104を任意のタイムウインドウΔθで周波数解析FFT分析し、振動レベルパワースペクトル106図5参照を算出する。このとき、タイムウインドウΔθを適切な範囲に設定することで、計測対象であるシリンダ2の振動レベルだけを抽出して、他のシリンダ2の振動の影響を排除した振動レベルパワースペクトル106を各シリンダ2について取得できる。この振動レベルパワースペクトル106の周波数帯域Δνにおける振動レベルPOAをPOA算出部34によって求める。このようにして得られた振動レベルPOAは、筒内圧パワースペクトル102から求めた部分オーバーオール値筒内圧POAと同様にシリンダ2のノック強度を反映したものであり、異常燃焼判定部36における燃焼状態の診断に用いることができる。なお、加速度センサは、シリンダ2のノック強度と相関のある振動レベルを検出可能であればシリンダ2又はシリンダヘッドの任意の位置に設けてもよい。例えば、シリンダ2にシリンダヘッドを取り付けるためのヘッドボルト上、シリンダヘッドの側面又は上面に加速度センサを設けてもよい。なかでも、ヘッドボルト上は、シリンダヘッドの内部構造冷却水通路等による影響を受けずに、シリンダ2の燃焼状態と相関の強い振動レベルを検出可能であることが経験的に分かっており、加速度センサを設ける位置として好適である。[第2実施形態]図6は、第2実施形態に係るガスエンジンの構成例を示す図である。なお、本実施形態に係るガスエンジンは、POA算出部で求めたサイクル毎の部分オーバーオール値POAと、シビアリティとを併用して燃焼状態を診断するようにした点を除けば、上述のガスエンジン1と同様である。したがって、ここではガスエンジン1と異なる点を中心に説明し、ガスエンジン1と共通する部分には図1と同一の符号を付してその説明を省略する。図6に示すガスエンジン50の燃焼診断装置30は、シビアリティ算出部38を有する。シビアリティ算出部38は、POA算出部34で算出した部分オーバーオール値POA筒内圧POA又は振動レベルPOAが所定の閾値を超えた頻度であるシビアリティSevを算出する。シビアリティSevは、部分オーバーオール値が所定の閾値を超える頻度を示す指標であり、具体的には、過去Nサイクル中、部分オーバーオール値が所定の閾値を超えたサイクル数をn回とすると以下の式1で表される。数1シビアリティ=n/N×1001異常燃焼判定部36は、POA算出部34により取得した現サイクルの部分オーバーオール値POAと、シビアリティ算出部38により取得した過去NサイクルのシビアリティSevとに基づいてシリンダ2の燃焼状態を診断する。そして、異常燃焼判定部36の判定結果に基づいて、コントローラ40がガスエンジン50の各部を制御する。ここで、ガスエンジン50における燃焼診断装置30による燃焼診断およびコントローラ40による各部制御の手順について説明する。図7は、ガスエンジン50における燃焼診断および各部制御の手順を示すフローチャートである。まず、ステップS2において、POA算出部34によって部分オーバーオール値POA筒内圧POA又は振動レベルPOAを算出し、シビアリティ算出部38によって過去NサイクルのシビアリティSevを算出する。図8は、部分オーバーオール値POAとシビアリティSevの経時変化の一例を示したグラフであり、図8aがサイクル毎の部分オーバーオール値を示し、図8bが過去NサイクルのシビアリティSevを示している。ステップS4では、現サイクルの部分オーバーオール値POAが第1閾値図8aの閾値A。以下、閾値Aと称する。よりも大きいか否かを判定する。この判定は各シリンダ2について行う。そして、いずれかのシリンダ2に関する部分オーバーオール値が閾値Aを超えている場合ステップS4のYES判定、ステップS5に進んでガスエンジン1の運転を停止する。一方、いずれのシリンダ2についても部分オーバーオール値POAが閾値A以下である場合ステップS4のNO判定、ステップS6に進んで、現サイクルの部分オーバーオール値POAが第2閾値図8aの閾値B<閾値A。以下、閾値Bと称する。よりも大きいか否かを判定する。この判定は各シリンダ2について行う。なお、ステップS4及び6で用いる閾値A及びBは、シビアリティSevを算出する際に用いる閾値所定の閾値に相当。よりも大きく設定されている。そして、いずれかのシリンダ2に関する部分オーバーオール値が閾値Bを超えている場合ステップS6のYES判定、ステップS7に進んでそのシリンダ2の負荷を下げた後、ステップS2に戻る。具体的には、コントローラ40の制御下において、主室ガス圧力調整弁23が吸気ポート4内に供給される燃料ガスの流量及び/又は圧力を低減し、シリンダ2の負荷を下げる。なお、シリンダ2の負荷の下げ代例えば10~20や下げ速度は任意に設定可能である。一方、いずれのシリンダ2についても部分オーバーオール値が閾値B以下である場合ステップS6のNO判定、ステップS8に進んで過去NサイクルのシビアリティSevが第3閾値図8bの閾値C。以下、閾値Cと称する。よりも大きいか否かを判定する。この判定は各シリンダ2について行う。ステップS8において、いずれかのシリンダ2に関するシビアリティSevが閾値Cよりも大きいと判定された場合YES判定、ステップS9に進んでそのシリンダ2の負荷を下げた後、ステップS2に戻る。具体的には、コントローラ40の制御下において、主室ガス圧力調整弁23が吸気ポート4内に供給される燃料ガスの流量及び/又は圧力を低減し、シリンダ2の負荷を下げる。なお、シリンダ2の負荷の下げ代例えば10~20や下げ速度は任意に設定可能である。一方、いずれのシリンダ2についてもシビアリティSevが閾値C以下である場合ステップS8のNO判定、ステップ10に進んで過去NサイクルのシビアリティSevが第4閾値図8bの閾値D<閾値C。以下、閾値Dと称する。よりも大きいか否かを判定する。この判定は各シリンダ2について行う。ステップS10において、いずれかのシリンダ2に関するシビアリティSevが閾値Dよりも大きいと判定された場合YES判定、ステップS11に進んでそのシリンダ2についてタイミングリタード着火タイミングを遅らせることを行った後、ステップS2に戻る。具体的には、コントローラ40の制御下において、点火プラグ9の点火タイミングを遅らせて、タイミングリタードを行う。なお、着火タイミングのリタード量例えばクランク角度換算で2度は任意に設定可能である。一方、いずれのシリンダ2についてもシビアリティSevが閾値D以下である場合ステップS10のNO判定、ステップS12に進んで、各シリンダ2の負荷及び着火タイミングが初期設定と異なるか否かを判断する。いずれかのシリンダ2の負荷及び着火タイミングの少なくとも一方が初期設定と異なる場合ステップS12のYES判定、ステップS13に進んで、そのシリンダ2の負荷又は着火タイミングを変更してから所定期間Xを経過したか否かを判断する。そして、ステップS13において所定期間Xを経過していると判断された場合YES判定、ステップS14に進んで、そのシリンダ2の負荷又は着火タイミングを初期設定に戻した後、ステップS2に戻る。このとき、負荷又は着火タイミングの初期設定への戻し速度は任意に決定可能である。一方、ステップS12において、いずれのシリンダ2の負荷及び着火タイミングも初期設定と同一である場合NO判定、ステップS2に戻る。同様に、ステップS13において所定期間Xを経過していないと判定された場合NO判定、ステップS2に戻る。このような手順で燃焼診断およびガスエンジン1各部の制御を行うことで、突発的に非常に大きい強度のノッキングが発生した時点ですなわちPOA>閾値Aの場合、ガスエンジン1の運転を即座に停止して安全を確保できる。また、いずれかのシリンダ2について、比較的大きな強度の振動が突発的に発生した時点すなわちPOA>閾値Bの場合、あるいは、所定の強度を超える振動が頻発したすなわちSev>閾値Cの場合、そのシリンダ2の負荷を下げて安定した運転を行うことができる。さらに、比較的大きな強度の振動がやや高い頻度で発生した場合すなわちSev>閾値Dの場合、そのシリンダ2の着火タイミングを遅らせて、負荷を下げることなく、該シリンダ2のノッキングの発生を抑制しつつ様子を見ることができる。なお、本実施形態では、部分オーバーオール値POAおよびシビアリティSevという二つのパラメータを指標として燃焼状態を判定する例について説明したが、シビアリティを部分オーバーオール値で重み付けたノック評価値に基づいて異常燃焼の有無を判定するようにしてもよい。例えば、以下の式2で表されるノック評価値を用いてもよい。数2ただし、i=1、2、...、Nであり、POAiはiサイクル目の部分オーバーオール値であり、XiはPOAiが所定の閾値を超えた場合は1、そうでない場合はゼロである。シビアリティを部分オーバーオール値で重み付けたノック評価値に基づいて内燃機関を制御することで、ノック評価値という一つのパラメータのみによって、突発的に生じる燃焼状態の大きな変化と緩やかな燃焼状態の変化との両方を把握して、燃焼状態に応じた細やかな制御を行うことができる。以上、本発明の実施形態について詳細に説明したが、本発明はこれに限定されず、本発明の要旨を逸脱しない範囲において、各種の改良や変形を行ってもよいのはいうまでもない。