(57)【要約】【課題】  ＬＮＧ貯蔵タンク内の貯残ＬＮＧのサンプリング及び分析作業を行なわずに計算のみにより貯残ＬＮＧ密度を求めるＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法を提供する。【解決手段】  貯蔵しているＬＮＧから蒸発したメタンガスを外部に払い出しているＬＮＧ貯蔵タンク１１内の貯残ＬＮＧの密度の算出方法において、ＬＮＧ貯蔵タンク１１を含むＬＮＧ貯蔵設備１０内のＬＮＧからの理論メタン蒸発量を求め、理論メタン蒸発量からＬＮＧ貯蔵タンク１１内のＬＮＧの組成割合を求めてＬＮＧ平均分子量及びモル容積をそれぞれ算出し、ＬＮＧ平均分子量をモル容積で除して貯残ＬＮＧの密度を求める。【選択図】        図１
2006308384.tif 000002
【特許請求の範囲】
【請求項１】
貯蔵しているＬＮＧから蒸発したメタンガスを外部に払い出しているＬＮＧ貯蔵タンク内の貯残ＬＮＧの密度の算出方法において、
前記ＬＮＧ貯蔵タンクを含むＬＮＧ貯蔵設備内のＬＮＧからの理論メタン蒸発量を求め、該理論メタン蒸発量から前記ＬＮＧ貯蔵タンク内のＬＮＧの組成割合を求めてＬＮＧ平均分子量及びモル容積をそれぞれ算出し、前記ＬＮＧ平均分子量を前記モル容積で除して前記貯残ＬＮＧの密度を求めることを特徴とするＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＬＮＧを貯蔵しているＬＮＧ貯蔵タンクに新たにＬＮＧを受入れる際に、ＬＮＧ貯蔵タンク内でのＬＮＧの層化を防止するためにＬＮＧ貯蔵タンク内の貯残ＬＮＧ（受入れ前ＬＮＧ）レベルを決定する根拠となるＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法に関する。
【背景技術】
【０００２】
ＬＮＧは、メタンを主成分としエタン、プロパン等の他の炭化水素を含む混合物であり、その組成は産地により異なっている。このため、ＬＮＧの密度も産地により異なることになる。また、ＬＮＧをＬＮＧ貯蔵タンクに受入れて払出しを行なっていると、ＬＮＧ貯蔵タンク内では沸点の低いメタンの蒸発が優先するためＬＮＧの組成が徐々に変化し、貯蔵中のＬＮＧの密度は徐々に増加する。従って、ＬＮＧが残留しているＬＮＧ貯蔵タンクに新たにＬＮＧを受入れる場合、貯残ＬＮＧの密度と受入れＬＮＧの密度の間には差が生じることになる。その結果、例えば、密度の高いＬＮＧが残留しているＬＮＧ貯蔵タンクに密度の低いＬＮＧを受入れる場合、密度の高いＬＮＧが下層に、密度の低いＬＮＧが上層に存在するという層化が生じ、その状態が長期間継続するとロールオーバーと呼ばれる現象が発生し、ボイルオフガスが多量に発生してＬＮＧ貯蔵タンク内の圧力が急上昇するという問題が生じる。
【０００３】
そこで、受入れＬＮＧの密度とＬＮＧ貯蔵タンク内の貯残ＬＮＧの密度との間の密度差及びＬＮＧ貯蔵タンク内の貯残ＬＮＧの液面高さの関係が層化の発生に及ぼす影響を予め把握しておき、ＬＮＧをＬＮＧ貯蔵タンクに受入れる際にＬＮＧ貯蔵タンク内の貯残ＬＮＧの密度を算出し、密度差と貯残ＬＮＧの液面高さの関係から層化の発生を判定し、層化が発生すると判定された場合は、貯残ＬＮＧの液面高さを減少させて層化を防止するようにしていた。ここで、貯残ＬＮＧの密度は、ＬＮＧ貯蔵タンクから貯残ＬＮＧをサンプリングし貯残ＬＮＧの組成割合を求めてＬＮＧ平均分子量及びモル容積をそれぞれ算出し、ＬＮＧ平均分子量をモル容積で除して求めている（例えば、非特許文献１参照）。
【０００４】
【非特許文献１】「ＬＮＧ−６記念出版集  ＬＮＧのすべて」、ＬＮＧ便覧、日本ＬＮＧ会議、昭和５６年５月、ｐ．３３３
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、貯残ＬＮＧの密度を求めるには、貯残ＬＮＧを採取するサンプリング、サンプリングした貯残ＬＮＧの組成分析に手間がかかるという問題が生じる。更に、ＬＮＧ貯蔵タンクから貯残ＬＮＧをサンプリングするためには払い出し装置を使用せねばならず、払い出し装置を稼動させるエネルギーが必要になるという問題が生じる。更に、ＬＮＧサンプリング時に気化したメタンガスが大気中に放散され、地球環境の観点からも問題が生じる。
【０００６】
本発明はかかる事情に鑑みてなされたもので、ＬＮＧ貯蔵タンク内の貯残ＬＮＧのサンプリング及び分析作業を行なわずに計算のみにより貯残ＬＮＧ密度を求めるＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
前記目的に沿う本発明に係るＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法は、貯蔵しているＬＮＧから蒸発したメタンガスを外部に払い出しているＬＮＧ貯蔵タンクを含むＬＮＧ貯蔵設備内のＬＮＧからの理論メタン蒸発量を求め、該理論メタン蒸発量から前記ＬＮＧ貯蔵タンク内のＬＮＧの組成割合を求めてＬＮＧ平均分子量及びモル容積をそれぞれ算出し、前記ＬＮＧ平均分子量を前記モル容積で除して前記貯残ＬＮＧの密度を求める。
【発明の効果】
【０００８】
請求項１記載のＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法においては、貯残ＬＮＧ密度を理論メタン蒸発量から算出したＬＮＧ貯蔵タンク内のＬＮＧの平均分子量及びモル容積を用いて算出するので、ＬＮＧ貯蔵タンク内の貯残ＬＮＧのサンプリング及び分析作業といった受入れ前作業の軽減及びＬＮＧの払い出し装置の運転コストの削減ができ、ＬＮＧの受入れコストを低減することが可能になる。更に、ＬＮＧのサンプリングに伴うメタンガスの大気中への放散が防止でき、地球温暖化防止にも寄与できる。
【発明を実施するための最良の形態】
【０００９】
続いて、添付した図面を参照しつつ、本発明を具体化した実施の形態につき説明し、本発明の理解に供する。
ここで、図１は本発明の一実施の形態に係るＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法が適用されるＬＮＧ貯蔵設備の説明図、図２は同ＬＮＧ貯蔵設備におけるメタンガス処理量とメタンガス発生量の関係を示す模式図である。
【００１０】
図１に示すように、本発明の一実施の形態に係るＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法を適用するＬＮＧ貯蔵設備１０は、ＬＮＧを貯蔵するＬＮＧ貯蔵タンク１１と、ＬＮＧ貯蔵タンク１１にＬＮＧを受入れるＬＮＧ受入れ配管１２と、タンクローリー等の輸送タンクにＬＮＧを払い出すＬＮＧポンプ１３を備えた出荷設備配管１４と、ＬＮＧ貯蔵タンク１１内のメタンガスを外部に払い出すボイルオフガス配管系１５とを備えている。
【００１１】
また、ＬＮＧ貯蔵タンク１１内には、貯残ＬＮＧの高さ、すなわち、貯蔵しているＬＮＧの液位を検出する、例えばディスプレーサ１６を備えた液位検出計１７と、ＬＮＧ貯蔵タンク１１の高さ方向に間隔を置いて設けられた、例えば、熱電対等の第１の温度センサ１８備えた温度計１９が設けられている。更に、ＬＮＧ貯蔵タンク１１の底部には複数の第２の温度センサ２０が取付けられている。そして、液位検出計１７で検出されたＬＮＧの液位と、温度計１９及び温度センサ２０で検出された温度は、図示しない演算処理装置に入力されるようになっている。
【００１２】
その結果、温度センサ２０の検出値から、例えば、ＬＮＧ貯蔵タンク１１の底部の平均温度を求めることができる。また、温度計１９の各検出値と液位検出計１７の液位検出値を組み合わせることで、ＬＮＧ貯蔵タンク１１の側壁でＬＮＧと接触している部分（液接触側壁部）の平均温度、メタンガスと接触している部分（ガス接触側壁部）の平均温度、及び、ＬＮＧ貯蔵タンク１１の天井部でメタンガスと接触している部分（ガス接触天井部）の平均温度の時間変化をそれぞれ求めることができる。
【００１３】
ボイルオフガス配管系１５は、ＬＮＧ貯蔵タンク１１からメタンガスを図示しないガス送り出し本管に移送するもので、ＬＮＧ貯蔵タンク１１に接続するガス取り出し管２１と、ガス取り出し管２１の他端が接続されガス取り出し管２１を介してＬＮＧ貯蔵タンク１１からのメタンガスを冷却するサクションクーラー２２と、サクションクーラー２２の出口側に連結管２３を介して接続されるガス圧縮機２４と、ガス圧縮機２４の出口に一端が接続され他端がガス送り出し本管に接続されて図示しない流量計を備えたガス圧縮機出口管２５を有している。
ここで、サクションクーラー２２には、出荷設備配管１４と連通するＬＮＧ供給配管２６によりＬＮＧが供給され、ガス取り出し管２１を介して流入するメタンガスとの間で熱交換が行なわれる。そして、ＬＮＧから蒸発した天然ガスは、ガス圧縮機２４に供給されるようになっている。なお、符号２７はＬＮＧ受入れ側遮断弁を指す。
【００１４】
次に、本発明の一実施の形態に係るＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法について説明する。
貯残ＬＮＧ密度Ｄは、ＬＮＧ貯蔵設備１０内のＬＮＧからの理論メタン蒸発量Ｇを求め、理論メタン蒸発量ＧからＬＮＧ貯蔵タンク１１内の貯残ＬＮＧの組成割合を求めてＬＮＧ平均分子量及びモル容積をそれぞれ算出し、ＬＮＧ平均分子量をモル容積で除すことにより求めている。ここで、ＬＮＧ貯蔵設備１０内のＬＮＧは、ＬＮＧ貯蔵タンク１１内のＬＮＧと、ＬＮＧ受入れ配管１２の冷却用にＬＮＧ受入れ配管１２内に注入しているＬＮＧと、出荷設備配管１４の冷却用に出荷設備配管１４内に注入しているＬＮＧと、ＬＮＧ供給配管２６を介してサクションクーラー２２の図示しない熱交換配管内に注入しているＬＮＧを指す。
【００１５】
理論メタン蒸発量Ｇは、ＬＮＧ貯蔵設備１０に外部からの入熱により生じるメタンの熱蒸発量と、ガス圧縮機２４によりＬＮＧ貯蔵タンク１１内のメタンガスが吸引されて外部に払い出されることに伴うＬＮＧ貯蔵タンク１１内の圧力降下によるメタン強制蒸発量ＧF との総和であるメタンガス発生量に一致する。
【００１６】
ここで、メタンの熱蒸発量は、ＬＮＧ貯蔵タンク１１の底部からの入熱に伴うメタン蒸発量ＧB と、ＬＮＧ貯蔵タンク１１の液接触側壁部からの入熱に伴うメタン蒸発量ＧS と、ＬＮＧ貯蔵タンク１１のガス接触側壁部からの入熱に伴ってＬＮＧ貯蔵タンク１１内のＬＮＧの液面から蒸発するメタン蒸発量ＧW と、ＬＮＧ貯蔵タンク１１のガス接触天井部からの入熱に伴ってＬＮＧ貯蔵タンク１１内のＬＮＧの液面から蒸発するメタン蒸発量ＧC と、ＬＮＧ受入れ配管１２の外側面からの入熱によりＬＮＧ受入れ配管１２内のＬＮＧから生じるメタン蒸発量ＧR と、出荷設備配管１４の外側面からの入熱により出荷設備配管１４内のＬＮＧから生じるメタン蒸発量ＧY と、サクションクーラー２２での熱交換により発生するメタン蒸発量ＧM との総和となる。
【００１７】
一方、メタンの各蒸発量ＧF 、ＧB 、ＧS 、ＧW 、ＧC 、ＧR 、ＧY 、ＧM は、ＬＮＧ貯蔵タンク１１内のＬＮＧが蒸発源となる蒸発量ＧF 、ＧB 、ＧS 、ＧW 、ＧC と、ＬＮＧ貯蔵タンク１１を除いた部分のＬＮＧが蒸発源となる蒸発量ＧR 、ＧY 、ＧM に分けられる。このため、ＬＮＧ貯蔵タンク１１内のＬＮＧからの総蒸発量を第１の理論メタン蒸発量ＧT 、ＬＮＧ貯蔵タンク１１を除いた部分のＬＮＧからの総蒸発量を第２の理論メタン蒸発量ＧA とすると、図２に示すように、理論メタン蒸発量Ｇは、第１の理論メタン蒸発量ＧT と第２の理論メタン蒸発量ＧA の和となり
Ｇ＝ＧT ＋ＧA ・・・・・（１）
と表される。ここで、第１の理論メタン蒸発量ＧT と第２の理論メタン蒸発量ＧA はそれぞれ次式のようになる。
ＧT ＝ＧB ＋ＧS ＋ＧF ＋ＧW ＋ＧC ・・・・・（２）
ＧA ＝ＧR ＋ＧY ＋ＧM ・・・・・（３）
【００１８】
ＬＮＧ貯蔵タンク１１の底部からの入熱は、ＬＮＧ貯蔵設備１０の周囲の大気温度ＴA とＬＮＧ貯蔵タンク１１の底部の平均温度ＴB との差に比例するので、メタン蒸発量ＧB （単位：ｔ／ｈ）は
ＧB ＝Ｃ1 ・（ＴA −ＴB ）・・・・・（４）
と表せる。ここで、Ｃ1 はＬＮＧ貯蔵タンク１１の寸法、ＬＮＧ貯蔵タンク１１の底部の断熱構造の特性、及び熱伝達特性等を反映した設備定数である。また、平均温度ＴB は、第２の温度センサ２０の各検出値の平均値を求めることにより得られる。
【００１９】
また、ＬＮＧ貯蔵タンク１１の液接触側壁部からの入熱は、ＬＮＧ貯蔵タンク１１内に貯蔵されているＬＮＧの液位ＬV に比例すると共に、ＬＮＧ貯蔵タンク１１の周囲の大気温度ＴA とＬＮＧ貯蔵タンク１１の液接触側壁部の平均温度ＴS との差に比例するので、メタン蒸発量ＧS （単位：ｔ／ｈ）は、
ＧS ＝Ｃ2 ・ＬV ・（ＴA −ＴS ）・・・・・（５）
と表せる。ここで、Ｃ2 はＬＮＧ貯蔵タンク１１の寸法、ＬＮＧ貯蔵タンク１１の側部の断熱構造の特性、及び熱伝達特性等を反映した設備定数である。また、液位ＬV は液位検出計１７による測定から得られ、平均温度ＴS は、温度計１９の各検出値と液位検出計１７の液位検出値から判明する液接触側壁部の範囲に存在する第１の温度センサ１８の各検出値の平均値を求めることにより得られる。
【００２０】
更に、ＬＮＧ貯蔵タンク１１のガス接触側壁部からの入熱は、ＬＮＧ貯蔵タンク１１内のメタンガスが存在する空間部の高さに比例すると共に、ＬＮＧ貯蔵タンク１１の周囲の大気温度ＴA とＬＮＧ貯蔵タンク１１のガス接触側壁部の平均温度との差に比例する。そして、メタンガスが存在する空間部の高さは、ＬＮＧ貯蔵タンク１１の側壁の高さＨと液位ＬV との差で表され、接触側壁部の平均温度は、温度計１９の各検出値と液位検出計１７の液位検出値から判明するガス接触側壁部の範囲に存在する第１の温度センサ１８の各検出値の平均値ＴW となる。このため、メタン蒸発量ＧW （単位：ｔ／ｈ）は、
ＧW ＝Ｃ3 ・ＬV ・（Ｈ−ＬV ）・（ＴA −ＴW ）・・・・・（６）
と表せる。ここで、Ｃ3 はＬＮＧ貯蔵タンク１１の寸法、ＬＮＧ貯蔵タンク１１の側部の断熱構造の特性、及び熱伝達特性等を反映した設備定数である。
【００２１】
また、ＬＮＧ貯蔵タンク１１のガス接触天井部からの入熱は、ＬＮＧ貯蔵タンク１１の周囲の大気温度ＴA とＬＮＧ貯蔵タンク１１のガス接触天井部の最高温度との差に比例する。そして、ガス接触天井部の最高温度は、ガス接触側壁部の範囲に存在する第１の温度センサ１８の検出値の最高温度ＴC として求まる。このため、メタン蒸発量ＧC （単位：ｔ／ｈ）は、
ＧC ＝Ｃ4 ・（ＴA −ＴC ）・・・・・（７）
と表せる。ここで、Ｃ4 はＬＮＧ貯蔵タンク１１の天井部の寸法、天井部の断熱構造の特性、及び熱伝達特性等を反映した設備定数である。
【００２２】
また、サクションクーラー２２で熱交換される熱量は、ガス取り出し管２１を介してサクションクーラー２２に流入するまでに、ガス取り出し管２１の外側面を介して流入してメタンガスの加熱に使用された熱量に等価と考えられる。従って、各メタン蒸発量ＧR 、ＧY 、ＧM はいずれも配管の外側面からの入熱により生じると考えられるため、配管内のＬＮＧ温度は一定として各配管からのメタン蒸発量合計ＧA を大気温度の関数として表す実験式を構築した。
ＧA ＝ＣA ｅｘｐ（Ａ・ＴA ）・・・・・（８）
ここで、ＣA とＡは、各配管の寸法、断熱構造の特性、及び熱伝達特性を等を反映した設備定数である。
【００２３】
メタン強制蒸発量ＧF を除いた各メタン蒸発量は、それぞれ（４）〜（８）式により定量的に求められる。一方、図２に示すように、ＬＮＧ貯蔵設備１０におけるメタンガス処理量は理論メタン蒸発量Ｇに一致し、ＬＮＧ貯蔵タンク１１から外部に払い出された排出ガス量ＧE と、ＬＮＧ貯蔵タンク１１内のＬＮＧ液面降下代Ｄに相当する空間内を満たしている液面降下代吸収ガス量ＧP との総和
Ｇ＝ＧE ＋ＧP ・・・・・（９）
として表せる。
【００２４】
ここで、排出ガス量ＧE は、ガス圧縮機出口管２５に設けた流量計により測定することができる。また、液面降下代吸収ガス量ＧP は、ＬＮＧ液面降下代Ｄに相当する空間の体積に比例する。このため、液面降下代吸収ガス量ＧP （単位：ｔ／ｈ）は
ＧP ＝Ｃ5 ・Ｄ・・・・・（１０）
と表せる。なお、Ｃ5 は貯残ＬＮＧの高さから貯残ＬＮＧ重量を求める際の換算係数である。また、ＬＮＧ液面降下代Ｄは、液位検出計１７の検出値から得られる。
【００２５】
また、式（１）〜（３）より、
ＧF ＝ＧE ＋ＧP −（ＧB ＋ＧS ＋ＧW ＋ＧC ＋ＧR ＋ＧY ＋ＧM ）
・・・・・（１１）
が成立し、メタン強制蒸発量ＧF を定量的に決定することができ、理論メタン蒸発量Ｇを、ガス圧縮機出口管２５に設けた流量計による測定から得られる排出ガス量ＧE 、大気温度ＴA 、ＬＮＧ貯蔵タンク１１内に設置した液位検出計１７、温度計１９、及び第２の温度センサ２０の検出値から得られるＬＮＧ貯蔵タンク１１の底部の平均温度ＴB 、液接触側壁部の平均温度ＴS 、ガス接触側壁部の平均温度ＴW 、ガス接触天井部の最高温度ＴC 、及び液位ＬV を用いて算出することができる。
【００２６】
そして、ＬＮＧ貯蔵タンク１１に貯蔵しているＬＮＧのある時刻における組成を一度求めておくと、その時刻を基準としてＺ時間が経過したときまでの総メタン蒸発量は、式（４）〜（８）及び（１１）から得られる理論メタン蒸発量Ｇを用いてＧ・Ｚとして求まるので、Ｚ時間が経過したときのＬＮＧの組成を算出することができる。すなわち、基準時刻におけるＬＮＧの各成分のモル数が、例えば、メタンがＭＣ1 モル、エタンがＭＣ2 モル、プロパンがＭＣ3 モル、ブタンがＭＣ4 モル、ペンタンがＭＣ5 モル、及び窒素がＭＮ2 モルであり、Ｇ・Ｚトンに相当するメタンのモル数がＭBOG モルであるとすると、基準時刻からＺ時間経過後のＬＮＧの各成分のモル数は、メタンが（ＭＣ1 −ＭBOG ）モル、エタンがＭＣ2 モル、プロパンがＭＣ3 モル、ブタンがＭＣ4 モル、ペンタンがＭＣ5 モル、及び窒素がＭＮ2 モルとなる。
【００２７】
従って、基準時刻からＺ時間経過後のＬＮＧのモル分率で示した組成割合は、メタンモル分率ｘ1 が（ＭＣ1 −ＭBOG ）／（Ｍ0 −ＭBOG ）、エタンモル分率ｘ2 がＭＣ2 ／（Ｍ0 −ＭBOG ）、プロパンモル分率ｘ3 がＭＣ3 ／（Ｍ0 −ＭBOG ）、ブタンモル分率ｘ4 がＭＣ4 ／（Ｍ0 −ＭBOG ）、ペンタンモル分率ｘ5 がＭＣ5 ／（Ｍ0 −ＭBOG ）、窒素モル分率ｘN がＭＮ2 ／（Ｍ0 −ＭBOG ）となる。なお、Ｍ0 ＝ＭＣ1 ＋ＭＣ2 ＋ＭＣ3 ＋ＭＣ4 ＋ＭＣ5 ＋ＭＮ2 である。
そして、メタン、エタン、プロパン、ブタン、ペンタン、及び窒素の各分子量ｍ1 、ｍ2 、ｍ3 、ｍ4 、ｍ5 、ｍN を用いて、基準時刻からＺ時間経過後のＬＮＧの平均分子量ｍA を求めると、
ｍA ＝ｘ1 ・ｍ1 ＋ｘ2 ・ｍ2 ＋ｘ3 ・ｍ3 ＋ｘ4 ・ｍ4 ＋ｘ5 ・ｍ5 ＋ｘN ・ｍN 
・・・・・（１２）
となる。
【００２８】
以降の、ＬＮＧ密度計算は、第１回ＬＮＧ国際会議で提唱されたエアープロダクションアンドケミカル社のクローゼクとマッキンレーの方法による。
メタン、エタン、プロパン、ブタン、ペンタン、及び窒素の各モル容積ｖ1 、ｖ2 、ｖ3 、ｖ4 、ｖ5 、ｖN を用いて、基準時刻からＺ時間経過後のＬＮＧの平均モル容積ｖA を求めると、
ｖA ＝ｘ1 ・ｖ1 ＋ｘ2 ・ｖ2 ＋ｘ3 ・ｖ3 ＋ｘ4 ・ｖ4 ＋ｘ5 ・ｖ5 ＋ｘN ・ｖN 
・・・・・（１３）
となる。ここで、（１３）式の平均モル容積は、ＬＮＧを各組成の単なる混合状態と仮定した場合のモル容積を示したものである。
【００２９】
一方、実際のＬＮＧでは、各組成が混合すると容積減少が生じ、この容積減少量Ｓは、ＬＮＧの平均分子量ｍA 、主成分であるメタンモル分率ｘ1 、及びＬＮＧの温度に依存することが判明し、定量的に求まっている。このため、ＬＮＧの平均分子量ｍA 当たりの実平均モル容積ｖE は、
ｖE ＝ｖA −Ｓ・・・・・（１４）
と表され、従って、貯残ＬＮＧの密度Ｄは
Ｄ＝ｍA ／（ｖA −Ｓ）・・・・・（１５）
を用いて算出できる。
【実施例】
【００３０】
次に、本発明の作用効果を確認するために行った実施例について説明する。
ＬＮＧ貯蔵タンクにＬＮＧの受入れが完了した時点を基準時刻として、貯蔵しているＬＮＧの払い出しを行ないながら、ＬＮＧ貯蔵タンク内の各部位の温度の計測、液位検出計によるＬＮＧ貯蔵タンク内の貯残ＬＮＧの高さＬV の計測、及びＬＮＧ貯蔵設備の周囲の大気温度をそれぞれ検出すると共に、ＬＮＧ貯蔵タンクから払い出された排出ガス量ＧE を計測した。一方、一定期間毎にＬＮＧ貯蔵タンクからＬＮＧをサンプリングして、ＬＮＧの組成をガスクロマトグラフを用いて分析し、得られた組成に基づいてＬＮＧの密度を算出することで実測貯残ＬＮＧ密度ＤM を求めた。
【００３１】
次いで、実測貯残ＬＮＧ密度ＤM を求めるためにＬＮＧのサンプリングを行なった時刻に対応するＬＮＧ貯蔵タンク内の各部位の温度、貯残ＬＮＧの高さＬV 、及びＬＮＧ貯蔵設備の周囲の大気温度、及び排出ガス量ＧE を用いて式（４）〜（８）及び（１１）から理論メタン蒸発量Ｇを求め、式（１５）を用いて基準時刻以降の貯残ＬＮＧ密度Ｄを順次求めた。そして、貯残ＬＮＧ密度Ｄと実測貯残ＬＮＧ密度ＤM との差ｄｅｖを求め、差ｄｅｖと実測貯残ＬＮＧ密度ＤM の関係を求めた。その結果を図３に示す。図３から、差ｄｅｖは、実測貯残ＬＮＧ密度ＤM に対して±１ｋｇ／ｍ3 の範囲となり、貯残ＬＮＧ密度Ｄは実測貯残ＬＮＧ密度ＤM に対して±０．２％の範囲内で一致することが確認できた。
【００３２】
以上、本発明の実施の形態を説明したが、本発明は、この実施の形態に限定されるものではなく、発明の要旨を変更しない範囲での変更は可能であり、前記したそれぞれの実施の形態や変形例の一部又は全部を組み合わせて本発明のＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法を構成する場合も本発明の権利範囲に含まれる。
例えば、ＬＮＧ貯蔵タンク内のＬＮＧの液位をディスプレーサ式の液位検出計で測定したが、超音波等を使用した非接触式の液位検出計を使用することができる。また、複数の温度センサを用いた温度計を複数本ＬＮＧ貯蔵タンク内に設置して、より正確にＬＮＧ貯蔵タンク内の温度を測定するようにしてもよい。
【図面の簡単な説明】
【００３３】
【図１】本発明の一実施の形態に係るＬＮＧ貯蔵タンク内の貯残ＬＮＧ密度の算出方法を適用するＬＮＧ貯蔵設備の説明図である。
【図２】同ＬＮＧ貯蔵設備におけるメタンガス処理量とメタンガス発生量の関係を示す模式図である。
【図３】貯残ＬＮＧ密度と実測貯残ＬＮＧ密度との差に対する実測貯残ＬＮＧ密度の関係を示すグラフである。
【符号の説明】
【００３４】
１０：ＬＮＧ貯蔵設備、１１：ＬＮＧ貯蔵タンク、１２：ＬＮＧ受入れ配管、１３：ＬＮＧポンプ、１４：出荷設備配管、１５：ボイルオフガス配管系、１６：ディスプレーサ、１７：液位検出計、１８：第１の温度センサ、１９：温度計、２０：第２の温度センサ、２１：ガス取り出し管、２２：サクションクーラー、２３：連結管、２４：ガス圧縮機、２５：ガス圧縮機出口管、２６：ＬＮＧ供給配管、２７：ＬＮＧ受入れ側遮断弁
【図１】
2006308384.tif 000003
【図２】
2006308384.tif 000004
【図３】
2006308384.tif 000005