JP2022191833A - メタン価の推定装置、ガスエンジンの制御装置およびメタン価の推定方法 - Google Patents

Abstract

【課題】ガスエンジンの運転時に、特別なセンサを用いることなくメタン価を推定可能なメタン価の推定装置、ガスエンジンの制御装置およびメタン価の推定方法を提供する。
【解決手段】対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定装置であって、ガスエンジンにおける、点火時期とノッキング強度と供給されるガス燃料のメタン価とが予め関連付けられた第１の関連付け情報、を取得する第１の関連付け情報取得部と、対象ガスエンジンにおける点火時期である対象点火時期を取得する点火時期取得部と、対象ガスエンジンの対象点火時期を含むサイクルにおけるノッキング強度である対象ノッキング強度を取得するノッキング強度取得部と、第１の関連付け情報に基づいて、対象点火時期および対象ノッキング強度から、対象ガスエンジンにおける気筒に供給されたガス燃料のメタン価を推定するメタン価推定部と、を備える。
【選択図】 図３

Description

本開示は、ガスエンジンに供給されるガス燃料のメタン価の推定装置、該メタン価の推定装置を備えるガスエンジンの制御装置、およびメタン価の推定方法に関する。

例えば、天然ガスの産地ごとにガス性状（例えば、メタン価）が異なることが知られている。従来のガスエンジンの運転状態を制御する一般的な制御装置は、ガスエンジンに供給されるガス燃料のガス性状（例えば、メタン価）の変化を捉えることができなかった。このため、ガス燃料の切り替え時などのような、ガスエンジンに供給されるガス燃料のガス性状が変化したときに、ガスエンジンの適切な運転領域から外れて、ノッキングや失火などの異常燃焼が生じる虞がある。

上記ガスエンジンにおいて、メタン価の異なる複数のガス燃料を使用する場合には、複数のガス燃料が混同しないように、ガス燃料の供給系統に弁を設けて、複数のガス燃料のうちの１つのガス燃料のみをガスエンジンに供給することが行われていた。ガス燃料を切り替える際には、ガスエンジンを停止し、新たなガス燃料に対応するように各種パラメータを変更した後に、ガスエンジンを再起動させていた。ガスエンジンを断続運転することは効率上好ましくはないため、ガス燃料を切り替える際にもガスエンジンを連続運転させることが望まれている。

ガス燃料を切り替える際にもガスエンジンを連続運転させるための方策として、ガス燃料を切り替える前にガスエンジンにおけるガス燃料の点火時期を安全側に変更することが考えられる。この場合には、ガスエンジンの出力や燃費性能の低下を招く虞がある。よって、ガスエンジンに供給されるガス燃料のガス性状（例えば、メタン価）の変化を取得することが望まれている。

特許第６００２２３５号公報

特許文献１には、気筒に供給されるガス燃料の発熱量を検出するカロリーメータや、気筒に供給されるガス燃料の組成を検出するガスクロマトグラフィなどの検出値に基づいてメタン価を算出し、算出したメタン価に基づいて点火時期を変更することが開示されている。カロリーメータやガスクロマトグラフィなどの特別なセンサは、通常のガスエンジンに備えられた機器ではないため、特別なセンサの設置や管理のためにガスエンジンの製造コストや管理コストの増大化を招く虞がある。

上述した事情に鑑みて、本開示の少なくとも一実施形態の目的は、ガスエンジンの運転時に、特別なセンサを用いることなくメタン価を推定可能なメタン価の推定装置、ガスエンジンの制御装置およびメタン価の推定方法を提供することにある。

本開示の一実施形態にかかるメタン価の推定装置は、
少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定装置であって、
ガスエンジンにおける、点火時期とノッキング強度と供給されるガス燃料のメタン価とが予め関連付けられた第１の関連付け情報、を取得する第１の関連付け情報取得部と、
前記対象ガスエンジンにおける前記点火時期である対象点火時期を取得する点火時期取得部と、
前記対象ガスエンジンの前記対象点火時期を含むサイクルにおける前記ノッキング強度である対象ノッキング強度を取得するノッキング強度取得部と、
前記第１の関連付け情報に基づいて、前記対象点火時期および前記対象ノッキング強度から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定部と、を備える。

本開示の一実施形態にかかるガスエンジンの制御装置は、
前記メタン価の推定装置と、
前記メタン価の推定装置で推定された前記ガス燃料のメタン価に応じて、前記対象ガスエンジンの点火時期である対象点火時期を制御するように構成された点火時期制御部と、を備える。

本開示の一実施形態にかかるメタン価の推定方法は、
少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定方法であって、
ガスエンジンにおける、点火時期とノッキング強度と供給されるガス燃料のメタン価とが予め関連付けられた第１の関連付け情報、を取得する第１の関連付け情報取得ステップ部と、
前記対象ガスエンジンにおける前記点火時期である対象点火時期を取得する点火時期取得ステップと、
前記対象ガスエンジンの前記対象点火時期を含むサイクルにおける前記ノッキング強度である対象ノッキング強度を取得するノッキング強度取得ステップと、
前記第１の関連付け情報に基づいて、前記対象点火時期および前記対象ノッキング強度から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定ステップと、を備える。

本開示の少なくとも一実施形態によれば、ガスエンジンの運転時に、特別なセンサを用いることなくメタン価を推定可能なメタン価の推定装置、ガスエンジンの制御装置およびメタン価の推定方法が提供される。

本開示の一実施形態にかかるガスエンジンの制御装置を備えるエンジンシステムの概略構成図である。 本開示の一実施形態におけるガス燃料供給系統の概略構成図である。 本開示の第１の実施形態にかかるメタン価の推定方法のフロー図である。 一気筒におけるノッキング強度と点火時期との関係を説明するための説明図である。 特定のノッキング強度を用いたメタン価の推定方法を説明するための説明図である。 本開示の第２の実施形態にかかるメタン価の推定方法のフロー図である。 ガス燃料の産地毎のメタン価と発熱量との関係を示す図である。 ガス燃料毎のガスエンジンにおけるエンジン負荷とガスエンジンに供給されるガス燃料のガス供給量との関係を示す図である。 本開示の第３の実施形態にかかるメタン価の推定方法のフロー図である。 クランク角に対する筒内圧の変化を示す図である。 低位発熱量とガス燃料のメタン価との関係を示す図である。 ガスエンジンに供給されるガス燃料のメタン価とガスエンジンの点火時期との関係を説明するための説明図である。

以下、添付図面を参照して本開示の幾つかの実施形態について説明する。ただし、実施形態として記載されている又は図面に示されている構成部品の寸法、材質、形状、その相対的配置等は、本開示の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。
例えば、「ある方向に」、「ある方向に沿って」、「平行」、「直交」、「中心」、「同心」或いは「同軸」等の相対的或いは絶対的な配置を表す表現は、厳密にそのような配置を表すのみならず、公差、若しくは、同じ機能が得られる程度の角度や距離をもって相対的に変位している状態も表すものとする。
例えば、「同一」、「等しい」及び「均質」等の物事が等しい状態であることを表す表現は、厳密に等しい状態を表すのみならず、公差、若しくは、同じ機能が得られる程度の差が存在している状態も表すものとする。
例えば、四角形状や円筒形状等の形状を表す表現は、幾何学的に厳密な意味での四角形状や円筒形状等の形状を表すのみならず、同じ効果が得られる範囲で、凹凸部や面取り部等を含む形状も表すものとする。
一方、一の構成要素を「備える」、「含む」、又は、「有する」という表現は、他の構成要素の存在を除外する排他的な表現ではない。
なお、同様の構成については同じ符号を付し説明を省略することがある。

（エンジンシステム）
図１は、本開示の一実施形態にかかるガスエンジンの制御装置を備えるエンジンシステムの概略構成図である。
幾つかの実施形態にかかるメタン価の推定装置１（１Ａ、１Ｂ、１Ｃ）は、少なくとも１つの気筒２１を有するガスエンジン２である対象ガスエンジン２ｔに供給されるガス燃料のメタン価である対象メタン価ＭＮｔを推定する装置である。ここで、対象ガスエンジン２ｔは、メタン価の推定装置１により供給されるガス燃料のメタン価が推定される対象のガスエンジン２である。以下、対象ガスエンジン２ｔを対象とするパラメータなどについては、対象ガスエンジン２ｔ以外のガスエンジン２と区別するために「対象」の語や符号ｔを付すことがある。

図示される実施形態では、メタン価の推定装置１は、ガスエンジン２の運転制御や燃焼制御を行う制御装置３に搭載されている。エンジンシステム４は、図１に示されるように、対象ガスエンジン２ｔと、対象ガスエンジン２ｔの制御装置３である対象制御装置３ｔと、対象ガスエンジン２ｔにおける気筒２１の内部にガス燃料を噴射するように構成されたガス燃料噴射装置（図示例では、燃焼噴射弁）４１と、ガス燃料噴射装置４１にガス燃料を供給するためのガス燃料供給系統５（図２参照）と、対象ガスエンジン２ｔにおける気筒２１の内部のガス燃料を含む気体（混合気）に点火する点火装置（図示例では、点火プラグ）４２と、を備える。

図示される実施形態では、対象ガスエンジン２ｔは、複数の気筒２１を有する。ガス燃料噴射装置４１や点火装置４２は、気筒２１毎に個別に設けられる。

