JPWO2016104270A1 - メタン価算出方法およびメタン価測定装置 - Google Patents
メタン価算出方法およびメタン価測定装置 Download PDFInfo
- Publication number
- JPWO2016104270A1 JPWO2016104270A1 JP2016566149A JP2016566149A JPWO2016104270A1 JP WO2016104270 A1 JPWO2016104270 A1 JP WO2016104270A1 JP 2016566149 A JP2016566149 A JP 2016566149A JP 2016566149 A JP2016566149 A JP 2016566149A JP WO2016104270 A1 JPWO2016104270 A1 JP WO2016104270A1
- Authority
- JP
- Japan
- Prior art keywords
- value
- methane number
- gas
- measurement target
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 674
- 238000004364 calculation method Methods 0.000 title claims description 65
- 239000007789 gas Substances 0.000 claims abstract description 259
- 230000014509 gene expression Effects 0.000 claims abstract description 107
- 238000000034 method Methods 0.000 claims abstract description 65
- 239000003345 natural gas Substances 0.000 claims abstract description 40
- 238000005259 measurement Methods 0.000 claims description 140
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 64
- 229910001873 dinitrogen Inorganic materials 0.000 claims description 52
- 239000003949 liquefied natural gas Substances 0.000 claims description 37
- 238000006243 chemical reaction Methods 0.000 claims description 28
- 230000008016 vaporization Effects 0.000 claims description 7
- 239000000203 mixture Substances 0.000 abstract description 29
- 239000002737 fuel gas Substances 0.000 abstract description 16
- 239000000446 fuel Substances 0.000 abstract description 14
- 238000012937 correction Methods 0.000 description 12
- TVCSLGBWHLPONF-UHFFFAOYSA-N 6-(hydroxymethyl)-9-methoxyphenazine-1-carboxylic acid Chemical compound C1=CC(C(O)=O)=C2N=C3C(OC)=CC=C(CO)C3=NC2=C1 TVCSLGBWHLPONF-UHFFFAOYSA-N 0.000 description 11
- 238000002485 combustion reaction Methods 0.000 description 6
- 229910052757 nitrogen Inorganic materials 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000005352 clarification Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 239000003502 gasoline Substances 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- TVMXDCGIABBOFY-UHFFFAOYSA-N octane Chemical compound CCCCCCCC TVMXDCGIABBOFY-UHFFFAOYSA-N 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- -1 specifically Substances 0.000 description 1
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C10—PETROLEUM, GAS OR COKE INDUSTRIES; TECHNICAL GASES CONTAINING CARBON MONOXIDE; FUELS; LUBRICANTS; PEAT
- C10L—FUELS NOT OTHERWISE PROVIDED FOR; NATURAL GAS; SYNTHETIC NATURAL GAS OBTAINED BY PROCESSES NOT COVERED BY SUBCLASSES C10G, C10K; LIQUEFIED PETROLEUM GAS; ADDING MATERIALS TO FUELS OR FIRES TO REDUCE SMOKE OR UNDESIRABLE DEPOSITS OR TO FACILITATE SOOT REMOVAL; FIRELIGHTERS
- C10L3/00—Gaseous fuels; Natural gas; Synthetic natural gas obtained by processes not covered by subclass C10G, C10K; Liquefied petroleum gas
- C10L3/06—Natural gas; Synthetic natural gas obtained by processes not covered by C10G, C10K3/02 or C10K3/04
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/22—Fuels; Explosives
- G01N33/225—Gaseous fuels, e.g. natural gas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K17/00—Measuring quantity of heat
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Physics & Mathematics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Food Science & Technology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Combustion & Propulsion (AREA)
- Organic Chemistry (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本発明は、測定対象ガスである天然ガスについて、一応の信頼性を有するメタン価をガス組成に拘わらず容易に得ることができるメタン価算出方法および燃料ガスとして利用される天然ガスの燃料性状の監視を行うことのできるメタン価測定装置を提供することを目的とする。本発明においては、各々天然ガスよりなる互いにメタン価の値が異なる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式を予め取得しておき、測定対象ガスである天然ガスの基礎熱量を測定し、測定された当該測定対象ガスの基礎熱量の値と、特定の関係式とから測定対象ガスのメタン価を算出する。
Description
本発明は、メタン価算出方法およびメタン価測定装置に関する。
近年、窒素酸化物（ＮＯx ）の排出量の抑制やＣＯ2 排出量低減などを目的として、液化天然ガス（ＬＮＧ）を船舶用燃料として利用したガスエンジンの開発、導入が進んでいる。
ＬＮＧを利用した船舶の実用化に係る課題としては、ＬＮＧの組成が産出地により異なる点や、ガスエンジンの起動時や負荷変動時に燃料ガスの消費量の変動が生ずることによって燃料ガスの組成が変化する点などが挙げられる。燃料ガスの組成が変わると、燃料ガスの熱量やメタン価などの特性が変化し、エンジンのノッキングや失火などの異常燃焼を引き起こす原因になりうる。ここで、メタン価とは、ガソリンエンジンのオクタン価に対応するノッキングに対する抵抗値を示す指標であって、純メタンを１００、水素を０として評価した指標である。
ＬＮＧを利用した船舶の実用化に係る課題としては、ＬＮＧの組成が産出地により異なる点や、ガスエンジンの起動時や負荷変動時に燃料ガスの消費量の変動が生ずることによって燃料ガスの組成が変化する点などが挙げられる。燃料ガスの組成が変わると、燃料ガスの熱量やメタン価などの特性が変化し、エンジンのノッキングや失火などの異常燃焼を引き起こす原因になりうる。ここで、メタン価とは、ガソリンエンジンのオクタン価に対応するノッキングに対する抵抗値を示す指標であって、純メタンを１００、水素を０として評価した指標である。
このような異常燃焼を回避するためには、燃料ガスの熱量やメタン価などの燃料性状をリアルタイムに捉え、このデータに基づきガスエンジンの燃焼制御を行うことが有効な手段であると考えられる。
例えばＬＮＧなどの燃料ガスの熱量を測定する方法としては、例えば熱量と特定の対応関係を有する物性値を測定し、測定値に基づいて熱量の値（換算熱量）を求める方法などが、本件出願人によって提案されている（例えば特許文献１参照。）。
一方、燃料ガスのメタン価を算出する方法としては、
（ａ）ＡＶＬ社が提案する方式（以下、「ＡＶＬ基準」ともいう。）、
（ｂ）カリフォルニア州大気資源評議会で規定された特定の演算式により算出する方式（以下、「ＣＡＲＢ基準」ともいう。）、
（ｃ）ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出する方式（以下、「ＧＲＩ（Ｌｃ）基準」ともいう。）、
（ｄ）ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出する方式（以下、「ＧＲＩ（Ｈ／Ｃ）基準」ともいう。）
の４種類が主に利用されている。ここに、メタン価は、同一の燃料ガスであっても、算出方法によって異なる値を示し、例えば地域ごとに異なる基準によるメタン価が要請されている。
しかしながら、いずれの方法も、メタン価をガス組成に基づいて算出するものであることから、上述したように、ガス組成の変動が生じた場合には、メタン価の算出にあたってガス組成を測定することが必要となる。
（ａ）ＡＶＬ社が提案する方式（以下、「ＡＶＬ基準」ともいう。）、
（ｂ）カリフォルニア州大気資源評議会で規定された特定の演算式により算出する方式（以下、「ＣＡＲＢ基準」ともいう。）、
（ｃ）ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出する方式（以下、「ＧＲＩ（Ｌｃ）基準」ともいう。）、
（ｄ）ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出する方式（以下、「ＧＲＩ（Ｈ／Ｃ）基準」ともいう。）
の４種類が主に利用されている。ここに、メタン価は、同一の燃料ガスであっても、算出方法によって異なる値を示し、例えば地域ごとに異なる基準によるメタン価が要請されている。
しかしながら、いずれの方法も、メタン価をガス組成に基づいて算出するものであることから、上述したように、ガス組成の変動が生じた場合には、メタン価の算出にあたってガス組成を測定することが必要となる。
而して、例えばＬＮＧを気化させて得られたガス（以下、「ＬＮＧ気化ガス」という。）には、通常、不燃ガス成分も含まれており、当該不燃ガス成分による発熱量に対する影響の程度に法則性がないことから、ＬＮＧ気化ガスの発熱量（真の熱量）とメタン価の値との間には、特定の相関関係は成立しない。
然るに、本発明者らが鋭意研究を重ねた結果、燃料ガスとして利用される天然ガスの基礎熱量の値と、上記の各基準により算出されるメタン価の値との間に特定の相関関係が成立し、測定対象ガスである天然ガスの基礎熱量を測定することにより、各々の基準に対応するメタン価の近似解を得ることができることを見出した。ここに、「基礎熱量」とは、天然ガスから不燃ガス成分を除いたときの燃焼性ガス成分の燃焼熱量をいう。例えばＬＮＧ気化ガスの場合には、不燃ガス成分をＮ2 だけと見なすことができるため、ＬＮＧ気化ガスの基礎熱量は、ＬＮＧ気化ガスからＮ2 を除いた場合の燃焼熱量をいう。
然るに、本発明者らが鋭意研究を重ねた結果、燃料ガスとして利用される天然ガスの基礎熱量の値と、上記の各基準により算出されるメタン価の値との間に特定の相関関係が成立し、測定対象ガスである天然ガスの基礎熱量を測定することにより、各々の基準に対応するメタン価の近似解を得ることができることを見出した。ここに、「基礎熱量」とは、天然ガスから不燃ガス成分を除いたときの燃焼性ガス成分の燃焼熱量をいう。例えばＬＮＧ気化ガスの場合には、不燃ガス成分をＮ2 だけと見なすことができるため、ＬＮＧ気化ガスの基礎熱量は、ＬＮＧ気化ガスからＮ2 を除いた場合の燃焼熱量をいう。
本発明は、以上のような事情に基づいてなされたものであって、測定対象ガスである天然ガスについて、一応の信頼性を有するメタン価をガス組成に拘わらず容易に得ることのできるメタン価算出方法を提供することを目的とする。
また、本発明は、測定対象ガスである天然ガスについて、一応の信頼性を有するメタン価をガス組成に拘わらず容易に得ることができ、燃料ガスとして利用される天然ガスの燃料性状の監視を行うことのできるメタン価測定装置を提供することを目的とする。
また、本発明は、測定対象ガスである天然ガスについて、一応の信頼性を有するメタン価をガス組成に拘わらず容易に得ることができ、燃料ガスとして利用される天然ガスの燃料性状の監視を行うことのできるメタン価測定装置を提供することを目的とする。
本発明のメタン価算出方法は、各々天然ガスよりなる互いにメタン価の値が異なる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式を予め取得しておき、
測定対象ガスである天然ガスの基礎熱量を測定し、
測定された当該測定対象ガスの基礎熱量の値と、前記特定の関係式とから当該測定対象ガスのメタン価を算出することを特徴とする。
測定対象ガスである天然ガスの基礎熱量を測定し、
測定された当該測定対象ガスの基礎熱量の値と、前記特定の関係式とから当該測定対象ガスのメタン価を算出することを特徴とする。
本発明のメタン価算出方法においては、前記特定の関係式が下記式（１）で表されるものであることが好ましい。
また、本発明のメタン価算出方法においては、カリフォルニア州大気資源評議会で規定された特定の演算式により算出されるメタン価の近似値を取得するに際しては、前記特定の関係式として下記式（２）で表されるものが用いられる。また、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得するに際しては、前記特定の関係式として下記式（３）で表されるものが用いられる。さらにまた、ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出されるメタン価の近似値を取得するに際しては、前記特定の関係式として下記式（４）で表されるものが用いられる。
本発明のメタン価算出方法は、各々窒素ガスを含有する天然ガスよりなる互いにメタン価の値が異なる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式を予め取得しておき、
測定対象ガスである窒素ガスを含有する天然ガスの基礎熱量および当該測定対象ガスに含まれる窒素ガス濃度を測定し、
測定された当該測定対象ガスの基礎熱量の値および窒素ガスの濃度値と、前記特定の関係式とから当該測定対象ガスのメタン価を算出することを特徴とする。
測定対象ガスである窒素ガスを含有する天然ガスの基礎熱量および当該測定対象ガスに含まれる窒素ガス濃度を測定し、
測定された当該測定対象ガスの基礎熱量の値および窒素ガスの濃度値と、前記特定の関係式とから当該測定対象ガスのメタン価を算出することを特徴とする。
本発明のメタン価算出方法においては、前記特定の関係式が下記式（５）で表されるものであることが好ましい。
また、本発明のメタン価算出方法においては、窒素ガスを含有する天然ガスについて、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得するに際しては、前記特定の関係式として下記式（５）で表されるものが用いられることが好ましい。
以上の本発明のメタン価算出方法においては、前記測定対象ガスである天然ガスは、液化天然ガスを気化させて得られたものであることが好ましい。
さらにまた、本発明のメタン価算出方法においては、前記測定対象ガスの基礎熱量は、当該測定対象ガスの屈折率から得られる屈折率換算熱量と、当該測定対象ガスの音速から得られる音速換算熱量とに基づいて得られるものであることが好ましい。
本発明のメタン価測定装置は、測定対象ガスである天然ガスの基礎熱量を測定する熱量測定機構と、
予め取得しておいた、各々メタン価の値が異なる天然ガスよりなる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式と、前記熱量測定機構により測定された当該測定対象ガスの基礎熱量の値とから、当該測定対象ガスのメタン価を算出するメタン価算出機構と
を備えることを特徴とする。
予め取得しておいた、各々メタン価の値が異なる天然ガスよりなる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式と、前記熱量測定機構により測定された当該測定対象ガスの基礎熱量の値とから、当該測定対象ガスのメタン価を算出するメタン価算出機構と
を備えることを特徴とする。
本発明のメタン価測定装置においては、前記特定の関係式が上記式（１）で表されるものであることが好ましい。
さらにまた、本発明のメタン価測定装置においては、カリフォルニア州大気資源評議会で規定された特定の演算式により算出されるメタン価の近似値を取得する上記式（２）で表される特定の関係式、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得する上記式（３）で表される特定の関係式、および、ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出されるメタン価の近似値を取得する上記式（４）で表される特定の関係式がさらに設定されていることが好ましい。
本発明のメタン価測定装置は、測定対象ガスである窒素ガスを含有する天然ガスの基礎熱量を測定する熱量測定機構と、
当該測定対象ガスに含まれる窒素ガス濃度を算出する濃度測定機構と、
予め取得しておいた、各々メタン価の値が異なる窒素ガスを含有する天然ガスよりなる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式と、前記熱量測定機構により測定された当該測定対象ガスの基礎熱量の値および前記濃度測定機構により測定された窒素ガスの濃度値とから、当該測定対象ガスのメタン価を算出するメタン価算出機構と
を備えることを特徴とする。
当該測定対象ガスに含まれる窒素ガス濃度を算出する濃度測定機構と、
予め取得しておいた、各々メタン価の値が異なる窒素ガスを含有する天然ガスよりなる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式と、前記熱量測定機構により測定された当該測定対象ガスの基礎熱量の値および前記濃度測定機構により測定された窒素ガスの濃度値とから、当該測定対象ガスのメタン価を算出するメタン価算出機構と
を備えることを特徴とする。
本発明のメタン価測定装置においては、前記特定の関係式が上記式（５）で表されるものであることが好ましい。
さらにまた、本発明のメタン価測定装置においては、カリフォルニア州大気資源評議会で規定された特定の演算式により算出されるメタン価の近似値を取得する上記式（２）で表される特定の関係式、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得する上記式（６）で表される特定の関係式、および、ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出されるメタン価の近似値を取得する上記式（４）で表される特定の関係式がさらに設定されていることが好ましい。
以上の本発明のメタン価測定装置においては、前記測定対象ガスである天然ガスは、液化天然ガスを気化させて得られたものであることが好ましい。
さらにまた、本発明のメタン価測定装置においては、前記熱量測定機構が、測定対象ガスの屈折率の値から屈折率換算熱量を求める屈折率換算熱量測定手段と、当該測定対象ガスの音速の値から音速換算熱量を求める音速換算熱量測定手段と、当該測定対象ガスの基礎熱量を屈折率換算熱量と音速換算熱量とに基づいて算出する熱量算出手段とを備えた構成とされていることが好ましい。
本発明のメタン価算出方法によれば、予め取得された基礎熱量とメタン価との特定の関係式を利用することにより、測定対象ガスの基礎熱量の値を測定するだけで、当該測定対象ガスのメタン価を求めることができる。特定の関係式は、各々メタン価の値が互いに異なる天然ガスよりなる複数種の基準ガスについて、基礎熱量とメタン価との相関関係を実験による裏づけにより定量的に明らかにしたものであるため、得られるメタン価は一応の信頼性を有するものとなる。
さらに、測定対象ガスに含まれる窒素ガス濃度に基づいた補正がなされることにより、得られるメタン価は一層信頼性の高いものとなる。
さらに、測定対象ガスに含まれる窒素ガス濃度に基づいた補正がなされることにより、得られるメタン価は一層信頼性の高いものとなる。
上記のメタン価算出方法が実行される本発明のメタン価測定装置によれば、熱量測定機構によって測定対象ガスの基礎熱量を連続的に測定することにより、実際の状況に即した測定対象ガスのメタン価を連続的に取得することができるので、燃料ガスとしての天然ガスの実際の燃料性状の監視を行うことができる。従って、ガス組成の変動が生じた場合には、ガス組成の変動に伴うメタン価の変動を速やかに検出することができる。
以下、本発明の実施の形態について詳細に説明する。
本発明においては、天然ガス、具体的には例えばＬＮＧ気化ガスを測定対象ガスとし、窒素ガスなどの不燃ガス成分を含んだものであってもよく、また、ＬＮＧ気化ガスの精製過程においてメタンより重質の炭化水素成分が除去またはその含有量が調整されたものも含む。
本発明においては、天然ガス、具体的には例えばＬＮＧ気化ガスを測定対象ガスとし、窒素ガスなどの不燃ガス成分を含んだものであってもよく、また、ＬＮＧ気化ガスの精製過程においてメタンより重質の炭化水素成分が除去またはその含有量が調整されたものも含む。
図１は、本発明のメタン価測定装置の一例における構成の概略を示すブロック図である。
このメタン価測定装置は、測定対象ガスの熱量を測定する熱量測定機構２０と、当該測定対象ガスのメタン価を算出するメタン価算出機構４０と、測定対象ガスの熱量およびメタン価などの情報を表示する表示機構４５とが、例えば防爆性容器１０内に配設されて構成されている。
このメタン価測定装置は、測定対象ガスの熱量を測定する熱量測定機構２０と、当該測定対象ガスのメタン価を算出するメタン価算出機構４０と、測定対象ガスの熱量およびメタン価などの情報を表示する表示機構４５とが、例えば防爆性容器１０内に配設されて構成されている。
熱量測定機構２０は、例えば、測定対象ガスの音速の値から求められる音速換算熱量Ｑｓを得るための音速換算熱量測定機構２５と、当該測定対象ガスの屈折率の値から求められる屈折率換算熱量Ｑｎを得るための屈折率換算熱量測定機構２１と、当該測定対象ガスに含まれる窒素ガス濃度ＸN2 〔ｖｏｌ％〕を測定する窒素濃度測定機構３０と、当該測定対象ガスの熱量Ｑｎｓの値および基礎熱量Ｑ´の値を算出する熱量計算機構３５とを備えている。
音速換算熱量測定機構２５は、測定対象ガス中における音波の伝播速度（測定対象ガスの音速）を測定する音速測定手段２６と、音速測定手段２６によって測定された音速の値に基づいて音速換算熱量Ｑｓの値を求める機能を有する音速−熱量換算処理手段２７とを備えている。
音速−熱量換算処理手段２７は、測定対象ガスであるＬＮＧ気化ガスにおいて不燃ガス成分（Ｎ2 ）を含まない燃焼性ガス成分（パラフィン系炭化水素ガス）のみからなる特定ガスについて、例えばグラフ化することなどによって予め取得された音速と熱量との相関関係を利用し、当該相関関係に対して、測定対象ガスについて得られた音速の値が特定ガスの音速である仮定して対照することにより音速換算熱量Ｑｓを算出する。
音速−熱量換算処理手段２７は、測定対象ガスであるＬＮＧ気化ガスにおいて不燃ガス成分（Ｎ2 ）を含まない燃焼性ガス成分（パラフィン系炭化水素ガス）のみからなる特定ガスについて、例えばグラフ化することなどによって予め取得された音速と熱量との相関関係を利用し、当該相関関係に対して、測定対象ガスについて得られた音速の値が特定ガスの音速である仮定して対照することにより音速換算熱量Ｑｓを算出する。
屈折率換算熱量測定機構２１は、測定対象ガスの屈折率を測定する屈折率測定手段２２と、屈折率測定手段２２によって測定された屈折率の値に基づいて屈折率換算熱量Ｑｎを求める機能を有する屈折率−熱量換算処理手段２３とを備えている。
屈折率−熱量換算処理手段２３は、測定対象ガスであるＬＮＧ気化ガスにおいて不燃ガス成分（Ｎ2 ）を含まない燃焼性ガス成分（パラフィン系炭化水素ガス）のみからなる特定ガスについて、例えばグラフ化することなどによって予め取得された屈折率と熱量との相関関係を利用し、当該相関関係に対して、測定対象ガスについて得られた屈折率の値が特定ガスの屈折率であると仮定して対照することにより屈折率換算熱量Ｑｎを算出する。
屈折率−熱量換算処理手段２３は、測定対象ガスであるＬＮＧ気化ガスにおいて不燃ガス成分（Ｎ2 ）を含まない燃焼性ガス成分（パラフィン系炭化水素ガス）のみからなる特定ガスについて、例えばグラフ化することなどによって予め取得された屈折率と熱量との相関関係を利用し、当該相関関係に対して、測定対象ガスについて得られた屈折率の値が特定ガスの屈折率であると仮定して対照することにより屈折率換算熱量Ｑｎを算出する。
窒素濃度測定機構３０は、音速換算熱量測定機構２５によって得られた音速換算熱量Ｑｓの値と、屈折率換算熱量測定機構２１によって得られた屈折率換算熱量Ｑｎの値とに基づいて、下記式（７）により、補正因子αとして１．１〜４．２の範囲内、好ましくは２．４０〜２．６０の範囲内において選択される値を用いる条件にて測定対象ガスに含まれる窒素ガスの濃度を算出する。式（７）におけるＸN2 は、体積百分率で表される窒素ガス濃度〔ｖｏｌ％〕である。ＫN2 は誤差係数であって、雑ガス成分としてのＮ2 が及ぼす誤差の影響の大きさを表す。計算に供される音速換算熱量Ｑｓの値および屈折率換算熱量Ｑｎの値の単位は、〔ＭＪ／ｍ3 〕である。
熱量計算機構３５は、音速換算熱量測定機構２５によって得られた音速換算熱量Ｑｓの値と、屈折率換算熱量測定機構２１によって得られた屈折率換算熱量Ｑｎの値とに基づいて、下記式（８）により、補正因子αとして１．１〜４．２の範囲内、好ましくは２．４０〜２．６０の範囲内において選択される値を用いる条件にて測定対象ガスの熱量Ｑｎｓの値を算出する。このようにして得られた熱量Ｑｎｓの値と、窒素濃度測定機構３０によって得られた窒素ガスの濃度ＸN2 の値とに基づいて、下記式（９）により、基礎熱量Ｑ´の値〔ＭＪ／ｍ3 〕を算出する。
メタン価算出機構４０は、測定対象ガスであるＬＮＧ気化ガスについて、上述した４つの基準（ａ）〜（ｄ）のうちから選択された基準に基づく方法により得られるメタン価の値の近似解を算出する。以下、例えば窒素ガスを含まないＬＮＧ気化ガスについて、（ａ）ＡＶＬ基準に基づく方法により得られるメタン価の値（以下、「ＡＶＬ値」ともいう。）の近似解を算出する場合を例に挙げて具体的に説明する。
メタン価算出機構４０は、熱量測定機構２０により測定された当該測定対象ガスの基礎熱量Ｑ´の値と、予め取得しておいた、各々メタン価の値が異なる天然ガスよりなる複数種の基準ガスについて、ＡＶＬ基準に基づく方法により得られるメタン価の値（ＡＶＬ値）と基礎熱量Ｑ´の値との特定の関係式とから、測定対象ガスのメタン価を算出する。
特定の関係式は、例えば図２に示すように、横軸を基礎熱量とし、縦軸をメタン価とする座標系において、前記複数種の基準ガスの各々について、基礎熱量Ｑ´の値とＡＶＬ値との関係を示す実測値を取得し、得られた実測値を例えば多項式で曲線近似することにより取得することができる。具体的には、ＡＶＬ基準に係る特定の関係式は、上記式（１）で表されるものであることが好ましい。
上記式（１）において、ＭＮはメタン価、具体的にはＡＶＬ値の近似解、ｆ(Q´) は、測定対象ガスの基礎熱量Ｑ´の値に応じて選択される上記式（ａ）〜上記式（ｄ）で表されるいずれかの関数である。
上記式（１）におけるＡは、−２．０〜２．０の範囲から選択される値である。このＡについて設定される数値範囲は、式（１）においてＡ＝０とした式で表される基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ａの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出される近似解のＡＶＬ値に対する誤差率が５．０％以内となり、高い信頼度が得られる。
上記式（１）におけるＡの値の具体的な設定方法としては、例えば、組成が既知の基準ガスについてメタン価を測定し、理論値（ＡＶＬ値）との差分を「Ａ」として設定（オフセット調整）することができる。
上記式（１）におけるＡは、−２．０〜２．０の範囲から選択される値である。このＡについて設定される数値範囲は、式（１）においてＡ＝０とした式で表される基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ａの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出される近似解のＡＶＬ値に対する誤差率が５．０％以内となり、高い信頼度が得られる。
上記式（１）におけるＡの値の具体的な設定方法としては、例えば、組成が既知の基準ガスについてメタン価を測定し、理論値（ＡＶＬ値）との差分を「Ａ」として設定（オフセット調整）することができる。
図２において実線で示す、基準となるメタン価算出曲線は、上記式（ａ）〜上記式（ｄ）で表される曲線部分の各々が変曲点を生じることなく連続するものである。この基準となるメタン価算出曲線は、例えば実在し得るＬＮＧ気化ガスのガス組成を考慮して、設定されている。
図２において破線で示す曲線は、式（１）においてＡ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（１）においてＡ＝−２．０としたメタン価算出曲線である。
このようなメタン価算出曲線を利用することにより、算出されるメタン価の信頼性を一層高いものとすることができる。
図２において破線で示す曲線は、式（１）においてＡ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（１）においてＡ＝−２．０としたメタン価算出曲線である。
このようなメタン価算出曲線を利用することにより、算出されるメタン価の信頼性を一層高いものとすることができる。
以上において、図１における１１は、測定対象ガスを音速測定手段２６および屈折率測定手段２２の各々に供給するための測定対象ガス導入部、１２は、屈折率測定手段２２において検知原理上必要とされる参照ガスを導入するための参照ガス導入部、１３はガス排出部である。また、図１における二点鎖線は、ガス配管を示す。
上記のメタン価測定装置は、例えば、適宜のガスサンプリング装置を介してガスパイプラインに配管接続され、ガスパイプライン内を流通するＬＮＧ気化ガスが測定対象ガスとして測定対象ガス導入部１１から音速換算熱量測定機構２５の音速測定手段２６および屈折率換算熱量測定機構２１の屈折率測定手段２２の各々に順次に供給される。一方、例えば空気などの参照ガスが参照ガス導入部１２から屈折率換算熱量測定機構２１の屈折率測定手段２２に供給される。これにより、音速換算熱量測定機構２５においては音速換算熱量Ｑｓが求められる共に、屈折率換算熱量測定機構２１においては屈折率換算熱量Ｑｎが求められる。そして、音速換算熱量Ｑｓの値と、屈折率換算熱量Ｑｎの値とに基づいて、上記式（７）および上記式（８）により、補正因子αとして特定の範囲内において選択された値を用いて、窒素ガス濃度ＸN2〔ｖｏｌ％〕および熱量Ｑｎｓが算出される。このようにして得られた熱量Ｑｎｓの値と、窒素ガス濃度ＸN2〔ｖｏｌ％〕の値とに基づいて、上記式（９）により、測定対象ガスの基礎熱量Ｑ´が算出される。
次いで、メタン価算出機構４０によって、熱量測定機構２０によって得られた基礎熱量Ｑ´の値と、上記の特定の関係式、例えばＡ＝０とした基準となるメタン価算出曲線を表す式（１）とから、ＡＶＬ値の近似解としてのメタン価が算出される。
以上のようにして得られた測定対象ガスのメタン価の値および熱量Ｑｎｓの値が表示機構４５に表示される。
なお、測定対象ガスおよび参照ガスは、ガス排出部１３を介して装置外部に排出される。
次いで、メタン価算出機構４０によって、熱量測定機構２０によって得られた基礎熱量Ｑ´の値と、上記の特定の関係式、例えばＡ＝０とした基準となるメタン価算出曲線を表す式（１）とから、ＡＶＬ値の近似解としてのメタン価が算出される。
以上のようにして得られた測定対象ガスのメタン価の値および熱量Ｑｎｓの値が表示機構４５に表示される。
なお、測定対象ガスおよび参照ガスは、ガス排出部１３を介して装置外部に排出される。
而して、上記のメタン価算出方法によれば、予め取得された基礎熱量Ｑ´とメタン価との特定の関係式を利用することにより、測定対象ガスの基礎熱量Ｑ´を測定するだけで、当該測定対象ガスのメタン価を求めることができる。特定の関係式は、各々ＡＶＬ基準に基づくメタン価の値（ＡＶＬ値）が互いに異なる天然ガスよりなる複数種の基準ガスについて、基礎熱量の値とＡＶＬ値との相関関係を実験による裏づけにより定量的に明らかにしたものであるため、得られるメタン価は一応の信頼性を有するものとなる。
従って、このようなメタン価算出方法が実行される上記構成のメタン価測定装置によれば、熱量測定機構２０によって測定対象ガスの基礎熱量Ｑ´を連続的に測定することにより、実際の状況に即した測定対象ガスのＡＶＬ値の近似解としてのメタン価を連続的に取得することができるので、燃料ガスとしての天然ガスの実際の燃料性状の監視を行うことができる。従って、ガス組成の変動が生じた場合には、ガス組成の変動に伴うメタン価の変動を速やかに検出することができる。
従って、このようなメタン価算出方法が実行される上記構成のメタン価測定装置によれば、熱量測定機構２０によって測定対象ガスの基礎熱量Ｑ´を連続的に測定することにより、実際の状況に即した測定対象ガスのＡＶＬ値の近似解としてのメタン価を連続的に取得することができるので、燃料ガスとしての天然ガスの実際の燃料性状の監視を行うことができる。従って、ガス組成の変動が生じた場合には、ガス組成の変動に伴うメタン価の変動を速やかに検出することができる。
また、上記のメタン価測定装置においては、熱量測定機構２０およびメタン価算出機構４０が防爆性容器１０内に配設されてなるものであることから、測定システムの構築および操作が簡便となる。しかも、測定に際しては、相当の時間を要することがなく、しかも基礎熱量Ｑ´の算出処理とメタン価の算出処理との間にタイムラグが生じることがないため、メタン価をリアルタイムに測定することができる。
さらにまた、熱量測定機構２０が、測定対象ガスの熱量を音速換算熱量Ｑｓおよび屈折率換算熱量Ｑｎの２つの換算熱量に基づいて算出する構成のものであることにより、得られる熱量Ｑｎｓは、測定対象ガスのガス組成によらずに当該測定対象ガスの熱量Ｑの真値との差の小さい値となるので、算出されるメタン価の値の信頼性が一層高いものなる。
上述したように、天然ガスを燃料ガスとして利用する場合には、実際には、例えば地域ごとに異なる基準に基づく算出方法によって算出されるメタン価が要請されるため、例えば、メタン価算出機構は、上記（ａ）ＡＶＬ基準、（ｂ）ＣＡＲＢ基準、（ｃ）ＧＲＩ（Ｌｃ）基準および（ｄ）ＧＲＩ（Ｈ／Ｃ）基準の各々に係るメタン価の近似解を算出する機能を有する構成とされていることが好ましい。このような構成のものにおいては、各基準に基づく方法により算出されるメタン価の値と、基礎熱量の値との特定の関係式を取得しておけばよい。
（ｂ）ＣＡＲＢ基準に係る特定の関係式は、上記式（２）で表されるものであることが好ましい。
上記式（２）におけるＢについての数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｂ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｂの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＣＡＲＢ基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。ＣＡＲＢ基準に係るメタン価算出曲線の一例を図３に示す。図３における実線で示す曲線が、基準となるメタン価算出曲線である。この基準となるメタン価算出曲線は、上記式（ｅ）および上記式（ｆ）で表される曲線部分の各々が変曲点を生じることなく連続するものであって、例えば実在し得るＬＮＧ気化ガスのガス組成を考慮して、設定されている。図３における破線で示す曲線は、式（２）においてＢ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（２）においてＢ＝−２．０としたメタン価算出曲線である。
上記式（２）におけるＢの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
上記式（２）におけるＢについての数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｂ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｂの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＣＡＲＢ基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。ＣＡＲＢ基準に係るメタン価算出曲線の一例を図３に示す。図３における実線で示す曲線が、基準となるメタン価算出曲線である。この基準となるメタン価算出曲線は、上記式（ｅ）および上記式（ｆ）で表される曲線部分の各々が変曲点を生じることなく連続するものであって、例えば実在し得るＬＮＧ気化ガスのガス組成を考慮して、設定されている。図３における破線で示す曲線は、式（２）においてＢ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（２）においてＢ＝−２．０としたメタン価算出曲線である。
上記式（２）におけるＢの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
（ｃ）ＧＲＩ（Ｌｃ）基準に係る特定の関係式は、上記式（３）で表されるものであることが好ましい。
上記式（３）におけるＣは、−２．０〜２．０の範囲から選択される値であり、この数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｃ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｃの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＧＲＩ（Ｌｃ）基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。ＧＲＩ（Ｌｃ）基準に係るメタン価算出曲線の一例を図４に示す。図４における実線で示す曲線が、式（３）においてＣ＝０とした基準となるメタン価算出曲線である。また、破線で示す曲線は、式（３）においてＣ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（３）においてＣ＝−２．０としたメタン価算出曲線である。
上記式（３）におけるＣの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
上記式（３）におけるＣは、−２．０〜２．０の範囲から選択される値であり、この数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｃ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｃの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＧＲＩ（Ｌｃ）基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。ＧＲＩ（Ｌｃ）基準に係るメタン価算出曲線の一例を図４に示す。図４における実線で示す曲線が、式（３）においてＣ＝０とした基準となるメタン価算出曲線である。また、破線で示す曲線は、式（３）においてＣ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（３）においてＣ＝−２．０としたメタン価算出曲線である。
上記式（３）におけるＣの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
（ｄ）ＧＲＩ（Ｈ／Ｃ）基準に係る特定の関係式は、上記式（４）で表されるものであることが好ましい。
上記式（４）におけるＤは、−２．０〜２．０の範囲から選択される値であり、この数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｄ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｄの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＧＲＩ（Ｈ／Ｃ）基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。ＧＲＩ（Ｈ／Ｃ）基準に係るメタン価算出曲線の一例を図５に示す。図５における実線で示す曲線が、式（４）においてＤ＝０とした基準となるメタン価算出曲線である。また、破線で示す曲線は、式（４）においてＤ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（４）においてＤ＝２．０としたメタン価算出曲線である。
上記式（４）におけるＤの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
上記式（４）におけるＤは、−２．０〜２．０の範囲から選択される値であり、この数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｄ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｄの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＧＲＩ（Ｈ／Ｃ）基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。ＧＲＩ（Ｈ／Ｃ）基準に係るメタン価算出曲線の一例を図５に示す。図５における実線で示す曲線が、式（４）においてＤ＝０とした基準となるメタン価算出曲線である。また、破線で示す曲線は、式（４）においてＤ＝２．０としたメタン価算出曲線であり、一点鎖線で示す曲線は、式（４）においてＤ＝２．０としたメタン価算出曲線である。
上記式（４）におけるＤの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
このような構成のメタン価測定装置においては、表示機構４５は、同一の測定対象ガスについて、上記複数の特定の関係式の各々と、熱量測定機構２０により測定される基礎熱量とから算出される各基準に係る複数のメタン価を同時に表示する機能を有するものであっても、選択された基準に係るメタン価を他のものと切り替え可能に表示する機能を有するものであってもいずれでもよい。このような構成とされていることにより、必要とされる基準に係るメタン価を得ることができるので、高い利便性が得られる。
以上においては、例えば窒素ガスを含まないＬＮＧ気化ガスである測定対象ガスのメタン価を算出する場合について説明したが、測定対象ガスが窒素ガスを含むものである場合には、実際上は問題にはならないが、窒素ガスが含有されていることに起因して測定誤差が生じてしまう。然るに、本発明者らが鋭意研究を重ねた結果、例えばＬＮＧ気化ガスに含まれる窒素ガス濃度と、当該窒素ガス濃度に起因するメタン価の変動量（誤差）との間に特定の相関関係が成立し、窒素ガス濃度に応じた補正量で補正することにより、各々の基準に対応するメタン価の近似解を一層高い信頼性をもって得ることができることを見出した。
以下、例えば窒素ガスを含むＬＮＧ気化ガスについて、（ａ）ＡＶＬ基準に基づく方法により得られるメタン価の値（ＡＶＬ値）の近似解を算出する場合を例に挙げて具体的に説明する。
メタン価の算出に際して、測定対象ガスに含まれる窒素ガス濃度に基づいた補正を行う場合には、メタン価算出機構４０は、熱量測定機構２０により測定された当該測定対象ガスの基礎熱量Ｑ´の値および窒素測定計算機構３５によって得られた窒素ガス濃度ＸN2〔ｖｏｌ％〕の値と、ＡＶＬ基準に基づく方法により得られるメタン価の値（ＡＶＬ値）と基礎熱量Ｑ´の値との特定の関係式とから、測定対象ガスのメタン価を算出する機能を有するものとされる。
窒素ガスを含むＬＮＧ気化ガスについてのＡＶＬ基準に係る特定の関係式は、上記式（５）で表されるものであることが好ましい。この特定の関係式は、上記のＡＶＬ基準に係る式（１）で表される特定の関係式と同様の方法により取得されたものである。
上記式（５）において、ＭＮはメタン価、具体的にはＡＶＬ値の近似解、ｆ(Q´) は、測定対象ガスの基礎熱量Ｑ´の値に応じて選択される上記式（ｇ）〜上記式（ｊ）で表されるいずれかの関数であり、Ｅは−２．０〜２．０の範囲から選択される値である。
上記式（ｇ）〜上記式（ｊ）における０．３２０ＸN2 の項は、窒素ガス濃度に基づくメタン価の補正量を示す。メタン価の補正量は、例えば図６に示すように、横軸を測定対象ガスに含まれる窒素ガスの、体積百分率で示される濃度〔ｖｏｌ％〕とし、縦軸をメタン価とする座標系において、メタン価の値が互いに異なる複数種の基準ガスの各々について、窒素ガスの濃度値ＸN2〔ｖｏｌ％〕とＡＶＬ値との関係を示す実測値を取得し、得られた実測値を例えば線形近似することにより取得された近似直線に基づいて設定されたものである。図６から明らかなように、各々の基準ガスについての近似曲線は、互いに同一の大きさの傾きを有していることが理解される。
上記式（５）におけるＥについて設定される数値範囲は、式（５）においてＥ＝０とした式で表される基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｅの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＡＶＬ値に対する誤差が５．０％以内となり、高い信頼度が得られる。
上記式（５）におけるＥの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
上記式（５）において、ＭＮはメタン価、具体的にはＡＶＬ値の近似解、ｆ(Q´) は、測定対象ガスの基礎熱量Ｑ´の値に応じて選択される上記式（ｇ）〜上記式（ｊ）で表されるいずれかの関数であり、Ｅは−２．０〜２．０の範囲から選択される値である。
上記式（ｇ）〜上記式（ｊ）における０．３２０ＸN2 の項は、窒素ガス濃度に基づくメタン価の補正量を示す。メタン価の補正量は、例えば図６に示すように、横軸を測定対象ガスに含まれる窒素ガスの、体積百分率で示される濃度〔ｖｏｌ％〕とし、縦軸をメタン価とする座標系において、メタン価の値が互いに異なる複数種の基準ガスの各々について、窒素ガスの濃度値ＸN2〔ｖｏｌ％〕とＡＶＬ値との関係を示す実測値を取得し、得られた実測値を例えば線形近似することにより取得された近似直線に基づいて設定されたものである。図６から明らかなように、各々の基準ガスについての近似曲線は、互いに同一の大きさの傾きを有していることが理解される。
上記式（５）におけるＥについて設定される数値範囲は、式（５）においてＥ＝０とした式で表される基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｅの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＡＶＬ値に対する誤差が５．０％以内となり、高い信頼度が得られる。
上記式（５）におけるＥの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
このようなメタン価測定装置においては、音速換算熱量測定機構２５において求められた音速換算熱量Ｑｓの値と、屈折率換算熱量測定機構２１において求められた屈折率換算熱量Ｑｎの値とに基づいて、上記式（７）および上記式（８）により、補正因子αとして特定の範囲内において選択された値を用いて、窒素ガスの濃度ＸN2〔ｖｏｌ％〕および熱量Ｑｎｓが算出される。そして、このようにして得られた熱量Ｑｎｓの値と、窒素ガス濃度ＸN2〔ｖｏｌ％〕の値とに基づいて、上記式（９）によって、測定対象ガスの基礎熱量Ｑ´が算出される。
次いで、メタン価算出機構４０によって、熱量測定機構２０によって得られた基礎熱量Ｑ´の値および窒素ガスの濃度値ＸN2〔ｖｏｌ％〕と、上記式（５）において例えばＥ＝０とした特定の関係式とから、ＡＶＬ値の近似解としてのメタン価が算出される。
次いで、メタン価算出機構４０によって、熱量測定機構２０によって得られた基礎熱量Ｑ´の値および窒素ガスの濃度値ＸN2〔ｖｏｌ％〕と、上記式（５）において例えばＥ＝０とした特定の関係式とから、ＡＶＬ値の近似解としてのメタン価が算出される。
而して、上記のメタン価算出方法によれば、実験による裏づけによって定量的に明らかにした、窒素ガス濃度ＸN2 とメタン価との相関関係を利用することにより、測定対象ガスに含まれる窒素ガスによって生ずる誤差が補償されるので、当該測定対象ガスのメタン価を一層高い信頼性をもって求めることができる。
従って、このようなメタン価算出方法が実行される上記構成のメタン価測定装置によれば、実際の状況に即した測定対象ガスのＡＶＬ値の近似解としてのメタン価を、一層高い信頼性をもって連続的に取得することができるので、燃料ガスとしての天然ガスの実際の燃料性状の監視を一層確実に行うことができる。
従って、このようなメタン価算出方法が実行される上記構成のメタン価測定装置によれば、実際の状況に即した測定対象ガスのＡＶＬ値の近似解としてのメタン価を、一層高い信頼性をもって連続的に取得することができるので、燃料ガスとしての天然ガスの実際の燃料性状の監視を一層確実に行うことができる。
また、窒素ガスを含む例えばＬＮＧ気化ガスについての、（ｃ）ＧＲＩ（Ｌｃ）基準に係る特定の関係式は、上記式（６）で表されるものであることが好ましい。
上記式（６）におけるＦは、−２．０〜２．０の範囲から選択される値であり、この数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｆ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｆの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＧＲＩ（Ｌｃ）基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。
上記式（６）におけるＦの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
上記式（６）におけるＦは、−２．０〜２．０の範囲から選択される値であり、この数値範囲は、上記式（１）におけるＡについての数値範囲と同様に、Ｆ＝０とした基準となるメタン価算出曲線自体の、実際のＬＮＧ気化ガスの燃料性状に即した補正が行われる実用的な許容範囲を示す。Ｆの値が上記数値範囲内であれば、後述する実験例の結果に示されるように、算出されるメタン価（近似解）の、ＧＲＩ（Ｌｃ）基準のメタン価に対する誤差率が５．０％以内となり、高い信頼度が得られる。
上記式（６）におけるＦの値は、例えば、上記式（１）におけるＡの値の設定方法と同様の方法により設定することができる。
（ｂ）ＣＡＲＢ基準に係るメタン価および（ｄ）ＧＲＩ（Ｈ／Ｃ）基準に係るメタン価の算出に際しては、測定対象ガスに含まれる窒素ガス濃度に拘わらず、上記式（２）で表される特定の関係式および上記式（４）で表される特定の関係式が用いられる。
以上、本発明の実施の形態について説明したが、本発明は上記の実施の形態に限定されるものではなく、種々の変更を加えることができる。
本発明においては、熱量測定機構は、上記構成のものに限定されるものではなく、熱伝導率換算熱量の値と、屈折率換算熱量の値とに基づいて熱量の値を求める構成を有する装置が用いられていてもよい。また、熱量と特定の対応関係を有する物性値の１つ、例えば屈折率、熱伝導率、音速のうちから選ばれる１つを測定し、その測定値に基づいて熱量を求める構成のものであってもよい。このようにして得られた測定対象ガスの熱量に基づいて算出される基礎熱量の値を利用しても、各基準に基づくメタン価の値に対して所定の許容範囲内のメタン価を得ることができる。
本発明においては、熱量測定機構は、上記構成のものに限定されるものではなく、熱伝導率換算熱量の値と、屈折率換算熱量の値とに基づいて熱量の値を求める構成を有する装置が用いられていてもよい。また、熱量と特定の対応関係を有する物性値の１つ、例えば屈折率、熱伝導率、音速のうちから選ばれる１つを測定し、その測定値に基づいて熱量を求める構成のものであってもよい。このようにして得られた測定対象ガスの熱量に基づいて算出される基礎熱量の値を利用しても、各基準に基づくメタン価の値に対して所定の許容範囲内のメタン価を得ることができる。
以下、本発明の実験例について説明する。
〔実験例１〕
下記表１に示すガス組成を有する１１種類の試料ガスを用意し、図１に示す構成のメタン価測定装置を用いて、各々の試料ガスＡ〜Ｋのメタン価を測定した。メタン価の算出は、上記式（１）においてＡ＝０とした特定の関係式を利用した。そして、試料ガスについて測定される基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ａ）を選択し、基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕より大きく、５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｂ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きく、６３．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｃ）を選択し、基礎熱量Ｑ´の値が６３．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｄ）を選択した。また、上記式（１）においてＡ＝−２．０とした特定の関係式、およびＡ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表２に示す。
下記表１に示すガス組成を有する１１種類の試料ガスを用意し、図１に示す構成のメタン価測定装置を用いて、各々の試料ガスＡ〜Ｋのメタン価を測定した。メタン価の算出は、上記式（１）においてＡ＝０とした特定の関係式を利用した。そして、試料ガスについて測定される基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ａ）を選択し、基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕より大きく、５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｂ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きく、６３．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｃ）を選択し、基礎熱量Ｑ´の値が６３．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｄ）を選択した。また、上記式（１）においてＡ＝−２．０とした特定の関係式、およびＡ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表２に示す。
各々の試料ガスＡ〜Ｋについて、ＡＶＬ基準に基づく方法（ＡＶＬ社製のメタン価算出ソフト）によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率を算出した。誤差率は、〔（誤差の絶対値）／真値〕×１００〔％〕で得られる値である。上記式（１）におけるＡの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表２に示す。
〔実験例２〕
上記式（２）においてＢ＝０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。メタン価の算出においては、試料ガスについて測定される基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｅ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｆ）を選択した。また、上記式（２）においてＢ＝−２．０とした特定の関係式、および、上記式（２）においてＢ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表３に示す。
各々の試料ガスＡ〜Ｋについて、ＣＡＲＢ基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（２）におけるＢの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表３に示す。
上記式（２）においてＢ＝０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。メタン価の算出においては、試料ガスについて測定される基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｅ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｆ）を選択した。また、上記式（２）においてＢ＝−２．０とした特定の関係式、および、上記式（２）においてＢ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表３に示す。
各々の試料ガスＡ〜Ｋについて、ＣＡＲＢ基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（２）におけるＢの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表３に示す。
〔実験例３〕
上記式（３）においてＣ＝０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。また、上記式（３）においてＣ＝−２．０とした特定の関係式、および、上記式（３）においてＣ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表４に示す。
各々の試料ガスＡ〜Ｋについて、ＧＲＩ（Ｌｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（３）におけるＣの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表４に示す。
上記式（３）においてＣ＝０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。また、上記式（３）においてＣ＝−２．０とした特定の関係式、および、上記式（３）においてＣ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表４に示す。
各々の試料ガスＡ〜Ｋについて、ＧＲＩ（Ｌｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（３）におけるＣの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表４に示す。
〔実験例４〕
上記式（４）においてＤ＝０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。また、上記式（４）においてＤ＝−２．０とした特定の関係式、および、上記式（４）においてＤ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表５に示す。
各々の試料ガスＡ〜Ｋについて、ＧＲＩ（Ｈ／Ｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（４）におけるＤの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表５に示す。
上記式（４）においてＤ＝０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。また、上記式（４）においてＤ＝−２．０とした特定の関係式、および、上記式（４）においてＤ＝２．０とした特定の関係式を利用したことの他は、実験例１と同様の方法により、各々の試料ガスＡ〜Ｋのメタン価を算出した。結果を下記表５に示す。
各々の試料ガスＡ〜Ｋについて、ＧＲＩ（Ｈ／Ｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（４）におけるＤの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表５に示す。
〔実験例５〕
下記表６に示すガス組成を有する１１種類の試料ガスを用意し、図１に示す構成のメタン価測定装置を用いて、各々の試料ガスａ〜ｋのメタン価を測定した。メタン価の算出は、上記式（５）においてＥ＝０とした特定の関係式を利用した。そして、試料ガスについて測定される基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｇ）を選択し、基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕より大きく、５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｈ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きく、６３．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｉ）を選択し、基礎熱量Ｑ´の値が６３．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｊ）を選択した。また、上記式（５）においてＥ＝−２．０とした特定の関係式、およびＥ＝２．０とした特定の関係式を利用したことの他は、同様の方法により、各々の各々の試料ガスａ〜ｋのメタン価を測定した。結果を下記表７に示す。
下記表６に示すガス組成を有する１１種類の試料ガスを用意し、図１に示す構成のメタン価測定装置を用いて、各々の試料ガスａ〜ｋのメタン価を測定した。メタン価の算出は、上記式（５）においてＥ＝０とした特定の関係式を利用した。そして、試料ガスについて測定される基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｇ）を選択し、基礎熱量Ｑ´の値が４２．０〔ＭＪ／ｍ3 〕より大きく、５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｈ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きく、６３．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｉ）を選択し、基礎熱量Ｑ´の値が６３．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｊ）を選択した。また、上記式（５）においてＥ＝−２．０とした特定の関係式、およびＥ＝２．０とした特定の関係式を利用したことの他は、同様の方法により、各々の各々の試料ガスａ〜ｋのメタン価を測定した。結果を下記表７に示す。
試料ガスａ〜ｋの各々について、ＡＶＬ基準に基づく方法（ＡＶＬ社製のメタン価算出ソフト）によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。誤差率は、〔（誤差の絶対値）／真値〕×１００〔％〕で得られる値である。上記式（５）におけるＥの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表７に示す。
〔実験例６〕
上記式（２）においてＢ＝０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。メタン価の算出においては、試料ガスについて測定される基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｅ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｆ）を選択した。また、上記式（２）においてＢ＝−２．０とした特定の関係式、および、上記式（２）においてＢ＝２．０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。結果を下記表８に示す。
試料ガスａ〜ｋの各々について、ＣＡＲＢ基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（２）におけるＢの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表８に示す。
上記式（２）においてＢ＝０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。メタン価の算出においては、試料ガスについて測定される基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕以下である場合には、上記式（ｅ）を選択し、基礎熱量Ｑ´の値が５５．０〔ＭＪ／ｍ3 〕より大きい場合には、上記式（ｆ）を選択した。また、上記式（２）においてＢ＝−２．０とした特定の関係式、および、上記式（２）においてＢ＝２．０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。結果を下記表８に示す。
試料ガスａ〜ｋの各々について、ＣＡＲＢ基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（２）におけるＢの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表８に示す。
〔実験例７〕
上記式（６）においてＦ＝０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。また、上記式（６）においてＦ＝−２．０とした特定の関係式、および、上記式（６）においてＦ＝２．０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。結果を下記表９に示す。
試料ガスａ〜ｋの各々について、ＧＲＩ（Ｌｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（６）におけるＦの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表９に示す。
上記式（６）においてＦ＝０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。また、上記式（６）においてＦ＝−２．０とした特定の関係式、および、上記式（６）においてＦ＝２．０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。結果を下記表９に示す。
試料ガスａ〜ｋの各々について、ＧＲＩ（Ｌｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（６）におけるＦの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表９に示す。
〔実験例８〕
上記式（４）においてＤ＝０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。また、上記式（４）においてＤ＝−２．０とした特定の関係式、および、上記式（４）においてＤ＝２．０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。結果を下記表１０に示す。
試料ガスａ〜ｋの各々について、ＧＲＩ（Ｈ／Ｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（４）におけるＤの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表１０に示す。
上記式（４）においてＤ＝０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。また、上記式（４）においてＤ＝−２．０とした特定の関係式、および、上記式（４）においてＤ＝２．０とした特定の関係式を利用したことの他は、実験例５と同様の方法により、各々の試料ガスａ〜ｋのメタン価を算出した。結果を下記表１０に示す。
試料ガスａ〜ｋの各々について、ＧＲＩ（Ｈ／Ｃ）基準に準拠した方法によって、ガス組成に基づいて算出されるメタン価を真値としたときの、基礎熱量に基づいて算出されたメタン価の値（熱量換算値）の当該真値に対する誤差率〔％〕を算出した。上記式（４）におけるＤの値を−２．０〜２．０の範囲内で設定したときに算出されるメタン価の値（熱量換算値）の真値に対する誤差率の最大値を下記表１０に示す。
以上の結果から、本発明に係るメタン価算出方法によれば、ＡＶＬ基準、ＣＡＲＢ基準、ＧＲＩ（Ｌｃ）基準およびＧＲＩ（Ｈ／Ｃ）基準のいずれの基準であっても、試料ガスの組成に拘わらず、これらの基準に係るメタン価に対して一定の誤差範囲内の値のメタン価（近似解）を算出することができることが確認された。ここに、誤差率が５．０％以内であれば、実用上問題ない程度の誤差であるといえる。
本発明は、燃料ガスとしての天然ガスのガス組成の変動に伴うメタン価の変動や熱量の変動といった燃料性状の変動をリアルタイムに検出することができるものであるので、ＬＮＧ燃料エンジンの燃焼制御を行うに際して、極めて有用なものとなることが期待される。
１０ 防爆性容器
１１ 測定対象ガス導入部
１２ 参照ガス導入部
１３ ガス排出部
２０ 熱量測定機構
２１ 屈折率換算熱量測定機構
２２ 屈折率測定手段
２３ 屈折率−熱量換算処理手段
２５ 音速換算熱量測定機構
２６ 音速測定手段
２７ 音速−熱量換算処理手段
３０ 窒素濃度測定機構
３５ 熱量計算機構
４０ メタン価算出機構
４５ 表示機構
１１ 測定対象ガス導入部
１２ 参照ガス導入部
１３ ガス排出部
２０ 熱量測定機構
２１ 屈折率換算熱量測定機構
２２ 屈折率測定手段
２３ 屈折率−熱量換算処理手段
２５ 音速換算熱量測定機構
２６ 音速測定手段
２７ 音速−熱量換算処理手段
３０ 窒素濃度測定機構
３５ 熱量計算機構
４０ メタン価算出機構
４５ 表示機構
Claims (20)
- 各々天然ガスよりなる互いにメタン価の値が異なる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式を予め取得しておき、
測定対象ガスである天然ガスの基礎熱量を測定し、
測定された当該測定対象ガスの基礎熱量の値と、前記特定の関係式とから当該測定対象ガスのメタン価を算出することを特徴とするメタン価算出方法。 - 前記測定対象ガスである天然ガスは、液化天然ガスを気化させて得られたものであることを特徴とする請求項１に記載のメタン価算出方法。
- 前記特定の関係式が下記式（１）で表されるものであることを特徴とする請求項２に記載のメタン価算出方法。
- 前記特定の関係式が下記式（２）で表されるものであり、カリフォルニア州大気資源評議会で規定された特定の演算式により算出されるメタン価の近似値を取得することを特徴とする請求項２に記載のメタン価算出方法。
- 前記特定の関係式が下記式（３）で表されるものであり、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得することを特徴とする請求項２に記載のメタン価算出方法。
- 前記特定の関係式が下記式（４）で表されるものであり、ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出されるメタン価の近似値を取得することを特徴とする請求項２に記載のメタン価算出方法。
- 各々窒素ガスを含有する天然ガスよりなる互いにメタン価の値が異なる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式を予め取得しておき、
測定対象ガスである窒素ガスを含有する天然ガスの基礎熱量および当該測定対象ガスに含まれる窒素ガス濃度を測定し、
測定された当該測定対象ガスの基礎熱量の値および窒素ガスの濃度値と、前記特定の関係式とから当該測定対象ガスのメタン価を算出することを特徴とするメタン価算出方法。 - 前記測定対象ガスである天然ガスは、液化天然ガスを気化させて得られたものであることを特徴とする請求項７に記載のメタン価算出方法。
- 前記特定の関係式が下記式（５）で表されるものであることを特徴とする請求項８に記載のメタン価算出方法。
- 前記特定の関係式が下記式（６）で表されるものであり、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得することを特徴とする請求項８に記載のメタン価算出方法。
- 前記測定対象ガスの基礎熱量は、当該測定対象ガスの屈折率から得られる屈折率換算熱量と、当該測定対象ガスの音速から得られる音速換算熱量とに基づいて得られるものであることを特徴とする請求項１乃至請求項１０のいずれかに記載のメタン価算出方法。
- 測定対象ガスである天然ガスの基礎熱量を測定する熱量測定機構と、
予め取得しておいた、各々メタン価の値が異なる天然ガスよりなる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式と、前記熱量測定機構により測定された当該測定対象ガスの基礎熱量の値とから、当該測定対象ガスのメタン価を算出するメタン価算出機構と
を備えることを特徴とするメタン価測定装置。 - 前記測定対象ガスである天然ガスは、液化天然ガスを気化させて得られたものであることを特徴とする請求項１２に記載のメタン価測定装置。
- 前記特定の関係式が下記式（１）で表されるものであることを特徴とする請求項１３に記載のメタン価測定装置。
- カリフォルニア州大気資源評議会で規定された特定の演算式により算出されるメタン価の近似値を取得する下記式（２）で表される特定の関係式、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得する下記式（３）で表される特定の関係式、および、ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出されるメタン価の近似値を取得する下記式（４）で表される特定の関係式がさらに設定されていることを特徴とする請求項１４に記載のメタン価測定装置。
- 測定対象ガスである窒素ガスを含有する天然ガスの基礎熱量を測定する熱量測定機構と、
当該測定対象ガスに含まれる窒素ガス濃度を測定する濃度測定機構と、
予め取得しておいた、各々メタン価の値が異なる窒素ガスを含有する天然ガスよりなる複数種の基準ガスについてのメタン価と基礎熱量との特定の関係式と、前記熱量測定機構により測定された当該測定対象ガスの基礎熱量の値および前記濃度測定機構により測定された窒素ガス濃度の値とから、当該測定対象ガスのメタン価を算出するメタン価算出機構と
を備えることを特徴とするメタン価測定装置。 - 前記測定対象ガスである天然ガスは、液化天然ガスを気化させて得られたものであることを特徴とする請求項１６に記載のメタン価測定装置。
- 前記特定の関係式が下記式（５）で表されるものであることを特徴とする請求項１７に記載のメタン価測定装置。
- カリフォルニア州大気資源評議会で規定された特定の演算式により算出されるメタン価の近似値を取得する下記式（２）で表される特定の関係式、ＩＳＯ／ＴＲ ２２３０２ ３．１．１に準拠した方法により算出されるメタン価の近似値を取得する下記式（６）で表される特定の関係式、および、ＩＳＯ／ＴＲ ２２３０２ ３．１．２に準拠した方法により算出されるメタン価の近似値を取得する下記式（４）で表される特定の関係式がさらに設定されていることを特徴とする請求項１８に記載のメタン価測定装置。
- 前記熱量測定機構が、測定対象ガスの屈折率の値から屈折率換算熱量を求める屈折率換算熱量測定手段と、当該測定対象ガスの音速の値から音速換算熱量を求める音速換算熱量測定手段と、当該測定対象ガスの基礎熱量を屈折率換算熱量と音速換算熱量とに基づいて算出する熱量算出手段とを備えることを特徴とする請求項１２乃至請求項１９のいずれか一項に記載のメタン価測定装置。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2014264811 | 2014-12-26 | ||
| JP2014264811 | 2014-12-26 | ||
| PCT/JP2015/085177 WO2016104270A1 (ja) | 2014-12-26 | 2015-12-16 | メタン価算出方法およびメタン価測定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2016104270A1 true JPWO2016104270A1 (ja) | 2017-10-05 |
| JP6702885B2 JP6702885B2 (ja) | 2020-06-03 |
Family
ID=56150295
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016566149A Active JP6702885B2 (ja) | 2014-12-26 | 2015-12-16 | メタン価算出方法およびメタン価測定装置 |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP6702885B2 (ja) |
| WO (1) | WO2016104270A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115790908A (zh) * | 2023-02-08 | 2023-03-14 | 成都千嘉科技股份有限公司 | 基于热量计量的天然气计量方法和装置 |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR3060796B1 (fr) * | 2016-12-15 | 2019-07-05 | Engie | Procede pour calculer en temps reel l'indice de methane mn d'un gaz naturel liquefie |
| US10976295B2 (en) * | 2017-06-26 | 2021-04-13 | Mustang Sampling Llc | System and methods for methane number generation |
| JP6959078B2 (ja) * | 2017-09-07 | 2021-11-02 | 理研計器株式会社 | ガス分析方法およびガス分析装置 |
| WO2019187710A1 (ja) * | 2018-03-28 | 2019-10-03 | 理研計器株式会社 | ガス検出方法およびガス検出装置 |
| WO2019186794A1 (ja) * | 2018-03-28 | 2019-10-03 | 理研計器株式会社 | ガス検出方法およびガス検出装置 |
| JP7445531B2 (ja) | 2020-06-15 | 2024-03-07 | 理研計器株式会社 | メタン価測定装置およびメタン価算出方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002509249A (ja) * | 1998-01-16 | 2002-03-26 | ビージー インテレクチュアル プロパティー リミテッド | ガス熱量測定方法及び装置 |
| JP2004162649A (ja) * | 2002-11-14 | 2004-06-10 | Tokyo Gas Co Ltd | 燃料ガスの燃焼性検知方法及び装置及びガス燃料内燃機関 |
| WO2012066828A1 (ja) * | 2010-11-15 | 2012-05-24 | 理研計器株式会社 | 燃焼性ガスの比重測定方法および比重測定装置並びにウォッベ指数測定装置 |
| WO2014054081A1 (ja) * | 2012-10-05 | 2014-04-10 | 川崎重工業株式会社 | ガスエンジン用の燃焼安定化装置 |
-
2015
- 2015-12-16 JP JP2016566149A patent/JP6702885B2/ja active Active
- 2015-12-16 WO PCT/JP2015/085177 patent/WO2016104270A1/ja active Application Filing
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002509249A (ja) * | 1998-01-16 | 2002-03-26 | ビージー インテレクチュアル プロパティー リミテッド | ガス熱量測定方法及び装置 |
| JP2004162649A (ja) * | 2002-11-14 | 2004-06-10 | Tokyo Gas Co Ltd | 燃料ガスの燃焼性検知方法及び装置及びガス燃料内燃機関 |
| WO2012066828A1 (ja) * | 2010-11-15 | 2012-05-24 | 理研計器株式会社 | 燃焼性ガスの比重測定方法および比重測定装置並びにウォッベ指数測定装置 |
| WO2014054081A1 (ja) * | 2012-10-05 | 2014-04-10 | 川崎重工業株式会社 | ガスエンジン用の燃焼安定化装置 |
Non-Patent Citations (1)
| Title |
|---|
| SCHLEY P,FRIELING H: "Bestimmung der Methanzahl aus Gasbeschaffenheitskenngroessen", GWF GAS ERDGAS, vol. 141, no. 1, JPN6016000764, 2000, pages 28 - 33, XP009506258, ISSN: 0004126526 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115790908A (zh) * | 2023-02-08 | 2023-03-14 | 成都千嘉科技股份有限公司 | 基于热量计量的天然气计量方法和装置 |
| CN115790908B (zh) * | 2023-02-08 | 2023-06-23 | 成都千嘉科技股份有限公司 | 基于热量计量的天然气计量方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2016104270A1 (ja) | 2016-06-30 |
| JP6702885B2 (ja) | 2020-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6759206B2 (ja) | メタン価算出方法およびメタン価測定装置 | |
| WO2016104270A1 (ja) | メタン価算出方法およびメタン価測定装置 | |
| JP6402387B2 (ja) | 熱量測定装置および熱量測定方法 | |
| JP5782454B2 (ja) | 燃焼性ガスの比重測定方法および比重測定装置並びにウォッベ指数測定装置 | |
| CN108195591A (zh) | 一种进气质量流量的修正方法及*** | |
| Koczak et al. | Particulate emissions in GDI vehicle transients: An examination of FTP, HWFET, and US06 measurements | |
| JP2000338093A (ja) | メタンガス中の窒素ガス含有率測定方法および装置 | |
| Bai-Gang et al. | A new equation of state for hydrogen gas | |
| Guido et al. | Assessment of closed-loop combustion control capability for biodiesel blending detection and combustion impact mitigation for an Euro5 automotive diesel engine | |
| JP5308842B2 (ja) | 熱量測定方法および熱量測定装置 | |
| CN110431415A (zh) | 确定含烃的燃烧气体混合物的甲烷指数的方法 | |
| JP6858859B2 (ja) | 液化天然ガスの液相中のメタン価ｍｎをリアルタイムで計算する方法 | |
| Dakkach et al. | Liquid densities and excess volumes of biofuel mixtures:(2-butanol+ di-isopropyl ether) system at pressures up to 140 MPa and temperatures from 293.15 K to 393.28 K | |
| Zannikou et al. | Preliminary investigation of experimental conditions and precision of an alternative method to determine high boiling point components in motor gasoline | |
| JP2021196263A (ja) | メタン価測定装置およびメタン価算出方法 | |
| JP2006047071A (ja) | 液化天然ガスの熱量計測方法及び装置 | |
| JP5184983B2 (ja) | 熱量測定方法および熱量測定装置 | |
| JP2021004891A (ja) | メタン濃度算出方法およびメタン濃度測定装置 | |
| JP5844556B2 (ja) | 熱伝導式ガスセンサの出力補正方法およびガス検知器 | |
| JP2012251813A (ja) | ガスの燃焼熱量測定方法およびウォッベ指数測定方法並びにガス測定装置 | |
| Alaoui et al. | Excess enthalpies of ternary mixtures of (oxygenated additives+ aromatic hydrocarbon) mixtures in fuels and bio-fuels:(Dibutyl-ether+ 1-propanol+ benzene), or toluene, at T=(298.15 and 313.15) K | |
| JP2004162649A (ja) | 燃料ガスの燃焼性検知方法及び装置及びガス燃料内燃機関 | |
| JP4504320B2 (ja) | 熱量測定方法 | |
| KR20180123438A (ko) | 가스 엔진 작동 방법 및 그 제어 장치 | |
| Lozano-Martín et al. | Thermodynamic (p, ρ, T) characterization of a reference high-calorific natural gas mixture when hydrogen is added up to 20%(mol/mol) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180913 |
|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20180913 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20191008 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191120 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200428 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200507 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6702885 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |