(57)【特許請求の範囲】
【請求項１】
液化ガスの組成を推定する組成推定装置であって、
組成が異なる種々の前記液化ガスについて、液密度又は熱量と組成とを実際に計測して得られるデータであって、前記液化ガスの液密度又は熱量の計測値と前記液化ガスの組成の計測値とが含まれる実績データを用いて前記液化ガスの液密度又は熱量と前記液化ガスの組成との相関関係を示す近似式を算出する近似式算出部と、
前記近似式算出部で算出された前記近似式を用い、前記液化ガスの液密度又は熱量の計測値から前記液化ガスの組成を推定する組成推定部と、
を備える組成推定装置。
【請求項２】
前記近似式算出部は、前記実績データのうち、前記近似式に対する距離が予め設定された基準距離以上のものを除外したものを用いて、新たな前記近似式を算出する、請求項１記載の組成推定装置。
【請求項３】
前記近似式算出部は、前記実績データを複数のグループに分け、前記グループ毎に前記近似式を算出する、請求項１又は請求項２記載の組成推定装置。
【請求項４】
前記近似式算出部は、前記液化ガスの液密度又は熱量の大きさに応じて前記実績データを複数の領域に分割し、分割された前記領域毎に前記近似式を算出する、請求項１から請求項３の何れか一項に記載の組成推定装置。
【請求項５】
前記液化ガスの液密度又は熱量の計測値と、液化ガスプラントを模擬するシミュレータによって得られる前記液化ガスの液密度又は熱量の推定値との差が予め規定された閾値を超えた場合に、前記組成推定部に前記液化ガスの組成を推定させる実行部を備える、請求項１から請求項４の何れか一項に記載の組成推定装置。
【請求項６】
前記実行部は、前記液化ガスの液密度又は熱量の計測値と、前記液化ガスの液密度又は熱量の推定値との差が予め規定された閾値を超えた場合に警告を発し、前記警告を発した後に前記液化ガスの組成を推定させる指示があったときに、前記組成推定部に前記液化ガスの組成を推定させる、請求項５記載の組成推定装置。
【請求項７】
液化ガスの組成を推定する組成推定方法であって、
組成が異なる種々の前記液化ガスについて、液密度又は熱量と組成とを実際に計測して得られるデータであって、前記液化ガスの液密度又は熱量の計測値と前記液化ガスの組成の計測値とが含まれる実績データを用いて前記液化ガスの液密度又は熱量と前記液化ガスの組成との相関関係を示す近似式を算出する近似式算出ステップと、
前記近似式算出ステップで算出された前記近似式を用い、前記液化ガスの液密度又は熱量の計測値から前記液化ガスの組成を推定する組成推定ステップと、
を有する組成推定方法。
【請求項８】
コンピュータを、液化ガスの組成を推定する組成推定装置として機能させる組成推定プログラムであって、
前記コンピュータを、
組成が異なる種々の前記液化ガスについて、液密度又は熱量と組成とを実際に計測して得られるデータであって、前記液化ガスの液密度又は熱量の計測値と前記液化ガスの組成の計測値とが含まれる実績データを用いて前記液化ガスの液密度又は熱量と前記液化ガスの組成との相関関係を示す近似式を算出する近似式算出手段と、
前記近似式算出手段で算出された前記近似式を用い、前記液化ガスの液密度又は熱量の計測値から前記液化ガスの組成を推定する組成推定手段と、
して機能させる組成推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、組成推定装置、組成推定方法、及び組成推定プログラムに関する。
【背景技術】
【０００２】
ＬＮＧ（Liquefied Natural Gas：液化天然ガス）は、メタンを主成分とし、エタン、プロパン、その他の炭化水素を含む混合物であり、その組成は産地によって異なっている。ＬＮＧの組成分析を行う場合には、例えば、液体クロマトグラフが用いられる。ここで、液体クロマトグラフは、分析及び精製装置の一種であり、物質の化学的相互作用や分子の大きさの違い等によって、混合成分を分離・定量する装置である。この液体クロマトグラフは、ＬＮＧの組成を求めることができるものの、分析に時間を要するため、即時性及び連続性が求められる用途には不向きである。
【０００３】
以下の特許文献１には、計測したＬＮＧの液密度を温度及び圧力の計測値を用いて補正し、予め求めておいた液密度と発熱量との相関関係から、ＬＮＧの発熱量を算出する方法が開示されている。また、以下の特許文献２には、ＬＮＧタンクのＢＯＧ（Boil off gas：自然入熱による気化ガス）量を推定し、ＢＯＧ量の分だけメタンが減少したとみなした場合のＬＮＧの組成割合を求めて、最終的にＬＮＧの液密度を算出する方法が開示されている。尚、以下の特許文献３には、現実のプラントを模擬するプラント・シミュレータの一種であるトラッキング・シミュレータが開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００６−４７０７１号公報
【特許文献２】特開２００６−３０８３８４号公報
【特許文献３】特開２００５−３３２３６０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、近年、コンピュータ上でプラントを再現するデジタルツインという考え方が盛んになっている。例えば、上述した特許文献３に開示されたトラッキング・シミュレータ等のプラント・シミュレータを実際のＬＮＧプラントに繋げ、ＬＮＧプラントを模擬することにより、オンラインでＬＮＧプラントの運転支援に役立てようとする試みがなされている。
【０００６】
ここで、プラント・シミュレータは、物理化学法則に基づいた計算によって成り立っているため、計算精度を保つためには、ＬＮＧプラントで扱われる物質の組成が正しく定義されている必要がある。しかしながら、前述した液体クロマトグラフは即時性に欠けるため、時々刻々変化する組成の変化を正確にプラント・シミュレータに取り込むことはできない。また、上述した特許文献１，２に開示された方法でも、時々刻々変化する組成の変化を正確に求めることはできない。
【０００７】
本発明は上記事情に鑑みてなされたものであり、時々刻々変化する液化ガスの組成を推定することができる組成推定装置、組成推定方法、及び組成推定プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するために、本発明の一態様による組成推定装置は、液化ガスの組成を推定する組成推定装置（１０）であって、前記液化ガスの液密度又は熱量の計測値と前記液化ガスの組成の計測値とが含まれる実績データを用いて前記液化ガスの液密度又は熱量と前記液化ガスの組成との相関関係を示す近似式を算出する近似式算出部（１３）と、前記近似式算出部で算出された前記近似式を用い、前記液化ガスの液密度又は熱量の計測値から前記液化ガスの組成を推定する組成推定部（１４）と、を備える。
【０００９】
また、本発明の一態様による組成推定装置は、前記近似式算出部が、前記実績データのうち、前記近似式に対する距離が予め設定された基準距離以上のものを除外したものを用いて、新たな前記近似式を算出する。
【００１０】
また、本発明の一態様による組成推定装置は、前記近似式算出部が、前記実績データを複数のグループに分け、前記グループ毎に前記近似式を算出する。
【００１１】
また、本発明の一態様による組成推定装置は、前記近似式算出部が、前記液化ガスの液密度又は熱量の大きさに応じて前記実績データを複数の領域に分割し、分割された前記領域毎に前記近似式を算出する。
【００１２】
また、本発明の一態様による組成推定装置は、前記液化ガスの液密度又は熱量の計測値と、液化ガスプラントを模擬するシミュレータによって得られる前記液化ガスの液密度又は熱量の推定値との差が予め規定された閾値を超えた場合に、前記組成推定部に前記液化ガスの組成を推定させる実行部（１２）を備える。
【００１３】
また、本発明の一態様による組成推定装置は、前記実行部が、前記液化ガスの液密度又は熱量の計測値と、前記液化ガスの液密度又は熱量の推定値との差が予め規定された閾値を超えた場合に警告を発し、前記警告を発した後に前記液化ガスの組成を推定させる指示があったときに、前記組成推定部に前記液化ガスの組成を推定させる。
【００１４】
本発明の一態様による組成推定方法は、液化ガスの組成を推定する組成推定方法であって、前記液化ガスの液密度又は熱量の計測値と前記液化ガスの組成の計測値とが含まれる実績データを用いて前記液化ガスの液密度又は熱量と前記液化ガスの組成との相関関係を示す近似式を算出する近似式算出ステップ（Ｓ１２、Ｓ１４）と、前記近似式算出ステップで算出された前記近似式を用い、前記液化ガスの液密度又は熱量の計測値から前記液化ガスの組成を推定する組成推定ステップ（Ｓ２４）と、を有する。
【００１５】
本発明の一態様による組成推定プログラムは、コンピュータを、液化ガスの組成を推定する組成推定装置（１０）として機能させる組成推定プログラムであって、前記コンピュータを、前記液化ガスの液密度又は熱量の計測値と前記液化ガスの組成の計測値とが含まれる実績データを用いて前記液化ガスの液密度又は熱量と前記液化ガスの組成との相関関係を示す近似式を算出する近似式算出手段（１３）と、前記近似式算出手段で算出された前記近似式を用い、前記液化ガスの液密度又は熱量の計測値から前記液化ガスの組成を推定する組成推定手段（１４）と、して機能させる。
【発明の効果】
【００１６】
本発明によれば、時々刻々変化する液化ガスの組成を推定することができるという効果がある。
【図面の簡単な説明】
【００１７】
【図１】本発明の一実施形態による組成推定装置の要部構成を示すブロック図である。
【図２】本発明の一実施形態における領域分割法を説明するための図である。
【図３】ＬＮＧに含まれる成分の一例を示す図である。
【図４】本発明の一実施形態による組成推定装置の近似式算出時の動作例を示すフローチャートである。
【図５】本発明の一実施形態による組成推定装置の組成推定時の動作例を示すフローチャートである。
【発明を実施するための形態】
【００１８】
以下、図面を参照して本発明の一実施形態による組成推定装置、組成推定方法、及び組成推定プログラムについて詳細に説明する。以下では、まず本発明の実施形態の概要について説明し、続いて本発明の実施形態の詳細について説明する。
【００１９】
〔概要〕
本発明の実施形態は、時々刻々変化する液化ガス（例えば、ＬＮＧ）の組成を推定することができるようにするものである。このように、液化ガスの組成を推定するのは、例えば、プラント・シミュレータによるオンラインでの液化ガスプラント（例えば、ＬＮＧプラント）の模擬を高い精度で行えるようにするためである。
【００２０】
従来、液化ガスプラントを模擬するプラント・シミュレータは、オフラインでの訓練等に使用されていた。このようなプラント・シミュレータでは、例えば、液体クロマトグラフ等によって測定された液化ガスの組成が入力されて用いられていた。ここで、プラント・シミュレータを実際の液化ガスプラントに繋げてオンラインで使用しようとした場合に、時々刻々変化する液化ガスの組成をプラント・シミュレータに取り込む必要がある。
【００２１】
しかしながら、前述した液体クロマトグラフは即時性に欠けるため、時々刻々変化する液化ガスの組成を取り込むことはできない。また、液体クロマトグラフは高価であるため、液化ガスを貯蔵するタンクの各々に液体クロマトグラフを設置すると飛躍的にコストが増大する。このため、現状では、タンクに入れる前の液化ガスの組成を液体クロマトグラフで計測し、タンク内に収容した後の液化ガスの組成変化をシミュレータで推定する方法が主流である。
【００２２】
但し、長時間に亘ってシミュレータを稼働すると、液化ガスの実際の組成とシミュレータで推定される液化ガスの組成とが、外乱の影響によって乖離することがある。尚、ここにいう外乱とは、タンクに貯蔵された液化ガスが送出されること、タンクに液化ガスが戻されること、外気温や圧力の変化、組成が異なる液化ガスの受け入れ（タンクに対する受け入れ）が幾度も行われること等である。
【００２３】
従って、プラント・シミュレータによるオンラインでの液化ガスプラントの模擬を高い精度で行うには、タンクに貯蔵されている液化ガスの組成と、プラント・シミュレータで推定される液化ガスの組成とのずれを補正する仕組みが必要になる。本発明の実施形態では、必ずといっていいほど計測されている液化ガスの液密度又は熱量を使用して、タンクに貯蔵された液化ガスの組成を推定するものである。ここで、液化ガスの液密度又は熱量と液化ガスの組成の相関関係は高いものも低いものもあり、液化ガスの産地によるばらつきもあるため、多量の実績データから液化ガスの液密度又は熱量と液化ガスの組成との相関関係を示す近似式を作成する。そして、上記の近似式を用いて、液化ガスの計測値から液化ガスの組成を推定するようにしている。
【００２４】
つまり、本発明の実施形態では、まず、液化ガスの液密度又は熱量の計測値と液化ガスの組成の計測値とが含まれる実績データを用いて液化ガスの液密度又は熱量と液化ガスの組成との相関関係を示す近似式を求める。次に、求められた近似式を用い、液化ガスの液密度又は熱量の計測値から液化ガスの組成を推定する。これにより、時々刻々変化する液化ガスの組成を推定することができる。
【００２５】
〔実施形態〕
〈組成推定装置の構成〉
図１は、本発明の一実施形態による組成推定装置の要部構成を示すブロック図である。図１に示す通り、本実施形態の組成推定装置１０は、ＬＮＧプラント１００を模擬するプラント・シミュレータ２００に設けられ、プラント・シミュレータ２００がＬＮＧプラント１００の模擬を高い精度で行うために必要となるＬＮＧの組成を推定する。
【００２６】
ここで、ＬＮＧプラント１００は、例えば、ガス田で採取された天然ガスから不純物、環境汚染物質、水分等を除去し、不純物等が除去された天然ガスの冷却を行って液化することによりＬＮＧを得るものである。ＬＮＧプラント１００には、ＬＮＧを貯蔵するＬＮＧタンク１０１が設けられている。ＬＮＧタンク１０１には、ＬＮＧタンク１０１に貯蔵されているＬＮＧの液密度を計測する液密度計１０２が設置されている。
【００２７】
プラント・シミュレータ２００は、ネットワーク（図示省略）を介してＬＮＧプラント１００に接続されており、ＬＮＧプラント１００で得られた測定データや制御データ等を取得してＬＮＧプラント１００を模擬する。プラント・シミュレータ２００は、組成推定装置１０及びシミュレータモデル２０を備える。
【００２８】
シミュレータモデル２０は、模擬の対象であるＬＮＧプラント１００をモデル化したものである。このシミュレータモデル２０には、ＬＮＧタンク１０１をモデル化したＬＮＧタンクモデル２１、液密度計１０２をモデル化した液密度計モデル２２が設けられている。ＬＮＧプラント１００から取得された測定データや制御データ等がシミュレータモデル２０に入力されることで、ＬＮＧプラント１００が模擬される。
【００２９】
図１に示す通り、組成推定装置１０は、データ入力部１１、実行部１２、近似式算出部１３（近似式算出手段）、及び組成推定部１４（組成推定手段）を備える。
【００３０】
データ入力部１１は、ＬＮＧの液密度の計測値とＬＮＧの組成の計測値とが含まれる実績データを入力する。上記の実績データは、産地、組成等が異なる種々のＬＮＧについて液密度と組成とを実際に計測して得られたデータである。ＬＮＧの液密度は、例えば、液密度計を用いて計測され、ＬＮＧの組成は、例えば、液体クロマトグラフを用いて計測される。尚、ＬＮＧの液密度は、ＬＮＧタンク１０１に設置された液密度計１０２を用いて計測されても良い。
【００３１】
実行部１２は、組成推定部１４に対して、ＬＮＧの推定される組成を求める演算（以下、「組成推定演算」という）を実行させる。具体的に、実行部１２は、ＬＮＧタンク１０１に設置された液密度計１０２で計測された液密度（以下、「実測値」という）と、ＬＮＧタンクモデル２１に設置された液密度計モデル２２で推定された液密度（以下、「推定値」という）とを比較する。そして、実行部１２は、液密度の実測値と推定値との差が予め規定された閾値を超えた場合に、組成推定演算を組成推定部１４に実行させる。
【００３２】
ここで、実行部１２は、組成推定演算を組成推定部１４に実行させる前に、液密度の実測値と推定値との差が予め規定された閾値を超えた旨を示す警告を発するようにしても良い。例えば、実行部１２は、上記の警告を表示装置３００に表示するようにしてもよい。そして、実行部１２は、上記の警告を発した後に、ＬＮＧプラント１００の運転員からＬＮＧの組成を推定させる指示があったときに、組成推定演算を組成推定部１４に実行させるようにしても良い。或いは、実行部１２は、上記の警告を発した後に、自動的に、組成推定演算を組成推定部１４に実行させるようにしても良い。
【００３３】
近似式算出部１３は、データ入力部１１から入力される実績データを用いて、ＬＮＧの液密度と組成との相関関係を示す近似式を算出して保持する。近似式算出部１３は、例えば、プラント・シミュレータ２００が動作する前（ＬＮＧプラント１００の模擬を行う前）に上記の近似式を算出して保持する。近似式算出部１３は、プラント・シミュレータ２００の動作中にデータ入力部１１から入力される実績データを用いて新たな近似式を算出し、保持している近似式を新たに算出された近似式に更新するようにしても良い。
【００３４】
近似式算出部１３は、近似式を算出した後に、異なる方法で近似式を再度算出することが可能である。近似式算出部１３は、例えば、以下の３つの方法で近似式を再度算出することが可能である。近似式算出部１３に近似式を再度算出させるか否かの指示、及び近似式算出部１３に近似式を再度算出させる際に用いる方法の指定は、近似式算出部１３に対するパラメータ設定により行う。
・外れ値除外法
・グループ化法
・領域分割法
【００３５】
「外れ値除外法」は、実績データのうち、既に算出されている近似式に対する距離が予め設定された基準距離以上のものを除外したものを用いて、新たな近似式を算出する方法である。この方法で近似式を算出することで、近似式の精度を向上させることができる。
【００３６】
「グループ化法」は、実績データを複数のグループに分け、グループ毎に近似式を算出する方法である。ＬＮＧは、産地毎の特性があり、産地毎に組成の配分がある程度偏っているため、実績データをグループ分けすることにより、近似式の精度を向上させることができる。
【００３７】
「領域分割法」は、ＬＮＧの液密度に応じて実績データを複数の領域に分割し、分割された領域毎に近似式を算出する方法である。図２は、本発明の一実施形態における領域分割法を説明するための図である。尚、図２に示すグラフは、横軸に液密度［ｋｇ／ｍ3］をとり、縦軸に熱量［ＭＪ］をとってある。
【００３８】
図２に示す通り、液密度が４２０〜４７０［ｋｇ／ｍ3］の範囲のデータが、実績データとしてデータ入力部１１から入力されるとする。領域分割法では、このような実績データを、実績データの分布状況に応じて複数の領域に分割する。図２に示す例では、例えば、４２０〜４３５［ｋｇ／ｍ3］、４３５〜４４６［ｋｇ／ｍ3］、４４６〜４７０［ｋｇ／ｍ3］の３つの領域に分割する。そして、分割した３つの領域毎に近似式を算出する。このように、分割された３つの領域毎に近似式を算出することで、全領域で１つの近似式を算出する場合よりも、近似式の精度を向上させることができる。尚、領域の分割数は３に制限されることはなく、任意の数に分割することが可能である。
【００３９】
ここで、ＬＮＧは、メタン（ＣＨ4）、エタン（Ｃ2Ｈ6）、プロパン（Ｃ3Ｈ8）、ブタン（Ｃ4Ｈ10）、ペンタン（Ｃ5Ｈ12）、窒素（Ｎ2）等を含んでおり、図３に示す通り、メタンが８５％以上占めるのが特徴である。図３は、ＬＮＧに含まれる成分の一例を示す図である。図３（ａ）は、ＬＮＧに含まれるメタンの一例を示す図であり、図３（ｂ）は、ＬＮＧに含まれるエタンを示す図である。尚、図３に示すグラフは、横軸に液密度［ｋｇ／ｍ3］をとり、縦軸に組成割合［％］をとってある。
【００４０】
図３（ａ）に示す通り、ＬＮＧに含まれるメタンの組成割合は、ＬＮＧの液密度が低くなるにつれて大きくなり、ＬＮＧの液密度が高くなるにつれて小さくなる傾向がある。これに対し、図３（ｂ）に示す通り、ＬＮＧに含まれるエタンの組成割合は、ＬＮＧの液密度が高くなるにつれて大きくなり、ＬＮＧの液密度が低くなるにつれて小さくなる傾向がある。
【００４１】
図３（ａ）に示す例において、ＬＮＧの液密度が４２０［ｋｇ／ｍ3］近辺である場合には、ＬＮＧに含まれるメタンの組成割合は９９［％］以上になる。このため、ＬＮＧの組成は、ほぼメタンであると特定できており、その熱量は、図２から、３９．８［ＭＪ／ｍ3］程度であると推定される。
【００４２】
つまり、図２，３を参照すると、ＬＮＧの液密度が低いほどメタンの組成割合が大きくなって他の組成の割合が小さくなり、ＬＮＧの液密度が高いほどメタンの組成割合が小さくなって他の組成の割合が大きくなる。言い換えると、ＬＮＧの液密度が高くなるほどＬＮＧの組成の推定精度が低くなってしまうため、統計処理により妥当な組成割合を計算する必要がある。
【００４３】
尚、上述した、外れ値除外法、グループ化法、及び領域分割法のうちの任意の２つの方法又は全ての方法を全て用いて新たな近似式を算出することも可能である。例えば、全ての方法を用いる場合には、実績データのうち、既に算出されている近似式に対する距離が基準距離以上のものを除外し（外れ値除外法）、産地毎に実績データをグループ分けし（グループ化法）、液密度に応じて実績データを複数の領域に分割し、各々の領域毎に近似式を算出する（領域分割法）。
【００４４】
組成推定部１４は、実行部１２の指示に基づいて、前述した組成推定演算を実行する。具体的に、組成推定部１４は、近似式算出部１３で算出された近似式を用い、ＬＮＧタンク１０１に設置された液密度計１０２で計測された液密度（実測値）に合うようにシミュレータモデル２０の液密度を変更したときのＬＮＧの組成を推定する。尚、組成推定部１４で推定されたＬＮＧの組成は、シミュレータモデル２０に入力される。ＬＮＧの組成が入力されたシミュレータモデル２０は、プラント・シミュレータ２００がＬＮＧプラント１００を模擬するために用いられる。
【００４５】
このような組成推定装置１０は、例えば、パーソナルコンピュータやワークステーション等のコンピュータにより実現される。組成推定装置１０がコンピュータにより実現される場合において、組成推定装置１０に設けられる各ブロック（データ入力部１１、実行部１２、近似式算出部１３、及び組成推定部１４）は、各々の機能を実現するためのプログラムが、コンピュータに設けられたＣＰＵ（中央処理装置）で実行されることによって実現される。つまり、組成推定装置１０に設けられる各ブロックは、ソフトウェアとハードウェア資源とが協働することによって実現される。尚、組成推定装置１０は、コンピュータによって実現されるプラント・シミュレータ２００の一機能として実現されても良い。
【００４６】
〈組成推定装置の動作〉
次に、組成推定装置の動作について説明する。組成推定装置１０の動作は、前述した近似式を算出する際の動作と、ＬＮＧの組成を推定する際の動作とに大別される。以下では、前述した近似式算出時の動作と、ＬＮＧの組成推定時の動作とを順に説明する。
【００４７】
《近似式算出時の動作》
図４は、本発明の一実施形態による組成推定装置の近似式算出時の動作例を示すフローチャートである。尚、図４に示すフローチャートの処理は、例えば、近似式の算出に用いられる実績データの準備が完了した後に、ＬＮＧプラント１００の運転員から近似式の作成開始指示がなされることによって開始される。
【００４８】
処理が開始されると、まず、実績データを組成推定装置１０に入力する処理が、データ入力部１１によって行われる（ステップＳ１１）。次に、組成推定装置１０に入力された実績データを用いて、予め入力されたパラメータ設定によるＬＮＧの液密度とＬＮＧの組成との相関関係を示す近似式を作成する処理が近似式算出部１３によって行われる（ステップＳ１２：近似式算出ステップ）。
【００４９】
次いで、近似式算出部１３に対するパラメータ設定変更がなされているか否かが近似式算出部１３によって判断される（ステップＳ１３）。つまり、近似式算出部１３に近似式を再度算出させるか否かの指示、及び近似式算出部１３に近似式を再度算出させる際に用いる方法の指定がなされているか否かが近似式算出部１３によって判断される。
【００５０】
近似式算出部１３に対するパラメータ設定変更がなされていると近似式算出部１３が判断した場合（ステップＳ１３の判断結果が「ＹＥＳ」の場合）には、実績データを用いて近似式を作成する処理が再び近似式算出部１３によって行われる（ステップＳ１２）。但し、ここでは、パラメータ設定変更によって指示された方法（前述した、外れ値除外法、グループ化法、及び領域分割法の少なくとも１つ）を用いた近似式の作成が行われる。
【００５１】
これに対し、近似式算出部１３に対するパラメータ設定変更がなされていないと近似式算出部１３が判断した場合（ステップＳ１３の判断結果が「ＮＯ」の場合）には、作成した近似式を保持する処理が近似式算出部１３によって行われる（ステップＳ１４：近似式算出ステップ）。以上の処理が終了すると、図４に示す一連の処理が終了する。このようにして近似式の作成が行われる。
【００５２】
《組成推定時の動作》
図５は、本発明の一実施形態による組成推定装置の組成推定時の動作例を示すフローチャートである。尚、図５に示すフローチャートの処理は、例えば、一定の時間間隔で繰り返し行われる。
【００５３】
処理が開始されると、まず、ＬＮＧタンク１０１に設置された液密度計１０２で計測された液密度の実測値と、ＬＮＧタンクモデル２１に設置された液密度計モデル２２で推定された液密度の推定値とを取得する処理が実行部１２によって行われる（ステップＳ２１）。次に、取得した液密度の実測値と推定値との差を求め、この差が予め規定された閾値を超えたか否かを判断する処理が実行部１２によって行われる（ステップＳ２２）。
【００５４】
液密度の実測値と推定値との差が閾値を超えていないと実行部１２が判断した場合（ステップＳ２２の判断結果が「ＮＯ」の場合）には、図５に示す一連の処理が終了する。これに対し、液密度の実測値と推定値との差が閾値を超えたと実行部１２が判断した場合（ステップＳ２２の判断結果が「ＹＥＳ」の場合）には、その旨を示す警告を発する処理が実行部１２で行われる。例えば、警告を表示装置３００に表示させる処理が実行部１２で行われる（ステップＳ２３）。
【００５５】
次いで、組成推定演算を組成推定部１４に実行させる処理が実行部１２によって行われる。ここで、実行部１２は、上記の警告を発した後に、ＬＮＧプラント１００の運転員からＬＮＧの組成を推定させる指示があったときに、組成推定演算を組成推定部１４に実行させるようにしても良い。或いは、実行部１２は、上記の警告を発した後に、自動的に、組成推定演算を組成推定部１４に実行させるようにしても良い。
【００５６】
このとき、ＬＮＧタンク１０１に貯蔵されているＬＮＧの産地が特定できるのであれば、ＬＮＧプラント１００の運転員は、組成推定装置１０に対してＬＮＧの産地を特定するようにしても良い。例えば、ＬＮＧプラント１００の運転員が、予め複数用意された産地の候補から特定の産地を選択することによってＬＮＧの産地を特定するようにしても良い。このような産地の特定を行うことで、組成推定部１４で組成推定演算が実行されるときに、産地に適した近似式を用いることで、高い精度でＬＮＧの組成を推定することが可能になる。
【００５７】
続いて、組成推定演算が組成推定部１４で実行される（ステップＳ２４：組成推定ステップ）。具体的には、近似式算出部１３で算出された近似式（図４のステップＳ１４で保持された近似式）を用い、ＬＮＧタンク１０１に設置された液密度計１０２で計測された液密度の実測値に合うようにシミュレータモデル２０の液密度を変更したときのＬＮＧの組成を推定する処理が行われる。
【００５８】
ここで、近似式算出部１３で算出された近似式を用いて推定されたＬＮＧの組成が１００％にならない場合がある。このような場合には、推定されたＬＮＧの組成の合計が１００％になるように調整する処理が組成推定部１４によって行われる。この処理では、例えば、各組成の上下限を設定して組成の合計が１００％になるように調整する方法や、各組成に重み付けを行って組成の合計が１００％になるように調整する方法が用いられる。
【００５９】
前者の方法では、例えば、メタンについて上下限値として８５〜９９．９％が設定される。そして、推定されたメタンの組成が設定された上下限値を逸脱している場合には、メタンの組成が上下限値に収まり、且つ組成の合計が１００％になるようにメタンの組成が調整される。
【００６０】
後者の方法では、例えば、メタンについては１００％の重みが設定され、他の成分については０％の重みが設定されたとする。このときには、他の成分の調整を行うことなく、組成の合計が１００％になるようにメタンの組成のみが調整される。また、例えば、メタンについては８０％の重みが設定され、エタンについては２０％の重みが設定されている場合において、近似式算出部１３で算出された近似式を用いて推定されたＬＮＧの組成が１０５％になったとする。このときには、メタンの組成を４％減ずるとともにエタンの組成を１％減ずることによって、組成の合計が１００％になるように調整される。
【００６１】
以上の処理が終了すると、推定された組成を出力して、シミュレータモデル２０に入力させる処理が組成推定部１４によって行われる（ステップＳ２５）。これにより、シミュレータモデル２０に新たに入力されたＬＮＧの組成を用いたＬＮＧプラント１００の模擬が、プラント・シミュレータ２００で行われる。
【００６２】
以上の通り、本実施形態では、近似式算出部１３が、ＬＮＧの液密度の計測値とＬＮＧの組成の計測値とが含まれる実績データを用いてＬＮＧの液密度と組成との相関関係を示す近似式を算出し、組成推定部１４が、算出された近似式を用い、ＬＮＧの液密度の計測値からＬＮＧの組成を推定するようにしている。これにより、時々刻々変化するＬＮＧの組成を推定することができる。
【００６３】
ここで、本実施形態では、ＬＮＧの液密度とＬＮＧの組成との相関関係を示す近似式を用いて、ＬＮＧの液密度の計測値からＬＮＧの組成を推定するようにしている。このため、液密度計１０２が設置されたＬＮＧタンク１０１に貯蔵されたＬＮＧの組成はもちろんのこと、ＬＮＧプラント内において液密度を計測可能な任意の箇所におけるＬＮＧの組成を推定することが可能である。
【００６４】
以上、本発明の一実施形態による組成推定装置、組成推定方法、及び組成推定プログラムについて説明したが、本発明は上記実施形態に制限される訳ではなく、本発明の範囲内で自由に変更が可能である。例えば、ＬＮＧタンク１０１に受け入れる前のＬＮＧの組成を液体クロマトグラフ等で計測できた場合（例えば、船上でＬＮＧの組成を計測した場合）に、得られた組成を組成推定装置１０に入力してＬＮＧの組成の推定精度を高めるようにしても良い。
【００６５】
また、上述した実施形態では、ＬＮＧの液密度とＬＮＧの組成との相関関係を示す近似式を算出し、算出した近似式を用いてＬＮＧの液密度の計測値からＬＮＧの組成を推定するようにしていた。しかしながら、一般的に、液密度と熱量とは相関が高いことが知られている。このため、ＬＮＧの熱量とＬＮＧの組成との相関関係を示す近似式を算出し、算出した近似式を用いてＬＮＧの熱量の計測値からＬＮＧの組成を推定するようにしても良い。
【００６６】
また、上述した実施形態では、ＬＮＧプラント１００のＬＮＧタンク１０１に貯蔵されたＬＮＧの組成を推定する例について説明した。しかしながら、本発明は、ＬＮＧの組成の推定に制限される訳ではなく、炭化水素を主成分とする液化ガスの組成の推定に適用することができる。このような液化ガスとしては、ＬＮＧ以外に、例えば、ＬＰＧ（Liquefied Petroleum Gas：液化石油ガス）が挙げられる。
【符号の説明】
【００６７】
１０    組成推定装置
１２    実行部
１３    近似式算出部
１４    組成推定部
【要約】
【課題】時々刻々変化する液化ガスの組成を推定することができる組成推定装置、組成推定方法、及び組成推定プログラムを提供する。
【解決手段】組成推定装置は、液化ガスの液密度又は熱量の計測値と液化ガスの組成の計測値とが含まれる実績データを用いて液化ガスの液密度又は熱量と液化ガスの組成との相関関係を示す近似式を算出する近似式算出部と、近似式算出部で算出された近似式を用い、液化ガスの液密度又は熱量の計測値から液化ガスの組成を推定する組成推定部と、を備える。
【選択図】図１【図１】
0006936378.tif 000002
【図２】
0006936378.tif 000003
【図３】
0006936378.tif 000004
【図４】
0006936378.tif 000005
【図５】
0006936378.tif 000006