JP3789324B2

JP3789324B2 - 水素化精製のシミュレーション方法およびそれを用いた水素化精製方法

Info

Publication number: JP3789324B2
Application number: JP2001219364A
Authority: JP
Inventors: 秀晶熊谷; 中村　　憲治; 博紀小山
Original assignee: Japan Energy Corp
Current assignee: Eneos Corp
Priority date: 2001-07-19
Filing date: 2001-07-19
Publication date: 2006-06-21
Anticipated expiration: 2021-07-19
Also published as: JP2003027069A

Description

【０００１】
【発明の属する技術分野】
石油留分などの炭化水素を水素化精製するために用いられる触媒の寿命を予測するためのシミュレーション方法およびそれを用いた水素化精製方法に関する。特に、重質油、すなわち原油、タールサンド、シェールオイルあいるは石炭液化油などを常圧蒸留または減圧蒸留することにより得られる各種の重質留分または残さ油などを水素化精製するための触媒の寿命予測方法に関する。
【０００２】
【従来の技術】
常圧蒸留または減圧蒸留の残さ油などの重質油には、バナジウム、ニッケルなどの金属分が多く含まれており、この金属分は脱硫、脱窒素、分解および脱金属などを目的とした水素化精製用触媒を被毒し、活性を低下させる。そこで、こうした原料油を、予め脱金属を目的とした水素化精製用触媒（以下、脱金属触媒ともいう）で適切な脱金属を行った後に、脱硫、脱窒素および分解などを目的とした水素化精製用触媒で目的の反応を行うことにより、長期間触媒の活性を維持させている。
【０００３】
水素化精製用触媒、特に脱硫触媒は、精製処理により活性が低下して、寿命が尽きると使用できなくなる。そのため、通常、重質油の水素化精製用触媒は半年〜４年程度の間隔で取り替えられている。この触媒の取り替えが、水素化精製コストに占める割合が大きく、運転条件や用いる触媒の組み合わせ、配置などに最適化が求められている。
【０００４】
こうした重質油の水素化精製を、運転方法および触媒組合せを含め、計画的に実行するために、触媒反応のシミュレーションが用いられる。しかしながら、シミュレーションの予測精度を高めるためには、個々の触媒の性能、特に寿命を正確に測定する必要がある。
【０００５】
【発明が解決しようとする課題】
従来は、個々の触媒の性能、特に寿命を正確に求める方法は提供されていなかった。ニッケル、バナジウムなどの金属分の触媒上への蓄積により、触媒の活性が低下することを説明するいくつかのモデルが提案されているが、寿命の予測には十分なものではなかった。
【０００６】
例えば、Ｄａｕｔｚｅｎｂｅｒｇ，Ｆ．Ｍ．ら（Ａ．Ｃ．Ｓ．Ｓｙｍｐ．Ｓｅｒ．ｖｏｌ．６５，Ｐ．２５４（１９７８））は、活性劣化を金属堆積による触媒細孔入口の閉塞によるものと仮定したモデルを提案した。ここで、脱金属反応を１次反応とみなし、触媒自体の脱金属活性と堆積した金属が有する脱金属活性は同じであると仮定した。ところが、発明者らの研究結果から、脱金属反応は１次反応ではなく、金属が堆積していない初期触媒細孔内の脱金属活性と堆積金属上の脱金属活性は異なり、活性の低下も金属による細孔入口閉塞だけでなく触媒活性点の被毒も重要であることがわかった。また、Ｄａｕｔｚｅｎｂｅｒｇらのモデルは、マクロポアを有するいわゆるバイモーダル型触媒には適用できない。Ｄａｕｔｚｅｎｂｅｒｇらは、上記の仮定に基づき、リサイクル型反応器を用いた実験で、生成油金属濃度に対する原料油金属濃度の比の２乗と反応時間との関係が直線であると見なし、前者の値（濃度比の２乗）が１となる時の反応時間を外挿により求め、さらにこの値をモデル式に代入して脱金属活性がなくなる金属堆積量である最大金属堆積量を求める方法を提案している。
【０００７】
もう一つの例は、Ｔａｍｍ，Ｐ．Ｗ．ら（Ｉｎｄ．Ｅｎｇ．Ｃｈｅｍ．Ｐｒｏｃ．Ｄｅｖ．，ｖｏｌ．２０，ｐ．２６２（１９８１））が提案している金属の分配因子（ＤｉｓｔｒｉｂｕｔｉｏｎＦａｃｔｏｒ）を用いる方法である。これは、水素化精製実験終了後に回収した触媒粒子について、エックス線マイクロアナライザ（ＸＭＡ）を用いて金属の触媒粒子内部の堆積分布測定を行い、触媒粒子内に堆積した平均の金属堆積濃度を触媒粒子外周付近の最大金属濃度で割った値を分散因子と定義している。しかしながら、触媒層内の反応流体の流れが完全に均一ではないこと、充填され触媒粒子表面の流体の流れは均一ではなく、局所的に反応量が異なり、従って金属の堆積の仕方が異なることとなる。このため、この方法は、触媒粒子のサンプリングの仕方によってＸＭＡ分析の結果が異なり、寿命を正確には予測できないという欠点がある。
【０００８】
本発明の目的は、このような課題を解決するものであり、水素化精製触媒の金属堆積による活性劣化を高い精度で、特殊な触媒評価方法を用いることなく、予測することのできるシミュレーション方法およびそれを用いた水素化精製方法を提供するものである。
【０００９】
【課題を解決するための手段】
本発明者らは、これまでの研究から、重質油の処理における水素化精製触媒の活性劣化の主な原因は、最初の数週間における大きな活性の低下を含め、堆積した金属による触媒活性点の被毒と触媒細孔入口の閉塞によるものであることを発見し、これらの触媒活性の変化を金属堆積量の関数として記述できることを発見した。
【００１０】
また、本発明者らは、主な堆積金属成分であるニッケル、バナジウム化合物などが多数の種類の分子からなるにも関わらず、これらの化合物をそれぞれ２つの疑似成分に分割することで、どのような重質油においても、一般に用いられる反応温度において、それぞれの反応速度が１次反応として表現できることも発見した。
【００１１】
これらの発見から、本発明者らは、積分型反応器を用いた水素化精製実験により触媒を評価し、金属堆積量を変数とする関数から触媒活性の変化を求める方法を見出し、本発明に想到した。
【００１２】
本発明による水素化精製のシミュレーション方法は、水素の存在下で触媒と、精製対象金属成分を含有する炭化水素油とを接触させる水素化精製をシミュレーションする方法において、
まず、触媒と炭化水素油を接触させて水素化精製を行い触媒を評価する際に、
（ａ）炭化水素油中の精製対象成分の濃度変化から求まる各精製対象成分の反応速度定数、および、
（ｂ）少なくとも２つの時点における少なくとも２つの異なる反応器内位置での触媒に含まれる精製対象金属成分の濃度を測定し、
その測定結果を、触媒に含まれる精製対象金属成分の濃度に基づき各精製対象成分の反応速度が変化する数式モデルにあてはめて劣化関数パラメータを抽出し、
次に、その劣化関数パラメータを用いた数式モデルを用いて、水素化精製時の各精製対象成分の反応速度および触媒に含まれる精製対象金属成分の経時変化を算出するものである。
【００１３】
このシミュレーション方法において、各精製対象成分の反応速度が、ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶの劣化関数ａ_ｍを含む１次反応（ここで、ΔＣは微少空間での精製対象成分の濃度変化、ｋ_０は触媒に含まれる精製対象金属成分の濃度が０のときの反応速度定数、Ｃは精製対象成分の濃度、ＬＨＳＶは微少空間での液空間速度である。）として表され、その劣化関数ａ_ｍが触媒に含まれる精製対象金属成分の濃度ｍおよびその最大金属堆積濃度ｍ_０を引数とする関数ｆ（ｍ／ｍ_０）で表される数式モデルを用いることが好ましい。特に、この劣化関数ａ_ｍが、
ａ_ｍ＝（ｚ／ｈ）×（１／ｈ）ｔａｎｈ（２ｈ_ｚ（１−ｚ）＋１）／（（１／ｈ）＋ｔａｎｈ（２ｈ_ｚ（１−ｚ）））
で定義さる劣化関数（ここで、ｚ＝（１−ｍ／ｍ_０）^０．５であり、ｍは触媒に含まれる精製対象金属成分の濃度、ｍ_０は最大金属堆積濃度、ｈおよびｈ_ｚは劣化関数パラメータである。）として表される数式モデルを用い、かつ、劣化関数パラメータとして、最大金属堆積濃度ｍ_０、劣化関数パラメータｈおよびｈ_ｚを用いることが好ましい。なお、本明細書において、最大金属堆積量は新触媒重量あたりの金属堆積量で示し、最大金属堆積濃度ｍ_０に等しい。さらには、各精製対象成分として硫黄含有化合物、バナジウム含有化合物およびニッケル含有化合物を含むことが、特には、硫黄含有化合物、バナジウム含有化合物およびニッケル含有化合物を難反応性化合物に対する反応速度式と易反応性化合物に対する反応速度式に分けた数式モデルを用いることが好ましい。
【００１４】
また、本発明の水素化精製のシミュレーション方法を用いて、複数の触媒を評価してその反応速度定数および劣化関数パラメータを抽出し、その複数の触媒を組み合わせた水素化精製をシミュレーションすることが好ましい。
【００１５】
すでに評価されている触媒を用いる場合には、炭化水素油中の精製対象成分の濃度変化を上述の反応速度式ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶで表し、かつ、上記精製対象成分である硫黄含有化合物、バナジウム含有化合物およびニッケル含有化合物をそれぞれ難反応性化合物と易反応性化合物に分割する反応モデルを用い、触媒に対して抽出された、難反応性化合物および易反応性化合物に対応するそれぞれの反応速度定数ｋ_０並びに最大金属堆積濃度ｍ_０を触媒パラメータとして用い、ある時点の触媒に含まれる精製対象金属成分の濃度ｍを計算し、次に、その結果に基づき次の時点における精製対象成分の濃度変化を計算する。複数の種類の触媒を用い、それぞれの触媒に対応する触媒パラメータを用いることが好ましい。
【００１６】
本発明による水素化精製方法は、これらの水素化精製のシミュレーション方法により、将来における水素化精製時の各精製対象成分の反応速度および触媒に含まれる精製対象金属成分の濃度の経時変化を算出し、その結果に基づき水素化精製の運転条件を決定するものである。
【００１７】
【発明の作用・効果】
本発明の方法を用いることによって、リサイクル型反応器のような特殊な微分型反応器を用いることなく、一般の積分型反応器を用いて抽出することができる劣化関数パラメータを基にシミュレーションすることができる。したがって、触媒の評価が容易であり、短時間の評価によって、長期間の触媒劣化をシミュレーションできるため、より正確な触媒寿命の予測や水素化精製の運転条件の最適化が可能となる。特に、複数の種類の触媒を組み合わせた場合の水素化精製を精度よくシミュレーションすることができる。
【００１８】
【好ましい実施の形態】
〔水素化精製〕原料油を水素の存在下で触媒と接触させて精製油を得る。原料油は、精製対象成分を含む炭化水素油、特には重質油、具体的には５０％留出温度３５０℃以上の炭化水素油が好ましい。反応形式としては、単塔または複数塔の固定床もしくは沸騰床での反応に適用できる。反応条件としては、反応温度：３００〜４５０℃、特には３５０〜４２０℃、液空間速度：０．１〜２．０／時、特には０．１〜１．０／時、水素／油比：２００〜２０００、特には４００〜１５００が好ましい。
【００１９】
〔精製対象成分〕精製対象成分としては、硫黄含有化合物、窒素含有化合物、金属含有化合物があげられる。具体的な精製対象金属成分としては、ニッケル、バナジウムなどの遷移金属があげられ、これらが金属元素重量として１０重量ｐｐｍ以上、特には１００重量ｐｐｍ以上含まれる原料油を精製する場合に適する。また、硫黄化合物の含有量である硫黄分が硫黄重量として１％以上、特には２％以上含まれる原料油を精製する場合に適する。
【００２０】
〔触媒〕代表的には、アルミナ、シリカ−アルミナなどの無機酸化物の多孔性担体に、タングステン、モリブデン、ニッケル、コバルトなどの水素化活性成分を担持した触媒が好ましく用いられる。複数種類を組み合わせて用いる場合にもシミュレーションすることができる。マクロポアを有するバイモーダル型触媒についても本発明の方法は適用できる。
【００２１】
〔触媒の評価〕本発明では、まず目的とする水素化精製とほぼ同等の条件で触媒の評価を行う。具体的には、評価時の水素化精製運転条件において、反応温度は任意に選択できるが、３４０〜４２０℃、特には３６０〜４００℃の範囲で行うことが好ましい。全触媒層の液空間速度は、０．１〜４／時、特には０．３〜２／時が好ましい。反応圧力および水素/油比は、最大金属堆積濃度ｍ_０および劣化関数パラメータに影響を与えるため、実際に予想しようとする反応条件に近い値を設定することが好ましい。
【００２２】
評価時の反応容器は、積分型固定床式等温反応器を用いても、微分型固定床式等温反応器を用いてもよいが、積分型固定床式等温反応器を用いることが好ましい。反応器は、１塔式よりも複数直列で各反応器出口の生成油を採取できる方が、最大金属堆積濃度ｍ_０および劣化関数パラメータを求めやすい。複数種の触媒は、別々に、または、組み合わせて評価することもできる。
【００２３】
劣化関数パラメータを抽出するためには、触媒硫化後、一定の原料で、一定の反応条件で触媒の寿命評価実験を行うことが好ましい。反応終了の条件として、反応終了時に触媒層の最上流部（最上部）における触媒に含まれる精製対象金属成分の濃度ｍが最大金属堆積濃度ｍ_０に近い値であること、特には、ｍ／ｍ_０が０．８〜１．０であることが劣化関数パラメータを正確に求める上で好ましい。実験終了後に、少なくとも触媒を上流から下流（上部から下部）にかけて５つに分けて、回収することが望ましい。回収した各触媒は、溶媒で洗浄後、元素分析により、新触媒重量当たりの金属堆積重量を求める。この評価では触媒が完全に失活するまでの実験を必要としないため、従来よりも短い期間で水素化精製用触媒の活性の変化を金属堆積濃度を変数とする関数として求めることができる。
【００２４】
〔反応速度定数の抽出〕精製対象成分として硫黄含有化合物およびバナジウム・ニッケル含有化合物が含まれ、さらには、窒素含有化合物などが含まれる。これらの各精製対象成分を難反応性化合物に対する反応速度定数と易反応性化合物に対する反応速度定数の２つにそれぞれ分けることで、シミュレーションの精度を向上できる。液空間速度ＬＨＳＶ及び反応温度は、それぞれ複数の条件で、好ましくは３以上の条件で評価を行うことが好ましい。
【００２５】
精製対象成分の除去は、１次反応式：ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶ（ここで、ΔＣは微少空間での精製対象成分の濃度変化、Ｃは精製対象成分の濃度、ＬＨＳＶは微少空間での液空間速度である。）により、１次反応の反応速度定数が定義され、劣化関数ａ_ｍは、反応速度定数の抽出時においてａ_ｍ＝１として評価する。
【００２６】
実運転に近い反応温度において、少なくとも３水準の液空間速度について精製対象成分、例えばバナジウム、ニッケルなどの転換率を測定する。精製対象成分について、２つの疑似成分の比率と２つの疑似成分の１次反応速度定数を最小自乗法により求める。この２つの疑似成分は、上述の難反応性化合物と易反応性化合物に対応する。さらに、これらの値を少なくとも３水準の反応温度について求め、アレニウスの式により、各成分の１次反応の活性エネルギーを最小自乗法により求めることができる。２つの疑似成分の比率が、反応温度により多少異なる場合は、それらの平均値を用いる。
【００２７】
〔劣化関数パラメータの抽出〕精製対象成分の除去は、１次反応式：ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶで表される。（ここで、ΔＣは微少空間での精製対象成分の濃度変化、ｋ_０は触媒に含まれる精製対象金属成分の濃度が０のときの反応速度定数、Ｃは精製対象成分の濃度、ＬＨＳＶは微少空間での液空間速度である。）劣化関数ａ_ｍは、触媒に含まれる精製対象金属成分の濃度ｍにより変化し、ｍ＝０のときａ_ｍ＝１、ｍ＝ｍ_０のときａ_ｍ＝０となる関数が劣化関数として選ばれる。触媒細孔入口が金属の堆積により閉塞して失活することを表現するパラメータである最大金属堆積濃度ｍ_０と、金属の堆積により活性点を被毒する効果および／または反応物の細孔内拡散抵抗の増大を表す少なくとも１つのパラメータが劣化関数に含まれることが好ましい。
【００２８】
通常、精製対象金属成分としては、バナジウムのみ、または、バナジウムおよびニッケルを対象とする。特に、劣化関数として、
ａ_ｍ＝（ｚ／ｈ）×（（１／ｈ）ｔａｎｈ（２ｈ_ｚ（１−ｚ）＋１）／（（１／ｈ）＋ｔａｎｈ（２ｈ_ｚ（１−ｚ）））
で定義される劣化関数（ここで、ｚ＝（１−ｍ／ｍ_０）^０．５であり、ｍは触媒に含まれる精製対象金属成分の濃度、ｍ_０は最大金属堆積量、ｈおよびｈ_ｚは劣化関数パラメータである、）を用いることが好ましい。劣化関数パラメータｈおよびｈ_ｚは、金属の堆積により活性点を被毒する効果および／または反応物の細孔内拡散抵抗の増大を表している。
【００２９】
水素化精製による他の主反応、例えば脱硫反応、脱窒素反応、分解反応についても、それぞれ精製成分、または、必要に応じて２成分以上の疑似成分に分割して、任意の次数の差分反応速度式で記述し、同様の劣化関数を乗じて反応速度の変化を表すことができる。通常の水素化精製の脱硫反応、脱窒素反応、脱金属反応は、１次の反応速度式で表すことができる。
【００３０】
このような劣化関数を配慮して、反応容器内の触媒を原料油の流れ方向に複数に、好ましくは５段階以上に分割し、それぞれの区分での金属堆積量を、劣化関数を含んだ反応速度式から算出することができる。さらに、これらの区分を組み合わせて、積分して金属堆積量を求めることで、精製対象成分の転換率の経時変化、および、反応終了後の触媒層内の金属堆積分布を計算することができる。この金属堆積分布が、上述の評価による金属堆積量と一致するように劣化関数の劣化関数パラメータを選ぶことで、任意の触媒の組み合わせ、運転条件での水素化精製をシミュレーションすることができる。同時に、精製対象成分の転換率の経時変化が上述の評価における値と一致するように劣化関数パラメータを選ぶことで、さらにシミュレーションの精度を上げることができる。
【００３１】
〔シミュレーション方法〕上述のシミュレーション方法を利用すれば、複数の触媒を組み合わせた水素化精製において、将来における水素化精製時の各精製対象成分の反応速度定数および触媒に含まれる精製対象金属成分の濃度の経時変化を算出することができるので、その結果に基づき水素化精製の運転条件を最適化することもできる。特に、脱金属触媒の下流に他の水素化精製触媒を充填するような水素化精製プロセスについても、その上流に脱金属触媒を充填することなく、直接原料油で水素化精製実験を行い、触媒活性の変化を金属堆積量を変数とする関数を求めることができる、
【００３２】
【実施例】
〔反応速度定数と活性化エネルギーの評価〕
表１の性状の市販重油脱硫触媒Ａについて、表２の性状のラタウィ重油を用い、表３の条件下で水素化精製を行った。温度の異なる４つの条件および液空間速度の異なる４つの条件で、脱硫反応並びにバナジウム化合物およびニッケル化合物の脱金属反応の反応速度定数を求めた。この際、硫黄、バナジウムおよびニッケルを含有する化合物をそれぞれ２つの疑似成分に分け、１次反応速度式の和で記述した。１次反応速度式は、ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶで定義され、ａ_ｍは劣化関数、ｋ_０は触媒に含まれる精製対象金属成分の濃度が０のときの反応速度定数、ΔＣは微少空間での精製対象成分の濃度変化、Ｃは精製対象成分の濃度、ＬＨＳＶは微少空間での液空間速度である。なお、この評価では劣化関数ａ_ｍ＝１として評価する。
【００３３】
【表１】