対象ガスエンジン２ｔ（２）は、内部でガス燃料を燃焼させることで動力を発生させるように構成されている。図示される実施形態では、エンジンシステム４は、対象ガスエンジン２ｔの駆動シャフト２２に接続された発電機４３をさらに備える。発電機４３は、駆動シャフト２２を介して対象ガスエンジン２ｔが発生させた動力が伝達され、対象ガスエンジン２ｔから伝達された動力により発電するように構成されている。

（エンジンシステムに搭載される測定機器）
エンジンシステム４は、対象ガスエンジン２ｔ（２）の気筒２１毎に設けられるノッキングセンサ４４と、対象ガスエンジン２ｔ（２）の気筒２１毎に設けられる筒内圧センサ４５と、を備える。ノッキングセンサ４４や筒内圧センサ４５は、通常のガスエンジンに搭載される測定機器である。エンジンシステム４には、カロリーメータやガスクロマトグラフィなどの特別なセンサは搭載されていない。

ノッキングセンサ４４は、装着された気筒２１の振動に応じた電気を発生させるように構成されている。或る実施形態では、ノッキングセンサ４４は、装着された気筒２１の振動に応じてセンサに内蔵された重りが振動し、重りの振動により圧電素子（圧電セラミックス）に力が加わることで電気を発生させる。ノッキング発生時には通常時よりも振動が大きくなるため、圧電素子が発生させる電気力が増大する。ノッキングセンサ４４で生じる電気量を対象制御装置３ｔ（３）でモニタリングすることで、ノッキングセンサ４４が装着された気筒２１におけるノッキングの発生の有無を検知できる。

筒内圧センサ４５は、装着された気筒２１の内部の圧力に応じた電気を発生させるように構成されている。或る実施形態では、筒内圧センサ４５は、圧電素子や歪ゲージなどを圧力検出素子として用いる。筒内圧センサ４５で生じる電気量を対象制御装置３ｔ（３）でモニタリングすることで、筒内圧センサ４５が装着された気筒２１の内部の圧力を監視できる。

（ガスエンジンの制御装置）
対象制御装置３ｔ（３）は、対象ガスエンジン２ｔやガス燃料噴射装置４１、点火装置４２などのエンジンシステム４が備える各装置の運転を制御するためのエンジンコントロールユニットからなる。対象制御装置３ｔ（３）は、図１に示されるように、エンジン制御部３１と、燃焼制御部３２と、燃焼診断部３３と、ガス燃料供給量制御部３４と、点火時期制御部３５と、上記メタン価の推定装置１と、を備える。

例えば、メタン価の推定装置１などの対象制御装置３ｔ（３）の各部には、対象ガスエンジン２ｔなどのエンジンシステム４が備える各装置から対象ガスエンジン２ｔの運転に関する情報が送られるようになっている。上記対象ガスエンジン２ｔの運転に関する情報には、例えば、対象ガスエンジン２ｔの回転数、エンジン（発電機４３の）出力、エンジン（発電機４３の）負荷、入口圧力、出口圧力、入口温度、出口温度、若しくは圧力比、気筒２１のクランク角、ノッキングセンサ４４の検出値、又は筒内圧センサ４５の検出値などが含まれる。

燃焼診断部３３は、ノッキングセンサ４４の検出値、筒内圧センサ４５の検出値などの対象ガスエンジン２ｔの運転に関する情報に基づいて、対象ガスエンジン２ｔにおける各気筒２１の燃焼診断を実行可能に構成されている。燃焼制御部３２は、燃焼診断部３３から送られた燃焼診断の結果や対象制御装置３ｔ（３）の各部から送られる情報に基づいて、対象ガスエンジン２ｔにおける各気筒２１に供給されるガス燃料の供給量や対象ガスエンジン２ｔにおける各気筒２１の点火時期を決定するように構成されている。

ガス燃料供給量制御部３４は、燃焼制御部３２において決定されたガス燃料の供給量が、対象ガスエンジン２ｔにおける各気筒２１に供給されるように、ガス燃料供給系統５に設けられた弁やガス燃料噴射装置４１の開閉を制御するように構成されている。点火時期制御部３５は、燃焼制御部３２において決定された点火時期に点火装置４２が点火するように、点火装置４２の点火時期を制御するように構成されている。

エンジン制御部３１は、対象ガスエンジン２ｔの運転に関する情報や対象制御装置３ｔ（３）の各部から送られる情報に基づいて、対象ガスエンジン２ｔの運転や停止を制御するように構成されている。

（ガス燃料供給系統）
図２は、本開示の一実施形態におけるガス燃料供給系統の概略構成図である。図２に示されるように、ガス燃料供給系統５は、第１のガス燃料を貯留する第１のガス燃料貯留装置（図示例では、ガス燃料貯留タンク）５１と、第１のガス燃料とはメタン価が異なる第２のガス燃料を貯留する第２のガス燃料貯留装置（図示例では、ガス燃料貯留タンク）５２と、第１のガス燃料貯留装置５１から対象ガスエンジン２ｔ（具体的には、対象ガスエンジン２ｔにおけるガス燃料噴射装置４１）に第１のガス燃料を送るための第１のガス燃料供給ライン５３と、第２のガス燃料貯留装置５２から対象ガスエンジン２ｔ（具体的には、対象ガスエンジン２ｔにおけるガス燃料噴射装置４１）に第２のガス燃料を送るための第２のガス燃料供給ライン５４と、対象ガスエンジン２ｔに送られるガス燃料の供給元を、第１のガス燃料貯留装置５１又は第２のガス燃料貯留装置５２に切り替え可能に構成され切り替え装置５５と、を備える。

図示される実施形態では、第２のガス燃料供給ライン５４は、合流部５６において第１のガス燃料供給ライン５３に合流している。切り替え装置５５は、第１のガス燃料供給ライン５３の合流部５６よりも上流側（第１のガス燃料貯留装置５１側）に設けられる第１の弁５５１と、第２のガス燃料供給ライン５４の合流部５６よりも上流側（第２のガス燃料貯留装置５２側）に設けられる第２の弁５５２と、を含む。第１の弁５５１および第２の弁５５２の夫々は、ガス燃料の供給路を開閉する弁体を可動させることで、上記弁体よりも下流側（対象ガスエンジン２ｔ側）に供給されるガス燃料の流量を調整可能である。

図２に示される実施形態では、第１のガス燃料貯留装置５１に貯留されるガス燃料は、第１の液化ガスからなる。第２のガス燃料貯留装置５２に貯留されるガス燃料は、第１の液化ガスとはメタン価が異なる第２の液化ガスからなる。図２中の実線部において、ガス燃料は液体状であり、図２中の点線部において、ガス燃料は気体状である。

第１のガス燃料供給ライン５３は、第１のポンプ５３１により第１のガス燃料貯留装置５１から抜き出した第１の液化ガスを、第１の気化装置５３２で気化させた後に第１の弁５５１に導くための第１の液化ガス供給ライン５３３と、第１のガス燃料貯留装置５１内で第１の液化ガスが気化した第１のボイルオフガスを第１の弁５５１に導くための第１のボイルオフガス供給ライン５３４と、を含む。

第２のガス燃料供給ライン５４は、第２のポンプ５４１により第２のガス燃料貯留装置５２から抜き出した第２の液化ガスを、第２の気化装置５４２で気化させた後に第２の弁５５２に導くための第２の液化ガス供給ライン５４３と、第２のガス燃料貯留装置５２内で第２の液化ガスが気化した第２のボイルオフガスを第２の弁５５２に導くための第２のボイルオフガス供給ライン５４４と、を含む。

（第１のメタン価の推定方法）
図３は、本開示の第１の実施形態にかかるメタン価の推定方法のフロー図である。
幾つかの実施形態にかかるメタン価の推定方法１００は、上述した対象メタン価ＭＮｔを推定するための方法である。メタン価の推定方法１００は、図３に示されるように、第１の関連付け情報取得ステップＳ１０１と、点火時期取得ステップＳ１０２と、ノッキング強度取得ステップＳ１０３と、メタン価推定ステップＳ１０４と、を備える。

図示される実施形態では、メタン価の推定方法１００における幾つかのステップ（第１の関連付け情報取得ステップＳ１０１、点火時期取得ステップＳ１０２、ノッキング強度取得ステップＳ１０３、メタン価推定ステップＳ１０４）は、メタン価の推定装置１（１Ａ）により行われる。メタン価の推定装置１（１Ａ）は、第１の関連付け情報取得ステップＳ１０１、点火時期取得ステップＳ１０２、ノッキング強度取得ステップＳ１０３およびメタン価推定ステップＳ１０４を実行可能に構成されており、これらのステップを行うようになっている。なお、メタン価の推定方法１００における幾つかのステップは、メタン価の推定装置１（１Ａ）以外の装置や機器により行われてもよいし、手動により行うようにしてもよい。

第１の関連付け情報取得ステップＳ１０１では、ガスエンジン２における、点火時期ＩＴとノッキング強度ＫＩと供給されるガス燃料のメタン価ＭＮとが予め関連付けられた第１の関連付け情報ＡＩ１、を取得することが行われる。図示される実施形態では、第１の関連付け情報取得部１１が、第１の関連付け情報取得ステップＳ１０１を実行する。第１の関連付け情報ＡＩ１は、第１の関連付け情報取得ステップＳ１０１よりも前に予め作成され、データベース部４６に記憶されている。第１の関連付け情報取得部１１は、データベース部４６から第１の関連付け情報ＡＩ１を取得する。

