特開2023-142721情報処理方法、反応温度算出装置、反応温度算出プログラム、及びコンピュータの非一時的可読記録媒体
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023142721
(43)【公開日】2023-10-05
(54)【発明の名称】情報処理方法、反応温度算出装置、反応温度算出プログラム、及びコンピュータの非一時的可読記録媒体
(51)【国際特許分類】
C10G 45/02 20060101AFI20230928BHJP
【ＦＩ】
C10G45/02
【審査請求】未請求
【請求項の数】10
【出願形態】ＯＬ
(21)【出願番号】P 2022049754
(22)【出願日】2022-03-25
(71)【出願人】
【識別番号】000105567
【氏名又は名称】コスモ石油株式会社
(74)【代理人】
【識別番号】100141139
【弁理士】
【氏名又は名称】及川 周
(74)【代理人】
【識別番号】100126882
【弁理士】
【氏名又は名称】五十嵐 光永
(74)【代理人】
【識別番号】100162868
【弁理士】
【氏名又は名称】伊藤 英輔
(74)【代理人】
【識別番号】100147267
【弁理士】
【氏名又は名称】大槻 真紀子
(74)【代理人】
【識別番号】100209347
【弁理士】
【氏名又は名称】内田 洋平
(72)【発明者】
【氏名】辻 浩二
(72)【発明者】
【氏名】薄井 大輔
(72)【発明者】
【氏名】関根 翼
(72)【発明者】
【氏名】廣瀬 敏之
(72)【発明者】
【氏名】出井 一夫
【テーマコード（参考）】
4H129
【Ｆターム（参考）】
4H129AA02
4H129CA07
4H129DA15
4H129KA06
4H129KB03
4H129LA14
4H129NA01
4H129NA45
(57)【要約】
【課題】重質コーカー分解軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる情報処理方法、及び前記反応温度を推定することができる反応温度算出装置、コンピュータを、前記反応温度算出装置として機能させるための反応温度算出プログラム、及び前記プログラムを記憶したコンピュータの非一時的可読記録媒体の提供。
【解決手段】反応を開始してから所定時間経過した際の重質コーカー分解軽油を含む原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、取得した前記情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
【選択図】なし
【解決手段】反応を開始してから所定時間経過した際の重質コーカー分解軽油を含む原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、取得した前記情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
【選択図】なし
【特許請求の範囲】
【請求項1】
重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、
取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、
前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、
前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
【請求項2】
前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成される、請求項１に記載の情報処理方法。
【請求項3】
前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、請求項１又は２に記載の情報処理方法。
【請求項4】
前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、請求項１～３のいずれか一項に記載の情報処理方法。
【請求項5】
重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む、反応温度算出装置。
【請求項6】
前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成される、請求項５に記載の反応温度算出装置。
【請求項7】
前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、請求項５又は６に記載の反応温度算出装置。
【請求項8】
前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、請求項５～７のいずれか一項に記載の反応温度算出装置。
【請求項9】
コンピュータを、請求項５～８のいずれか一項に記載の反応温度算出装置として機能させるための反応温度算出プログラム。
【請求項10】
請求項９に記載のプログラムを記憶したコンピュータの非一時的可読記録媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、情報処理方法、反応温度算出装置、反応温度算出プログラム、及びコンピュータの非一時的可読記録媒体に関する。
【背景技術】
【0002】
分解油水添脱硫装置（Coker Distillate HDS Unit, ＤＨＵ）のHigh Coker Gas Oil Section（ＨＳ）は、コーカー分解油の重油留分である重質コーカー分解軽油を主な原料油とする水素化装置であり、流動接触分解装置の原料油となる脱硫重質コーカー分解軽油の硫黄分を低減することを目的として運転している。
【0003】
重質コーカー分解軽油の水素化処理においては、副生成物としてコークが生成し、コークが水素化処理触媒に堆積することにより、水素化処理触媒の活性は経時的に低下する。
また、重質コーカー分解軽油には、ＤＨＵ－ＨＳ装置の前段にあるコーカーに投入されている泡消剤由来のシリコンが含まれており、これらの化合物が水素化処理触媒に堆積することによっても、水素化処理触媒の活性は経時的に低下する。したがって、コーク及びシリコン堆積に伴う水素化処理触媒の活性低下に対して、生成油中の硫黄含有量を一定レベル以下に保つため、反応温度を上げて運転を行わなければならない。
また、重質コーカー分解軽油には、ＤＨＵ－ＨＳ装置の前段にあるコーカーに投入されている泡消剤由来のシリコンが含まれており、これらの化合物が水素化処理触媒に堆積することによっても、水素化処理触媒の活性は経時的に低下する。したがって、コーク及びシリコン堆積に伴う水素化処理触媒の活性低下に対して、生成油中の硫黄含有量を一定レベル以下に保つため、反応温度を上げて運転を行わなければならない。
【0004】
この反応温度の設定において、反応温度が高すぎると、触媒の活性低下が進行するため、所定の運転時間（日数）を達成できなくなり生産性が低下する。一方、反応温度が低すぎると、触媒の活性低下が緩和されるため所定の運転時間までに余力が発生し、余力を活用できないと生産性が低下する。したがって、所定の反応条件（温度、処理量）を精度よく推定するための方法が望まれている。
【0005】
原油精製分野において、最適な反応温度を精度よく推定するための方法については種々の検討が行われている。例えば、特許文献１には、常圧蒸留軽油を含む原料油の水素化処理において、原料油として、性状の異なる常圧蒸留軽油を含む原料油に切り替える場合に、切り替え前の常圧蒸留軽油の性状に関する情報及び運転条件と、切り替え後の常圧蒸留軽油の性状に関する情報及び反応温度以外の運転条件から、切り替え後に必要となる反応温度を推定する方法が開示されている。
【先行技術文献】
【特許文献】
【0006】
【特許文献1】特開平１０－６０４５５号公報
【発明の概要】
【発明が解決しようとする課題】
【0007】
特許文献１に記載の反応温度を推定する方法は、原料油の性状に着目した方法であり、水素化処理触媒の活性低下という観点では検討が行われていない。したがって、水素化処理触媒の活性低下を考慮した上で、最適な反応温度を推定する方法が望まれている。
【0008】
本発明は、上記事情に鑑みてなされたものであって、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる情報処理方法、及び前記反応温度を推定することができる反応温度算出装置、コンピュータを、前記反応温度算出装置として機能させるための反応温度算出プログラム、及び前記プログラムを記憶したコンピュータの非一時的可読記録媒体を提供することを課題とする。
【課題を解決するための手段】
【0009】
上記課題を解決するため、本発明は、以下の態様を有する。
［１］ 重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
［２］ 前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成される、［１］に記載の情報処理方法。
［３］ 前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、［１］又は［２］に記載の情報処理方法。
［４］ 前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、［１］～［３］のいずれか一項に記載の情報処理方法。
［５］ 重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む、反応温度算出装置。
［６］ 前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成される、［５］に記載の反応温度算出装置。
［７］ 前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、［５］又は［６］に記載の反応温度算出装置。
［８］ 前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、［５］～［７］のいずれか一項に記載の反応温度算出装置。
［９］ コンピュータを、［５］～［８］のいずれか一項に記載の反応温度算出装置として機能させるための反応温度算出プログラム。
［１０］ ［９］に記載のプログラムを記憶したコンピュータの非一時的可読記録媒体。
［１］ 重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
［２］ 前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成される、［１］に記載の情報処理方法。
［３］ 前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、［１］又は［２］に記載の情報処理方法。
［４］ 前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、［１］～［３］のいずれか一項に記載の情報処理方法。
［５］ 重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む、反応温度算出装置。
［６］ 前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成される、［５］に記載の反応温度算出装置。
［７］ 前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、［５］又は［６］に記載の反応温度算出装置。
［８］ 前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、［５］～［７］のいずれか一項に記載の反応温度算出装置。
［９］ コンピュータを、［５］～［８］のいずれか一項に記載の反応温度算出装置として機能させるための反応温度算出プログラム。
［１０］ ［９］に記載のプログラムを記憶したコンピュータの非一時的可読記録媒体。
【発明の効果】
【0010】
本発明によれば、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる情報処理方法、及び前記反応温度を推定することができる反応温度算出装置、コンピュータを、前記反応温度算出装置として機能させるための反応温度算出プログラム、及び前記プログラムを記憶したコンピュータの非一時的可読記録媒体を提供することができる。
【図面の簡単な説明】
【0011】
【図1】一実施形態に係る情報処理方法のフローチャートである。
【図2】一実施形態に係る情報処理方法のフローチャートである。
【図3】一実施形態に係る情報処理方法のフローチャートである。
【図4】一実施形態に係る反応温度算出装置の構成ブロック図である。
【発明を実施するための形態】
【0012】
以下、本発明の実施の形態について詳細に説明するが、以下の記載は本発明の実施態様の一例であり、本発明はこれらの内容に限定されず、その要旨の範囲内で変形して実施することができる。
【0013】
≪情報処理方法≫
本実施形態の情報処理方法は、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップ（図１のＳ１）と、取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップ（図１のＳ２）と、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップ（図１のＳ３）と、を含む。以下、各ステップについて説明を行う。なお、以下に示す各ステップは、例えば、本実施形態の反応温度算出装置１によって実行される。例えば、Ｓ１は取得部１１によって実行され、Ｓ２及びＳ３は計算機本体１２における演算部１３によって実行される。
本実施形態の情報処理方法は、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップ（図１のＳ１）と、取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップ（図１のＳ２）と、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップ（図１のＳ３）と、を含む。以下、各ステップについて説明を行う。なお、以下に示す各ステップは、例えば、本実施形態の反応温度算出装置１によって実行される。例えば、Ｓ１は取得部１１によって実行され、Ｓ２及びＳ３は計算機本体１２における演算部１３によって実行される。
【0014】
≪情報取得ステップ≫
本実施形態の情報取得ステップは、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得するステップである。反応を開始してから所定時間経過した際とは、例えば反応を開始してから任意のｔ日経過時である。ｔは整数でも小数でもよく、例えばｔが０．５の場合、反応を開始してから１２時間経過時を意味する。また、ｔ日経過時は、本実施形態の情報処理方法を実施する時からみて過去でも、現在でも、未来でもよい。例えば、本実施形態の情報処理方法を実施する時が、反応を開始してから２日経過時であり、ｔが４の場合、２日後（未来）における反応温度を推定することとなる。
本実施形態の情報取得ステップは、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得するステップである。反応を開始してから所定時間経過した際とは、例えば反応を開始してから任意のｔ日経過時である。ｔは整数でも小数でもよく、例えばｔが０．５の場合、反応を開始してから１２時間経過時を意味する。また、ｔ日経過時は、本実施形態の情報処理方法を実施する時からみて過去でも、現在でも、未来でもよい。例えば、本実施形態の情報処理方法を実施する時が、反応を開始してから２日経過時であり、ｔが４の場合、２日後（未来）における反応温度を推定することとなる。
【0015】
（原料油に関する情報）
原料油に関する情報としては、原料油の組成に関する情報が例として挙げられる。原料油の組成に関する情報としては、原料油中の硫黄濃度及びシリコン濃度に関する情報が例として挙げられる。
原料油に関する情報としては、原料油の組成に関する情報が例として挙げられる。原料油の組成に関する情報としては、原料油中の硫黄濃度及びシリコン濃度に関する情報が例として挙げられる。
【0016】
原料油中の硫黄濃度に関する情報は、本分野で公知の硫黄濃度の測定方法により得ることができ、例えば、紫外蛍光法、波長分散蛍光Ｘ線法等により求めることができる。また、原料油中の硫黄濃度は、原料油を変更することにより制御することができる。原料油中の硫黄濃度に関する情報は設定値であることが好ましい。すなわち、使用を予定している原料油の硫黄濃度を使用することができる。
【0017】
原料油中のシリコン濃度に関する情報は、本分野で公知のシリコン濃度の測定方法により得ることができ、例えば、ＩＣＰ－ＭＳ法、蛍光Ｘ線分析法等により求めることができる。また、原料油中のシリコン濃度は、原料油を変更することにより制御することができる。原料油中のシリコン濃度に関する情報は設定値であることが好ましい。すなわち、使用を予定している原料油のシリコン濃度を使用することができる。
【0018】
（生成油に関する情報）
生成油に関する情報としては、生成油の組成に関する情報が例として挙げられる。生成油の組成に関する情報としては、生成油中の硫黄濃度及びシリコン濃度に関する情報が例として挙げられる。
生成油に関する情報としては、生成油の組成に関する情報が例として挙げられる。生成油の組成に関する情報としては、生成油中の硫黄濃度及びシリコン濃度に関する情報が例として挙げられる。
【0019】
生成油中の硫黄濃度に関する情報は、上述の原料油の場合と同様に求めることができる。本実施形態においては、生成油中の硫黄濃度に関する情報は、設定値であることが好ましい。すなわち、目的とする生成油の硫黄濃度を使用することができる。
【0020】
生成油中のシリコン濃度に関する情報は、上述の原料油の場合と同様に求めることができる。本実施形態においては、生成油中のシリコン濃度に関する情報は、設定値であることが好ましい。すなわち、目的とする生成油のシリコン濃度を使用することができる。なお、上述の生成油中の硫黄濃度を一定とした場合、生成油中のシリコン濃度も一定となる。
【0021】
（運転条件に関する情報）
運転条件に関する情報としては、水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報が例として挙げられる。例えば、運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報の少なくとも何れかである。中でも、運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含むことが好ましく、水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報の全てを含むことがより好ましい。また、運転条件には、反応を開始してから任意のｔ日経過時等の時間の情報も含まれることが好ましい。また、運転条件には、反応温度の実測値に関する情報が含まれることが好ましい。
運転条件に関する情報としては、水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報が例として挙げられる。例えば、運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報の少なくとも何れかである。中でも、運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含むことが好ましく、水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報の全てを含むことがより好ましい。また、運転条件には、反応を開始してから任意のｔ日経過時等の時間の情報も含まれることが好ましい。また、運転条件には、反応温度の実測値に関する情報が含まれることが好ましい。
【0022】
水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報は、本分野において公知の方法により求めることができる。水素分圧に関する情報、触媒充填量に関する情報、原料油の供給量に関する情報、及び水素の供給量に関する情報は、重質コーカー分解軽油を含む原料油の水素化処理反応において制御することが可能である。運転条件に関する情報は設定値であることが好ましい。すなわち、予定の運転条件を使用する。
【0023】
≪劣化度算出ステップ≫
本実施形態の劣化度算出ステップは、取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出するステップである。
本実施形態の劣化度算出ステップは、取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出するステップである。
【0024】
＜劣化度＞
劣化度は下式１で表される。
Φ＝ｋｔ／ｋ０ 式１
前記式１中、ｋ０は反応０日経過時（すなわち、反応開始時）の触媒の反応速度定数であり、ｋｔは任意の反応ｔ日経過時の触媒の反応速度定数である。なお、ｋ０、ｋｔは後述の温度ＴＳＯＲにおける反応速度定数である。
劣化度は下式１で表される。
Φ＝ｋｔ／ｋ０ 式１
前記式１中、ｋ０は反応０日経過時（すなわち、反応開始時）の触媒の反応速度定数であり、ｋｔは任意の反応ｔ日経過時の触媒の反応速度定数である。なお、ｋ０、ｋｔは後述の温度ＴＳＯＲにおける反応速度定数である。
【0025】
本実施形態においては、取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出することができる。
【0026】
＜劣化関数＞
劣化関数は、触媒の劣化度を算出する関数である。本実施形態においては、劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成されることが好ましい。このような劣化関数としては、特に限定されないが、例えば下式２で表される劣化関数が例として挙げられる。
劣化関数は、触媒の劣化度を算出する関数である。本実施形態においては、劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数と、シリコンの堆積による触媒の劣化に関するシリコン劣化関数から構成されることが好ましい。このような劣化関数としては、特に限定されないが、例えば下式２で表される劣化関数が例として挙げられる。
【0027】
Φ＝ΦＣΦSi 式２
前記式２中、ΦＣはコーク劣化関数であり、ΦSiはシリコン劣化関数である。
前記式２中、ΦＣはコーク劣化関数であり、ΦSiはシリコン劣化関数である。
【0028】
以下、コーク劣化関数と、シリコン劣化関数について説明を行う。
【0029】
＜コーク劣化関数１＞
コーク劣化関数としては、触媒のコーク劣化に関する劣化度合いを算出可能な関数であれば、特に限定されないが、例えば下式３で表されるコーク劣化関数１が例として挙げられる。
コーク劣化関数としては、触媒のコーク劣化に関する劣化度合いを算出可能な関数であれば、特に限定されないが、例えば下式３で表されるコーク劣化関数１が例として挙げられる。
【0030】
ΦＣ＝ｅｘｐ（－Ｄｔ） 式３
前記式３中、Ｄは触媒の活性種のコークによる劣化係数であり、ｔは反応経過日数（日）である。
前記式３中、Ｄは触媒の活性種のコークによる劣化係数であり、ｔは反応経過日数（日）である。
【0031】
Ｄは下式４により、求めることができる。下式４は、特定のパラメータにより触媒の活性種のコークによる劣化係数を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
【0032】
【数1】
【0033】
本明細書において、「要求温度」とは、所定の反応条件を達成するために必要な反応温度を意味する。すなわち、０日目の要求温度とは、反応開始時において、上記ＳＦ、ＳＰ、ＬＨＳＶ、及びＰの反応条件を達成するために必要な反応温度を意味する。
【0034】
本明細書において、「基準水素分圧」とは、実際の反応条件の標準的な圧力を意味する。後述する水素分圧係数ａを決定する際に用いる反応圧力の平均値として求められる。
【0035】
本明細書において、「基準反応温度」とは、実際に運転される可能性のある標準的な運転条件の振れ幅で得られるＴＳＯＲの平均値を意味する。
【0036】
前記式４において、ＴＳＯＲは触媒の初期活性を表し、αは触媒の劣化速度を表す。すなわち、ＴＳＯＲとαの数値が大きいほど、触媒の劣化が大きいことを意味し、この劣化挙動はＤの値に反映されることになる。
【0037】
前記式４において、（１／ＳＰ
ｎ－１－１／ＳＦ
ｎ－１）ＬＨＳＶで表される項は、脱硫反応速度定数であり、上述した通り、ＳＰ、ＳＦ、ＬＨＳＶを設定値とし、一定の条件で運転する場合、定数となる。前記式４において、（ＰＢ／Ｐ）ａで表される項は、水素分圧依存性を示す項であり、上述した通り、Ｐを設定値とし、一定の条件で運転する場合、定数となる。前記式４においてｅｘｐ［Ｅｃ／Ｒ（１／ＴＢ－１／ＴＳＯＲ）］で表される項は、温度依存性を示す項であり、定数となる。
【0038】
前記式４中、ＳＦ、ＳＰ、ＬＨＳＶ、Ｐは、上述した情報取得ステップで取得した原料油に関する情報、生成油に関する情報、及び運転条件に関する情報に基づき代入される値である。なお、ＬＨＳＶは、原料油の供給量（体積／時間）を触媒充填量（体積）で除すことにより求めることができる。
【0039】
上述した通り、ＳＦ、ＬＨＳＶ、Ｐは制御することが可能なパラメータである。また、ＳＰは目的とする生成油の硫黄濃度である。すなわち、前記式４により、上記ＳＦ、ＳＰ、ＬＨＳＶ、及びＰの反応条件における、触媒の活性種のコークによる劣化係数を算出することができる。重質コーカー分解軽油を含む原料油の水素化処理反応の反応次数であるｎの求め方は後述する。
【0040】
（コークの基本劣化パラメータの求め方）
前記式４中、α、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは定数である。以下、これらのパラメータを総称して「コークの基本劣化パラメータ１」という。コークの基本劣化パラメータ１は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。コークの基本劣化パラメータを求める際の原料油として重質コーカー分解軽油を含む原料油を使用すると、コークの堆積による活性低下に加え、シリコンの堆積による活性低下も反映されることとなる。したがって、コークの基本劣化パラメータを求める際の原料油としては、シリコンを実質的に含まない原料油を使用することが好ましい。シリコンを実質的に含まないとは、原料油中のシリコン濃度が０．３重量ｐｐｍ未満であることを意味する。本願の発明者らが検討を行った所、重質コーカー分解軽油を含む原料油の水素化処理反応における触媒のコークの堆積による劣化は、常圧蒸留残渣を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応（間接脱硫反応）における触媒のコークの堆積による劣化とほぼ同じであることが判明した。すなわち、コークの基本劣化パラメータを求める上では、減圧蒸留軽油を原料油として用いることができる。以下、コーク劣化関数のパラメータを求める際に「想定する実機運転条件」と記す場合、シリコン濃度以外が想定する実機運転条件であり、原料油中のシリコン濃度は０．３ｐｐｍ未満である。
以下、コークの基本劣化パラメータＰＢ、ａ、Ｅｃ、ＴＢの求め方について例を示す。以下で示す例は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（ＰＢ、ａ、Ｅｃ、ＴＢの求め方）である。また、コークの基本劣化パラメータα、ＴＳＯＲの求め方については２例を示すが、本発明はこれらに限定されるものではない。１例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（α及びＴＳＯＲの求め方１）であり、２例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される反応温度プロファイルから求める方法（α及びＴＳＯＲの求め方２）である。
前記式４中、α、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは定数である。以下、これらのパラメータを総称して「コークの基本劣化パラメータ１」という。コークの基本劣化パラメータ１は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。コークの基本劣化パラメータを求める際の原料油として重質コーカー分解軽油を含む原料油を使用すると、コークの堆積による活性低下に加え、シリコンの堆積による活性低下も反映されることとなる。したがって、コークの基本劣化パラメータを求める際の原料油としては、シリコンを実質的に含まない原料油を使用することが好ましい。シリコンを実質的に含まないとは、原料油中のシリコン濃度が０．３重量ｐｐｍ未満であることを意味する。本願の発明者らが検討を行った所、重質コーカー分解軽油を含む原料油の水素化処理反応における触媒のコークの堆積による劣化は、常圧蒸留残渣を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応（間接脱硫反応）における触媒のコークの堆積による劣化とほぼ同じであることが判明した。すなわち、コークの基本劣化パラメータを求める上では、減圧蒸留軽油を原料油として用いることができる。以下、コーク劣化関数のパラメータを求める際に「想定する実機運転条件」と記す場合、シリコン濃度以外が想定する実機運転条件であり、原料油中のシリコン濃度は０．３ｐｐｍ未満である。
以下、コークの基本劣化パラメータＰＢ、ａ、Ｅｃ、ＴＢの求め方について例を示す。以下で示す例は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（ＰＢ、ａ、Ｅｃ、ＴＢの求め方）である。また、コークの基本劣化パラメータα、ＴＳＯＲの求め方については２例を示すが、本発明はこれらに限定されるものではない。１例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（α及びＴＳＯＲの求め方１）であり、２例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される反応温度プロファイルから求める方法（α及びＴＳＯＲの求め方２）である。
【0041】
（ＰＢ、ａ、Ｅｃ、ＴＢの求め方）
本実施形態のコークの基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。この触媒の劣化挙動（劣化度）は、前記式１と同じ考え方となり、下式５で表すことができる。
Φ’＝ｋｔ’／ｋ０’ 式５
前記式５中、ｋ０’は反応０日経過時（すなわち、反応開始時）の触媒の反応速度定数であり、ｋｔ’は任意の反応ｔ日経過時の触媒の反応速度定数である。なお、ｋ０’、ｋｔ’は後述の温度ＴＳＯＲ’における反応速度定数である。
本実施形態のコークの基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。この触媒の劣化挙動（劣化度）は、前記式１と同じ考え方となり、下式５で表すことができる。
Φ’＝ｋｔ’／ｋ０’ 式５
前記式５中、ｋ０’は反応０日経過時（すなわち、反応開始時）の触媒の反応速度定数であり、ｋｔ’は任意の反応ｔ日経過時の触媒の反応速度定数である。なお、ｋ０’、ｋｔ’は後述の温度ＴＳＯＲ’における反応速度定数である。
【0042】
前記式５は、前記式１と同様、反応速度定数に基づいた劣化関数である。反応速度定数は下式６のアレニウスの式で表される。
【0043】
【数2】
【0044】
反応開始時温度ＴＳＯＲ’における反応速度定数ｋ０’は、反応ｔ日経過時にはｋｔ’に低下する。反応ｔ日経過時に反応速度定数ｋ０’に相当する活性を得るためには、反応温度をＴｔ’としなければならないものとすると、前記式５及び前記式６より、下式７が導かれる。なお、本実施形態の反応は、重質コーカー分解軽油を含む原料油の水素化処理反応のため、活性化エネルギーＥを脱硫の活性化エネルギーＥａ（ｋＪ／ｍｏｌ）としている。
【0045】
【数3】
【0046】
（Ｅｃ及びＴＢの求め方）
ＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度を一定の条件とし、生成油中の硫黄濃度を一定の値ＳＰｎになるように、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度をＳＰｎとするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｎｘ＋ｂｎ（０＜ａｎである。）で表される直線が得られる。ａｎ及びｂｎは触媒の劣化挙動を反映した値である。この直線におけるｂｎが前記式７におけるＴＳＯＲ’となる。前記式７において、ＴＳＯＲ’にｂｎを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。なお、脱硫の活性化エネルギーＥａは、後述の方法により求められた値を使用することができる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｎ ’ｘ（|－ａｎ ’|＝ａｎ ’である。）で表される直線が得られる。ａｎ ’は、触媒の劣化速度を表している。
ＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度を一定の条件とし、生成油中の硫黄濃度を一定の値ＳＰｎになるように、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度をＳＰｎとするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｎｘ＋ｂｎ（０＜ａｎである。）で表される直線が得られる。ａｎ及びｂｎは触媒の劣化挙動を反映した値である。この直線におけるｂｎが前記式７におけるＴＳＯＲ’となる。前記式７において、ＴＳＯＲ’にｂｎを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。なお、脱硫の活性化エネルギーＥａは、後述の方法により求められた値を使用することができる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｎ ’ｘ（|－ａｎ ’|＝ａｎ ’である。）で表される直線が得られる。ａｎ ’は、触媒の劣化速度を表している。
【0047】
ｎ種類の硫黄濃度ＳＰｎについて同様の反応を行い、同様にして、ｎ個のａｎ、ｂｎを求め、上記と同様の方法により、ｎ個のａｎ
’を得る。ｎは３以上の整数である。ｎの数が多い方が、より精度の高いＥｃを得ることができる。一方、ｎの数が多すぎると、Ｅｃを得るために時間が掛かり効率的ではない。本実施形態においては、ｎは３～２０であることが好ましく、３～１０であることがより好ましい。
このようにして得られたｎ個のａｎ ’及びｂｎを下式８にそれぞれ代入する。下式８は、コークの活性化エネルギー、及び基準反応温度を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
このようにして得られたｎ個のａｎ ’及びｂｎを下式８にそれぞれ代入する。下式８は、コークの活性化エネルギー、及び基準反応温度を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
【0048】
ｌｎ（ａｎ
’）＝ｌｎ（Ａ）－（Ｅｃ／Ｒｂｎ） 式８
前記式８中、Ａは頻度因子であり、Ｅｃはコーク劣化の活性化エネルギー（ｋＪ／ｍｏｌ）であり、Ｒは気体定数：０．００８３１（ｋＪ／（ｍｏｌ・Ｋ））である。
前記式８中、Ａは頻度因子であり、Ｅｃはコーク劣化の活性化エネルギー（ｋＪ／ｍｏｌ）であり、Ｒは気体定数：０．００８３１（ｋＪ／（ｍｏｌ・Ｋ））である。
【0049】
ｎ個のａｎ
’及びｂｎの組み合わせに関し、ｌｎ（ａｎ
’）を縦軸に、１／ｂｎを横軸にプロットし、回帰直線を引き、その傾きを求める。この傾きはＥｃ／Ｒであるため、傾きからＲを除すことによりコーク劣化の活性化エネルギーＥｃを求めることができる。
【0050】
また、ｎ個のｂｎを平均することによりＴＢを求めることができる。
【0051】
Ｅｃ及びＴＢを求める上でのＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのようなＬＨＳＶとしては、例えば０．５～２．０時間－１であり、水素分圧としては、例えば１４～１７ＭＰａであり、水素／原料油比としては、例えば５００～１３００［Ｎｍ３／ｋＬ］であり、原料油中の硫黄濃度としては、例えば３～５質量％である。
ｎ種類の硫黄濃度ＳＰｎも同様に実機運転条件に即した条件とすることが好ましい。このようなＳＰｎとしては、例えば０．１質量％以下である。
反応期間としては、例えば３００～２５００日である。
そのようなＬＨＳＶとしては、例えば０．５～２．０時間－１であり、水素分圧としては、例えば１４～１７ＭＰａであり、水素／原料油比としては、例えば５００～１３００［Ｎｍ３／ｋＬ］であり、原料油中の硫黄濃度としては、例えば３～５質量％である。
ｎ種類の硫黄濃度ＳＰｎも同様に実機運転条件に即した条件とすることが好ましい。このようなＳＰｎとしては、例えば０．１質量％以下である。
反応期間としては、例えば３００～２５００日である。
【0052】
（ＰＢ及びａの求め方）
ＬＨＳＶ、水素／原料油比、原料油中の硫黄濃度、シリコン濃度、生成油中の硫黄濃度を一定の条件とし、水素分圧Ｐｍの条件で、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度を一定の値とするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｍｘ＋ｂｍ（０＜ａｍである。）で表される直線が得られる。この直線におけるｂｍが前記式７におけるＴＳＯＲ’となる。前記式７において、ＴＳＯＲ’にｂｍを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｍ ’ｘ（|―ａｍ ’|＝ａｍ ’である。）で表される直線が得られる。ａｍ ’は、触媒の劣化速度を表している。
ＬＨＳＶ、水素／原料油比、原料油中の硫黄濃度、シリコン濃度、生成油中の硫黄濃度を一定の条件とし、水素分圧Ｐｍの条件で、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度を一定の値とするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｍｘ＋ｂｍ（０＜ａｍである。）で表される直線が得られる。この直線におけるｂｍが前記式７におけるＴＳＯＲ’となる。前記式７において、ＴＳＯＲ’にｂｍを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｍ ’ｘ（|―ａｍ ’|＝ａｍ ’である。）で表される直線が得られる。ａｍ ’は、触媒の劣化速度を表している。
【0053】
ｍ種類の水素分圧Ｐｍについて同様の反応を行い、同様にして、ｍ個のａｍ、ｂｍを求め、上記と同様の方法により、ｍ個のａｍ
’を得る。ｍは３以上の整数である。ｍの数が多い方が、より精度の高いａを得ることができる。一方、ｍの数が多すぎると、ａを得るために時間が掛かり効率的ではない。本実施形態においては、ｍは３～２０であることが好ましく、３～１０であることがより好ましい。
このようにして得られたｍ個のａｍ’及びＰｍを下式９にそれぞれ代入する。下式９は、水素分圧係数、及び基準水素分圧を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
このようにして得られたｍ個のａｍ’及びＰｍを下式９にそれぞれ代入する。下式９は、水素分圧係数、及び基準水素分圧を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
【0054】
ｌｎ（ａｍ
’）＝－ａｌｎ（Ｐｍ）＋Ｂ１ 式９
前記式９中、Ｂ１は０とすることができる。
前記式９中、Ｂ１は０とすることができる。
【0055】
ｍ個のａｍ’及びＰｍの組み合わせに関し、ｌｎ（ａｍ
’）を縦軸に、ｌｎ（Ｐｍ）を横軸にプロットし、回帰直線を引き、その傾きを求める。この傾きが水素分圧係数ａである。
【0056】
また、上記ｍ個の水素分圧Ｐｍを平均することによりＰＢを求めることができる。
【0057】
ａ及びＰＢを求める上でのＬＨＳＶ、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのようなＬＨＳＶとしては、例えば０．５～２．０時間－１であり、水素／原料油比としては、例えば５００～１３００［Ｎｍ３／ｋＬ］であり、原料油中の硫黄濃度としては、例えば３～５質量％であり、生成油中の硫黄濃度としては、例えば０．１質量％以下である。
ｍ種類の水素分圧Ｐｍも同様に実機運転条件に即した条件とすることが好ましい。このようなＰｍとしては、例えば１４～１７ＭＰａである。反応期間としては、例えば３００～２５００日である。
そのようなＬＨＳＶとしては、例えば０．５～２．０時間－１であり、水素／原料油比としては、例えば５００～１３００［Ｎｍ３／ｋＬ］であり、原料油中の硫黄濃度としては、例えば３～５質量％であり、生成油中の硫黄濃度としては、例えば０．１質量％以下である。
ｍ種類の水素分圧Ｐｍも同様に実機運転条件に即した条件とすることが好ましい。このようなＰｍとしては、例えば１４～１７ＭＰａである。反応期間としては、例えば３００～２５００日である。
【0058】
（α及びＴＳＯＲの求め方１）
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。想定する実機運転条件とは、（Ｅｃ及びＴＢの求め方）、（ＰＢ及びａの求め方）で説明した運転条件が例として挙げられる。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａαｘ＋ｂαで表される直線が得られる。この直線におけるａαが前記式４におけるαに相関した値であり、ｂαが前記式４におけるＴＳＯＲとなる。前記式７において、ＴＳＯＲ’にｂαを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａα ’ｘ（|－ａα ’|＝ａα ’とする。）で表される直線が得られる。ａα ’は、触媒の劣化速度を表している。
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。想定する実機運転条件とは、（Ｅｃ及びＴＢの求め方）、（ＰＢ及びａの求め方）で説明した運転条件が例として挙げられる。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａαｘ＋ｂαで表される直線が得られる。この直線におけるａαが前記式４におけるαに相関した値であり、ｂαが前記式４におけるＴＳＯＲとなる。前記式７において、ＴＳＯＲ’にｂαを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａα ’ｘ（|－ａα ’|＝ａα ’とする。）で表される直線が得られる。ａα ’は、触媒の劣化速度を表している。
【0059】
運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（すなわち、ｂα）を前記式４に代入してＤを求める。求められたＤを前記式３に代入するとΦＣが得られる。この場合、ΦＣはαの関数となる。前記式２において、ΦＳｉ＝１とすると、Φ＝ΦＣとなり、Φはαの関数となる。反応時間を横軸に、Φの対数を縦軸にプロットし、αを０＜αとなるように変えて回帰直線を引くと、ｙ＝－ａα
”ｘ（|－ａα
”|＝ａα
”である。）で表される複数の直線がαの値ごとに得られる。各直線におけるａα
”と、上述のａα
’が等しい値となるときのαを前記式４におけるαとすることができる。
【0060】
（α及びＴＳＯＲの求め方２）
（α及びＴＳＯＲの求め方１）と同様の反応を行い、ｙ＝ａαｘ＋ｂαで表される直線を得、ｂαを前記式４におけるＴＳＯＲとする。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（すなわち、ｂα）を前記式４に代入してＤを求める。求められたＤを前記式３に代入するとΦＣが得られる。この場合、ΦＣはαの関数となる。前記式２において、ΦSi＝１とすると、Φ＝ΦＣとなり、Φはαの関数となる。得られたΦを前記式７のΦ’に、ＴＳＯＲ（すなわち、ｂα）を前記式７のＴＳＯＲ ’に代入して、Ｔｔ’について整理すると、Ｔｔ’はαの関数となる。実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのαを前記式４におけるαとすることができる。同様にＮ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのαを前記式４におけるαとすることが好ましい。Ｎは１０以上の整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
（α及びＴＳＯＲの求め方１）と同様の反応を行い、ｙ＝ａαｘ＋ｂαで表される直線を得、ｂαを前記式４におけるＴＳＯＲとする。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（すなわち、ｂα）を前記式４に代入してＤを求める。求められたＤを前記式３に代入するとΦＣが得られる。この場合、ΦＣはαの関数となる。前記式２において、ΦSi＝１とすると、Φ＝ΦＣとなり、Φはαの関数となる。得られたΦを前記式７のΦ’に、ＴＳＯＲ（すなわち、ｂα）を前記式７のＴＳＯＲ ’に代入して、Ｔｔ’について整理すると、Ｔｔ’はαの関数となる。実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのαを前記式４におけるαとすることができる。同様にＮ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのαを前記式４におけるαとすることが好ましい。Ｎは１０以上の整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
【0061】
（コーク劣化関数１の変形例）
コーク劣化関数１の変形例を以下に説明する。コーク劣化関数１としては下式１０で表されるコーク劣化関数１－１を使用してもよい。
コーク劣化関数１の変形例を以下に説明する。コーク劣化関数１としては下式１０で表されるコーク劣化関数１－１を使用してもよい。
【0062】
ΦＣ＝ｅｘｐ（－Ｄ’ｔ） 式１０
前記式１０中、Ｄ’は触媒の活性種のコークの劣化係数であり、ｔは反応経過日数（日）である。
前記式１０中、Ｄ’は触媒の活性種のコークの劣化係数であり、ｔは反応経過日数（日）である。
【0063】
Ｄ’は下式１１により、求めることができる。
【0064】
【数4】
【0065】
前記式１１中のＳＦ、ＳＰ、ＬＨＳＶ、Ｐは、前記式４の説明と同様である。前記式１１中、Ｇは、上述した情報取得ステップで取得した原料油に関する情報及び運転条件に関する情報に基づき代入される値である。具体的には、Ｇは、水素の供給量（Ｎｍ３／時間）を原料油の供給量（ｋＬ／時間）で除すことにより求めることができる。
【0066】
本明細書において、「基準水素／原料油比」とは、実際の反応条件の標準的な水素／原料油比を意味する。後述する水素／原料油比係数ｂを決定する際に用いる水素／原料油比の平均値として求められる。
【0067】
前記式１１において、（ＧＢ／Ｇ）ｂで表される項は、水素／原料油比依存性を示す項であり、上述した通り、Ｇを設定値とし、一定の条件で運転する場合、定数となる。
【0068】
ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式４で説明した方法と同じ方法で求めることができる。以下、ＧＢ、ｂ、α’の求め方を説明する。
【0069】
ＧＢ、ｂはＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲと同様に、使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。以下、ＧＢ、ｂの求め方について説明する。
【0070】
（ＧＢ、ｂの求め方）
本実施形態の基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。
本実施形態の基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。
【0071】
ＬＨＳＶ、水素分圧、原料油中の硫黄濃度、生成油中の硫黄濃度を一定の条件とし、水素／原料油比Ｇｈの条件で、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度を一定の値とするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｈｘ＋ｂｈ（０＜ａｈである。）で表される直線が得られる。この直線におけるｂｈが前記式７におけるＴＳＯＲ’となる。前記式７において、ＴＳＯＲ’にｂｈを代入し、Ｔｔ
’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｈ
’ｘ（|―ａｈ
’|＝ａｈ
’である。）で表される直線が得られる。ａｈ
’は、触媒の劣化速度を表している。
【0072】
ｈ種類の水素／原料油比Ｇｈについて同様の反応を行い、同様にして、ｈ個のａｈ、ｂｈを求め、上記と同様の方法により、ｈ個のａｈ
’を得る。ｈは３以上の整数である。ｈの数が多い方が、より精度の高いｂを得ることができる。一方、ｈの数が多すぎると、ｂを得るために時間が掛かり効率的ではない。本実施形態においては、ｈは３～２０であることが好ましく、３～１０であることがより好ましい。
このようにして得られたｈ個のａｈ’及びＧｈを下式１２にそれぞれ代入する。下式１２は、水素／原料油比係数、及び基準水素／原料油比を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
このようにして得られたｈ個のａｈ’及びＧｈを下式１２にそれぞれ代入する。下式１２は、水素／原料油比係数、及び基準水素／原料油比を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
【0073】
ｌｎ（ａｈ’）＝－ｂｌｎ（Ｇｈ）＋Ｂ３ 式１２
前記式１２中、Ｂ３は０とすることができる。
前記式１２中、Ｂ３は０とすることができる。
【0074】
ｈ個のａｈ’及びＧｈの組み合わせに関し、ｌｎ（ａｈ’）を縦軸に、ｌｎ（Ｇｈ）を横軸にプロットし、回帰直線を引き、その傾きを求める。この傾きが水素／原料油比係数ｂである。
【0075】
また、上記ｈ種類の水素／原料油比Ｇｈを平均することにより、ＧＢを求めることができる。
【0076】
ｂ及びＧＢを求める上でのＬＨＳＶ、水素分圧、原料油中の硫黄濃度、生成油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのようなＬＨＳＶとしては、例えば０．５～２．０時間－１であり、水素分圧としては、例えば１４～１７ＭＰａであり、原料油中の硫黄濃度としては、例えば３～５質量％であり、生成油中の硫黄濃度としては、例えば０．１質量％以下である。
ｈ種類の水素／原料油比Ｇｈも同様に実機運転条件に即した条件とすることが好ましい。このようなＧｈとしては、例えば５００～１３００［Ｎｍ３／ｋＬ］である。反応期間としては、例えば３００～２５００日である。
そのようなＬＨＳＶとしては、例えば０．５～２．０時間－１であり、水素分圧としては、例えば１４～１７ＭＰａであり、原料油中の硫黄濃度としては、例えば３～５質量％であり、生成油中の硫黄濃度としては、例えば０．１質量％以下である。
ｈ種類の水素／原料油比Ｇｈも同様に実機運転条件に即した条件とすることが好ましい。このようなＧｈとしては、例えば５００～１３００［Ｎｍ３／ｋＬ］である。反応期間としては、例えば３００～２５００日である。
【0077】
また、α’は、前記式４の代わりに前記式１１を使用する以外は、コーク劣化関数１において説明したαの求め方（（α及びＴＳＯＲの求め方１）及び（α及びＴＳＯＲの求め方２））と同じ方法により求めることができる。
【0078】
＜コーク劣化関数２＞
本実施形態においては、コークの堆積による触媒の劣化に関するコーク劣化関数は、触媒の易失活活性種の劣化に関する易失活活性種劣化関数と、触媒の難失活活性種の劣化に関する難失活活性種劣化関数から構成される下式１３で表されるコーク劣化関数２が好ましい。
本実施形態においては、コークの堆積による触媒の劣化に関するコーク劣化関数は、触媒の易失活活性種の劣化に関する易失活活性種劣化関数と、触媒の難失活活性種の劣化に関する難失活活性種劣化関数から構成される下式１３で表されるコーク劣化関数２が好ましい。
【0079】
ΦＣ＝ｋ1×ｅｘｐ（－Ｄ１ｔ）＋ｋ２×ｅｘｐ（－Ｄ２ｔ） 式１３
前記式１３中、ｋ１は触媒の易失活活性種の活性点係数であり、ｋ２は触媒の難失活活性種の活性点係数であり、活性点係数は両活性種の相対反応速度定数を表している。Ｄ１は触媒の易失活活性種のコークによる劣化係数であり、Ｄ２は触媒の難失活活性種のコークによる劣化係数であり、ｔは反応経過日数（日）であり、ｋ１＋ｋ２＝１である。
前記式１３中、ｋ１は触媒の易失活活性種の活性点係数であり、ｋ２は触媒の難失活活性種の活性点係数であり、活性点係数は両活性種の相対反応速度定数を表している。Ｄ１は触媒の易失活活性種のコークによる劣化係数であり、Ｄ２は触媒の難失活活性種のコークによる劣化係数であり、ｔは反応経過日数（日）であり、ｋ１＋ｋ２＝１である。
【0080】
重質コーカー分解軽油の水素化処理反応においては、上述した通り、コークの堆積により触媒が劣化するため、生成油中の硫黄含有量を一定レベル以下に保つため、反応温度を上げて運転を行う必要がある。重質コーカー分解軽油の水素化処理反応において、反応開始初期には、反応温度が急激に上昇する。この反応温度の急激な上昇は、反応開始初期の触媒の急激な劣化を意味する。一方、反応中期以降においては、反応温度は緩やかに上昇する。この反応温度の緩やかな上昇は、反応中期以降の触媒の緩やかな劣化を意味する。
【0081】
すなわち、重質コーカー分解軽油の水素化処理反応においては、反応開始初期の触媒の急激な劣化と、反応中期以降の触媒の緩やかな劣化が起こっていることが、反応時間に対する反応温度のプロファイルから示唆される。
【0082】
本願の発明者らは、上記反応時間に対する反応温度のプロファイルに基づき、触媒には、反応開始初期にコークにより失活する易失活活性種と、反応中期以降にコークにより失活する難失活活性種が存在するとの仮定に基づき、コーク劣化関数１をさらに改良し、コーク劣化関数２を見出した。その結果、コーク劣化関数２によれば、コーク劣化関数１よりも、より精度よく触媒のコークによる劣化度を算出可能となることを見出した。易失活活性種とは、主に反応初期に活性が失われる活性種であり、難失活活性種とは、反応中期以降に活性が失われる活性種を意味する。
【0083】
前記式１３において、ｋ１は、触媒の易失活活性種の活性点係数を示し、ｋ２は、触媒の難失活活性種の活性点係数を示す。ｋ１及びｋ２は触媒固有の定数であり、その求め方は後述する。
【0084】
Ｄ１は下式１４により、Ｄ２は下式１５により求めることができる。
【0085】
【数5】
【0086】
【数6】
【0087】
前記式１４及び前記式１５中、ＳＦは任意の反応ｔ日経過時の原料油中の硫黄濃度（質量％）であり、ＳＰは任意の反応ｔ日経過時の生成油中の硫黄濃度（質量％）であり、ｎは、重質コーカー分解軽油を含む原料油の水素化処理反応の反応次数であり、ＬＨＳＶは任意の反応ｔ日経過時の液空間速度（時間－１）であり、ＰＢは基準水素分圧（ＭＰａ）であり、Ｐは任意の反応ｔ日経過時の水素分圧（ＭＰａ）であり、ａは水素分圧係数であり、Ｅｃはコーク劣化の活性化エネルギー（ｋＪ／ｍｏｌ）であり、Ｒは気体定数：０．００８３１（ｋＪ／（ｍｏｌ・Ｋ））であり、ＴＢは基準反応温度（Ｋ）であり、ＴＳＯＲは０日目の要求温度（Ｋ）である。
前記式１４中、α１は易失活活性点の触媒定数（触媒のコークによる劣化速度を表す定数）であり、前記式１５中、α２は難失活活性点の触媒定数（触媒のコークによる劣化速度を表す定数）である。
前記式１４中、α１は易失活活性点の触媒定数（触媒のコークによる劣化速度を表す定数）であり、前記式１５中、α２は難失活活性点の触媒定数（触媒のコークによる劣化速度を表す定数）である。
【0088】
前記式１４及び前記式１５中、ＳＦ、ＳＰ、ＬＨＳＶ、Ｐは、前記式４の説明と同様に、上述した情報取得ステップで取得した原料油に関する情報、生成油に関する情報、及び運転条件に関する情報に基づき代入される値である。なお、ＬＨＳＶは原料油の供給量（体積／時間）を触媒充填量（体積）で除すことにより求めることができる。
【0089】
上述した通り、ＳＦ、ＬＨＳＶ、Ｐは制御することが可能なパラメータである。また、ＳＰは目的とする生成油の硫黄濃度である。すなわち、前記式１４及び前記式１５により、上記ＳＦ、ＳＰ、ＬＨＳＶ、及びＰの反応条件における、触媒の易失活活性種のコークによる劣化係数、及び触媒の難失活活性種のコークによる劣化係数をそれぞれ算出することができる。重質コーカー分解軽油を含む原料油の水素化処理反応の反応次数であるｎの求め方は後述する。
【0090】
（コークの基本劣化パラメータの求め方）
前記式１４及び前記式１５中、α１、α２、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式４と同様、定数であり、これらのパラメータを総称して「コークの基本劣化パラメータ２」という。コークの基本劣化パラメータ２は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。前記式１４及び前記式１５中、ＰＢ、ａ、Ｅｃ、ＴＢは前記式４と同じ方法により求めることができる。
一方、α１、α２、ＴＳＯＲは例えば以下の２種類の方法により求めることができる。また、前記式１３における触媒の易失活活性種の活性点数係数ｋ１、触媒の難失活活性種の活性点数係数ｋ２も同時に以下のように求めることができる。
前記式１４及び前記式１５中、α１、α２、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式４と同様、定数であり、これらのパラメータを総称して「コークの基本劣化パラメータ２」という。コークの基本劣化パラメータ２は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。前記式１４及び前記式１５中、ＰＢ、ａ、Ｅｃ、ＴＢは前記式４と同じ方法により求めることができる。
一方、α１、α２、ＴＳＯＲは例えば以下の２種類の方法により求めることができる。また、前記式１３における触媒の易失活活性種の活性点数係数ｋ１、触媒の難失活活性種の活性点数係数ｋ２も同時に以下のように求めることができる。
【0091】
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方１）
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。反応時間を横軸に、反応温度を縦軸にプロットする。これらのプロットの回帰直線を引く場合、上述した通り、反応開始初期（ｘ１～ｘｎ）の、急激な反応温度の上昇に相関するｙ＝ａ１ｘ＋ｂ１で表される直線と、反応中期以降（ｘｎ＋１～ｘｍ）の、緩やかな反応温度の上昇に相関するｙ＝ａ２ｘ＋ｂ２で表される２本の直線が得られる。前記式において、ａ１＞ａ２＞０であり、ｂ１＜ｂ２であり、ｘ１＜ｘｎ＜ｘｎ＋１＜ｘｍである。ｘｎ、ｘｍは反応開始からのｎ（ｍ）番目のプロットに相関した反応時間を意味する。
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。反応時間を横軸に、反応温度を縦軸にプロットする。これらのプロットの回帰直線を引く場合、上述した通り、反応開始初期（ｘ１～ｘｎ）の、急激な反応温度の上昇に相関するｙ＝ａ１ｘ＋ｂ１で表される直線と、反応中期以降（ｘｎ＋１～ｘｍ）の、緩やかな反応温度の上昇に相関するｙ＝ａ２ｘ＋ｂ２で表される２本の直線が得られる。前記式において、ａ１＞ａ２＞０であり、ｂ１＜ｂ２であり、ｘ１＜ｘｎ＜ｘｎ＋１＜ｘｍである。ｘｎ、ｘｍは反応開始からのｎ（ｍ）番目のプロットに相関した反応時間を意味する。
【0092】
上述のａ１がα１に相関する値であり、ｂ１がｋ１＋ｋ２に相関する値であり、ＴＳＯＲである。また、ａ２がα２に相関する値であり、ｂ２がｋ２に相関する値である。次に、ｙ＝ａ１ｘ＋ｂ１とｙ＝ａ２ｘ＋ｂ２での交点（ｘｉｐ、ｙｉｐ）を計算する。この交点はｙ＝ａ１ｘ＋ｂ１とｙ＝ａ２ｘ＋ｂ２の変曲点を意味する。すなわち、（ｘｉｐ、ｙｉｐ）までは、触媒の易失活活性種のみが存在するとみなし、（ｘｉｐ、ｙｉｐ）からは触媒の易失活活性種及び触媒の難失活活性種の２つの活性種が存在するとみなす。
【0093】
前記式７において、ＴＳＯＲ’にｂ１を代入し、Ｔｔ’に反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、ｘ１～ｘｉｐまでの回帰直線を引くと、ｙ＝－ａ１
’ｘ（|－ａ１
’|＝ａ１
’である。）で表される直線が得られる。また、ｘiｐ～ｘｍまでの回帰直線を引くと、ｙ＝－ａ２
’ｘ－ｂ２
’（|－ａ２
’|＝ａ２
’であり、|－ｂ２
’|＝ｂ２
’である。）で表される直線が得られる。ａ１
’はα１に相関した値であり、ａ２
’はα２に相関した値であり、ｂ２
’はｋ２に相関した値である。
【0094】
ｋ２は、上述の回帰直線で得られたｂ２
’を下式１６に代入して求めることができる。ｋ１はｋ１＋ｋ２＝１であるから、ｋ１＝１－ｋ２より求めることができる。
ｋ２＝ｅｘｐ（－ｂ２ ’） 式１６
ｋ２＝ｅｘｐ（－ｂ２ ’） 式１６
【0095】
運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（ｂ１）を前記式１４及び前記式１５に代入してＤ１及びＤ２を得る。求められたＤ１、Ｄ２、ｋ１、ｋ２を前記式１３に代入するとΦＣが得られる。この場合、ΦＣはα１とα２の関数となる。前記式２において、ΦＳｉ＝１とすると、Φ＝ΦＣとなり、Φはα１とα２の関数となる。反応時間を横軸に、Φの対数を縦軸にプロットし、α２＝０としてα１を０＜α１になるように変えてｘ１～ｘｉｐまでの回帰直線を引くと、ｙ＝－ａα１
”ｘ（|－ａα１
“|＝ａα１
”である。）で表される複数の直線がα１の値ごとに得られる。各直線におけるａα１
”と、上述のａ１
’が等しい値となるときのα１を前記式１４におけるα１とすることができる。
次に、上述の方法で得られたα１とα２の関数であるΦＣ（Φ）に、得られたα１を代入し、α２を０＜α２になるように変えてｘiｐ～ｘｍまでの回帰直線を引くと、ｙ＝－ａα２ ”ｘ－ｂα２ “（|－ａα２ ”|＝ａα２ “である。）で表される複数の直線がα２の値ごとに得られる。各直線におけるａα２と、上述のａ２ ’が等しい値となるときのα２を前記式１５におけるα２とすることができる。
次に、上述の方法で得られたα１とα２の関数であるΦＣ（Φ）に、得られたα１を代入し、α２を０＜α２になるように変えてｘiｐ～ｘｍまでの回帰直線を引くと、ｙ＝－ａα２ ”ｘ－ｂα２ “（|－ａα２ ”|＝ａα２ “である。）で表される複数の直線がα２の値ごとに得られる。各直線におけるａα２と、上述のａ２ ’が等しい値となるときのα２を前記式１５におけるα２とすることができる。
【0096】
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方２）
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方１）と同様の反応を行い、ＴＳＯＲ（ｂ１）を得る。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（ｂ１）を前記式１４及び前記式１５に代入してＤ１及びＤ２を得る。求められたＤ１、Ｄ２を前記式１３に代入するとΦＣが得られる。この場合、ΦＣはα１、α２、ｋ１、ｋ２の関数となる。前記式２において、ΦＳｉ＝１とすると、Φ＝ΦＣとなり、Φはα１、α２、ｋ１、ｋ２の関数となる。得られたΦを前記式７のΦ’に、ＴＳＯＲ（すなわち、ｂ１）を前記式７のＴＳＯＲ ’に代入して、Ｔｔ’について整理すると、Ｔｔ’はα１、α２、ｋ１、ｋ２の関数となる。α２を０とし、ｋ２＝１－ｋ１とするとＴｔ’はα１、ｋ１の関数となる。ｘ１～ｘｉｐにおける実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのα１、ｋ１の組み合わせを求め、このα１を前記式１４におけるα１とすることができる。なお、この時のｋ１は仮値である。同様にＭ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα１を前記式１４におけるα１とすることが好ましい。Ｍは１０以上の整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
得られたα１を使用し、ｋ１＝１－ｋ２とするとＴｔ’はα２、ｋ２の関数となる。ｘiｐ～ｘｍにおける実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのα２、ｋ２の組み合わせを求め、このα２、ｋ２を前記式１５におけるα２、ｋ２とすることができる。得られたｋ２をｋ１＝１－ｋ２に代入することによりｋ１を求め、このｋ１を前記式１４におけるｋ１とすることができる。同様にＬ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα２、ｋ２を前記式１５におけるα２、ｋ２とすることが好ましい。また、得られたｋ２から求められたｋ１を前記式１４におけるｋ１とすることが好ましい。Ｌは１０以上の整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方１）と同様の反応を行い、ＴＳＯＲ（ｂ１）を得る。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（ｂ１）を前記式１４及び前記式１５に代入してＤ１及びＤ２を得る。求められたＤ１、Ｄ２を前記式１３に代入するとΦＣが得られる。この場合、ΦＣはα１、α２、ｋ１、ｋ２の関数となる。前記式２において、ΦＳｉ＝１とすると、Φ＝ΦＣとなり、Φはα１、α２、ｋ１、ｋ２の関数となる。得られたΦを前記式７のΦ’に、ＴＳＯＲ（すなわち、ｂ１）を前記式７のＴＳＯＲ ’に代入して、Ｔｔ’について整理すると、Ｔｔ’はα１、α２、ｋ１、ｋ２の関数となる。α２を０とし、ｋ２＝１－ｋ１とするとＴｔ’はα１、ｋ１の関数となる。ｘ１～ｘｉｐにおける実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのα１、ｋ１の組み合わせを求め、このα１を前記式１４におけるα１とすることができる。なお、この時のｋ１は仮値である。同様にＭ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα１を前記式１４におけるα１とすることが好ましい。Ｍは１０以上の整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
得られたα１を使用し、ｋ１＝１－ｋ２とするとＴｔ’はα２、ｋ２の関数となる。ｘiｐ～ｘｍにおける実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのα２、ｋ２の組み合わせを求め、このα２、ｋ２を前記式１５におけるα２、ｋ２とすることができる。得られたｋ２をｋ１＝１－ｋ２に代入することによりｋ１を求め、このｋ１を前記式１４におけるｋ１とすることができる。同様にＬ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα２、ｋ２を前記式１５におけるα２、ｋ２とすることが好ましい。また、得られたｋ２から求められたｋ１を前記式１４におけるｋ１とすることが好ましい。Ｌは１０以上の整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
【0097】
このように、α１、α２、ＴＳＯＲ、ｋ１、及びｋ２を求めるためには、ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２を得る必要がある。ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２は例えば、以下のように求めることができる。
【0098】
上述の方法で、反応時間を横軸に、反応温度を縦軸にプロットする。反応開始から反応終了までの回帰直線を引くと、ｙ＝ａ’ｘ＋ｂ’で表される直線が得られる。この直線は、変曲点を考慮していない直線である。反応終了からプロットを順に削除し、相関係数が１に近づくように前記ｙ＝ａ’ｘ＋ｂ’を調整し、ｙ＝ａ１’ｘ＋ｂ１’を得る。同様に、反応開始からプロットを順に削除し、相関係数が１に近づくように前記ｙ＝ａ’ｘ＋ｂ’を調整し、ｙ＝ａ２’ｘ＋ｂ２’を得る。ｙ＝ａ１’ｘ＋ｂ１’とｙ＝ａ２’ｘ＋ｂ２’の相関係数の平均が最も１に近づく場合の各直線がｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２となる。なお、全プロットは、ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２のどちらか一方に属するようにする。
【0099】
ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２を得るために必要な反応時間は通常１００日以上である。また、一般に、上記、ｙ＝ａ’ｘ＋ｂ’の相関関数が０．５以上となるまで反応を行えば充分である。
【0100】
（コーク劣化関数２の変形例）
コーク劣化関数２の変形例を以下に説明する。コーク劣化関数２としては下式１７で表されるコーク劣化関数２－１を使用してもよい。
コーク劣化関数２の変形例を以下に説明する。コーク劣化関数２としては下式１７で表されるコーク劣化関数２－１を使用してもよい。
【0101】
ΦＣ＝ｋ1×ｅｘｐ（－Ｄ１’ｔ）＋ｋ２×ｅｘｐ（－Ｄ２’ｔ） 式１７
前記式１７中、ｋ１、ｋ２、ｔは前記式１３と同じであり、Ｄ１’は触媒の易失活活性種のコークの劣化係数であり、Ｄ２’は触媒の難失活活性種のコークの劣化係数である。
前記式１７中、ｋ１、ｋ２、ｔは前記式１３と同じであり、Ｄ１’は触媒の易失活活性種のコークの劣化係数であり、Ｄ２’は触媒の難失活活性種のコークの劣化係数である。
【0102】
Ｄ１’は下式１８により、Ｄ２’は下式１９により求めることができる。
【0103】
【数7】
【0104】
【数8】
【0105】
前記式１８及び前記式１９中、ＳＦ、ＳＰ、ｎ、ＬＨＳＶ、ＰＢ、Ｐ、ａ、Ｅｃ、Ｒ、ＴＢ、ＴＳＯＲは前記式１４及び前記式１５と同じであり、ＧＢ、Ｇ、ｂは前記式１１と同じであり、前記式１８中、α１’は易失活活性点の触媒定数（触媒のコークによる劣化速度を表す定数）であり、前記式１９中、α２’は難失活活性点の触媒定数（触媒のコークによる劣化速度を表す定数）である。
【0106】
ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式１４及び前記式１５で説明した方法と同じ方法で求めることができ、ＧＢ、ｂは、前記式１１で説明した方法と同じ方法で求めることができる。
【0107】
また、α１’、α２’は、前記式１３の代わりに前記式１７を、前記式１４の代わりに前記式１８を、前記式１５の代わりに前記式１９を使用する以外は、劣化関数２において説明したα１、α２の求め方と同じ方法により求めることができる。
【0108】
＜シリコン劣化関数＞
シリコン劣化関数としては、触媒のシリコン劣化に関する劣化度合いを算出可能な関数であれば、特に限定されないが、例えば下式２０で表される劣化関数が例として挙げられる。下式２０で表される劣化関数は、非特許文献（化学工学論文集、第２４巻、第４号（１９９８）、ｐ６５６）に記載されている重質油の拡散性劣化モデルにおけるメタル劣化をシリコン劣化に応用して本願の発明者らが初めて見出した関数である。
シリコン劣化関数としては、触媒のシリコン劣化に関する劣化度合いを算出可能な関数であれば、特に限定されないが、例えば下式２０で表される劣化関数が例として挙げられる。下式２０で表される劣化関数は、非特許文献（化学工学論文集、第２４巻、第４号（１９９８）、ｐ６５６）に記載されている重質油の拡散性劣化モデルにおけるメタル劣化をシリコン劣化に応用して本願の発明者らが初めて見出した関数である。
【0109】
【数9】
【0110】
前記式２０中、ηは、触媒の細孔内での反応物、すなわち硫黄化合物の拡散状況を表す。前記式２０は、シリコンの堆積によって生じる細孔閉塞に伴う触媒の劣化挙動を触媒有効係数ηで表すシリコン劣化関数である。η＝１のときは、拡散阻害がなく細孔内部まで有効に反応に使われることを示しており、小さくなるほど触媒の外表面近傍での反応になることを表している。前記式２０中のＳｉＯＣは下式２１により求めることができる。η、ＳｉＯＣ∞の求め方は後述する。
【0111】
【数10】
【0112】
前記式２１中、ＷＨＳＶ、Ｓｉｆ、Ｓｉｐは、上述した情報取得ステップで取得した原料油に関する情報、生成油に関する情報、及び運転条件に関する情報に基づき代入される値である。なお、ＷＨＳＶは、原料油の供給量（重量／時間）を触媒充填量（重量）で除すことにより求めることができる。
【0113】
上述した通り、ＷＨＳＶ、Ｓｉｆは制御することが可能なパラメータである。また、Ｓｉｐは目的とする生成油のシリコン濃度である。すなわち、前記式２０及び２１により、上記ＷＨＳＶ、Ｓｉｆ、Ｓｉｐの反応条件における、シリコン劣化関数を算出することができる。
【0114】
重質コーカー分解軽油を含む原料油の水素化処理反応においては、上述した通り、シリコンの堆積により触媒が劣化するため、生成油中の硫黄含有量を一定レベル以下に保つため、反応温度を上げて運転を行う必要がある。このシリコンの堆積による触媒の劣化の現象は、原料油の触媒細孔への拡散速度の減少により説明することができる。すなわち、重質コーカー分解軽油の水素化処理反応においては、原料油が触媒細孔内に拡散し、活性点にアクセスする必要がある。このとき、反応速度より拡散速度が充分に大きい場合には、原料油は全ての活性点にアクセスできるが、反応速度に対して拡散速度が小さい場合は、原料油は全ての活性点にアクセスできず、全ての活性点が有効に機能しない。拡散速度は、シリコンの堆積により触媒の細孔が閉塞することにより経時的に減少する。
【0115】
触媒細孔内の活性点がどの程度有効に使用されるかは、反応速度と拡散速度との関係により定まることになる。この関係は、触媒有効係数ηで理論的に表すことができる。拡散阻害の影響がなく、触媒細孔内の活性点が実質的に全て使用されている場合、ηは１となる。一方、触媒細孔内の活性点が実質的に全て使用されない場合、ηは１よりも小さくなる。そして、触媒の外表面近傍の活性点のみが使用されている場合は、ηは０に近づく。
【0116】
上述した通り、本願の発明者らは、反応開始初期にコークにより失活する易失活活性種と、反応中期以降にコークにより失活する難失活活性種が存在するとの仮定に基づき、コーク劣化関数を見出した。ここで、反応期間全般に渡ってシリコン堆積による活性低下も起こっていると考えられる。本願の発明者らは、コーク劣化関数とシリコン劣化関数を組み合わせることにより、より精度よく触媒の劣化度を算出可能となることを見出した。
【0117】
前記式２０中、η、ＳｉＯＣ∞は、定数である。以下、これらのパラメータを総称して「シリコンの基本劣化パラメータ」という。シリコンの基本劣化パラメータは、使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。
以下、基本劣化パラメータη、ＳｉＯＣ∞の求め方について例を説明するが、本発明はこれらに限定されるものではない。
以下、基本劣化パラメータη、ＳｉＯＣ∞の求め方について例を説明するが、本発明はこれらに限定されるものではない。
【0118】
（ＳｉＯＣ∞の求め方）
ＳｉＯＣ∞は触媒固有の値であり本分野で公知の方法に求めることができる。ＳｉＯＣ∞の求め方の例としては、例えば、実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、シリコン濃度、生成油中の硫黄濃度、シリコン濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。実機運転条件とは、後述の（ηの求め方１）に記載の劣化触媒を得るための反応条件が例として挙げられる。で説明する運転条件が例として挙げられる。長時間反応を継続すると、反応温度を機器使用最高温度としても（例えばＢＡＳＥ（全反応期間をｔｔｏｔａｌ（日）とした時に反応開始～０．３４×ｔｔｏｔａｌ（日）における平均反応温度）＋９０℃）、生成油中のシリコン濃度が所定の値とならなくなる。このまま反応を継続し、最終的に、原料油中のシリコン濃度＝生成油中のシリコン濃度となったときに反応を停止し、触媒を抜出す。抜き出した触媒のシリコン濃度を測定し、得られた値をＳｉＯＣ∞とすることができる。
ＳｉＯＣ∞は触媒固有の値であり本分野で公知の方法に求めることができる。ＳｉＯＣ∞の求め方の例としては、例えば、実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、シリコン濃度、生成油中の硫黄濃度、シリコン濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。実機運転条件とは、後述の（ηの求め方１）に記載の劣化触媒を得るための反応条件が例として挙げられる。で説明する運転条件が例として挙げられる。長時間反応を継続すると、反応温度を機器使用最高温度としても（例えばＢＡＳＥ（全反応期間をｔｔｏｔａｌ（日）とした時に反応開始～０．３４×ｔｔｏｔａｌ（日）における平均反応温度）＋９０℃）、生成油中のシリコン濃度が所定の値とならなくなる。このまま反応を継続し、最終的に、原料油中のシリコン濃度＝生成油中のシリコン濃度となったときに反応を停止し、触媒を抜出す。抜き出した触媒のシリコン濃度を測定し、得られた値をＳｉＯＣ∞とすることができる。
【0119】
（ηの求め方１）
ベンチスケールにおいて、新触媒及び前記新触媒に対応する劣化触媒を用いて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、シリコン濃度、生成油中の硫黄濃度、シリコン濃度となるように、一定期間反応を行う。劣化触媒としては、例えば、ＬＨＳＶ０．５～２時間－１、水素分圧１１～１７ＭＰａ、水素／原料油比５００～１３００［Ｎｍ３／ｋＬ］、原料油中の硫黄濃度３～５質量％、シリコン濃度０．３～７重量ｐｐｍ、生成油中の硫黄濃度０．１質量％未満、シリコン濃度０．３重量ｐｐｍ未満で、１００日以上反応を行った触媒を使用することができる。反応は実機で行っても、ベンチスケールで行ってもよい。
生成油の硫黄濃度が安定した数日間について、劣化触媒および新触媒の脱硫反応速度定数をそれぞれ後述の式２４に基づき求める。劣化触媒の脱硫反応速度定数を新触媒の触媒重量平均温度（反応温度）で補正した脱硫反応速度定数を算出する。具体的には、下式Ａのｋ１として算出する。下式Ａは、前記式６で表されるアレニウスの式に基づいた式である。
ベンチスケールにおいて、新触媒及び前記新触媒に対応する劣化触媒を用いて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、シリコン濃度、生成油中の硫黄濃度、シリコン濃度となるように、一定期間反応を行う。劣化触媒としては、例えば、ＬＨＳＶ０．５～２時間－１、水素分圧１１～１７ＭＰａ、水素／原料油比５００～１３００［Ｎｍ３／ｋＬ］、原料油中の硫黄濃度３～５質量％、シリコン濃度０．３～７重量ｐｐｍ、生成油中の硫黄濃度０．１質量％未満、シリコン濃度０．３重量ｐｐｍ未満で、１００日以上反応を行った触媒を使用することができる。反応は実機で行っても、ベンチスケールで行ってもよい。
生成油の硫黄濃度が安定した数日間について、劣化触媒および新触媒の脱硫反応速度定数をそれぞれ後述の式２４に基づき求める。劣化触媒の脱硫反応速度定数を新触媒の触媒重量平均温度（反応温度）で補正した脱硫反応速度定数を算出する。具体的には、下式Ａのｋ１として算出する。下式Ａは、前記式６で表されるアレニウスの式に基づいた式である。
【0120】
【数11】
すなわち、前記式ＡにＴ１、ｋ２、Ｔ２、Ｅ（脱硫の活性化エネルギーの測定方法については後述する）、Ｒを代入することにより、前記新触媒の脱硫反応速度定数を測定したときの反応温度で補正した前記劣化触媒の脱硫反応速度定数ｋ１を求めることができる。
【0121】
新触媒の脱硫反応速度定数（ｋｆｒｅｓｈ）に対する前記補正した劣化触媒の脱硫反応速度定数（ｋｕｓｅｄ）である相対活性（ｋｕｓｅｄ／ｋｆｒｅｓｈ）を求める。
ｋｆｒｅｓｈ／ｋｆｒｅｓｈとｋｕｓｅｄ／ｋｆｒｅｓｈの差には、シリコンの堆積に加えてコークの堆積による活性低下分も含まれているので、その低下分をｋｕｓｅｄ／ｋｆｒｅｓｈに補正（加算）する必要がある。発明者らが検討を行った所、上述したように重質コーカー分解軽油を含む原料油の水素化処理反応における触媒のコークの堆積による劣化は、常圧蒸留残渣を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応（間接脱硫反応）における触媒のコークの堆積による劣化とほぼ同じであり、劣化挙動についてもカーボン堆積量に対して直線的であることが判明した。すなわち、（ｋｆｒｅｓｈ／ｋｆｒｅｓｈ－ｋｕｓｅｄ／ｋｆｒｅｓｈ）におけるコークの堆積による活性低下分は、シリコンを含まない減圧蒸留軽油を含む原料油の水素化処理反応により求めることができる。
ｋｆｒｅｓｈ／ｋｆｒｅｓｈとｋｕｓｅｄ／ｋｆｒｅｓｈの差には、シリコンの堆積に加えてコークの堆積による活性低下分も含まれているので、その低下分をｋｕｓｅｄ／ｋｆｒｅｓｈに補正（加算）する必要がある。発明者らが検討を行った所、上述したように重質コーカー分解軽油を含む原料油の水素化処理反応における触媒のコークの堆積による劣化は、常圧蒸留残渣を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応（間接脱硫反応）における触媒のコークの堆積による劣化とほぼ同じであり、劣化挙動についてもカーボン堆積量に対して直線的であることが判明した。すなわち、（ｋｆｒｅｓｈ／ｋｆｒｅｓｈ－ｋｕｓｅｄ／ｋｆｒｅｓｈ）におけるコークの堆積による活性低下分は、シリコンを含まない減圧蒸留軽油を含む原料油の水素化処理反応により求めることができる。
【0122】
前記新触媒を用いて減圧蒸留軽油を含む原料油の水素化処理反応を行い、間接脱硫劣化触媒を製造する。間接脱硫劣化触媒としては、例えば、ＬＨＳＶ０．５～２．５時間－１、水素分圧３～７ＭＰａ、水素／原料油比１００～４００［Ｎｍ３／ｋＬ］、原料油中の硫黄濃度０．５～３質量％、シリコン濃度０．３重量ｐｐｍ未満、生成油中の硫黄濃度０．３質量％未満、シリコン濃度０．３重量ｐｐｍ未満で、１００日以上反応を行った触媒を使用することができる。なお、反応時間は、上記劣化触媒と同程度の反応時間とすることが好ましい。反応は実機で行っても、ベンチスケールで行ってもよい。得られた間接脱硫劣化触媒を用いて上述の要領で、重質コーカー分解軽油を含む原料油の水素化処理反応における、脱硫反応速度定数（ｋｕｓｅｄ（ｃｏｋｅ））を求める。
コーク堆積量を横軸、コーク堆積による新触媒からの活性低下量（ｋｆｒｅｓｈ／ｋｆｒｅｓｈ（=１）－ｋｕｓｅｄ（ｃｏｋｅ）／ｋｆｒｅｓｈ）を縦軸にプロットし、プロットした点と原点とを結ぶ直線の傾きからカーボン堆積量１質量％あたりの活性低下量を算出する。炭素含有量が異なる劣化触媒としては、例えば同じ反応塔の上部と下部からそれぞれ抜き出したものや別の時期に運転した触媒などを用いることが考えられる。同じ銘柄の触媒であっても反応塔での充填箇所や運転期間や原料油性状によって触媒の炭素含有量が異なってくるためである。算出したカーボン堆積量１質量％あたりの活性低下量に劣化触媒に堆積したカーボン量を乗じることで、劣化触媒に堆積したカーボンによる活性低下量を算出する。
算出した堆積カーボンによる活性低下量を元の劣化触媒の相対活性に加算することで、シリコン堆積のみによる劣化触媒の相対活性を算出する。同じ新触媒を用いて（ＳｉＯＣ／ＳｉＯＣ∞）の異なる（すなわちＳｉＯＣの異なる）複数の劣化触媒に関し、上述の方法で、シリコン堆積のみによる劣化触媒の相対活性を算出する。ＳｉＯＣの異なる複数の劣化触媒としては、例えば同じ反応等塔の上部と下部からそれぞれ抜き出したものや別の時期に運転した触媒などを用いることが考えられる。同じ銘柄の触媒であっても反応塔での充填箇所や運転期間や原料油性状によって触媒のシリコン含有量が異なってくるためである。
横軸に劣化触媒のシリコン許容量に対するシリコン許容割合（ＳｉＯＣ／ＳｉＯＣ∞）、縦軸に新触媒の活性を１とした場合のシリコン堆積のみによる相対活性（ΦＳｉ）をプロットし、それらの点が最も近似する前記式２０で表されるシリコン劣化関数で用いられているηの値を求める。プロットの数（ＳｉＯＣ／ＳｉＯＣ∞の異なる劣化触媒の数）としては、２～５であることが好ましい。
コーク堆積量を横軸、コーク堆積による新触媒からの活性低下量（ｋｆｒｅｓｈ／ｋｆｒｅｓｈ（=１）－ｋｕｓｅｄ（ｃｏｋｅ）／ｋｆｒｅｓｈ）を縦軸にプロットし、プロットした点と原点とを結ぶ直線の傾きからカーボン堆積量１質量％あたりの活性低下量を算出する。炭素含有量が異なる劣化触媒としては、例えば同じ反応塔の上部と下部からそれぞれ抜き出したものや別の時期に運転した触媒などを用いることが考えられる。同じ銘柄の触媒であっても反応塔での充填箇所や運転期間や原料油性状によって触媒の炭素含有量が異なってくるためである。算出したカーボン堆積量１質量％あたりの活性低下量に劣化触媒に堆積したカーボン量を乗じることで、劣化触媒に堆積したカーボンによる活性低下量を算出する。
算出した堆積カーボンによる活性低下量を元の劣化触媒の相対活性に加算することで、シリコン堆積のみによる劣化触媒の相対活性を算出する。同じ新触媒を用いて（ＳｉＯＣ／ＳｉＯＣ∞）の異なる（すなわちＳｉＯＣの異なる）複数の劣化触媒に関し、上述の方法で、シリコン堆積のみによる劣化触媒の相対活性を算出する。ＳｉＯＣの異なる複数の劣化触媒としては、例えば同じ反応等塔の上部と下部からそれぞれ抜き出したものや別の時期に運転した触媒などを用いることが考えられる。同じ銘柄の触媒であっても反応塔での充填箇所や運転期間や原料油性状によって触媒のシリコン含有量が異なってくるためである。
横軸に劣化触媒のシリコン許容量に対するシリコン許容割合（ＳｉＯＣ／ＳｉＯＣ∞）、縦軸に新触媒の活性を１とした場合のシリコン堆積のみによる相対活性（ΦＳｉ）をプロットし、それらの点が最も近似する前記式２０で表されるシリコン劣化関数で用いられているηの値を求める。プロットの数（ＳｉＯＣ／ＳｉＯＣ∞の異なる劣化触媒の数）としては、２～５であることが好ましい。
【0123】
以上、コークの劣化関数、シリコン劣化関数、及びこれらの関数に用いられるパラメータ（コークの基本劣化パラメータ、シリコンの基本劣化パラメータ、易失活活性種の活性点数係数ｋ１、触媒の難失活活性種の活性点数係数ｋ２等）の求め方を説明した。以上の説明では、コークの劣化関数で用いられるパラメータと、シリコン劣化関数で用いられるパラメータを別々に求めたが、これらを同時にフィッティングにより求めてもよい。具体的には、想定する実機運転条件で、実機又はベンチスケールで反応を行い、反応温度の経時変化のデータを取得する。例えば、コーク劣化関数１を使用する場合、後述の式２３、前記式２、前記式３、前記式４、前記式２０を用いる。前記式４、前記式２０に反応条件を代入する。そして、後述の式２３のＴｔが、前記実測の反応温度に一致するような、前記式４におけるコーク劣化の基本劣化パラメータ及び前記式２０におけるシリコンの基本劣化パラメータを、フィッティングで求めてもよい。
【0124】
≪反応温度算出ステップ≫
本実施形態の反応温度算出ステップは、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出するステップである。反応温度は、アレニウスの式に基づいた劣化速度式により算出することが好ましい。
本実施形態の反応温度算出ステップは、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出するステップである。反応温度は、アレニウスの式に基づいた劣化速度式により算出することが好ましい。
【0125】
＜劣化速度式＞
劣化速度式は前記式６で表されるアレニウスの式に基づいた式である。前記式７の算出方法と同様に、前記式１及び前記式６より下式２２が、導かれる。
劣化速度式は前記式６で表されるアレニウスの式に基づいた式である。前記式７の算出方法と同様に、前記式１及び前記式６より下式２２が、導かれる。
【0126】
【数12】
前記式２２を整理すると下式２３が得られる。
【0127】
【数13】
【0128】
前記式２３に、上述の方法で得られたＴＳＯＲ、Φを代入することにより、任意の反応ｔ日経過時におけるＴｔ（Ｋ）を得ることができる。なお、前記式２３における脱硫反応の活性化エネルギーは以下のように求めることができる。
【0129】
（脱硫反応の活性化エネルギーの求め方）
脱硫反応の活性化エネルギーは、前記式６で表されるアレニウスの式に基づいて、本分野で公知の方法により求めることができる。以下、一例を説明する。
脱硫反応の活性化エネルギーは、前記式６で表されるアレニウスの式に基づいて、本分野で公知の方法により求めることができる。以下、一例を説明する。
【0130】
はじめに、重質コーカー分解軽油を含む原料油の脱硫反応の反応次数を求める。反応温度、水素分圧、水素／原料油比、原料油中の硫黄濃度は一定の条件とし、ＬＨＳＶ（ｘ）の条件で、反応を行い、生成油中の硫黄濃度を測定する。下式２４で表される脱硫反応速度式のＳＦに原料油中の硫黄濃度を、ＳＰに得られた生成油中の硫黄濃度を、ＬＨＳＶにＬＨＳＶ（ｘ）を代入する。得られた左辺の結果を縦軸に、１／ＬＨＳＶを横軸にプロットする。なお、この場合、縦軸はｎの関数となる。
【0131】
【数14】
【0132】
ｘ種類のＬＨＳＶ（ｘ）について同様の反応を行い、ｘ個の上記プロットを得る。得られたプロットを基に原点を通るように回帰直線を引き、ｙ＝ｃｘで表される直線を得る。ｙは（１／ｎ－１）（（１／ＳＰ
ｎ－１）－（１／ＳＦ
ｎ－１））であり、ｘは１／ＬＨＳＶであり、ｃはｋとなる。ｅｘｃｅｌ等で相関関数を求め、相関係数が最も１に近づくｎを求める。得られたｎが反応次数となる。なお、ｎは小数点第１位のオーダーで求めればよい。
【0133】
上記ｘは３以上の整数である。ｘの数が多い方が、より精度の高いｎを得ることができる。一方、ｘの数が多すぎると、ｎを得るために時間が掛かり効率的ではない。本実施形態においては、ｘは３～２０であることが好ましく、３～１０であることがより好ましい。
【0134】
ｎを求める上での反応温度、水素分圧、水素／原料油比、原料油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのような反応温度としては、例えば３３０～４１０℃であり、水素分圧は１４～１７ＭＰａであり、水素／原料油比は５００～１３００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度は３～５質量％である。
ｘ種類のＬＨＳＶ（ｘ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＬＨＳＶ（ｘ）としては、０．５～２．０時間－１である。
そのような反応温度としては、例えば３３０～４１０℃であり、水素分圧は１４～１７ＭＰａであり、水素／原料油比は５００～１３００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度は３～５質量％である。
ｘ種類のＬＨＳＶ（ｘ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＬＨＳＶ（ｘ）としては、０．５～２．０時間－１である。
【0135】
前記式６で表されるアレニウスの式中の活性化エネルギーＥを脱硫の活性化エネルギーＥａとし、両辺の自然対数を取ると、下式２５で表される式となる。
【数15】
【0136】
水素分圧、水素／原料油比、ＬＨＳＶ、原料油中の硫黄濃度は一定の条件とし、反応温度Ｔ（ｙ）の条件で反応を行い、生成油中の硫黄濃度を測定する。前記式２４で表される脱硫反応速度式のＳＦに原料油中の硫黄濃度を、ＳＰに得られた生成油中の硫黄濃度を、ＬＨＳＶにＬＨＳＶを、上記得られたｎを代入して反応速度定数ｋを求める。得られた反応速度定数を前記式２５に代入し、得られた左辺の結果（ｌｎｋ）を縦軸に、１／Ｔ（１／Ｔ（ｙ））を横軸にプロットする。
【0137】
ｙ種類の反応温度Ｔ（ｙ）について同様の反応を行い、ｙ個の上記プロットを得る。得られたプロットより、回帰直線を引き、その傾きを求める。この傾きはＥａ／Ｒであるため、傾きからＲを除すことにより脱硫の活性化エネルギーＥａを求めることができる。
【0138】
上記ｙは３以上の整数である。ｙの数が多い方が、より精度の高いＥａを得ることができる。一方、ｙの数が多すぎると、Ｅａを得るために時間が掛かり効率的ではない。本実施形態においては、ｙは３～２０であることが好ましく、３～１０であることがより好ましい。
【0139】
Ｅａを求める上での水素分圧、水素／原料油比、ＬＨＳＶ、原料油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのような水素分圧としては、例えば１４～１７ＭＰａであり、水素／原料油比は５００～１３００［Ｎｍ３/ｋＬ］であり、ＬＨＳＶは０．５～２．０時間－１であり、原料油中の硫黄濃度は３～５質量％である。
ｙ種類の反応温度Ｔ（ｙ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＴ（ｙ）としては、３３０～４１０℃である。
そのような水素分圧としては、例えば１４～１７ＭＰａであり、水素／原料油比は５００～１３００［Ｎｍ３/ｋＬ］であり、ＬＨＳＶは０．５～２．０時間－１であり、原料油中の硫黄濃度は３～５質量％である。
ｙ種類の反応温度Ｔ（ｙ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＴ（ｙ）としては、３３０～４１０℃である。
【0140】
このように得られた各パラメータを前記式２３に代入することにより、所定の反応条件を達成するために必要な反応温度Ｔｔを求めることができる。
【0141】
実測の反応温度Ｔｏｂｓに対する、本実施形態の情報処理方法により求められる反応温度Ｔｔの割合であるＴｔ／Ｔｏｂｓは、℃換算で、０．９７～１．０３であることが好ましく、０．９８５～１．０１５であることがより好ましい。Ｔｔ／Ｔｏｂｓが前記範囲内である場合、精度よく反応温度を推定できていると判断することができる。
【0142】
≪情報出力ステップ≫
このようにして得られた前記反応温度を示す情報を出力する情報出力ステップ（図１のＳ４）をさらに有してもよい。例えば、Ｓ４は出力部１４によって実行される。
このようにして得られた前記反応温度を示す情報を出力する情報出力ステップ（図１のＳ４）をさらに有してもよい。例えば、Ｓ４は出力部１４によって実行される。
【0143】
≪重質コーカー分解軽油を含む原料油の水素化処理反応≫
重質コーカー分解軽油を含む原料油の水素化処理反応について概要を説明する。
重質コーカー分解軽油は、コーカー分解油を蒸留して得られる、沸点範囲が常圧換算で３５０～５２０℃程度の留分である。重質コーカー分解軽油の密度は０．９６～１．００ｇ／ｍＬである。原料油中の重質コーカー分解軽油の含有量としては、例えば５０体積％以上であってもよく、７０体積％以上であってもよい。
なお、原料油に含まれる重質コーカー分解軽油以外の油種としては、流動接触分解装置から得られる重質サイクル油が例として挙げられる。原料油中の重質サイクル油の含有量としては、例えば０～５０体積％が例示される。
重質コーカー分解軽油を含む原料油の水素化処理反応について概要を説明する。
重質コーカー分解軽油は、コーカー分解油を蒸留して得られる、沸点範囲が常圧換算で３５０～５２０℃程度の留分である。重質コーカー分解軽油の密度は０．９６～１．００ｇ／ｍＬである。原料油中の重質コーカー分解軽油の含有量としては、例えば５０体積％以上であってもよく、７０体積％以上であってもよい。
なお、原料油に含まれる重質コーカー分解軽油以外の油種としては、流動接触分解装置から得られる重質サイクル油が例として挙げられる。原料油中の重質サイクル油の含有量としては、例えば０～５０体積％が例示される。
【0144】
重質コーカー分解軽油を含む原料油の水素化処理反応は、水素存在下、重質コーカー分解軽油を含む原料油と水素化処理触媒とを接触処理することにより行うことができる。
水素化処理触媒は、特に限定されるものではなく、本分野において公知の水素化処理触媒を使用することができる。触媒の担体として、種々のものが使用でき、例えばシリカ、アルミナ、ボリア、マグネシア、チタニア、シリカ－アルミナ、シリカ－マグネシア、シリカ－ジルコニア、シリカ－トリア、シリカ－ベリリア、シリカ－チタニア、シリカ－ボリア、アルミナ－ジルコニア、アルミナ－チタニア、アルミナ－ボリア、アルミナ－クロミア、チタニア－ジルコニア、シリカ－アルミナ－トリア、シリカ－アルミナ－ジルコニア、シリカ－アルミナ－マグネシア、シリカ－マグネシア－ジルコニアなど、又はこれらの２種以上の混合物が挙げられる。これらの無機酸化物のうち、好ましいものとしては、アルミナ、シリカ－アルミナ、アルミナ－チタニア、アルミナ－ボリア、アルミナ－ジルコニアが挙げられ、特に好ましくは、アルミナが挙げられ、アルミナの中でもγアルミナが特に好ましい。これらの無機酸化物は、１種単体で用いてもよいし、２種以上を組み合わせて用いてもよい。
水素化処理触媒は、特に限定されるものではなく、本分野において公知の水素化処理触媒を使用することができる。触媒の担体として、種々のものが使用でき、例えばシリカ、アルミナ、ボリア、マグネシア、チタニア、シリカ－アルミナ、シリカ－マグネシア、シリカ－ジルコニア、シリカ－トリア、シリカ－ベリリア、シリカ－チタニア、シリカ－ボリア、アルミナ－ジルコニア、アルミナ－チタニア、アルミナ－ボリア、アルミナ－クロミア、チタニア－ジルコニア、シリカ－アルミナ－トリア、シリカ－アルミナ－ジルコニア、シリカ－アルミナ－マグネシア、シリカ－マグネシア－ジルコニアなど、又はこれらの２種以上の混合物が挙げられる。これらの無機酸化物のうち、好ましいものとしては、アルミナ、シリカ－アルミナ、アルミナ－チタニア、アルミナ－ボリア、アルミナ－ジルコニアが挙げられ、特に好ましくは、アルミナが挙げられ、アルミナの中でもγアルミナが特に好ましい。これらの無機酸化物は、１種単体で用いてもよいし、２種以上を組み合わせて用いてもよい。
【0145】
前記担体に活性成分として含有させる金属は、周期表第６族金属及び第８～１０族金属の中から選ばれる少なくとも１種類以上の金属であり、好ましくはモリブデン、タングステン、コバルト及びニッケルの金属である。これらの金属は、金属状態又は金属酸化物、金属硫化物の何れの形態でも有効であり、また、イオン交換などにより金属が触媒担体と結合した形態で存在してもよい。この金属成分の含有量は、通常、触媒基準かつ酸化物換算で、約１～２５質量％の範囲内である。金属含有量が１質量％より少ないと、活性点として働く金属の絶対量が少ないために、脱硫活性を始めとする水素化処理活性（以下、簡単に水素化処理活性と言う）が発現せず、逆に担持される金属の含有量が２５質量％より多すぎると、金属の凝集が起こり活性点の数が減少し、その結果、水素化処理活性が却って低下するからである。更に、必要に応じて、元素周期律表第６族金属及び第８族金属からなる活性金属に加えて、リン、ホウ素、亜鉛、ジルコニア等を含ませることができる。本発明方法を適用するに当たり、触媒層の形態には制約はなく、例えば固定床、移動床、流動床等の触媒層の反応器に適用できる。
【0146】
重質コーカー分解軽油を含む原料油の水素化処理反応における条件としては、一般に反応温度が３３０～４１０℃、好ましくは３６０～４００℃であり、水素分圧が１０～２０ＭＰａ、好ましくは１４～１７ＭＰａであり、ＬＨＳＶが０．１～２．０ｈｒ－１、好ましくは０．５～１．２ｈｒ－１であり、水素／原料油比が１７０～１４００［Ｎｍ３/ｋＬ］、好ましくは５００～１３００［Ｎｍ３/ｋＬ］である。
【0147】
重質コーカー分解軽油を含む原料油中の硫黄濃度は、通常３～５質量％である。また、生成油中の硫黄濃度は、通常０．１質量％以下である。重質コーカー分解軽油を含む原料油中のシリコン濃度は、通常０．３重量ｐｐｍ～７重量ｐｐｍである。また、生成油中のシリコン濃度は、通常０．３重量ｐｐｍ未満である。
【0148】
≪反応温度算出装置≫
本実施形態の反応温度算出装置は、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む。本実施形態の反応温度算出装置は、算出された反応温度を示す情報を出力する出力部を有していてもよい。
本実施形態の反応温度算出装置は、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む。本実施形態の反応温度算出装置は、算出された反応温度を示す情報を出力する出力部を有していてもよい。
【0149】
本実施形態の反応温度算出装置１は、パーソナルコンピューターやサーバー装置や専用装置等の情報処理装置を用いて構成される。反応温度算出装置１は、１台又は複数台の情報処理装置を用いて構成されてもよい。例えば、反応温度算出装置１は、クラスタマシンとして構築されてもよいし、クラウドとして構築されてもよいし、どのような態様で構築されてもよい。反応温度算出装置１は、例えば、図４に示されるように取得部１１と、取得部からの情報を処理する計算機本体１２とを有する。反応温度算出装置１は、計算機本体１２において処理された情報を外部に出力する出力部１４を有していてもよい。これらの構成要素は、例えば、ＣＰＵ（Central Processing Unit）等のハードウェアプロセッサがプログラム（ソフトウェア）を実行することにより実現される。また、これらの構成要素のうち一部または全部は、ＬＳＩ（Large Scale Integrated circuit）やＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field-Programmable Gate Array）、ＧＰＵ（Graphics Processing Unit）等のハードウェア（回路部；circuitryを含む）によって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。プログラムは、予めＨＤＤ（Hard Disk Drive）やフラッシュメモリ等の記憶装置（非一過性の記憶媒体を備える記憶装置）に格納されていてもよいし、ＤＶＤやＣＤ－ＲＯＭ等の着脱可能な記憶媒体（非一過性の記憶媒体）に格納されており、記憶媒体がドライブ装置に装着されることで記憶装置にインストールされてもよい。記憶装置は、例えば、ＨＤＤ、フラッシュメモリ、ＥＥＰＲＯＭ（Electrically Erasable Programmable Read Only Memory）、ＲＯＭ（Read Only Memory）、またはＲＡＭ（Random Access Memory）等により構成される。
【0150】
取得部１１は、反応のオペレーターによって所定の情報が入力され、この入力により取得した情報を計算機本体１２に送信するものである。本実施形態の取得部１１が取得する情報は、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報である。反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報としては上述した通りである。例えば、取得部１１は、上述した情報取得ステップを実行する。取得部１１は、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得すればよく、その取得方法には特に限定されない。
【0151】
本実施形態では、取得部１１は単一のキーボードによって構成されている。取得部１１の具体的構成は、限定されず、本実施形態ではキーボードであるが、タッチパネル等であってもよい。なお、各種情報を取得する取得部が別々に構成され、それぞれが独立して計算機本体１２に接続されてもよい。また、取得部１１は、反応器の制御等に用いられる計算機等から有線又は無線により前記の各情報を直接取得するように構成されてもよい。
【0152】
計算機本体１２は、例えば、種々の情報を処理可能な、いわゆるコンピュータである。計算機本体１２は、演算部１３を備える。例えば、この計算機本体１２には、所定のプログラムが組み込まれ、このプログラムの実行によって機能的に演算部１３が構成される。具体的には、この演算部１３において、前記取得部１１で取得した反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する。劣化関数としては上述した通りである。上述したように反応温度は、例えば劣化速度式から求めることができる。例えば、演算部１３は、上述した劣化度算出ステップ、及び反応温度算出ステップを実行する。演算部１３は、例えば、ＣＰＵ（Central Processing Unit)又はＭＰＵ（Micro Processing Unit）などのプロセッサ及び不揮発性又は揮発性の半導体メモリ（例えば、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、フラッシュメモリ、ＥＰＲＯＭ（Erasable Programmable Read Only Memory）、ＥＥＰＲＯＭ（Electrically Erasable Programmable Read Only Memory））を備えてもよい。例えば、演算部１３は、ＭＣＵなどのマイクロコントローラであってもよい。
【0153】
演算部１３は、上記のように求めた前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を示す情報を出力部１４に出力してもよい。
【0154】
出力部１４は、計算機本体１２（詳しくは、演算部１３）が出力した計算結果（反応温度）を受信し、受信した計算結果を外部に出力するものである。本実施形態の出力部１４は、ＣＲＴディスプレイや液晶ディスプレイ、ＰＤＰ等の表示部によって構成されているが、これに限定されず、プリンタ等の印刷部や、他の装置（例えば、重質コーカー分解軽油を含む原料油の水素化処理反応の制御等に用いられる計算機等）等へ出力するように構成されてもよい。また、出力部１４は、これらを組み合わせたものでもよい。例えば、出力部１４は、上述した情報出力ステップを実行する。
【0155】
また、本実施形態においては、コンピュータを反応温度算出装置として機能させるための反応温度算出プログラム及び当該プログラムを記憶したコンピュータの非一時的可読記録媒体が提供される。コンピュータの非一時的可読記録媒体としては、例えば、磁気テープ（デジタルデータストレージ（ＤＳＳ）など）、磁気ディスク（ハードディスクドライブ（ＨＤＤ）、フレキシブルディスク（ＦＤ）など）、光ディスク（コンパクトディスク（ＣＤ）、デジタルバーサタイルディスク（ＤＶＤ）、ブルーレイディスク（ＢＤ）など）、光磁気ディスク（ＭＯ）、フラッシュメモリ（ＳＳＤ（Solid State Drive）、メモリーカード、ＵＳＢメモリなど）が挙げられる。
【0156】
＜情報処理方法及び反応温度算出装置の活用方法＞
本実施形態の情報処理方法及び反応温度算出装置によると、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる。本実施形態の情報処理方法及び反応温度算出装置によると、反応温度の推定値の経時的なプロットを得ることができる。前記プロットと、機器使用最高温度等との関係から、以下のような活用法が考えられる。
本実施形態の情報処理方法及び反応温度算出装置によると、重質コーカー分解軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる。本実施形態の情報処理方法及び反応温度算出装置によると、反応温度の推定値の経時的なプロットを得ることができる。前記プロットと、機器使用最高温度等との関係から、以下のような活用法が考えられる。
【0157】
第１の活用方法としては、所定の反応条件における触媒を交換するのに必要な運転時間（日数）を推定することである。すなわち、前記プロットにより、機器使用最高温度となる温度となる時間を推定することができ、当該時間から触媒を交換するのに必要な運転時間（日数）を推定することができる。また、推定された反応温度が機器設定温度以下でない場合、図２のＳ４－１～Ｓ４－３で示される情報処理方法が行われてもよい。図２のＳ４－１～Ｓ４－３で示される情報処理方法は、目的とする反応温度に関する情報を取得する情報取得ステップ（図２のＳ４－１）と、取得した目的とする反応温度に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップ（図２のＳ４－２）と、劣化度に基づき、目的とする反応温度を満たすために必要な、原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を算出する反応条件算出ステップ（図２のＳ４－３）と、を含む。このようにして得られた前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報を出力する情報出力ステップ（図２のＳ４－４）をさらに有してもよい。具体的には、例えば、コーク劣化関数１及びシリコン劣化関数を使用する場合、前記式２３のＴｔに目的とする反応温度を代入し、Φを求める。得られたΦを前記式２に代入し、ΦＣ、ΦSiを求める。得られたΦＣを前記式３に代入し、Ｄを求める。得られたＤを前記式４に代入して、前記式４の等式が成り立つような、ＳＰ、ＳＦ、ＬＨＳＶ、Ｐの組み合わせを求めればよい。また、上記組み合わせにおいて、例えばＳＰ、ＳＦを設定値として、Ｐ及びＬＨＳＶを求めてもよい。また、コーク劣化関数としては、コーク劣化関数１、コーク劣化関数１－１、コーク劣化関数２、コーク劣化関数２－１のいずれを使用してもよい。また、上記はΦSiを変更しない前提で説明を行ったが、ΦＣを変更せず、ΦＳｉを変更してもよい。具体的には、シリコン劣化関数において、Ｓｉｆ、Ｓｉｐ、ＷＨＳＶの条件を変更してもよい。なお、以上の各ステップは、例えば、本実施形態の反応温度算出装置１によって実行される。例えば、Ｓ４－１は取得部１１によって実行され、Ｓ４－２及びＳ４－３は計算機本体１２における演算部１３によって実行され、Ｓ４－４は出力部１４によって実行される。なお、目的温度は予め計算機本体１２における演算部１３に記憶させておいてもよい。
【0158】
第２の活用方法としては、所定の運転時間を達成するための反応条件（処理量（ＬＨＳＶ）等）を推定することである。このような情報処理方法は、図３のＳ１Ａ～３Ａで表される。図３のステップ１Ａ～３Ａで示される情報処理方法は、所定の運転時間における目的とする反応温度に関する情報を取得する情報取得ステップ（図３のＳ１Ａ）と、取得した目的とする反応温度に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップ（図３のＳ２Ａ）と、劣化度に基づき、目的とする反応温度を満たすために必要な、原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を算出する反応条件算出ステップ（図３のＳ３Ａ）と、を含む。このようにして得られた前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報を出力する情報出力ステップ（図３のＳ４Ａ）をさらに有してもよい。具体的には、例えば、コーク劣化関数１及びシリコン劣化関数を使用する場合、前記式２３のＴｔに目的とする反応温度を代入し、Φを求める。得られたΦを前記式２に代入し、ΦＣ、ΦＳｉを求める。得られたΦＣを前記式３に代入し、Ｄを求める。得られたＤを前記式４に代入して、前記式４の等式が成り立つような、ＳＰ、ＳＦ、ＬＨＳＶ、Ｐの組み合わせを求めればよい。また、上記組み合わせにおいて、例えばＳＰ、ＳＦを設定値として、Ｐ及びＬＨＳＶを求めてもよい。また、コーク劣化関数としては、コーク劣化関数１、コーク劣化関数１－１、コーク劣化関数２、コーク劣化関数２－１のいずれを使用してもよい。また、上記はΦSiを変更しない前提で説明を行ったが、ΦＣを変更せず、ΦＳｉを変更してもよい。具体的には、シリコン劣化関数において、Ｓｉｆ、Ｓｉｐ、ＷＨＳＶの条件を変更してもよい。なお、以上の各ステップは、例えば、本実施形態の反応温度算出装置１によって実行される。例えば、Ｓ１Ａは取得部１１によって実行され、Ｓ２Ａ及びＳ３Ａは計算機本体１２における演算部１３によって実行され、Ｓ４Ａは出力部１４によって実行される。なお、目的温度は予め計算機本体１２における演算部１３に記憶させておいてもよい。
【実施例0159】
以下、実施例により本発明をさらに具体的に説明するが、本発明は以下の実施例に限定されるものではない。
【0160】
［実施例］
ベンチスケールで、重質コーカー分解軽油を１００体積％含む原料油を水素化処理触媒に接触処理させることにより水素化処理反応を行った。得られた結果を基に、劣化関数を算出した。なお、本実施例では、前記式１０で表されるコーク劣化関数１－１を使用した。前記式１１及び前記式２０におけるパラメータを上述の方法により求めた所、α’＝０．０１３５、ＰＢ＝１５．９（ＭＰａ）、ａ＝０．６、Ｅｃ＝１５６（ｋＪ／ｍｏｌ）、ＴＢ＝６４３．１５（Ｋ）、ＴＳＯＲ=６２８．９５（Ｋ）、ＧＢ=５５０（Ｎｍ３/ｋＬ）、ｂ＝０．２、ｎ＝１．５、η＝０．６、ＳｉＯＣ∞＝９．０（重量％）であった。
ベンチスケールで、重質コーカー分解軽油を１００体積％含む原料油を水素化処理触媒に接触処理させることにより水素化処理反応を行った。得られた結果を基に、劣化関数を算出した。なお、本実施例では、前記式１０で表されるコーク劣化関数１－１を使用した。前記式１１及び前記式２０におけるパラメータを上述の方法により求めた所、α’＝０．０１３５、ＰＢ＝１５．９（ＭＰａ）、ａ＝０．６、Ｅｃ＝１５６（ｋＪ／ｍｏｌ）、ＴＢ＝６４３．１５（Ｋ）、ＴＳＯＲ=６２８．９５（Ｋ）、ＧＢ=５５０（Ｎｍ３/ｋＬ）、ｂ＝０．２、ｎ＝１．５、η＝０．６、ＳｉＯＣ∞＝９．０（重量％）であった。
【0161】
ベンチスケールで行った水素化処理反応と同様の原料油、水素化処理触媒を用いて、実機で反応を行った。前記式１１および式２１における運転条件として、ＳＦ＝３．５４（質量％）、Ｓｉｆ＝１．５（重量ｐｐｍ）、Ｐ＝１５．６（ＭＰａ）、Ｇ＝６１４（Ｎｍ３/ｋＬ）、ＬＨＳＶ＝０．７８７（時間－１）、ＳＰ＝０．１（質量％）、Ｓｉｐ＝０（重量ｐｐｍ）とした。
また、事前に上述の方法で、脱硫の活性化エネルギーを求めた所、Ｅａ＝１２０（ｋＪ／ｍｏｌ）であった。
また、事前に上述の方法で、脱硫の活性化エネルギーを求めた所、Ｅａ＝１２０（ｋＪ／ｍｏｌ）であった。
【0162】
これらの基本劣化パラメータ及び反応条件を前記式１１に代入し、任意の反応ｔ日経過時のＤ’を得た。得られたＤ’及び反応経過日数ｔを前記式１０に代入しΦＣを得た。また、Ｓｉｆ、Ｓｉｐ、反応経過日数ｔ、ＬＨＳＶと原油通油量、触媒充填量から算出したＷＨＳＶを前記式２１に代入してＳｉＯＣを得、得られたＳｉＯＣ、上記のη及びＳｉＯＣ∞を前記式２０に代入しΦＳｉを得た。得られたΦＣとΦＳｉから前記式２によりΦを求めた。このようにして得られたΦ、Ｅａ、ＴＳＯＲを前記式２３に代入して、任意の反応ｔ日経過時の要求温度Ｔｔを求めた。表１にｔ＝５１３日、１００８日、２００２日、２４７４日の反応温度の実測値（但し、ＳＰが０．１（質量％）になるように補正した後の反応温度）、及び上記方法で求めた要求温度Ｔｔ、及び要求温度Ｔｔ／反応温度の実測値の割合を示す。
【0163】
【表1】
【0164】
表１に示された通り、本発明によって求められた要求温度Ｔｔは、反応温度の実測値とほぼ同等となることがわかった。
【符号の説明】
【0165】
１・・・反応温度算出装置
１１・・・取得部
１２・・・計算機本体
１３・・・演算部
１４・・・出力部
１１・・・取得部
１２・・・計算機本体
１３・・・演算部
１４・・・出力部
【図1】
【図2】
【図3】
【図4】