(57)【要約】【課題】任意組成の多成分系液化ガスから気化した気化ガスの組成を定量的且つ容易に推算する気体組成推算プログラム及び気体組成推算装置を提供する。【解決手段】本発明の気体組成推算装置は、入力手段，数値仮定手段，モル分率算出手段，モル密度算出手段，フガシティー算出手段，分配率算出手段，条件判定手段，繰り返し手段，及び出力手段を備えており、多成分系の液化ガスを所定温度下で所定容積のボンベに仕込んだ場合の、ボンベ内の気相の組成を推算する装置である。気体組成推算装置に温度Ｔ，密閉空間の容積Ｖ，液化ガスの仕込物質量Ｆ，及び液化ガスを構成する各成分のモル分率ｚを入力すると、圧力Ｐ，物質量比ξ，分配率Ｋが仮定され、気相及び液相の各成分の組成ｘi ，ｙi 及び圧力Ｐが試行錯誤法により算出され、出力される。【選択図】図１３
2011126965.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも１種類の炭化水素化合物を含有する多成分系の液化ガスを、所定温度下で所定容積の密閉空間内に仕込んだ場合の、前記密閉空間内の気相の組成を推算するプログラムであって、
前記温度Ｔ，前記密閉空間の容積Ｖ，前記液化ガスの仕込物質量Ｆ，及び仕込んだ前記液化ガスを構成する各成分のモル分率ｚの入力を待つ入力ステップと、
前記入力ステップで入力された前記温度Ｔ，前記密閉空間の容積Ｖ，前記仕込物質量Ｆ，及び前記各成分のモル分率ｚに基づいて前記密閉空間内の圧力Ｐを仮定するとともに、前記密閉空間内の気相の物質量Ｄと前記仕込物質量Ｆとの比Ｄ／Ｆ、及び、前記密閉空間内の気相の各成分のモル分率ｙと液相の各成分のモル分率ｘとの比ｙ／ｘを任意の数値に仮定する数値仮定ステップと、
前記数値仮定ステップで仮定した前記比Ｄ／Ｆを物質量比ξ、前記比ｙ／ｘを分配率Ｋとし、式ｘ＝ｚ／｛１＋（Ｋ−１）ξ｝を用いて前記各成分の前記モル分率ｘをそれぞれ算出するとともに、算出した前記モル分率ｘを式Ｋ＝ｙ／ｘに代入して前記各成分の前記モル分率ｙをそれぞれ算出するモル分率算出ステップと、
前記モル分率算出ステップで算出した前記モル分率ｘ，ｙを、Benedict-Webb-Rubin 状態方程式に代入して、前記密閉空間内の液相のモル密度ρL 及び気相のモル密度ρV をそれぞれ算出するモル密度算出ステップと、
前記密閉空間内の液相のフガシティーをｆL 、気相のフガシティーをｆV とし、下記式１（ただし、該式１中のｎは物質量、Ｚは圧縮係数であり、ｉは成分の種類を表す）を用いて前記各成分のｆL ／（ｘＰ）、ｆV ／（ｙＰ）をそれぞれ算出するフガシティー算出ステップと、
【数１】
2011126965.tif 000035
前記フガシティー算出ステップで算出したｆL ／（ｘＰ）及びｆV ／（ｙＰ）から前記各成分の｛ｆL ／（ｘＰ）｝／｛ｆV ／（ｙＰ）｝を算出し、Ｋ’とする分配率算出ステップと、
前記分配率算出ステップで算出した前記各成分のＫ’が前記分配率Ｋとそれぞれ等しいという第一条件、前記モル分率算出ステップで算出した前記各成分の前記モル分率ｘの総和が１に等しいという第二条件、前記モル分率算出ステップで算出した前記各成分の前記モル分率ｙの総和が１に等しいという第三条件、並びに、前記モル密度算出ステップで算出した前記モル密度ρL ，ρV が下記式２を満足するという第四条件がそれぞれ満たされているかを判定する条件判定ステップと、
【数２】
2011126965.tif 000036
前記条件判定ステップにおいて前記第一条件から前記第四条件までの４つの条件のうち少なくとも１つが満たされていなかった場合は、前記数値仮定ステップに戻って前記圧力Ｐ，前記物質量比ξ，及び前記分配率Ｋを別の数値に変更し、前記モル分率算出ステップ，前記モル密度算出ステップ，前記フガシティー算出ステップ，及び前記分配率算出ステップの計算を再度行い、前記４つの条件が全て満たされるまで前記数値仮定ステップ，前記モル分率算出ステップ，前記モル密度算出ステップ，前記フガシティー算出ステップ，及び前記分配率算出ステップの操作を繰り返す繰り返しステップと、
前記条件判定ステップにおいて前記４つの条件が全て満たされた場合は、前記モル分率算出ステップで算出した前記各成分の前記モル分率ｘをそれぞれ出力する出力ステップと、
を有する処理をコンピュータに実行させることを特徴とする気体組成推算プログラム。
【請求項２】
少なくとも１種類の炭化水素化合物を含有する多成分系の液化ガスを、所定温度下で所定容積の密閉空間内に仕込んだ場合の、前記密閉空間内の気相の組成を推算する装置であって、
前記温度Ｔ，前記密閉空間の容積Ｖ，前記液化ガスの仕込物質量Ｆ，及び仕込んだ前記液化ガスを構成する各成分のモル分率ｚを入力する入力手段と、
前記入力手段により入力された前記温度Ｔ，前記密閉空間の容積Ｖ，前記仕込物質量Ｆ，及び前記各成分のモル分率ｚに基づいて前記密閉空間内の圧力Ｐを仮定するとともに、前記密閉空間内の気相の物質量Ｄと前記仕込物質量Ｆとの比Ｄ／Ｆ、及び、前記密閉空間内の気相の各成分のモル分率ｙと液相の各成分のモル分率ｘとの比ｙ／ｘを任意の数値に仮定する数値仮定手段と、
前記数値仮定手段により仮定した前記比Ｄ／Ｆを物質量比ξ、前記比ｙ／ｘを分配率Ｋとし、式ｘ＝ｚ／｛１＋（Ｋ−１）ξ｝を用いて前記各成分の前記モル分率ｘをそれぞれ算出するとともに、算出した前記モル分率ｘを式Ｋ＝ｙ／ｘに代入して前記各成分の前記モル分率ｙをそれぞれ算出するモル分率算出手段と、
前記モル分率算出手段により算出した前記モル分率ｘ，ｙを、Benedict-Webb-Rubin 状態方程式に代入して、前記密閉空間内の液相のモル密度ρL 及び気相のモル密度ρV をそれぞれ算出するモル密度算出手段と、
前記密閉空間内の液相のフガシティーをｆL 、気相のフガシティーをｆV とし、下記式１（ただし、該式１中のｎは物質量、Ｚは圧縮係数であり、ｉは成分の種類を表す）を用いて前記各成分のｆL ／（ｘＰ）、ｆV ／（ｙＰ）をそれぞれ算出するフガシティー算出手段と、
【数３】
2011126965.tif 000037
前記フガシティー算出手段により算出したｆL ／（ｘＰ）及びｆV ／（ｙＰ）から前記各成分の｛ｆL ／（ｘＰ）｝／｛ｆV ／（ｙＰ）｝を算出し、Ｋ’とする分配率算出手段と、
前記分配率算出手段により算出した前記各成分のＫ’が前記分配率Ｋとそれぞれ等しいという第一条件、前記モル分率算出手段により算出した前記各成分の前記モル分率ｘの総和が１に等しいという第二条件、前記モル分率算出手段により算出した前記各成分の前記モル分率ｙの総和が１に等しいという第三条件、並びに、前記モル密度算出手段により算出した前記モル密度ρL ，ρV が下記式２を満足するという第四条件がそれぞれ満たされているかを判定する条件判定手段と、
【数４】
2011126965.tif 000038
前記条件判定手段による判定で前記第一条件から前記第四条件までの４つの条件のうち少なくとも１つが満たされていなかった場合は、前記数値仮定手段により前記圧力Ｐ，前記物質量比ξ，及び前記分配率Ｋを別の数値に変更し、前記モル分率算出手段，前記モル密度算出手段，前記フガシティー算出手段，及び前記分配率算出手段による計算を再度行い、前記４つの条件が全て満たされるまで前記数値仮定手段，前記モル分率算出手段，前記モル密度算出手段，前記フガシティー算出手段，及び前記分配率算出手段による操作を繰り返す繰り返し手段と、
前記条件判定手段による判定で前記４つの条件が全て満たされた場合は、前記モル分率算出手段により算出した前記各成分の前記モル分率ｘをそれぞれ出力する出力手段と、
を備えることを特徴とする気体組成推算装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、任意組成の多成分系液化ガスから気化した気化ガスの組成を定量的に推算するプログラム及び装置に関する。
【背景技術】
【０００２】
プロパン，ブタン等の炭化水素化合物を主成分とする液化石油ガス（ＬＰＧ）と他の燃料成分とを混合した液化混合燃料ガスをボンベに充填した場合は、揮発性の高い成分を多く含む気化ガスが生成してボンベ外に供給されるため、ボンベからの気化ガスの消費に伴って気化ガスの組成が徐々に変化することとなる。よって、上記のような液化混合燃料ガスから発生した気化ガスを燃焼させるガス燃焼器（例えば内燃機関，ボイラー，ガスコンロ）は、気化ガスの組成の変化に応じて燃焼条件の設定を自動的に最適化する機構や、気化ガスの組成が変化しても問題なく気化ガスを燃焼させることが可能な構造を備えている必要がある。
【０００３】
前記のような機構や構造を備えるガス燃焼器を開発する場合には、気化ガスの組成がどのように変化するのか、そのデータが予め必要である。例えば、加藤らや大井らは、１９７０年代に、ＬＰＧボンベからの自然気化モデルとして、ＬＰＧボンベの内圧が７０ｋＰａに到達する時間を外気温，風速，及び湿度の影響を考慮して予測する方法を提案している。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Benedict et al., J. Chem. Phys., 8, 334(1940)
【非特許文献２】Cooper et al., Hydrocarbon Processing, 46, 142(1967)
【非特許文献３】Tsuji et al., Proceeding of AIChE Annual Meeting, San Francisco (2006)
【非特許文献４】Joffe, Chemical Engineering Process, 45,160(1949)
【非特許文献５】Reid et al., Properties of Gases & Liquids 4th Ed., McGraw Hill(1986)
【非特許文献６】Younglove et al., J. Phys. Chem. Ref. Data, 16, 577(1987)
【非特許文献７】Stotler et al., Chem. Eng. Progr. Symp. Ser. No. 6, 49, 25(1953)
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、液化混合燃料ガスから自然気化して発生した気化ガスの組成を推算する方法は従来知られていたものの、共沸組成のような気液組成が合致するようなごく限られた場合にしか適用できないものであり、任意組成の多成分系液化混合燃料ガスから自然気化して発生した気化ガスの組成を、定量的且つ容易に推算する方法は知られていなかった。そのため、前記のような機構や構造を備えるガス燃焼器を開発する場合には、気化ガスのあらゆる組成を想定して実験する必要があるため、多大な時間と労力を要していた。
そこで、本発明は上記のような従来技術が有する問題点を解決し、任意組成の多成分系液化ガスから気化した気化ガスの組成を定量的且つ容易に推算するプログラム及び装置を提供することを課題とする。
【課題を解決するための手段】
【０００６】
前記課題を解決するため、本発明は次のような構成からなる。すなわち、本発明に係る気体組成推算プログラムは、少なくとも１種類の炭化水素化合物を含有する多成分系の液化ガスを、所定温度下で所定容積の密閉空間内に仕込んだ場合の、前記密閉空間内の気相の組成を推算するプログラムであって、前記温度Ｔ，前記密閉空間の容積Ｖ，前記液化ガスの仕込物質量Ｆ，及び仕込んだ前記液化ガスを構成する各成分のモル分率ｚの入力を待つ入力ステップと、前記入力ステップで入力された前記温度Ｔ，前記密閉空間の容積Ｖ，前記仕込物質量Ｆ，及び前記各成分のモル分率ｚに基づいて前記密閉空間内の圧力Ｐを仮定するとともに、前記密閉空間内の気相の物質量Ｄと前記仕込物質量Ｆとの比Ｄ／Ｆ、及び、前記密閉空間内の気相の各成分のモル分率ｙと液相の各成分のモル分率ｘとの比ｙ／ｘを任意の数値に仮定する数値仮定ステップと、前記数値仮定ステップで仮定した前記比Ｄ／Ｆを物質量比ξ、前記比ｙ／ｘを分配率Ｋとし、式ｘ＝ｚ／｛１＋（Ｋ−１）ξ｝を用いて前記各成分の前記モル分率ｘをそれぞれ算出するとともに、算出した前記モル分率ｘを式Ｋ＝ｙ／ｘに代入して前記各成分の前記モル分率ｙをそれぞれ算出するモル分率算出ステップと、前記モル分率算出ステップで算出した前記モル分率ｘ，ｙを、Benedict-Webb-Rubin 状態方程式に代入して、前記密閉空間内の液相のモル密度ρL 及び気相のモル密度ρV をそれぞれ算出するモル密度算出ステップと、前記密閉空間内の液相のフガシティーをｆL 、気相のフガシティーをｆV とし、下記式１（ただし、該式１中のｎは物質量、Ｚは圧縮係数であり、ｉは成分の種類を表す）を用いて前記各成分のｆL ／（ｘＰ）、ｆV ／（ｙＰ）をそれぞれ算出するフガシティー算出ステップと、前記フガシティー算出ステップで算出したｆL ／（ｘＰ）及びｆV ／（ｙＰ）から前記各成分の｛ｆL ／（ｘＰ）｝／｛ｆV ／（ｙＰ）｝を算出し、Ｋ’とする分配率算出ステップと、前記分配率算出ステップで算出した前記各成分のＫ’が前記分配率Ｋとそれぞれ等しいという第一条件、前記モル分率算出ステップで算出した前記各成分の前記モル分率ｘの総和が１に等しいという第二条件、前記モル分率算出ステップで算出した前記各成分の前記モル分率ｙの総和が１に等しいという第三条件、並びに、前記モル密度算出ステップで算出した前記モル密度ρL ，ρV が下記式２を満足するという第四条件がそれぞれ満たされているかを判定する条件判定ステップと、前記条件判定ステップにおいて前記第一条件から前記第四条件までの４つの条件のうち少なくとも１つが満たされていなかった場合は、前記数値仮定ステップに戻って前記圧力Ｐ，前記物質量比ξ，及び前記分配率Ｋを別の数値に変更し、前記モル分率算出ステップ，前記モル密度算出ステップ，前記フガシティー算出ステップ，及び前記分配率算出ステップの計算を再度行い、前記４つの条件が全て満たされるまで前記数値仮定ステップ，前記モル分率算出ステップ，前記モル密度算出ステップ，前記フガシティー算出ステップ，及び前記分配率算出ステップの操作を繰り返す繰り返しステップと、前記条件判定ステップにおいて前記４つの条件が全て満たされた場合は、前記モル分率算出ステップで算出した前記各成��の前記モル分率ｘをそれぞれ出力する出力ステップと、を有する処理をコンピュータに実行させることを特徴とする。
【０００７】
【数１】
2011126965.tif 000003
【０００８】
【数２】
2011126965.tif 000004
【０００９】
また、本発明に係る気体組成推算装置は、少なくとも１種類の炭化水素化合物を含有する多成分系の液化ガスを、所定温度下で所定容積の密閉空間内に仕込んだ場合の、前記密閉空間内の気相の組成を推算する装置であって、前記温度Ｔ，前記密閉空間の容積Ｖ，前記液化ガスの仕込物質量Ｆ，及び仕込んだ前記液化ガスを構成する各成分のモル分率ｚを入力する入力手段と、前記入力手段により入力された前記温度Ｔ，前記密閉空間の容積Ｖ，前記仕込物質量Ｆ，及び前記各成分のモル分率ｚに基づいて前記密閉空間内の圧力Ｐを仮定するとともに、前記密閉空間内の気相の物質量Ｄと前記仕込物質量Ｆとの比Ｄ／Ｆ、及び、前記密閉空間内の気相の各成分のモル分率ｙと液相の各成分のモル分率ｘとの比ｙ／ｘを任意の数値に仮定する数値仮定手段と、前記数値仮定手段により仮定した前記比Ｄ／Ｆを物質量比ξ、前記比ｙ／ｘを分配率Ｋとし、式ｘ＝ｚ／｛１＋（Ｋ−１）ξ｝を用いて前記各成分の前記モル分率ｘをそれぞれ算出するとともに、算出した前記モル分率ｘを式Ｋ＝ｙ／ｘに代入して前記各成分の前記モル分率ｙをそれぞれ算出するモル分率算出手段と、前記モル分率算出手段により算出した前記モル分率ｘ，ｙを、Benedict-Webb-Rubin 状態方程式に代入して、前記密閉空間内の液相のモル密度ρL 及び気相のモル密度ρV をそれぞれ算出するモル密度算出手段と、前記密閉空間内の液相のフガシティーをｆL 、気相のフガシティーをｆV とし、下記式１（ただし、該式１中のｎは物質量、Ｚは圧縮係数であり、ｉは成分の種類を表す）を用いて前記各成分のｆL ／（ｘＰ）、ｆV ／（ｙＰ）をそれぞれ算出するフガシティー算出手段と、前記フガシティー算出手段により算出したｆL ／（ｘＰ）及びｆV ／（ｙＰ）から前記各成分の｛ｆL ／（ｘＰ）｝／｛ｆV ／（ｙＰ）｝を算出し、Ｋ’とする分配率算出手段と、前記分配率算出手段により算出した前記各成分のＫ’が前記分配率Ｋとそれぞれ等しいという第一条件、前記モル分率算出手段により算出した前記各成分の前記モル分率ｘの総和が１に等しいという第二条件、前記モル分率算出手段により算出した前記各成分の前記モル分率ｙの総和が１に等しいという第三条件、並びに、前記モル密度算出手段により算出した前記モル密度ρL ，ρV が下記式２を満足するという第四条件がそれぞれ満たされているかを判定する条件判定手段と、前記条件判定手段による判定で前記第一条件から前記第四条件までの４つの条件のうち少なくとも１つが満たされていなかった場合は、前記数値仮定手段により前記圧力Ｐ，前記物質量比ξ，及び前記分配率Ｋを別の数値に変更し、前記モル分率算出手段，前記モル密度算出手段，前記フガシティー算出手段，及び前記分配率算出手段による計算を再度行い、前記４つの条件が全て満たされるまで前記数値仮定手段，前記モル分率算出手段，前記モル密度算出手段，前記フガシティー算出手段，及び前記分配率算出手段による操作を繰り返す繰り返し手段と、前記条件判定手段による判定で前記４つの条件が全て満たされた場合は、前記モル分率算出手段により算出した前記各成分の前記モル分率ｘをそれぞれ出力する出力手段と、を備えることを特徴とする。
【００１０】
【数３】
2011126965.tif 000005
【００１１】
【数４】
2011126965.tif 000006
【発明の効果】
【００１２】
本発明の気体組成推算プログラム及び気体組成推算装置によれば、任意組成の多成分系液化ガスから気化した気化ガスの組成を定量的且つ容易に推算することができる。
【図面の簡単な説明】
【００１３】
【図１】プロパンのＰ−Ｖ−Ｔ関係を示すグラフである。
【図２】ジメチルエーテルのＰ−Ｖ−Ｔ関係を示すグラフである。
【図３】純物質の飽和蒸気圧を計算するフローチャートである。
【図４】プロパン及びブタンの飽和蒸気圧を示すグラフである。
【図５】ジメチルエーテルの飽和蒸気圧を示すグラフである。
【図６】２成分系気液平衡関係を計算するフローチャートである。
【図７】プロパン−ジメチルエーテル混合系の気液平衡関係を示すグラフである。
【図８】３成分系気液平衡関係の気液共存領域を示す三角線図である。
【図９】３成分系気液平衡関係のフラッシュ計算の概念を示す図である。
【図１０】３成分系気液平衡関係を計算するフローチャートである。
【図１１】プロパン−ジメチルエーテル−ブタン混合系の気液平衡関係を示すグラフである。
【図１２】体積収支を考慮した場合の３成分系気液平衡関係のフラッシュ計算の概念を示す図である。
【図１３】体積収支を考慮した場合の３成分系気液平衡関係を計算するフローチャートである。
【図１４】実施例１の計算結果を示すグラフである。
【図１５】実施例２の計算結果を示すグラフである。
【図１６】実施例３の計算結果を示すグラフである。
【発明を実施するための形態】
【００１４】
本発明に係る気体組成推算プログラム及び気体組成推算装置の実施の形態を、図面を参照しながら詳細に説明する。
ボンベ等の密閉空間内に多成分系の液化ガスを仕込むと、気相と液相とが形成される。本実施形態の気体組成推算プログラム及び気体組成推算装置は、少なくとも１種類の炭化水素化合物を含有する多成分系の液化ガスを、所定温度下で所定容積の密閉空間内に仕込んだ場合に形成される気相及び液相のそれぞれについて、組成及び物質量を定量的且つ容易に推算することができる上、密閉空間の内圧も推算することができる。
液化ガスを構成する成分の数は特に限定されるものではなく、２成分，３成分でも、あるいは４成分以上でも、問題なく推算を行うことができる。ただし、本発明に係る気体組成推算プログラム及び気体組成推算装置は、３成分以上の液化ガスに対して特に好適である。
【００１５】
また、複数の成分を含有する液化ガスの組成は特に限定されるものではなく（例えば共沸組成に限定されるものではなく）、任意組成の液化ガスを密閉空間内に仕込んだ場合の気相の組成等を推算することができる。
さらに、成分の種類は、複数の成分のうち少なくとも１種類が炭化水素化合物であり、残りの成分の種類は特に限定されない。つまり、液化ガスは、複数種の炭化水素化合物の混合物であってもよいし、１種又は複数種の炭化水素化合物と他種の化合物との混合物であってもよい。炭化水素化合物の種類は特に限定されるものではなく、エタン，プロパン，ブタン等を例示することができる。また、他種の化合物としては、例えば、エーテル類（ジメチルエーテル等），アルコール類（メチルアルコール等）があげられる。
【００１６】
密閉空間内に仕込む液化ガスの代表例としては、液化石油ガス（ＬＰＧ），液化天然ガス（ＬＮＧ）のような液化燃料ガスや、該液化燃料ガスにエーテル類，アルコール類を添加した液化混合燃料ガスがあげられる。ただし、本発明において使用可能な液化ガスは、燃料用の液化ガスに限定されるものではなく、他の用途の液化ガスを問題なく使用可能である。
本実施形態の気体組成推算装置は、後に詳述する入力手段，数値仮定手段，モル分率算出手段，モル密度算出手段，フガシティー算出手段，分配率算出手段，条件判定手段，繰り返し手段，及び出力手段を備えている。
【００１７】
入力手段により、温度Ｔ，密閉空間の容積Ｖ，液化ガスの仕込物質量Ｆ，及び仕込んだ液化ガスを構成する各成分のモル分率ｚを気体組成推算装置に入力すると、これらの入力値に基づいて、数値仮定手段が密閉空間内の圧力Ｐを仮定するとともに、密閉空間内の気相の物質量Ｄと仕込物質量Ｆとの比Ｄ／Ｆで定義される物質量比ξ、及び、密閉空間内の気相の各成分のモル分率ｙと液相の各成分のモル分率ｘとの比ｙ／ｘで定義される分配率Ｋを任意の数値に仮定する。
すると、モル分率算出手段が、入力手段により入力された各成分のモル分率ｚと、数値仮定手段により仮定された物質量比ξ，分配率Ｋと、式ｘ＝ｚ／｛１＋（Ｋ−１）ξ｝とを用いて気相の各成分のモル分率ｘをそれぞれ算出するとともに、算出したモル分率ｘを式Ｋ＝ｙ／ｘに代入して液相の各成分のモル分率ｙをそれぞれ算出する。
【００１８】
次に、モル密度算出手段が、モル分率算出手段により算出されたモル分率ｘ，ｙと、Benedict-Webb-Rubin 状態方程式（以降においては、ＢＷＲ式と記すこともある）とを用いて、密閉空間内の液相のモル密度ρL 及び気相のモル密度ρV をそれぞれ算出する。ＢＷＲ式の定数は物質によって異なる固有の数値であるが、計算の都度入力手段により入力するようにしてもよいし、気体組成推算プログラム又は気体組成推算装置にデータとして格納しておき、計算時にデータを取り出して使用するようにしてもよい。
モル分率算出手段及びモル密度算出手段の計算とは別に、フガシティー算出手段が前記式１を用いて各成分のｆL ／（ｘＰ）、ｆV ／（ｙＰ）をそれぞれ算出する。ｆL は密閉空間内の液相のフガシティーであり、ｆV は気相のフガシティーである。
【００１９】
次に、分配率算出手段が、フガシティー算出手段により算出されたｆL ／（ｘＰ）及びｆV ／（ｙＰ）から各成分の｛ｆL ／（ｘＰ）｝／｛ｆV ／（ｙＰ）｝（＝Ｋ’）を算出する。
なお、モル分率算出手段及びモル密度算出手段の計算と、フガシティー算出手段及び分配率算出手段の計算との時間的順序は特に問わない。すなわち、いずれか一方の計算を先に行ってもよいし、両計算を平行して行ってもよい。
【００２０】
次に、これらの計算が終了したら、条件判定手段が下記の４つの条件について判定を行う。
第一条件：分配率算出手段により算出した各成分のＫ’が分配率Ｋとそれぞれ等しい。
第二条件：モル分率算出手段により算出した各成分のモル分率ｘの総和が１に等しい。
第三条件：モル分率算出手段により算出した各成分のモル分率ｙの総和が１に等しい。
第四条件：モル密度算出手段により算出したモル密度ρL ，ρV が前記式２を満足する。
【００２１】
これらの判定を行って、前記４つの条件のうち１つでも満たされていない条件があった場合は、繰り返し手段が作動し、数値仮定手段が圧力Ｐ，物質量比ξ，及び分配率Ｋを別の数値に変更し、変更された数値を用いてモル分率算出手段，モル密度算出手段，フガシティー算出手段，及び分配率算出手段が上記計算を再度行う。そして、再計算が終了したら、条件判定手段が前記４つの条件について再判定を行う。前記４つの条件が全て満たされないと繰り返し手段が作動するようになっているので、数値仮定手段，モル分率算出手段，モル密度算出手段，フガシティー算出手段，及び分配率算出手段による操作が繰り返される。
【００２２】
条件判定手段による判定で前記４つの条件が全て満たされた場合は、繰り返し手段が作動しないようになっているので、モル分率算出手段により算出した気相の各成分のモル分率ｘが出力手段により出力される。
なお、圧力Ｐ，物質量比ξ，及び分配率Ｋを仮定して計算を行い、４つの条件が満たされなかった場合には、前記３つの数値の全てを別の数値に変更して再計算を行ってもよいが、数値の変更を段階的に行って再計算を行ってもよい。すなわち、前記３つの数値のうち２つは変更せず１つを順次変更して再計算を行い、４つの条件を満たす組み合わせが発見できなかった場合は別の１つの数値を順次変更して再計算を行うという方法（試行錯誤法）である。
【００２３】
このとき、圧力Ｐは、取り得る最高圧力と最低圧力が分かっているし（詳細は後述する）、物質量比ξは０と１の間の数値であるから、最初に圧力Ｐ及び物質量比ξを仮定する際には、上記上限値と下限値の中間の値を設定することが好ましい。そして、計算の結果、４つの条件が満たされなかったら、圧力Ｐや物質量比ξを別の数値に変更するが、その際の変更値は、前記上限値と前記中間の値との間の中間値、及び、前記下限値と前記中間の値との間の中間値のうち、前記計算結果から考えて４つの条件が一致する可能性が高いと思われる方とすることが好ましい。再計算の結果、４つの条件が満たされなかったら、再度同様にして数値の変更を行い、再計算を続けて４つの条件を満たす数値を探すとよい。このようにして数値を仮定すれば、全体の計算量を少なくすることができる。
【００２４】
また、条件判定手段による判定においては、誤差を考慮することが好ましい。例えば、Ｋ’がＫと等しいという第一条件を満足するか否かの判定に際しては、Ｋ’とＫとが完全に一致する可能性は低いので、一定の誤差範囲内で一致すれば第一条件を満足すると判定することが好ましい。この誤差範囲は、計算の都度所望の値に設定してもよいし、全ての計算において固定してもよい。
【００２５】
この気体組成推算装置は、入力手段，数値仮定手段，モル分率算出手段，モル密度算出手段，フガシティー算出手段，分配率算出手段，条件判定手段，繰り返し手段，及び出力手段などをソフトウェア上で実現するためのコンピュータシステムを備えている。このコンピュータシステムのハードウェア構成は、例えば、各種制御や演算処理を担う中央演算処理装置であるＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）と、主記憶装置（Ｍａｉｎ  Ｓｔｏｒａｇｅ）を構成するＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）と、読み出し専用の記憶装置であるＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）と、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）等の外部記憶装置（Ｓｅｃｏｎｄａｒｙ  Ｓｔｏｒａｇｅ）と、印刷装置，ディスプレイ装置等の出力装置と、操作パネル，マウス，キーボード，スキャナ等の入力装置と、を接続したものである。
【００２６】
そして、電源を投入すると、ＲＯＭに予め記憶されたコンピュータプログラム（気体組成推算プログラム）、又は、ＣＤ−ＲＯＭ，ＤＶＤ−ＲＯＭ，フレキシブルディスク（ＦＤ）等の記憶媒体を介して若しくはインターネット等の通信回線を介して記憶装置にインストールされたコンピュータプログラム（気体組成推算プログラム）がＲＡＭにロードされる。そのＲＡＭにロードされたプログラムに記述された命令に従って、ＣＰＵが所定の制御及び演算処理を行うことによって、前述のような各手段の各機能をソフトウェア上で実現できるようになっている。
【００２７】
なお、液化ガスの仕込物質量Ｆ及び仕込んだ液化ガスを構成する各成分のモル分率ｚを入力する代わりに、液化ガスの仕込質量及び仕込んだ液化ガスを構成する各成分の質量分率を入力し、これらの仕込質量及び質量分率から仕込物質量Ｆ及び各成分のモル分率ｚを算出してもよい。この算出に際しては、液化ガスを構成する各成分の分子量が必要となるが、この分子量の数値は入力手段により入力してもよいし、気体組成推算プログラム又は気体組成推算装置にデータとして格納しておき、計算時にデータを取り出して使用してもよい。
【００２８】
〔Benedict-Webb-Rubin 状態方程式について〕
まず、本発明において用いるBenedict-Webb-Rubin 状態方程式について説明する。ＢＷＲ式（下記の式３及び非特許文献１を参照）は、ＬＰＧの主成分である炭化水素化合物（例えばプロパン，ブタン）や、ＬＰＧに燃料成分として添加されるエーテル類（例えばジメチルエーテル），アルコール類（例えばメチルアルコール）の何れに対しても十分な精度を有することが確認されている、ビリアル展開型の状態方程式である。
ここで、Ｐ，Ｒ，Ｔ，ρはそれぞれ圧力，気体定数，温度，及びモル密度である。また、Ａ0 ，Ｂ0 ，Ｃ0 ，a ，b ，c ，α，γの８定数は物質固有の値であり、プロパン及びブタンについてはCooperらの文献値（非特許文献２を参照）により容易に入手することができる。
【００２９】
【数５】
2011126965.tif 000007
【００３０】
一方、ジメチルエーテルの値については報告されていなかったが、近年Tsuji らがJoffe の対応状態原理（非特許文献４を参照）とLee-Keslarの方法（非特許文献５を参照）とを併用し、さらに４０２．００ＫにおけるＰ−Ｖ−Ｔ関係の実測値を用いて、前記８定数を決定している（非特許文献３を参照）。プロパン，ジメチルエーテル，及びブタンに対するＢＷＲ式の８定数を表１に示す。
【００３１】
【表１】
2011126965.tif 000008
【００３２】
また、プロパン及びジメチルエーテルに対するＰ−Ｖ−Ｔ関係の文献値と計算値の比較を、図１，２のグラフに示す。実線が、ＢＷＲ式を用いた計算値である。また、白丸印（○印）がTsuji らの実測値であり、黒丸印（●印）がYounglove らのプロパン標準値（非特許文献６を参照）である。図１，２のグラフから、点線で示した臨界点の近傍を含めて、気相のＰ−Ｖ−Ｔ関係について推算結果は良好であることが分かる。
【００３３】
〔純物質の飽和蒸気圧について〕
一般に、純物質（１成分系）において気液両相が平衡状態にあるとき、次式４が成立する。すなわち、液相のフガシティーｆL と気相のフガシティーｆV とが一致する。なお、本発明においては、フガシティーの場合に限らず、添え字のＬ，Ｖはそれぞれ液相，気相を表す。
【００３４】
【数６】
2011126965.tif 000009
【００３５】
ここでフガシティーｆは、液相，気相ともに、ＢＷＲ式が圧力Ｐに対して与えられる場合は、次式５を適用する。ここで、Ｖは体積、ｎは物質量、Ｚ（＝ＰＶ／ｎＲＴ）は圧縮係数である。
【００３６】
【数７】
2011126965.tif 000010
【００３７】
すなわち、任意の温度に対して、式４，５を同時に満たすρL ，ρV が、それぞれ液相，気相の飽和密度に対応する。計算のフローチャートを図３に示す。はじめに温度Ｔを入力した後に圧力Ｐを仮定し、液相のモル密度ρL 及び気相のモル密度ρV を求める。これは、同一圧力において正の大きい値は液相の比容積１／ρL である。そして、式５を用いてｆL とｆV を算出し、ρL ，ρV に対して式４の成立を確認したら、平衡圧力が決定する。
【００３８】
式４が不成立の場合は、圧力Ｐを別の値に変更して、前記と同様にｆL とｆV を再度算出する。そして、式４が成立するまで圧力Ｐの変更を繰り返す。なお、温度Ｔを変化させた際の圧力軌跡は蒸気圧曲線となる。
ここで、プロパン，ブタン，及びジメチルエーテルのそれぞれの飽和蒸気圧の計算結果を、図４，５に示す。実線が、上記のようにして求めた計算値である。また、図４の白丸印（○印）がYounglove らの実測値である。さらに、図５の白丸印（○印）がTsuji らの実測値であり、黒丸印（●印）がHolldorff らの実測値である。図４，５から、飽和蒸気圧についても推算結果は良好であることが分かる。
【００３９】
〔２成分系気液平衡関係の計算について〕
多成分系の平衡条件は、純物質の場合と同様に各相の各成分のフガシティーが等しいことに帰着する（式６を参照）。
【００４０】
【数８】
2011126965.tif 000011
【００４１】
純物質と異なるのは、式６が各成分すべてに対して成立しなければならず、２成分系では２元、３成分系では３元の連立方程式を満たしていなければならない。また、多成分系に対するフガシティーも純物質とは多少異なる次式１で与えられる。
【００４２】
【数９】
2011126965.tif 000012
【００４３】
一方、多成分系にＢＷＲ式を適用するには、８つの定数を混合則と呼ばれる組成関数化する必要がある。ここでは、混合則として次式７〜１４に示すStotler-Benedict式（非特許文献７を参照）を用いた。
【００４４】
【数１０】
2011126965.tif 000013
【００４５】
【数１１】
2011126965.tif 000014
【００４６】
【数１２】
2011126965.tif 000015
【００４７】
【数１３】
2011126965.tif 000016
【００４８】
【数１４】
2011126965.tif 000017
【００４９】
【数１５】
2011126965.tif 000018
【００５０】
【数１６】
2011126965.tif 000019
【００５１】
【数１７】
2011126965.tif 000020
【００５２】
ここで、ｍijは異種分子間パラメータである。添え字ｉ，ｊは成分の種類を表す。プロパン−ジメチルエーテル混合系、ジメチルエーテル−ブタン混合系、及びブタン−プロパン混合系の２成分系の異種分子間パラメータを表２に示す。これは、２９８．１５Ｋにおけるデータを用いて最適化したものである。
【００５３】
【表２】
2011126965.tif 000021
【００５４】
気液平衡関係を計算するフローチャートを、図６に示す。はじめに温度Ｔ及び圧力Ｐを入力し、次に気相及び液相の組成（モル分率）ｘi ，ｙi を仮定する。ここで、各成分の分配率Ｋi （気相の各成分のモル分率ｙi と液相の各成分のモル分率ｘi との比ｙi ／ｘi ）も予め計算しておく。さらに、気相及び液相の個々に混合則を適用し、それぞれ異なる定数から気相及び液相のモル密度ρL ，ρV を得る。
【００５５】
フガシティーが等しくなることを利用し、式１から気液相の分配率Ｋ’i を算出し、はじめに仮定した分配率Ｋi と誤差範囲内において同じになる組成を見出す。すなわち、Ｋ’i ＝Ｋi が成立する組成ｘi ，ｙi を見出し、それらを出力する。Ｋ’i ＝Ｋi が不成立の場合は、仮定した組成ｘi ，ｙi を別の値に変更して、前記と同様にＫ’i を再度算出する。そして、Ｋ’i ＝Ｋi が成立するまで組成ｘi ，ｙi の変更を繰り返す。さらに、このとき純物質と同様に飽和密度も同時に求めていることになる。
計算結果の一例として、２９８．１５Ｋにおけるプロパン−ジメチルエーテル混合系の気液平衡関係の実測値と計算値を、図７のグラフに示す。図７のグラフより、実線で示した計算値は、Tsuji らの実測値（○印）を良く再現していることがわかる。
【００５６】
〔３成分系気液平衡関係の計算について〕
Gibbs の相律によると３成分３相系では自由度は３となり、３成分系気液平衡関係を求めるには、温度及び圧力の他に組成を与える必要がある。３成分系気液平衡関係の気液共存領域を、図８の三角線図に示した。
一般に３成分系では温度及び圧力を与えると、いわゆる露点線と沸点線が決まる。すなわち、気液平衡は、平衡な露点と沸点を結ぶタイラインが無数に存在し、その軌跡が露点線と沸点線になる。例えば温度及び圧力の他に気相における第１成分の組成ｙ1 、又は液相における第１成分の組成ｘ1 を与えると、沸点線又は露点線と交差して、その点を通るタイラインが得られる。しかし、ＢＷＲ式を用いた解法では、２成分系で示したように温度及び圧力を与えるのみで組成を与えない。そこで、仕込点Ｆを通るタイラインは必ず１本のみ定まることを利用して、物質収支を考えてフラッシュ計算を行うものである。
図９にフラッシュ計算の概念を示した。フラッシュ計算では、仕込物質量Ｆ，組成（モル分率）ｚ1 ，ｚ2 ，ｚ3 の３成分系液化ガスを温度Ｔのボンベにフラッシュさせて、圧力がＰになるようにする。ここで、成分ｉの物質収支を考えると、下記の式１５が成立する。
【００５７】
【数１８】
2011126965.tif 000022
【００５８】
一方、成分ｉの気相及び液相の分配率Ｋi は、次式１６で定義される。
【００５９】
【数１９】
2011126965.tif 000023
【００６０】
したがって、気相及び液相のモル分率は次式１７，１８となる。
【００６１】
【数２０】
2011126965.tif 000024
【００６２】
【数２１】
2011126965.tif 000025
【００６３】
ここで、ξ（＝Ｄ／Ｆ）は、全体に対する気相の物質量比となる。これを三角線図で考えると、次式１９に対応している。
【００６４】
【数２２】
2011126965.tif 000026
【００６５】
実際の計算のフローチャートを図１０に示した。はじめにＴ，Ｐ，ｚi を入力して、次にξを仮定し、さらにＫi を仮定する。なお、ξは０〜１の数値である。前述の２成分系の場合とほぼ同様に、式１，１７，１８，２０から気液相の分配率Ｋ’i を算出し、はじめに仮定した分配率Ｋi と誤差範囲内において同じになるか否か判定する。Ｋ’i ＝Ｋi が不成立の場合は、先ほど仮定したＫi を別の値に変更し、Ｋi ＝Ｋ’i を満たすＫi を試行錯誤法により求める（前述の２成分系の場合とほぼ同様に、数値の変更と再計算を繰り返す）。
【００６６】
【数２３】
2011126965.tif 000027
【００６７】
ただし、Ｋ’i ＝Ｋi が成立しても、Ｔ，Ｐ，ｚi ，ξより得た液相及び気相のモル分率ｘi ，ｙi は、各成分の総和が必ずしも１とはならないので、次式２１を用いて確認を行う。なお、式２１中の「ｉ」は、多成分系液化ガス中の第ｉ成分を意味する。また、式２１の総和記号は、気相又は液相の各成分の組成ｘ又は組成ｙ（モル分率）を、第１成分から第ｉ成分まで総和することを意味している。他の式の総和記号についても、ほぼ同様の意味である。
【００６８】
【数２４】
2011126965.tif 000028
【００６９】
図１０のフローチャートに示すように、式２１が成立するまで、先ほど仮定したξを別の値に変更して、前記と同様にＫ’i を再度算出する操作を繰り返す。そして、Ｋ’i ＝Ｋi が成立し且つ式２１が成立したら、組成ｘi ，ｙi を出力する。
計算結果の一例として、２９８．１５Ｋにおけるプロパン−ジメチルエーテル−ブタン混合系の気液平衡関係を、図１１に示した。
【００７０】
〔体積収支を考慮した場合の３成分系気液平衡関係の計算について〕
前述のフラッシュ計算では、ある組成の仕込物質量に対して温度及び圧力を与えることにより３成分気液平衡関係が求まることを示した。しかしながら、実際の運用を考えると、ボンベ等の密閉空間への液化ガスの充填前に圧力が規定されることは少なく、温度と内容積が既知の密閉空間に組成と質量を定めた液化ガスを仕込むことが多い。
ここでは、温度と密閉空間の内容積が既知の場合に、気相組成，液相組成，及び圧力を求める方法を考える。ＢＷＲ式による推算では、媒介変数として気相及び液相の比容積が算出されるので、これを考慮して体積収支を考えると次式２となる。
【００７１】
【数２５】
2011126965.tif 000029
【００７２】
体積収支を考慮したフラッシュ計算の概念を図１２に示すともに、計算フローチャートを図１３に示した。計算方法は基本的には前述の３成分系の場合と同様であるが、圧力が定まっていないので、その値を仮定する。仮定する範囲としては、仕込点が露点上にある場合の圧力（最低圧力Ｐmin ）と、沸点線上にある場合の圧力（最高圧力Ｐmax ）との範囲内である。この間で圧力Ｐの値を少しづつ変化させて、式２を満たすξを求めることになる。実際には、Ｐ，ξ，Ｋi を少しづつ変化させて試行錯誤的に求めることになるので、計算量はさらに多くなる。
【００７３】
実際の計算においては、はじめにＴ，Ｖ，ｚi を入力して、次にＰ，ξ，Ｋi を仮定する。なお、ξは０〜１の数値である。前述の３成分系の場合とほぼ同様に、式１，１７，１８，２０から気液の分配率Ｋ’i を算出し、はじめに仮定した分配率Ｋi と誤差範囲内において同じになるか否か判定する。Ｋ’i ＝Ｋi が不成立の場合は、先ほど仮定したＫi を別の値に変更し、Ｋi ＝Ｋ’i を満たすＫi を試行錯誤法により求める（前述の３成分系の場合とほぼ同様に、数値の変更と再計算を繰り返す）。
さらに、前述の３成分系の場合と同様に、Ｋ’i ＝Ｋi が成立してもＴ，Ｖ，Ｐ，ｚi ，ξより得た液相及び気相のモル分率ｘi ，ｙi は、各成分の総和が必ずしも１とはならないので、式２１を用いて確認を行う。
【００７４】
図１３のフローチャートに示すように、式２１が成立するまで、ξを別の値に変更して前記と同様にＫ’i を再度算出する操作を繰り返す。そして、Ｋ’i ＝Ｋi が成立し且つ式２１が成立したら、さらに式２の成立を確認する。式２が不成立の場合は、Ｐを別の値に変更し、式２が成立するように数値の変更と再計算を繰り返す。そして、Ｋ’i ＝Ｋi が成立し且つ式２１，式２が成立したら、組成ｘi ，ｙi を出力する。なお、圧力Ｐ等の他の数値を同時に出力してもよい。
【実施例】
【００７５】
次に、プロパン−ジメチルエーテル−ブタンの３成分系液化ガスを、温度２９８．１５Ｋ（２５℃）下で容積１１７．５Ｌのボンベに仕込んだ場合に形成される気相及び液相のそれぞれについて、組成及び物質量を推算した。
まず、前述した本実施形態の気体組成推算装置を起動し、コンピュータシステム上で本実施形態の気体組成推算プログラムを実行させる。すると、入力ステップにより数値入力待ち状態となるので、温度Ｔ（２９８．１５Ｋ）、ボンベの容積Ｖ（１１７．５Ｌ）、仕込んだ液化ガスの質量ｗ（５０ｋｇ）、及び液化ガスの各成分（プロパン，ジメチルエーテル，ブタン）の質量分率ｚwiを入力する。
【００７６】
プロパン，ジメチルエーテル，ブタンの分子量をそれぞれＭｗ1 ，Ｍｗ2 ，Ｍｗ3 とすると、仕込物質量Ｆと液化ガスの各成分のモル分率ｚi は次式２２，２３となるので、これらの式に従って仕込物質量Ｆと各成分のモル分率ｚi を算出する。なお、分子量Ｍｗ1 ，Ｍｗ2 ，Ｍｗ3 は、計算の都度入力してもよいが、気体組成推算装置のＲＯＭ，ＨＤＤ等に予め格納されたデータを読み出して計算に使用してもよい。また、仕込物質量Ｆと液化ガスの各成分のモル分率ｚi は、このように算出せずに入力するようにしてもよい。
【００７７】
【数２６】
2011126965.tif 000030
【００７８】
【数２７】
2011126965.tif 000031
【００７９】
次に、入力された各数値や算出された各成分のモル分率ｚi を数値仮定ステップに与えて、均一気相の圧力を求める。このとき、ξ＝１であり、得られた均一気相の圧力は最低圧力Ｐmin になる。そこで、数値仮定ステップにより、Ｐmin よりも僅かに大きい値を圧力Ｐとして仮定するとともに、さらにξを０から１の間の任意の数値に、Ｋi を任意の数値にそれぞれ仮定して、これらの仮定した数値をモル分率算出ステップに送る。
【００８０】
モル分率算出ステップでは、これらの仮定した数値から前述と同様にして気相及び液相のモル分率ｘi ，ｙi を算出する。さらに、入力された各数値や仮定した各数値や算出されたモル分率ｘi ，ｙi を、モル密度算出ステップやフガシティー算出ステップに送り、前述と同様にして気相及び液相のモル密度ρL ，ρV を算出するとともに、各成分のｆL ／（ｘＰ）、ｆV ／（ｙＰ）をそれぞれ算出する。さらに、フガシティー算出ステップで算出したｆL ／（ｘＰ）及びｆV ／（ｙＰ）を分配率算出ステップに送り、前述と同様にして分配率Ｋ’を算出する。
【００８１】
次に、得られた分配率Ｋ’i を条件判定ステップに送り、はじめに仮定した分配率Ｋi と誤差範囲内において同じになるか判定する。Ｋ’i ＝Ｋi が不成立の場合は、繰り返しステップによって処理が数値仮定ステップに戻る。そして、数値仮定ステップによりＫi を別の値に変更して、Ｋi ＝Ｋ’i が成立するまでＫi の変更を繰り返す。
Ｋ’i ＝Ｋi が成立した場合は、条件判定ステップにより、液相及び気相の各成分のモル分率ｘi ，ｙi の総和がそれぞれ１になるか否か判定する（前述の式２１を参照）。式２１が不成立の場合は、繰り返しステップによって処理が数値仮定ステップに戻る。そして、数値仮定ステップによりξを別の値に変更して、式２１が成立するまでξの変更を繰り返す。
【００８２】
さらに、Ｋ’i ＝Ｋi が成立し且つ式２１が成立したら、条件判定ステップにより前述の式２の成立を確認する。式２が不成立の場合は、繰り返しステップによって処理が数値仮定ステップに戻る。そして、数値仮定ステップにより圧力Ｐを別の値に変更し、式２が成立するまで圧力Ｐの変更を繰り返す。
そして、Ｋ’i ＝Ｋi が成立し且つ式２１，式２が成立したら計算を終了し、組成ｘi ，ｙi や圧力Ｐ等を出力ステップに送って、印刷装置，ディスプレイ装置等の出力装置に出力する。
【００８３】
次に、気相の一部をボンベから抜き出した状態を仮定する。気相の一部を抜き出すと、ボンベ内の液化ガスの質量及び液化ガスの各成分（プロパン，ジメチルエーテル，ブタン）の質量分率が変化するが、抜き出し後のボンベ内に残る液化ガスの新しい質量ｗnew 及び液化ガスの各成分の新しい質量分率ｚwinew は、次式２４，２５により計算できる。なお、式２４，２５中のΔＤは、抜き出した気化ガスの質量（消費量）である。
【００８４】
【数２８】
2011126965.tif 000032
【００８５】
【数２９】
2011126965.tif 000033
【００８６】
そして、式２４，２５の計算結果を当てはめることにより、気相の一部を抜き出した後の組成ｘi ，ｙi や圧力Ｐを、前述と同様にして推算することができる。この抜き出し及び推算を繰り返せば、液化ガスを仕込んだボンベから気相の気化ガスを徐々に抜き出していった場合のボンベ内の液化ガスの質量、圧力、組成ｘi ，ｙi 等を推算することができる。なお、計算量が膨大になるため、気化速度は熱移動に追随して影響を及ぼさない範囲で計算する。また、計算は、ボンベ内の圧力が７０ｋＰａゲージ圧となるまで、又は、液相が消滅までとする。
【００８７】
表３の実施例１のような条件において上記のような推算を行った結果を、図１４のグラフに実線で示す。グラフの横軸は、抜き出した気化ガスの量（ガス消費量）ΔＤであり、縦軸は、気相のジメチルエーテルの質量分率である。なお、液化ガスがプロパン９０質量％とジメチルエーテル１０質量％の２成分系である場合の推算結果を、図１４のグラフに破線で示してある。
【００８８】
【表３】
2011126965.tif 000034
【００８９】
図１４のグラフから、ジメチルエーテルの濃縮の過程が再現されていることが分かる。また、ブタンの混入により、気相におけるジメチルエーテルの濃縮がやや緩慢になる傾向にあることが分かる。ただし、現在の計算プログラムでは、ボンベ内の液化ガスの残量が２．８９ｋｇ付近になると、気液共存領域の追跡が困難となるため、最終的にボンベ内が気相のみになる場合の気化ガスの組成は得られていない。
【００９０】
また、表３の実施例２，３のような条件において上記のような推算を行った結果を、図１５，１６のグラフに実線で示す。なお、図１５，１６のグラフにおいても、液化ガスがプロパン９０質量％とジメチルエーテル１０質量％の２成分系である場合の推算結果を、図１４のグラフと同様に破線で示してある。また、図１５，１６のグラフには、アタム技研株式会社で行われた実測データをプロット（○印）してある。
【００９１】
実施例２では、３成分系の推算結果よりも２成分系の推算結果の方が実測値を良く再現していると言える。一方、実施例３では全く異なる傾向となり、３成分系の推算結果の方が実測値を良く再現している。これらの結果から、ブタンの質量分率が１質量％程度と低い場合は、２成分系として計算しても差し支えないが、ブタンの質量分率が５質量％程度と多くなると、３成分系の計算が必要と思われる。
【産業上の利用可能性】
【００９２】
本発明の気体組成推算プログラム及び気体組成推算装置は、液化ガスから発生した気化ガスを燃焼させるガス燃焼器（例えば、自動車等のエンジン，ボイラー，湯沸かし器，ガスコンロ，ガスストーブ）の設計、ボイラーの運転制御、ガスエンジンへのガス噴射タイミングの設定等に好適に利用可能である。
例えば、ガス燃焼器においては、ボンベからの消費に伴って組成が徐々に変化する気化ガスを最適な条件で燃焼させる必要があるので、気化ガスの組成の変化に応じて空気混合比等の燃焼条件の設定を自動的に最適化する機構や、気化ガスの組成が変化しても問題なく気化ガスを燃焼させることが可能な構造を開発する必要があるが、本発明の気体組成推算プログラム及び気体組成推算装置は、上記のような機構や構造の開発に好適に利用可能である。
【００９３】
また、家庭用，業務用ガスボンベに使用される燃料ガス、都市ガス、ガスライター用ガス、モータガス（タクシー用燃料）、燃料電池用ガス（燃料電池の水素源）の製造にも好適に利用可能である。
例えば、多成分系の液化ガスをボンベに充填して気化ガスを使用し、ボンベ内の液化ガスが少量になったら液化ガスを再充填して再使用する場合は、ボンベ内に残存している液化ガスと再充填される液化ガスの組成が異なるので、再充填を繰り返すうちにボンベ内の液化ガスの組成が元の液化ガスとは大きく異なってしまい、燃料ガスとして使用困難な組成になってしまうおそれがある。このような場合に、本発明の気体組成推算プログラム及び気体組成推算装置を利用すれば、再充填を何回繰り返すとボンベ内の液化ガスの組成が燃料ガスとして使用困難な組成になるのか、実験しなくても予測することができる。
【図１】
2011126965.tif 000039
【図２】
2011126965.tif 000040
【図３】
2011126965.tif 000041
【図４】
2011126965.tif 000042
【図５】
2011126965.tif 000043
【図６】
2011126965.tif 000044
【図７】
2011126965.tif 000045
【図８】
2011126965.tif 000046
【図９】
2011126965.tif 000047
【図１０】
2011126965.tif 000048
【図１１】
2011126965.tif 000049
【図１２】
2011126965.tif 000050
【図１３】
2011126965.tif 000051
【図１４】
2011126965.tif 000052
【図１５】
2011126965.tif 000053
【図１６】
2011126965.tif 000054