データベース部４６は、図３に示される実施形態では、メタン価の推定装置１の外部に設けられているが、他の幾つかの実施形態では、メタン価の推定装置１の内部に設けられていてもよい。また、データベース部４６は、制御装置３の内部に設けられていてもよいし、ガスエンジン２や制御装置３から遠隔地に設けられて、ネットワーク回線などを介して、メタン価の推定装置１と通信可能に設けられていてもよい。

第１の関連付け情報ＡＩ１は、ガスエンジン２における点火時期ＩＴと、ガスエンジン２におけるノッキング強度ＫＩと、ガスエンジン２に供給されるガス燃料のメタン価ＭＮと、の間の対応関係を示すものであり、点火時期ＩＴおよびノッキング強度ＫＩを入力情報とした際に、入力情報である点火時期ＩＴおよびノッキング強度ＫＩに対応するメタン価ＭＮを出力情報として取得できるものであればよい。

点火時期取得ステップＳ１０２では、対象ガスエンジン２ｔにおける点火時期ＩＴである対象点火時期ＩＴｔを取得することが行われる。図示される実施形態では、点火時期取得部１２が、点火時期取得ステップＳ１０２を実行する。点火時期取得部１２は、燃焼制御部３２や点火時期制御部３５などから対象点火時期ＩＴｔを取得する。

ノッキング強度取得ステップＳ１０３では、対象ガスエンジン２ｔの対象点火時期ＩＴｔを含むサイクルにおけるノッキング強度ＫＩである対象ノッキング強度ＫＩｔを取得することが行われる。図示される実施形態では、ノッキング強度取得部１３が、ノッキング強度取得ステップＳ１０３を実行する。

対象ノッキング強度ＫＩｔは、エンジンシステム４において一般的に測定されるパラメータ（例えば、ノッキングセンサ４４の検出値および筒内圧センサ４５の検出値）を用いて公知の手法により算出される。或る実施形態では、ノッキング強度取得部１３は、ノッキングセンサ４４の検出値および筒内圧センサ４５の検出値を取得し、ノッキングセンサ４４の検出値から横軸をクランク角とし、縦軸を圧力とする波形を形成する。ノッキング強度取得部１３は、上記波形から筒内圧センサ４５の検出値を用いて低周波成分を除去した後の、各サイクルにおける最大圧力振幅を対象ノッキング強度ＫＩｔとしている。なお、燃焼診断部３３が対象ノッキング強度ＫＩｔを算出している場合には、ノッキング強度取得部１３は、燃焼診断部３３が算出した対象ノッキング強度ＫＩｔを取得してもよい。

メタン価推定ステップＳ１０４では、上述した第１の関連付け情報ＡＩ１に基づいて、対象点火時期ＩＴｔおよび対象ノッキング強度ＫＩｔから、対象メタン価ＭＮｔを推定することが行われる。図示される実施形態では、メタン価推定部１４が、メタン価推定ステップＳ１０４を実行する。具体的には、メタン価推定部１４が、第１の関連付け情報取得部１１が取得した第１の関連付け情報ＡＩ１に基づいて、点火時期取得部１２が取得した対象点火時期ＩＴｔおよびノッキング強度取得部１３が取得した対象ノッキング強度ＫＩｔから、対象メタン価ＭＮｔを推定する。

上記の方法によれば、ガスエンジン２における、点火時期ＩＴとノッキング強度ＫＩと供給されるガス燃料のメタン価ＭＮとの関係性を示す第１の関連付け情報ＡＩ１を利用することで、対象点火時期ＩＴｔおよび対象ノッキング強度ＫＩｔから、対象メタン価ＭＮｔを推定できる。点火時期ＩＴやノッキング強度ＫＩは、対象ガスエンジン２ｔの運転時において対象ガスエンジン２ｔに通常装着されるセンサなどから取得可能である。よって、上記の方法によれば、対象ガスエンジン２ｔの運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価ＭＮｔを推定可能である。

幾つかの実施形態では、上述したメタン価の推定装置１（１Ａ）は、図３に示されるように、上述した第１の関連付け情報ＡＩ１を取得する第１の関連付け情報取得部１１と、上述した対象点火時期ＩＴｔを取得する点火時期取得部１２と、上述した対象点火時期ＩＴｔを含むサイクルにおける対象ノッキング強度ＫＩｔを取得するノッキング強度取得部１３と、上述した第１の関連付け情報ＡＩ１に基づいて、対象点火時期ＩＴｔおよび対象ノッキング強度ＫＩｔから、対象メタン価ＭＮｔを推定するメタン価推定部１４と、を備える。

上記の構成によれば、ガスエンジン２における、点火時期ＩＴとノッキング強度ＫＩと供給されるガス燃料のメタン価ＭＮとの関係性を示す第１の関連付け情報ＡＩ１を利用することで、対象点火時期ＩＴｔおよび対象ノッキング強度ＫＩｔから、対象メタン価ＭＮｔを推定できる。点火時期ＩＴやノッキング強度ＫＩは、対象ガスエンジン２ｔの運転時において対象ガスエンジン２ｔに通常装着されるセンサなどから取得可能である。よって、上記の構成によれば、対象ガスエンジン２ｔの運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価ＭＮｔを推定可能である。

図４は、一気筒におけるノッキング強度を点火時期との関係を説明するための説明図である。図５は、特定のノッキング強度を用いたメタン価の推定方法を説明するための説明図である。
幾つかの実施形態では、上述したメタン価推定部１４（メタン価推定ステップＳ１０４）は、図５に示されるように、対象ガスエンジン２ｔの対象エンジン負荷ＥＬｔが一定である一定負荷条件において、特定のノッキング強度ＫＩｓに到達したときの対象ノッキング強度ＫＩｔ１、および特定のノッキング強度ＫＩｓに到達したときの対象ノッキング強度ＫＩｔ１を含むサイクルにおける対象点火時期ＩＴｔ１から、第１の関連付け情報ＡＩ１に基づいて、対象ガスエンジン２ｔにおける気筒２１に供給されたガス燃料のメタン価（対象メタン価ＭＮｔ）を推定する。

図４は、対象ガスエンジン２ｔの対象エンジン負荷ＥＬｔが一定である一定負荷条件における、一つの気筒２１の時間Ｔの経過に対する、ノッキング強度ＫＩと点火時期ＩＴの変化が示されている。図４に示されるように、点火時期ＩＴを徐々に進角させると、ノッキング強度ＫＩ（振幅幅）が徐々に増大し、特定のノッキング強度ＫＩｓを超える。

図５では、対象ガスエンジン２ｔの対象エンジン負荷ＥＬｔが一定である一定負荷条件において、対象ガスエンジン２ｔの一つの気筒２１の点火時期ＩＴを横軸とし、上記一つの気筒２１のノッキング強度ＫＩを縦軸とするグラフが示されている。上記グラフには、第１のメタン価を有するガス燃料における点火時期ＩＴとノッキング強度ＫＩとの関係を示す第１の曲線Ｃ１と、第１のメタン価よりも高い第２のメタン価を有するガス燃料における点火時期ＩＴとノッキング強度ＫＩとの関係を示す第２の曲線Ｃ２と、が示されている。第１の曲線Ｃ１や第２の曲線Ｃ２は、上述した第１の関連付け情報ＡＩ１に含まれる。

上述したメタン価推定部１４（メタン価推定ステップＳ１０４）は、図５に示されるように、対象ガスエンジン２ｔの一つの気筒２１において、対象ノッキング強度ＫＩｔが特定のノッキング強度ＫＩｓに到達するまで、予め設定された設定点火時期ＩＴｓから対象点火時期ＩＴｔを進角させるように、点火時期制御部３５に対して指示する。

上述したメタン価推定部１４（メタン価推定ステップＳ１０４）は、図５に示されるように、対象ノッキング強度ＫＩｔが特定のノッキング強度ＫＩｓに到達した（特定のノッキング強度ＫＩｓを上回った）サイクルにおける対象ノッキング強度ＫＩｔ１および対象点火時期ＩＴｔ１を求める。上述したメタン価推定部１４（メタン価推定ステップＳ１０４）は、第１の曲線Ｃ１の特定のノッキング強度ＫＩｓにおける点火時期ＩＴｃ１および第２の曲線Ｃ２の特定のノッキング強度ＫＩｓにおける点火時期ＩＴｃ２から、対象点火時期ＩＴｔ１の割合を求め、対象点火時期ＩＴｔ１の割合と第１のメタン価と第２のメタン価から対象メタン価ＭＮｔを算出できる。

特定のノッキング強度ＫＩｓを高く設定するほど、ガス燃料毎の対象点火時期ＩＴｔの差を大きくできるため、対象メタン価ＭＮｔの推定精度の向上が図れるが、ノッキングが生じる可能性が高くなる虞がある。上記の構成によれば、特定のノッキング強度ＫＩｓに到達したサイクルにおける対象ノッキング強度ＫＩｔ１や対象点火時期ＩＴｔ１を用いて対象メタン価ＭＮｔを推定することで、対象ガスエンジン２ｔのノッキングを抑制しつつ、対象メタン価ＭＮｔの推定精度を高いものにすることができる。

