特開2022-154058情報処理方法、反応温度算出装置、反応温度算出プログラム、及びコンピュータの非一時的可読記録媒体
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2022154058
(43)【公開日】2022-10-13
(54)【発明の名称】情報処理方法、反応温度算出装置、反応温度算出プログラム、及びコンピュータの非一時的可読記録媒体
(51)【国際特許分類】
C10G 45/02 20060101AFI20221005BHJP
G06Q 50/04 20120101ALI20221005BHJP
【ＦＩ】
C10G45/02
G06Q50/04
【審査請求】未請求
【請求項の数】10
【出願形態】ＯＬ
(21)【出願番号】P 2021056915
(22)【出願日】2021-03-30
(71)【出願人】
【識別番号】000105567
【氏名又は名称】コスモ石油株式会社
(74)【代理人】
【識別番号】100141139
【弁理士】
【氏名又は名称】及川 周
(74)【代理人】
【識別番号】100126882
【弁理士】
【氏名又は名称】五十嵐 光永
(74)【代理人】
【識別番号】100162868
【弁理士】
【氏名又は名称】伊藤 英輔
(74)【代理人】
【識別番号】100147267
【弁理士】
【氏名又は名称】大槻 真紀子
(74)【代理人】
【識別番号】100209347
【弁理士】
【氏名又は名称】内田 洋平
(72)【発明者】
【氏名】大崎 貴之
(72)【発明者】
【氏名】出井 一夫
【テーマコード（参考）】
4H129
5L049
【Ｆターム（参考）】
4H129AA02
4H129CA09
4H129DA15
4H129DA16
4H129DA18
4H129KA06
4H129KA07
4H129KA08
4H129KB03
4H129LA09
4H129LA14
4H129NA01
4H129NA37
4H129NA45
5L049CC04
(57)【要約】
【課題】減圧蒸留軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる情報処理方法、及び前記反応温度を推定することができる反応温度算出装置、コンピュータを、前記反応温度算出装置として機能させるための反応温度算出プログラム、及び前記プログラムを記憶したコンピュータの非一時的可読記録媒体の提供。
【解決手段】反応を開始してから所定時間経過した際の減圧蒸留軽油を含む原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、取得した前記情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
【選択図】なし
【解決手段】反応を開始してから所定時間経過した際の減圧蒸留軽油を含む原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、取得した前記情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
【選択図】なし
【特許請求の範囲】
【請求項1】
常圧蒸留残渣油を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する情報取得ステップと、
取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、
前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップと、
前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する反応温度算出ステップと、を含む、情報処理方法。
【請求項2】
前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数であり、触媒の易失活活性種の劣化に関する易失活活性種劣化関数と、触媒の難失活活性種の劣化に関する難失活活性種劣化関数から構成される、請求項１に記載の情報処理方法。
【請求項3】
前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、請求項１又は２に記載の情報処理方法。
【請求項4】
前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、請求項１～３のいずれか一項に記載の情報処理方法。
【請求項5】
常圧蒸留残渣油を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、
前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む、反応温度算出装置。
前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む、反応温度算出装置。
【請求項6】
前記劣化関数は、コークの堆積による触媒の劣化に関するコーク劣化関数であり、触媒の易失活活性種の劣化に関する易失活活性種劣化関数と、触媒の難失活活性種の劣化に関する難失活活性種劣化関数から構成される、請求項５に記載の反応温度算出装置。
【請求項7】
前記原料油に関する情報は、原料油中の硫黄濃度に関する情報を含み、前記生成油に関する情報は、生成油中の硫黄濃度に関する情報を含む、請求項５又は６に記載の反応温度算出装置。
【請求項8】
前記運転条件に関する情報は、水素分圧に関する情報、触媒充填量に関する情報、及び原料油の供給量に関する情報を含む、請求項５～７のいずれか一項に記載の反応温度算出装置。
【請求項9】
コンピュータを、請求項５～８のいずれか一項に記載の反応温度算出装置として機能させるための反応温度算出プログラム。
【請求項10】
請求項９に記載のプログラムを記憶したコンピュータの非一時的可読記録媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、情報処理方法、反応温度算出装置、反応温度算出プログラム、及びコンピュータの非一時的可読記録媒体に関する。
【背景技術】
【0002】
原油を常圧蒸留して得られる常圧蒸留残渣油には、高濃度の硫黄が含有している。常圧蒸留残渣油を水素化処理することにより硫黄含有量を下げる方式は、直接脱硫として知られている。一方、常圧蒸留残渣油を減圧蒸留して、減圧蒸留軽油と減圧蒸留残渣油とに分け、減圧蒸留軽油を水素化処理することにより硫黄含有量を下げ、水素化処理された減圧蒸留軽油と減圧蒸留残渣とを混合する方式は、間接脱硫として知られている。
【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】
【数1】
【0028】
本明細書において、「要求温度」とは、所定の反応条件を達成するために必要な反応温度を意味する。すなわち、０日目の要求温度とは、反応開始時において、上記ＳＦ、ＳＰ、ＬＨＳＶ、及びＰの反応条件を達成するために必要な反応温度を意味する。
【0029】
本明細書において、「基準水素分圧」とは、実際の反応条件の標準的な圧力を意味する。後述する水素分圧係数ａを決定する際に用いる反応圧力の平均値として求められる。
【0030】
本明細書において、「基準反応温度」とは、実際に運転される可能性のある標準的な運転条件の振れ幅で得られるＴＳＯＲの平均値を意味する。
【0031】
前記式３において、ＴＳＯＲは触媒の初期活性を表し、αは触媒の劣化速度を表す。すなわち、ＴＳＯＲとαの数値が大きいほど、触媒の劣化が大きいことを意味し、この劣化挙動はＤの値に反映されることになる。
【0032】
前記式３において、（１／ＳＰ
ｎ－１－１／ＳＦ
ｎ－１）ＬＨＳＶで表される項は、脱硫反応速度定数であり、上述した通り、ＳＰ、ＳＦ、ＬＨＳＶを設定値とし、一定の条件で運転する場合、定数となる。前記式３において、（ＰＢ／Ｐ）ａで表される項は、水素分圧依存性を示す項であり、上述した通り、Ｐを設定値とし、一定の条件で運転する場合、定数となる。前記式３においてｅｘｐ［Ｅｃ／Ｒ（１／ＴＢ－１／ＴＳＯＲ）］で表される項は、温度依存性を示す項であり、定数となる。
【0033】
前記式３中、ＳＦ、ＳＰ、ＬＨＳＶ、Ｐは、上述した情報取得ステップで取得した原料油に関する情報、生成油に関する情報、及び運転条件に関する情報に基づき代入される値である。なお、ＬＨＳＶは、原料油の供給量（体積／ｈ）を触媒充填量（体積）で除すことにより求めることができる。
【0034】
上述した通り、ＳＦ、ＬＨＳＶ、Ｐは制御することが可能なパラメータである。また、ＳＰは目的とする生成油の硫黄濃度である。すなわち、前記式３により、上記ＳＦ、ＳＰ、ＬＨＳＶ、及びＰの反応条件における、触媒の活性種の劣化係数を算出することができる。減圧蒸留軽油を含む原料油の水素化処理反応の反応次数であるｎの求め方は後述する。
【0035】
（基本劣化パラメータの求め方）
前記式３中、α、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは定数である。以下、これらのパラメータを総称して「基本劣化パラメータ１」という。基本劣化パラメータ１は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。
以下、基本劣化パラメータＰＢ、ａ、Ｅｃ、ＴＢの求め方について例を示す。上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（ＰＢ、ａ、Ｅｃ、ＴＢの求め方）である。また、α、ＴＳＯＲの求め方については２例を示すが、本発明はこれらに限定されるものではない。１例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（α及びＴＳＯＲの求め方１）であり、２例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される反応温度プロファイルから求める方法（α及びＴＳＯＲの求め方２）である。
前記式３中、α、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは定数である。以下、これらのパラメータを総称して「基本劣化パラメータ１」という。基本劣化パラメータ１は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。
以下、基本劣化パラメータＰＢ、ａ、Ｅｃ、ＴＢの求め方について例を示す。上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（ＰＢ、ａ、Ｅｃ、ＴＢの求め方）である。また、α、ＴＳＯＲの求め方については２例を示すが、本発明はこれらに限定されるものではない。１例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動（反応速度定数の変化）から求める方法（α及びＴＳＯＲの求め方１）であり、２例目は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される反応温度プロファイルから求める方法（α及びＴＳＯＲの求め方２）である。
【0036】
（ＰＢ、ａ、Ｅｃ、ＴＢの求め方）
本実施形態の基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。この触媒の劣化挙動（劣化度）は、前記式１と同じ考え方となり、下式４で表すことができる。
Φ’＝ｋｔ’／ｋ０’ 式４
前記式４中、ｋ０’は反応０日経過時（すなわち、反応開始時）の触媒の反応速度定数であり、ｋｔ’は任意の反応ｔ日経過時の触媒の反応速度定数である。なお、ｋ０’、ｋｔ’は後述の温度ＴＳＯＲ’における反応速度定数である。
本実施形態の基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。この触媒の劣化挙動（劣化度）は、前記式１と同じ考え方となり、下式４で表すことができる。
Φ’＝ｋｔ’／ｋ０’ 式４
前記式４中、ｋ０’は反応０日経過時（すなわち、反応開始時）の触媒の反応速度定数であり、ｋｔ’は任意の反応ｔ日経過時の触媒の反応速度定数である。なお、ｋ０’、ｋｔ’は後述の温度ＴＳＯＲ’における反応速度定数である。
【0037】
前記式４は、前記式１と同様、反応速度定数に基づいた劣化関数である。反応速度定数は下式５のアレニウスの式で表される。
【0038】
【数2】
【0039】
反応開始時温度ＴＳＯＲ’における反応速度定数ｋ０’は、反応ｔ日経過時にはｋｔ’に低下する。反応ｔ日経過時に反応速度定数ｋ０’に相当する活性を得るためには、反応温度をＴｔ’としなければならないものとすると、前記式４及び前記式５より、下式６が導かれる。なお、本実施形態の反応は、減圧蒸留軽油を含む原料油の水素化処理反応のため、活性化エネルギーＥを脱硫の活性化エネルギーＥａ（ｋＪ／ｍｏｌ）としている。
【0040】
【数3】
【0041】
（Ｅｃ及びＴＢの求め方）
ＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度を一定の条件とし、生成油中の硫黄濃度を一定の値ＳＰｎになるように、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度をＳＰｎとするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｎｘ＋ｂｎ（０＜ａｎである。）で表される直線が得られる。ａｎ及びｂｎは触媒の劣化挙動を反映した値である。この直線におけるｂｎが前記式６におけるＴＳＯＲ’となる。前記式６において、ＴＳＯＲ’にｂｎを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。なお、脱硫の活性化エネルギーＥａは、後述の方法により求められた値を使用することができる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｎ ’ｘ（|－ａｎ ’|＝ａｎ ’である。）で表される直線が得られる。ａｎ ’は、触媒の劣化速度を表している。
ＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度を一定の条件とし、生成油中の硫黄濃度を一定の値ＳＰｎになるように、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度をＳＰｎとするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｎｘ＋ｂｎ（０＜ａｎである。）で表される直線が得られる。ａｎ及びｂｎは触媒の劣化挙動を反映した値である。この直線におけるｂｎが前記式６におけるＴＳＯＲ’となる。前記式６において、ＴＳＯＲ’にｂｎを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。なお、脱硫の活性化エネルギーＥａは、後述の方法により求められた値を使用することができる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｎ ’ｘ（|－ａｎ ’|＝ａｎ ’である。）で表される直線が得られる。ａｎ ’は、触媒の劣化速度を表している。
【0042】
ｎ種類の硫黄濃度ＳＰｎについて同様の反応を行い、同様にして、ｎ個のａｎ、ｂｎを求め、上記と同様の方法により、ｎ個のａｎ
’を得る。ｎは３以上の整数である。ｎの数が多い方が、より精度の高いＥｃを得ることができる。一方、ｎの数が多すぎると、Ｅｃを得るために時間が掛かり効率的ではない。本実施形態においては、ｎは３～２０であることが好ましく、３～１０であることがより好ましい。
このようにして得られたｎ個のａｎ ’及びｂｎを下式７にそれぞれ代入する。下式７は、コークの活性化エネルギー、及び基準反応温度を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
このようにして得られたｎ個のａｎ ’及びｂｎを下式７にそれぞれ代入する。下式７は、コークの活性化エネルギー、及び基準反応温度を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
【0043】
ｌｎ（ａｎ
’）＝ｌｎ（Ａ）－（Ｅｃ／Ｒｂｎ） 式７
前記式７中、Ａは頻度因子であり、Ｅｃはコーク劣化の活性化エネルギー（ｋＪ／ｍｏｌ）であり、Ｒは気体定数：０．００８３１（ｋＪ／（ｍｏｌ・Ｋ））である。
前記式７中、Ａは頻度因子であり、Ｅｃはコーク劣化の活性化エネルギー（ｋＪ／ｍｏｌ）であり、Ｒは気体定数：０．００８３１（ｋＪ／（ｍｏｌ・Ｋ））である。
【0044】
ｎ個のａｎ
’及びｂｎの組み合わせに関し、ｌｎ（ａｎ
’）を縦軸に、１／ｂｎを横軸にプロットし、回帰直線を引き、その傾きを求める。この傾きはＥｃ／Ｒであるため、傾きからＲを除すことによりコーク劣化の活性化エネルギーＥｃを求めることができる。
【0045】
また、ｎ個のｂｎを平均することによりＴＢを求めることができる。
【0046】
Ｅｃ及びＴＢを求める上でのＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのようなＬＨＳＶとしては、例えば０．５～２．５ｈ－１であり、水素分圧としては、例えば３～７ＭＰａであり、水素／原料油比としては、例えば１００～４００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度としては、例えば０．５～３．０質量％である。
ｎ種類の硫黄濃度ＳＰｎも同様に実機運転条件に即した条件とすることが好ましい。このようなＳＰｎとしては、例えば０．３質量％以下である。
反応期間としては、例えば１００～１５００日である。
そのようなＬＨＳＶとしては、例えば０．５～２．５ｈ－１であり、水素分圧としては、例えば３～７ＭＰａであり、水素／原料油比としては、例えば１００～４００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度としては、例えば０．５～３．０質量％である。
ｎ種類の硫黄濃度ＳＰｎも同様に実機運転条件に即した条件とすることが好ましい。このようなＳＰｎとしては、例えば０．３質量％以下である。
反応期間としては、例えば１００～１５００日である。
【0047】
（ＰＢ及びａの求め方）
ＬＨＳＶ、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度を一定の条件とし、水素分圧Ｐｍの条件で、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度を一定の値とするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｍｘ＋ｂｍ（０＜ａｍである。）で表される直線が得られる。この直線におけるｂｍが前記式６におけるＴＳＯＲ’となる。前記式６において、ＴＳＯＲ’にｂｍを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｍ ’ｘ（|―ａｍ ’|＝ａｍ ’である。）で表される直線が得られる。ａｍ ’は、触媒の劣化速度を表している。
ＬＨＳＶ、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度を一定の条件とし、水素分圧Ｐｍの条件で、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度を一定の値とするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｍｘ＋ｂｍ（０＜ａｍである。）で表される直線が得られる。この直線におけるｂｍが前記式６におけるＴＳＯＲ’となる。前記式６において、ＴＳＯＲ’にｂｍを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｍ ’ｘ（|―ａｍ ’|＝ａｍ ’である。）で表される直線が得られる。ａｍ ’は、触媒の劣化速度を表している。
【0048】
ｍ種類の水素分圧Ｐｍについて同様の反応を行い、同様にして、ｍ個のａｍ、ｂｍを求め、上記と同様の方法により、ｍ個のａｍ
’を得る。ｍは３以上の整数である。ｍの数が多い方が、より精度の高いａを得ることができる。一方、ｍの数が多すぎると、ａを得るために時間が掛かり効率的ではない。本実施形態においては、ｍは３～２０であることが好ましく、３～１０であることがより好ましい。
このようにして得られたｍ個のａｍ’及びＰｍを下式８にそれぞれ代入する。下式８は、水素分圧係数、及び基準水素分圧を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
このようにして得られたｍ個のａｍ’及びＰｍを下式８にそれぞれ代入する。下式８は、水素分圧係数、及び基準水素分圧を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
【0049】
ｌｎ（ａｍ
’）＝－ａｌｎ（Ｐｍ）＋Ｂ１ 式８
前記式８中、Ｂ１は０とすることができる。
前記式８中、Ｂ１は０とすることができる。
【0050】
ｍ個のａｍ’及びＰｍの組み合わせに関し、ｌｎ（ａｍ
’）を縦軸に、ｌｎ（Ｐｍ）を横軸にプロットし、回帰直線を引き、その傾きを求める。この傾きが水素分圧係数ａである。
【0051】
また、上記ｍ個の水素分圧Ｐｍを平均することによりＰＢを求めることができる。
【0052】
ａ及びＰＢを求める上でのＬＨＳＶ、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのようなＬＨＳＶとしては、例えば０．５～２．５ｈ－１であり、水素／原料油比としては、例えば１００～４００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度としては、例えば０．５～３．０質量％であり、生成油中の硫黄濃度としては、例えば０．３質量％以下である。
ｍ種類の水素分圧Ｐｍも同様に実機運転条件に即した条件とすることが好ましい。このようなＰｍとしては、例えば３～７ＭＰａである。反応期間としては、例えば５０～６００日である。
そのようなＬＨＳＶとしては、例えば０．５～２．５ｈ－１であり、水素／原料油比としては、例えば１００～４００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度としては、例えば０．５～３．０質量％であり、生成油中の硫黄濃度としては、例えば０．３質量％以下である。
ｍ種類の水素分圧Ｐｍも同様に実機運転条件に即した条件とすることが好ましい。このようなＰｍとしては、例えば３～７ＭＰａである。反応期間としては、例えば５０～６００日である。
【0053】
（α及びＴＳＯＲの求め方１）
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。想定する実機運転条件とは、（Ｅｃ及びＴＢの求め方）、（ＰＢ及びａの求め方）で説明した運転条件が例として挙げられる。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａαｘ＋ｂαで表される直線が得られる。この直線におけるａαが前記式３におけるαに相関した値であり、ｂαが前記式３におけるＴＳＯＲとなる。前記式６において、ＴＳＯＲ’にｂαを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａα ’ｘ（|－ａα ’|＝ａα ’とする。）で表される直線が得られる。ａα ’は、触媒の劣化速度を表している。
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。想定する実機運転条件とは、（Ｅｃ及びＴＢの求め方）、（ＰＢ及びａの求め方）で説明した運転条件が例として挙げられる。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａαｘ＋ｂαで表される直線が得られる。この直線におけるａαが前記式３におけるαに相関した値であり、ｂαが前記式３におけるＴＳＯＲとなる。前記式６において、ＴＳＯＲ’にｂαを代入し、Ｔｔ’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａα ’ｘ（|－ａα ’|＝ａα ’とする。）で表される直線が得られる。ａα ’は、触媒の劣化速度を表している。
【0054】
運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（すなわち、ｂα）を前記式３に代入してＤを求める。求められたＤを前記式２に代入するとΦが得られる。この場合、Φはαの関数となる。反応時間を横軸に、Φの対数を縦軸にプロットし、αを０＜αとなるように変えて回帰直線を引くと、ｙ＝－ａα
”ｘ（|－ａα
”|＝ａα
”である。）で表される複数の直線がαの値ごとに得られる。各直線におけるａα
”と、上述のａα
’が等しい値となるときのαを前記式３におけるαとすることができる。
【0055】
（α及びＴＳＯＲの求め方２）
（α及びＴＳＯＲの求め方１）と同様の反応を行い、ｙ＝ａαｘ＋ｂαで表される直線を得、ｂαを前記式３におけるＴＳＯＲとする。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（すなわち、ｂα）を前記式３に代入してＤを求める。求められたＤを前記式２に代入するとΦが得られる。この場合、Φはαの関数となる。得られたΦを前記式６のΦ’に、ＴＳＯＲ（すなわち、ｂα）を前記式６のＴＳＯＲ ’に代入して、Ｔｔ’について整理すると、Ｔｔ’はαの関数となる。実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのαを前記式３におけるαとすることができる。同様にＮ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのαを前記式３におけるαとすることが好ましい。Ｎは１０以上整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
（α及びＴＳＯＲの求め方１）と同様の反応を行い、ｙ＝ａαｘ＋ｂαで表される直線を得、ｂαを前記式３におけるＴＳＯＲとする。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（すなわち、ｂα）を前記式３に代入してＤを求める。求められたＤを前記式２に代入するとΦが得られる。この場合、Φはαの関数となる。得られたΦを前記式６のΦ’に、ＴＳＯＲ（すなわち、ｂα）を前記式６のＴＳＯＲ ’に代入して、Ｔｔ’について整理すると、Ｔｔ’はαの関数となる。実測の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）が１となるときのαを前記式３におけるαとすることができる。同様にＮ個の反応温度Ｔｏｂｓに対するＴｔ’の比（Ｔｔ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのαを前記式３におけるαとすることが好ましい。Ｎは１０以上整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
【0056】
（劣化関数１の変形例）
劣化関数１の変形例を以下に説明する。劣化関数１としては下式９で表される劣化関数１－１を使用してもよい。
劣化関数１の変形例を以下に説明する。劣化関数１としては下式９で表される劣化関数１－１を使用してもよい。
【0057】
Φ＝ｅｘｐ（－Ｄ’ｔ） 式９
前記式９中、Ｄ’は触媒の活性種の劣化係数であり、ｔは反応経過日数（日）である。
前記式９中、Ｄ’は触媒の活性種の劣化係数であり、ｔは反応経過日数（日）である。
【0058】
Ｄ’は下式１０により、求めることができる。
【0059】
【数4】
【0060】
前記式１０中のＳＦ、ＳＰ、ＬＨＳＶ、Ｐは、前記式３の説明と同様である。前記式１０中、Ｇは、上述した情報取得ステップで取得した原料油に関する情報及び運転条件に関する情報に基づき代入される値である。具体的には、Ｇは、水素の供給量（Ｎｍ３／時間）を原料油の供給量（ｋＬ／時間）で除すことにより求めることができる。
【0061】
本明細書において、「基準水素／原料油比」とは、実際の反応条件の標準的な水素／原料油比を意味する。後述する水素／原料油比係数ｂを決定する際に用いる水素／原料油比の平均値として求められる。
【0062】
前記式１０において、（ＧＢ／Ｇ）ｂで表される項は、水素／原料油比依存性を示す項であり、上述した通り、Ｇを設定値とし、一定の条件で運転する場合、定数となる。
【0063】
ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式３で説明した方法と同じ方法で求めることができる。以下、ＧＢ、ｂ、α’の求め方を説明する。
【0064】
ＧＢ、ｂはＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲと同様に、使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。以下、ＧＢ、ｂの求め方について説明する。
【0065】
（ＧＢ、ｂの求め方）
本実施形態の基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。
本実施形態の基本劣化パラメータの求め方は、上述の実機での反応、又は実機運転条件に基づいてベンチスケールでの反応で得られたデータから解析される触媒の劣化挙動に基づく。
【0066】
ＬＨＳＶ、水素分圧、原料油中の硫黄濃度、生成油中の硫黄濃度を一定の条件とし、水素／原料油比Ｇｈの条件で、一定期間反応を行う。触媒は反応により劣化するため、生成油中の硫黄濃度を一定の値とするために、反応温度を上げながら運転を行う。反応時間を横軸に、実測の反応温度を縦軸にプロットし、回帰直線を引くと、ｙ＝ａｈｘ＋ｂｈ（０＜ａｈである。）で表される直線が得られる。この直線におけるｂｈが前記式６におけるＴＳＯＲ’となる。前記式６において、ＴＳＯＲ’にｂｈを代入し、Ｔｔ
’に実測の反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、回帰直線を引くと、ｙ＝－ａｈ
’ｘ（|―ａｈ
’|＝ａｈ
’である。）で表される直線が得られる。ａｈ
’は、触媒の劣化速度を表している。
【0067】
ｈ種類の水素／原料油比Ｇｈについて同様の反応を行い、同様にして、ｈ個のａｈ、ｂｈを求め、上記と同様の方法により、ｈ個のａｈ
’を得る。ｈは３以上の整数である。ｈの数が多い方が、より精度の高いｂを得ることができる。一方、ｈの数が多すぎると、ｂを得るために時間が掛かり効率的ではない。本実施形態においては、ｈは３～２０であることが好ましく、３～１０であることがより好ましい。
このようにして得られたｈ個のａｈ’及びＧｈを下式１１にそれぞれ代入する。下式１１は、水素／原料油比係数、及び基準水素／原料油比を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
このようにして得られたｈ個のａｈ’及びＧｈを下式１１にそれぞれ代入する。下式１１は、水素／原料油比係数、及び基準水素／原料油比を算出可能な式であり、本願の発明者らが、実機の運転結果等を基に初めて見出した式である。
【0068】
ｌｎ（ａｈ’）＝－ｂｌｎ（Ｇｈ）＋Ｂ３ 式１１
前記式１１中、Ｂ３は０とすることができる。
前記式１１中、Ｂ３は０とすることができる。
【0069】
ｈ個のａｈ’及びＧｈの組み合わせに関し、ｌｎ（ａｈ’）を縦軸に、ｌｎ（Ｇｈ）を横軸にプロットし、回帰直線を引き、その傾きを求める。この傾きが水素／原料油比係数ｂである。
【0070】
また、上記ｈ種類の水素／原料油比Ｇｈを平均することにより、ＧＢを求めることができる。
【0071】
ｂ及びＧＢを求める上でのＬＨＳＶ、水素分圧、原料油中の硫黄濃度、生成油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのようなＬＨＳＶとしては、例えば０．５～２．５ｈ－１であり、水素分圧としては、例えば３～７ＭＰａであり、原料油中の硫黄濃度としては、例えば０．５～３．０質量％であり、生成油中の硫黄濃度としては、例えば０．３質量％以下である。
ｈ種類の水素／原料油比Ｇｈも同様に実機運転条件に即した条件とすることが好ましい。このようなＧｈとしては、例えば１００～４００［Ｎｍ３／ｋＬ］である。反応期間としては、例えば５０～６００日である。
そのようなＬＨＳＶとしては、例えば０．５～２．５ｈ－１であり、水素分圧としては、例えば３～７ＭＰａであり、原料油中の硫黄濃度としては、例えば０．５～３．０質量％であり、生成油中の硫黄濃度としては、例えば０．３質量％以下である。
ｈ種類の水素／原料油比Ｇｈも同様に実機運転条件に即した条件とすることが好ましい。このようなＧｈとしては、例えば１００～４００［Ｎｍ３／ｋＬ］である。反応期間としては、例えば５０～６００日である。
【0072】
また、α’は、前記式３の代わりに前記式１０を使用する以外は、劣化関数１において説明したαの求め方（（α及びＴＳＯＲの求め方１）及び（α及びＴＳＯＲの求め方２））と同じ方法により求めることができる。
【0073】
＜劣化関数２＞
本実施形態においては、コークの堆積による触媒の劣化に関するコーク劣化関数は、触媒の易失活活性種の劣化に関する易失活活性種劣化関数と、触媒の難失活活性種の劣化に関する難失活活性種劣化関数から構成される下式１２で表される劣化関数２が好ましい。
本実施形態においては、コークの堆積による触媒の劣化に関するコーク劣化関数は、触媒の易失活活性種の劣化に関する易失活活性種劣化関数と、触媒の難失活活性種の劣化に関する難失活活性種劣化関数から構成される下式１２で表される劣化関数２が好ましい。
【0074】
Φ＝ｋ1×ｅｘｐ（－Ｄ１ｔ）＋ｋ２×ｅｘｐ（－Ｄ２ｔ） 式１２
前記式１２中、ｋ１は触媒の易失活活性種の活性点係数であり、ｋ２は触媒の難失活活性種の活性点係数であり、活性点係数は両活性種の相対反応速度定数を表している。Ｄ１は触媒の易失活活性種の劣化係数であり、Ｄ２は触媒の難失活活性種の劣化係数であり、ｔは反応経過日数（日）であり、ｋ１＋ｋ２＝１である。
前記式１２中、ｋ１は触媒の易失活活性種の活性点係数であり、ｋ２は触媒の難失活活性種の活性点係数であり、活性点係数は両活性種の相対反応速度定数を表している。Ｄ１は触媒の易失活活性種の劣化係数であり、Ｄ２は触媒の難失活活性種の劣化係数であり、ｔは反応経過日数（日）であり、ｋ１＋ｋ２＝１である。
【0075】
減圧蒸留軽油の水素化処理反応においては、上述した通り、コークの堆積により触媒が劣化するため、生成油中の硫黄含有量を一定レベル以下に保つため、反応温度を上げて運転を行う必要がある。減圧蒸留軽油の水素化処理反応において、反応開始初期には、反応温度が急激に上昇する。この反応温度の急激な上昇は、反応開始初期の触媒の急激な劣化を意味する。一方、反応中期以降においては、反応温度は緩やかに上昇する。この反応温度の緩やかな上昇は、反応中期以降の触媒の緩やかな劣化を意味する。
【0076】
すなわち、減圧蒸留軽油の水素化処理反応においては、反応開始初期の触媒の急激な劣化と、反応中期以降の触媒のゆるやかな劣化が起こっていることが、反応時間に対する反応温度のプロファイルから示唆される。
【0077】
本願の発明者らは、上記反応時間に対する反応温度のプロファイルに基づき、触媒には、反応開始初期に失活する易失活活性種と、反応中期以降に失活する難失活活性種が存在するとの仮定に基づき、劣化関数１をさらに改良し、劣化関数２を見出した。その結果、劣化関数２によれば、劣化関数１よりも、より精度よく触媒の劣化度を算出可能となることを見出した。易失活活性種とは、主に反応初期に活性が失われる活性種であり、難失活活性種とは、反応中期以降に活性が失われる活性種を意味する。
【0078】
前記式１２において、ｋ１は、触媒の易失活活性種の活性点係数を示し、ｋ２は、触媒の難失活活性種の活性点係数を示す。ｋ１及びｋ２は触媒固有の定数であり、その求め方は後述する。
【0079】
Ｄ１は下式１３により、Ｄ２は下式１４により求めることができる。
【0080】
【数5】
【0081】
【数6】
【0082】
前記式１３及び前記式１４中、ＳＦは任意の反応ｔ日経過時の原料油中の硫黄濃度（質量％）であり、ＳＰは任意の反応ｔ日経過時の生成油中の硫黄濃度（質量％）であり、ｎは、減圧蒸留軽油を含む原料油の水素化処理反応の反応次数であり、ＬＨＳＶは任意の反応ｔ日経過時の液空間速度（ｈ－１）であり、ＰＢは基準水素分圧（ＭＰａ）であり、Ｐは任意の反応ｔ日経過時の水素分圧（ＭＰａ）であり、ａは水素分圧係数であり、Ｅｃはコーク劣化の活性化エネルギー（ｋＪ／ｍｏｌ）であり、Ｒは気体定数：０．００８３１（ｋＪ／（ｍｏｌ・Ｋ））であり、ＴＢは基準反応温度（Ｋ）であり、ＴＳＯＲは０日目の要求温度（Ｋ）である。
前記式１３中、α１は易失活活性点の触媒定数（触媒の劣化速度を表す定数）であり、前記式１４中、α２は難失活活性点の触媒定数（触媒の劣化速度を表す定数）である。
前記式１３中、α１は易失活活性点の触媒定数（触媒の劣化速度を表す定数）であり、前記式１４中、α２は難失活活性点の触媒定数（触媒の劣化速度を表す定数）である。
【0083】
前記式１３及び前記式１４中、ＳＦ、ＳＰ、ＬＨＳＶ、Ｐは、前記式３の説明と同様に、上述した情報取得ステップで取得した原料油に関する情報、生成油に関する情報、及び運転条件に関する情報に基づき代入される値である。なお、ＬＨＳＶは原料油の供給量（体積／ｈ）を触媒充填量（体積）で除すことにより求めることができる。
【0084】
上述した通り、ＳＦ、ＬＨＳＶ、Ｐは制御することが可能なパラメータである。また、ＳＰは目的とする生成油の硫黄濃度である。すなわち、前記式１３及び前記式１４により、上記ＳＦ、ＳＰ、ＬＨＳＶ、及びＰの反応条件における、触媒の易失活活性種の劣化係数、及び触媒の難失活活性種の劣化係数をそれぞれ算出することができる。減圧蒸留軽油を含む原料油の水素化処理反応の反応次数であるｎの求め方は後述する。
【0085】
（基本劣化パラメータの求め方）
前記式１３及び前記式１４中、α１、α２、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式３と同様、定数であり、これらのパラメータを総称して「基本劣化パラメータ２」という。基本劣化パラメータ２は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。前記式１３及び前記式１４中、ＰＢ、ａ、Ｅｃ、ＴＢは前記式３と同じ方法により求めることができる。
一方、α１、α２、ＴＳＯＲは例えば以下の２種類の方法により求めることができる。また、前記式１２における触媒の易失活活性種の活性点数係数ｋ１、触媒の難失活活性種の活性点数係数ｋ２も同時に以下のように求めることができる。
前記式１３及び前記式１４中、α１、α２、ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式３と同様、定数であり、これらのパラメータを総称して「基本劣化パラメータ２」という。基本劣化パラメータ２は使用する触媒に応じて定められるパラメータであり、実機で反応を行いながら求めてもよいし、実機運転条件に基づいてベンチスケールにおいて事前に求めてもよい。本実施形態においては、実機運転条件に基づいてベンチスケールにおいて事前に求めることが好ましい。前記式１３及び前記式１４中、ＰＢ、ａ、Ｅｃ、ＴＢは前記式３と同じ方法により求めることができる。
一方、α１、α２、ＴＳＯＲは例えば以下の２種類の方法により求めることができる。また、前記式１２における触媒の易失活活性種の活性点数係数ｋ１、触媒の難失活活性種の活性点数係数ｋ２も同時に以下のように求めることができる。
【0086】
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方１）
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。反応時間を横軸に、反応温度を縦軸にプロットする。これらのプロットの回帰直線を引く場合、上述した通り、反応開始初期（ｘ１～ｘｎ）の、急激な反応温度の上昇に相関するｙ＝ａ１ｘ＋ｂ１で表される直線と、反応中期以降（ｘｎ＋１～ｘｍ）の、緩やかな反応温度の上昇に相関するｙ＝ａ２ｘ＋ｂ２で表される２本の直線が得られる。前記式において、ａ１＞ａ２＞０であり、ｂ１＜ｂ２であり、ｘ１＜ｘｎ＜ｘｎ＋１＜ｘｍである。ｘｎ、ｘｍは反応開始からのｎ（ｍ）番目のプロットに相関した反応時間を意味する。
実機又はベンチスケールにおいて、想定する実機運転条件のＬＨＳＶ、水素分圧、水素／原料油比、原料油中の硫黄濃度、生成油中の硫黄濃度となるように、一定期間反応を行う。触媒は反応により劣化するため、反応温度を上げながら運転を行う。反応時間を横軸に、反応温度を縦軸にプロットする。これらのプロットの回帰直線を引く場合、上述した通り、反応開始初期（ｘ１～ｘｎ）の、急激な反応温度の上昇に相関するｙ＝ａ１ｘ＋ｂ１で表される直線と、反応中期以降（ｘｎ＋１～ｘｍ）の、緩やかな反応温度の上昇に相関するｙ＝ａ２ｘ＋ｂ２で表される２本の直線が得られる。前記式において、ａ１＞ａ２＞０であり、ｂ１＜ｂ２であり、ｘ１＜ｘｎ＜ｘｎ＋１＜ｘｍである。ｘｎ、ｘｍは反応開始からのｎ（ｍ）番目のプロットに相関した反応時間を意味する。
【0087】
上述のａ１がα１に相関する値であり、ｂ１がｋ１＋ｋ２に相関する値であり、ＴＳＯＲである。また、ａ２がα２に相関する値であり、ｂ２がｋ２に相関する値である。次に、ｙ＝ａ１ｘ＋ｂ１とｙ＝ａ２ｘ＋ｂ２での交点（ｘｉｐ、ｙｉｐ）を計算する。この交点はｙ＝ａ１ｘ＋ｂ１とｙ＝ａ２ｘ＋ｂ２の変曲点を意味する。すなわち、（ｘｉｐ、ｙｉｐ）までは、触媒の易失活活性種のみが存在するとみなし、（ｘｉｐ、ｙｉｐ）からは触媒の易失活活性種及び触媒の難失活活性種の２つの活性種が存在するとみなす。
【0088】
前記式６において、ＴＳＯＲ’にｂ１を代入し、Ｔｔ’に反応温度を代入すると、任意の反応ｔ日経過時における劣化度Φ’が得られる。反応時間を横軸に、Φ’の対数を縦軸にプロットし、ｘ１～ｘｉｐまでの回帰直線を引くと、ｙ＝－ａ１
’ｘ （|－ａ１
’|＝ａ１
’である。）で表される直線が得られる。また、ｘiｐ～ｘｍまでの回帰直線を引くと、ｙ＝－ａ２
’ｘ－ｂ２
’ （|－ａ２
’|＝ａ２
’であり、|－ｂ２
’|＝ｂ２
’である。）で表される直線が得られる。ａ１
’はα１に相関した値であり、ａ２
’はα２に相関した値であり、ｂ２
’はｋ２に相関した値である。
【0089】
ｋ２は、上述の回帰直線で得られたｂ２
’を下式１５に代入して求めることができる。ｋ１はｋ１＋ｋ２＝１であるから、ｋ１＝１－ｋ２より求めることができる。
ｋ２＝ｅｘｐ（－ｂ２ ’） 式１５
ｋ２＝ｅｘｐ（－ｂ２ ’） 式１５
【0090】
運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（ｂ１）を前記式１３及び前記式１４に代入してＤ１及びＤ２を得る。求められたＤ１、Ｄ２、ｋ１、ｋ２を前記式１２に代入するとΦが得られる。この場合、Φはα１とα２の関数となる。反応時間を横軸に、Φの対数を縦軸にプロットし、α２＝０としてα１を０＜α１になるように変えてｘ１～ｘｉｐまでの回帰直線を引くと、ｙ＝－ａα1
”ｘ（|－ａα１
“|＝ａα１
”である。）で表される複数の直線がα１の値ごとに得られる。各直線におけるａα１
”と、上述のａ１
’が等しい値となるときのα１を前記式１３におけるα１とすることができる。
次に、上述の方法で得られたα１とα２の関数であるΦに、得られたα１を代入し、α２を０＜α２になるように変えてｘiｐ～ｘｍまでの回帰直線を引くと、ｙ＝－ａα２ ”ｘ－ｂα２ “（|－ａα２ ”|＝ａα２ “である。）で表される複数の直線がα２の値ごとに得られる。各直線におけるａα２と、上述のａ２ ’が等しい値となるときのα２を前記式１４におけるα２とすることができる。
次に、上述の方法で得られたα１とα２の関数であるΦに、得られたα１を代入し、α２を０＜α２になるように変えてｘiｐ～ｘｍまでの回帰直線を引くと、ｙ＝－ａα２ ”ｘ－ｂα２ “（|－ａα２ ”|＝ａα２ “である。）で表される複数の直線がα２の値ごとに得られる。各直線におけるａα２と、上述のａ２ ’が等しい値となるときのα２を前記式１４におけるα２とすることができる。
【0091】
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方２）
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方１）と同様の反応を行い、ＴＳＯＲ（ｂ１）を得る。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（ｂ１）を前記式１３及び前記式１４に代入してＤ１及びＤ２を得る。求められたＤ１、Ｄ２を前記式１２に代入するとΦが得られる。この場合、Φはα１、α２、ｋ１、ｋ２の関数となる。得られたΦを前記式６のΦ’に、ＴＳＯＲ（すなわち、ｂ１）を前記式６のＴＳＯＲ ’に代入して、Ｔｔ ’について整理すると、Ｔｔ ’はα１、α２、ｋ１、ｋ２の関数となる。α２を０とし、ｋ２＝１－ｋ１とするとＴｔ ’はα１、ｋ１の関数となる。ｘ１～ｘｉｐにおける実測の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）が１となるときのα１、ｋ１の組み合わせを求め、このα１を前記式１３におけるα１とすることができる。なお、この時のｋ１は仮値である。同様にＭ個の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα１を前記式１３におけるα１とすることが好ましい。Ｍは１０以上整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
得られたα１を使用し、ｋ１＝１－ｋ２とするとＴｔ ’はα２、ｋ２の関数となる。ｘiｐ～ｘｍにおける実測の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）が１となるときのα２、ｋ２の組み合わせを求め、このα２、ｋ２を前記式１４におけるα２、ｋ２とすることができる。得られたｋ２をｋ１＝１－ｋ２に代入することによりｋ１を求め、このｋ１を前記式１３におけるｋ１とすることができる。同様にＬ個の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα２、ｋ２を前記式１４におけるα２、ｋ２とすることが好ましい。また、得られたｋ２から求められたｋ１を前記式１３におけるｋ１とすることが好ましい。Ｌは１０以上整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
（α１、α２、ＴＳＯＲ、ｋ１、及びｋ２の求め方１）と同様の反応を行い、ＴＳＯＲ（ｂ１）を得る。運転条件であるＳＰ、ＳＦ、ＬＨＳＶ、Ｐ、及び前記の方法で求めたＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲ（ｂ１）を前記式１３及び前記式１４に代入してＤ１及びＤ２を得る。求められたＤ１、Ｄ２を前記式１２に代入するとΦが得られる。この場合、Φはα１、α２、ｋ１、ｋ２の関数となる。得られたΦを前記式６のΦ’に、ＴＳＯＲ（すなわち、ｂ１）を前記式６のＴＳＯＲ ’に代入して、Ｔｔ ’について整理すると、Ｔｔ ’はα１、α２、ｋ１、ｋ２の関数となる。α２を０とし、ｋ２＝１－ｋ１とするとＴｔ ’はα１、ｋ１の関数となる。ｘ１～ｘｉｐにおける実測の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）が１となるときのα１、ｋ１の組み合わせを求め、このα１を前記式１３におけるα１とすることができる。なお、この時のｋ１は仮値である。同様にＭ個の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα１を前記式１３におけるα１とすることが好ましい。Ｍは１０以上整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
得られたα１を使用し、ｋ１＝１－ｋ２とするとＴｔ ’はα２、ｋ２の関数となる。ｘiｐ～ｘｍにおける実測の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）が１となるときのα２、ｋ２の組み合わせを求め、このα２、ｋ２を前記式１４におけるα２、ｋ２とすることができる。得られたｋ２をｋ１＝１－ｋ２に代入することによりｋ１を求め、このｋ１を前記式１３におけるｋ１とすることができる。同様にＬ個の反応温度Ｔｏｂｓに対するＴｔ ’の比（Ｔｔ ’／Ｔｏｂｓ）を計算し、これらの平均が最も１に近づくときのα２、ｋ２を前記式１４におけるα２、ｋ２とすることが好ましい。また、得られたｋ２から求められたｋ１を前記式１３におけるｋ１とすることが好ましい。Ｌは１０以上整数であり、１０～５００であることが好ましく、５０～２００であることがより好ましい。
【0092】
このように、α１、α２、ＴＳＯＲ、ｋ１、及びｋ２を求めるためには、ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２を得る必要がある。ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２は例えば、以下のように求めることができる。
【0093】
上述の方法で、反応時間を横軸に、反応温度を縦軸にプロットする。反応開始から反応終了までの回帰直線を引くと、ｙ＝ａ’ｘ＋ｂ’で表される直線が得られる。この直線は、変曲点を考慮していない直線である。反応終了からプロットを順に削除し、相関係数が１に近づくように前記ｙ＝ａ’ｘ＋ｂ’を調整し、ｙ＝ａ１’ｘ＋ｂ１’を得る。同様に、反応開始からプロットを順に削除し、相関係数が１に近づくように前記ｙ＝ａ’ｘ＋ｂ’を調整し、ｙ＝ａ２’ｘ＋ｂ２’を得る。ｙ＝ａ１’ｘ＋ｂ１’とｙ＝ａ２’ｘ＋ｂ２’の相関係数の平均が最も１に近づく場合の各直線がｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２となる。なお、全プロットは、ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２のどちらか一方に属するようにする。
【0094】
ｙ＝ａ１ｘ＋ｂ１、ｙ＝ａ２ｘ＋ｂ２を得るために必要な反応時間は通常１００日以上である。また、一般に、上記、ｙ＝ａ’ｘ＋ｂ’の相関関数が０．５以上となるまで反応を行えば充分である。
【0095】
（劣化関数２の変形例）
劣化関数２の変形例を以下に説明する。劣化関数１としては下式１６で表される劣化関数２－１を使用してもよい。
劣化関数２の変形例を以下に説明する。劣化関数１としては下式１６で表される劣化関数２－１を使用してもよい。
【0096】
Φ＝ｋ1×ｅｘｐ（－Ｄ１’ｔ）＋ｋ２×ｅｘｐ（－Ｄ２’ｔ） 式１６
前記式１６中、ｋ１、ｋ２、ｔは前記式１２と同じであり、Ｄ１’は触媒の易失活活性種の劣化係数であり、Ｄ２’は触媒の難失活活性種の劣化係数である。
前記式１６中、ｋ１、ｋ２、ｔは前記式１２と同じであり、Ｄ１’は触媒の易失活活性種の劣化係数であり、Ｄ２’は触媒の難失活活性種の劣化係数である。
【0097】
Ｄ１’は下式１７により、Ｄ２’は下式１８により求めることができる。
【0098】
【数7】
【0099】
【数8】
【0100】
前記式１７及び前記式１８中、ＳＦ、ＳＰ、ｎ、ＬＨＳＶ、ＰＢ、Ｐ、ａ、Ｅｃ、Ｒ、ＴＢ、ＴＳＯＲは前記式１３及び前記式１４と同じであり、ＧＢ、Ｇ、ｂは前記式１０と同じであり、前記式１７中、α１’は易失活活性点の触媒定数（触媒の劣化速度を表す定数）であり、前記式１８中、α２’は難失活活性点の触媒定数（触媒の劣化速度を表す定数）である。
【0101】
ＰＢ、ａ、Ｅｃ、ＴＢ、ＴＳＯＲは前記式１３及び前記式１４で説明した方法と同じ方法で求めることができ、ＧＢ、ｂは、前記式１０で説明した方法と同じ方法で求めることができる。
【0102】
また、α１’、α２’は、前記式１２の代わりに前記式１６を、前記式１３の代わりに前記式１７を、前記式１４の代わりに前記式１８を使用する以外は、劣化関数２において説明したα１、α２の求め方と同じ方法により求めることができる。
【0103】
≪反応温度算出ステップ≫
本実施形態の反応温度算出ステップは、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出するステップである。反応温度は、アレニウスの式に基づいた劣化速度式により算出することが好ましい。
本実施形態の反応温度算出ステップは、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出するステップである。反応温度は、アレニウスの式に基づいた劣化速度式により算出することが好ましい。
【0104】
＜劣化速度式＞
劣化速度式は前記式５で表されるアレニウスの式に基づいた式である。前記式６の算出方法と同様に、前記式１及び前記式５より下式１９が、導かれる。
劣化速度式は前記式５で表されるアレニウスの式に基づいた式である。前記式６の算出方法と同様に、前記式１及び前記式５より下式１９が、導かれる。
【0105】
【数9】
前記式１９を整理すると下式２０が得られる。
【0106】
【数10】
【0107】
前記式２０に、上述の方法で得られたＴＳＯＲ、Φを代入することにより、任意の反応ｔ日経過時におけるＴｔ（Ｋ）を得ることができる。なお、前記式２０における脱硫反応の活性化エネルギーは以下のように求めることができる。
【0108】
（脱硫反応の活性化エネルギーの求め方）
脱硫反応の活性化エネルギーは、前記式５で表されるアレニウスの式に基づいて、本分野で公知の方法により求めることができる。以下、一例を説明する。
脱硫反応の活性化エネルギーは、前記式５で表されるアレニウスの式に基づいて、本分野で公知の方法により求めることができる。以下、一例を説明する。
【0109】
はじめに、減圧蒸留軽油を含む原料油の脱硫反応の反応次数を求める。反応温度、水素分圧、水素／原料油比、原料油中の硫黄濃度は一定の条件とし、ＬＨＳＶ（ｘ）の条件で、反応を行い、生成油中の硫黄濃度を測定する。下式２１で表される脱硫反応速度式のＳＦに原料油中の硫黄濃度を、ＳＰに得られた生成油中の硫黄濃度を、ＬＨＳＶにＬＨＳＶ（ｘ）を代入する。得られた左辺の結果を縦軸に、１／ＬＨＳＶを横軸にプロットする。なお、この場合、縦軸はｎの関数となる。
【0110】
【数11】
【0111】
ｘ種類のＬＨＳＶ（ｘ）について同様の反応を行い、ｘ個の上記プロットを得る。得られたプロットを基に原点を通るように回帰直線を引き、ｙ＝ｃｘで表される直線を得る。ｙは（１／ｎ－１（（１／ＳＰ
ｎ－１）－（１／ＳＦ
ｎ－１）））であり、ｘは１／ＬＨＳＶであり、ｃはｋとなる。ｅｘｃｅｌ等で相関関数を求め、相関係数が最も１に近づくｎを求める。得られたｎが反応次数となる。なお、ｎは小数点第１位のオーダーで求めればよい。
【0112】
上記ｘは３以上の整数である。ｘの数が多い方が、より精度の高いｎを得ることができる。一方、ｘの数が多すぎると、ｎを得るために時間が掛かり効率的ではない。本実施形態においては、ｘは３～２０であることが好ましく、３～１０であることがより好ましい。
【0113】
ｎを求める上での反応温度、水素分圧、水素／原料油比、原料油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのような反応温度としては、例えば３００～４２０℃であり、水素分圧は３～７ＭＰａであり、水素／原料油比は１００～４００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度は０．５～３．０質量％である。
ｘ種類のＬＨＳＶ（ｘ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＬＨＳＶ（ｘ）としては、０．５～２．５ｈ－１である。
そのような反応温度としては、例えば３００～４２０℃であり、水素分圧は３～７ＭＰａであり、水素／原料油比は１００～４００［Ｎｍ３/ｋＬ］であり、原料油中の硫黄濃度は０．５～３．０質量％である。
ｘ種類のＬＨＳＶ（ｘ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＬＨＳＶ（ｘ）としては、０．５～２．５ｈ－１である。
【0114】
前記式５で表されるアレニウスの式中の活性化エネルギーＥを脱硫の活性化エネルギーＥａとし、両辺の自然対数を取ると、下式２２で表される式となる。
【数12】
【0115】
水素分圧、水素／原料油比、ＬＨＳＶ、原料油中の硫黄濃度は一定の条件とし、反応温度Ｔ（ｙ）の条件で反応を行い、生成油中の硫黄濃度を測定する。前記式２１で表される脱硫反応速度式のＳＦに原料油中の硫黄濃度を、ＳＰに得られた生成油中の硫黄濃度を、ＬＨＳＶにＬＨＳＶを、上記得られたｎを代入して反応速度定数ｋを求める。得られた反応速度定数を前記式２２に代入し、得られた左辺の結果（ｌｎｋ）を縦軸に、１／Ｔ（１／Ｔ（ｙ））を横軸にプロットする。
【0116】
ｙ種類の反応温度Ｔ（ｙ）について同様の反応を行い、ｙ個の上記プロットを得る。得られたプロットより、回帰直線を引き、その傾きを求める。この傾きはＥａ／Ｒであるため、傾きからＲを除すことにより脱硫の活性化エネルギーＥａを求めることができる。
【0117】
上記ｙは３以上の整数である。ｙの数が多い方が、より精度の高いＥａを得ることができる。一方、ｙの数が多すぎると、Ｅａを得るために時間が掛かり効率的ではない。本実施形態においては、ｙは３～２０であることが好ましく、３～１０であることがより好ましい。
【0118】
Ｅａを求める上での水素分圧、水素／原料油比、ＬＨＳＶ、原料油中の硫黄濃度は、実機運転条件に即した条件とすることが好ましい。
そのような水素分圧としては、例えば３～７ＭＰａであり、水素／原料油比は１００～４００［Ｎｍ３/ｋＬ］であり、ＬＨＳＶは０．５～２．５ｈ－１であり、原料油中の硫黄濃度は０．５～３．０質量％である。
ｙ種類の反応温度Ｔ（ｙ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＴ（ｙ）としては、３００～４２０℃である。
そのような水素分圧としては、例えば３～７ＭＰａであり、水素／原料油比は１００～４００［Ｎｍ３/ｋＬ］であり、ＬＨＳＶは０．５～２．５ｈ－１であり、原料油中の硫黄濃度は０．５～３．０質量％である。
ｙ種類の反応温度Ｔ（ｙ）も同様に実機運転条件に即した条件とすることが好ましい。このようなＴ（ｙ）としては、３００～４２０℃である。
【0119】
このように得られた各パラメータを前記式２０に代入することにより、所定の反応条件を達成するために必要な反応温度Ｔｔを求めることができる。
【0120】
実測の反応温度Ｔｏｂｓに対する、本実施形態の情報処理方法により求められる反応温度Ｔｔの割合であるＴｔ／Ｔｏｂｓは、℃換算で、０．９７～１．０３であることが好ましく、０．９８５～１．０１５であることがより好ましい。Ｔｔ／Ｔｏｂｓが前記範囲内である場合、精度よく反応温度を推定できていると判断することができる。
【0121】
≪情報出力ステップ≫
このようにして得られた前記反応温度を示す情報を出力する情報出力ステップ（図１のＳ４）をさらに有してもよい。例えば、Ｓ４は出力部１４によって実行される。
このようにして得られた前記反応温度を示す情報を出力する情報出力ステップ（図１のＳ４）をさらに有してもよい。例えば、Ｓ４は出力部１４によって実行される。
【0122】
≪減圧蒸留軽油を含む原料油の水素化処理反応≫
減圧蒸留軽油を含む原料油の水素化処理反応について概要を説明する。
減圧蒸留軽油は、常圧蒸留残渣油を減圧蒸留して得られる、沸点範囲が２００～５５０℃の留分である。減圧蒸留軽油の密度は０．９１～０．９５ｇ／ｍＬである。原料油中の減圧蒸留軽油の含有量としては、例えば３０体積％以上であってもよく、５０体積％以上等であってもよい。
なお、原料油に含まれる減圧蒸留軽油以外の油種としては、原油を常圧蒸留装置で常圧蒸留して得られる常圧重油、水素化分解重油等潤滑油基油の溶剤抽出により抽出除去される油分の中で特に重質な油分である重質エキストラクト、常圧蒸留残渣、流動接触分解残油、脱礫油等が挙げられる。
減圧蒸留軽油を含む原料油の水素化処理反応について概要を説明する。
減圧蒸留軽油は、常圧蒸留残渣油を減圧蒸留して得られる、沸点範囲が２００～５５０℃の留分である。減圧蒸留軽油の密度は０．９１～０．９５ｇ／ｍＬである。原料油中の減圧蒸留軽油の含有量としては、例えば３０体積％以上であってもよく、５０体積％以上等であってもよい。
なお、原料油に含まれる減圧蒸留軽油以外の油種としては、原油を常圧蒸留装置で常圧蒸留して得られる常圧重油、水素化分解重油等潤滑油基油の溶剤抽出により抽出除去される油分の中で特に重質な油分である重質エキストラクト、常圧蒸留残渣、流動接触分解残油、脱礫油等が挙げられる。
【0123】
減圧蒸留軽油を含む原料油の水素化処理反応は、水素存在下、減圧蒸留軽油を含む原料油と水素化処理触媒とを接触処理することにより行うことができる。
水素化処理触媒は、特に限定されるものではなく、本分野において公知の水素化処理触媒を使用することができる。触媒の担体として、種々のものが使用でき、例えばシリカ、アルミナ、ボリア、マグネシア、チタニア、シリカ－アルミナ、シリカ－マグネシア、シリカ－ジルコニア、シリカ－トリア、シリカ－ベリリア、シリカ－チタニア、シリカ－ボリア、アルミナ－ジルコニア、アルミナ－チタニア、アルミナ－ボリア、アルミナ－クロミア、チタニア－ジルコニア、シリカ－アルミナ－トリア、シリカ－アルミナ－ジルコニア、シリカ－アルミナ－マグネシア、シリカ－マグネシア－ジルコニアなど、又はこれらの２種以上の混合物が挙げられる。これらの無機酸化物のうち、好ましいものとしては、アルミナ、シリカ－アルミナ、アルミナ－チタニア、アルミナ－ボリア、アルミナ－ジルコニアが挙げられ、特に好ましくは、アルミナが挙げられ、アルミナの中でもγアルミナが特に好ましい。これらの無機酸化物は、１種単体で用いてもよいし、２種以上を組み合わせて用いてもよい。
水素化処理触媒は、特に限定されるものではなく、本分野において公知の水素化処理触媒を使用することができる。触媒の担体として、種々のものが使用でき、例えばシリカ、アルミナ、ボリア、マグネシア、チタニア、シリカ－アルミナ、シリカ－マグネシア、シリカ－ジルコニア、シリカ－トリア、シリカ－ベリリア、シリカ－チタニア、シリカ－ボリア、アルミナ－ジルコニア、アルミナ－チタニア、アルミナ－ボリア、アルミナ－クロミア、チタニア－ジルコニア、シリカ－アルミナ－トリア、シリカ－アルミナ－ジルコニア、シリカ－アルミナ－マグネシア、シリカ－マグネシア－ジルコニアなど、又はこれらの２種以上の混合物が挙げられる。これらの無機酸化物のうち、好ましいものとしては、アルミナ、シリカ－アルミナ、アルミナ－チタニア、アルミナ－ボリア、アルミナ－ジルコニアが挙げられ、特に好ましくは、アルミナが挙げられ、アルミナの中でもγアルミナが特に好ましい。これらの無機酸化物は、１種単体で用いてもよいし、２種以上を組み合わせて用いてもよい。
【0124】
前記担体に活性成分として含有させる金属は、周期表第６族金属及び第８～１０族金属の中から選ばれる少なくとも１種類以上の金属であり、好ましくはモリブデン、タングステン、コバルト及びニッケルの金属である。これらの金属は、金属状態又は金属酸化物、金属硫化物の何れの形態でも有効であり、また、イオン交換などにより金属が触媒担体と結合した形態で存在してもよい。この金属成分の含有量は、通常、触媒基準かつ酸化物換算で、約１０～２５質量％の範囲内である。金属含有量が１０質量％より少ないと、活性点として働く金属の絶対量が少ないために、脱硫活性を始めとする水素化処理活性（以下、簡単に水素化処理活性と言う）が発現せず、逆に担持される金属の含有量が２５質量％より多すぎると、金属の凝集が起こり活性点の数が減少し、その結果、水素化処理活性が却って低下するからである。更に、必要に応じて、元素周期律表第６族金属及び第８族金属からなる活性金属に加えて、リン、ホウ素、亜鉛、ジルコニア等を含ませることができる。本発明方法を適用するに当たり、触媒層の形態には制約はなく、例えば固定床、移動床、流動床等の触媒層の反応器に適用できる。
【0125】
減圧蒸留軽油を含む原料油の水素化処理反応における条件としては、一般に反応温度が３００～４２０℃、好ましくは３４０～４００℃であり、水素分圧が３～７ＭＰａ、好ましくは４～６ＭＰａであり、ＬＨＳＶが０．５～２．５ｈｒ－１、好ましくは０．６～１．５ｈｒ－１であり、水素／原料油比が１００～４００［Ｎｍ３/ｋＬ］、好ましくは２００～３００［Ｎｍ３/ｋＬ］である。
【0126】
減圧蒸留軽油を含む原料油中の硫黄濃度は、通常０．５～３．０質量％である。また、生成油中の硫黄濃度は、通常０．０１～０．３質量％である。
【0127】
≪反応温度算出装置≫
本実施形態の反応温度算出装置は、常圧蒸留残渣油を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む。本実施形態の反応温度算出装置は、算出された反応温度を示す情報を出力する出力部を有していてもよい。
本実施形態の反応温度算出装置は、常圧蒸留残渣油を減圧蒸留して得られる減圧蒸留軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得する取得部と、前記取得部で取得した前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、算出した前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する演算部と、を含む。本実施形態の反応温度算出装置は、算出された反応温度を示す情報を出力する出力部を有していてもよい。
【0128】
本実施形態の反応温度算出装置１は、パーソナルコンピューターやサーバー装置や専用装置等の情報処理装置を用いて構成される。反応温度算出装置１は、１台又は複数台の情報処理装置を用いて構成されてもよい。例えば、反応温度算出装置１は、クラスタマシンとして構築されてもよいし、クラウドとして構築されてもよいし、どのような態様で構築されてもよい。反応温度算出装置１は、例えば、図４に示されるように取得部１１と、取得部からの情報を処理する計算機本体１２とを有する。反応温度算出装置１は、計算機本体１２において処理された情報を外部に出力する出力部１４を有していてもよい。これらの構成要素は、例えば、ＣＰＵ（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）等により構成される。
【0129】
取得部１１は、反応のオペレーターによって所定の情報が入力され、この入力により取得した情報を計算機本体１２に送信するものである。本実施形態の取得部１１が取得する情報は、減圧蒸留軽油を含む原料油の水素化処理反応に関して、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報である。反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報としては上述した通りである。例えば、取得部１１は、上述した情報取得ステップを実行する。取得部１１は、反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を取得すればよく、その取得方法には特に限定されない。
【0130】
本実施形態では、取得部１１は単一のキーボードによって構成されている。取得部１１の具体的構成は、限定されず、本実施形態ではキーボードであるが、タッチパネル等であってもよい。なお、各種情報を取得する取得部が別々に構成され、それぞれが独立して計算機本体１２に接続されてもよい。また、取得部１１は、反応器の制御等に用いられる計算機等から有線又は無線により前記の各情報を直接取得するように構成されてもよい。
【0131】
計算機本体１２は、例えば、種々の情報を処理可能な、いわゆるコンピュータである。計算機本体１２は、演算部１３を備える。例えば、この計算機本体１２には、所定のプログラムが組み込まれ、このプログラムの実行によって機能的に演算部１３が構成される。具体的には、この演算部１３において、前記取得部１１で取得した反応を開始してから所定時間経過した際の原料油に関する情報、生成油に関する情報、及び運転条件に関する情報に基づき、劣化関数により、触媒の劣化度を算出し、前記触媒の劣化度に基づき、前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を算出する。劣化関数としては上述した通りである。上述したように反応温度は、例えば劣化速度式から求めることができる。例えば、演算部１３は、上述した劣化度算出ステップ、及び反応温度算出ステップを実行する。演算部１３は、例えば、ＣＰＵ（Central Processing Unit)又はＭＰＵ（Micro Processing Unit）などのプロセッサ及び不揮発性又は揮発性の半導体メモリ（例えば、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、フラッシュメモリ、ＥＰＲＯＭ（Erasable Programmable Read Only Memory）、ＥＥＰＲＯＭ（Electrically Erasable Programmable Read Only Memory））を備えてもよい。例えば、演算部１３は、ＭＣＵなどのマイクロコントローラであってもよい。
【0132】
演算部１３は、上記のように求めた前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件を満たすために必要な反応温度を示す情報を出力部１４に出力してもよい。
【0133】
出力部１４は、計算機本体１２（詳しくは、演算部１３）が出力した計算結果（反応温度）を受信し、受信した計算結果を外部に出力するものである。本実施形態の出力部１４は、ＣＲＴディスプレイや液晶ディスプレイ、ＰＤＰ等の表示部によって構成されているが、これに限定されず、プリンタ等の印刷部や、他の装置（例えば、減圧蒸留軽油を含む原料油の水素化処理反応の制御等に用いられる計算機等）等へ出力するように構成されてもよい。また、出力部１４は、これらを組み合わせたものでもよい。例えば、出力部１４は、上述した情報出力ステップを実行する。
【0134】
また、本実施形態においては、コンピュータを反応温度算出装置として機能させるための反応温度算出プログラム及び当該プログラムを記憶したコンピュータの非一時的可読記録媒体が提供される。コンピュータの非一時的可読記録媒体としては、例えば、磁気テープ（デジタルデータストレージ（ＤＳＳ）など）、磁気ディスク（ハードディスクドライブ（ＨＤＤ）、フレキシブルディスク（ＦＤ）など）、光ディスク（コンパクトディスク（ＣＤ）、デジタルバーサタイルディスク（ＤＶＤ）、ブルーレイディスク（ＢＤ）など）、光磁気ディスク（ＭＯ）、フラッシュメモリ（ＳＳＤ（Solid State Drive）、メモリーカード、ＵＳＢメモリなど）が挙げられる。
【0135】
＜情報処理方法及び反応温度算出装置の活用方法＞
本実施形態の情報処理方法及び反応温度算出装置によると、減圧蒸留軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる。本実施形態の情報処理方法及び反応温度算出装置によると、反応温度の推定値の経時的なプロットを得ることができる。前記プロットと、機器使用最高温度等との関係から、以下のような活用法が考えられる。
本実施形態の情報処理方法及び反応温度算出装置によると、減圧蒸留軽油を含む原料油の水素化処理反応に関して、所定の反応条件を達成するために必要な反応温度を推定することができる。本実施形態の情報処理方法及び反応温度算出装置によると、反応温度の推定値の経時的なプロットを得ることができる。前記プロットと、機器使用最高温度等との関係から、以下のような活用法が考えられる。
【0136】
第１の活用方法としては、所定の反応条件における触媒を交換するのに必要な運転時間（日数）を推定することである。すなわち、前記プロットにより、機器使用最高温度となる温度となる時間を推定することができ、当該時間から触媒を交換するのに必要な運転時間（日数）を推定することができる。また、推定された反応温度が機器設定温度以下でない場合、図２のＳ４－１～４－３で示される情報処理方法が行われてもよい。図２のＳ４－１～４－３で示される情報処理方法は、目的とする反応温度に関する情報を取得する情報取得ステップ（図２のＳ４－１）と、取得した目的とする反応温度に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップ（図２のＳ４－２）と、劣化度に基づき、目的とする反応温度を満たすために必要な、原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を算出する反応条件算出ステップ（図２のＳ４－３）と、を含む。このようにして得られた前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報を出力する情報出力ステップ（図２のＳ４－４）をさらに有してもよい。具体的には、例えば、劣化関数１を使用する場合、前記式１９のＴｔに目的とする反応温度を代入し、Φを求める。得られたΦを前記式２に代入し、Ｄを求める。得られたＤを前記式３に代入して、前記式３の等式が成り立つような、ＳＰ、ＳＦ、ＬＨＳＶ、Ｐの組み合わせを求めればよい。また、上記組み合わせにおいて、例えばＳＰ、ＳＦを設定値として、Ｐ及びＬＨＳＶを求めてもよい。また、劣化関数としては、劣化関数１、劣化関数１－１、劣化関数２、劣化関数２－１のいずれを使用してもよい。なお、以上の各ステップは、例えば、本実施形態の反応温度算出装置１によって実行される。例えば、Ｓ４－１は取得部１１によって実行され、Ｓ４－２及びＳ４－３は計算機本体１２における演算部１３によって実行され、Ｓ４－４は出力部１４によって実行される。なお、目的温度は予め計算機本体１２における演算部１３に記憶させておいてもよい。
【0137】
第２の活用方法としては、所定の運転時間を達成するための反応条件（処理量（ＬＨＳＶ）等）を推定することである。このような情報処理方法は、図３のＳ１Ａ～３Ａで表される。図３のＳ１Ａ～３Ａで示される情報処理方法は、所定の運転時間における目的とする反応温度に関する情報を取得する情報取得ステップ（図３のＳ１Ａ）と、取得した目的とする反応温度に基づき、劣化関数により、触媒の劣化度を算出する劣化度算出ステップ（図３のＳ２Ａ）と、劣化度に基づき、目的とする反応温度を満たすために必要な、原料油に関する情報、生成油に関する情報、及び運転条件に関する情報を算出する反応条件算出ステップ（図３のＳ３Ａ）と、を含む。このようにして得られた前記原料油に関する情報、前記生成油に関する情報、及び前記運転条件に関する情報を出力する情報出力ステップ（図３のＳ４Ａ）をさらに有してもよい。具体的には、例えば、劣化関数１を使用する場合、前記式１９のＴｔに所定の運転時間ｔ日経過時における目的とする反応温度を代入し、Φを求める。得られたΦを前記式２に代入し、Ｄを求める。得られたＤを前記式３に代入して、前記式３の等式が成り立つような、ＳＰ、ＳＦ、ＬＨＳＶ、Ｐの組み合わせを求めればよい。また、上記組み合わせにおいて、例えばＳＰ、ＳＦを設定値として、Ｐ及びＬＨＳＶを求めてもよい。また、劣化関数としては、劣化関数１、劣化関数１－１、劣化関数２、劣化関数２－１のいずれを使用してもよい。なお、以上の各ステップは、例えば、本実施形態の反応温度算出装置１によって実行される。例えば、Ｓ１Ａは取得部１１によって実行され、Ｓ２Ａ及びＳ３Ａは計算機本体１２における演算部１３によって実行され、Ｓ４Ａは出力部１４によって実行される。なお、所定の運転時間における目的とする温度は予め計算機本体１２における演算部１３に記憶させておいてもよい。
【実施例0138】
以下、実施例により本発明をさらに具体的に説明するが、本発明は以下の実施例に限定されるものではない。
【0139】
［実施例］
ベンチスケールで、減圧蒸留軽油を７０体積％含む原料油を水素化処理触媒に接触処理させることにより水素化処理反応を行った。得られた結果を基に、劣化関数を算出した。なお、本実施例では、前記式１２で表される劣化関数２を使用した。前記式１２～１４におけるパラメータを上述の方法により求めた所、ｋ１＝０．９、ｋ２＝０．１、α１＝０．００１、α２＝０．０００５、ＰＢ＝４．９（ＭＰａ）、ａ＝０．６、Ｅｃ＝１５０（ｋＪ／ｍｏｌ）、ＴＢ＝６７３（Ｋ）、ＴＳＯＲ＝６３０（Ｋ）、ｎ＝１．５であった。
ベンチスケールで、減圧蒸留軽油を７０体積％含む原料油を水素化処理触媒に接触処理させることにより水素化処理反応を行った。得られた結果を基に、劣化関数を算出した。なお、本実施例では、前記式１２で表される劣化関数２を使用した。前記式１２～１４におけるパラメータを上述の方法により求めた所、ｋ１＝０．９、ｋ２＝０．１、α１＝０．００１、α２＝０．０００５、ＰＢ＝４．９（ＭＰａ）、ａ＝０．６、Ｅｃ＝１５０（ｋＪ／ｍｏｌ）、ＴＢ＝６７３（Ｋ）、ＴＳＯＲ＝６３０（Ｋ）、ｎ＝１．５であった。
【0140】
ベンチスケールで行った水素化処理反応と同じ原料油、水素化処理触媒を用いて、実機で反応を行った。前記式１３、１４における運転条件として、ＳＦ＝２．６２（質量％）、Ｐ＝６．０（ＭＰａ）、ＬＨＳＶ＝０．４８５（ｈ－１）、ＳＰ＝０．０７（質量％）とした。
また、事前に上述の方法で、脱硫の活性化エネルギーを求めた所、Ｅａ＝１２５（ｋＪ／ｍｏｌ）であった。
これらの基本劣化パラメータ及び反応条件を前記式１３、１４に代入し、任意の反応ｔ日経過時のＤ１、Ｄ２を得た。得られたＤ１、Ｄ２及び反応経過日数ｔを前記式１２に代入しΦを得た。得られたΦ、Ｅａ、ＴＳＯＲを前記式２０に代入して、任意の反応ｔ日経過時の要求温度Ｔｔを求めた。表１にｔ＝２５１日、５１４日、１１６４日の反応温度の実測値、及び上記方法で求めた要求温度Ｔｔ、及び要求温度Ｔｔ／反応温度の実測値の割合を示す。
また、事前に上述の方法で、脱硫の活性化エネルギーを求めた所、Ｅａ＝１２５（ｋＪ／ｍｏｌ）であった。
これらの基本劣化パラメータ及び反応条件を前記式１３、１４に代入し、任意の反応ｔ日経過時のＤ１、Ｄ２を得た。得られたＤ１、Ｄ２及び反応経過日数ｔを前記式１２に代入しΦを得た。得られたΦ、Ｅａ、ＴＳＯＲを前記式２０に代入して、任意の反応ｔ日経過時の要求温度Ｔｔを求めた。表１にｔ＝２５１日、５１４日、１１６４日の反応温度の実測値、及び上記方法で求めた要求温度Ｔｔ、及び要求温度Ｔｔ／反応温度の実測値の割合を示す。
【0141】
【表1】
【0142】
表１に示された通り、本発明によって求められた要求温度Ｔｔは、反応温度の実測値とほぼ同等となることがわかった。
【符号の説明】
【0143】
１・・・反応温度算出装置
１１・・・取得部
１２・・・計算機本体
１３・・・演算部
１４・・・出力部
１１・・・取得部
１２・・・計算機本体
１３・・・演算部
１４・・・出力部
【図1】
【図2】
【図3】
【図4】