(57)【特許請求の範囲】
【請求項１】
天然ガスを含むタンク（１）に含まれる液化天然ガスのメタン価ＭＮをリアルタイムで計算するためにコンピュータによって実装される方法であって、
前記天然ガスは、
−液体状態（１）の天然ガスの層（１１）であって、所与の時点ｔにおいて、その温度Ｔ（ｔ）およびその密度ρ（ｔ）によって定義され、気体状態（ｇ）の天然ガスの層（１２）とバランスがとれている、液体状態（１）の天然ガスの層（１１）と、
−気体状態（ｇ）の天然ガスの層（１２）であって、所与の時点ｔにおいて、その圧力Ｐ（ｔ）によって定義される、気体状態（ｇ）の天然ガスの層（１２）と
に分配され、
前記方法は、所与の時点ｔで、
Ａ．前記液体状態（１）の前記天然ガスの層（１１）の温度Ｔ０（ｔ）および密度ρ０（ｔ）、ならびに前記気体状態（ｇ）の前記天然ガスの層（１２）の圧力Ｐｏ（ｔ）を測定により決定するステップと、
Ｂ．応力最小値計算アルゴリズムによって、前記タンク（１）に含まれる液化天然ガス（１１）の組成を近似計算するステップであって、前記アルゴリズムが、以下のサブステップ、
ｂ１）応力最小値を計算することによって、圧力Ｐｏ（ｔ）で密度ρ０（ｔ）の第１の組成、または圧力Ｐｏ（ｔ）で同じ密度ρ０（ｔ）を有する組成の第１の包絡線を決定するステップと、
ｂ２）応力最小値を計算することによって、温度Ｔ０（ｔ）および圧力Ｐｏ（ｔ）で同じ密度ρ０（ｔ）を有し、かつ圧力Ｐｏ（ｔ）で平衡温度Ｔｅｑ（ｔ）がＴ０（ｔ）に等しい組成物の第２の包絡線を決定するステップと、
ｂ３）ＬＮＧの実際の組成に近似する組成を提供する前記第２の包絡線の非特異点を決定するステップと、
ｂ４）前記近似の組成から前記メタン価ＭＮを計算するステップと
を含む、ステップと
を含むアルゴリズムからなることを特徴とする、方法。
【請求項２】
組成物の前記第１の包絡線を決定するためのステップｂ１）は、以下の方程式を解くことによって、以下のように実行され、
0006858859.tif 000012
前記方程式において、
ρはｘおよびＰｏから計算される密度を意味し、
ρ０は測定された密度を意味し、
Ｐｏは測定された圧力を意味し、
ｘはｘｉからなる合成ベクトルを意味し、
ｘｉは成分ｉのモル分率を意味し、
ｌおよびｕは、それぞれモル分率ｘｉの下限および上限であることを示す添字である、請求項１に記載の方法。
【請求項３】
組成物の前記第１の包絡線を決定するためのステップｂ１）は、以下の方程式を解くことによって、以下のように実行され、
0006858859.tif 000013
前記方程式において、
ρはｘおよびＰｏから計算される密度を意味し、
ρ０は測定された密度を意味し、
Ｐｏは測定された圧力を意味し、
ｘはｘｉからなる合成ベクトルを意味し、
ｘｉは成分ｉのモル分率を意味し、
ｌおよびｕは、それぞれモル分率ｘｉの下限および上限であることを示す添字である、請求項１に記載の方法。
【請求項４】
組成物の前記第２の包絡線を決定するためのステップｂ２）は、以下の方程式を解くことによって、以下のように実行され、
0006858859.tif 000014
前記方程式において、
Ｔ０は測定された温度を意味し、
Ｔｅｑは圧力ＰｏでＬＮＧの平衡温度を意味し、
ρはｘおよびＰｏから計算される密度を意味し、
ρ０は測定された密度を意味し、
Ｐｏは測定された圧力を意味し、
ｘはｘｉからなる合成ベクトルを意味し、
ｘｉは成分ｉのモル分率を意味し、
ｌおよびｕは、それぞれモル分率ｘｉの下限および上限であることを示す添字である、請求項１〜３のいずれか一項に記載の方法。
【請求項５】
組成物の前記第２の包絡線を決定するためのステップｂ２）は、以下の方程式を解くことによって、以下のように実行され、
0006858859.tif 000015
前記方程式において、
Ｔ０は測定された温度を意味し、
Ｔｅｑは圧力ＰｏでＬＮＧの平衡温度を意味し、
ρはｘおよびＰｏから計算される密度を意味し、
ρ０は測定された密度を意味し、
Ｐｏは測定された圧力を意味し、
ｘはｘｉからなる合成ベクトルを意味し、
ｘｉは成分ｉのモル分率を意味し、
ｌおよびｕは、それぞれモル分率ｘｉの下限および上限であることを示す添字である、請求項１〜３のいずれか一項に記載の方法。
【請求項６】
前記第２の包絡線の前記非特異点が、前記第２の包絡線の重心、または前記第２の包絡線の前記重心に最も近い前記第２の包絡線の点である、請求項１〜５のいずれか一項に記載の方法。
【請求項７】
前記第２の包絡線の前記非特異点が、前記第２の包絡線の標準化組成物の再標準化平均として計算される、請求項１〜５のいずれか一項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般に、タンク内、特に船上タンク内に含まれる液化天然ガスのメタン価をリアルタイムで計算する方法に関する。
【背景技術】
【０００２】
燃料ＬＮＧは、ＣＯ２排出量、汚染粒子、エネルギー密度のすべての観点から、従来の燃料に代わる簡易で効率的な代替燃料である。道路輸送、海上輸送、鉄道輸送のように、ますます多くの関係者がその利用に目を向け始めている。
【０００３】
しかしながら、従来の精製された燃料とは異なり、ＬＮＧは標準的な組成を有しない。加えて、堆積物の所在地によっては、ＬＮＧは少なくとも重化合物（Ｃ２、Ｃ３、Ｃ４など）を含み得る。
【０００４】
この組成の違いは難点を伴う。すなわち、ノッキングの現象は、重化合物の割合が高い場合にエンジンで発生する可能性があり、不適切な場所で不適切な時間に重炭化水素によって引き起こされる不適切な燃焼によるものであって、燃焼による第２の面、したがって、性能、効率および機械的強度の問題の発生、さらにはエンジン停止などにつながる。メタン価は、この現象が出現し、必要な修正措置をとることのきっかけとなり得るパラメータである。
【０００５】
ノッキング現象の出現というこの問題を解決するために、メタン価はできるだけ安定していなければならない。しかし、メタン価は今日ではほとんど知られていない大きさであり、その決定方法は結果に関して非常に多様である。したがって、端末で受け取ったＬＮＧ組成からメタン価を計算することからなる現在使用されている方法は、依然として理論的見解であり、加えて、このメタン価が使用されるタンクを燃料補給するための方法および選択された計算規則に従って経時変化することが分かっている。したがって、メタン価を決定可能にすることは、小売ＬＮＧ燃料市場の発展を民主化し、促進するための不可欠な課題である。
【０００６】
現在の最先端技術では、最初のノッキングを検出するだけでメタン価の微分を検出することが可能になり、次いで必要な補正を実行することが可能になり、現在使用される技術は、音響センサによってノッキングを検出することであり、それによってエンジンシステムの空気入口を適合させ、その後のノッキングを防ぐためのエンジンシステムの修正を始動する。
【０００７】
ＬＮＧのメタン価を決定するために、気相の組成を決定することに基づく解決策が当業者に知られており、それによってメタン価に近似することが可能になる。気相のこの決定は、特に以下のように実施され得る。
−分光測定による光学的測定によるものであって、これにより、重炭化水素に対しては十分な精度を可能にするが、しかしメタン価の制御の観点からガスを遮断するために使用される窒素に対しては不可能になる（分光測定はＣＨ結合のみのスペクトル応答を測定する）。さらに、組成を決定するためのそのような方法は、今日ガソリンエンジンにのみ使用されている（したがって、ＬＮＧの要求を満たすためのその能力についての不確実性もある）。最後に、そのような方法は、光学センサのコストが約１００ドルであるので非常に高価である。
−温度および圧力センサとの経験的相関による計算によるが、そのような解決策のコストは、現在、ＬＮＧタンクのコストの３分の１以上であり、これはかなり思いとどまらせ
る要因である。さらに、現在知られている温度および圧力センサは、通常、ガソリンエンジンを対象としているが、必要に応じてＬＮＧにも適合されるべきである。
−クロマトグラフィーによる測定によるが、装置の体積およびコストのために、この解決法は船上タンクでは考えられないものになっている。
【発明の概要】
【０００８】
上記の欠点を克服するために、本出願人は、船上タンクのようなタンクに含まれる液化天然ガス（一般に頭字語、ＬＮＧで呼ばれる）のメタン価ＭＮをリアルタイムで計算する方法を実装した。メタン価を計算するためのこの方法は、多成分の液相中のＬＮＧの組成を決定することに基づいており、これは３種類の異なるセンサを有利に使用するに過ぎない。したがって、この決定は経済的であり、いかなる気相測定も含まない。
【０００９】
先行技術（上記で引用）から知られている解決策は、気相の組成を決定することからメタン価の値に近似することを目指しており、したがって、得られる値はエンジンに噴射されるガスの実際の品質には対応しない。実際、気相の組成は、成分の揮発性の相違により、液相の組成とは異なる。
【００１０】
それにもかかわらず、現在の噴射技術は液体ＬＮＧの気化を促進する。したがって、近似しなければならないのはこの液相の組成である。一旦組成が分かれば、選択された規格に応じて、計算によってメタン価に近似することが容易になる。
【００１１】
より具体的には、本発明は、したがって、
天然ガスを含むタンク（１）に含まれる液化天然ガスのメタン価ＭＮをリアルタイムで計算するためにコンピュータによって実装される方法であって、
前記天然ガスは、
−液体状態（１）の天然ガスの層（１１）であって、所与の時点ｔにおいてその温度Ｔ（ｔ）およびその密度ρ（ｔ）で定義され、気体状態（ｇ）の天然ガスの層（１２）とバランスがとれている、液体状態（１）の天然ガスの層（１１）と、
−気体状態（ｇ）の天然ガスの層（１２）と
に分配され、
前記方法は、所与の時点ｔで、以下のステップ、
Ａ．液体状態の前記天然ガス層の温度Ｔ０（ｔ）および密度ρ０（ｔ）、ならびに気体状態（ｇ）の天然ガス層の圧力Ｐｏ（ｔ）を測定により決定するステップと、
Ｂ．応力最小値計算アルゴリズムによって、タンクに含まれる液化天然ガスの組成を近似計算するステップであって、前記アルゴリズムが、以下のサブステップ
ｂ１）応力最小値を計算することによって、圧力Ｐｏ（ｔ）で密度ρ０（ｔ）の第１の組成、または圧力Ｐｏ（ｔ）で同じ密度ρ０（ｔ）を有する組成の第１の包絡線を決定するステップと、
ｂ２）応力最小値を計算することによって、温度Ｔ０（ｔ）および圧力Ｐｏ（ｔ）で同じ密度ρ０（ｔ）を有し、かつ圧力Ｐｏ（ｔ）で平衡温度Ｔｅｑ（ｔ）がＴ０（ｔ）に等しい組成物の第２の包絡線を決定するステップと、
ｂ３）ＬＮＧの実際の組成に近似する組成を提供する第２の包絡線の非特異点を決定するステップと、
ｂ４）前記近似の組成からメタン価ＭＮを計算するステップと
を有する、アルゴリズムからなることを特徴する方法に関する。
【００１２】
本発明による計算方法は、コンピュータによって実装される。
タンクとは、本発明の意味において、加圧タンク、または周囲圧力でのメタンのタンク、またはガスターミナルタイプのタンクを意味する。
平衡温度Ｔｅｑ（または気泡温度）とは、本発明の意味において、膨張空間の圧力にお
いて、加熱された場合に最初の気泡が現れる温度を意味する。
【００１３】
本発明による方法の範囲、特にステップＡでは、温度Ｔ０（ｔ）の測定は少なくとも１つの温度センサを用いて実行される。しかしながら、ＬＮＧタンク内での測定の均一性を保証するために、有利には、図１に示す（実施例１参照）ように、いくつかの温度センサ、例えば膨張空間内の１つのセンサおよび液体状態の天然ガスの層内の３つのセンサを使用することが可能である。
【００１４】
本発明による方法のステップＡの範囲では、液体状態の天然ガス層の密度ρ０（ｔ）を測定することは、少なくとも１つの密度センサを用いて実行される。単一のセンサが密度を測定するために十分である可能性があり、ＬＮＧが貯蔵タンク内の組成および温度に関して比較的均一であるからである。
【００１５】
本発明による方法のステップＡの範囲内では、気体状態（ｇ）の天然ガス層の圧力Ｐｏ（ｔ）を測定することは、少なくとも１つの圧力センサを用いて実行される。単一のセンサが圧力を測定するのに十分である可能性があり、圧力は気相中で均一であるからである。
【００１６】
本発明による方法は、３種類の異なるセンサによって行われた測定値（特に、天然ガスの液相の温度および密度ならびに天然ガスの気相の圧力）から、非常に低い誤差（５％未満）でＬＮＧの組成に近似することができるという利点を有する。
【００１７】
実際、理論が共有されるならば、その組成、したがってＭＮ価を決定しようとしているＬＮＧは、最高でＣ５および窒素の化合物を含み、それを含む通常９つの要素ＣＨ４、Ｃ２Ｈ６、Ｃ３Ｈ８、ｉＣ４Ｈ１０、ｎＣ４Ｈ１０、ｉＣ５Ｈ１２、ｎＣ５Ｈ１２、Ｃ６Ｈ１４、およびＮ２が列挙され得る。物理的な大きさが各成分のモル分率を関連させる方程式を提供すると仮定することによって、８つの異なる情報項目（最後の成分の分率は他の８つから推定される）を得ることが必要であろう。
【００１８】
この目的のために、本発明による方法は、ＬＮＧの液相の温度、密度、およびガスの圧力の３つの重要な測定可能なパラメータを共有し、それに理想的な「統計的」条件付き理論（ＬＮＧ業界およびＧＩＩＧＮＬ（液化天然ガス輸入業者の国際グループ［１］）による出版物で定義される）を加え、メタン価を求める点から考慮して、方程式（各成分のモル分率を関連させる）を解くために計算機（例えば、マイクロプロセッサまたはＰＣ）によって実装される非線形最適化アルゴリズムを使用する。
【００１９】
本発明による方法の主なステップは、計算機によって実装される応力最小値計算アルゴリズムによるタンク内に含まれる液化天然ガスの組成の近似計算のステップＢである。ユーザがＬＮＧの特性（組成、ＭＮ価）を知りたいと思うとすぐに、計算機は、温度、圧力および密度からＬＮＧの組成を決定するためにアルゴリズムを実行する。
【００２０】
この計算は、主な理論に基づいており、ＬＮＧはその気相とバランスがとれており、すなわち液体はその気泡曲線上で熱力学的に位置するということである。この仮定は、計算の基礎である。確かに、バランスのとれたＬＮＧ（ＬＮＧ貯蔵庫ではほとんどの場合）を仮定することは、ＬＮＧの温度、圧力、組成を達成することを可能にする。したがって、この式および密度の値（組成に応じて）は、組成を決定するために解決するシステムを形成する。
【００２１】
本発明による方法のアルゴリズムは、３つのステップ（ステップｂ１からｂ３）で、ＬＮＧの１組の可能な組成を提供することを可能にし、
−最初に、測定された圧力Ｐｏ（ｔ）で測定された密度ρ０（ｔ）を有するＬＮＧ組成物またはＬＮＧ組成物の包絡線が決定される（ステップｂ１）。次いで、
−ステップｂ１で計算された組成、または組成の包絡線から、圧力Ｐｏ（ｔ）での平衡温度Ｔｅｑ（ｔ）が、測定された温度Ｔ０（ｔ）に等しい、組成の第２の包絡線が決定される（ステップｂ２）。
−それにより、非特異点が決定され、ＬＮＧの実際の組成に近似する組成を提供する（ステップｂ３）。
【００２２】
ステップｂ１およびｂ２のそれぞれから得られる解の空間（すなわち可能な組成）を減らすために、これらは勾配型の応力最小値計算ステップ（共役勾配、推定ニュートン、内点など）である。
【００２３】
次に、ステップｂ３で得られた近似の組成からメタン価ＭＮが算出される（ステップｂ４）。この計算は簡単であり、ＬＮＧ業界の標準的な計算の１つ（会社のタイトルが付いたＡＶＬ法、ＳＡＥインターナショナルによって提案されたＣＡＲＢ法、ガス研究所のＧＲＩ法など）を使用する。
【００２４】
どの計算方法が使用されようとも、メタン価は０から１００の間の無次元数であり、値「０」は理論的に１００％の水素を含む組成物に相当し、値「１００」は１００％のメタンを含む組成物に相当する。
組成の計算誤差が約５〜１０％に近づくと、物理的な温度および圧力条件に従って、０．２％未満のメタン価の値の精度を達成することが可能である。
【００２５】
本発明による方法のステップｂ１）の第１の有利な実施形態によると、組成物の第１の包絡線を決定することが、以下の方程式を解くことによって、以下のように実行されることができ、
【数１】
0006858859.tif 000002
上記の方程式において、
ρはｘおよびＰｏから計算される密度を意味し、
ρ０は測定された密度を意味し、
Ｐｏは測定された圧力を意味し、
ｘはｘｉからなる合成ベクトルを意味し、
ｘｉは成分ｉのモル分率を意味し、
ｌおよびｕは、それぞれモル分率ｘｉの下限および上限であることを示す添字である。
【００２６】
本発明による方法のステップｂ１）の第２の有利な実施形態によると、組成物の第１の包絡線を決定することが、以下の方程式を解くことによって、以下のように実行されることができ、
【数２】
0006858859.tif 000003
上記の方程式において、
ρはｘおよびＰｏから計算される密度を意味し、
ρ０は測定された密度を意味し、
Ｐｏは測定された圧力を意味し、
ｘはｘｉからなる合成ベクトルを意味し、
ｘｉは成分ｉのモル分率を意味し、
ｌおよびｕは、それぞれモル分率ｘｉの下限および上限であることを示す添字である。
【００２７】
当業者にとって、他の既知の種類の応力最小値計算を使用することによって、上述の２つの有利な実施形態からステップｂ１）の他の実施形態を見つけることは明らかである。
【００２８】
本発明による方法のステップｂ２）の第１の有利な実施形態によると、組成物の第２の包絡線を決定することが、以下の方程式を解くことによって、以下のように実行されることができ、
【数３】
0006858859.tif 000004
上記の方程式において、
Ｔ０は測定温度を意味し、
Ｔｅｑは圧力ＰｏでのＬＮＧの平衡温度を意味し、Ｐｏ、ρ、ρ０、ｘ、ｘｉは上記に定義されたようなものである。
【００２９】
本発明による方法のステップｂ２）の第２の有利な実施形態によれば、組成の第２の包絡線を決定することは、以下の式
【数４】
0006858859.tif 000005
を解くことによって以下のように実行されることができ、上記の式において、
Ｔ０は測定温度を意味し、
Ｔｅｑは圧力ＰｏでのＬＮＧの平衡温度を意味し、Ｐｏ、ρ、ρ０、ｘ、ｘｉは上記に定義されたようなものである。
【００３０】
当業者にとって、他の既知の種類の応力最小値計算を使用することによって、上述の２つの有利な実施形態からステップｂ２）の他の実施形態を見つけることは明らかである。
【００３１】
本発明による方法のステップｂ３）の第１の有利な実施形態によれば、第２の包絡線の非特異点は、第２の包絡線の重心、または第２の包絡線の重心に最も近い第２の包絡線の点である。
【００３２】
本発明による方法のステップｂ３）の第２の有利な実施形態によれば、第２の包絡線の非特異点は、第２の包絡線の標準化組成物の再標準化平均として計算され得る。
ステップｂ３）の他の実施形態は、他の既知の種類の特異点を用いることによって、上述の２つの有利な実施形態から実施可能である。
【００３３】
本発明の他の有利な点および詳細は、限定しない例として提供され、異なる温度、密度および圧力センサを有する天然ガスを含むタンクの例を概略的に表す単一の添付の図を参照してなされる以下の説明から生じる。
より具体的には、この図は、
−所与の時点ｔにその温度Ｔ（ｔ）およびその密度ρ（ｔ）によって定義される液体状態１の天然ガスの層１１であって、天然ガスの層１２とバランスがとれている、天然ガスの層１１と、
−所与の時点ｔにおいて圧力Ｐ（ｔ）によって定義される、気体状態ｇの天然ガスの層１２と
に分配される天然ガスを含む水平タンク１を示す。
タンク１は、
−ＬＮＧの温度Ｔ０（ｔ）を決定するための複数の温度センサ２１、２２、２３、２４であって、これらのセンサは支持部２５上に配置され、これらのうちの１つのセンサ２１は、ＬＮＧの気相（または膨張空間）に配置され、これらのうちの３つのセンサは、液相に配置されている（ガスの温度は計算で考慮されない）、複数の温度センサ２１、２２、２３、２４と、
−液体状態のＬＮＧの密度ρ０（ｔ）を決定するための密度センサ２６（液体状態のＬＮＧが均一な密度を有する場合には単一のセンサで十分である）と、
−膨張空間の圧力Ｐｏ（ｔ）を決定するための圧力センサ３１と
を備える。
【００３４】
一旦各時点で測定されると、温度値Ｔ０（ｔ）、密度ρ０（ｔ）、および圧力Ｐｏ（ｔ）は、本発明による方法を実装する計算機に導入され、ＬＮＧに近似する組成を定義する
。これらのシミュレーションの結果は、以下の実施例１に提示されている。
本発明による計算方法は、以下の実施例においてより詳細に説明される。
【実施例】
【００３５】
以下の実施例は、当業者によるその読解および議論の理解を容易にするために、教育目的のために限定された数の点（６つの組成の包絡線）で作製されている。簡単に推定することによって、測定点の数を増やすことによって、その数の結果が得られることを当業者は容易に理解することができ、その精度は測定点の数と共に増すであろう。
【００３６】
本発明による方法の入力パラメータは以下の通りである。
−Ｔ０（ｔ）＝−１６０．７６℃
−ρ０（ｔ）＝４４８．１１ｇ／ｃｍ３
−Ｐｏ（ｔ）＝１１９０ｍｂａｒ
【００３７】
これらの値から、したがって、測定された圧力Ｐｏ（ｔ）で測定された密度ρ０（ｔ）を有するＬＮＧ組成が決定される（ステップｂ１）。この第１の組成は、以下の表１に提供されている。
【００３８】
【表１】
0006858859.tif 000006
【００３９】
この第１の組成では、平衡温度は従来の位相包絡線計算方法（例えば、Ｒａｃｈｆｏｒｄ−Ｒｉｃｅ方式）［２］によって計算されており、これは−１５８．３６℃である。検証のために密度ρ（ｔ）を再計算すると、Ｔ０（ｔ）＝−１６０．７６℃に対してρ（ｔ）＝４４８．１１ｇ／ｃｍ３である。
【００４０】
次に、ステップｂ１で計算されたこの第１の組成から、圧力Ｐｏ（ｔ）における平衡温度Ｔｅｑ（ｔ）が、測定された温度Ｔ０（ｔ）、すなわち−１６０．７６℃に等しい組成の包絡線が決定される（ステップｂ２）。この包絡線は６つの組成物Ｃ１〜Ｃ６を含み、それらは以下の表２に詳述されている。
【００４１】
【表２】
0006858859.tif 000007
【００４２】
組成物Ｃ１〜Ｃ６のこの包絡線について、（同様にＲａｃｈｆｏｒｄ−Ｒｉｃｅ方式のような従来の位相包絡線計算方法［２］によって）平衡温度が計算されている：Ｔｅｑ（ｔ）＝−１６０．７６℃であり、測定された温度Ｔ０であることから、計算が正しく行われていることを示す。
【００４３】
組成物Ｃ１からＣ６のこの包絡線から、ＬＮＧの実際の組成に近似の組成を提供する第２の包絡線から非特異点が決定される。次に、ステップｂ３で得られた近似の組成からメタン価ＭＮを算出する（ステップｂ４）。
【００４４】
第１の場合では、下記の表３に示されるように、重心が非特異点として取られるようにＬＮＧの実際の組成に近似の組成を提供する（ステップｂ３）。
【００４５】
【表３】
0006858859.tif 000008
【００４６】
次に、ステップｂ３で得られた近似の組成からメタン価ＭＮが算出される（ステップｂ４）。この計算は簡単であり、ＬＮＧ業界の標準的な計算方法（ＣＡＲＢ、ＧＲＩ  Ｈ／Ｃ、線形係数ＧＲＩ、およびＭＯＮ／ＫＩＷＡ、ＭＷＭ、ＡＶＬ）を使用し、これらは様々な経験およびエンジンの種類に基づく当業者に周知の経験的方法である。ＭＷＭ方式およびＡＶＬ方式は、当業界で最も一般的に使用されている周知の方式である。
【００４７】
これらの計算の結果は、表３の組成からメタン価を計算するための異なる方法のそれぞれについて下記の表４に示されている。
【００４８】
【表４】
0006858859.tif 000009
【００４９】
異なる方法は、異なる結果を提供することに留意されたい。メタン価を計算するために使用される方法について、今日当業界では合意がないので、これは通常のことである。様々な方法があり、各社は支持する傾向の方法を有する。それほど拡張的ではない他の計算方法が存在し、上記で決定された組成からすべてが使用される可能性もある。
【００５０】
第２の場合では、上記で得られた組成物Ｃ１〜Ｃ６の包絡線から、第２の包絡線の標準化された組成の再標準化された平均を非特異点と見なし、その組成が以下の表５に提供される。
【００５１】
【表５】
0006858859.tif 000010
【００５２】
次いで、メタン価ＭＮが、上記のＬＮＧ業界の標準的な計算方法（ＭＷＭ、ＡＶＬ）を使用することによって、ステップｂ３において得られたこの近似の組成（表５に提供された）から計算される（ステップｂ４）。メタン価計算の結果は、表５の組成から、使用された異なる方法のそれぞれについて下記の表６に示されている。
【００５３】
【表６】
0006858859.tif 000011
表６で計算された値は、表４の値に近いことが観察される。
【００５４】
表４と表６との間の差はすべて、１．３％であるＭＷＭ方式に関するものを除いて、０．２％未満である。ＭＷＭ方式は、安定した方法ではなく、非常に近い組成に対して絶対値１の差を提供することが当業界で知られているので、したがってこの結果は予想され得る。
【００５５】
（参考文献一覧）
［１］http://www.giignl.org/publications
［２］Procedure for use of electronic digital computers in calculating flash vaporization hydrocarbon equilibrium. RACHFORD Jr, H. H., & RICE, J. D. s.l.: Journal of Petroleum Technology, 1952, Vol. 4(10), 19-3.
【図１】
0006858859.tif 000016