メタン価推定部１４が推定した対象メタン価ＭＮｔは、燃焼制御部３２に送られる。燃焼制御部３２は、メタン価推定部１４が推定した対象メタン価ＭＮｔに基づいて、対象ガスエンジン２ｔにおける各気筒２１の点火時期を決定する。点火時期制御部３５は、メタン価推定部１４が推定した対象メタン価ＭＮｔに基づいて、燃焼制御部３２において決定された点火時期に点火装置４２が点火するように、点火装置４２の点火時期を制御するように構成されている。つまり、点火時期制御部３５は、メタン価推定部１４が推定した対象メタン価ＭＮｔに応じて、対象点火時期ＩＴｔを制御するように構成されている。

幾つかの実施形態では、メタン価の推定装置１又は制御装置３の少なくとも一方は、対象ガスエンジン２ｔの一つの気筒２１においてメタン価推定部１４が推定した対象メタン価ＭＮｔを、対象ガスエンジン２ｔの上記一つの気筒２１以外の他の気筒２１に反映させてもよい。この場合には、対象メタン価ＭＮｔを推定する際の対象ガスエンジン２ｔの性能変化の影響を小さくできる。

幾つかの実施形態では、メタン価の推定装置１は、対象ガスエンジン２ｔの気筒２１毎に対象メタン価ＭＮｔを推定してもよい。この場合には、気筒２１毎に適切な対象メタン価ＭＮｔを推定できるため、気筒２１毎の燃焼のばらつきを抑制できる。

（第２のメタン価の推定方法）
図６は、本開示の第２の実施形態にかかるメタン価の推定方法のフロー図である。
幾つかの実施形態にかかるメタン価の推定方法２００は、上述した対象メタン価ＭＮｔを推定するための方法である。メタン価の推定方法２００は、図６に示されるように、第２の関連付け情報取得ステップＳ２０１と、エンジン負荷取得ステップＳ２０２と、ガス供給量取得ステップＳ２０３と、メタン価推定ステップＳ２０４と、を備える。

図示される実施形態では、メタン価の推定方法２００における幾つかのステップ（第２の関連付け情報取得ステップＳ２０１、エンジン負荷取得ステップＳ２０２、ガス供給量取得ステップＳ２０３、メタン価推定ステップＳ２０４）は、メタン価の推定装置１（１Ｂ）により行われる。メタン価の推定装置１（１Ｂ）は、第２の関連付け情報取得ステップＳ２０１、エンジン負荷取得ステップＳ２０２、ガス供給量取得ステップＳ２０３およびメタン価推定ステップＳ２０４を実行可能に構成されており、これらのステップを行うようになっている。なお、メタン価の推定方法２００における幾つかのステップは、メタン価の推定装置１（１Ｂ）以外の装置や機器により行われてもよいし、手動により行うようにしてもよい。

第２の関連付け情報取得ステップＳ２０１では、ガスエンジン２における、エンジン負荷ＥＬと供給されるガス燃料のガス供給量ＦＳとガス燃料のメタン価ＭＮとが予め関連付けられた第２の関連付け情報ＡＩ２、を取得することが行われる。図示される実施形態では、第２の関連付け情報取得部６１が、第２の関連付け情報取得ステップＳ２０１を実行する。第２の関連付け情報ＡＩ２は、第２の関連付け情報取得ステップＳ２０１よりも前に予め作成され、データベース部４６に記憶されている。第２の関連付け情報取得部６１は、データベース部４６から第２の関連付け情報ＡＩ２を取得する。

第２の関連付け情報ＡＩ２は、ガスエンジン２におけるエンジン負荷ＥＬと、ガスエンジン２に供給されるガス燃料のガス供給量ＦＳと、ガスエンジン２に供給されるガス燃料のメタン価ＭＮと、の間の対応関係を示すものであり、エンジン負荷ＥＬおよびガス供給量ＦＳを入力情報とした際に、入力情報であるエンジン負荷ＥＬおよびガス供給量ＦＳに対応するメタン価ＭＮを出力情報として取得できるものであればよい。

図７は、ガス燃料の産地毎のメタン価と発熱量との関係を示す図である。図７では、発熱量を横軸とし、メタン価を縦軸とするグラフが示されており、このグラフ上に産地毎にガス燃料がプロットされている。図７に示されるように、発熱量の増加に伴い、メタン価が低減する傾向がある。

図８は、ガス燃料毎のガスエンジンにおけるエンジン負荷とガスエンジンに供給されるガス燃料のガス供給量との関係を示す図である。図８では、ガスエンジンにおけるエンジン負荷を横軸とし、ガスエンジンに供給されるガス燃料のガス供給量を縦軸とするグラフが示されている。上記グラフには、第１のメタン価を有するガス燃料における上記エンジン負荷と上記ガス供給量との関係を示す第１の直線ＳＬ１と、第１のメタン価よりも高い第２のメタン価を有するガス燃料における上記エンジン負荷と上記ガス供給量との関係を示す第２の直線ＳＬ２と、が示されている。図８に示されるように、エンジン負荷の増加に伴い、ガス供給量が増加する傾向がある。

一般的に、発熱量の増加に伴い、ガス供給量が低減する傾向がある。以上から、ガス供給量と発熱量とメタン価との間には、一定の相関がある。発熱量の増加に伴い、ガス供給量とメタン価とが低減する傾向がある。第２の関連付け情報ＡＩ２には、ガス供給量と発熱量とメタン価との間の相関を示す情報が含まれる。

エンジン負荷取得ステップＳ２０２では、対象ガスエンジン２ｔにおけるエンジン負荷ＥＬである対象エンジン負荷ＥＬｔを取得することが行われる。図示される実施形態では、エンジン負荷取得部６２が、エンジン負荷取得ステップＳ２０２を実行する。エンジン負荷取得部６２は、対象ガスエンジン２ｔや対象ガスエンジン２ｔから動力を伝達可能に接続された発電機４３、対象ガスエンジン２ｔに通常装着されるセンサ、エンジン制御部３１、燃焼制御部３２などから対象エンジン負荷ＥＬｔを取得する。なお、上記エンジン負荷ＥＬや対象エンジン負荷ＥＬｔは、ガスエンジン２に動力を伝達可能に接続された発電機の負荷であってもよい。エンジン負荷取得ステップＳ２０２では、上述した発電機４３の負荷を対象エンジン負荷ＥＬｔとして取得してもよい。

ガス供給量取得ステップＳ２０３では、対象ガスエンジン２ｔの対象エンジン負荷ＥＬｔを取得した期間におけるガス供給量ＦＳである対象ガス供給量ＦＳｔを取得することが行われる。図示される実施形態では、ガス供給量取得部６３が、ガス供給量取得ステップＳ２０３を実行する。ガス供給量取得部６３は、対象ガスエンジン２ｔや対象ガスエンジン２ｔに通常装着されるセンサ、燃焼制御部３２、ガス燃料供給量制御部３４などから対象ガス供給量ＦＳｔを取得する。

メタン価推定ステップＳ２０４では、上述した第２の関連付け情報ＡＩ２に基づいて、対象エンジン負荷ＥＬｔおよび対象ガス供給量ＦＳｔから、対象メタン価ＭＮｔを推定することが行われる。図示される実施形態では、メタン価推定部６４が、メタン価推定ステップＳ２０４を実行する。具体的には、メタン価推定部６４が、第２の関連付け情報取得部６１が取得した第２の関連付け情報ＡＩ２に基づいて、エンジン負荷取得部６２が取得した対象エンジン負荷ＥＬｔおよびガス供給量取得部６３が取得した対象ガス供給量ＦＳｔから、対象メタン価ＭＮｔを推定する。

上記の方法によれば、ガスエンジン２における、エンジン負荷ＥＬと供給されるガス燃料のガス供給量ＦＳとメタン価ＭＮとの関係性を示す第２の関連付け情報ＡＩ２を利用することで、対象エンジン負荷ＥＬｔおよび対象ガス供給量ＦＳｔから、対象メタン価ＭＮｔを推定できる。エンジン負荷ＥＬやガス供給量ＦＳは、対象ガスエンジン２ｔの運転時において対象ガスエンジン２ｔに通常装着されるセンサなどから取得可能である。よって、上記の方法によれば、対象ガスエンジン２ｔの運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価ＭＮｔを推定可能である。

幾つかの実施形態では、上述したメタン価の推定装置１（１Ｂ）は、図６に示されるように、上述した第２の関連付け情報ＡＩ２を取得する第２の関連付け情報取得部６１と、上述した対象エンジン負荷ＥＬｔを取得するエンジン負荷取得部６２と、上述した対象エンジン負荷ＥＬｔを取得した期間における対象ガス供給量ＦＳｔを取得するガス供給量取得部６３と、上述した第２の関連付け情報ＡＩ２に基づいて、対象エンジン負荷ＥＬｔおよび対象ガス供給量ＦＳｔから、対象メタン価ＭＮｔを推定するメタン価推定部６４と、を備える。

上記の構成によれば、ガスエンジン２における、エンジン負荷ＥＬと供給されるガス燃料のガス供給量ＦＳとメタン価ＭＮとの関係性を示す第２の関連付け情報ＡＩ２を利用することで、対象エンジン負荷ＥＬｔおよび対象ガス供給量ＦＳｔから、対象メタン価ＭＮｔを推定できる。エンジン負荷ＥＬやガス供給量ＦＳは、対象ガスエンジン２ｔの運転時において対象ガスエンジン２ｔに通常装着されるセンサなどから取得可能である。よって、上記の構成によれば、対象ガスエンジン２ｔの運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価ＭＮｔを推定可能である。