【００３４】
【表２】

【００３５】
【表３】

【００３６】
実験結果に合致するように、難反応性と易反応性の成分に硫黄、バナジウムおよびニッケルを含有する化合物を分けて、難反応性成分と易反応性成分の割合および６つの反応速度定数と活性化エネルギーを定めた。この結果を表４に示す。
【００３７】
【表４】

【００３８】
バナジウム化合物の脱金属反応を２つの疑似成分に分けた計算結果と測定結果を図１、図２に示す。図１は、液空間速度によるバナジウム化合物の転換率の変化を示した図であり、破線は易反応性成分を、一点破線は難反応性成分の転換率を表４の反応速度定数により計算した結果をそれぞれ示している。この両者の合計を実線で示しているが、○印で示した実験結果と一致していることがわかる。
【００３９】
図２は、バナジウム化合物のアレニウスプロット（反応温度による反応速度定数の変化）を示した図であり、実線は易反応性成分を、破線は難反応性成分の反応速度定数を表４の活性化エネルギーにより計算した結果をそれぞれ示している。各温度での実験結果から表４の難反応性成分と易反応性成分の割合を用いて計算した反応速度定数を○印と△印で示しているが、これらが計算結果と一致していることがわかる。
【００４０】
〔触媒の劣化関数パラメータの評価〕
表１の性状の市販重油脱硫触媒Ａについて、表２の性状のラタウィ重油を用い、２段の反応器を用いて表５の条件で水素化精製を行った。
【００４１】
【表５】

