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