（第３のメタン価の推定方法）
図９は、本開示の第３の実施形態にかかるメタン価の推定方法のフロー図である。
幾つかの実施形態にかかるメタン価の推定方法３００は、上述した対象メタン価ＭＮｔを推定するための方法である。メタン価の推定方法３００は、図９に示されるように、第３の関連付け情報取得ステップＳ３０１と、熱発生量取得ステップＳ３０２と、供給量取得ステップＳ３０３と、低位発熱量算出ステップＳ３０４と、メタン価推定ステップＳ３０５と、を備える。

図示される実施形態では、メタン価の推定方法３００における幾つかのステップ（第３の関連付け情報取得ステップＳ３０１、熱発生量取得ステップＳ３０２、供給量取得ステップＳ３０３、低位発熱量算出ステップＳ３０４、メタン価推定ステップＳ３０５）は、メタン価の推定装置１（１Ｃ）により行われる。メタン価の推定装置１（１Ｃ）は、第３の関連付け情報取得ステップＳ３０１、熱発生量取得ステップＳ３０２、供給量取得ステップＳ３０３、低位発熱量算出ステップＳ３０４およびメタン価推定ステップＳ３０５を実行可能に構成されており、これらのステップを行うようになっている。なお、メタン価の推定方法３００における幾つかのステップは、メタン価の推定装置１（１Ｃ）以外の装置や機器により行われてもよいし、手動により行うようにしてもよい。

第３の関連付け情報取得ステップＳ３０１では、低位発熱量ＬＨＶと供給されるガス燃料のメタン価ＭＮとが予め関連付けられた第３の関連付け情報ＡＩ３、を取得することが行われる。図示される実施形態では、第３の関連付け情報取得部７１が、第３の関連付け情報取得ステップＳ３０１を実行する。第３の関連付け情報ＡＩ３は、第３の関連付け情報取得ステップＳ３０１よりも前に予め作成され、データベース部４６に記憶されている。第３の関連付け情報取得部７１は、データベース部４６から第３の関連付け情報ＡＩ３を取得する。

第３の関連付け情報ＡＩ３は、ガスエンジン２における低位発熱量ＬＨＶと、ガスエンジン２に供給されるガス燃料のメタン価ＭＮと、の間の対応関係を示すものであり、低位発熱量ＬＨＶを入力情報とした際に、入力情報である低位発熱量ＬＨＶに対応するメタン価ＭＮを出力情報として取得できるものであればよい。

熱発生量取得ステップＳ３０２では、対象ガスエンジン２ｔにおける１サイクルあたりの熱発生量ＱＣである対象熱発生量ＱＣｔを取得することが行われる。図示される実施形態では、熱発生量取得部７２が、熱発生量取得ステップＳ３０２を実行する。

図１０は、クランク角に対する筒内圧の変化を示す図である。図１０では、クランク角θを横軸とし、筒内圧を縦軸とするグラフが示されている。このグラフには、対象ガスエンジン２ｔにおいてガス燃料を燃焼させた場合における複数の筒内圧センサ４５の検出値の平均値から形成された波形ＰＣＡと、対象ガスエンジン２ｔにおいてガス燃料を燃焼させない場合における筒内圧センサ４５の検出値から形成された波形ＰＣＭとが示されている。

熱発生量取得部７２は、筒内圧センサ４５の検出値から形成された波形から１サイクルあたりの対象熱発生量ＱＣｔを導出する。１サイクルあたりの対象熱発生量ＱＣｔは、１サイクルあたりの筒内圧の変化や気筒２１内容積の変化から、公知の手法により導出可能である。或る実施形態では、１サイクルあたりの筒内圧の変化や気筒２１内容積の変化から、クランク角θ毎の熱発生速度を導出し、クランク角θ毎の熱発生速度を積分することで、１サイクルあたりの対象熱発生量ＱＣｔを算出する。

供給量取得ステップＳ３０３では、対象熱発生量ＱＣｔを取得した期間における対象ガスエンジン２ｔに供給されるガス燃料の１サイクルあたりの供給量ＭＦである対象供給量ＭＦｔを取得することが行われる。図示される実施形態では、供給量取得部７３が、供給量取得ステップＳ３０３を実行する。供給量取得部７３は、対象ガスエンジン２ｔや対象ガスエンジン２ｔに通常装着されるセンサ、燃焼制御部３２、ガス燃料供給量制御部３４などから対象供給量ＭＦｔを取得する。

低位発熱量算出ステップＳ３０４では、対象熱発生量ＱＣｔと対象供給量ＭＦｔから、対象ガスエンジン２ｔにおける低位発熱量ＬＨＶである対象低位発熱量ＬＨＶｔを算出することが行われる。図示される実施形態では、低位発熱量算出部７４が、低位発熱量算出ステップＳ３０４を実行する。

図示される実施形態では、低位発熱量算出部７４は、以下の式（１）を用いて、熱発生量取得部７２で取得した対象熱発生量ＱＣｔと、供給量取得部７３で取得した対象供給量ＭＦｔから、対象低位発熱量ＬＨＶｔを算出する。
ＬＨＶｔ ＝ （ＱＣｔ＋Ｑｈｌ）／ＭＦｔ ・・・（１）
なお、上記式（１）におけるＱｈｌは、熱損失であり、定数としてもよい。

なお、低位発熱量算出部７４は、上記式（１）で算出した対象低位発熱量ＬＨＶｔにカルマンフィルターなどのフィルタを用いて推定した対象低位発熱量ＬＨＶｔを、低位発熱量算出部７４で取得する対象低位発熱量ＬＨＶｔとしてもよい。

図１１は、低位発熱量とガス燃料のメタン価との関係を示す図である。図１１では、低位発熱量ＬＨＶを横軸とし、ガス燃料のメタン価ＭＮを縦軸とするグラフである。このグラフには、産地ごとのガス燃料がプロットされており、プロットから形成された回帰直線ＲＬが示されている。図１１に示されるように、低位発熱量ＬＨＶの増大に伴い、メタン価が低減する傾向がある。

メタン価推定ステップＳ３０５では、上述した第３の関連付け情報ＡＩ３に基づいて、対象低位発熱量ＬＨＶｔから、対象メタン価ＭＮｔを推定することが行われる。図示される実施形態では、メタン価推定部７５が、メタン価推定ステップＳ３０５を実行する。具体的には、メタン価推定部７５が、第３の関連付け情報取得部７１が取得した第３の関連付け情報ＡＩ３に基づいて、低位発熱量算出部７４が算出した対象低位発熱量ＬＨＶｔから、対象メタン価ＭＮｔを推定する。

上記の方法によれば、ガスエンジン２における、低位発熱量ＬＨＶと供給されるガス燃料のメタン価ＭＮとの関係性を示す第３の関連付け情報ＡＩ３を利用することで、対象低位発熱量ＬＨＶｔから、対象メタン価ＭＮｔを推定できる。低位発熱量ＬＨＶは、対象ガスエンジン２ｔの運転時において対象ガスエンジン２ｔに通常装着されるセンサなどから取得できる、上記熱発生量ＱＣおよび上記供給量ＭＦから算出できる。よって、上記の方法によれば、対象ガスエンジン２ｔの運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価ＭＮｔを推定可能である。

幾つかの実施形態では、上述したメタン価の推定装置１（１Ｃ）は、図９に示されるように、上述した第３の関連付け情報ＡＩ３を取得する第３の関連付け情報取得部７１と、上述した対象熱発生量ＱＣｔを取得する熱発生量取得部７２と、上述した対象熱発生量ＱＣｔを取得した期間における対象供給量ＭＦｔを取得する供給量取得部７３と、上述した対象熱発生量ＱＣｔと上述した対象供給量ＭＦｔから、対象低位発熱量ＬＨＶｔを算出する低位発熱量算出部７４と、上述した第３の関連付け情報ＡＩ３に基づいて、対象低位発熱量ＬＨＶｔから、対象メタン価ＭＮｔを推定するメタン価推定部７５と、を備える。

上記の構成によれば、ガスエンジン２における、低位発熱量ＬＨＶと供給されるガス燃料のメタン価ＭＮとの関係性を示す第３の関連付け情報ＡＩ３を利用することで、対象低位発熱量ＬＨＶｔから、対象メタン価ＭＮｔを推定できる。低位発熱量ＬＨＶは、対象ガスエンジン２ｔの運転時において対象ガスエンジン２ｔに通常装着されるセンサなどから取得できる、熱発生量ＱＣおよび供給量ＭＦから算出できる。よって、上記の構成によれば、対象ガスエンジン２ｔの運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価ＭＮｔを推定可能である。

上述した第１の関連付け情報ＡＩ１、第２の関連付け情報ＡＩ２および第３の関連付け情報ＡＩ３の夫々には、上記入力情報と上記出力情報との対応関係を示すリストや表、マップ、関数、機械学習のモデルなどが含まれる。第１の関連付け情報ＡＩ１、第２の関連付け情報ＡＩ２および第３の関連付け情報ＡＩ３の夫々は、定常試験データを基に作成してもよいし、定常試験データ以外の過去の実績値や実験値、数値解析結果などを基に作成してもよい。

上述した第１の関連付け情報ＡＩ１、第２の関連付け情報ＡＩ２および第３の関連付け情報ＡＩ３の夫々は、対象ガスエンジン２ｔから取得される情報だけでなく、対象ガスエンジン２ｔ以外のガスエンジン２から取得される情報を含んでいてもよい。また、第１の関連付け情報ＡＩ１、第２の関連付け情報ＡＩ２および第３の関連付け情報ＡＩ３の夫々は、対象ガスエンジン２ｔから取得される情報を含まずに対象ガスエンジン２ｔ以外のガスエンジン２から取得される情報のみを含んでいてもよい。これらの場合には、情報を取得する対象ガスエンジン２ｔ以外のガスエンジン２は、対象ガスエンジン２ｔと同一機種や類似機種であることが望ましい。

