(57)【特許請求の範囲】
【請求項１】
実ノッキング出現率の遅れ演算値が目標出現率となるように点火時期を調整する、リーンバーンで運転されるガスエンジンの制御方法において、
前記ガスエンジンのシリンダ内の最大圧力の遅れ演算値が基準値以上のときは、前記シリンダ内の最大圧力の遅れ演算値を前記目標出現率よりも大きな値以上の仮想ノッキング出現率に換算し、前記実ノッキング出現率の遅れ演算値を前記目標出現率と比較する前に前記仮想ノッキング出現率を前記実ノッキング出現率の遅れ演算値に加算するか、前記実ノッキング出現率の遅れ演算値に代えて前記仮想ノッキング出現率を前記目標出現率と比較する、ガスエンジンの制御方法。
【請求項２】
前記基準値よりも小さな下側設定値と前記基準値の間では、前記シリンダ内の最大圧力の遅れ演算値を、ゼロと前記目標出現率よりも大きな値との間での比例計算により仮想ノッキング出現率に換算し、前記実ノッキング出現率の遅れ演算値を前記目標出現率と比較する前に前記仮想ノッキング出現率を前記実ノッキング出現率の遅れ演算値に加算するか、前記実ノッキング出現率の遅れ演算値に代えて前記仮想ノッキング出現率を前記目標出現率と比較する、請求項１に記載のガスエンジンの制御方法。
【請求項３】
前記シリンダ内の最大圧力の遅れ演算値が前記基準値以上のときも前記基準値未満のときも、前記実ノッキング出現率の遅れ演算値を前記目標出現率と比較する前に前記仮想ノッキング出現率を前記実ノッキング出現率の遅れ演算値に加算し、
前記シリンダ内の最大圧力の遅れ演算値が前記下側設定値以下のときは、前記シリンダ内の最大圧力の遅れ演算値をゼロの仮想ノッキング出現率に換算する、請求項２に記載のガスエンジンの制御方法。
【請求項４】
前記シリンダ内の最大圧力の遅れ演算値が前記基準値よりも大きな上側設定値以上のときは、前記シリンダ内の最大圧力の遅れ演算値を前記目標出現率の１．５倍以上の値の仮想ノッキング出現率に換算する、請求項１〜３のいずれか一項に記載のガスエンジンの制御方法。
【請求項５】
前記シリンダ内の最大圧力の遅れ演算値は、最大圧力の移動平均値である、請求項１〜４のいずれか一項に記載のガスエンジンの制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ガスエンジンの制御方法に関する。
【背景技術】
【０００２】
従来から、ＮＯＸの排出量が少なくかつ高効率、高出力での運転が実現可能なリーンバーンで運転されるガスエンジンが知られている。ガスエンジンの運転領域には、空燃比（空気過剰率）と正味平均有効圧（ＢＭＥＰ）との関係でノッキング領域および失火領域が存在する。このため、ノッキングおよび失火を防止しながら高効率、高出力で運転できるようにガスエンジンを制御することが望まれる。
【０００３】
例えば、特許文献１には、実ノッキング出現率の遅れ演算値が目標出現率となるように点火時期を調整するガスエンジンの制御方法が開示されている。この制御方法では、シリンダ内の圧力から１サイクルごとに許容ノッキングが出現したか否かが判定される。前記の実ノッキング出現率は、所定のサイクル数に対する許容ノッキングが出現したサイクル数の割合である。そして、この制御方法では、実ノッキング出現率の遅れ演算値（例えば、移動平均値）が目標出現率よりも低い場合はそれらの偏差に基づいて点火時期が進角され、実ノッキング出現率の遅れ演算値が目標出現率よりも高い場合はそれらの偏差に基づいて点火時期が遅角される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−８４６８１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、燃料ガスとしてメタン価が高い燃料ガスが用いられる場合は、ノッキングが殆ど発生しない。従って、上記のように実ノッキング出現率の遅れ演算値に基づいて点火時期を制御した場合には、実ノッキング出現率の遅れ演算値がほぼゼロとなり、点火時期が進角され過ぎることになる。
【０００６】
そこで、本発明は、燃料ガスとしてメタン価が高い燃料ガスが用いられる場合でも、適切な点火時期での運転が可能なガスエンジンの制御方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
前記課題を解決するために、本発明の発明者らは、鋭意研究の結果、点火時期が進角され過ぎるとシリンダ内の最大圧力が過度に高くなるため、シリンダ内の最大圧力を仮想的にノッキング出現率に換算することを思い付いた。このような仮想的なノッキング出現率を用いれば、上記の実ノッキング出現率の遅れ演算値に基づいて点火時期を制御する制御方法を殆ど変更することなく使用することができる。本発明は、このような観点からなされたものである。
【０００８】
すなわち、本発明のガスエンジンの制御方法は、実ノッキング出現率の遅れ演算値が目標出現率となるように点火時期を調整する、リーンバーンで運転されるガスエンジンの制御方法において、前記ガスエンジンのシリンダ内の最大圧力の遅れ演算値が基準値以上のときは、前記シリンダ内の最大圧力の遅れ演算値を前記目標出現率よりも大きな値以上の仮想ノッキング出現率に換算し、前記実ノッキング出現率の遅れ演算値を前記目標出現率と比較する前に前記仮想ノッキング出現率を前記実ノッキング出現率の遅れ演算値に加算するか、前記実ノッキング出現率の遅れ演算値に代えて前記仮想ノッキング出現率を前記目標出現率と比較する、ことを特徴とする。
【０００９】
上記の構成によれば、シリンダ内でノッキングが殆ど発生せずに、点火時期が進角され過ぎてシリンダ内の最大圧力が過度に高くなったときには、シリンダ内の最大圧力が目標出現率よりも大きな値以上の仮想ノッキング出現率に換算され、点火時期が遅角される。従って、燃料ガスとしてメタン価が高い燃料ガスが用いられる場合でも、適切な点火時期でガスエンジンを運転することができる。
【００１０】
前記基準値よりも小さな下側設定値と前記基準値の間では、前記シリンダ内の最大圧力の遅れ演算値を、ゼロと前記目標出現率よりも大きな値との間での比例計算により仮想ノッキング出現率に換算し、前記実ノッキング出現率の遅れ演算値を前記目標出現率と比較する前に前記仮想ノッキング出現率を前記実ノッキング出現率の遅れ演算値に加算するか、前記実ノッキング出現率の遅れ演算値に代えて前記仮想ノッキング出現率を前記目標出現率と比較してもよい。この構成によれば、仮想ノッキング出現率がゼロから徐々に大きくなる。従って、仮想ノッキング出現率をゼロから急激に立ち上げたときに比べ、安定性を向上させることができる。
【００１１】
前記シリンダ内の最大圧力の遅れ演算値が前記基準値以上のときも前記基準値未満のときも、前記実ノッキング出現率の遅れ演算値を前記目標出現率と比較する前に前記仮想ノッキング出現率を前記実ノッキング出現率の遅れ演算値に加算し、前記シリンダ内の最大圧力の遅れ演算値が前記下側設定値以下のときは、前記シリンダ内の最大圧力の遅れ演算値をゼロの仮想ノッキング出現率に換算してもよい。この構成によれば、常に仮想ノッキング出現率を実ノッキング出現率に加算することができ、制御ロジックを簡単にすることができる。
【００１２】
前記シリンダ内の最大圧力の遅れ演算値が前記基準値よりも大きな上側設定値以上のときは、前記シリンダ内の最大圧力の遅れ演算値を前記目標出現率の１．５倍以上の値の仮想ノッキング出現率に換算してもよい。この構成によれば、シリンダ内の最大圧力が上側設定値（例えば、機械的な上限）近くになったときは、点火時期が大きく遅角される。従って、たとえ好ましくない状況になったとしてもそれを早急に改善することができる。
【００１３】
例えば、前記シリンダ内の最大圧力の遅れ演算値は、最大圧力の移動平均値であってもよい。
【発明の効果】
【００１４】
本発明によれば、燃料ガスとしてメタン価が高い燃料ガスが用いられる場合でも、適切な点火時期でガスエンジンを運転することができる。
【図面の簡単な説明】
【００１５】
【図１】ガスエンジンの一部の断面図である。
【図２】本発明の一実施形態に係るガスエンジンの制御方法におけるフローチャートである。
【図３】シリンダ内の圧力波形の一例である。
【図４】実ノッキング出現率とその移動平均値を説明するための図である。
【図５】シリンダ内の最大圧力の移動平均値を仮想ノッキング出現率に換算するためのグラフである。
【図６】変形例のガスエンジンの制御方法におけるフローチャートである。
【図７】別の変形例のガスエンジンの制御方法におけるフローチャートである。
【発明を実施するための形態】
【００１６】
図１〜図５を参照して、本発明の一実施形態に係るガスエンジンの制御方法を説明する。図１に、本実施形態で用いられるガスエンジン１を示す。このガスエンジン１は、リーンバーンで運転される。
【００１７】
ガスエンジン１は、燃料ガス（例えば、天然ガス）のみを燃焼させるガス専焼エンジンであってもよいし、燃料ガスと燃料油の一方または双方を燃焼させる二元燃料エンジンであってもよい。本実施形態では、ガスエンジン１が４ストロークエンジンであるが、ガスエンジン１は２ストロークエンジンであってもよい。
【００１８】
ガスエンジン１は、複数のシリンダ２１を有する（図１では１つのシリンダ２１のみを図示）。全てのシリンダ２１の構成は同じであるため、以下では１つのシリンダ２１のみについて説明する。
【００１９】
シリンダ２１内にはピストン２２が往復動自在に配設されており、シリンダ２１およびピストン２２によって燃焼室２が形成されている。ピストン２２は、図略の連接棒により図略のクランク軸と連結されている。
【００２０】
４ストロークエンジンの場合、シリンダ２１において、ピストン２２が二往復することにより、ガスエンジン１の１サイクル（給気、圧縮、膨張、排気）が行われる。シリンダ２１における１サイクルの間のガスエンジン１の位相角（０〜７２０度）は、位相角検出器６２により検出される。位相角としては、クランク軸の回転角（クランク角）やピストン２２の位置などを用いることができる。例えば、位相角検出器６２は、電磁ピックアップ、近接スイッチまたはロータリーエンコーダである。また、位相角検出器６２からは、ガスエンジン１の実回転数も検出される。
【００２１】
燃焼室２へは、給気路４１を通じて空気が供給され、燃焼室２からは、排気路４２を通じて排気が排出される。給気路４１には、燃焼室２へ供給される空気中に燃料ガスを噴射する主燃料噴射弁５１が設けられている。
【００２２】
シリンダ２１には、給気路４１の燃焼室２に対する開口である給気ポートを開閉する給気弁２３と、排気路４２の燃焼室２に対する開口である排気ポートを開閉する排気弁２４が設けられている。また、シリンダ２１には、燃焼室２内で燃料ガスと空気の混合気に点火するための点火プラグ（点火装置の一例）５５が設けられている。
【００２３】
本実施形態では、燃焼室２が、給気路４１および排気路４２と連通する主燃焼室２Ａと、連通孔が設けられた隔壁２５により主燃焼室２Ａと隔てられた副燃焼室２Ｂとからなる。点火プラグ５５は副燃焼室２Ｂに配置されており、副燃焼室２Ｂ内には副燃料噴射弁５２から燃料ガスが噴射される。副燃料噴射弁５２からの燃料ガスの噴射により副燃焼室２Ｂ内にはリッチな混合気が形成され、この混合気が点火プラグ５５により点火される。これにより副燃焼室２Ｂ内に火炎が発生し、その火炎が隔壁２５の連通孔を通じて主燃焼室２Ａ内に伝播することにより主燃焼室２Ａ内のリーンな混合気にも点火される。シリンダ２１には、当該シリンダ２１内の圧力Ｐ（主燃焼室２Ａ内の圧力）を検出する圧力計６１が設けられている。
【００２４】
ただし、点火装置は、副燃焼室２Ｂ内の混合気に点火する点火プラグ５５に限られない。例えば、点火装置としては、主燃焼室２Ａ内に高圧のパイロット燃料（オイルや燃料ガス）を直接的に噴射することにより当該パイロット燃料を自己発火させるパイロット燃料噴射弁を採用することも可能である。
【００２５】
燃料噴射弁５１，５２および点火プラグ５５は、位相角検出器６２で検出される位相角などに基づいて制御装置３により制御される。制御装置３は、実ノッキング出現率Ｋｒの遅れ演算値が目標出現率Ｋｔ（例えば、７．５％）となるように、１サイクルごとに点火時期を調整する。
【００２６】
実ノッキング出現率Ｋｒは、所定のサイクル数Ｃｙ（例えば、５０サイクル）に対する許容ノッキングが出現したサイクル数Ｃｋの割合である（Ｋｒ＝Ｃｋ／Ｃｙ）。本実施形態では、実ノッキング出現率Ｋｒの遅れ演算値が、実ノッキング出現率Ｋｒの移動平均値である。ただし、遅れ演算値は、一次遅れ演算値であってもよい。
【００２７】
以下、図２のフローチャートを参照して、制御装置３により実行される制御方法について詳しく説明する。
【００２８】
制御装置３は、圧力計６１で検出されるシリンダ２１内の圧力Ｐに基づいて、１サイクルごとに、燃焼状態が、許容不可ノック（Ｄ）、許容ノック（Ｃ）、通常（Ａ）、失火（Ｂ）のいずれであったかを判定する（ステップＳ１）。
【００２９】
燃焼状態の判定の仕方としては、例えば次のような方法がある。制御装置３は、図３に示すシリンダ２１内の圧力波形をフィルタに通して、ピストン２２が上死点に達してから所定時間ΔＴ内の高周波成分を抽出する。さらに、制御装置３は、抽出した高周波成分のうちから複数個の高周波成分をサンプリングし、サンプリングされた高周波成分の平均値ＰＡを算出する。平均値ＰＡが第１閾値γ１以上の場合（γ１≦ＰＡ）は燃焼状態が許容不可ノックであったと判定され、平均値ＰＡが第１閾値γ１未満かつ第２閾値γ２以上の場合（γ２≦ＰＡ＜γ１）は燃焼状態が許容ノックであったと判定される。また、制御装置３は、上死点前後のシリンダ２１内の圧力Ｐの偏差ΔＰを算出する。偏差ΔＰが第３閾値γ３未満の場合（ΔＰ＜γ３）には燃焼状態が失火であったと判定される。偏差ΔＰが第３閾値γ３以上であって上述した平均値ＰＡが第２閾値γ２未満の場合（γ３≦ΔＰ、ＰＡ＜γ２）には燃焼状態が通常であったと判定される。
【００３０】
図４に示すように、現在のサイクルをｎ回目とすると、制御装置３は、現在のサイクルからＣｙ回前までの燃焼状態を記憶している。制御装置３は、現在のサイクルからＣｙ回前までで許容ノッキング（Ｃ）が出現したサイクル数Ｃｋを求め、ＣｋをＣｙで割ることによって、現在のサイクルにおける実ノッキング出現率Ｋｒnを算出する（ステップＳ２）。なお、算出された実ノッキング出現率Ｋｒは、図４に示すように、現在のサイクルから所定サイクル数Ｎ分だけ制御装置３に格納される。
【００３１】
次に、制御装置３は、実ノッキング出現率Ｋｒの移動平均値Ｋｒ(AVE)を算出する（ステップＳ３）。より詳しくは、図４に示すように、現在のサイクルをｎ回目とすると、制御装置は、現在のサイクルとＮ回前のサイクルとの間で実ノッキング出現率Ｋｒを平均化する（Ｋｒ(AVE)＝（Ｋｒn＋Ｋｒn-1＋・・・＋Ｋｒn-N+1）／Ｎ）。
【００３２】
ガスエンジン１が陸上での発電に用いられる場合などでガスエンジン１の負荷が一定もしくはその変動が小さい場合は、Ｎは比較的に多いことが望ましい（例えば、５０サイクル）。一方、ガスエンジン１が船舶の主機に用いられる場合などでガスエンジン１の負荷が頻繁に変動する場合は、Ｎは比較的に少ないことが望ましい（例えば、１０サイクル）。
【００３３】
次に、制御装置３は、シリンダ２１内の最大圧力Ｐｍの遅れ演算値として、最大圧力Ｐｍの移動平均値Ｐｍ(AVE)を算出する（ステップＳ４）。ただし、遅れ演算値は、一次遅れ演算値であってもよい。制御装置３には、現在のサイクルからＭ回前までの最大圧力Ｐｍが格納されている（例えば、Ｍは５サイクル）。制御装置３は、現在のサイクルをｎ回目とすると、現在のサイクルとＭ回前のサイクルとの間で最大圧力Ｐｍを平均化する（Ｐｍ(AVE)＝（Ｐｍn＋Ｐｍn-1＋・・・＋Ｐｍn-M+1）／Ｍ）。
【００３４】
その後、制御装置３は、シリンダ２１内の最大圧力Ｐｍの移動平均値Ｐｍ(AVE)を仮想ノッキング出現率Ｋｖに換算する（ステップＳ５）。具体的に、制御装置３は、図５に示すように、移動平均値Ｐｍ(AVE)が基準値Ｐ０（例えば、２０ＭＰａ）以上のときは、移動平均値Ｐｍ(AVE)を目標出現率Ｋｔよりも大きな値Ｒ１（例えば、９．０％）以上の仮想ノッキング出現率Ｋｖに換算する。
【００３５】
本実施形態では、制御装置３が、基準値Ｐ０よりも小さな下側設定値Ｐ１と基準値Ｐ０の間では、移動平均値Ｐｍ(AVE)を、ゼロとＲ１との間での比例計算により仮想ノッキング出現率Ｋｖに換算する。換言すれば、仮想ノッキング出現率Ｋｖは、Ｐ１とＰ０の間ではゼロからＲ１まで一定の傾きで立ち上がる。例えば、下側設定値Ｐ１は、基準値Ｐ０よりも１〜２ＭＰａ小さい（例えば、１８．５ＭＰａ）。また、制御装置３は、移動平均値Ｐｍ(AVE)が下側設定値Ｐ１以下のときは、移動平均値Ｐｍ(AVE)をゼロの仮想ノッキング出現率Ｋｖに換算する。
【００３６】
さらに、本実施形態では、制御装置３が、移動平均値Ｐｍ(AVE)が基準値Ｐ０よりも大きな上側設定値Ｐ２以上のときは、移動平均値Ｐｍ(AVE)を目標出現率Ｋｔの１．５倍以上の値Ｒ２（例えば、２０％）の仮想ノッキング出現率Ｋｖに換算する。例えば、上側設定値Ｐ２は、基準値Ｐ０よりも２〜４ＭＰａ大きい（例えば、２３ＭＰａ）。
【００３７】
次に、制御装置３は、シリンダ２１内の最大圧力Ｐｍの移動平均値Ｐｍ(AVE)が基準値Ｐ０以上のときも基準値Ｐ０未満のときも、移動平均値Ｋｒ(AVE)に仮想ノッキング出現率Ｋｖを加算して算出出現率Ｋｃを算出する（ステップＳ６）。その後、制御装置３は、算出出現率Ｋｃ（仮想ノッキング出現率Ｋｖが加算された移動平均値Ｋｒ(AVE)）を目標出現率Ｋｔと比較する（ステップＳ７）。換言すれば、実ノッキング出現率Ｋｒの移動平均値Ｋｒ(AVE)が目標出現率Ｋｔと比較される前に、移動平均値Ｋｒ(AVE)に仮想ノッキング出現率Ｋｖが加算される。
【００３８】
算出出現率Ｋｃが目標出現率Ｋｔよりも低い場合（ステップＳ７でＹＥＳ）、制御装置３は、それらの偏差ΔＫ（＝Ｋｔ−Ｋｃ）に基づいて点火タイミングを進角させる（ステップＳ８）。一方、算出出現率Ｋｃが目標出現率Ｋｔよりも高い場合（ステップＳ７でＮＯ）、制御装置３は、それらの偏差ΔＫ（＝Ｋｃ−Ｋｔ）に基づいて点火タイミングを遅角させる（ステップＳ９）。なお、図２では、Ｋｃ＝Ｋｔの場合はステップＳ８に進むようになっているが、Ｋｃ＝Ｋｔの場合はステップＳ９に進んでもよい。その後は、ステップＳ１に戻る。
【００３９】
以上説明したように、本実施形態の制御方法では、シリンダ２１内でノッキングが殆ど発生せずに、点火時期が進角され過ぎてシリンダ２１内の最大圧力Ｐｍが過度に高くなったときには、シリンダ２１内の最大圧力Ｐｍが目標出現率Ｋｔよりも大きな値Ｒ１以上の仮想ノッキング出現率Ｋｖに換算され、点火時期が遅角される。従って、燃料ガスとしてメタン価が高い燃料ガスが用いられる場合でも、適切な点火時期でガスエンジン１を運転することができる。
【００４０】
また、本実施形態では、シリンダ２１内の最大圧力Ｐｍの移動平均値Ｐｍ(AVE)が基準値Ｐ０以上のときも基準値Ｐ０未満のときも算出出現率Ｋｃが算出されるため、常に仮想ノッキング出現率Ｋｖを実ノッキング出現率Ｋｒに加算することができ、制御ロジックを簡単にすることができる。
【００４１】
さらに、本実施形態では、シリンダ２１内の最大圧力Ｐｍの移動平均値Ｐｍ(AVE)が上側設定値Ｐ２以上のときは、移動平均値Ｐｍ(AVE)が目標出現率Ｋｔの１．５倍以上の値Ｒ２の仮想ノッキング出現率Ｋｖに換算されるため、シリンダ２１内の最大圧力Ｐｍが上側設定値Ｐ２（例えば、機械的な上限）近くになったときは、点火時期が大きく遅角される。従って、たとえ好ましくない状況になったとしてもそれを早急に改善することができる。
【００４２】
（変形例）
本発明は上述した実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々の変形が可能である。
【００４３】
例えば、前記実施形態では、図５に示すように下側設定値Ｐ１と基準値Ｐ０の間では、仮想ノッキング出現率ＫｖがゼロからＲ１まで一定の傾きで立ち上がっている。しかし、仮想ノッキング出現率Ｋｖは、シリンダ２１内の最大圧力Ｐｍの移動平均値Ｐｍ(AVE)が基準値Ｐ０となったときに、ゼロからＲ１まで垂直に立ち上がってもよい。ただし、前記実施形態のように、仮想ノッキング出現率Ｋｖがゼロから徐々に大きくなれば、仮想ノッキング出現率Ｋ０をゼロから急激に立ち上げたときに比べ、安定性を向上させることができる。
【００４４】
また、前記実施形態では、シリンダ２１内の最大圧力Ｐｍの移動平均値Ｐｍ(AVE)が下側設定値Ｐ１以下のときは、移動平均値Ｐｍ(AVE)をゼロの仮想ノッキング出現率Ｋｖに換算しているが、図６に示すように、ステップＳ５の前に、移動平均値Ｐｍ(AVE)を下側設定値Ｐ１と比較するステップＳ１０を設け、移動平均値Ｐｍ(AVE)が下側設定値Ｐ１以上のときにのみ移動平均値Ｐｍ(AVE)を仮想ノッキング出現率Ｋｖに換算してもよい。
【００４５】
この場合、移動平均値Ｐｍ(AVE)が下側設定値Ｐ１未満のときは（ステップＳ１０でＮＯ）、実ノッキング出現率Ｋｒの移動平均値Ｋｒ(AVE)を目標出現率Ｋｔと直接的に比較すればよい（ステップＳ１１）。
【００４６】
あるいは、移動平均値Ｐｍ(AVE)が下側設定値Ｐ１以上のときは、実ノッキング出現率Ｋｒの移動平均値Ｋｒ(AVE)がほぼゼロであると想定されるので、図６のステップＳ６およびステップＳ７に代えて、図７に示すように、仮想ノッキング出現率Ｋｖを目標出現率Ｋｔと比較するステップＳ１２を採用してもよい。換言すれば、実ノッキング出現率Ｋｒの移動平均値Ｋｒ(AVE)に代えて、仮想ノッキング出現率Ｋｖが目標出現率Ｋｔと比較されてもよい。そして、仮想ノッキング出現率Ｋｖが目標出現率Ｋｔよりも低い場合はそれらの偏差ΔＫ（＝Ｋｔ−Ｋｖ）に基づいて点火時期を進角し、仮想ノッキング出現率Ｋｖが目標出現率Ｋｔよりも高い場合はそれらの偏差Δｋ（＝Ｋｖ−Ｋｔ）に基づいて点火時期を遅角する。
【００４７】
なお、シリンダ２１内の最大圧力Ｐｍの移動平均値Ｐｍ(AVE)が基準値Ｐ０となったときに仮想ノッキング出現率ＫｖをゼロからＲ１まで垂直に立ち上げる場合は、図６および図７に示すステップＳ１０において、移動平均値Ｐｍ(AVE)を、下側設定値Ｐ１に代えて基準値Ｐ０と比較してもよい。
【符号の説明】
【００４８】
１    ガスエンジン
２１  シリンダ
Ｐｍ  最大圧力
Ｐｒ(AVE)  移動平均値
Ｋｒ  実ノッキング出現率
Ｋｒ(AVE)  移動平均値（遅れ演算値）
Ｋｖ  仮想ノッキング出現率
Ｋｃ  算出出現率
ΔＫ  偏差
【要約】      （修正有）【課題】燃料ガスとしてメタン価が高い燃料ガスが用いられる場合でも、適切な点火時期での運転が可能なガスエンジンの制御方法を提供する。【解決手段】実ノッキング出現率の遅れ演算値が目標出現率となるように点火時期を調整する、リーンバーンで運転されるガスエンジンの制御方法において、ガスエンジンのシリンダ内の最大圧力の遅れ演算値が基準値以上のときは、遅れ演算値を目標出現率よりも大きな値以上の仮想ノッキング出現率に換算しＳ５、実ノッキング出現率の遅れ演算値を目標出現率と比較する前に仮想ノッキング出現率を実ノッキング出現率の遅れ演算値に加算するＳ６ことにより、メタン価が高い燃料ガスが用いられシリンダ内でノッキングが殆ど発生しない場合でも、点火時期が進角され過ぎてシリンダ内の最大圧力Ｐｍが過度に高くなったときには、点火時期が遅角される。【選択図】図２【図１】
0005990316.tif 000002
【図２】
0005990316.tif 000003
【図３】
0005990316.tif 000004
【図４】
0005990316.tif 000005
【図５】
0005990316.tif 000006
【図６】
0005990316.tif 000007
【図７】
0005990316.tif 000008