複数の気筒を備えたガスエンジンの燃焼制御装置において、前記複数の気筒の中から選択された代表気筒の筒内圧を検出する筒内圧センサーと、前記筒内圧センサーで検出された代表気筒の筒内圧波形から、代表気筒のノッキング限界点火時期を算出する第1の算出手段と、筒内圧波形以外の運転状態量、エンジン基本諸元及び補正係数とを含む算式から、代表気筒のノッキング限界点火時期を算出する第2の算出手段と、前記第1の算出手段で算出された代表気筒のノッキング限界点火時期と、前記第2の算出手段で算出された代表気筒のノッキング限界点火時期とを比較し、前記補正係数を求める第3の算出手段と、求めた補正係数と、前記運転状態量及びエンジン基本諸元とを前記算式に代入し、代表気筒以外の気筒のノッキング限界点火時期を算出する第4の算出手段と、前記第1の算出手段及び第4の算出手段で算出された全気筒のノッキング限界点火時期に基づいて各気筒のノッキング限界点火時期を設定する点火時期設定手段とを備えていることを特徴とするガスエンジンの燃焼制御装置。前記運転状態量が各気筒に供給される給気の温度、空気過剰率及びエンジン回転数であり、前記エンジン基本諸元がシリンダボア、ストローク、圧縮比及びコンロッド長であることを特徴とする請求項1に記載のガスエンジンの燃焼制御装置。前記第4の算出手段において、予め取得したガスエンジンの運転試験データに基づいて、前記補正係数を気筒毎に修正するようにしたことを特徴とする請求項1又は2に記載のガスエンジンの燃焼制御装置。各気筒のうち最もノッキングが起きやすい気筒を代表気筒に選択することを特徴とする請求項1に記載のガスエンジンの燃焼制御装置。しかし、筒内圧センサーは高価であり、気筒数が多数ある場合、筒内圧センサーを全気筒に付設すると、高コストになる。また、筒内圧センサーの故障や誤動作の確率も高くなるという問題がある。本発明は、かかる従来技術の課題に鑑み、複数の気筒を備えたガスエンジンにおいて、簡易かつ低コストな手段で、全気筒のノッキングを抑制しつつ、点火時期を気筒毎に最大限進角することで、ガスエンジンの出力及び熱効率を向上させることを目的とする。かかる目的を達成するため、本発明のガスエンジンの燃焼制御装置は、複数の気筒を備えたガスエンジンの燃焼制御装置において、複数の気筒の中から選択された代表気筒の筒内圧を検出する筒内圧センサーと、筒内圧センサーで検出された代表気筒の筒内圧波形から、代表気筒のノッキング限界点火時期を算出する第1の算出手段と、筒内圧波形以外の運転状態量、エンジン基本諸元及び補正係数とを含む算式から、代表気筒のノッキング限界点火時期を算出する第2の算出手段と、第1の算出手段で算出された代表気筒のノッキング限界点火時期と、第2の算出手段で算出された代表気筒のノッキング限界点火時期とを比較し、補正係数を求める第3の算出手段と、求めた補正係数と、運転状態量及びエンジン基本諸元とを前記算式に代入し、代表気筒以外の気筒のノッキング限界点火時期を算出する第4の算出手段と、第1の算出手段及び第4の算出手段で算出された全気筒のノッキング限界点火時期に基づいて各気筒のノッキング限界点火時期を設定する点火時期設定手段とを備えていることを特徴とする。本発明では、全気筒に筒内圧センサーを設ける代わりに、代表気筒を選択し、代表気筒のみに筒内圧センサーを設ける。他の気筒では、筒内圧センサーを設けることなく、他の運転状態量及びエンジン基本諸元からノッキング限界点火時期を算出するので、低コスト化できる。そして、筒内圧センサーで検出した筒内圧波形と、各気筒で検出した運転状態量とから、気筒毎にノッキング限界点火時期を算出できる。そのため、低コストで出力、燃費及び熱効率が良好な運転を実現できる。前記運転状態量として、各気筒に供給される給気の温度、空気過剰率及びエンジン回転数を選択するとよい。これによって、ガスエンジンに装着するセンサー類を削減でき、低コスト化できる。また、前記エンジン基本諸元として、シリンダボア、ストローク、圧縮比及びコンロッド長を選択するとよい。基本諸元としてこれらのファクタを選択することで、各気筒のノッキング限界点火時期を正確に算出できる。第4の算出手段において、予め取得したガスエンジンの運転試験データに基づいて、補正係数を気筒毎に修正するようにするとよい。これによって、気筒毎にさらに正確なノッキング限界点火時期を算出できる。本発明では、代表気筒のノッキング限界点火時期を元に、他の気筒のノッキング限界点火時期を算出している。そのため、最もノッキングが起きやすい気筒を代表気筒に選択し、その気筒のノッキング限界点火時期に基づいて、他の気筒のノッキング限界点火時期を算出することで、ノッキングを確実に抑制できる安全サイドの運転状態から点火時期を進角し、ノッキング限界点火時期に調整することが可能である。本発明によれば、全気筒に筒内圧センサーを設ける代わりに、筒内圧センサーを代表気筒に1個のみ設けただけの低コストな手段で、代表気筒の筒内圧波形と各気筒の運転状態量及びエンジン基本諸元とから、全気筒のノッキング限界点火時期を気筒毎に求めることができる。そのため、全気筒のノッキングを抑制しつつ、点火時期を最大限進角でき、ガスエンジンの出力及び熱効率を向上できる。以下、本発明を図に示した実施形態を用いて詳細に説明する。但し、この実施形態に記載されている構成部品の寸法、材質、形状、その相対配置などは特に特定的な記載がない限り、この発明の範囲をそれのみに限定する趣旨ではない。本発明をガスエンジンの燃焼制御に適用した一実施形態を図1~図4に基づいて説明する。図1において、本実施形態に係るガスエンジン10は、エンジンブロック12に一列に6個の気筒14a~fが設けられている。各気筒14a~fのうち左端に配置された気筒14aを代表気筒に選択している。気筒14aには筒内圧センサー16が設けられている。各気筒14a~fには、エアクーラ18から給気マニホールド20を介して給気aが送られる。各気筒14a~fで燃焼した後の排気eは、夫々排気マニホールド22を介して排気管24に集合し、排気管24から排出される。エンジンブロック12に設けられたクランク軸図示省略に、上死点及び下死点を含むクランク角を検出する回転センサー26が設けられている。筒内圧センサー16及び回転センサー26から、クランク角に対応した筒内圧波形Pが検出される。サイクル毎の筒内圧の変化より、ノッキング発生頻度とノッキング限界点火時期が導出される。各気筒14a~fの入口に給気温度センサー28a~fが設けられている。図1に示すように、各気筒14a~fに吸入される給気温度Tsは夫々異なっている。気筒列の両端にある気筒14a及び14fの給気温度Tsが一番高く、中央の気筒14c及び14dの給気温度Tsが一番低い。エアクーラ18や給気マニホールド20は両端側で、給気が滞留する傾向があるため、壁温が高くなる傾向がある。そのため、エンジンブロック12の両端に配置された気筒14a及び14fに流入する給気の温度Tsが最も高くなる傾向にある。従って、エンジンブロック12の両端に配置された気筒でノッキングが起きやすくなる。本実施形態では、左端に配置された気筒14aを代表気筒に選択している。図2にガスエンジン10の燃焼制御装置30を示す。図2において、第1算出部32には、筒内圧センサー16で検出された代表気筒14aの筒内圧波形が入力される。第1算出部32では、この筒内圧波形から代表気筒14aのノッキング限界点火時期を算出する。第2算出部34では、運転状態量、エンジン基本諸元及び補正係数Cを含む次の算式1から、代表気筒14aのノッキング限界点火時期を算出する。Ti=fC、B、S、ε、CR、Ne、λ、Ts1ここで、Cは補正係数である。Bはシリンダボア、Sはストローク、εは圧縮比、CRはコンロッド長であり、これらはエンジン基本諸元である。Neはエンジン回転数、λは空気過剰率、Tsは給気温度である。これらは運転状態量である。エンジン回転数Neは回転センサー26で検出され、給気温度Tsは、給気温度センサー28a~fで気筒毎に検出される。空気過剰率λは、配置通路に酸素濃度センサーを設け、酸素濃度センサーにより設定値どおりの空気過剰率λになっているかどうか確認するのが望ましい。コンプレッサ前でガス燃料の供給を行う場合は、気筒間ばらつきが生じないので、酸素濃度センサーは排気管24に1個設ければよい。第3算出部36では、第1算出部32で算出された代表気筒14aのノッキング限界点火時期と、第2算出部34で算出された代表気筒14aのノッキング限界点火時期とを比較し、補正係数Cを求める。第4算出部38では、第3算出部36で求められた補正係数Cに基づいて、各気筒毎に運転状態量及びエンジン基本諸元を算式1に入力することで、代表気筒14a以外の各気筒のノッキング限界点火時期を算出する。点火時期設定部40では、算出された各気筒のノッキング限界点火時期に基づいて、点火装置図示省略における各気筒のノッキング限界点火時期を設定する。図3は、制御装置30によるノッキング限界点火時期設定手順を示すフロー図である。図3において、まず、筒内圧センサー16で代表気筒14aの筒内圧波形を検出するS10。第1算出部32で、検出された筒内圧波形から代表気筒14aのノッキング限界点火時期Tiaを算出するS12。一方、前述のように、ガスエンジン10に付設された各センサーで運転状態量として、エンジン回転数Ne、空気過剰率λ及び気筒毎の給気温度Tsを選択するS14。これらの運転状態量及び前記エンジン基本諸元が第2算出部34に入力される。第2算出部34では、算式1から代表気筒14aのノッキング限界点火時期Tia’が算出されるS16。次に、第3算出部36で、第1算出部32で算出したノッキング限界点火時期Tiaと、第2算出部34で算出したノッキング限界点火時期Tia’とを比較し、補正係数Cを算出するS18。即ち、算式1にノッキング限界点火時期Tiaと、前記運転状態量及びエンジン基本諸元を代入し、補正係数Cを求める。次に、求めた補正係数C及び各気筒の運転状態量及びエンジン基本諸元に基づいて、算式1から、各気筒のノッキング限界点火時期Tib~Tifを算出するS20。なお、ガスエンジン10の同型の試験機を運転したときの試験データを予め取得しておく。この試験データに基づいて、補正係数Cを各気筒毎に修正する。修正した補正係数Cに基づいて、各気筒のノッキング限界点火時期Tib~Tifを算出する。図4は、試験機の試験データに基づいて作成した、ノッキング頻度の気筒間ばらつきを示す線図であり、給気温度を一定としたときのデータである。例えば、給気温度Tsの差によってノッキング頻度の気筒間ばらつきが発生するが、実際には各気筒間で給気温度Tsが同一であっても、各気筒の燃焼室壁面温度等の違いによって、ノッキング頻度の気筒間ばらつきが発生する。そのため、気筒毎に正確なノッキング限界点火時期を算出するためには、気筒毎に補正係数Cを修正するとよい。こうして算出された各気筒14a~fのノッキング限界点火時期に基づいて、点火時期設定部40で点火装置における各気筒14a~fの点火時期を設定する。こうして、時々刻々代表気筒14aの筒内圧波形及びその他の運転状態量を検出しながら、各気筒の点火時期をリアルタイムに変更する。本実施形態によれば、選択した代表気筒14aのみに筒内圧センサー16を設ければよく、他の気筒14b~fでは、筒内圧センサーを設けることなく、ノッキング限界点火時期を算出するので、低コスト化できる。そして、筒内圧センサー16で検出した代表気筒14aの筒内圧波形と、各気筒で検出した他の運転状態量及びエンジン基本諸元とから、各気筒毎にノッキング限界点火時期を算出するので、熱効率が良好な運転を実現できる。また、ガスエンジン10と同型の試験機から予め取得した運転試験データに基づいて、補正係数Cを気筒毎に修正するようにしているので、各気筒の正確なノッキング限界点火時期を算出できる。また、代表気筒として、給気温度Tsが最も高く、最もノッキングが起きやすい気筒14aを選択し、気筒14aのノッキング限界点火時期に基づいて、他の気筒のノッキング限界点火時期を算出するようにしているので、全気筒において、安全サイドの運転状態から点火時期を進角し、ノッキング限界点火時期まで調整することが可能であり、ノッキングを確実に抑制できる。