幾つかの実施形態にかかるガスエンジン２の制御装置３は、図１に示されるように、上述したメタン価の推定装置１（１Ａ、１Ｂ、１Ｃ）と、メタン価の推定装置１で推定されたガス燃料の対象メタン価ＭＮｔに応じて、対象ガスエンジン２ｔの点火時期ＩＴである対象点火時期ＩＴｔを制御するように構成された点火時期制御部３５と、を備える。

図１２は、ガスエンジンに供給されるガス燃料のメタン価とガスエンジンの点火時期との関係を説明するための説明図である。図１２に示されるように、メタン価ＭＮが高ければ、その分だけ点火時期ＩＴを進角させることができる。点火時期制御部３５は、対象メタン価ＭＮｔが高くなった場合には、その分だけ対象点火時期ＩＴｔを進角させ、対象メタン価ＭＮｔが低くなった場合には、その分だけ対象点火時期ＩＴｔを遅角させることで、対象ガスエンジン２ｔの高効率での運転が可能となる。

上記の構成によれば、ガスエンジン２の制御装置３は、点火時期制御部３５において、メタン価の推定装置１で推定されたガス燃料のメタン価ＭＮｔに応じて、対象ガスエンジン２ｔの点火時期ＩＴである対象点火時期ＩＴｔを制御することで、対象ガスエンジン２ｔにおける気筒２１に供給されるガス燃料のメタン価（対象メタン価ＭＮｔ）の変化に応じた対象ガスエンジン２ｔの燃焼制御が可能となる。この場合には、ガス燃料の切り替え時に連続してガス燃料のメタン価ＭＮｔを推定し、連続して対象点火時期ＩＴｔを調整できるため、ガス燃料の切り替え時に対象ガスエンジン２ｔを停止させずに高効率での連続運転が可能となる。

メタン価の推定装置１は、供給されるガス燃料が、第１の液化ガスを気化した気化ガスと第１のボイルオフガスの混合ガス、第２の液化ガスを気化した気化ガスと第２のボイルオフガスの混合ガス、又は第１の液化ガスを気化した気化ガスと第２の液化ガスを気化した気化ガスの混合ガスの何れにおいても対象メタン価を推定可能である。このため、上記制御装置３を備える対象ガスエンジン２ｔは、上述した混合ガスがガス燃料として供給された場合においても、高効率での連続運転が可能である。

本開示は上述した実施形態に限定されることはなく、上述した実施形態に変形を加えた形態や、これらの形態を適宜組み合わせた形態も含む。

上述した幾つかの実施形態に記載の内容は、例えば以下のように把握されるものである。

１）本開示の少なくとも一実施形態にかかるメタン価の推定装置（１Ａ（１））は、
少なくとも１つの気筒（２１）を備えるガスエンジン（２）である対象ガスエンジン（２ｔ）に供給されるガス燃料のメタン価（ＭＮｔ）を推定するメタン価の推定装置（１Ａ）であって、
ガスエンジン（２）における、点火時期（ＩＴ）とノッキング強度（ＫＩ）と供給されるガス燃料のメタン価（ＭＮ）とが予め関連付けられた第１の関連付け情報（ＡＩ１）、を取得する第１の関連付け情報取得部（１１）と、
前記対象ガスエンジン（２ｔ）における前記点火時期（ＩＴ）である対象点火時期（ＩＴｔ）を取得する点火時期取得部（１２）と、
前記対象ガスエンジン（２ｔ）の前記対象点火時期（ＩＴｔ）を含むサイクルにおける前記ノッキング強度（ＫＩ）である対象ノッキング強度（ＫＩｔ）を取得するノッキング強度取得部（１３）と、
前記第１の関連付け情報（ＡＩ１）に基づいて、前記対象点火時期（ＩＴｔ）および前記対象ノッキング強度（ＫＩｔ）から、前記対象ガスエンジン（２ｔ）における前記気筒（２１）に供給された前記ガス燃料のメタン価（ＭＮｔ）を推定するメタン価推定部（１４）と、を備える。

上記１）の構成によれば、ガスエンジン（２）における、点火時期（ＩＴ）とノッキング強度（ＫＩ）と供給されるガス燃料のメタン価（ＭＮ）との関係性を示す第１の関連付け情報（ＡＩ１）を利用することで、対象点火時期（ＩＴｔ）および対象ノッキング強度（ＫＩｔ）から、対象メタン価ＭＮｔ（対象ガスエンジン２ｔにおける気筒２１に供給されたガス燃料のメタン価）を推定できる。点火時期（ＩＴ）やノッキング強度（ＫＩ）は、対象ガスエンジン（２ｔ）の運転時において対象ガスエンジン（２ｔ）に通常装着されるセンサなどから取得可能である。よって、上記１）の構成によれば、対象ガスエンジン（２ｔ）の運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価（ＭＮｔ）を推定可能である。

２）幾つかの実施形態では、上記１）に記載のメタン価の推定装置（１Ａ）であって、
前記メタン価推定部（１４）は、
前記対象ガスエンジン（２ｔ）のエンジン負荷（ＥＬｔ）が一定である一定負荷条件において、特定のノッキング強度に到達したときの前記対象ノッキング強度（ＫＩｔ１）、および前記特定のノッキング強度に到達したときの前記対象ノッキング強度（ＫＩｔ１）を含むサイクルにおける前記対象点火時期（ＩＴｔ１）から、前記第１の関連付け情報（ＡＩ１）に基づいて、前記対象ガスエンジン（２ｔ）における前記気筒（２１）に供給された前記ガス燃料のメタン価（ＭＮｔ）を推定する。

特定のノッキング強度を高く設定するほど、ガス燃料毎の対象点火時期（ＩＴｔ１）の差を大きくできるため、対象メタン価（ＭＮｔ）の推定精度の向上が図れるが、ノッキングが生じる可能性が高くなる虞がある。上記２）の構成によれば、特定のノッキング強度に到達したサイクルにおける対象ノッキング強度（ＫＩｔ１）や対象点火時期（ＩＴｔ１）を用いて対象メタン価（ＭＮｔ）を推定することで、対象ガスエンジン（２ｔ）のノッキングを抑制しつつ、対象メタン価（ＭＮｔ）の推定精度を高いものにすることができる。

３）本開示の少なくとも一実施形態にかかるメタン価の推定装置（１Ｂ（１））は、
少なくとも１つの気筒（２１）を備えるガスエンジン（２）である対象ガスエンジン（２ｔ）に供給されるガス燃料のメタン価（ＭＮｔ）を推定するメタン価の推定装置（１Ｂ）であって、
ガスエンジン（２）における、エンジン負荷（ＥＬ）と供給されるガス燃料のガス供給量（ＦＳ）と前記ガス燃料のメタン価（ＭＮ）とが予め関連付けられた第２の関連付け情報（ＡＩ２）、を取得する第２の関連付け情報取得部（６１）と、
前記対象ガスエンジン（２ｔ）における前記エンジン負荷（ＥＬ）である対象エンジン負荷（ＥＬｔ）を取得するエンジン負荷取得部（６２）と、
前記対象ガスエンジン（２ｔ）の前記対象エンジン負荷（ＥＬｔ）を取得した期間における前記ガス供給量（ＦＳ）である対象ガス供給量（ＦＳｔ）を取得するガス供給量取得部（６３）と、
前記第２の関連付け情報（ＡＩ２）に基づいて、前記対象エンジン負荷（ＥＬｔ）および前記対象ガス供給量（ＦＳｔ）から、前記対象ガスエンジン（２ｔ）における前記気筒（２１）に供給された前記ガス燃料のメタン価（ＭＮｔ）を推定するメタン価推定部（６４）と、を備える。

上記３）の構成によれば、ガスエンジン（２）における、エンジン負荷（ＥＬ）と供給されるガス燃料のガス供給量（ＦＳ）とメタン価（ＭＮ）との関係性を示す第２の関連付け情報（ＡＩ２）を利用することで、対象エンジン負荷（ＥＬｔ）および対象ガス供給量（ＦＳｔ）から、対象メタン価（ＭＮｔ、対象ガスエンジン２ｔにおける気筒２１に供給されたガス燃料のメタン価）を推定できる。エンジン負荷（ＥＬ）やガス供給量（ＦＳ）は、対象ガスエンジン（２ｔ）の運転時において対象ガスエンジン（２ｔ）に通常装着されるセンサなどから取得可能である。よって、上記３）の構成によれば、対象ガスエンジン（２ｔ）の運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価（ＭＮｔ）を推定可能である。

