(57)【要約】本発明は、液化天然ガス（LNG）の組成を特定する必要なしに、LNGを収容する非冷却式加圧タンク内の残存化学エネルギの量をリアルタイムに計算する方法およびシステムに関する。
2019519740.tif 000002
【特許請求の範囲】
【請求項１】
形状および寸法によって画定され、液化天然ガス（LNG）の層を収容する加圧タンク（1）内に含まれる残存化学エネルギEをリアルタイムに計算する方法であって、前記LNG層は、所与の瞬間tにおいて、前記タンク内でのLNG層の温度T、LNG層の密度ρ、およびLNG層
の高さhによって定義され、
前記方法は、所与の瞬間tにおいて、以下のステップ：
A．以下の測定によって、前記LNG層の特性パラメータを取得するステップと、
−  高さセンサ（2）を使用する前記タンク内の前記LNG層の前記高さhの測定、
−  温度センサ（3）を使用する前記温度Tの測定、
−  密度センサ（4）を使用する前記密度ρの測定、
B．前記タンク（1）内に収容された前記LNGの全質量mtを求めるステップと、
を含むアルゴリズムからなる、方法において、
前記アルゴリズムは、各瞬間tに対して、以下のステップ：
C．以下の式に従い、前記液体の前記温度および前記密度をパラメータとして用いる関
数fを使用して、前記LNGの質量総熱量GCVmassを計算するステップと、
【数１】
2019519740.tif 000012
D．以下の式に従って、前記残存化学エネルギEを計算するステップと、
【数２】
2019519740.tif 000013
をさらに含むことを特徴とする方法。
【請求項２】
・前記アルゴリズムは、前記タンク（1）を使用するオペレータ（7）から要求されたときに繰り返して行われるか、
・または、前記アルゴリズムは、所与の時間間隔Δtが経過した直後に自動的に行われ
るか、
のいずれかである、請求項1に記載の方法。
【請求項３】
前記タンク（1）に収容されたLNGの前記全質量mtの導出は、バランスまたは応力ゲージを使用する直接測定によって行われる、請求項1または2に記載の方法。
【請求項４】
前記タンク（1）に収容されたLNGの前記全質量mtの導出は、以下の式に従った計算によって行われる、請求項1または2に記載の方法。
【数３】
2019519740.tif 000014
式中、
−  hは、前記タンク内の前記LNG層の前記高さ
−  ρは、前記LNGの前記密度
−  gは、前記タンクの前記形状に関係する関数
【請求項５】
前記質量総熱量GCVmassを前記パラメータTおよびρに関係付ける関数fは、以下の形を
取る、請求項1〜4のいずれか一項に記載の方法。
【数４】
2019519740.tif 000015
式中、
−  Aは、所与の温度に対して一定値
−  Bは、組成とは無関係の定数
【請求項６】
形状および寸法によって画定され、液化天然ガス（LNG）の層を収容する加圧タンク（1
）内に含まれる残存化学エネルギEを請求項1〜5のいずれか一項に記載の方法に従って、
リアルタイムに計算するためのシステムであって、前記LNG層は、所与の瞬間tにおいて、前記タンク内でのLNG層の温度T、LNG層の密度ρ、およびLNG層の高さhによって定義され
、
前記システムは、
−  前記タンク（1）に設けられた高さセンサ（2）、温度センサ（3）、および密度セ
ンサ（4）に接続されることを意図され、請求項1〜5のいずれか一項に記載の方法の前記
アルゴリズムを実行することができる計算機（5）と、
−  請求項1〜5のいずれか一項に記載の方法の前記アルゴリズムによって得られた残存化学エネルギの量をオペレータに知らせるために、前記計算機（5）と相互作用するMMIインターフェイス（6）と、
を含むことを特徴とする、システム。
【請求項７】
車載システムであり、
・前記計算機（5）は、前記高さセンサ（2）、温度センサ（3）、および密度センサ（4）に接続された車載計算機であり、前記計算機（5）は、特に、本発明に記載の方法の前
記アルゴリズムを実行するように設計され、
・前記MMIインターフェイス（6）は、車載ダッシュボードタイプの車載インターフェイス、またはオフセットインターフェイスである、
請求項6に記載のシステム。
【請求項８】
液化天然ガスの層を収容し、高さセンサ（2）、温度センサ（3）、および密度センサ（4）を設けられた加圧タンク（1）を含む車両であって、請求項6または7に記載のシステムをさらに含むことを特徴とする車両。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、概略的には、液化天然ガス（LNG）の組成を特定する必要なしに、LNGを収容する非冷却式加圧タンク内の残存化学エネルギの量をリアルタイムに計算する方法およびシステムに関する。
【背景技術】
【０００２】
LNG燃料は、CO2排出物、汚染粒子、およびエネルギ密度の観点から、従来の燃料に対する単純かつ有効な代替物である。道路上、海上、またはレール上の運搬機など、LNG燃料
の使用に切り換える関係者が増えている。
【０００３】
しかし、従来の燃料とは異なり、LNGの体積エネルギ密度、すなわち、LNGの単位体積当たりに含まれるエネルギは一定ではない。これは、2つの別の現象から説明することがで
きる。第1に、LNGの温度は、残存熱入力により、非冷却式加圧タンクでの貯蔵期間全体にわたって上昇する。この場合に、この温度上昇により、（体積が最大で20％を超える範囲で増大することがある）流体の熱膨張が生じ、ひいては、流体のエネルギ密度が落ちる。
【０００４】
LNGのエネルギ密度の変化を説明する第2の現象はLNGの組成の変化である。LNGは精製された生産物ではなく、したがって、炭化水素の組成が、使用される保管所に応じて変わることがある。
【０００５】
非冷却式の貯蔵器に保管されたLNGの体積エネルギ密度の変動は、燃料消費の細かい監
視を必要とするシステムにおいて問題であり得る。通常、LNGで動く貨物自動車において
、600LのLNGを収容する同じ貯蔵器の場合に、同一組成のLNGに対して、LNGが重くかつ冷
たいかどうか、またはLNGが軽くかつ温かいかどうかに応じて、約15-20％のLNGの体積エ
ネルギ密度の違いが認められることがある。これは、実際上、比較例に示すように、始動時に同じ量のLNGが導入された場合に、走行キロメートル数に関して100キロメートル程度の差をもたらす。
【０００６】
現在のところ、LNGのタンク内に含まれる残存エネルギについて、加圧タンクのオペレ
ータにリアルタイムで知らせる問題解決策はない。オペレータにとって入手可能な情報は、ガス組成物の圧力、（最良の場合の）LNGの温度、およびタンクの充填高さのみである
。
【０００７】
一般的に、燃料供給装置によるタンクへの充填時に、（供給装置から供給される）LNG
とLNGの転送質量との既知の最新組成を使用して、国際規格ISO 6976.1995に準拠したエネルギ計算が行われる。この計算は、金融取引のための基準として使用される。このように、LNGが実質的にメタン、エタン、プロパン、イソブタン、n−ブタン、イソ−ペンタン、n-ペンタン、および窒素からなることを前提として、式（1）による、LNGの燃焼温度におけるこの計算から、LNGの総熱量GCVmassが求められ、
【数１】
2019519740.tif 000003
式中、
−  GCVmassは、LNGの熱量を表し、
−  Tc  GCVが計算される燃焼温度、
−  xj  混合物の成分jのモル分率、
−  Mj  成分jのモル質量、
−  M  標準規格NF EN ISO 976から得られるLNGのモル質量、
−  GCVmass＿j  ISO 6976.1995の図表から得られる成分jの総熱量、
である。
【０００８】
ただし、この計算は、LNGの組成に依存する。しかし、この組成は、特定するのに手間
がかかることがある。実際上、クロマトグラフの導入が必要である。
【０００９】
タンク内に含まれるエネルギに関するリアルタイムでの情報がないのは、以下のいくつかの理由から問題である。
−  供給管理：現在のところ、特定のタンク（特に、貨物自動車のタンク）に対するLNGの供給管理は、貯蔵器に残った液体の体積のみに基づいている。しかし、タンクに接続
されたユニットによって要求されたエネルギに基づく管理の方が一貫性がある。それが、例えば、引き続き走行できるキロメートル数を評価するために必要とされるデータだからである。
−  燃料不足および燃料切れの回避：LNGのエネルギ密度に応じて、ユニットの体積消
費量が急激に多くなることがあり、その理由は、この場合に、同じエネルギ量を得るために、より多量のLNGが必要とされるからである。オペレータが予定しないこの変動は、燃
料の予期しない不足、ひいては燃料切れを引き起こすことがある。
−  オペレータのトレーニング：LNG燃料市場は、比較的規模が小さい。市場関係者は
、ほとんどが、LNGを取り扱うのに適したトレーニングを受け、よい練習をした専門家で
ある。しかし、市場が急速に成長した場合に、トレーニング不足の関係者が、LNGの消費
に対処し、かつ／またはLNGの消費を管理する必要がある。タンク内に含まれるエネルギ
量が分かることで、これらのオペレータが容易に理解できる大きさ（例えば、残りのキロメートル数）を簡単に計算することが可能になる。
【発明の概要】
【００１０】
これを踏まえて、LNG燃料の発展を確実にするために、本出願人は、タンク内で測定さ
れる熱力学的パラメータ（タンク内のLNGの密度、LNGの層の温度および高さ）だけを使用して、タンクに収容されたLNGの組成を知ることなく、LNGのエネルギ含量をリアルタイムに、より良好に計画することを可能にする問題解決策を用意した。
【００１１】
特に、本発明は、目的を達成するために、形状および寸法によって画定され、液化天然ガス（LNG）の層を収容する非冷却式加圧タンク内に含まれる残存化学エネルギEをリアルタイムに計算する方法を有し、前記LNG層は、所与の瞬間tにおいて、前記タンク内での温度T、密度ρ、および高さhによって定義され、
前記方法は、所与の瞬間tにおいて、以下のステップ：
A．以下の測定によって、LNG層の特性パラメータを取得するステップと、
−  高さセンサを使用する、タンク内のLNG層の高さhの測定、
−  温度センサを使用する温度Tの測定、
−  密度センサを使用する密度ρの測定、
B．タンク内に収容されたLNGの全質量mtを求めるステップと、
を含むアルゴリズムからなり、
前記方法は、前記アルゴリズムが、各瞬間tに対して、以下のステップ：
C．以下の式に従い、液体の温度および密度をパラメータとして用いる関数fを使用して、LNGの質量総熱量GCVmassを計算するステップと、
【数２】
2019519740.tif 000004
D．以下の式に従って残存化学エネルギEを計算するステップと、
【数３】
2019519740.tif 000005
をさらに含むことを特徴とする。
【００１２】
天然ガスの質量総熱量（mass gross caloric value）という用語は、本発明において、一定の圧力および所与の温度で空気中に含まれる、対象となる単位質量の天然ガスの完全燃焼によって供給される熱量を意味する。質量総熱量は、単位質量の燃料当たりの熱量として表される（本発明の枠組みではkWh/m3）。
【００１３】
タンクの形状および寸法、LNGの温度、LNGの層の高さ、ならびにLNGの密度などの入力
情報を用いて、本発明による方法のアルゴリズムは、任意のタンク内に含まれる残存化学エネルギの実際の量を即座に計算することを可能にする。
【００１４】
さらに、この方法は、LNGのGCVmassを求めるために、クロマトグラフまたは熱量計の使用を必要とするLNGの組成の特定を必要としないので、この方法を用意するのは簡単であ
る。実際上、通常では、LNGの質量GCVは、一般に、LNGが単にメタン、エタン、プロパン
、イソブタン、n-ブタン、イソ−ペンタン、n-ペンタン、および窒素で構成されるという近似を行うことで、LNGの組成に従って計算される。
【００１５】
本発明による方法の場合、基本的にLNGの正確な組成を用いないことで生じる誤差は、
最大で約3％であり、これは、対象とする組成の温度と同じ温度での、重LNG（メタン以外の炭化水素を10％を超えて含有）のGCVmassと軽LNG（純粋なメタンを99％を超えて含有）との間に認められる差である。
【００１６】
比較として、LNGのGCVmassを求めるための、本発明とは異なる方法で生じる誤差は、LNGのGCVmassが不正確な温度で求められる場合に、組成が正確な場合を含めても、約20％の値に急伸することがある。
【００１７】
有利にも、前記アルゴリズムは、前記タンクを使用するオペレータの要求時に繰り返して行うか、または所与の時間間隔Δtが経過した直後に自動で行うことができ、この間隔
は、例えば、約1秒とすることも、または適用可能な場合に、使用されるセンサ技術に応
じて、潜在的遅延を考慮して最適に設定することもできる。
【００１８】
LNGの全質量の導出は、様々な方法で行うことができる。
【００１９】
導出するための第1の方法によれば、タンクに収容されたLNGの全質量mtは、有利にも、バランスまたは応力ゲージを使用する直接測定によって求めることができる。
【００２０】
別の有益な実施形態によれば、タンクに収容されたLNGの全質量mtの導出は、以下の式
に従った計算により行うことができる。
【数４】
2019519740.tif 000006
式中、
−  hは、タンク内のLNGの層の高さ
−  ρはLNGの密度
−  gは、タンクの形状に関係する関数であり、体積に対して一定の値を付与する
【００２１】
全質量mtを求めるこの方法は、特に、質量の直接測定が、タンクで実施するのに煩雑な場合、例えば、タンクが測定時に動く場合に使用することができる。
【００２２】
有利にも、質量総熱量GCVmassをパラメータTおよびρに関係付ける関数fは、以下の形
を取ることができる。
【数５】
2019519740.tif 000007
式中、
−  Aは所与の温度に対して一定値
−  Bは、組成とは無関係の定数
【００２３】
関数fに存在する2つの定数値は、LNG Industry magazine 2014などの業界出版物、または科学文献で規定されている。
【００２４】
本発明はまた、目的を達成するために、形状および寸法によって画定され、液化天然ガス（LNG）の層を収容する加圧タンク内に含まれる残存化学エネルギEを本発明の方法に従ってリアルタイムに計算するシステムを有し、LNGの前記層は、所与の瞬間tにおいて、前記タンク内での温度T、密度ρ、および高さhによって定義され、
前記システムは、
−  前記タンクに設けられた高さセンサ、温度センサ、および密度センサに接続されることを意図され、本発明による方法のアルゴリズムを実行することができる計算機と、
−  接続された計算機を用いて本発明による方法のアルゴリズムが実行された場合に、このアルゴリズムによって得られた残存化学エネルギの量をオペレータに知らせるために、前記計算機と相互作用するMMIインターフェイスと、
を含むことを特徴とする。
【００２５】
MMIインターフェイスという用語は、本発明において、活動に関する適切な決定を行う
ために、使用者が、残存エネルギの量に関する情報の任意の可聴信号または機械信号によって、見る、または通知を受けることを可能にするマン−マシンインターフェイスを意味する。
【００２６】
本発明の枠組み内で使用できるMMIインターフェイスとして、車両のダッシュボード、
コンピュータキーボード、LED表示灯、タッチスクリーンおよびタブレット、拡声器など
を挙げることができる。
【００２７】
本発明によるシステムの有益な実施形態によれば、このシステムは車載システムとすることができ、この車載システムでは、
・計算機は、前記高さセンサ、温度センサ、および密度センサに接続された車載計算機とすることができ、計算機は、特に、本発明による方法のアルゴリズムを実行するように設計され、
・MMIインターフェイスも車載式とするか、あるいは（例えば、車両が制御センタに接
続されている場合にオフセットとすることができる。
【００２８】
このMMIインターフェイスは、車載式の場合に、本発明の方法に従って計算された自律
性（autonomy）の持続時間をオペレータ（ここでは運転者）に知らせるために、特に、前記車載計算機と相互作用する車載ダッシュボードタイプとすることができる。
【００２９】
特に、本発明による方法のアルゴリズムを実行するように設計された計算機という用語は、本発明において、専用の記憶装置およびインターフェイスマザーボードと組み合わせたプロセッサを含む車載コンピュータであって、これらの要素のすべてが、機械的、熱力学的、および電磁的耐性の観点から見た「車載コンピュータ」ユニットの頑強性を付与し、したがって、LNG車両で使用するために「車載コンピュータ」ユニットが適合するのを
可能にするように組み立てられた車載コンピュータを意味する。
【００３０】
本発明によるシステムは、オペレータが、LNGの取扱いに適する任意のトレーニングを
受けていない場合でさえ、タンク内内に含まれる残存化学エネルギ量の値をオペレータにとって容易に利用可能にする。本発明によるシステムはまた、車載コンピュータなどのサードパーティシステムにこの値を供給するのを可能にする。
【００３１】
有利にも、システムは、タンク内に収容されたLNGの全質量を直接測定するために、バ
ランスまたは応力ゲージをさらに含むことができる。
【００３２】
最後に、本発明は、目的を達成するために、液化天然ガスの層を収容し、高さセンサ、温度センサ、および密度センサを設けられた加圧タンクを含む車両（陸上用、海上用、または空中用）をさらに有し、前記車両は、本発明によるシステムをさらに含むことを特徴とする。
【００３３】
本発明によるシステムにより、この車両は、LNGの取扱いに関するきめ細かなトレーニ
ングを受けていないオペレータが容易に使用することができる。実際上、このシステムは、タンク内の残存エネルギの値を表示するか、またはコンピュータに残存エネルギの値を送ることを可能にし、次いで、コンピュータは、タンクにもう一度補充することが必要になる前に、残りのキロメートル数を残存エネルギから推定することができる。
【図面の簡単な説明】
【００３４】
非限定的な例として提示され、添付図を参照する以下の説明を読むことで、本発明の他の利点および詳細が明らかになるであろう。
【００３５】
【図１】所与の温度および組成において、液体天然ガスの密度に対応したLNGの熱量のいくつかの測定結果を示している。
【図２】本発明による測定システムの特定の実施形態の図を示している。
【図３】本発明の枠組み内で使用できる非冷却式加圧タンクの例の図面を示し（円筒形の水平タンクの場合）、このタンクに収容されたLNGの質量を計算するのを可能にする関数g(h)を定めるのを可能にする様々なパラメータを示している。
【図４】本発明の枠組み内で使用できる非冷却式加圧タンクの例の図を示し（球形タンクの場合）、このタンクに収容されたLNGの質量を計算するのを可能にする関数g(h)を定めるのを可能にする様々なパラメータを示している。
【図５−７】それぞれ、円筒形であり、水平であるLNGのタンクを移送する車両のダッシュボードのスクリーン取込み図であり、本発明の方法に従った残存化学エネルギEの計算に使用される入力データと、さらにはこの計算の結果とを示している。
【発明を実施するための形態】
【００３６】
図1は、所与の温度（−160℃）において、LNGの様々な密度値に対して行われた総熱量
の一連の測定の結果を示している。これらの測定点は、−160℃でのこの特定の場合にお
いて、式f（ρ）＝0．0283ρ−0．7791で表される回帰直線によって、十分に関係付ける
ことができる（相関係数R2＝0．957）。したがって、この式fは、LNGが−160℃の温度に
あるときに、LNGのGCVmassを求めるための相関関数として使用することができる。
【００３７】
図2は、タンク1が円筒形であり、垂直である場合の本発明の特定の実施形態の間略図を示している。連続して行うことができる測定が行われる場合に、時間間隔Δtが経過した
後、またはオペレータ7からの命令が届いた後、タンク内に存在する密度センサ4、温度センサ3、および高さセンサ2は、タンク内の液体の温度と、この液体の密度と、さらには高さとの値を読み込む。次いで、この情報は、オペレータ7が、マン−マシンインターフェ
イス（MMI）6を介して、タンク1の形状、およびタンク1の特徴的な寸法、この特定の場合
には半径を前もって入力している計算機5に送られる。これは、計算機5が、タンク内に収容されたLNGの全質量mtを求めるために使用される関数g(h)を定めることを可能にする。
【００３８】
図3は、水平に置かれた円筒形タンクの図を示している。この場合に、このタンク内のLNG層の体積の計算は、円盤の部分面積を計算するのと同様である。この場合に、関数g(h)は以下である。
【数６】
2019519740.tif 000008
【００３９】
タンクが垂直に置かれた場合、g(h)は、単にg(h)＝π×R2×hである。
【００４０】
図4は、球形タンクを示している。この場合に、このタンク内のLNG層の体積の計算は、球形キャップを計算するのと同様である。この場合に、関数g(h)は以下である。
【数７】
2019519740.tif 000009
【００４１】
この情報を使用して、次いで、計算機5は、タンク1内に収容されたLNGの全質量mtとLNGの総熱量GCVmassとを計算し、この場合に、これらの値により、計算機が、測定時にタン
ク内内に含まれる残存エネルギEの値を得るのが可能になる。次いで、残存エネルギEの値は、MMI6を介してオペレータに供給することも、または残りのキロメートル数などの容易に理解できる情報を得るために再処理することもできる。
【００４２】
本発明は、下記に実施例でより詳細に示される。
【実施例１】
【００４３】
［例1（比較）］
この実施例は、非冷却式の貯蔵器に貯蔵されたLNGの体積エネルギ密度の変動を示して
いる。
【００４４】
この場合に、標準規格ISO 6976:1995の式（1）を使用する計算によって、600L（すなわち、0．6m3）のLNGを収容した貯蔵器内で、重く冷たいLNG（ケースa）：3barで平衡）の
場合と、同じ組成であるが軽く温かいLNG（ケースb）：14barで平衡）の場合とにおける
残存化学エネルギEが求められる。
【００４５】
ケースa）：重く冷たいLNG（3barで平衡）
前提として、LNGは、表1で下記に示す以下の組成を有する。
【００４６】
【表１】
2019519740.tif 000010
【００４７】
−  燃焼条件
○燃焼温度Tc＝0℃
○圧力：1．01325bar
−  質量GCV（Tc）＝14．99kWh／kg、標準規格ISO 6976:1995の式に準拠して計算
−  LNGの温度T＝−147．07℃
−  密度＝443．7153kg／m3
E＝0．6×密度×GCVmass＝3990kWh
【００４８】
ケースb）：軽く高温のLNG（14barで平衡）
LNGは、下記に表2で示す組成と同じ組成を有する。
【００４９】
【表２】
2019519740.tif 000011
【００５０】
−  燃焼条件
○燃焼温度Tc＝0℃
○圧力：1．01325bar
−  質量GCV（Tc）＝15．37kWh／kg、標準規格ISO 6976:1995の式に準拠して計算
−  LNGの温度T＝−112．5℃
−  密度＝355．65kg／m3
E＝0．6×密度×GCVmass＝3279kWh
【００５１】
結果的に、ケースa）およびケースb）でそれぞれ計算されたエネルギ値E間に17％を超
える差が認められる。言い換えると、600リットルの同じ初期体積のLNGに対するエネルギのこの差は、貯蔵器に導入されるLNGが、冷たく重い場合に（ケースa）、ケースb）で走
行するキロメートル数と比較して、さらに100キロメートル走行できるようにする。
【００５２】
［例2（本発明）］
図5〜7は、それぞれ、円筒形であり、水平であるLNGのタンクを移送する車両のダッシ
ュボードのスクリーン取込み図であり、本発明の方法に従った残存化学エネルギEの計算
に使用される入力データと、さらにはこの計算の結果とを示している。
【００５３】
特に、図5は、タンクに特有の入力データを示すダッシュボードのスクリーン取込み図
である。
−  形状：タンクを担持する車両内に水平に配置されたシリンダ、
−  寸法：
○長さ：1．2m、
○直径：0．7m、
【００５４】
図6は、LNGの層に特有の入力データを示すダッシュボードのスクリーン取込み図である。
−  温度T：−152．2℃、
−  密度ρ：420．2kg／m3、
−  高さh：0．501m、
【図１】
2019519740.tif 000016
【図２】
2019519740.tif 000017
【図３】
2019519740.tif 000018
【図４】
2019519740.tif 000019
【図５】
2019519740.tif 000020
【図６】
2019519740.tif 000021
【図７】
2019519740.tif 000022
【国際調査報告】
2019519740.tif 000023
2019519740.tif 000024
2019519740.tif 000025
2019519740.tif 000026