【００４２】
第１段の反応器の出口と第２段の反応器の出口における硫黄、バナジウムおよびニッケルの含有量の経時変化を測定した。所定の時間反応を行った後、各反応器の触媒をそれぞれ流れ方向に対して５等分に分割して回収した。各使用済み触媒をソックスレー抽出器で洗浄し、ＩＣＰ発光分析法で元素分析を行い、新触媒の単位重量当たりに含まれるバナジウムおよびニッケルの重量を堆積量として算出した。
【００４３】
触媒の劣化を考慮するため、反応速度定数と活性化エネルギーを求める際に用いた１次反応式：
ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶにおいて、劣化関数ａ_ｍは、Ｋｈａｎｇ，Ｓ−Ｊら（Ｉｎｄ．Ｅｎｇ．Ｃｈｅｍ．Ｐｒｏｃ．Ｄｅｖ．，ｖｏｌ．２５，ｐ．４３７（１９８６））が導いた式：
ａ_ｍ＝（ｚ／ｈ）×（１／ｈ）ｔａｎｈ（２ｈ_ｚ（１−ｚ）＋１）／（（１／ｈ）＋ｔａｎｈ（２ｈ_ｚ（１−ｚ）））
を用いた。ここで、ｚ＝（１−ｍ／ｍ_０）^０．５であり、ｍは触媒に含まれるバナジウムおよびニッケルの濃度、ｍ_０は最大金属堆積量、ｈおよびｈ_ｚは劣化関数パラメータであり、ｈ×ｈ_ｚ＝２０とした。
【００４４】
この一次反応式を用い、反応容器を１０に区分し、それぞれの区分での金属堆積量から反応を予測し、この１０区分を積分した結果が、実験結果と合致するように最大金属堆積量ｍ_０および劣化関数パラメータｈを決定した。この結果を表４に併せて示す。
【００４５】
反応器内のバナジウムおよびニッケルの濃度分布を図３に示す。●印および実線は、表４のパラメータを用いて計算した値であり、○印で示した実験結果と一致していることがわかる。また、第１段の反応器の出口と第２段の反応器の出口におけるバナジウム含有量の経時変化を図４に示す。△印は第１段の反応器の出口、○印は第２段の反応器の出口におけるバナジウム濃度の測定値を示しており、表４の値から計算された濃度変化を実線で示しており、両者はよく一致していることがわかる。
【００４６】
〔商業装置における検証〕
重油脱硫の商業装置によりシミュレーション方法を検証した。５種類の市販触媒を上述の方法で評価し、得られた結果を用いて、重油脱硫のシミュレーションを行った。商業生産装置からの精製油の硫黄濃度、バナジウム濃度の経時変化を図５、図６に示す。○印で運転実績を示しており、●印で示したシミュレーション結果とよく一致していることがわかる。
【図面の簡単な説明】
【図１】触媒評価における液空間速度によるバナジウム化合物の転換率の変化を示した図である。
【図２】触媒評価におけるバナジウム化合物のアレニウスプロットを示した図である。
【図３】触媒寿命評価における反応器内のバナジウムおよびニッケルの濃度分布を示した図である。
【図４】触媒寿命評価における反応器出口のバナジウム濃度の経時変化を示した図である。
【図５】商業生産装置における精製油の硫黄濃度の経時変化をシミュレーション結果とともに示した図である。
【図６】商業生産装置における精製油のバナジウム濃度の経時変化をシミュレーション結果とともに示した図である。