４）本開示の少なくとも一実施形態にかかるメタン価の推定装置（１Ｃ（１））は、
少なくとも１つの気筒（２１）を備えるガスエンジン（２）である対象ガスエンジン（２ｔ）に供給されるガス燃料のメタン価（ＭＮｔ）を推定するメタン価の推定装置（１Ｃ）であって、
ガスエンジン（２）における、低位発熱量（ＬＨＶ）と供給されるガス燃料のメタン価（ＭＮ）とが予め関連付けられた第３の関連付け情報（ＡＩ３）、を取得する第３の関連付け情報取得部（７１）と、
前記対象ガスエンジン（２ｔ）における１サイクルあたりの熱発生量（ＱＣ）である対象熱発生量（ＱＣｔ）を取得する熱発生量取得部（７２）と、
前記対象熱発生量（ＱＣｔ）を取得した期間における前記対象ガスエンジン（２ｔ）に供給されるガス燃料の１サイクルあたりの供給量（ＭＦ）である対象供給量（ＭＦｔ）を取得する供給量取得部（７３）と、
前記対象熱発生量（ＱＣｔ）と前記対象供給量（ＭＦｔ）から、前記対象ガスエンジン（２ｔ）における低位発熱量（ＬＨＶ）である対象低位発熱量（ＬＨＶｔ）を算出する低位発熱量算出部（７４）と、
前記第３の関連付け情報（ＡＩ３）に基づいて、前記対象低位発熱量（ＬＨＶｔ）から、前記対象ガスエンジン（２ｔ）における前記気筒（２１）に供給された前記ガス燃料のメタン価（ＭＮｔ）を推定するメタン価推定部（７５）と、を備える。

上記４）の構成によれば、ガスエンジン（２）における、低位発熱量（ＬＨＶ）と供給されるガス燃料のメタン価（ＭＮ）との関係性を示す第３の関連付け情報（ＡＩ３）を利用することで、対象低位発熱量（ＬＨＶｔ）から、対象メタン価（ＭＮｔ、対象ガスエンジン２ｔにおける気筒２１に供給されたガス燃料のメタン価）を推定できる。低位発熱量（ＬＨＶ）は、対象ガスエンジン（２ｔ）の運転時において対象ガスエンジン（２ｔ）に通常装着されるセンサなどから取得できる、上記熱発生量（ＱＣ）および上記供給量（ＭＦ）から算出できる。よって、上記４）の構成によれば、対象ガスエンジン（２ｔ）の運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価（ＭＮｔ）を推定可能である。

５）本開示の少なくとも一実施形態にかかるガスエンジン（２）の制御装置（３）は、
上記１）から４）までの何れかに記載のメタン価の推定装置（１（１Ａ、１Ｂ、１Ｃ））と、
前記メタン価の推定装置（１）で推定された前記ガス燃料のメタン価（ＭＮｔ）に応じて、前記対象ガスエンジン（２ｔ）の点火時期（ＩＴ）である対象点火時期（ＩＴｔ）を制御するように構成された点火時期制御部（３５）と、を備える。

上記５）の構成によれば、ガスエンジン（２）の制御装置（３）は、点火時期制御部（３５）において、メタン価の推定装置（１）で推定されたガス燃料のメタン価（ＭＮｔ）に応じて、対象ガスエンジン（２ｔ）の点火時期（ＩＴ）である対象点火時期（ＩＴｔ）を制御することで、対象ガスエンジン（２ｔ）における気筒（２１）に供給されるガス燃料のメタン価（ＭＮｔ）の変化に応じた対象ガスエンジン（２ｔ）の燃焼制御が可能となる。この場合には、ガス燃料の切り替え時に連続してガス燃料のメタン価（ＭＮｔ）を推定し、連続して対象点火時期（ＩＴｔ）を調整できるため、ガス燃料の切り替え時に対象ガスエンジン（２ｔ）を停止させずに高効率での連続運転が可能となる。

６）本開示の少なくとも一実施形態にかかるメタン価の推定方法（１００）は、
少なくとも１つの気筒（２１）を備えるガスエンジン（２）である対象ガスエンジン（２ｔ）に供給されるガス燃料のメタン価（ＭＮｔ）を推定するメタン価の推定方法（１００）であって、
ガスエンジン（２）における、点火時期（ＩＴ）とノッキング強度（ＫＩ）と供給されるガス燃料のメタン価（ＭＮ）とが予め関連付けられた第１の関連付け情報（ＡＩ１）、を取得する第１の関連付け情報取得ステップ（Ｓ１０１）と、
前記対象ガスエンジン（２ｔ）における前記点火時期（ＩＴ）である対象点火時期（ＩＴｔ）を取得する点火時期取得ステップ（Ｓ１０２）と、
前記対象ガスエンジン（２ｔ）の前記対象点火時期（ＩＴｔ）を含むサイクルにおける前記ノッキング強度（ＫＩ）である対象ノッキング強度（ＫＩｔ）を取得するノッキング強度取得ステップ（Ｓ１０３）と、
前記第１の関連付け情報（ＡＩ１）に基づいて、前記対象点火時期（ＩＴｔ）および前記対象ノッキング強度（ＫＩｔ）から、前記対象ガスエンジン（２ｔ）における前記気筒（２１）に供給された前記ガス燃料のメタン価（ＭＮｔ）を推定するメタン価推定ステップ（Ｓ１０４）と、を備える。

上記６）の方法によれば、ガスエンジン（２）における、点火時期（ＩＴｔ）とノッキング強度（ＫＩ）と供給されるガス燃料のメタン価（ＭＮ）との関係性を示す第１の関連付け情報（ＡＩ１）を利用することで、対象点火時期（ＩＴｔ）および対象ノッキング強度（ＫＩｔ）から、対象ガスエンジン（２ｔ）における気筒（２１）に供給されたガス燃料のメタン価（対象メタン価ＭＮｔ）を推定できる。点火時期（ＩＴ）やノッキング強度（ＫＩ）は、対象ガスエンジン（２ｔ）の運転時において対象ガスエンジン（２ｔ）に通常装着されるセンサなどから取得可能である。よって、上記６）の方法によれば、対象ガスエンジン（２ｔ）の運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価（ＭＮｔ）を推定可能である。

７）本開示の少なくとも一実施形態にかかるメタン価の推定方法（２００）は、
少なくとも１つの気筒（２１）を備えるガスエンジン（２）である対象ガスエンジン（２ｔ）に供給されるガス燃料のメタン価（ＭＮｔ）を推定するメタン価の推定方法（２００）であって、
ガスエンジン（２）における、エンジン負荷（ＥＬ）と供給されるガス燃料のガス供給量（ＦＳ）と前記ガス燃料のメタン価（ＭＮ）とが予め関連付けられた第２の関連付け情報（ＡＩ２）、を取得する第２の関連付け情報取得ステップ（Ｓ２０１）と、
前記対象ガスエンジン（２ｔ）における前記エンジン負荷（ＥＬ）である対象エンジン負荷（ＥＬｔ）を取得するエンジン負荷取得ステップ（Ｓ２０２）と、
前記対象ガスエンジン（２ｔ）の前記対象エンジン負荷（ＥＬｔ）を取得した期間における前記ガス供給量（ＦＳ）である対象ガス供給量（ＦＳｔ）を取得するガス供給量取得ステップ（Ｓ２０３）と、
前記第２の関連付け情報（ＡＩ２）に基づいて、前記対象エンジン負荷（ＥＬｔ）および前記対象ガス供給量（ＦＳｔ）から、前記対象ガスエンジン（２ｔ）における前記気筒（２１）に供給された前記ガス燃料のメタン価（ＭＮｔ）を推定するメタン価推定ステップ（Ｓ２０４）と、を備える。

上記７）の方法によれば、ガスエンジン（２）における、エンジン負荷（ＥＬ）と供給されるガス燃料のガス供給量（ＦＳ）とメタン価（ＭＮ）との関係性を示す第２の関連付け情報（ＡＩ２）を利用することで、対象エンジン負荷（ＥＬｔ）および対象ガス供給量（ＦＳｔ）から、対象メタン価（ＭＮｔ、対象ガスエンジン２ｔにおける気筒２１に供給されたガス燃料のメタン価）を推定できる。エンジン負荷（ＥＬ）やガス供給量（ＦＳ）は、対象ガスエンジン（２ｔ）の運転時において対象ガスエンジン（２ｔ）に通常装着されるセンサなどから取得可能である。よって、上記７）の方法によれば、対象ガスエンジン（２ｔ）の運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価（ＭＮｔ）を推定可能である。

８）本開示の少なくとも一実施形態にかかるメタン価の推定方法（３００）は、
少なくとも１つの気筒（２１）を備えるガスエンジン（２）である対象ガスエンジン（２ｔ）に供給されるガス燃料のメタン価（ＭＮｔ）を推定するメタン価の推定方法（３００）であって、
ガスエンジン（２）における、低位発熱量（ＬＨＶ）と供給されるガス燃料のメタン価（ＭＮ）とが予め関連付けられた第３の関連付け情報（ＡＩ３）、を取得する第３の関連付け情報取得ステップ（Ｓ３０１）と、
前記対象ガスエンジン（２ｔ）における１サイクルあたりの熱発生量（ＱＣ）である対象熱発生量（ＱＣｔ）を取得する熱発生量取得ステップ（Ｓ３０２）と、
前記対象熱発生量（ＱＣｔ）を取得した期間における前記対象ガスエンジン（２ｔ）に供給されるガス燃料の１サイクルあたりの供給量（ＭＦ）である対象供給量（ＭＦｔ）を取得する供給量取得ステップ（Ｓ３０３）と、
前記対象熱発生量（ＱＣｔ）と前記対象供給量（ＭＦｔ）から、前記対象ガスエンジン（２ｔ）における低位発熱量（ＬＨＶ）である対象低位発熱量（ＬＨＶｔ）を算出する低位発熱量算出ステップ（Ｓ３０４）と、
前記第３の関連付け情報（ＡＩ１）に基づいて、前記対象低位発熱量（ＬＨＶｔ）から、前記対象ガスエンジン（２ｔ）における前記気筒（２１）に供給された前記ガス燃料のメタン価（ＭＮｔ）を推定するメタン価推定ステップ（Ｓ３０５）と、を備える。

