JP6725112B6 - 多成分系溶液の性状推定方法、装置、及びプログラム - Google Patents

多成分系溶液の性状推定方法、装置、及びプログラム Download PDF

Info

Publication number
JP6725112B6
JP6725112B6 JP2016150586A JP2016150586A JP6725112B6 JP 6725112 B6 JP6725112 B6 JP 6725112B6 JP 2016150586 A JP2016150586 A JP 2016150586A JP 2016150586 A JP2016150586 A JP 2016150586A JP 6725112 B6 JP6725112 B6 JP 6725112B6
Authority
JP
Japan
Prior art keywords
component
liquid phase
solution
aggregation
fraction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016150586A
Other languages
English (en)
Other versions
JP6725112B2 (ja
JP2018018448A (ja
Inventor
正人 森本
正人 森本
隆三 田中
隆三 田中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chiyoda Corp
Cosmo Oil Co Ltd
Idemitsu Kosan Co Ltd
National Institute of Advanced Industrial Science and Technology AIST
Japan Petroleum Energy Center JPEC
Eneos Corp
Original Assignee
Chiyoda Corp
Cosmo Oil Co Ltd
Idemitsu Kosan Co Ltd
National Institute of Advanced Industrial Science and Technology AIST
Japan Petroleum Energy Center JPEC
Eneos Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chiyoda Corp, Cosmo Oil Co Ltd, Idemitsu Kosan Co Ltd, National Institute of Advanced Industrial Science and Technology AIST, Japan Petroleum Energy Center JPEC, Eneos Corp filed Critical Chiyoda Corp
Priority to JP2016150586A priority Critical patent/JP6725112B6/ja
Publication of JP2018018448A publication Critical patent/JP2018018448A/ja
Application granted granted Critical
Publication of JP6725112B2 publication Critical patent/JP6725112B2/ja
Publication of JP6725112B6 publication Critical patent/JP6725112B6/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Description

本発明は、多成分系溶液の性状推定方法、装置及びプログラムに係り、より詳細には、新たに想定した多成分系凝集モデルに基づいて、重質油をはじめとする種々の多成分系溶液における各成分の溶解、凝集及び析出等の性状を定量的に推定するための方法、装置、及びその方法をコンピュータに実行させるコンピュータプログラムに関する。
石油精製の過程においては、重質油の分解処理をはじめとする様々な場面で、コークの前駆体であるアスファルテンの凝集を緩和、抑制することが必要である。そのため、アスファルテンの凝集を緩和、抑制する溶媒（溶剤）を原料油に予め混合しておくことが通常行われているが、混合される溶媒の種類及び量は、永年蓄積されてきた知見に基づいて経験的に選定されているのが現状である。
一方、本出願に係る発明者らは、重質油に含まれるアスファルテンの凝集を抑制、緩和する技術について鋭意研究した結果、アスファルテンと溶媒とのハンセン溶解度指数値の差（Δδ）と、アスファルテンの凝集度との間に一定の関係があることを見出し、その差（Δδ）を指標として溶媒を選定する方法を開発した。かかる選定方法は、下記の特許文献１に記載されている。
特開２０１４−２１８６４３号公報
ところで、溶媒と混合した重質油においては、液相のハンセン溶解度指数値は、溶媒の値そのものではなく、溶媒の値と重質油中の液相成分の値とに基づいて決まると考えられる。このため、重質油に含まれる各成分のうち、どの成分が溶解して液相をなすのか、それにより、液相のハンセン溶解度指数値はどう変化するのか、その結果、どの成分が凝集又は析出するのかについて、把握することが重要となる。
しかしながら、複数の成分からなる溶液である多成分系溶液の液相、凝集相及び固相の性状、その相変化のメカニズム、特に重質油中のアスファルテンの溶解、凝集及び析出の挙動については、産業上大きなテーマであるにもかかわらず、その解明はほとんどなされていない。
本発明は、かかる事情に鑑みてなされたものであり、多成分系溶液中の各成分の性状を推定することができる方法、装置及びプログラム、さらには、かかる方法を用いた多成分系溶液の処理方法を提供することを目的としている。
本出願に係る発明者らは、多成分系溶液中における溶質の凝集のメカニズム、とりわけ重質油中におけるアスファルテンの凝集のメカニズムについて鋭意検討した結果、アスファルテンの凝集性状を適切に表しているであろうと考えられる多成分系凝集モデル（Multi-Component Aggregation Model：ＭＣＡＭ）を想定した。その結果、かかるモデルに基づけば、重質油中のアスファルテンの溶解、凝集及び析出といった性状を定量的に推定することができることを見出し、本発明を完成させた。
本発明の多成分系溶液の性状推定方法は、コンピュータによる多成分系溶液の性状推定方法であって、
（１）前記多成分系溶液を構成する各成分の分率、融点及びハンセン溶解度指数値を取得するステップと、
（２）前記多成分系溶液を構成する各成分のうちの所望の温度未満の融点を有する成分を液相成分として分類し、前記所望の温度以上の融点を有する成分を非液相成分として分類するステップと、
（３）前記液相成分として分類された各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、液相全体の平均ハンセン溶解度指数値を算出するステップと、
（４）液相全体の前記平均ハンセン溶解度指数値と、非液相成分における各成分のハンセン溶解度指数値との差を算出するステップと、
（５）非液相成分における各成分を、前記差に基づいて、液相成分又は非液相成分として再分類し、液相成分として再分類された各成分を非液相成分から液相成分に編入して、液相成分及び非液相成分を更新するステップと、
（６）更新後の液相成分における各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、更新後の液相全体の平均ハンセン溶解度指数値を算出するステップと、
（７）前記ステップ（４）〜（６）を、前記ステップ（５）において液相成分として再分類される非液相成分がなくなる最終段階まで繰り返すステップと
を含むことを特徴としている。
前記「多成分系溶液を構成する各成分」において、「成分」とは、「溶液をある特定の物理的又は化学的性状を基準として括った塊」、言い換えれば、「ある特定の物理的又は化学的性状を基準として分画された分画物(フラクション)」という意味と捉えることができる。特定の物理的又は化学的性状を基準として括る方法としては、例えば、蒸留試験における沸点範囲を特定して、その温度範囲にあるものを一つの成分として分画する方法等が考えられる。この場合、溶液は「分画物(フラクション)の集合体」ということになる。
或いは、「成分」を「同一の分子種に属すると認められる分子の集合体」と捉えることも可能である。即ち、「同一の分子種に属すると認められる分子を集めて一括りにしたもの」と捉えることである。ここで「同一の分子種」とは、「分子構造を完全に特定し、その上で同一であると認められる分子の一括り」という意味と考えてもよく、また、「ある任意に設定した特定の条件を満たす分子の一括り」という意味で捉えてもよい。「ある任意に設定した特定の条件」とは、例えば、構造上同じ官能基を有していることのみを条件とする場合や、ある構造の芳香環を有していることのみを条件とする場合など、条件は任意に設定してもよいという意味である。
「溶液を構成する各成分」という文言において、「成分」という語を、前記後者の「同一の分子種に属すると認められる分子の集合体」という意味で用いる場合、「構成する各成分」とは、厳密に当該溶液中に存在するすべての分子種が構成の対象となるというわけではなく、溶液中において一定の存在量(存在割合)以上を持つ分子種のみを対象とすると考えても良い。当該溶液中に存在しているできる限り多くの分子種を対象とすることが望ましいが、例えば、重質油には何十万という種類の分子種が存在しているとも言われているので、その場合は、微量しか存在していないような分子種は無視してもよい。試料とする溶液を前もって成分分析し、各分子種の存在量(存在割合)を以て、対象とする分子種の選定基準にしてもよい。
さらに、「分率」については、重量分率、容量分率又はモル分率等、存在割合を示すものであれば何でもよく、いずれをも含む概念である。液相全体の平均ハンセン溶解度指数値を算出する場合は、好ましくは容量分率が用いられ、各成分の当該液相における容量分率で重み付けした加重平均値として算出される。
本発明によれば、実施形態の欄で後述するように、多成分系凝集モデルに基づいて、多成分系溶液の性状を推定することができる。具体的には、多成分系溶液の性状として、液相の成分、及び液相の平均ハンセン溶解度指数値を推定することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、最終段階での更新後の液相成分における各成分の分率の合計を液相分率として算出するステップを更に含む。
これにより、多成分系溶液の性状として、液相の量（溶液全体に対する分率）を推定することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、所望の温度における前記最終段階での更新後の非液相成分における各成分の凝集度を、液相全体の前記平均ハンセン溶解度指数値と前記非液相成分における各成分のハンセン溶解度指数値との差及び前記最終段階での更新後の非液相成分における各成分の濃度に基づいて算出するステップを更に有する。
これにより、多成分系溶液の性状として、多成分系溶液中における非液相の各成分の凝集度をそれぞれ推定することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、凝集度（Ｄ）を下記の式（Ａ）により算出する。
Ｄ(ｐ，ｑ) = ＭＡＳ (Ｋ０ +Ｋ１ p +Ｋ２ ｑ +Ｋ３ ｐ２ +Ｋ４ ｐｑ +Ｋ５ ｑ２+Ｋ６ｐ３ +Ｋ７ｐ２ｑ +Ｋ８ｐｑ２ +Ｋ９ ｑ３) ・・・（Ａ）
ここで、ｐは、前記所望の温度Ｔが、Ｔ≦１５０℃のときに、ｐ = (Ｌ０(Ｔ - 25) + Ｌ１)ＲＥＤｇ、前記所望の温度Ｔが、１５０℃＜Ｔ≦２００℃のときに、ｐ = (Ｌ０(150 - 25) + Ｌ１)ＲＥＤｇ、前記所望の温度Ｔが、２００℃＜Ｔのときに、ｐ = (Ｌ０(Ｔ - 25) + Ｌ２)ＲＥＤｇで表され、
Ｌ０、Ｌ１及びＬ２は、係数であり、
ＲＥＤｇは、ＲＥＤ≧０．３のときに、ＲＥＤｇ＝ＲＥＤ、ＲＥＤ＜０．３のときに、ＲＥＤｇ＝０．３と表され、
ＲＥＤは、ＲＥＤ＝Δδ／Ｒ０で表され、
Δδは、液相全体の前記平均ハンセン溶解度指数値と前記非液相成分における各成分のハンセン溶解度指数値との差であり、
Ｒ０は、非液相成分における各成分の定数であり、
ｑは、ｑ＝ｌｏｇＣで表され、
Ｃは、非液相成分における各成分の濃度であり、
ＭＡＳ、及びＫ０〜Ｋ９は、係数である。
これにより、上記の（Ａ）式を用いて、多成分系溶液中における非液相の各成分の凝集度をそれぞれ推定することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、最終段階での更新後の非液相成分における成分のうち、前記凝集度が所定の閾値未満の成分を凝集相成分として分類し、前記凝集度が前記所定の閾値以上の成分を固相成分として分類するステップを更に有する。
これにより、多成分系溶液の性状として、多成分系溶液中の凝集相の成分、固相の成分を推定することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、凝集相成分として分類された各成分の分率の合計を凝集相分率として算出するステップと、前記固相成分として分類された各成分の分率の合計を固相分率として算出するステップとを更に有する。
これにより、多成分系溶液の性状として、多成分系溶液中の凝集相の量（溶液全体に対する分率）と、固相の量（溶液全体に対する分率）を推定することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、前記凝集相成分として分類された各成分の凝集度の和を当該成分の数で除した値として、凝集相全体の平均凝集度を算出するステップを更に有する。
これにより、多成分系溶液の性状として、多成分系溶液中の凝集相の平均凝集度を推定することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、多成分系溶液には、所望の種類の溶媒が所望の分率で混合されている。
本発明は、溶剤等の溶媒が混合された多成分系溶液にも適用することができる。
また、本発明の多成分系溶液の性状推定方法において、好ましくは、多成分系溶液は、重質油、又は重質油と溶媒との混合液である。
これにより、重質油の性状を推定し、アスファルテンの凝集を制御して、重出油の処理に寄与することができる。
また、本発明の多成分系溶液の性状推定装置は、多成分系溶液の性状推定装置であって、
多成分系溶液を構成する各成分の分率、融点及びハンセン溶解度指数値を取得する成分情報取得部と、多成分系溶液を構成する各成分のうちの所望の温度未満の融点を有する成分を液相成分として分類し、前記所望の温度以上の融点を有する成分を非液相成分として分類する初期分類部と、液相の性状を推定する液相演算部とを備え、
液相演算部は、前記液相成分として分類された各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、液相全体の平均ハンセン溶解度指数値を算出し、液相全体の前記平均ハンセン溶解度指数値と、非液相成分における各成分のハンセン溶解度指数値との差を算出し、非液相成分における各成分を、前記差に基づいて、液相成分又は非液相成分として再分類し、液相成分として再分類された各成分を非液相成分から液相成分に編入して、液相成分及び非液相成分を更新し、更新後の液相成分における各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、更新後の液相全体の平均ハンセン溶解度指数値を算出し、液相成分として再分類される非液相成分がなくなる最終段階まで、平均ハンセン溶解度指数値、液相成分、及び非液相成分の更新を繰り返すことを特徴としている。
本発明によれば、実施形態の欄で後述する多成分系凝集モデルに基づいて、多成分系溶液の性状を推定することができる。具体的には、多成分系溶液の性状として、液相の成分、及び液相の平均ハンセン溶解度指数値を推定することができる。
また、本発明の多成分系溶液の性状推定装置において、好ましくは、非液相の性状を推定する非液相演算部を更に備え、前記非液相演算部は、所望の温度における前記最終段階での更新後の非液相成分における各成分の凝集度を、液相全体の前記平均ハンセン溶解度指数値と前記非液相成分における各成分のハンセン溶解度指数値との差及び前記最終段階での更新後の非液相成分における各成分の濃度に基づいて算出する。
これにより、多成分系溶液の性状として、多成分系溶液中における非液相の各成分の凝集度をそれぞれ推定することができる。
また、本発明のプログラムは、上記の多成分系溶液の性状推定方法を実行させることを特徴としている。
本発明のプログラムによれば、コンピュータを上記の多成分系溶液の性状推定装置として機能させることができる。これにより、多成分系溶液中の各成分の性状を推定することができる。
また、本発明の多成分系溶液の処理方法は、上記の多成分系溶液の性状推定方法により、原多成分系溶液の性状を推定し、さらに当該性状推定結果に基づいて、原多成分系溶液に所望の種類の溶媒を所望の分率で混合した場合、原多成分系溶液の温度を所望の温度に変更した場合、又は原多成分系溶液に所望の種類の溶媒を所望の分率で混合しかつ原多成分系溶液の温度を所望の温度に変更した場合における新多成分系溶液の性状を予測し、その予測を基に、その所望の種類の溶媒をその所望の分率で原多成分系溶液に混合する、又は原多成分系溶液の温度をその所望の温度に変更する、又はその所望の種類の溶媒をその所望の分率で原多成分系溶液に混合しかつ原多成分系溶液の温度をその所望の温度に変更することを特徴としている。
本発明の多成分溶液系の処理方法によれば、多成分系溶液に溶媒を混合したり、多成分系溶液の温度を変化させたりしたときの多成分系溶液の性状を予想することができ、かかる予想に基づいて、多成分系溶液の処理を行うことが可能となる。
また、本発明の多成分系溶液の処理方法において、好ましくは、その所望の種類の溶媒及びその所望の分率、及び／又はその所望の温度は、予測された新多成分系溶液の性状が、推定された原多成分系溶液における凝集相成分及び固相成分の全部又は一部が溶解し、又は、凝集相成分及び固相成分の凝集又は析出が抑制された性状となる溶媒及び分率、及び／又は温度である。
これにより、溶媒種を適切に選択する、新たな溶媒を追加する、又は、溶液の温度を変更する措置を講じて、溶液中における凝集相及び固相の全部又は一部を溶解し、又は溶質の凝集若しくは析出を抑制することができる。
また、本発明の多成分系溶液の処理方法において、好ましくは、その所望の種類の溶媒及びその所望の分率、及び／又はその所望の温度は、予測された新多成分系溶液の性状が、推定された原多成分系溶液における非固相成分のうちの一種類以上の成分が析出した性状となる溶媒及び分率、及び／又は温度である。
これにより、溶媒種を適切に選択する、新たな溶媒を追加する、又は、溶液の温度を変更する措置を講じて、一種以上の成分を析出させることができる。
本発明によれば、多成分系溶液中の各成分の性状を推定することができる。例えば、本発明によれば、多成分系溶液中における液相、凝集相及び固相に関し、各相の成分及びその量（分率）を推定し、また、凝集相における各成分の凝集度及び凝集相の平均凝集度を算出することもできる。このため、例えば、溶媒種を適切に選択する、或いは新たな溶媒を追加する、溶液の温度を変更する等の措置を適切に講ずれば、溶液中における凝集相及び固相の全部又は一部を溶解し、又は溶質の凝集若しくは析出を抑制すること、或いは一種以上の成分を析出させることが可能となる。
本発明の実施形態による多成分系溶液の性状推定装置の機能ブロック図である。 本発明の実施形態による多成分系溶液の性状推定方法のフローチャートである。 ＨＳＰ値の模式図である。 無溶媒でのモデル重質油の液相、凝集相及び固相の量と凝集相の平均凝集度を表すグラフである。 トルエン溶媒でのモデル重質油の液相、凝集相及び固相の量と凝集相の平均凝集度を表すグラフである。 トルエンとブロモベンゼンの等重量被混合溶媒でのモデル重質油の液相、凝集相及び固相の量と凝集相の平均凝集度を表すグラフである。
（１．多成分凝集モデル）
本発明の実施形態の説明に先立ち、本発明が立脚する多成分凝集モデルについて説明する。以下、「ハンセン溶解度指数値」を「ＨＳＰ値」と記すことがある。
重質油を例にとって、多成分系凝集モデル（Multi-Component Aggregation Model：ＭＣＡＭ）を前記ステップ(1)〜(7)に合わせて説明する。ここで、重質油というものを、説明の簡便上、Ａ分子、Ｂ分子、Ｃ分子、Ｄ分子及びＥ分子から構成されているものと仮定する。
(1)まず、各成分の分率、融点及びハンセン溶解度指数値を取得する。(2)次にある所望の温度Ｔを設定したとする。Ａ分子、Ｂ分子、Ｃ分子、Ｄ分子及びＥ分子のうち、温度Ｔよりも融点が低いＡ分子及びＢ分子は液体となるので液相を形成し、Ｃ分子、Ｄ分子及びＥ分子は非液相を形成する。(3)これにより、液相のＨＳＰ値は、Ａ分子のＨＳＰ値とＢ分子のＨＳＰ値をそれぞれの存在割合(好ましくは容量分率)により加重平均した値となる。(4)この時点までは、Ｃ分子、Ｄ分子及びＥ分子は非液相成分であったが、液相のＨＳＰ値とＣ分子、Ｄ分子及びＥ分子の各々のＨＳＰ値との差を見ると、(5)液相のＨＳＰ値との差が極めて小さいＣ分子は液相に移ることになる。(6)その結果、液相のＨＳＰ値は、Ａ分子、Ｂ分子及びＣ分子の各々のＨＳＰ値の各々の存在割合(好ましくは容量分率) により加重平均した値に変化する。(7)その結果、Ｄ分子は、当該更新後の液相とのＨＳＰ値差が小さくなったため、固相を維持できなくなり、凝集相となる。一方、Ｅ分子は、液相のＨＳＰ値が変化したにもかかわらずＨＳＰ値差が大きいため、依然として、固相を維持している。
このようにして、液相の成分の再分類が行われ、最終的に、全成分（Ａ分子、Ｂ分子、Ｃ分子、Ｄ分子及びＥ分子）が、液相成分、凝集した凝集成分、析出した固相成分に分けられた状態となる。
現実の重質油というものは、何十万種という成分から成っているものであるので、液相のＨＳＰ値及び各成分のＨＳＰ値に基づいて液相成分と非液相成分の再分類が次々と行われ、最終的に、液相、凝集相及び固相に分類された状態となる。
因みに、重質油中に存在するアスファルテンというのは、アスファルテン分子がアボガドロ数(6×10の23乗)以上塊集すると、目に見える程度の大きさになり固体として析出するが、塊集度合がそれよりはるかに少ない場合は、重質油中で分散・浮遊した状態にあると推定できる。この状態にあるものを「凝集体」と名付け、その凝集度合を小さくすることを「凝集の緩和、抑制」と呼ぶ。
（２．装置）
次に、図1を参照して、本発明の多成分系溶液の性状推定装置の実施形態を説明する。図１は、実施形態の多成分系溶液の性状推定装置の機能ブロック図である。コンピュータに本発明のプログラムを実行させることにより、コンピュータが多成分系溶液の性状推定装置として機能する。
なお、図１では、情報の入力及び出力を行うインタフェースの図示を省略している。
本装置は、演算装置１と記憶部２とを備えている。演算装置１は、１つのＣＰＵで構成してもよいし、通信回線を介して互いに接続された複数の演算装置で構成されてもよい。また、記憶部２は、演算装置１に内蔵されていてもよいし、外部装置であってもよいし、通信回線を介して接続された記憶装置であってもよい。
本演算装置１は、成分情報取得部１０と、初期分類部２０と、液相演算部３０と、非液相演算部４０とを有している。
（２−１．成分情報取得部）
成分情報取得部１０は、対象とする多成分系溶液を構成する各成分について、その分率、融点、及びＨＳＰ値を取得する。これらの成分の情報は、多成分系溶液についての情報がデータベースとして格納された記憶部２から取得するとよい。
データベースにこれらの成分の情報が格納されていない場合には、成分情報算出部１１によって、各成分の必要なパラメータを推算するとよい。
多成分系溶液を構成している成分の融点とＨＳＰ値を推算する方法の一例として、「分子組成(分子構造)に関する情報を基に行う方法」を挙げることができる。
（１）この方法では、先ず、試料である溶液を構成している各分子種につき、各々の分子種の分子組成(分子構造)に関する情報を得ることが必要である。ここで、溶液を構成している分子種とは、当該溶液中に存在している厳密にすべての分子種を指すというわけではなく、溶液中において一定の存在量(存在割合)以上を持つ分子種を指すと考えても良い。当該溶液中に存在しているできる限り多くの分子種を対象とすることが望ましいが、微量しか存在していないような分子種は無視してもよい。試料とする溶液を前もって成分分析し、各分子種の存在量(存在割合)を以て、対象とする分子種の選定基準にしてもよい。
あるいはまた、前述のように、「成分」を「溶液をある特定の物理的又は化学的性状を基準として括った塊」、言い換えれば、「ある特定の物理的又は化学的性状を基準として分画された分画物(フラクション)」という意味で用いる場合には、この「分子組成(分子構造)に関する情報を基に行う方法」は、次のようにして適用することが可能である。
即ち、「ある特定の物理的又は化学的性状を基準として分画された分画物(フラクション)」の各々について、ＮＭＲ、元素分析、質量スペクトル等を測定することにより、公知の方法を用いて、その分画物(フラクション)の「平均分子構造」を得ることができる。こうして得られた「平均分子構造」を用いれば、この方法を適用することができる。
（２）次に、得られた各々の分子種の分子組成(分子構造)に関する情報に基づいて、各々の分子種の融点を推定する。これには、原子団寄与法(Marrero Gani法)を用いることができる。即ち、分子種について、その分子構造を形成している「基」に分解し、各々の基が持つ固有のパラメータ値からその分子種の融点を算出する。
前記のように「平均分子構造」を用いる場合も、当該構造を基に原子団寄与法を用いて、当該「ある特定の物理的又は化学的性状を基準として分画された分画物(フラクション)」の融点を算出する。
（３）また、各々の分子種の分子組成(分子構造)に関する情報に基づいて、各々の分子種のＨＳＰ値を推定する。前記のように「平均分子構造」を用いる場合も、同様である。ＨＳＰ値の推定も、原子団寄与法にて行うことができる。即ち、分子種についてその分子構造を形成している基に分解し、各々の基が持つ固有のパラメータ値(Fd値, Fp値, Eh値, Vc値)を用いて、その分子種のＨＳＰ値を算出する。
原子団寄与法としては、例えば、D.W.van Krevelen , K.te Nijenhuis著「Properties of Polymers(4ed.2009)」に記載のKrevelen & Hoftyzerの方法やEmmanuel Stefanis, Costas Panayiotou著「Prediction of Hansen Solubility Parameters with a New Group-Contribution Method」に記載のStefanis & Panayiotouの方法等を用いることができる。
これらについては、分子構造から推算するプログラムを利用することもできる。このようなプログラムとして、例えば、ＨＳＰiP（http：//www.hansen-solubility．Com/）がある。
さらには、これらの諸方法で得た値を参考にして、適宜修正を加えて得た値を用いることもできる。
ＨＳＰ値は、文献（Hansen，C.M.，Hansen solubility parameters:a user’s handbook. ２nd ed．; CRC： Boca Raton； London，2007）に、1200種類を超える物質のＨＳＰ値が報告されているので、その文献値を使用することもできる。
Krevelen & Hoftyzerの方法の概略は、次のとおりである。
ｉ．まず、物質のＨＳＰ値(δt)は、次の式で求められることが広く知られている。
図３に、ＨＳＰ値の模式図を示す。ＨＳＰ値とは、ある物質がある溶媒にどのくらい溶けるのかを示す溶解性の指標であり、溶解性を［分散項（δｄ）、分極項（δｐ）、水素結合項（δｈ）］のベクトルで表すものである。
ii．次に、δd, δp, δhは、Krevelen & Hoftyzerの原子団寄与法によると、次の式で求められる。
ここで、Fd値((MJ/m3)1/2・mol-1)とは、分散力に起因する各基(原子団)のモル引力定数であり、Fp値((MJ/m3)1/2・mol-1)とは、双極子間力に起因する各基(原子団)のモル引力定数であり、Eh値(J・mol-1)とは、各基(原子団)の水素結合エネルギーである。
また、Vはモル体積(cm3・mol-1)であり、RheineckおよびLinが提案した以下の式(3)により求めることができる。
ここでＶｃは各基(原子団)のモル体積である。
D.W.van Krevelen , K.te Nijenhuis著「Properties of Polymers(4ed.2009)」の文献によれば、多くの基について、Fd値, Fp値, Eh値, Vc値が示されているので、当該文献にて値が記載されている基については、その値を用いればよい。(D.W.van Krevelen , K.te Nijenhuis著「Properties of Polymers(4ed.2009)」195〜197ページ及び215ページ)。値が記載されていない基については、構造的に近似する他の基の情報を用いて推定した値を用いるなどを行えばよい。
iii．以上により、Krevelen & Hoftyzer法によりＨＳＰ値(δt)を算出することができる。
（２−２．初期分類部）
初期分類部２０は、多成分系溶液を構成する各成分のうちの所望の温度未満の融点を有する成分を液相成分として分類し、所望の温度以上の融点を有する成分を非液相成分として分類する。すなわち、溶媒の融点以上のある任意の温度以上において、その温度における「液相」の量及び組成を求める。融点がその温度より低い成分は、液相に存在する成分となる。このときの「液相」の量及び成分が求まる。
（２−３．液相演算部）
液相演算部３０は、液相の性状を推定するために、平均ＨＳＰ算出部３１と、Δδ（ＨＳＰ値差）算出部３２と、再分類部３３と、液相成分情報算出部３４とを備えている。
平均ＨＳＰ算出部３１は、液相全体の平均ＨＳＰ値を算出する。ここで、液相全体の平均ＨＳＰ値は、当該液相成分における各成分のＨＳＰ値を各成分の当該液相における分率、好ましくは容量分率で重み付けした加重平均値として算出されるものである。
ＨＳＰ値差（Δδ）算出部３２は、液相全体の平均ＨＳＰ値と、非液相成分における各成分のＨＳＰ値との差（Δδ）を算出する。
再分類部３３は、非液相成分における各成分を、差（Δδ）に基づいて、液相成分と非液相成分とに再分類し、液相成分として再分類された各成分を非液相成分から液相成分に編入して、液相成分及び非液相成分を更新する。
再分類部３３は、溶解する成分があればそれを液相に加えて液相全体のＨＳＰ値を再計算する。
平均ＨＳＰ算出部３１は、更新後の液相全体の平均ＨＳＰ値を算出する。ここで、更新後の液相全体の平均ＨＳＰ値は、更新後の液相成分における各成分のＨＳＰ値を各成分の当該液相における分率、好ましくは容量分率で重み付けした加重平均値として算出されるものである。
そして、液相成分に再分類される非液相成分がなくなる最終段階まで、平均ＨＳＰ値、液相成分及び非液相成分（凝集相、固相）の更新を繰り返す。
さらに、液相情報算出部３４は、最終段階での更新後の液相成分の分率の合計を液相分率として算出する。
（２−４．非液相演算部）
非液相演算部４０は、非液相の性状を推定するために、凝集度算出部４１、凝集相、固相分類部４２、凝集相情報算出部４３、及び固相情報算出部４４を有する。非液相演算部４０は、非液相の性状として、例えば、凝集相の量、成分、凝集している成分それぞれの凝集度及び凝集相の平均凝集度並びに固相の量及び組成を決定する。
凝集度算出部４１は、所望の温度における最終段階での更新後の非液相成分における各成分の凝集度を、液相全体の平均ＨＳＰ値と前記非液相成分における各成分のＨＳＰ値との差及び最終段階での更新後の非液相成分における各成分の濃度Ｃに基づいて算出する。具体的には、以下のようにして分類する。
最終的に液相に溶解しなかった非液相成分における各成分について、そのＨＳＰ値と液相全体のＨＳＰ値に基づいてそれぞれの凝集度を決定する。凝集している成分それぞれの凝集度Ｄは、液相のＨＳＰ値、凝集している成分のＨＳＰ値、凝集している成分の濃度及び場の温度を変数とする関数式(Ａ)により、算出することができる。
Ｄ(ｐ，ｑ) = ＭＡＳ (Ｋ０ +Ｋ１ p +Ｋ２ ｑ +Ｋ３ ｐ２ +Ｋ４ ｐｑ +Ｋ５ ｑ２+Ｋ６ｐ３ +Ｋ７ｐ２ｑ +Ｋ８ｐｑ２ +Ｋ９ ｑ３) ・・・（Ａ）
ここで、ｐは、前記所望の温度Ｔが、Ｔ≦１５０℃のときに、ｐ = (Ｌ０(Ｔ - 25) + Ｌ１)ＲＥＤｇ、前記所望の温度Ｔが、１５０℃＜Ｔ≦２００℃のときに、ｐ = (Ｌ０(150 - 25) + Ｌ１)ＲＥＤｇ、前記所望の温度Ｔが、２００℃＜Ｔのときに、ｐ = (Ｌ０(Ｔ - 25) + Ｌ２)ＲＥＤｇで表される。
ＲＥＤｇは、ＲＥＤ≧０．３のときに、ＲＥＤｇ＝ＲＥＤ、ＲＥＤ＜０．３のときに、ＲＥＤｇ＝０．３と表され、ＲＥＤは、ＲＥＤ＝Δδ／Ｒ０で表され、Δδは、液相全体の前記平均ＨＳＰ値と前記非液相成分における各成分のＨＳＰ値との差であり、Ｒ０は、非液相成分における各成分ごとの定数である。
Ｌ０、Ｌ１及びＬ２は、経験的に得た係数であり、下記の定数値を有する。
L0 ＝−0.0031262、 L1 ＝ 1.07815、 L2 ＝ 1.15631
ｑは、ｑ＝ｌｏｇＣで表され、Ｃは、非液相成分における凝集している当該成分の濃度である。
ＭＡＳは、成分種により定まった定数であり、例えば、多成分系溶液の凝集相成分及び固相成分がアスファルテンの場合、以下のとおりである。カナダ産オイルサンド系アスファルテン(CaAs)：1.319、中東産アスファルテン1 (ArAs1)：1.000、中東産アスファルテン2 (ArAs2)：1.136である。
Ｋ０〜Ｋ９は、経験的に得た係数であり、以下の定数値を有する。
K0＝−1.26929、K1＝ 9.42231、K2＝ 0.363439、K3＝−11.1925、K4＝ 0.093622、K5＝−0.15436、K6＝ 5.337433、K7＝−0.20868、K8＝ 0.077223、K9＝ 0.019492である。
以上より、ある温度において、ある溶液中においてある成分が凝集している場合、その凝集している成分の凝集度Ｄの値を算出することができる。
尚、上記において数値で示したＬ０、Ｌ１、Ｌ２、ＭＡＳ及びＫ０〜Ｋ９等の値は、対象により種々の数値を採り得るものであり、上記の数値に限定されるものではない。
凝集相、固相分類部４２は、最終段階での更新後の非液相成分のうち、凝集度が所定の閾値未満の成分を凝集相成分に分類し、凝集度が所定の閾値以上の成分を固相成分に分類する。すなわち、凝集度が凝集レベルにある成分を凝集相成分とし、析出レベルにある成分を固相成分とする。ここで、「凝集レベルにある」とは、概念的には、凝集粒子の大きさが数百nm以下で液中に分散していることをいい、「析出レベルにある」とは、凝集粒子の大きさがサブミクロン以上で液中に分散できず沈殿していることと考えられる。凝集度Ｄ≧５であるとき、おおむね、その成分種は析出すると判断できるが、この閾値は、成分種により変化しうるものである。
凝集相情報算出部４３は、凝集相成分として分類された各成分の量（溶液全体に対する分率）の合計を凝集相分率として算出する。さらに、凝集相情報算出部４３は、凝集相成分として分類された各成分の凝集度の和を当該成分の数で除した値を凝集相全体の平均凝集度として算出する。
固相情報算出部４４は、固相成分として分類された各成分の量（溶液全体に対する分率）の合計を固相分率として算出する。
（３．方法）
次に、図２を参照して、本発明の多成分系溶液の性状推定方法の実施形態を説明する。図２は、実施形態の多成分系溶液の性状推定方法のフローチャートである。本方法は、本発明のプログラムによりコンピュータによって実行される。
本実施形態では、多成分系溶液として、以下に示す１１種類の成分（分子種Ｏ−０１〜Ｏ−１１）を含むモデル重質油を用いて、このモデル重質油の性状（溶解、凝集及び析出性状）をコンピューターにより推定する例を説明する。
（Ｓ１：成分情報取得）
まず、モデル重質油の性状推定にあたり、先ず、モデル重質油を構成する各成分の分率、融点及びＨＳＰ値(δt)を取得する（Ｓ１）。
表１に、モデル重質油の分子組成(分率)、推算した各分子の融点及びＨＳＰ値(δt)を示す。
上記の表１に示した各分子の融点及びＨＳＰ値(δt)は、データベースから既知のものを取得してもよいし、原子団寄与法により、コンピュータを用いて算出してもよい。
以下に、ＨＳＰ値(δt)及び融点を原子団寄与法により算出する一例を説明する。
ｉ．ＨＳＰ値の求め方
（１）「Ｏ−０１」分子の場合
前出のKrevelen & Hoftyzerの文献において、分子を形成している基について示されているFd値, Fp値, Eh値及びVc値の数値を用いて、Krevelen & Hoftyzerの方法により算出できる。
即ち、「Ｏ−０１」分子の場合、基は、「ＣＨ3−」が２個と「−ＣＨ2−」が28個とからなっている。前出のKrevelen & Hoftyzerの文献によると、「ＣＨ3−」のFd値は420、Vc値は33.5と記載されており、「−ＣＨ2−」のFd値は270、Vc値は16.1と記載されているのでこれらの値を用い、前出の式(3)及び式(2)より、「O-01」のδdは、{(420×2)＋(270×28)}／{(33.5×2)＋(16.1×28)}＝16.22となる。δp 及びδhについても同様に計算して、ＨＳＰ値(δt)＝16.22 と算出することができる。
「Ｏ−０２」〜「Ｏ−１１」分子の場合も同様にして算出する。分子を形成している基に関し、Krevelen & Hoftyzerの文献に記載されている基については、Fd値, Fp値, Eh値及びVc値の記載値を用いればよい。記載されていない基については、構造的に近似する他の基の情報を用いて推定した値を用いるなどを行えばよい。
ii．融点の求め方
融点の推算は、原子団寄与法の一つであり、「Joback, K. G., Reid, R. C., Chem. Eng. Comm., 57, 233 (1987).」に記載されているJoback法を用いればよい。
次に、多成分系溶液を構成する各成分のうちの所望の温度未満の融点を有する成分を液相成分として分類し、前記所望の温度以上の融点を有する成分を非液相成分として分類する（Ｓ２）。
例えば、当該モデル重質油の場合、その液温が２５０℃である場合、この温度よりも低い融点を有する「Ｏ−０１」〜「Ｏ−０３」、「Ｏ−０９」及び「Ｏ−１０」分子が液相成分として分類され、一方、この温度よりも高い融点を有する「Ｏ−０４」〜「Ｏ−０８」及び「Ｏ−１１」分子が非液相成分として分類される。
次に、液相成分として分類された各成分のＨＳＰ値について、各成分の当該液相における容量分率で重み付けした加重平均値を、液相全体の平均ＨＳＰ値として算出する（Ｓ３）。各成分について、密度、分子量等の物性に関する諸情報を予め記憶部２に格納しておくことにより、容量分率は成分情報算出部１１において、容易に算出することができる。
次に、液相全体の前記平均ＨＳＰ値と、非液相成分における各成分のＨＳＰ値との差を算出する（Ｓ４）。
次に、非液相成分における各成分を、差（Δδ）に基づいて、液相成分又は非液相成分として再分類し、液相成分として再分類された各成分を非液相成分から液相成分に編入して、液相成分及び非液相成分を更新する（Ｓ５）。
なお、この再分類における更新は、非液相成分における各成分について、１つずつ順番に行ってもよし、複数の成分ごとに行ってもよい。
次に、更新後の液相成分における各成分のＨＳＰ値について、各成分の当該更新後の液相における容量分率で重み付けした加重平均値を、更新後の液相全体の平均ＨＳＰ値として算出する（Ｓ６）。
そして、ステップＳ４〜Ｓ５を、ステップＳ５において液相成分として再分類される非液相成分がなくなる最終段階まで繰り返す（Ｓ６）。
このようにして、最終段階での更新後の液相成分をモデル重質油のその温度における液相の成分とし、また、最終段階での更新後の平均ＨＳＰ値をモデル重質油のその温度での液相全体の平均ＨＳＰ値として決定する。また、最終段階での更新後の液相成分における分類された各成分の分率の合計を液相分率として算出する。
続いて、所望の温度における最終段階での更新後の非液相成分の凝集度Ｄを算出する（Ｓ７）。
次に、最終段階での更新後の非液相成分における各成分を、凝集度Ｄに基づいて、凝集相成分と固相成分とに分類する（Ｓ８）
さらに、これらの結果に基づいて、モデル重質油の種々の性状を表すパラメータを算出する。例えば、凝集相成分に分類された各成分の分率の合計を凝集相分率として算出し、また、固相成分に分類された各成分の分率の合計を固相分率として算出する。さらに、凝集相成分に分類された各成分の凝集度の和を当該成分の数で除した値を凝集相全体の平均凝集度として算出する。
上記の方法により、このモデル重質油について、先ず無溶媒での液相、凝集相及び固相の量及び組成並びに凝集相における各分子の凝集度 及び凝集相の平均凝集度を算出、推定した。実験は-50℃から350℃まで100℃間隔で実施した。
表２に、各温度において、各分子がいかなる相に存在するかを示す。これはまた、各相の分子組成を示しているものである。
さらに、図４に、各温度における各相の量と凝集相の平均凝集度を示す。図４中、「Liquid」は「液相」、「Aggregate」は「凝集相」、「Solid」は「固相」の意味である。また、「Phase wt ratio」は、「各相の重量分率」の意味である。さらに、「Averaged Dagg」は、「平均凝集度」の意味である。「Temperature」は「温度」の意味である。本明細書における他の図においても同様である。
表２及び図４に示すように、無溶媒の場合、150℃までは全て固相であり、250℃以上では全て凝集相又は液相となり、350℃では液相の割合が増加する。凝集相の平均凝集度は250℃と300℃でそれぞれ1.53と1.25である。
（４．処理方法）
次に、本発明の多成分系溶液の処理方法の実施形態を説明する。
処理にあっては、上記のように、原多成分系溶液としてモデル重質油の性状を推定し、さらに、モデル重質油に、溶媒を混合したり、温度を変更したりしたときに新モデル重質油の性状を、上記の方法により予測する。そして、予測と同じ条件で、溶媒をその分率でモデル重質油に混合したり、モデル重質油の温度の予測時の温度に変更することによって、溶液の処理を行う。
これにより、本発明においては、上記の方法により算出又は推定された溶液中における液相、凝集相及び固相の各々の量及び組成並びに凝集相における各成分の凝集度及び凝集相の平均凝集度に基づいて、凝集相及び固相の全部又は一部を溶解し、又は成分の凝集若しくは析出を抑制するための措置を施すことにより溶液を処理することができる。
（４−１．凝集緩和処理）
まず、処理例として、溶液における凝集相及び固相の溶解、又は成分の凝集若しくは析出抑制のための処理を説明する。溶液中における液相、凝集相及び固相の各々の量及び組成並びに凝集相における各成分の凝集度及び凝集相の平均凝集度が判ったならば、凝集相及び固相の全部又は一部を溶解し、又は成分の凝集若しくは析出を抑制するためには、いかなる措置をとればよいのかという指針を得ることが可能となる。
ここで、モデル重質油を30wt%トルエン溶液とした場合の液相、凝集相及び固相の量及び組成、並びに凝集相の平均凝集度を算出、推定した。下記の表３に、各温度において、各分子がいかなる相に存在するかを示す。これはまた、各相の分子組成を示しているものである。
さらに、図５に、各温度における各相の量と凝集相の平均凝集度を示す。
表３及び図５に示すように、トルエン溶液では、-50℃においても固相は10%に過ぎず、90%は凝集相で溶媒中に分散している。50℃では100%が凝集相となり、150℃以上では昇温と共に凝集相が減少し液相が増加する。凝集相の平均凝集度は-50℃の1.81から350℃の1.27に向かって昇温と共に減少することがわかる。
さらに、モデル重質油のトルエン溶液について、凝集相及び固相を溶解するには、ＨＳＰ値が最も大きいＯ−１１分子を基準にして、Ｏ−１１分子が溶解するよう、添加する溶剤の種類及び量を決めることになる。具体的には、溶媒をトルエンから、トルエンとブロモベンゼン(ＨＳＰ値20.4（δｄ＝21.0、δｐ＝2.5、δｈ＝2.0、δｔ＝20.4）)の等重量比混合溶媒に変更することで、３５０℃の液温下で、すべての凝集相及び固相を溶解させることができる。
さらに、図６に、各温度における各相の量と凝集相の平均凝集度を示す。
このように、モデル重質油における凝集相及び固相の析出量、析出性状等により、
とりうる措置は様々であるが、具体的には、例えば、現状と同じ溶媒を添加する、或いは、新たに別種の溶媒(溶剤)を添加する、溶液の温度を上げる等が考えられる。特に、溶媒(溶剤)の添加については、現状の溶液、追加する溶媒(溶剤)、凝集相及び固相の各々のＨＳＰ値に基づいて溶媒(溶剤)追加後の性状を予測することにより、いかなるＨＳＰ値を有する溶媒(溶剤)を、いかなる量を添加すればよいのかを決めることができる。
（４−２．析出促進処理）
次に、処理例として、溶液における析出促進のための処理を説明する。
溶液中における液相、凝集相及び固相の各々の量及び組成並びに凝集相における各成分の凝集度及び凝集相の平均凝集度が判ったならば、推定された溶液中における液相、凝集相及び固相の各々の量及び組成並びに凝集相における各成分の凝集度及び凝集相の平均凝集度に基づいて、一種以上の成分を析出させるためには、いかなる措置をとればよいのかという指針を得ることが可能となる。
このように、当該溶液における液相の性状等により、とりうる措置は様々であるが、具体的には、例えば、新たに別種の溶媒(溶剤)を添加する、溶液の温度を下げる等が考えられる。特に、溶媒(溶剤)の添加については、現状の溶液、追加する溶媒(溶剤)、凝集相及び固相の各々のＨＳＰ値に基づいて溶媒(溶剤)追加後の性状を予測することにより、いかなるＨＳＰ値を有する溶媒(溶剤)を、いかなる量を添加すればよいのかを決めることができる。
このように、溶液中における液相、凝集相及び固相に関し、それらの量及び組成を算出、推定し、また凝集相における各成分の凝集度及び凝集相の平均凝集度を算出することができるので、例えば、溶媒種を適切に選択する、或いは新たな溶媒を追加する、溶液の温度を変更する等の措置を適切に講ずれば、溶液中における凝集相及び固相の全部又は一部を溶解し、又は成分の凝集若しくは析出を抑制すること、或いは一種以上の成分を析出させることが可能となる。
なお、上記の実施形態では、多成分系溶液として、モデル重質油の性状を推定した例を説明したが、本発明では、多成分系溶液は、重質油に限定されない。
本発明により、コンピュータを用いて、溶液中における液相、凝集相及び固相に関し、それらの量並びに凝集相における各成分の凝集度及び凝集相の平均凝集度を算出し、これらの各相の組成を推定する方法を特にアスファルテンの凝集緩和、抑制に応用することは、石油産業にとって非常に大きな意味を持っている。
すなわち、アスファルテンの凝集・析出を効果的に制御できれば、原油採掘における坑井の閉塞、原油輸送のパイプラインにおけるアスファルテンの析出による閉塞、タンク底部でのスラッジ(析出固体)の滞留、さらには、石油精製過程における多くの課題、例えば、加熱炉パイプの詰まり、熱交換器中へのスラッジ析出による熱効率の低下、重質油分解装置における触媒上にコークが析出することによる触媒性能の劣化等、諸課題に対し、その解決策を与えることができるものであり、また溶剤の添加による脱アスファルテン処理等にも応用が可能である。

Claims (15)

  1. コンピュータによる多成分系溶液の性状推定方法であって、
    （１）前記多成分系溶液を構成する各成分の分率、融点及びハンセン溶解度指数値を取得するステップと、
    （２）前記多成分系溶液を構成する各成分のうち、所望の温度未満の融点を有する成分を液相成分として分類し、前記所望の温度以上の融点を有する成分を非液相成分として分類するステップと、
    （３）前記液相成分として分類された各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、液相全体の平均ハンセン溶解度指数値を算出するステップと、
    （４）液相全体の前記平均ハンセン溶解度指数値と、非液相成分における各成分のハンセン溶解度指数値との差を算出するステップと、
    （５）非液相成分における各成分を、前記差に基づいて、液相成分又は非液相成分として再分類し、液相成分として再分類された各成分を非液相成分から液相成分に編入して、液相成分及び非液相成分を更新するステップと、
    （６）更新後の液相成分における各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、更新後の液相全体の平均ハンセン溶解度指数値を算出するステップと、
    （７）前記ステップ（４）〜（６）を、前記ステップ（５）において液相成分として再分類される非液相成分がなくなる最終段階まで繰り返すステップとを含むことを特徴とする、多成分系溶液の性状推定方法。
  2. 前記最終段階での更新後の液相成分における各成分の分率の合計を液相分率として算出するステップを更に含む
    ことを特徴とする、請求項１記載の多成分系溶液の性状推定方法。
  3. 前記所望の温度における前記最終段階での更新後の非液相成分における各成分の凝集度を、前記最終段階での更新後の液相全体の前記平均ハンセン溶解度指数値と前記最終段階での更新後の非液相成分における各成分のハンセン溶解度指数値との差及び前記最終段階での更新後の非液相成分における各成分の濃度に基づいて算出するステップを更に有する
    ことを特徴とする、請求項１又は２記載の多成分系溶液の性状推定方法。
  4. 前記非液相成分における各成分の凝集度（Ｄ）を下記の式（Ａ）により算出する
    ことを特徴とする、請求項３記載の多成分系溶液の性状推定方法。
    Ｄ(ｐ，ｑ) = ＭＡＳ (Ｋ０ +Ｋ１ p +Ｋ２ ｑ +Ｋ３ ｐ２ +Ｋ４ ｐｑ +Ｋ５ ｑ２+Ｋ６ｐ３ +Ｋ７ｐ２ｑ +Ｋ８ｐｑ２ +Ｋ９ ｑ３) ・・・（Ａ）
    ここで、ｐは、前記所望の温度Ｔが、Ｔ≦１５０℃のときに、ｐ = (Ｌ０(Ｔ - 25) + Ｌ１)ＲＥＤｇ、前記所望の温度Ｔが、１５０℃＜Ｔ≦２００℃のときに、ｐ = (Ｌ０(150 - 25) + Ｌ１)ＲＥＤｇ、前記所望の温度Ｔが、２００℃＜Ｔのときに、ｐ = (Ｌ０(Ｔ - 25) + Ｌ２)ＲＥＤｇで表され、
    Ｌ０、Ｌ１及びＬ２は、係数であり、
    ＲＥＤｇは、ＲＥＤ≧０．３のときに、ＲＥＤｇ＝ＲＥＤ、ＲＥＤ＜０．３のときに、ＲＥＤｇ＝０．３と表され、
    ＲＥＤは、ＲＥＤ＝Δδ／Ｒ０で表され、
    Δδは、液相全体の前記平均ハンセン溶解度指数値と前記非液相成分における各成分のハンセン溶解度指数値との差であり、
    Ｒ０は、非液相成分における各成分の定数であり、
    ｑは、ｑ＝ｌｏｇＣで表され、
    Ｃは、非液相成分における各成分の濃度であり、
    ＭＡＳ、及びＫ０〜Ｋ９は、係数である。
  5. 前記最終段階での更新後の非液相成分における各成分のうち、前記凝集度が所定の閾値未満の成分を凝集相成分として分類し、前記凝集度が前記所定の閾値以上の成分を固相成分として分類するステップを更に有する
    ことを特徴とする、請求項３又は４記載の多成分系溶液の性状推定方法。
  6. 前記凝集相成分として分類された各成分の分率の合計を凝集相分率として算出するステップと、
    前記固相成分として分類された各成分の分率の合計を固相分率として算出するステップとを更に有する
    ことを特徴とする、請求項５記載の多成分系溶液の性状推定方法。
  7. 前記凝集相成分として分類された各成分の凝集度の和を当該成分の数で除した値として、凝集相全体の平均凝集度を算出するステップを更に有する
    ことを特徴とする、請求項５又は６記載の多成分系溶液の性状推定方法。
  8. 前記多成分系溶液には、所望の種類の溶媒が所望の分率で混合されている
    ことを特徴とする、請求項１〜７のいずれか一項に記載の多成分系溶液の性状推定方法。
  9. 前記多成分系溶液は、重質油、又は重質油と溶媒との混合液である
    ことを特徴とする、請求項１〜８のいずれか一項に記載の多成分系溶液の性状推定方法。
  10. 多成分系溶液の性状推定装置であって、
    前記多成分系溶液を構成する各成分の分率、融点及びハンセン溶解度指数値を取得する成分情報取得部と、
    前記多成分系溶液を構成する各成分のうちの所望の温度未満の融点を有する成分を液相成分として分類し、前記所望の温度以上の融点を有する成分を非液相成分として分類する初期分類部と、
    液相の性状を推定する液相演算部と
    を備え、
    前記液相演算部は、
    前記液相成分として分類された各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、液相全体の平均ハンセン溶解度指数値を算出し、
    液相全体の前記平均ハンセン溶解度指数値と、前記非液相成分における各成分のハンセン溶解度指数値との差を算出し、
    非液相成分における各成分を、前記差に基づいて、液相成分又は非液相成分として再分類し、液相成分として再分類された各成分を非液相成分から液相成分に編入して、液相成分及び非液相成分を更新し、
    更新後の液相成分における各成分のハンセン溶解度指数値を各成分の当該液相における分率で重み付けした加重平均値として、更新後の液相全体の平均ハンセン溶解度指数値を算出し、
    液相成分として再分類される非液相成分がなくなる最終段階まで、平均ハンセン溶解度指数値、液相成分、及び非液相成分の更新を繰り返すことを特徴とする多成分系溶液の性状推定装置。
  11. 非液相の性状を推定する非液相演算部を更に備え、
    前記非液相演算部は、前記所望の温度における前記最終段階での更新後の非液相成分における各成分の凝集度を、液相全体の前記平均ハンセン溶解度指数値と前記非液相成分における各成分のハンセン溶解度指数値との差と、前記最終段階での更新後の非液相成分における各成分の濃度とに基づいて算出することを特徴とする請求項１０記載の多成分系溶液の性状推定装置。
  12. コンピュータに、請求項１〜９のいずれか一項に記載の多成分系溶液の性状推定方法を実行させるためのプログラム。
  13. 請求項１〜８記載の多成分系溶液の性状推定方法により、原多成分系溶液の性状を推定し、さらに当該性状推定結果に基づいて、前記原多成分系溶液に所望の種類の溶媒を所望の分率で混合した場合、前記原多成分系溶液の温度を所望の温度に変更した場合、又は前記原多成分系溶液に所望の種類の溶媒を所望の分率で混合しかつ前記原多成分系溶液の温度を所望の温度に変更した場合における新多成分系溶液の性状を予測し、その予測を基に、前記所望の種類の溶媒を前記所望の分率で前記原多成分系溶液に混合する、又は前記原多成分系溶液の温度を前記所望の温度に変更する、又は前記所望の種類の溶媒を前記所望の分率で前記原多成分系溶液に混合しかつ前記原多成分系溶液の温度を前記所望の温度に変更する
    ことを特徴とする、多成分系溶液の処理方法。
  14. 前記所望の種類の溶媒及び前記所望の分率、及び／又は前記所望の温度は、予測された前記新多成分系溶液の性状が、推定された前記原多成分系溶液における凝集相成分及び固相成分の全部又は一部が溶解し、又は、凝集相成分及び固相成分の凝集又は析出が抑制された性状となる溶媒及び分率、及び／又は温度である
    ことを特徴とする、請求項１３記載の多成分系溶液の処理方法。
  15. 前記所望の種類の溶媒及び前記所望の分率、及び／又は前記所望の温度は、予測された前記新多成分系溶液の性状が、推定された前記原多成分系溶液における非固相成分のうちの一種類以上の成分が析出した性状となる溶媒及び分率、及び／又は温度である
    ことを特徴とする、請求項１３記載の多成分系溶液の処理方法。
JP2016150586A 2016-07-29 2016-07-29 多成分系溶液の性状推定方法、装置、及びプログラム Active JP6725112B6 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016150586A JP6725112B6 (ja) 2016-07-29 2016-07-29 多成分系溶液の性状推定方法、装置、及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016150586A JP6725112B6 (ja) 2016-07-29 2016-07-29 多成分系溶液の性状推定方法、装置、及びプログラム

Publications (3)

Publication Number Publication Date
JP2018018448A JP2018018448A (ja) 2018-02-01
JP6725112B2 JP6725112B2 (ja) 2020-07-15
JP6725112B6 true JP6725112B6 (ja) 2020-08-19

Family

ID=61076236

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016150586A Active JP6725112B6 (ja) 2016-07-29 2016-07-29 多成分系溶液の性状推定方法、装置、及びプログラム

Country Status (1)

Country Link
JP (1) JP6725112B6 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7320971B2 (ja) * 2019-03-29 2023-08-04 一般財団法人石油エネルギー技術センター 多成分混合物の物性値推算方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4572777A (en) * 1982-12-14 1986-02-25 Standard Oil Company (Indiana) Recovery of a carbonaceous liquid with a low fines content
JP4523458B2 (ja) * 2005-03-03 2010-08-11 株式会社神戸製鋼所 石油系重質油の水素化分解方法
TWI374150B (en) * 2005-03-17 2012-10-11 Dow Global Technologies Llc Filled polymer compositions made from interpolymers of ethylene α-olefins and uses thereof
BR112014017618B1 (pt) * 2012-01-18 2021-03-30 Prad Research And Development Limited Método para a caracterização de fluido de petróleo em um reservatório atravessado por, pelo menos, uma boca de poço
JP6395178B2 (ja) * 2013-04-09 2018-09-26 国立研究開発法人産業技術総合研究所 アスファルテンの処理方法及びアスファルテン凝集緩和用材料

Also Published As

Publication number Publication date
JP6725112B2 (ja) 2020-07-15
JP2018018448A (ja) 2018-02-01

Similar Documents

Publication Publication Date Title
CA2640914C (en) A multi-stage process for treating bitumen froth using a paraffinic diluent
Mansoori A unified perspective on the phase behaviour of petroleum fluids
AU2016303647B2 (en) System and method of predictive analytics for dynamic control of a hydrocarbon refining process
Stratiev et al. Investigation of relationships between petroleum properties and their impact on crude oil compatibility
Jamei et al. Estimating the density of hybrid nanofluids for thermal energy application: Application of non-parametric and evolutionary polynomial regression data-intelligent techniques
Duran et al. Kinetics of asphaltene precipitation/aggregation from diluted crude oil
JP6725112B6 (ja) 多成分系溶液の性状推定方法、装置、及びプログラム
Alade et al. Development of models to predict the viscosity of a compressed Nigerian bitumen and rheological property of its emulsions
Vatti et al. Asphaltene aggregation in aqueous solution using different water models: A classical molecular dynamics study
KR101542292B1 (ko) 탄화수소 스트림 내에서 풀라들을 분산하는 방법
Xu et al. Descriptor-based methodology for designing heterogeneous microstructural materials system
Aswath et al. Probing the effect of aliphatic ionic liquids on asphaltene aggregation using classical molecular dynamics simulations
Islam et al. Aggregation thermodynamics of asphaltenes: Prediction of asphaltene precipitation in petroleum fluids with NRTL-SAC
Wang et al. Predicting wax appearance temperature and precipitation profile of normal alkane systems: An explicit co-crystal model
Li et al. Effect of temperature on asphaltene precipitation in crude oils from xinjiang oilfield
Yu et al. Impact of Centrifuge-Separated Fractions from Light Tight Oils on the Propensity and Mechanisms of Fouling, and Improved Crude Oil Fouling Potential Prediction Model
Panuganti et al. Harness AI and machine learning in de-emulsifier chemical selection
Durga Prasada Rao et al. Multi-objective optimization of mechanical properties of aluminium 7075-based hybrid metal matrix composite using genetic algorithm
Del Bianco et al. Tailoring hydrocarbon streams for asphaltene removal
Mansoori Phase behavior in petroleum fluids (A Detailed Descriptive and Illustrative Account)
Alhajri et al. Predicting asphaltene precipitation in CO2 flooding processes using artificial intelligence techniques
Agostinho Assessment of Asphaltene Stability at Different Temperatures Using Salicylic and ISO-Phthalic Acid as Inhibitors
Shadman et al. Developing a novel colloidal model for predicting asphaltene precipitation from crude oil by alkane dilution
Adeyanju et al. Experimental Study of Effects of Water-in-Oil Emulsions Flow on Wax Deposition in Subsea Pipelines
Zeidani et al. Analysis of crude oil electrostatic desalters performance

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20160801

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190606

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200519

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200609

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200618

R150 Certificate of patent or registration of utility model

Ref document number: 6725112

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250