Claims

水素の存在下で触媒と、精製対象金属成分を含有する炭化水素油とを接触させる水素化精製をシミュレーションする方法において、
まず、触媒と炭化水素油を接触させて水素化精製を行い触媒を評価する際に、
（ａ）炭化水素油中の精製対象成分の濃度変化から求まる各精製対象成分の反応速度定数、および、
（ｂ）少なくとも２つの時点における少なくとも２つの異なる反応器内位置での触媒に含まれる精製対象金属成分の濃度を測定し、
その測定結果を、触媒に含まれる精製対象金属成分の濃度に基づき各精製対象成分の反応速度が変化する数式モデルにあてはめて劣化関数パラメータを抽出し、
次に、その劣化関数パラメータを用いた数式モデルを用いて、水素化精製時の各精製対象成分の反応速度および触媒に含まれる精製対象金属成分の経時変化を算出する水素化精製のシミュレーション方法。
請求項１において、
各精製対象成分の反応速度が、ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶの劣化関数ａ_ｍを含む１次反応（ここで、ΔＣは微少空間での精製対象成分の濃度変化、ｋ_０は触媒に含まれる精製対象金属成分の濃度が０のときの反応速度定数、Ｃは精製対象成分の濃度、ＬＨＳＶは微少空間での液空間速度である。）として表され、その劣化関数ａ_ｍが触媒に含まれる精製対象金属成分の濃度ｍおよびその最大金属堆積濃度ｍ_０を引数とする関数ｆ（ｍ／ｍ_０）で表される数式モデルを用いる請求項１記載の水素化精製のシミュレーション方法。
請求項１において、各精製対象成分として硫黄含有化合物、バナジウム含有化合物およびニッケル含有化合物を含む請求項１または２記載の水素化精製のシミュレーション方法。
請求項３において、硫黄含有化合物、バナジウム含有化合物およびニッケル含有化合物を難反応性化合物に対する反応速度式と易反応性化合物に対する反応速度式に分けた数式モデルを用いる請求項３記載の水素化精製のシミュレーション方法。
水素の存在下で触媒と、精製対象金属成分を含有する炭化水素油とを接触させる水素化精製の経時変化をシミュレーションする方法において、
炭化水素油中の精製対象成分の濃度変化を反応速度式ΔＣ＝ａ_ｍ×ｋ_０×Ｃ／ＬＨＳＶ（ここで、ΔＣは微少空間での精製対象成分の濃度変化、ａ_ｍは劣化関数であり、ｋ_０は触媒に含まれる精製対象金属成分の濃度が０のときの反応速度定数、Ｃは精製対象成分の濃度、ＬＨＳＶは微少空間での液空間速度である。劣化関数ａ_ｍは触媒に含まれる精製対象金属成分の濃度ｍおよびその最大金属堆積濃度ｍ_０を引数とする関数ｆ（ｍ／ｍ_０）で表される。）で表し、かつ、上記精製対象成分である硫黄含有化合物、バナジウム含有化合物およびニッケル含有化合物をそれぞれ難反応性化合物と易反応性化合物に分割する反応モデルを用い、
触媒に対して抽出された、難反応性化合物および易反応性化合物に対応するそれぞれの反応速度定数ｋ_０並びに最大金属堆積濃度ｍ_０を触媒パラメータとして用い、
ある時点の触媒に含まれる精製対象金属成分の濃度ｍを計算し、次に、
その結果に基づき次の時点における精製対象成分の濃度変化を計算する水素化精製のシミュレーション方法。
請求項５において、複数の種類の触媒を用い、それぞれの触媒に対応する触媒パラメータを用いる請求項５記載の水素化精製のシミュレーション方法。
請求項１〜６記載の水素化精製のシミュレーション方法により、将来における水素化精製時の各精製対象成分の反応速度および触媒に含まれる精製対象金属成分の濃度の経時変化を算出し、その結果に基づき水素化精製の運転条件を決定する水素化精製方法。