上記８）の方法によれば、ガスエンジン（２）における、低位発熱量（ＬＨＶ）と供給されるガス燃料のメタン価（ＭＮ）との関係性を示す第３の関連付け情報（ＡＩ３）を利用することで、対象低位発熱量（ＬＨＶｔ）から、対象メタン価（ＭＮｔ、対象ガスエンジン２ｔにおける気筒２１に供給されたガス燃料のメタン価）を推定できる。低位発熱量（ＬＨＶ）は、対象ガスエンジン（２ｔ）の運転時において対象ガスエンジン（２ｔ）に通常装着されるセンサなどから取得できる、上記熱発生量（ＱＣ）および上記供給量（ＭＦ）から算出できる。よって、上記８）の方法によれば、対象ガスエンジン（２ｔ）の運転時に、カロリーメータなどの特別なセンサを用いることなく対象メタン価（ＭＮｔ）を推定可能である。

１，１Ａ～１Ｃ メタン価の推定装置
２ ガスエンジン
２ｔ 対象ガスエンジン
３ 制御装置
３ｔ 対象制御装置
４ エンジンシステム
５ ガス燃料供給系統
１１ 第１の関連付け情報取得部
１２ 点火時期取得部
１３ ノッキング強度取得部
１４ メタン価推定部
２１ 気筒
２２ 駆動シャフト
３１ エンジン制御部
３２ 燃焼制御部
３３ 燃焼診断部
３４ ガス燃料供給量制御部
３５ 点火時期制御部
４１ 燃料噴射弁
４１ ガス燃料噴射装置
４２ 点火装置
４３ 発電機
４４ ノッキングセンサ
４５ 筒内圧センサ
４６ データベース部
５１ 第１のガス燃料貯留装置
５２ 第２のガス燃料貯留装置
５３ 第１のガス燃料供給ライン
５４ 第２のガス燃料供給ライン
５５ 切り替え装置
５６ 合流部
６１ 第２の関連付け情報取得部
６２ エンジン負荷取得部
６３ ガス供給量取得部
６４ メタン価推定部
７１ 第３の関連付け情報取得部
７２ 熱発生量取得部
７３ 供給量取得部
７４ 低位発熱量算出部
７５ メタン価推定部
１００，２００，３００ メタン価の推定方法
ＡＩ１ 第１の関連付け情報
ＡＩ２ 第２の関連付け情報
ＡＩ３ 第３の関連付け情報
ＥＬ エンジン負荷
ＥＬｔ 対象エンジン負荷
ＦＳ ガス供給量
ＦＳｔ 対象ガス供給量
ＩＴ，ＩＴｃ１，ＩＴｃ２ 点火時期
ＩＴｓ 設定点火時期
ＩＴｔ，ＩＴｔ１ 対象点火時期
ＫＩ，ＫＩｓ ノッキング強度
ＫＩｔ，ＫＩｔ１ 対象ノッキング強度
ＬＨＶ 低位発熱量
ＬＨＶｔ 対象低位発熱量
ＭＦ 供給量
ＭＦｔ 対象供給量
ＭＮ メタン価
ＭＮｔ 対象メタン価
ＱＣ 熱発生量
ＱＣｔ 対象熱発生量
Ｓ１０１ 第１の関連付け情報取得ステップ
Ｓ１０２ 点火時期取得ステップ
Ｓ１０３ ノッキング強度取得ステップ
Ｓ１０４ メタン価推定ステップ
Ｓ２０１ 第２の関連付け情報取得ステップ
Ｓ２０２ エンジン負荷取得ステップ
Ｓ２０３ ガス供給量取得ステップ
Ｓ２０４ メタン価推定ステップ
Ｓ３０１ 第３の関連付け情報取得ステップ
Ｓ３０２ 熱発生量取得ステップ
Ｓ３０３ 供給量取得ステップ
Ｓ３０４ 低位発熱量算出ステップ
Ｓ３０５ メタン価推定ステップ

Claims (8)

1. 少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定装置であって、
   ガスエンジンにおける、点火時期とノッキング強度と供給されるガス燃料のメタン価とが予め関連付けられた第１の関連付け情報、を取得する第１の関連付け情報取得部と、
   前記対象ガスエンジンにおける前記点火時期である対象点火時期を取得する点火時期取得部と、
   前記対象ガスエンジンの前記対象点火時期を含むサイクルにおける前記ノッキング強度である対象ノッキング強度を取得するノッキング強度取得部と、
   前記第１の関連付け情報に基づいて、前記対象点火時期および前記対象ノッキング強度から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定部と、を備える、
   メタン価の推定装置。
2. 前記メタン価推定部は、
   前記対象ガスエンジンのエンジン負荷が一定である一定負荷条件において、特定のノッキング強度に到達したときの前記対象ノッキング強度、および前記特定のノッキング強度に到達したときの前記対象ノッキング強度を含むサイクルにおける前記対象点火時期から、前記第１の関連付け情報に基づいて、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定する、
   請求項１に記載のメタン価の推定装置。
3. 少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定装置であって、
   ガスエンジンにおける、エンジン負荷と供給されるガス燃料のガス供給量と前記ガス燃料のメタン価とが予め関連付けられた第２の関連付け情報、を取得する第２の関連付け情報取得部と、
   前記対象ガスエンジンにおける前記エンジン負荷である対象エンジン負荷を取得するエンジン負荷取得部と、
   前記対象ガスエンジンの前記対象エンジン負荷を取得した期間における前記ガス供給量である対象ガス供給量を取得するガス供給量取得部と、
   前記第２の関連付け情報に基づいて、前記対象エンジン負荷および前記対象ガス供給量から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定部と、を備える、
   メタン価の推定装置。
4. 少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定装置であって、
   ガスエンジンにおける、低位発熱量と供給されるガス燃料のメタン価とが予め関連付けられた第３の関連付け情報、を取得する第３の関連付け情報取得部と、
   前記対象ガスエンジンにおける１サイクルあたりの熱発生量である対象熱発生量を取得する熱発生量取得部と、
   前記対象熱発生量を取得した期間における前記対象ガスエンジンに供給されるガス燃料の１サイクルあたりの供給量である対象供給量を取得する供給量取得部と、
   前記対象熱発生量と前記対象供給量から、前記対象ガスエンジンにおける低位発熱量である対象低位発熱量を算出する低位発熱量算出部と、
   前記第３の関連付け情報に基づいて、前記対象低位発熱量から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定部と、を備える、
   メタン価の推定装置。
5. 請求項１乃至４の何れか１項に記載のメタン価の推定装置と、
   前記メタン価の推定装置で推定された前記ガス燃料のメタン価に応じて、前記対象ガスエンジンの点火時期である対象点火時期を制御するように構成された点火時期制御部と、を備える、
   ガスエンジンの制御装置。
6. 少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定方法であって、
   ガスエンジンにおける、点火時期とノッキング強度と供給されるガス燃料のメタン価とが予め関連付けられた第１の関連付け情報、を取得する第１の関連付け情報取得ステップと、
   前記対象ガスエンジンにおける前記点火時期である対象点火時期を取得する点火時期取得ステップと、
   前記対象ガスエンジンの前記対象点火時期を含むサイクルにおける前記ノッキング強度である対象ノッキング強度を取得するノッキング強度取得ステップと、
   前記第１の関連付け情報に基づいて、前記対象点火時期および前記対象ノッキング強度から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定ステップと、を備える、
   メタン価の推定方法。
7. 少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定方法であって、
   ガスエンジンにおける、エンジン負荷と供給されるガス燃料のガス供給量と前記ガス燃料のメタン価とが予め関連付けられた第２の関連付け情報、を取得する第２の関連付け情報取得ステップと、
   前記対象ガスエンジンにおける前記エンジン負荷である対象エンジン負荷を取得するエンジン負荷取得ステップと、
   前記対象ガスエンジンの前記対象エンジン負荷を取得した期間における前記ガス供給量である対象ガス供給量を取得するガス供給量取得ステップと、
   前記第２の関連付け情報に基づいて、前記対象エンジン負荷および前記対象ガス供給量から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定ステップと、を備える、
   メタン価の推定方法。
8. 少なくとも１つの気筒を備えるガスエンジンである対象ガスエンジンに供給されるガス燃料のメタン価を推定するメタン価の推定方法であって、
   ガスエンジンにおける、低位発熱量と供給されるガス燃料のメタン価とが予め関連付けられた第３の関連付け情報、を取得する第３の関連付け情報取得ステップと、
   前記対象ガスエンジンにおける１サイクルあたりの熱発生量である対象熱発生量を取得する熱発生量取得ステップと、
   前記対象熱発生量を取得した期間における前記対象ガスエンジンに供給されるガス燃料の１サイクルあたりの供給量である対象供給量を取得する供給量取得ステップと、
   前記対象熱発生量と前記対象供給量から、前記対象ガスエンジンにおける低位発熱量である対象低位発熱量を算出する低位発熱量算出ステップと、
   前記第３の関連付け情報に基づいて、前記対象低位発熱量から、前記対象ガスエンジンにおける前記気筒に供給された前記ガス燃料のメタン価を推定するメタン価推定ステップと、を備える、
   メタン価の推定方法。

Images