(57)【要約】【課題】未知の物質の飽和蒸気圧および活量係数を容易に算出できるようにする。【解決手段】シミュレーション装置は、実験データに基づいて混合溶液中の各物質の気相および液相のモル分率を決定し（Ｓ１０）、適宜決定した定数を利用しＮＲＴＬ式およびＡｎｔｏｉｎｅ式に従って各物質についての活量係数と飽和蒸気圧を暫定的に算出し（Ｓ２０〜Ｓ５０）、そして、算出結果である活量係数および飽和蒸気圧と上記の液相のモル分率を利用して各物質の気相のモル分率（ｙｉ（ｃａｌ））を算出する（Ｓ６０〜Ｓ８０）。そして、実験データに基づいて決定した気相のモル分率（ｙｉ（ｅｘｐ））とｙｉ（ｃａｌ）から各物質についての△ｙｉを算出し、混合溶液内のすべての物質の△ｙｉの総和が所定の条件を満足すれば、そのときに利用した各物質についての定数を、各物質の飽和蒸気圧および活量係数を算出するための式の定数として出力する（Ｓ１１０）。【選択図】図２
2007212356.tif 000002
【特許請求の範囲】
【請求項１】
物質の飽和蒸気圧および活量係数を算出するための物性値算出方法であって、
特定の物性値が既知の基準物質を含み密閉容器内で気液平衡状態にされた混合溶液について、前記基準物質の前記特定の物性値に基づいて、前記混合溶液に含まれる各物質の液相および気相のモル分率を決定するモル分率決定ステップと、
前記モル分率決定ステップで決定した前記混合溶液に含まれる各物質の液相および気相のモル分率を、前記混合溶液の温度と関連付けて記憶するステップと、
第１の定数群を含み温度と飽和蒸気圧の相関式である第１の相関式に基づいて前記混合溶液内の各物質の暫定的な飽和蒸気圧を算出する飽和蒸気圧算出ステップと、
第２の定数群を含み温度と液相のモル分率の相関式である第２の相関式に基づいて前記混合溶液内の各物質の暫定的な活量係数を算出する活量係数算出ステップと、
前記活量係数算出ステップで算出した暫定的な活量係数、前記飽和蒸気圧算出ステップで算出した暫定的な飽和蒸気圧、および、前記モル分率決定ステップで決定した前記混合溶液内の各物質の液相のモル分率に基づいて、前記密閉容器内の気液平衡状態での前記混合溶液内の各物質の分圧を算出する分圧算出ステップと、
前記分圧算出ステップで算出した分圧に基づいて、前記密閉容器内の気相の全圧を算出する全圧算出ステップと、
前記分圧算出ステップで算出した各物質の分圧と前記全圧算出ステップで算出した気相の全圧に基づいて、前記所定の温度での前記密閉容器内の前記混合溶液内の各物質の気相のモル分率を算出する気相モル分率算出ステップと、
前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで算出した気相のモル分率が所定の条件を満たしているか否かを判断する判断ステップと、
前記判断ステップで前記所定の条件を満たしていると判断した場合に、前記暫定的な活量係数の算出に利用した第１の定数群を前記所定の物質に関する前記所定の温度での前記第１の相関式に利用される定数群として出力し、かつ、前記暫定的な飽和蒸気圧の算出に利用した第２の定数群を前記所定の物質についての第２の相関式に利用される定数群として出力する、定数群出力ステップとを含む、物性値算出方法。
【請求項２】
前記活量係数算出ステップおよび前記飽和蒸気圧算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記第１の定数群および前記第２の定数群の少なくとも一方に代入する値を変更して、前記混合溶液内の各物質の前記暫定的な飽和蒸気圧および前記暫定的な活量係数を再度算出し、
前記分圧算出ステップ、前記全圧算出ステップ、および、前記気相モル分率算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、再度算出された前記暫定的な飽和蒸気圧および前記暫定的な活量係数に基づいて、前記混合溶液内の各物質の分圧、前記混合溶液内の気相の全圧、および、前記混合溶液内の各物質の気相のモル分率を算出し、
前記判断ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで再度算出された前記気相のモル分率が所定の条件を満たしているか否かを判断する、請求項１に記載の物性値算出方法。
【請求項３】
前記モル分率決定ステップは、
前記密閉容器内の気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフの入力を受付けるステップと、
前記密閉容器内で気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフにおいて各ピークの面積を算出するステップと、
前記液相および気相のガスクロマトグラフの前記基準物質に対応するピーク以外のピークに対応する成分のマススペクトルに基づいて、前記混合溶液内の各物質の分子量を決定するステップと、
前記ガスクロマトグラフにおける各ピークの面積と前記決定された前記混合溶液内の各物質の分子量とに基づいて、前記所定の物質の液相のモル分率および気相のモル分率を決定するステップとを含む、請求項１または請求項２に記載の物性値算出方法。
【請求項４】
前記第１の定数群を含む温度と飽和蒸気圧との相関式は、飽和蒸気圧をｐｓａｔとし、第１の定数群をＡ、Ｂ、および、Ｃとし、前記所定の温度の絶対温度をＴとした場合に、式（１）で表される、
【数１】
2007212356.tif 000026
請求項１〜請求項３のいずれかに記載の物性値算出方法。
【請求項５】
前記第２の定数群を含む温度と液相のモル分率との相関式は、前記所定の物質を含む溶液がＮ個の物質を含み、前記所定の物質ｉの活量係数をγｉとし、前記Ｎ個の物質の中の前記所定の物質以外の物質を物質ｊとし、Ｔを気体定数とし、ｘを各物質の前記所定の物質を含む溶液の中のモル分率とし、物質ｊと物質ｉに関する前記第２の定数群であるＧｉｊ、Ｇｊｉ、τｉｊ、τｊｉ、αｉｊ、αｊｉ、ｇｉｊ、ｇｊｊ、ｇｊｉ、ｇｉｉをそれぞれ次の式（２）〜（６）で表されるものとした場合に、式（７）で表される、
【数２】
2007212356.tif 000027
請求項１〜請求項４のいずれかに記載の物性値算出方法。
【請求項６】
物質の飽和蒸気圧および活量係数を算出するための物性値算出プログラムであって、
コンピュータに、
特定の物性値が既知の基準物質を含み密閉容器内で気液平衡状態にされた混合溶液について、前記基準物質の前記特定の物性値に基づいて、前記混合溶液に含まれる各物質の液相および気相のモル分率を決定するモル分率決定ステップと、
前記モル分率決定ステップで決定した前記混合溶液に含まれる各物質の液相および気相のモル分率を、前記混合溶液の温度と関連付けて記憶するステップと、
第１の定数群を含み温度と飽和蒸気圧の相関式である第１の相関式に基づいて前記混合溶液内の各物質の暫定的な飽和蒸気圧を算出する飽和蒸気圧算出ステップと、
第２の定数群を含み温度と液相のモル分率の相関式である第２の相関式に基づいて前記混合溶液内の各物質の暫定的な活量係数を算出する活量係数算出ステップと、
前記活量係数算出ステップで算出した暫定的な活量係数、前記飽和蒸気圧算出ステップで算出した暫定的な飽和蒸気圧、および、前記モル分率決定ステップで決定した前記混合溶液内の各物質の液相のモル分率に基づいて、前記密閉容器内の気液平衡状態での前記混合溶液内の各物質の分圧を算出する分圧算出ステップと、
前記分圧算出ステップで算出した分圧に基づいて、前記密閉容器内の気相の全圧を算出する全圧算出ステップと、
前記分圧算出ステップで算出した各物質の分圧と前記全圧算出ステップで算出した気相の全圧に基づいて、前記所定の温度での前記密閉容器内の前記混合溶液内の各物質の気相のモル分率を算出する気相モル分率算出ステップと、
前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで算出した気相のモル分率が所定の条件を満たしているか否かを判断する判断ステップと、
前記判断ステップで前記所定の条件を満たしていると判断した場合に、前記暫定的な活量係数の算出に利用した第１の定数群を前記所定の物質に関する前記所定の温度での前記第１の相関式に利用される定数群として出力し、かつ、前記暫定的な飽和蒸気圧の算出に利用した第２の定数群を前記所定の物質についての第２の相関式に利用される定数群として出力する、定数群出力ステップとを実行させる、物性値算出プログラム。
【請求項７】
前記活量係数算出ステップおよび前記飽和蒸気圧算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記第１の定数群および前記第２の定数群の少なくとも一方に代入する値を変更して、前記混合溶液内の各物質の前記暫定的な飽和蒸気圧および前記暫定的な活量係数を再度算出し、
前記分圧算出ステップ、前記全圧算出ステップ、および、前記気相モル分率算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、再度算出された前記暫定的な飽和蒸気圧および前記暫定的な活量係数に基づいて、前記混合溶液内の各物質の分圧、前記混合溶液内の気相の全圧、および、前記混合溶液内の各物質の気相のモル分率を算出し、
前記判断ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで再度算出された前記気相のモル分率が所定の条件を満たしているか否かを判断する、請求項６に記載の物性値算出プログラム。
【請求項８】
前記モル分率決定ステップは、
前記密閉容器内の気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフの入力を受付けるステップと、
前記密閉容器内で気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフにおいて各ピークの面積を算出するステップと、
前記液相および気相のガスクロマトグラフの前記基準物質に対応するピーク以外のピークに対応する成分のマススペクトルに基づいて、前記混合溶液内の各物質の分子量を決定するステップと、
前記ガスクロマトグラフにおける各ピークの面積と前記決定された前記混合溶液内の各物質の分子量とに基づいて、前記所定の物質の液相のモル分率および気相のモル分率を決定するステップとを含む、請求項６または請求項７に記載の物性値算出プログラム。
【請求項９】
前記第１の定数群を含む温度と飽和蒸気圧との相関式は、飽和蒸気圧をｐｓａｔとし、第１の定数群をＡ、Ｂ、および、Ｃとし、前記所定の温度の絶対温度をＴとした場合に、式（１）で表される、
【数３】
2007212356.tif 000028
請求項６〜請求項８のいずれかに記載の物性値算出プログラム。
【請求項１０】
前記第２の定数群を含む温度と液相のモル分率との相関式は、前記所定の物質を含む溶液がＮ個の物質を含み、前記所定の物質ｉの活量係数をγｉとし、前記Ｎ個の物質の中の前記所定の物質以外の物質を物質ｊとし、Ｔを気体定数とし、ｘを各物質の前記所定の物質を含む溶液の中のモル分率とし、物質ｊと物質ｉに関する前記第２の定数群であるＧｉｊ、Ｇｊｉ、τｉｊ、τｊｉ、αｉｊ、αｊｉ、ｇｉｊ、ｇｊｊ、ｇｊｉ、ｇｉｉをそれぞれ次の式（２）〜（６）で表されるものとした場合に、式（７）で表される、
【数４】
2007212356.tif 000029
請求項６〜請求項９のいずれかに記載の物性値算出プログラム。
【請求項１１】
請求項６〜請求項１０のいずれかに記載の物性値算出プログラムを記録した、コンピュータ読取可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、物性値算出方法に関し、特に、飽和蒸気圧および活量係数を算出する物性値算出方法、物性値算出プログラム、および、コンピュータ読取可能な記録媒体に関する。
【背景技術】
【０００２】
物質の飽和蒸気圧とは、当該物質の液相と蒸気とが一定温度において平衡状態にあるときに得られる当該物質の一定の蒸気圧力である。そして、このような物質の飽和蒸気圧は、たとえば、非特許文献１にも記載されるように、当該物質を蒸留操作によって分離および／または精製する際に、気相と液相の平衡状態での関係を考慮するために必要とされる、重要な物性値である。
【０００３】
従来から、混合液における気相と液相の平衡関係は、非特許文献２等に示されるように、以下に示す低圧気液平衡の式に基づいて求めることができる。
【０００４】
ｙｉｐＴ＝γｉｘｉｐｉｓａｔ（Ｔ）
なお、上記の式中で、ｉは混合液中の各物質を特定するための変数であり、ｙｉは気相中の物質ｉのモル分率であり、ｐＴは或る温度Ｔにおける全圧であり、γｉは物質ｉの活量係数であり、ｘｉは液相中の物質ｉのモル分率であり、そして、ｐｉｓａｔ（Ｔ）は或る温度Ｔでの物質ｉの飽和蒸気圧である。
【０００５】
つまり、ユーザは、混合液の中から物質ｉを単離する際には、当該物質ｉのγｉとｐｉｓａｔ（Ｔ）が得られれば、上記の式に従って、ある圧力、ある温度での物質ｉの液相組成（ｘｉ）および気相組成（ｙｉ）を得ることができる。
【０００６】
物質の飽和蒸気圧は、一般に、非特許文献３等に示されるように、Ａｎｔｏｉｎｅの式やＷａｇｎｅｒの式、Ｆｒｏｎｔ−Ｋａｌｋｗａｒｆの式に基づいて求めることができる。ここでは、一例として、Ａｎｔｏｉｎｅの式を式（１）として示す。
【０００７】
【数１】
2007212356.tif 000003
【０００８】
式（１）では、飽和蒸気圧はｐｓａｔで示され、また、Ａ，Ｂ，Ｃは物質ごとに決められる定数であり、Ｔは絶対温度（Ｋ）である。
【０００９】
また、活量係数は、一般に、非特許文献４等に示されるように、２成分系であれば、Ｍａｒｇｕｌｅｓ式、Ｖａｎ  Ｌａａｒ式、Ｗｉｌｓｏｎ式、ＮＲＴＬ（Ｎｏｎ  Ｒａｎｄｏｍ  Ｔｗｏ−Ｌｉｑｕｉｄ）式等に基づいて、また、多成分系であればＷｉｌｓｏｎ式、ＮＲＴＬ式、ＵＮＩＱＵＡＣ式等に基づいて、求めることができる。ここでは、一例として、多成分系のＮＲＴＬ式を式（７）として示す。式（７）では、活量係数はγｉで示されている。ｉとは多成分系の中で活量係数を求める対象となる物質に対応した記号であり、ｊは多成分系の中のそれ以外の物質に対応した記号である。また、式（２）〜式（６）には、式（７）で利用される定数であるＧｉｊ、Ｇｊｉ、τｉｊ、τｊｉ、αｉｊ、αｊｉ、ｇｉｊ、ｇｊｊ、ｇｊｉ、ｇｉｉの各値の定義が示されている。
【００１０】
【数２】
2007212356.tif 000004
【非特許文献１】平田光穂編著、「最新蒸留工学」、昭和４６年１月３０日発行、ｐ．７０−７４
【非特許文献２】社団法人化学工学会編、「化学工学便覧」、改訂六版、平成１１年２月２５日発行、ｐ．７４−７５
【非特許文献３】社団法人化学工学会編、「化学工学便覧」、改訂六版、平成１１年２月２５日発行、ｐ．５２
【非特許文献４】社団法人化学工学会編、「化学工学便覧」、改訂六版、平成１１年２月２５日発行、ｐ．６６−７４
【発明の開示】
【発明が解決しようとする課題】
【００１１】
上記した式（１）および式（７）等のような式を利用して活量係数を求めようとする場合、まず、各式で利用される定数を決定する必要がある。定数は、一般的には、まず、所定の物質（液相）を、密閉容器において所定の温度で平衡状態とし、当該容器中の液相および気相のそれぞれにおける所定の物質のモル分率を求め、その値に合うように活量係数の計算値を収束させることによって、決定されることが考えられる。
【００１２】
ただし、従来、多成分系の場合、混合溶液中の全ての物質の活量係数を求めるためには、２成分の組み合わせごとに液相や気相でのモル分率自体を求めた上で、推算や実験をすることによって定数および活量係数を決定しなければならず、実用的にこれらを求めることは困難であった。
【００１３】
また、混合溶液中の全ての物質の飽和蒸気圧を求めるには、まず、純物質を単離してその飽和蒸気圧を求めなければならず、多成分系の中に未知物質が存在する場合、通常、単離された当該未知物質は入手困難であることから、当該物質を含む混合溶液中のすべての物質の飽和蒸気圧を求めることは困難であった。
【００１４】
液相や気相でのモル分率は、対象となる物質の単離が可能であれば、容易に得られると考えられる。しかしながら、単離のために必要な飽和蒸気圧等の物性値が、単離を行なわないと得られないという問題があった。
【００１５】
本発明は、かかる実情に鑑み考えだされたものであり、その目的は、未知の物質を含む混合溶液において、当該未知の物質の飽和蒸気圧および活量係数を容易に算出できるようにすることである。
【課題を解決するための手段】
【００１６】
本発明に従った物性値算出方法は、物質の飽和蒸気圧および活量係数を算出するための物性値算出方法であって、特定の物性値が既知の基準物質を含み密閉容器内で気液平衡状態にされた混合溶液について、前記基準物質の前記特定の物性値に基づいて、前記混合溶液に含まれる各物質の液相および気相のモル分率を決定するモル分率決定ステップと、前記モル分率決定ステップで決定した前記混合溶液に含まれる各物質の液相および気相のモル分率を、前記混合溶液の温度と関連付けて記憶するステップと、第１の定数群を含み温度と飽和蒸気圧の相関式である第１の相関式に基づいて前記混合溶液内の各物質の暫定的な飽和蒸気圧を算出する飽和蒸気圧算出ステップと、第２の定数群を含み温度と液相のモル分率の相関式である第２の相関式に基づいて前記混合溶液内の各物質の暫定的な活量係数を算出する活量係数算出ステップと、前記活量係数算出ステップで算出した暫定的な活量係数、前記飽和蒸気圧算出ステップで算出した暫定的な飽和蒸気圧、および、前記モル分率決定ステップで決定した前記混合溶液内の各物質の液相のモル分率に基づいて、前記密閉容器内の気液平衡状態での前記混合溶液内の各物質の分圧を算出する分圧算出ステップと、前記分圧算出ステップで算出した分圧に基づいて、前記密閉容器内の気相の全圧を算出する全圧算出ステップと、前記分圧算出ステップで算出した各物質の分圧と前記全圧算出ステップで算出した気相の全圧に基づいて、前記所定の温度での前記密閉容器内の前記混合溶液内の各物質の気相のモル分率を算出する気相モル分率算出ステップと、前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで算出した気相のモル分率が所定の条件を満たしているか否かを判断する判断ステップと、前記判断ステップで前記所定の条件を満たしていると判断した場合に、前記暫定的な活量係数の算出に利用した第１の定数群を前記所定の物質に関する前記所定の温度での前記第１の相関式に利用される定数群として出力し、かつ、前記暫定的な飽和蒸気圧の算出に利用した第２の定数群を前記所定の物質についての第２の相関式に利用される定数群として出力する、定数群出力ステップとを含むことを特徴とする。
【００１７】
また、本発明に従った物性値算出方法では、前記活量係数算出ステップおよび前記飽和蒸気圧算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記第１の定数群および前記第２の定数群の少なくとも一方に代入する値を変更して、前記混合溶液内の各物質の前記暫定的な飽和蒸気圧および前記暫定的な活量係数を再度算出し、前記分圧算出ステップ、前記全圧算出ステップ、および、前記気相モル分率算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、再度算出された前記暫定的な飽和蒸気圧および前記暫定的な活量係数に基づいて、前記混合溶液内の各物質の分圧、前記混合溶液内の気相の全圧、および、前記混合溶液内の各物質の気相のモル分率を算出し、前記判断ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで再度算出された前記気相のモル分率が所定の条件を満たしているか否かを判断することが好ましい。
【００１８】
また、本発明に従った物性値算出方法では、前記モル分率決定ステップは、前記密閉容器内の気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフの入力を受付けるステップと、前記密閉容器内で気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフにおいて各ピークの面積を算出するステップと、前記液相および気相のガスクロマトグラフの前記基準物質に対応するピーク以外のピークに対応する成分のマススペクトルに基づいて、前記混合溶液内の各物質の分子量を決定するステップと、前記ガスクロマトグラフにおける各ピークの面積と前記決定された前記混合溶液内の各物質の分子量とに基づいて、前記所定の物質の液相のモル分率および気相のモル分率を決定するステップとを含むことが好ましい。
【００１９】
また、本発明に従った物性値算出方法では、前記第１の定数群を含む温度と飽和蒸気圧との相関式は、飽和蒸気圧をｐｓａｔとし、第１の定数群をＡ、Ｂ、および、Ｃとし、前記所定の温度の絶対温度をＴとした場合に、式（１）で表されることが好ましい。
【００２０】
【数３】
2007212356.tif 000005
【００２１】
また、本発明に従った物性値算出方法では、前記第２の定数群を含む温度と液相のモル分率との相関式は、前記所定の物質を含む溶液がＮ個の物質を含み、前記所定の物質ｉの活量係数をγｉとし、前記Ｎ個の物質の中の前記所定の物質以外の物質を物質ｊとし、Ｔを気体定数とし、ｘを各物質の前記所定の物質を含む溶液の中のモル分率とし、物質ｊと物質ｉに関する前記第２の定数群であるＧｉｊ、Ｇｊｉ、τｉｊ、τｊｉ、αｉｊ、αｊｉ、ｇｉｊ、ｇｊｊ、ｇｊｉ、ｇｉｉをそれぞれ次の式（２）〜（６）で表されるものとした場合に、式（７）で表されることが好ましい。
【００２２】
【数４】
2007212356.tif 000006
【００２３】
本発明に従った物性値算出プログラムは、物質の飽和蒸気圧および活量係数を算出するための物性値算出プログラムであって、コンピュータに、特定の物性値が既知の基準物質を含み密閉容器内で気液平衡状態にされた混合溶液について、前記基準物質の前記特定の物性値に基づいて、前記混合溶液に含まれる各物質の液相および気相のモル分率を決定するモル分率決定ステップと、前記モル分率決定ステップで決定した前記混合溶液に含まれる各物質の液相および気相のモル分率を、前記混合溶液の温度と関連付けて記憶するステップと、第１の定数群を含み温度と飽和蒸気圧の相関式である第１の相関式に基づいて前記混合溶液内の各物質の暫定的な飽和蒸気圧を算出する飽和蒸気圧算出ステップと、第２の定数群を含み温度と液相のモル分率の相関式である第２の相関式に基づいて前記混合溶液内の各物質の暫定的な活量係数を算出する活量係数算出ステップと、前記活量係数算出ステップで算出した暫定的な活量係数、前記飽和蒸気圧算出ステップで算出した暫定的な飽和蒸気圧、および、前記モル分率決定ステップで決定した前記混合溶液内の各物質の液相のモル分率に基づいて、前記密閉容器内の気液平衡状態での前記混合溶液内の各物質の分圧を算出する分圧算出ステップと、前記分圧算出ステップで算出した分圧に基づいて、前記密閉容器内の気相の全圧を算出する全圧算出ステップと、前記分圧算出ステップで算出した各物質の分圧と前記全圧算出ステップで算出した気相の全圧に基づいて、前記所定の温度での前記密閉容器内の前記混合溶液内の各物質の気相のモル分率を算出する気相モル分率算出ステップと、前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで算出した気相のモル分率が所定の条件を満たしているか否かを判断する判断ステップと、前記判断ステップで前記所定の条件を満たしていると判断した場合に、前記暫定的な活量係数の算出に利用した第１の定数群を前記所定の物質に関する前記所定の温度での前記第１の相関式に利用される定数群として出力し、かつ、前記暫定的な飽和蒸気圧の算出に利用した第２の定数群を前記所定の物質についての第２の相関式に利用される定数群として出力する、定数群出力ステップとを実行させることを特徴とする。
【００２４】
また、本発明に従った物性値算出プログラムでは、前記活量係数算出ステップおよび前記飽和蒸気圧算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記第１の定数群および前記第２の定数群の少なくとも一方に代入する値を変更して、前記混合溶液内の各物質の前記暫定的な飽和蒸気圧および前記暫定的な活量係数を再度算出し、前記分圧算出ステップ、前記全圧算出ステップ、および、前記気相モル分率算出ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、再度算出された前記暫定的な飽和蒸気圧および前記暫定的な活量係数に基づいて、前記混合溶液内の各物質の分圧、前記混合溶液内の気相の全圧、および、前記混合溶液内の各物質の気相のモル分率を算出し、前記判断ステップは、前記判断ステップで前記所定の条件を満たしていないと判断した場合に、前記モル分率決定ステップで決定した気相のモル分率と前記気相モル分率算出ステップで再度算出された前記気相のモル分率が所定の条件を満たしているか否かを判断することが好ましい。
【００２５】
また、本発明に従った物性値算出プログラムでは、前記モル分率決定ステップは、前記密閉容器内の気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフの入力を受付けるステップと、前記密閉容器内で気液平衡状態にある前記混合溶液の液相および気相のガスクロマトグラフにおいて各ピークの面積を算出するステップと、前記液相および気相のガスクロマトグラフの前記基準物質に対応するピーク以外のピークに対応する成分のマススペクトルに基づいて、前記混合溶液内の各物質の分子量を決定するステップと、前記ガスクロマトグラフにおける各ピークの面積と前記決定された前記混合溶液内の各物質の分子量とに基づいて、前記所定の物質の液相のモル分率および気相のモル分率を決定するステップとを含むことが好ましい。
【００２６】
また、本発明に従った物性値算出プログラムでは、前記第１の定数群を含む温度と飽和蒸気圧との相関式は、飽和蒸気圧をｐｓａｔとし、第１の定数群をＡ、Ｂ、および、Ｃとし、前記所定の温度の絶対温度をＴとした場合に、式（１）で表されることが好ましい。
【００２７】
【数５】
2007212356.tif 000007
【００２８】
また、本発明に従った物性値算出プログラムでは、前記第２の定数群を含む温度と液相のモル分率との相関式は、前記所定の物質を含む溶液がＮ個の物質を含み、前記所定の物質ｉの活量係数をγｉとし、前記Ｎ個の物質の中の前記所定の物質以外の物質を物質ｊとし、Ｔを気体定数とし、ｘを各物質の前記所定の物質を含む溶液の中のモル分率とし、物質ｊと物質ｉに関する前記第２の定数群であるＧｉｊ、Ｇｊｉ、τｉｊ、τｊｉ、αｉｊ、αｊｉ、ｇｉｊ、ｇｊｊ、ｇｊｉ、ｇｉｉをそれぞれ次の式（２）〜（６）で表されるものとした場合に、式（７）で表されることが好ましい。
【００２９】
【数６】
2007212356.tif 000008
【００３０】
本発明に従ったコンピュータ読取可能な記録媒体は、上記した物性値算出プログラムを記録していることを特徴とする。
【発明の効果】
【００３１】
本発明によれば、飽和蒸気圧が未知の物質を含む混合溶液についての物性の測定結果を利用することによって、種々の値の算出や種々の判断を行なうことによって、当該未知の物質についての飽和蒸気圧および活量係数を算出するための式の定数を得ることができる。これにより、未知の物質の飽和蒸気圧および活量係数を決定することができる。つまり、本発明によれば、未知の物質を含む混合溶液において、当該混合溶液から未知の物質を単離することなく、当該未知の物質の飽和蒸気圧および活量係数を算出できる。
【発明を実施するための最良の形態】
【００３２】
図１は、本実施の形態において物性値の算出を行なうための装置（シミュレーション装置）のブロック図である。シミュレーション装置は、当該装置の全体的な動作を制御するＣＰＵ（central processing unit）１１、ＣＰＵ１１の動作プログラムを含む種々の情報を記憶するハードディスク（ＨＤ）１３、外部から挿入されるＣＤ−ＲＯＭ（compact disc read only memory）やＤＶＤ（Digital Versatile Disk）等のディスク１００に記憶された情報を読取るディスクドライブ１４、ＣＰＵ１１のワークエリアとなるＲＡＭ（random access memory）１５、ユーザからの情報の入力を受付ける入力部１６、および、液晶ディスプレイ等で構成される表示部１８を備えている。ＣＰＵ１１は、入出力インターフェイス（Ｉ／Ｏ）１２を介して、ＨＤ１３、ＲＯＭ１４、ＲＡＭ１５、入力部１６、データベース１７、および、表示部１８と接続されている。
【００３３】
ＣＰＵ１１は、ＨＤ１３に記憶されているプログラムを実行することもできるし、ディスク１００等、シミュレーション装置から着脱可能な記録媒体に記録されたプログラムを実行することもできる。
【００３４】
図１に示されたシミュレーション装置は、混合溶液中の各物質について飽和蒸気圧と活量係数とを決定する処理（物性値算出処理）を実行する。ここで、このような処理について、当該処理のフローチャートである図２を参照しつつ説明する。なお、シミュレーション装置では、上記混合溶液中の各物質について物性値算出処理を実行することにより、各物質の飽和蒸気圧が算出される。そして、この混合溶液における各物質の活量係数を求めることができる。
【００３５】
物性値算出処理が実行される際には、分子量および飽和蒸気圧を求める式についての定数（本実施の形態では、Ａｎｔｏｉｎｅ式の定数Ａ〜Ｃ）が既知である物質（以下、基準物質と呼ぶ）を含む混合溶液が準備される。そして、当該混合溶液中の基準物質を含む各物質の液相および気相のモル分率を算出するために、当該混合溶液が気液平衡状態にされて、当該混合溶液の気相および液相のガスクロマトグラフが測定される。
【００３６】
また、物性値算出処理では、上記の混合溶液の各物質について、暫定的な飽和蒸気圧が算出される。この飽和蒸気圧の算出には、たとえば式（１）として示すＡｎｔｏｉｎｅの式を利用できる。
【００３７】
【数７】
2007212356.tif 000009
【００３８】
式（１）では、飽和蒸気圧はｐｓａｔで示されている。また、Ａ，Ｂ，Ｃは物質ごとに決められる定数であり、Ｔは絶対温度（Ｋ）である。
【００３９】
物性値算出処理では、上記の混合溶液中の各物質について、後述する所定の条件が成立するまで、パラメータＡ〜Ｃの値を変更して暫定的な飽和蒸気圧の算出が繰り返される。ここで、定数Ａ〜Ｃの初期値は特に限定されないが、Ａｎｔｏｉｎｅの式の場合、たとえば、Ａは０〜５０００程度、Ｂは−５０００〜０程度、Ｃは０〜５０００程度である。そして、所定の条件が成立すると、当該条件が成立したときに利用されたＡ〜Ｃが、対応する物質についてのＡ〜Ｃとされる。
【００４０】
なお、本発明では、飽和蒸気圧の算出に利用できる式はＡｎｔｏｉｎｅの式に限定されない。Ｗａｇｎｅｒの式や、Ｆｒｏｎｔ−Ｋａｌｋｗａｒｆの式等の、他の式に基づいて算出されても良い。本発明によれば、後述する所定の条件が成立するまで、混合溶液中の各物質について、飽和蒸気圧を算出するための式（Ａｎｔｏｉｎｅの式に限定されない）における定数に代入する値を変えることによって暫定的な飽和蒸気圧の算出が繰り返される。そして、所定の条件が成立すると、シミュレーション装置は、そのときに利用されていた定数の値を、その物質の飽和蒸気圧を算出するための式に利用される定数の値として出力する。
【００４１】
また、物性値算出処理では、上記の混合溶液中の各物質について、後述する所定の条件が成立するまで、算出の条件を変えて、暫定的な活量係数の算出が繰り返される。活量係数の算出には、たとえば、式（７）として示されるＮＲＴＬ式が利用される。
【００４２】
【数８】
2007212356.tif 000010
【００４３】
式（７）では、活量係数はγｉで示されている。ｉとは多成分系の中で活量係数を求める対象となる物質に対応した記号であり、ｊは多成分系の中のそれ以外の物質に対応した記号である。また、式（２）〜式（６）には、式（７）で利用される定数であるＧｉｊ、Ｇｊｉ、τｉｊ、τｊｉ、αｉｊ、αｊｉ、ｇｉｊ、ｇｊｊ、ｇｊｉ、ｇｉｉの各値の定義が示されている。
【００４４】
ここで、ｂｉｊ＝ｇｉｊ−ｇｊｊ、かつ、ｂｊｉ＝ｇｊｉ−ｇｉｉであり、さらに式（６）から、ＮＲＴＬ式では、定数は、ｂｉｊ、ｂｊｉ、αｉｊの３種類と考えることができる。上記したように暫定的な活量係数の算出が繰り返される際、算出の条件として、この３種類の定数の値が変更される。そして、所定の条件が成立すると、シミュレーション装置は、そのときに定数に代入されていた値を、その物質についての活量係数を算出するための定数として出力する。
【００４５】
また、式（７）では、ｘｉは混合溶液についての物質ｉの液相のモル分率を示す。また、本明細書では、ｙｉは混合溶液についての物質ｉの気相のモル分率を示す。
【００４６】
本発明では、活量係数の算出に利用する式はＮＲＴＬ式に限定されない。Ｍａｒｇｕｌｅｓ式、Ｖａｎ  Ｌａａｒ式、Ｗｉｌｓｏｎ式、Ｗｉｌｓｏｎ式、ＵＮＩＱＵＡＣ式等の他の式に基づいて算出されても良い。つまり、本発明によれば、所定の条件が成立するまで、活量係数を算出するための式（ＮＲＴＬ式に限定されない）が利用されて暫定的な活量係数の算出が繰り返される。なお、暫定的な活量係数の算出が繰り返される際には、当該式の定数として代入する値が変更される。そして、所定の条件が成立すると、シミュレーション装置は、そのときに定数に代入されていた値を、その物質についての活量係数を算出するための定数として出力する。
【００４７】
図２を参照して、物性値算出処理では、ＣＰＵ１１は、まずステップＳ１０で、上記した混合溶液中の各物質の液相のモル分率ｘｉおよび気相のモル分率ｙｉを算出する。ここでのモル分率は、たとえば、密閉容器内で気液平衡状態にある混合溶液の液相および気相のガスクロマトグラフから算出される。
【００４８】
具体的には、混合溶液の液相のガスクロマトグラフから、対象となる物質の面積百分率組成ｘＡＩを算出して、その値をｘｉに決定し、また、混合溶液の気相のガスクロマトグラフから、対象となる物質の面積百分率組成ｙＡＩを算出して、その値をｙｉに決定する。なお、ここで算出される気相のモル分率ｙｉは、混合溶液について実験的に求められる値である。一方、本実施の形態では、後述するように、Ａｎｔｏｉｎｅ式およびＮＲＴＬ式に基づいて暫定的に算出されるｙｉの計算値が得られる。そこで、本明細書では、２つのｙｉを区別するために、ステップＳ１０において実験的に求められるｙｉを適宜ｙｉ（ｅｘｐ）と記述し、暫定的に算出されるｙｉの計算値をｙｉ（ｃａｌ）と記述する。
【００４９】
液相のモル分率ｘｉおよび気相のモル分率ｙｉの算出について、さらに具体的に説明する。ＣＰＵ１１は、面積百分率組成ｘＡＩおよび面積百分率組成ｙＡＩをそれぞれ以下の式（８）および式（９）に基づいて算出する。
【００５０】
【数９】
2007212356.tif 000011
【００５１】
式（８）において、Ａｘｉ，Ａｘｊは、それぞれｉ，ｊで特定される物質の混合溶液の液相のガスクロマトグラフの面積百分率であり、Ｍｗｉ，Ｍｗｊは、それぞれｉ，ｊで特定される物質の分子量である。なお、物質の分子量は、未知の場合には、たとえば、ガスクロマトグラフ質量分析装置（ＧＣ−ＭＣ）において、ガスクロマトグラフの各ピークに対応するマススペクトルから、求めることができる。また、式（９）において、Ａｙｉ，Ａｙｊは、それぞれｉ，ｊで特定される物質の混合溶液の気相のガスクロマトグラフの面積百分率であり、Ｍｗｉ，Ｍｗｊは、それぞれｉ，ｊで特定される物質の分子量である。
【００５２】
また、本実施の形態では、測定対象となる混合溶液を、気液平衡状態として、気相と液相を導入することから、このような導入が可能なヘッドスペース試料導入装置を具備したガスクロマトグラフ装置を利用されることが好ましい。なお、本実施の形態のシミュレーション装置では、混合溶液の液相および気相についてのガスクロマトグラフおよびマススペクトル（の測定結果）が入力されればよく、必ずしも、ガスクロマトグラフ質量分析装置等の測定用の装置に直接的に接続されている必要はない。また、本実施の形態では、シミュレーション装置は、少なくとも混合溶液の液相および気相の各物質のモル分率を決定できればよく、決定の方法は、ガスクロマトグラフによるものに限定されない。
【００５３】
次に、ＣＰＵ１１は、ステップＳ２０で、ＮＲＴＬ式の定数（上記したｂｉｊ、ｂｊｉ、αｉｊ）を決定する。なお、この場合、定数は、初回には、予め定められた値またはユーザによって入力された値に決定され、それ以降（後述するステップＳ１００におけるＮＯ判断時に、処理が戻された際）には、たとえば収束関数等が利用されることによって、それまで決定された値から適宜変更されるように、決定される。
【００５４】
そして、ＣＰＵ１１は、ステップＳ３０で、直前に実行したステップＳ２０において決定した定数を利用して、物質ｉの活量係数γｉを算出する。
【００５５】
次に、ＣＰＵ１１は、ステップＳ４０で、Ａｎｔｏｉｎｅ式の定数（上記したＡ，Ｂ，Ｃ）を決定する。なお、この場合、定数は、初回には、予め定められた値またはユーザによって入力された値に決定され、それ以降（後述するステップＳ１００におけるＮＯ判断時に、処理が戻された際）には、たとえば収束関数等が利用されることによって、それまで決定された値から適宜変更されるように、決定される。
【００５６】
そして、ＣＰＵ１１は、ステップＳ５０で、直前に実行したステップＳ４０において決定した定数を利用して、物質ｉの飽和蒸気圧ｐｉｓａｔを算出する。
【００５７】
次に、ＣＰＵ１１は、ステップＳ６０で、ステップＳ３０で算出した活量係数γｉ、ステップＳ５０で算出した飽和蒸気圧ｐｉｓａｔ、および、ステップＳ１０で決定した液相のモル分率ｘｉを利用し、上記した低圧気液平衡の式に基づいて、物質ｉの分圧ｐｉを算出する。
【００５８】
次に、ＣＰＵ１１は、ステップＳ７０で、ステップＳ６０と同様に混合溶液中のすべての物質の気相での分圧を算出し、それらの総和（全圧ｐＴ）を算出する。
【００５９】
次に、ＣＰＵ１１は、ステップＳ８０で、ステップＳ６０で算出したｐｉとステップＳ７０で算出したｐＴを利用して、物質ｉの気相の分圧の計算値（ｙｉ（ｃａｌ））を算出する。
【００６０】
次に、偏差を算出する。具体的には、ＣＰＵ１１は、ステップＳ９０で、ステップＳ１０で決定したｙｉ（ｙｉ（ｅｘｐ））とステップＳ８０で算出したｙｉ（ｃａｌ）の差（絶対値）のｙｉ（ｃａｌ）に対する割合（△ｙｉ）を算出する。なお、本発明において、偏差の算出方法はこれに限定されない。たとえば、最小二乗法やニュートン法等の他の方法が利用されても良い。
【００６１】
そして、ＣＰＵ１１は、ステップＳ１００で、ステップＳ９０で算出した△ｙｉが予め定められた設定値以下であるか否かを判断し、そうであればステップＳ１１０でそのときに利用していた、各物質についてのＡｎｔｏｉｎｅ式の定数（Ａ，Ｂ，Ｃ）、および、ＮＲＴＬ式についての各２物質間の定数（ｂｉｊ、ｂｊｉ、αｉｊ）を出力して、処理を終了する。一方、△ｙｉが予め定められた設定値を越えている場合には、ステップＳ２０へ処理を戻す。これにより、ＮＲＴＬ式やＡｎｔｏｉｎｅ式に利用される定数が変更されて、再度、飽和蒸気圧ｐｉｓａｔと活量係数γｉの算出が行なわれる。
【００６２】
以上説明した本実施の形態では、まず、混合溶液中の各物質の気相および液相のモル分率が実験データ（ガスクロマトグラフ）に基づいて算出される（Ｓ１０）。
【００６３】
次に、各物質について、モル分率と温度と活量係数の相関式（ＮＲＴＬ式）および温度と飽和蒸気圧の相関式（Ａｎｔｏｉｎｅ式）に利用される定数が適宜決定され、決定された定数と温度やモル分率（実験データに基づいて算出されたモル分率）がこれらの式に代入されることにより、暫定的に活量係数と飽和蒸気圧が算出される（Ｓ２０〜Ｓ５０）。なお、暫定的な活量係数と飽和蒸気圧の算出に用いられる温度は、混合溶液の温度である。
【００６４】
そして、算出結果である暫定的な活量係数および飽和蒸気圧と、混合溶液についての実験データから求められる各物質の液相のモル分率とが利用されて、混合溶液の各物質についての気相のモル分率の計算値（上記したｙｉ（ｃａｌ））が算出される（Ｓ６０〜Ｓ８０）。
【００６５】
なお、ｙｉ（ｃａｌ）が算出される際には、まず、混合溶液中の各物質についての分圧が算出される（Ｓ６０）。分圧ｐｉは、以下の式（１０）で示す低圧気液平衡の式から導かれる式（１１）に基づき、上記の暫定的な活量係数および飽和蒸気圧と実験データから求められた液相のモル分率が利用されて、算出される。
【００６６】
ｙｉｐＴ＝γｉｘｉｐｉｓａｔ（Ｔ）  …（１０）
ｐｉ＝γｉｘｉｐｉｓａｔ（Ｔ）  …（１１）
そして、混合溶液内のすべての物質について算出された分圧ｐｉの和が算出されることによって、密閉容器内の全圧ｐＴが算出され（Ｓ７０）、各物質の全圧ｐＴに対する分圧ｐｉの比として、ｙｉ（ｃａｌ）が算出される（Ｓ８０）。
【００６７】
そして、上記した割合△ｙｉが算出されて（Ｓ１００）、割合△ｙｉに対する判断がなされる。なお、割合△ｙｉに対する判断は、混合溶液内のすべての物質のそれぞれに対して算出された割合△ｙｉの総和に対して行なわれる。
【００６８】
そして、所定の条件が成立した割合△ｙｉの算出に利用した各物質についての定数（Ａ，Ｂ，Ｃ，ｂｉｊ，ｂｊｉ，αｉｊ）が、各物質の飽和蒸気圧および活量係数を算出するための式の定数として、出力される（Ｓ１１０）。
【００６９】
以上説明した本実施の形態では、所定の条件が成立するまで上記した割合△ｙｉが算出されるが、本実施の形態のシミュレーション装置における割合△ｙｉの算出の終了条件はこれに限定されない。シミュレーション装置は、ユーザ等によって指定される回数だけ上記した割合△ｙｉを算出し、指定回数算出された割合△ｙｉをそのとき利用した各物質についての定数（Ａ，Ｂ，Ｃ，ｂｉｊ，ｂｊｉ，αｉｊ）の組とともにすべて記憶し、その中で所定の条件を満たす（たとえば、最小値）割合△ｙｉを算出したときに利用した定数（Ａ，Ｂ，Ｃ，ｂｉｊ，ｂｊｉ，αｉｊ）の組を、各物質についての定数として出力しても良い。
【００７０】
また、本実施の形態のシミュレーション装置は、定数を出力するとともに、当該定数を適宜記憶し、対応する物質の飽和蒸気圧や活量係数を算出しても良い。たとえば、出力される定数を利用し、対応する物質の温度が入力されることに応じて、Ａｎｔｏｉｎｅ式に基づいて当該入力された温度での飽和蒸気圧を算出する。また、出力される定数を利用し、２つの物質の混合溶液についてのモル分率と温度を入力されることにより、ＮＲＴＬ式に基づいて、これらの物質の活量係数を算出する。
【実施例】
【００７１】
本発明のシミュレーション装置による活量係数および飽和蒸気圧の算出の一実施例として、ｎ−ｈｅｘａｎｅ（ＨＡＮ），２−ｐｒｏｐａｎｏｌ（ＩＰＡ），ｔｏｌｕｅｎｅ（ＴＯＬ），ｎ−ｈｅｐｔａｎｅ（ＮＨＰ）の４物質を含む混合溶液を用い、ＮＨＰのみの飽和蒸気圧が既知（Ａｎｔｏｉｎｅ式の定数（Ａ，Ｂ，Ｃ）が既知）であるとして、他の物質のＡｎｔｏｉｎｅ式の定数およびＮＲＴＬ式の定数を求めた。
【００７２】
なお、表１に、４物質のＡｎｔｏｉｎｅ式の定数（Ａ，Ｂ，Ｃ）の文献値を示す。
【００７３】
【表１】
2007212356.tif 000012
【００７４】
また、表２に、４物質についてのＮＲＴＬ式についての２物質間の定数の文献値を示す。なお、表２では、ｂ１２，ｂ２１は、それぞれ上記したｂｉｊ，ｂｊｉを意味し、ｂ１２，ｂ２１の「１」とは表２中の（１）の欄に記載された物質を意味し、ｂ１２，ｂ２１の「２」とは表２中の（２）の欄に記載された物質を意味する。
【００７５】
【表２】
2007212356.tif 000013
【００７６】
まず、混合溶液として、表３〜表５に条件１〜条件９８として示したような、温度および上記の４種類の物質の濃度が異なる９８種類の混合溶液を準備し、これらのそれぞれについての液相および気相のクロマトグラフをシミュレーション装置に入力した。
【００７７】
これに応じて、シミュレーション装置は、この９８種類の混合溶液のそれぞれについて、表３〜表５に示すような液相と気相の組成（モル分率）を算出し、ＲＡＭ１５に記憶する。
【００７８】
【表３】
2007212356.tif 000014
【００７９】
【表４】
2007212356.tif 000015
【００８０】
【表５】
2007212356.tif 000016
【００８１】
表３〜表５において、Ｅは常用対数を示している。つまり、表３中の「２．３６Ｅ−０１」とは、「０．２３６」を意味する。
【００８２】
次に、シミュレーション装置は、ＲＡＭ１５から各条件の液相のモル分率と温度を読出し、ｂｉｊ、ｂｊｉ、αｉｊを適宜読出し、ＮＲＴＬ式に従って、それぞれの混合溶液内の各物質について、暫定的な活量係数を算出し、ＲＡＭ１５に記憶する。これにより、ＲＡＭ１５には、表６に示すようなγ１０１〜γ４９８が記憶される。
【００８３】
【表６】
2007212356.tif 000017
【００８４】
なお、γ１０１〜γ４９８は、それぞれ、活量係数の算出結果である数値を意味する。また、表６では、ＨＡＮについての活量係数がγ１０１〜γ１９８で示され、ＩＰＡについての活量係数がγ２０１〜γ２９８で示され、ＴＯＬについての活量係数がγ３０１〜γ３９８で示され、ＮＨＰについての活量係数がγ４０１〜γ４９８で示されている。
【００８５】
また、γ１０１〜γ４９８の算出の際に読出されるｂｉｊ、ｂｊｉ、αｉｊは、混合溶液内の物質の２物質の組ごとに設定されたものである。つまり、混合溶液４種類の物質が含まれる場合には、２物質の組合わせが６種類考えられるため、ｂｉｊ、ｂｊｉ、αｉｊの値の組は６種類の組合せに対して１組ずつ利用されることになる。
【００８６】
次に、シミュレーション装置は、ＲＡＭ１５から各条件の温度を読出し、Ａ、Ｂ、Ｃを適宜読出し、Ａｎｔｏｉｎｅ式に従って、それぞれの混合溶液内の各物質についての暫定的な飽和蒸気圧を算出し、ＲＡＭ１５に記憶する。これにより、ＲＡＭ１５には、表７に示すようなｐｓａｔ１０１〜ｐｓａｔ４９８が記憶される。
【００８７】
【表７】
2007212356.tif 000018
【００８８】
なお、ｐｓａｔ１０１〜ｐｓａｔ４９８は、それぞれ、飽和蒸気圧の算出結果である数値を意味する。また、表７では、ＨＡＮについての飽和蒸気圧がｐｓａｔ１０１〜ｐｓａｔ１９８で示され、ＩＰＡについての飽和蒸気圧がｐｓａｔ２０１〜ｐｓａｔ２９８で示され、ＴＯＬについての飽和蒸気圧がｐｓａｔ３０１〜ｐｓａｔ３９８で示され、ＮＨＰについての飽和蒸気圧がｐｓａｔ４０１〜ｐｓａｔ４９８で示されている。
【００８９】
次に、シミュレーション装置は、各条件について算出した暫定的な活量係数と飽和蒸気圧を読出し、クロマトグラフから求めた各条件での液相のモル分率を読出し、そして、上記した式（１１）に従って、各条件での各物質の分圧を算出し、ＲＡＭ１５に記憶する。これにより、ＲＡＭ１５には、表８に示すようなｐ１０１〜ｐ４９８が記憶される。
【００９０】
【表８】
2007212356.tif 000019
【００９１】
なお、ｐ１０１〜ｐ４９８は、それぞれ、分圧の算出結果である数値を意味する。また、表８では、ＨＡＮについての分圧がｐ１０１〜ｐ１９８で示され、ＩＰＡについての飽和蒸気圧がｐ２０１〜ｐ２９８で示され、ＴＯＬについての飽和蒸気圧がｐ３０１〜ｐ３９８で示され、ＮＨＰについての飽和蒸気圧がｐ４０１〜ｐ４９８で示されている。
【００９２】
次に、シミュレーション装置は、各条件について、算出した各物質の分圧の総和（全圧ｐＴ）を算出し、ＲＡＭ１５に記憶する。これにより、条件１〜条件９８のそれぞれについての全圧ｐＴとして、表９に示すようなｐＴ１〜ｐＴ９８がＲＡＭ１５に記憶される。
【００９３】
【表９】
2007212356.tif 000020
【００９４】
なお、ｐＴ１〜ｐＴ９８は、それぞれ、全圧の算出結果である数値を意味する。
次に、シミュレーション装置は、上記のｐ１０１〜ｐ４９８とｐＴ１〜ｐＴ９８を読出して、各条件での各物質の気相のモル分率を算出し、ＲＡＭ１５に記憶させる。これにより、ＲＡＭ１５には、表１０に示すようなｙ（ｃａｌ）１０１〜ｙ（ｃａｌ）４９８が記憶される。
【００９５】
【表１０】
2007212356.tif 000021
【００９６】
なお、ｙ（ｃａｌ）１０１〜ｙ（ｃａｌ）４９８は、それぞれ、気相のモル分率の算出結果である数値を意味する。また、表１０では、ＨＡＮについての気相のモル分率がｙ（ｃａｌ）１０１〜ｙ（ｃａｌ）１９８で示され、ＩＰＡについての気相のモル分率がｙ（ｃａｌ）２０１〜ｙ（ｃａｌ）２９８で示され、ＴＯＬについての気相のモル分率がｙ（ｃａｌ）３０１〜ｙ（ｃａｌ）３９８で示され、ＮＨＰについての気相のモル分率がｙ（ｃａｌ）４０１〜ｙ（ｃａｌ）４９８で示されている。
【００９７】
次に、シミュレーション装置は、クロマトグラフから算出した各物質の各条件での気相のモル分率（ｙ（ｅｘｐ）：表３〜表５参照）と上記のｙ（ｃａｌ）１０１〜ｙ（ｃａｌ）４９８とを読出し、以下の式（１２）に従って、各条件での各物質の△ｙを算出する。
【００９８】
【数１０】
2007212356.tif 000022
【００９９】
これにより、ＲＡＭ１５には、表１１に示すような△ｙ１０１〜△ｙ４９８が記憶される。
【０１００】
【表１１】
2007212356.tif 000023
【０１０１】
なお、△ｙ１０１〜△ｙ４９８は、それぞれ、算出結果である数値を意味する。また、表１１では、ＨＡＮについての算出結果が△ｙ１０１〜△ｙ１９８で示され、ＩＰＡについての算出結果が△ｙ２０１〜△ｙ２９８で示され、ＴＯＬについての算出結果が△ｙ３０１〜△ｙ３９８で示され、ＮＨＰについての算出結果が△ｙ４０１〜△ｙ４９８で示されている。
【０１０２】
そして、シミュレーション装置は、△ｙ１０１〜△ｙ４９８の総和Σ△ｙを算出し、所定の条件が成立したか否かを判断する。所定の条件は、適宜ユーザ等によって設定可能なものであると考えられる。なお、たとえば、所定の条件が、総和Σ△ｙが予め定められた設定値以下となることとされれば、シミュレーション装置では、そのような判断がなされる。そして、シミュレーション装置は、所定の条件が成立したΣ△ｙを算出したときに利用していた定数（ｂｉｊ、ｂｊｉ、αｉｊ、Ａ、Ｂ、Ｃ）を、混合溶液内の各物質の活量係数および飽和蒸気圧を求める式の定数であるとして出力する。
【０１０３】
表１２および表１３に、シミュレーション装置が出力した各物質についての定数（Ａ、Ｂ、Ｃ）および各２物質間の定数（ｂｉｊ、ｂｊｉ、αｉｊ）の一例を推算値として示す。
【０１０４】
【表１２】
2007212356.tif 000024
【０１０５】
【表１３】
2007212356.tif 000025
【０１０６】
表１と表１２の対応する各値を比較することにより、また、表２と表１３の対応する各値を比較することにより理解されるように、本発明に従ったシミュレーション装置は、文献値に対して精度の高い推算値を出力できる。
【０１０７】
本実施例では、気液平衡状態にある混合溶液についてその温度と各物質の気相および液相のモル分率が実験データに基づいて求められると、それらを用い、そして、適宜ＮＲＴＬ式の定数（ｂｉｊ、ｂｊｉ、αｉｊ）が定められることによって暫定的に活量係数が算出され、さらに、適宜Ａｎｔｏｉｎｅ式の定数（Ａ，Ｂ，Ｃ）が定められることによって飽和蒸気圧が算出される。そして、これらの算出結果と実験データに基づいて求められた気相および液相のモル分率とによってΣ△ｙが求められて、所定の条件が成立したかどうかの判断がなされる。そして、所定の条件が成立すれば、その時点で利用されていた定数が、各式の定数として出力される。なお、暫定的な活量係数の算出と暫定的な飽和蒸気圧の算出については、たとえば、マイクロソフト（登録商標）社のソフトウェアであるエクセル（登録商標）のソルバーという機能によって実現されるような、利用される定数が適宜変更される回帰計算が行なわれる。
【０１０８】
なお、暫定的な活量係数と飽和蒸気圧を算出することによって行なわれるΣ△ｙの算出は、当該Σ△ｙが所定の設定値以下となる回帰条件が成立するまで実行されても良いし、または、ユーザ等によって指定された回数だけ実行されても良い。指定された回数実行される場合には、シミュレーション装置は、各回で算出したΣ△ｙをＲＡＭ１５に蓄積し、そして、Σ△ｙの最小値が算出されたときに各式に利用した定数を、混合溶液の各物質についての定数として出力する。
【０１０９】
以上説明した実施例では、すべての条件のすべての物質について、活量係数が算出された後で飽和蒸気圧が算出されたが、本発明はこれに限定されない。これらの値は、△ｙｉを算出するために算出されるものであるため、算出される順番は問われない。つまり、たとえば、条件の番号順に、活量係数と飽和蒸気圧が算出されても良い。
【０１１０】
以上説明した実施例では、混合溶液を構成する複数の物質の中の、少なくとも１つの物質は、分子量および飽和蒸気圧を算出するための式に用いる定数（Ａ，Ｂ，Ｃ）が既知である必要がある。そして、上記の計算が行なわれる際には、この物質についての定数（Ａ，Ｂ，Ｃ）が固定値として代入されている。
【０１１１】
混合溶液中に、分子量および飽和蒸気圧が既知の物質が含まれていない場合には、そのような物質を混合溶液中に添加し、当該物質が含まれた状態で、液相および気相の各物質のモル分率を算出し、上記の計算を行なえば良い。なお、添加される物質は、混合溶液中の各物質が著しい非理想性を示さないために、添加前に混合溶液を構成する物質と化学構造が類似していることが好ましい。また、添加される物質は、計算の精度を上げるために、添加前に混合溶液を構成する物質と蒸気圧が近いことが好ましい。
【０１１２】
今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【図面の簡単な説明】
【０１１３】
【図１】本実施の形態において物性値の算出を行なうための装置のブロック図である。
【図２】図１のＣＰＵが実行する物性値算出処理のフローチャートである。
【符号の説明】
【０１１４】
１１  ＣＰＵ、１２  Ｉ／Ｏ、１３  ＨＤ、１４  ディスクドライブ、１５  ＲＡＭ、１６  入力部、１８  表示部、１００  ディスク。
【図１】
2007212356.tif 000030
【図２】
2007212356.tif 000031