JP2002372507A - 固定床触媒の寿命予測方法及び寿命予測プログラム - Google Patents
固定床触媒の寿命予測方法及び寿命予測プログラムInfo
- Publication number
- JP2002372507A JP2002372507A JP2001179615A JP2001179615A JP2002372507A JP 2002372507 A JP2002372507 A JP 2002372507A JP 2001179615 A JP2001179615 A JP 2001179615A JP 2001179615 A JP2001179615 A JP 2001179615A JP 2002372507 A JP2002372507 A JP 2002372507A
- Authority
- JP
- Japan
- Prior art keywords
- temperature distribution
- catalyst
- deterioration
- degree
- correlation coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Devices And Processes Conducted In The Presence Of Fluids And Solid Particles (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
(57)【要約】
【課題】 実現の困難な反応速度式モデル等の計算モデ
ルを構築すること無しに、精度のよい寿命予測を行うこ
とのできる固定床触媒の寿命予測方法及び寿命予測プロ
グラムを提供すること。 【解決手段】 反応器内の所定の箇所に固定床触媒が充
填された触媒充填層の温度分布を測定することにより前
記固定床触媒の寿命を予測する固定床触媒の寿命予測方
法であって、予め、触媒充填層の温度分布サンプルを劣
化度合いの進行程度ごとに多数取得して保存しておくス
テップと、寿命予測しようとする触媒充填層の温度分布
データを取得するステップと、温度分布データと多数の
温度分布サンプルとの相関係数を取得するステップと、
相関係数と劣化度合いとの組み合わせを取得するステッ
プと、組み合わせに基づいて相関係数が最大となる劣化
度合いを取得するステップとを有する。
ルを構築すること無しに、精度のよい寿命予測を行うこ
とのできる固定床触媒の寿命予測方法及び寿命予測プロ
グラムを提供すること。 【解決手段】 反応器内の所定の箇所に固定床触媒が充
填された触媒充填層の温度分布を測定することにより前
記固定床触媒の寿命を予測する固定床触媒の寿命予測方
法であって、予め、触媒充填層の温度分布サンプルを劣
化度合いの進行程度ごとに多数取得して保存しておくス
テップと、寿命予測しようとする触媒充填層の温度分布
データを取得するステップと、温度分布データと多数の
温度分布サンプルとの相関係数を取得するステップと、
相関係数と劣化度合いとの組み合わせを取得するステッ
プと、組み合わせに基づいて相関係数が最大となる劣化
度合いを取得するステップとを有する。
Description
【０００１】
【発明の属する技術分野】本発明は、反応器内の所定の
箇所に固定床触媒が充填された触媒充填層の温度分布を
測定することにより前記固定床触媒の寿命を予測する固
定床触媒の寿命予測方法及び寿命予測プログラムに関す
るものである。
箇所に固定床触媒が充填された触媒充填層の温度分布を
測定することにより前記固定床触媒の寿命を予測する固
定床触媒の寿命予測方法及び寿命予測プログラムに関す
るものである。
【０００２】
【従来の技術】固定床触媒は、化学産業において種々の
用途に広く用いられている。例えば、反応器内の所定の
箇所には、固定床触媒が充填された触媒充填層が設けら
れている。この反応器の入口から原料ガスを導入する
と、出口から反応ガスが取り出される。原料ガスを導入
することにより、固定床触媒上で化学反応が起こり、反
応に伴う発熱や吸熱が起こる。一般に触媒は使用するに
つれて、その活性を失う。経済的に成り立たない度合い
まで活性が下がると、 寿命が来たものとして、入れ替
えや再生が行われる。かかる産業分野において、固定床
触媒の寿命を予測することは重要である。その理由は、
製品在庫等の生産管理をバランスよく行うことができる
ようにするためである。そこで、固定床触媒の寿命予測
が行われており、触媒充填層の温度分布を測定して観察
することで予測をすることが試みられている。先ほど説
明したように原料ガスの導入により、固定床触媒上での
化学反応による発熱や吸熱現象が生じ、触媒充填層の深
さ方向に沿った温度分布は、劣化度合いが進行するにつ
れて変化してくる。したがって、この温度分布を測定し
て観測することにより、寿命予測を行うことができると
考えられる。
用途に広く用いられている。例えば、反応器内の所定の
箇所には、固定床触媒が充填された触媒充填層が設けら
れている。この反応器の入口から原料ガスを導入する
と、出口から反応ガスが取り出される。原料ガスを導入
することにより、固定床触媒上で化学反応が起こり、反
応に伴う発熱や吸熱が起こる。一般に触媒は使用するに
つれて、その活性を失う。経済的に成り立たない度合い
まで活性が下がると、 寿命が来たものとして、入れ替
えや再生が行われる。かかる産業分野において、固定床
触媒の寿命を予測することは重要である。その理由は、
製品在庫等の生産管理をバランスよく行うことができる
ようにするためである。そこで、固定床触媒の寿命予測
が行われており、触媒充填層の温度分布を測定して観察
することで予測をすることが試みられている。先ほど説
明したように原料ガスの導入により、固定床触媒上での
化学反応による発熱や吸熱現象が生じ、触媒充填層の深
さ方向に沿った温度分布は、劣化度合いが進行するにつ
れて変化してくる。したがって、この温度分布を測定し
て観測することにより、寿命予測を行うことができると
考えられる。
【０００３】しかし、定量的な評価が難しいため、寿命
予測の精度がよくないという問題がある。そこで、より
精度のよい寿命予測を行うために、反応速度式モデル、
触媒劣化速度式モデル、及び、伝熱計算モデルを構築し
た後に、シミュレーション計算により温度分布を計算
し、これと実際に測定された温度分布と比較することで
精度のよい予測が可能であると考えられる。
予測の精度がよくないという問題がある。そこで、より
精度のよい寿命予測を行うために、反応速度式モデル、
触媒劣化速度式モデル、及び、伝熱計算モデルを構築し
た後に、シミュレーション計算により温度分布を計算
し、これと実際に測定された温度分布と比較することで
精度のよい予測が可能であると考えられる。
【０００４】
【発明が解決しようとする課題】しかしながら、上記の
ような計算モデルを構築するためには、数多くの劣化試
験や化工計算により、計算用パラメータを特定する必要
がある。これに要する時間と労力は、寿命予測を行う上
で大きな妨げとなる。
ような計算モデルを構築するためには、数多くの劣化試
験や化工計算により、計算用パラメータを特定する必要
がある。これに要する時間と労力は、寿命予測を行う上
で大きな妨げとなる。
【０００５】本発明は上記実情に鑑みてなされたもので
あり、その課題は、実現の困難な反応速度式モデル等の
計算モデルを構築すること無しに、精度のよい寿命予測
を行うことのできる固定床触媒の寿命予測方法及び寿命
予測プログラムを提供することである。
あり、その課題は、実現の困難な反応速度式モデル等の
計算モデルを構築すること無しに、精度のよい寿命予測
を行うことのできる固定床触媒の寿命予測方法及び寿命
予測プログラムを提供することである。
【０００６】
【課題を解決するための手段】上記課題を解決するため
本発明に係る固定床触媒の寿命予測方法は、反応器内の
所定の箇所に固定床触媒が充填された触媒充填層の温度
分布を測定することにより前記固定床触媒の寿命を予測
する固定床触媒の寿命予測方法であって、予め、前記触
媒充填層の温度分布サンプルを劣化度合いの進行程度ご
とに多数取得して保存しておくステップと、寿命予測し
ようとする触媒充填層の温度分布データを取得するステ
ップと、前記温度分布データと多数の前記温度分布サン
プルとの相関係数を取得するステップと、前記相関係数
と前記劣化度合いとの組み合わせを取得するステップ
と、前記組み合わせに基づいて前記相関係数が最大とな
る前記劣化度合いを取得するステップとを有することを
特徴とするものである。
本発明に係る固定床触媒の寿命予測方法は、反応器内の
所定の箇所に固定床触媒が充填された触媒充填層の温度
分布を測定することにより前記固定床触媒の寿命を予測
する固定床触媒の寿命予測方法であって、予め、前記触
媒充填層の温度分布サンプルを劣化度合いの進行程度ご
とに多数取得して保存しておくステップと、寿命予測し
ようとする触媒充填層の温度分布データを取得するステ
ップと、前記温度分布データと多数の前記温度分布サン
プルとの相関係数を取得するステップと、前記相関係数
と前記劣化度合いとの組み合わせを取得するステップ
と、前記組み合わせに基づいて前記相関係数が最大とな
る前記劣化度合いを取得するステップとを有することを
特徴とするものである。
【０００７】この構成による固定床触媒の寿命予測方法
の作用・効果を、本発明の原理を説明する図１も利用し
ながら以下説明する。
の作用・効果を、本発明の原理を説明する図１も利用し
ながら以下説明する。
【０００８】（１）まず、あらかじめ、触媒充填層の温
度分布サンプルを劣化度合いの進行程度ごとに取得して
おく。図１の例では、１０％劣化のサンプルａ、４０％
劣化のサンプルｂ、６０％劣化のサンプルｃ、９０％劣
化のサンプルｄの４つの（「多数の」に相当する。）温
度分布サンプルが取得されている。このサンプルは、過
去の運転実績から抽出するか、予め実験をすることで得
ることができる。このサンプルは、横軸が温度を示し、
縦軸が触媒充填層の深さ（測定位置）を示すグラフとな
っている。
度分布サンプルを劣化度合いの進行程度ごとに取得して
おく。図１の例では、１０％劣化のサンプルａ、４０％
劣化のサンプルｂ、６０％劣化のサンプルｃ、９０％劣
化のサンプルｄの４つの（「多数の」に相当する。）温
度分布サンプルが取得されている。このサンプルは、過
去の運転実績から抽出するか、予め実験をすることで得
ることができる。このサンプルは、横軸が温度を示し、
縦軸が触媒充填層の深さ（測定位置）を示すグラフとな
っている。
【０００９】（２）次に、寿命予測をしようとする触媒
充填層の温度分布データを測定して取得する。これを図
１にＸで示す。温度分布サンプルも温度分布データも実
体は同じものであり、識別するための用語として「サン
プル」「データ」を付加する。
充填層の温度分布データを測定して取得する。これを図
１にＸで示す。温度分布サンプルも温度分布データも実
体は同じものであり、識別するための用語として「サン
プル」「データ」を付加する。
【００１０】（３）次に、温度分布データと多数の温度
分布サンプルとの相関係数を取得する。相関係数とは、
両グラフの近似度の高さを表わす指標であり、相関係数
が高いほど両グラフはよく近似していることを示す。図
１の例では、１０％劣化のサンプルとの相関係数ｒ2 ＝
０．４、４０％劣化のサンプルではｒ2 ＝０．９、６０
％劣化のサンプルではｒ2 ＝０．８、９０％劣化のサン
プルではｒ2 ＝０．２となっている。なお、相関係数は
公知の数学的手法により求められる。すなわち、ｘとｙ
の相関係数は、ｘ，ｙの共分散をｘの標準偏差とｙの標
準偏差の積で割ったものをいう。
分布サンプルとの相関係数を取得する。相関係数とは、
両グラフの近似度の高さを表わす指標であり、相関係数
が高いほど両グラフはよく近似していることを示す。図
１の例では、１０％劣化のサンプルとの相関係数ｒ2 ＝
０．４、４０％劣化のサンプルではｒ2 ＝０．９、６０
％劣化のサンプルではｒ2 ＝０．８、９０％劣化のサン
プルではｒ2 ＝０．２となっている。なお、相関係数は
公知の数学的手法により求められる。すなわち、ｘとｙ
の相関係数は、ｘ，ｙの共分散をｘの標準偏差とｙの標
準偏差の積で割ったものをいう。
【００１１】（４）次に、相関係数と劣化度合いとの組
み合わせを取得する。図１の例では、（１０％劣化：ｒ
2 ＝０．４）（４０％劣化：ｒ2 ＝０．９）（６０％劣
化：ｒ 2 ＝０．８）（９０％劣化：ｒ2 ＝０．２）であ
り、この組み合わせを図１に示すようなグラフに示すこ
ともできる。
み合わせを取得する。図１の例では、（１０％劣化：ｒ
2 ＝０．４）（４０％劣化：ｒ2 ＝０．９）（６０％劣
化：ｒ 2 ＝０．８）（９０％劣化：ｒ2 ＝０．２）であ
り、この組み合わせを図１に示すようなグラフに示すこ
ともできる。
【００１２】（５）次に、組み合わせに基づいて相関係
数が最大となる劣化度合いを取得する。例えば、上記求
められた組み合わせのうち（４０％劣化：ｒ2 ＝０．
９）の場合が最も相関係数が高い。従って、現在の劣化
度合いは４０％（又は、０．４）であると判断してもよ
い。しかし、より精度のよい劣化度合いを取得するに
は、次のような処理を実行する。
数が最大となる劣化度合いを取得する。例えば、上記求
められた組み合わせのうち（４０％劣化：ｒ2 ＝０．
９）の場合が最も相関係数が高い。従って、現在の劣化
度合いは４０％（又は、０．４）であると判断してもよ
い。しかし、より精度のよい劣化度合いを取得するに
は、次のような処理を実行する。
【００１３】すなわち、 上記取得された組み合わせ
は、グラフ上でプロットすることができ、このプロット
された各点を通る曲線（近似曲線）を求める。そして、
この曲線を微分することにより、相関係数が極値（最大
値）を取る場合の劣化度合いを演算により取得する。こ
の劣化度合いは、図１のグラフではおよそ４５％であ
り、より精度のよい劣化度合いを取得することが可能で
ある。
は、グラフ上でプロットすることができ、このプロット
された各点を通る曲線（近似曲線）を求める。そして、
この曲線を微分することにより、相関係数が極値（最大
値）を取る場合の劣化度合いを演算により取得する。こ
の劣化度合いは、図１のグラフではおよそ４５％であ
り、より精度のよい劣化度合いを取得することが可能で
ある。
【００１４】さて、以上の構成によれば、あらかじめ温
度分布サンプルを取得して（記録媒体等に）保存してお
けばよい。この温度分布サンプルと温度分布データとの
相関係数を求めることにより、劣化度合いを判断する手
法であるから、実現の困難な反応速度式モデル等の計算
モデルを構築する必要はない。これにより、精度のよい
寿命予測を行うことのできる固定床触媒の寿命予測方法
を提供することができる。
度分布サンプルを取得して（記録媒体等に）保存してお
けばよい。この温度分布サンプルと温度分布データとの
相関係数を求めることにより、劣化度合いを判断する手
法であるから、実現の困難な反応速度式モデル等の計算
モデルを構築する必要はない。これにより、精度のよい
寿命予測を行うことのできる固定床触媒の寿命予測方法
を提供することができる。
【００１５】本発明の好適な実施形態として、前記温度
分布データを取得するステップと、前記相関係数を取得
するステップと、前記組み合わせを取得するステップ
と、前記劣化度合いを取得するステップとを繰り返し実
行ことにより寿命を予測するステップを有するものがあ
げられる。
分布データを取得するステップと、前記相関係数を取得
するステップと、前記組み合わせを取得するステップ
と、前記劣化度合いを取得するステップとを繰り返し実
行ことにより寿命を予測するステップを有するものがあ
げられる。
【００１６】すでに説明したステップ（２）〜（５）を
繰り返して取得する。例えば、原料ガスの反応量が所定
値になるごとに、劣化度合いを取得する。原料ガスの累
計反応量が増えていけば、劣化度合いは進行する。よっ
て、上記ステップを繰り返して、劣化度合いを取得する
ことにより、劣化度合いが１００％（寿命時期）となる
時点を予測することが可能である。
繰り返して取得する。例えば、原料ガスの反応量が所定
値になるごとに、劣化度合いを取得する。原料ガスの累
計反応量が増えていけば、劣化度合いは進行する。よっ
て、上記ステップを繰り返して、劣化度合いを取得する
ことにより、劣化度合いが１００％（寿命時期）となる
時点を予測することが可能である。
【００１７】本発明の別の好適な実施形態として、前記
温度分布サンプルは、操業条件ごとに求めておくことも
のがあげられる。
温度分布サンプルは、操業条件ごとに求めておくことも
のがあげられる。
【００１８】操業条件、例えば、原料ガスの導入速度や
導入量が異なると、触媒充填層の温度分布も異なってく
ると考えられる。そこで、操業条件ごとに温度分布サン
プルを求めておくことにより、より精度の高い寿命予測
を行うことができる。
導入量が異なると、触媒充填層の温度分布も異なってく
ると考えられる。そこで、操業条件ごとに温度分布サン
プルを求めておくことにより、より精度の高い寿命予測
を行うことができる。
【００１９】上記課題を解決するため本発明に係る固定
床触媒の寿命予測プログラムは、反応器内の所定の箇所
に固定床触媒が充填された触媒充填層の温度分布を測定
することにより前記固定床触媒の寿命を予測する固定床
触媒の寿命予測プログラムであって、予め、前記触媒充
填層の温度分布サンプルを劣化度合いの進行程度ごとに
多数取得して保存しておく処理と、寿命予測しようとす
る触媒充填層の温度分布データを取得する処理と、前記
温度分布データと多数の前記温度分布サンプルとの相関
係数を取得する処理と、前記相関係数と前記劣化度合い
との組み合わせを取得する処理と、前記組み合わせに基
づいて前記相関係数が最大となる前記劣化度合いを取得
する処理とをコンピュータに実行させることを特徴とす
るものである。
床触媒の寿命予測プログラムは、反応器内の所定の箇所
に固定床触媒が充填された触媒充填層の温度分布を測定
することにより前記固定床触媒の寿命を予測する固定床
触媒の寿命予測プログラムであって、予め、前記触媒充
填層の温度分布サンプルを劣化度合いの進行程度ごとに
多数取得して保存しておく処理と、寿命予測しようとす
る触媒充填層の温度分布データを取得する処理と、前記
温度分布データと多数の前記温度分布サンプルとの相関
係数を取得する処理と、前記相関係数と前記劣化度合い
との組み合わせを取得する処理と、前記組み合わせに基
づいて前記相関係数が最大となる前記劣化度合いを取得
する処理とをコンピュータに実行させることを特徴とす
るものである。
【００２０】かかるプログラムをコンピュータにインス
トールすることにより、コンピュータに本発明による寿
命予測方法を実行させることができる。なお、プログラ
ムのコンピュータへのインストールは、適宜の方法によ
り行うことができる。
トールすることにより、コンピュータに本発明による寿
命予測方法を実行させることができる。なお、プログラ
ムのコンピュータへのインストールは、適宜の方法によ
り行うことができる。
【００２１】
【発明の実施の形態】本発明に係る固定床触媒の寿命予
測方法の好適な実施形態を図面を用いて説明する。
測方法の好適な実施形態を図面を用いて説明する。
【００２２】図２は、固定床触媒が充填された触媒充填
層を有する反応器Ａを示す図である。反応器Ａの上部に
は、原料ガスを導入する導入部１が、下部には反応ガス
が取り出される導出部２が設けられている。反応器Ａの
上下方向のほぼ中央部には、触媒充填層３が設けられ
る。触媒充填層３には、多数（数千本オーダー）のチュ
ーブが上下方向に沿って配置されている。このチューブ
の内部に固定床触媒が充填されている。
層を有する反応器Ａを示す図である。反応器Ａの上部に
は、原料ガスを導入する導入部１が、下部には反応ガス
が取り出される導出部２が設けられている。反応器Ａの
上下方向のほぼ中央部には、触媒充填層３が設けられ
る。触媒充填層３には、多数（数千本オーダー）のチュ
ーブが上下方向に沿って配置されている。このチューブ
の内部に固定床触媒が充填されている。
【００２３】導入部１から原料ガスが導入されると、固
定床触媒上での化学反応により、導出部２から反応ガス
が取り出される。上記固定床触媒上での化学反応に伴
い、発熱や吸熱現象が起こる。この発熱や吸熱に伴っ
て、触媒充填層には温度分布が生じる。その温度分布の
例を図２（ｂ）のグラフに示す。このグラフにおいて、
横軸は温度を表わし、縦軸は触媒充填層の上下（深さ）
方向の座標（すなわち、温度測定位置）を示す。固定床
触媒を充填した直後は、触媒充填層の上層の方で主に反
応が起こるために、（ｂ）のグラフに示すように上層の
方に温度分布のピークが見られる。そして、固定床触媒
を使用しているうちに、この温度分布のピークは徐々に
下方に下がってくる。すなわち、 劣化度合いが進行し
ていくにつれて、温度分布のピークは徐々に下方に移動
していく傾向にある。このように、温度分布曲線から劣
化度合いを知ることが可能である。
定床触媒上での化学反応により、導出部２から反応ガス
が取り出される。上記固定床触媒上での化学反応に伴
い、発熱や吸熱現象が起こる。この発熱や吸熱に伴っ
て、触媒充填層には温度分布が生じる。その温度分布の
例を図２（ｂ）のグラフに示す。このグラフにおいて、
横軸は温度を表わし、縦軸は触媒充填層の上下（深さ）
方向の座標（すなわち、温度測定位置）を示す。固定床
触媒を充填した直後は、触媒充填層の上層の方で主に反
応が起こるために、（ｂ）のグラフに示すように上層の
方に温度分布のピークが見られる。そして、固定床触媒
を使用しているうちに、この温度分布のピークは徐々に
下方に下がってくる。すなわち、 劣化度合いが進行し
ていくにつれて、温度分布のピークは徐々に下方に移動
していく傾向にある。このように、温度分布曲線から劣
化度合いを知ることが可能である。
【００２４】そこで、触媒充填層の温度分布を測定する
ために触媒充填層のなかに熱電対による温度計１８が設
けられている。本実施形態では、上下方向の１５箇所の
温度を測定できるようにしている。
ために触媒充填層のなかに熱電対による温度計１８が設
けられている。本実施形態では、上下方向の１５箇所の
温度を測定できるようにしている。
【００２５】＜固定床触媒の寿命予測装置＞図３は、本
発明に係る固定床触媒の寿命予測方法を実現するための
コンピュータ（寿命予測装置）の構成を示す図である。
発明に係る固定床触媒の寿命予測方法を実現するための
コンピュータ（寿命予測装置）の構成を示す図である。
【００２６】ＣＰＵ１０は、固定床触媒の寿命を予測す
るための寿命予測プログラムを実行するにあたり、各部
の作動制御を行う。モニター１１は、寿命予測を行うに
あたり、各種の演算結果や寿命予測結果をグラフ等を利
用してモニター画面に表示させる。キーボード１２及び
マウス１３は、各種データを入力したり、プログラムの
実行開始入力等のために用いられる。ＲＡＭ１４には、
寿命予測プログラムを実行するために、プログラムを格
納する。
るための寿命予測プログラムを実行するにあたり、各部
の作動制御を行う。モニター１１は、寿命予測を行うに
あたり、各種の演算結果や寿命予測結果をグラフ等を利
用してモニター画面に表示させる。キーボード１２及び
マウス１３は、各種データを入力したり、プログラムの
実行開始入力等のために用いられる。ＲＡＭ１４には、
寿命予測プログラムを実行するために、プログラムを格
納する。
【００２７】温度分布サンプルファイル１５には、寿命
予測を行おうとする固定床触媒の温度分布サンプルが多
数格納（保存）されている。このサンプルは劣化度合い
の進行程度ごとに格納される。本実施形態では、劣化度
合いが０．００，０．０１，０．０２，０．０３，０．
０５，０．１，０．２，０．３，０．４，０．５，０．
６，０．７，０．８，０．９，１．０の１５種類の温度
分布サンプルを予め測定しておく。
予測を行おうとする固定床触媒の温度分布サンプルが多
数格納（保存）されている。このサンプルは劣化度合い
の進行程度ごとに格納される。本実施形態では、劣化度
合いが０．００，０．０１，０．０２，０．０３，０．
０５，０．１，０．２，０．３，０．４，０．５，０．
６，０．７，０．８，０．９，１．０の１５種類の温度
分布サンプルを予め測定しておく。
【００２８】劣化度合いの定義について説明する。上記
の劣化度合いは、完全に劣化した状態を１．０（１００
％）で表現している。例えば、１００日で劣化するもの
と仮定すれば、１０日目は劣化度合い０．１（１０
％）、５０日目は０．５（５０％）である。なお、日に
ちで表現するのではなく、完全に劣化した状態の累計反
応量を１．０で表現してもよい。劣化度合いの定義は、
管理上の都合等を考慮して適宜に決めることができる。
の劣化度合いは、完全に劣化した状態を１．０（１００
％）で表現している。例えば、１００日で劣化するもの
と仮定すれば、１０日目は劣化度合い０．１（１０
％）、５０日目は０．５（５０％）である。なお、日に
ちで表現するのではなく、完全に劣化した状態の累計反
応量を１．０で表現してもよい。劣化度合いの定義は、
管理上の都合等を考慮して適宜に決めることができる。
【００２９】また、上記説明においては、劣化度合いが
初期の段階（上記例では０．１以下）のサンプルの取り
方は、劣化度合いが中期・末期の段階に比べてサンプル
の取り方に比べて、密度が高く（劣化度合いの間隔が小
さい）なっている。上記の例では、初期の段階では０．
０１間隔で温度分布サンプルを用意している。また、中
期・末期の段階では０．１間隔でサンプルを用意してい
る。これは、初期の段階の方が、中期・末期に比べて温
度分布の変化が大きいからである。このように、劣化度
合いの進行程度ごとに温度分布サンプルを取得しておく
場合に、どのような間隔で用意しておくかは、温度分布
の変化の度合いに応じて適切に設定することができる。
これにより、精度のよい寿命予測を行うことができる。
初期の段階（上記例では０．１以下）のサンプルの取り
方は、劣化度合いが中期・末期の段階に比べてサンプル
の取り方に比べて、密度が高く（劣化度合いの間隔が小
さい）なっている。上記の例では、初期の段階では０．
０１間隔で温度分布サンプルを用意している。また、中
期・末期の段階では０．１間隔でサンプルを用意してい
る。これは、初期の段階の方が、中期・末期に比べて温
度分布の変化が大きいからである。このように、劣化度
合いの進行程度ごとに温度分布サンプルを取得しておく
場合に、どのような間隔で用意しておくかは、温度分布
の変化の度合いに応じて適切に設定することができる。
これにより、精度のよい寿命予測を行うことができる。
【００３０】なお、固定床触媒の種類が異なれば温度分
布も異なってくるので、必要な種類の分だけ用意してお
く。また、操業条件が異なると、温度分布も異なってく
ると考えられるので、操業条件ごとに用意しておくこと
が好ましい。操業条件とは、例えば、原料ガスの導入速
度や導入量があげられる。
布も異なってくるので、必要な種類の分だけ用意してお
く。また、操業条件が異なると、温度分布も異なってく
ると考えられるので、操業条件ごとに用意しておくこと
が好ましい。操業条件とは、例えば、原料ガスの導入速
度や導入量があげられる。
【００３１】図３に戻り、温度分布データファイル１６
には、寿命予測を行おうとする触媒充填層の温度分布デ
ータが格納される。この温度分布データは、後述するよ
うに、温度分布サンプルと比較される。温度計１８は、
先ほど説明したように、反応器Ａの触媒充填層３の中に
挿入されて温度を測定する。この温度データは、インタ
ーフェイス１９を介してＡ−Ｄ変換された後に、所定の
ファイル形式で温度分布データファイル１６に格納され
る。
には、寿命予測を行おうとする触媒充填層の温度分布デ
ータが格納される。この温度分布データは、後述するよ
うに、温度分布サンプルと比較される。温度計１８は、
先ほど説明したように、反応器Ａの触媒充填層３の中に
挿入されて温度を測定する。この温度データは、インタ
ーフェイス１９を介してＡ−Ｄ変換された後に、所定の
ファイル形式で温度分布データファイル１６に格納され
る。
【００３２】寿命予測プログラム１７には、固定床触媒
の寿命予測を行うためのプログラムが格納されている。
このプログラムが有する機能の一部が図３に例示されて
いる。まず、相関係数取得手段１７ａは、温度分布サン
プルと温度分布データとを対比して相関係数を取得す
る。組み合わせ取得手段１７ｂは、相関係数と劣化度合
いとの組み合わせを取得する。劣化度合い取得手段１７
ｃは、相関係数と劣化度合いとの組み合わせのデータか
ら、現時点での劣化度合いを取得する。グラフ作成手段
１７ｄは、モニター１１に演算結果をグラフで表示させ
るためのグラフ化機能を備えている。
の寿命予測を行うためのプログラムが格納されている。
このプログラムが有する機能の一部が図３に例示されて
いる。まず、相関係数取得手段１７ａは、温度分布サン
プルと温度分布データとを対比して相関係数を取得す
る。組み合わせ取得手段１７ｂは、相関係数と劣化度合
いとの組み合わせを取得する。劣化度合い取得手段１７
ｃは、相関係数と劣化度合いとの組み合わせのデータか
ら、現時点での劣化度合いを取得する。グラフ作成手段
１７ｄは、モニター１１に演算結果をグラフで表示させ
るためのグラフ化機能を備えている。
【００３３】以上説明した、温度分布サンプルファイル
１５、温度分布データファイル１６、寿命予測プログラ
ム１７は、ハードディスク等の記録媒体に格納されてお
り、プログラムを実行するにあたり必要に応じてＲＡＭ
１４に格納される。また、図３にはデータバス２０が示
されている。
１５、温度分布データファイル１６、寿命予測プログラ
ム１７は、ハードディスク等の記録媒体に格納されてお
り、プログラムを実行するにあたり必要に応じてＲＡＭ
１４に格納される。また、図３にはデータバス２０が示
されている。
【００３４】＜寿命予測方法の手順＞次に、寿命予測を
行う場合の手順を説明する。図４は、寿命予測を行う場
合の手順を示すフローチャートである。
行う場合の手順を説明する。図４は、寿命予測を行う場
合の手順を示すフローチャートである。
【００３５】まずは、寿命予測を行おうとする固定床触
媒が充填された触媒充填層について、予め温度分布サン
プルを取得しておく（＃１）。この温度分布サンプルの
例を図５に示す。ここには、ある固定床触媒に関して４
通りの温度分布サンプルが示されている。各グラフにお
いて、横軸は温度データを取得した測定位置を示す。グ
ラフ中の原点（左端）は、触媒充填層の一番上の位置で
あり、右端は触媒充填層の一番下の位置である。グラフ
中に、○で示されるのが測定ポイントを示す。グラフの
縦軸は、温度（℃）を示す。図５において、（ａ）は劣
化度合い０．１（１０％）、（ｂ）は劣化度合い０．２
（２０％）、（ｃ）は劣化度合い０．５（５０％）、
（ｄ）は劣化度合い１．０（１００％）の場合のグラフ
を図示している。これらのデータは、温度分布サンプル
ファイル１５に格納されている。この図５からも分かる
ように、温度分布中に現れるピーク位置は、劣化度合い
が進行していくにつれて、右側に移行する。また、劣化
度合いが進行して中期・末期になると、ピーク位置が明
確に現れなくなるという傾向も示す。
媒が充填された触媒充填層について、予め温度分布サン
プルを取得しておく（＃１）。この温度分布サンプルの
例を図５に示す。ここには、ある固定床触媒に関して４
通りの温度分布サンプルが示されている。各グラフにお
いて、横軸は温度データを取得した測定位置を示す。グ
ラフ中の原点（左端）は、触媒充填層の一番上の位置で
あり、右端は触媒充填層の一番下の位置である。グラフ
中に、○で示されるのが測定ポイントを示す。グラフの
縦軸は、温度（℃）を示す。図５において、（ａ）は劣
化度合い０．１（１０％）、（ｂ）は劣化度合い０．２
（２０％）、（ｃ）は劣化度合い０．５（５０％）、
（ｄ）は劣化度合い１．０（１００％）の場合のグラフ
を図示している。これらのデータは、温度分布サンプル
ファイル１５に格納されている。この図５からも分かる
ように、温度分布中に現れるピーク位置は、劣化度合い
が進行していくにつれて、右側に移行する。また、劣化
度合いが進行して中期・末期になると、ピーク位置が明
確に現れなくなるという傾向も示す。
【００３６】次に、寿命予測を行おうとしている触媒充
填層の温度分布データを温度計１８により測定して、温
度分布データファイル１６に格納する（＃２）。この温
度分布データをグラフ化したものが図６に示される。こ
のグラフの縦軸及び横軸のパラメータは図５と同じであ
る。
填層の温度分布データを温度計１８により測定して、温
度分布データファイル１６に格納する（＃２）。この温
度分布データをグラフ化したものが図６に示される。こ
のグラフの縦軸及び横軸のパラメータは図５と同じであ
る。
【００３７】次に、取得された温度分布データと、予め
取得してある温度分布サンプルのそれぞれとの相関係数
を求める。温度分布サンプルと温度分布データとの類似
度が高ければ高いほど相関係数は大きくなる。相関係数
は、０〜１．０までの数値で表現することができ、相関
係数が１．０であれば、温度分布サンプルと温度分布デ
ータとは全く同じデータ分布（曲線）を示すことにな
る。なお、相関係数を求めるための演算式は周知であ
り、その説明は省略する（＃３）。
取得してある温度分布サンプルのそれぞれとの相関係数
を求める。温度分布サンプルと温度分布データとの類似
度が高ければ高いほど相関係数は大きくなる。相関係数
は、０〜１．０までの数値で表現することができ、相関
係数が１．０であれば、温度分布サンプルと温度分布デ
ータとは全く同じデータ分布（曲線）を示すことにな
る。なお、相関係数を求めるための演算式は周知であ
り、その説明は省略する（＃３）。
【００３８】次に、取得された相関係数と劣化度合いと
の組み合わせを取得する（＃４）。これを図７に示す。
図７（ａ）は表として組み合わせ結果を示したものであ
る。この表によると、例えば、劣化度合い０．０１の温
度分布サンプルと、図６の温度分布データとの相関係数
は、０．１９５となっている。図７（ｂ）は、（ａ）に
示される組み合わせ結果をグラフ化したものである。プ
ロットした点は、グラフ中において、□で示される。こ
のグラフの横軸は、劣化度合いを示し、縦軸は、相関係
数を示している。以上のように、組み合わせを表わすに
は、（ａ）のような表形式でもよいし、（ｂ）のような
グラフ形式でもよい。
の組み合わせを取得する（＃４）。これを図７に示す。
図７（ａ）は表として組み合わせ結果を示したものであ
る。この表によると、例えば、劣化度合い０．０１の温
度分布サンプルと、図６の温度分布データとの相関係数
は、０．１９５となっている。図７（ｂ）は、（ａ）に
示される組み合わせ結果をグラフ化したものである。プ
ロットした点は、グラフ中において、□で示される。こ
のグラフの横軸は、劣化度合いを示し、縦軸は、相関係
数を示している。以上のように、組み合わせを表わすに
は、（ａ）のような表形式でもよいし、（ｂ）のような
グラフ形式でもよい。
【００３９】次に、図７に示される組み合わせ結果に基
づいて、相関係数が最大となる場合の劣化度合を取得す
る（＃５）。図７（ａ）においては、相関係数が最も高
くなるのは劣化度合いが０．３の場合であり、このとき
の相関係数は０．９４５となっている。従って、この場
合の劣化度合いは０．３であると判断をしてもよく、こ
れも本発明の枠内に入る。しかし、以下のような手順を
行うことにより、より精度のよい劣化度合いを取得する
ことができる。
づいて、相関係数が最大となる場合の劣化度合を取得す
る（＃５）。図７（ａ）においては、相関係数が最も高
くなるのは劣化度合いが０．３の場合であり、このとき
の相関係数は０．９４５となっている。従って、この場
合の劣化度合いは０．３であると判断をしてもよく、こ
れも本発明の枠内に入る。しかし、以下のような手順を
行うことにより、より精度のよい劣化度合いを取得する
ことができる。
【００４０】まず、図７（ｂ）の各プロット点を多項式
により近似する。この図では１５のプロットがあるが、
これをすべて考慮すると、近似多項式が複雑化すると共
に、相関係数の低い部分はノイズ成分（この部分を考慮
しても精度は上がらない。）となる。そこで、近似式に
回帰する場合には、相関係数の低いプロットは考慮せ
ず、相関係数の高い部分のみを考慮して近似式に回帰す
る。
により近似する。この図では１５のプロットがあるが、
これをすべて考慮すると、近似多項式が複雑化すると共
に、相関係数の低い部分はノイズ成分（この部分を考慮
しても精度は上がらない。）となる。そこで、近似式に
回帰する場合には、相関係数の低いプロットは考慮せ
ず、相関係数の高い部分のみを考慮して近似式に回帰す
る。
【００４１】そこで、図７（ｂ）に黒い四角で示すよう
に、相関係数が最も高い点と、その点の前後の点の合計
３点を抽出し、この３点を通る多項式（二次曲線）を演
算して求める。この演算式は、 ｙ＝−３．７６９４ｘ2＋２．１０８４ｘ＋０．６５１
８ で示される。ただし、ｘは劣化度合い、ｙは相関係数を
示す。次に、この二次曲線を微分することで、相関係数
が極値（最大値）となる劣化度合いの値を求める。上記
式の場合、ｘ＝０．２８となる。この数値を元に、固定
床触媒の寿命予測を行うことができる。ただし、より精
度の高い寿命予測を行うためには、上記のステップ＃２
〜＃５を繰り返して行う。繰り返して温度分布データを
取得するときの、測定間隔は適宜設定することができ
る。
に、相関係数が最も高い点と、その点の前後の点の合計
３点を抽出し、この３点を通る多項式（二次曲線）を演
算して求める。この演算式は、 ｙ＝−３．７６９４ｘ2＋２．１０８４ｘ＋０．６５１
８ で示される。ただし、ｘは劣化度合い、ｙは相関係数を
示す。次に、この二次曲線を微分することで、相関係数
が極値（最大値）となる劣化度合いの値を求める。上記
式の場合、ｘ＝０．２８となる。この数値を元に、固定
床触媒の寿命予測を行うことができる。ただし、より精
度の高い寿命予測を行うためには、上記のステップ＃２
〜＃５を繰り返して行う。繰り返して温度分布データを
取得するときの、測定間隔は適宜設定することができ
る。
【００４２】図８は、固定床触媒を使用開始してからの
劣化度合いの推移を示すグラフである。上記ステップを
繰り返すことで、係るグラフを取得することができる
（＃６）。なお、このグラフにおいて、横軸は累計反応
量（ｋｇ）を示し、縦軸は劣化度合いを示す。図８から
も分かるように、累計反応量が増えていくほど劣化度合
いが進行していくことが理解される。そして、このグラ
フからも分かるように、累計反応量と劣化度合いとの関
係は、ほぼ直線関係にある。そこで、各プロット点を元
にして近似直線式を求めることができる。この直線か
ら、劣化度合いが１．０となる場合の累計反応量を演算
で求めることにより、固定床触媒の寿命予測を行うこと
ができる。
劣化度合いの推移を示すグラフである。上記ステップを
繰り返すことで、係るグラフを取得することができる
（＃６）。なお、このグラフにおいて、横軸は累計反応
量（ｋｇ）を示し、縦軸は劣化度合いを示す。図８から
も分かるように、累計反応量が増えていくほど劣化度合
いが進行していくことが理解される。そして、このグラ
フからも分かるように、累計反応量と劣化度合いとの関
係は、ほぼ直線関係にある。そこで、各プロット点を元
にして近似直線式を求めることができる。この直線か
ら、劣化度合いが１．０となる場合の累計反応量を演算
で求めることにより、固定床触媒の寿命予測を行うこと
ができる。
【００４３】＜別実施形態＞本発明が適用される固定床
触媒は特定の種類に限定されるものではなく、反応器の
構造についても特定のものに限定されるものではない。
温度分布サンプルをどの程度細かく取得しておくかにつ
いては、適宜設定できるものである。
触媒は特定の種類に限定されるものではなく、反応器の
構造についても特定のものに限定されるものではない。
温度分布サンプルをどの程度細かく取得しておくかにつ
いては、適宜設定できるものである。
【図１】本発明の原理を説明する図
【図２】反応器の構成及び温度分布を示す図
【図３】本発明に係る固定床触媒の寿命予測方法を実現
するためのコンピュータの構成を示す図
するためのコンピュータの構成を示す図
【図４】寿命予測を行う場合の手順を示すフローチャー
ト
ト
【図５】温度分布サンプルの例を示す図
【図６】温度分布データのグラフを示す図
【図７】相関係数と劣化度合いの組み合わせを取得した
結果を示す図
結果を示す図
【図８】固定床触媒を使用開始してからの劣化度合いの
推移を示すグラフ
推移を示すグラフ
Ａ 反応器 ３ 触媒充填層 １５ 温度分布サンプルファイル １６ 温度分布データファイル １７ 寿命予測プログラム １７ａ 相関係数取得手段 １７ｂ 組み合わせ取得手段 １７ｃ 劣化度合い取得手段 １８ 温度計
───────────────────────────────────────────────────── フロントページの続き Ｆターム(参考） 2G040 AB12 BB09 CA02 CB02 DA02 DA15 FA01 GA07 ZA05 2G050 BA01 CA02 DA01 EA05 EB10 EC01 4G070 AA01 AB04 BB01 CA01 CB30 CC03 DA30
Claims (4)
- 【請求項１】 反応器内の所定の箇所に固定床触媒が充
填された触媒充填層の温度分布を測定することにより前
記固定床触媒の寿命を予測する固定床触媒の寿命予測方
法であって、 予め、前記触媒充填層の温度分布サンプルを劣化度合い
の進行程度ごとに多数取得して保存しておくステップ
と、 寿命予測しようとする触媒充填層の温度分布データを取
得するステップと、 前記温度分布データと多数の前記温度分布サンプルとの
相関係数を取得するステップと、 前記相関係数と前記劣化度合いとの組み合わせを取得す
るステップと、 前記組み合わせに基づいて前記相関係数が最大となる前
記劣化度合いを取得するステップとを有することを特徴
とする固定床触媒の寿命予測方法。 - 【請求項２】 前記温度分布データを取得するステップ
と、前記相関係数を取得するステップと、前記組み合わ
せを取得するステップと、前記劣化度合いを取得するス
テップとを繰り返し実行ことにより寿命を予測するステ
ップを有することを特徴とする請求項１に記載の固定床
触媒の寿命予測方法。 - 【請求項３】 前記温度分布サンプルは、操業条件ごと
に求めておくことを特徴とする請求項１又は２に記載の
固定床触媒の寿命予測方法。 - 【請求項４】 反応器内の所定の箇所に固定床触媒が充
填された触媒充填層の温度分布を測定することにより前
記固定床触媒の寿命を予測する固定床触媒の寿命予測プ
ログラムであって、 予め、前記触媒充填層の温度分布サンプルを劣化度合い
の進行程度ごとに多数取得して保存しておく処理と、 寿命予測しようとする触媒充填層の温度分布データを取
得する処理と、 前記温度分布データと多数の前記温度分布サンプルとの
相関係数を取得する処理と、 前記相関係数と前記劣化度合いとの組み合わせを取得す
る処理と、 前記組み合わせに基づいて前記相関係数が最大となる前
記劣化度合いを取得する処理とをコンピュータに実行さ
せる特徴とする固定床触媒の寿命予測プログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001179615A JP3669573B2 (ja) | 2001-06-14 | 2001-06-14 | 固定床触媒の寿命予測方法及び寿命予測プログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001179615A JP3669573B2 (ja) | 2001-06-14 | 2001-06-14 | 固定床触媒の寿命予測方法及び寿命予測プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002372507A true JP2002372507A (ja) | 2002-12-26 |
| JP3669573B2 JP3669573B2 (ja) | 2005-07-06 |
Family
ID=19020160
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001179615A Expired - Fee Related JP3669573B2 (ja) | 2001-06-14 | 2001-06-14 | 固定床触媒の寿命予測方法及び寿命予測プログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3669573B2 (ja) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015208703A (ja) * | 2014-04-24 | 2015-11-24 | 東ソー株式会社 | 触媒寿命予測方法および触媒寿命解析装置 |
| JP2015208702A (ja) * | 2014-04-24 | 2015-11-24 | 東ソー株式会社 | 触媒性能の経時変化解析方法および触媒性能の経時変化解析装置 |
| JP2015208704A (ja) * | 2014-04-24 | 2015-11-24 | 東ソー株式会社 | 触媒寿命解析装置を備えた触媒反応製造システム |
| JP2016170117A (ja) * | 2015-03-13 | 2016-09-23 | 東ソー株式会社 | 触媒寿命予測方法と触媒寿命解析装置 |
| JP2016170116A (ja) * | 2015-03-13 | 2016-09-23 | 東ソー株式会社 | 触媒性能の経時変化解析方法および触媒性能の経時変化解析装置 |
| WO2018216746A1 (ja) * | 2017-05-25 | 2018-11-29 | コスモ石油株式会社 | プラントの推奨運転条件を提供する方法、サーバ、コンピュータ可読命令および記録媒体 |
| CN110210117A (zh) * | 2019-05-31 | 2019-09-06 | 西安工程大学 | 一种细纱机剩余运行寿命的预测方法及系统 |
| WO2019187962A1 (ja) * | 2018-03-28 | 2019-10-03 | コスモ石油株式会社 | Ｒｆ装置の推奨触媒を提案する装置、方法、プログラム、非一時的コンピュータ可読記録媒体 |
| WO2023048284A1 (ja) * | 2021-09-27 | 2023-03-30 | 日本化薬株式会社 | 多管式反応器の運転またはその準備行為をサポートする方法 |
-
2001
- 2001-06-14 JP JP2001179615A patent/JP3669573B2/ja not_active Expired - Fee Related
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015208703A (ja) * | 2014-04-24 | 2015-11-24 | 東ソー株式会社 | 触媒寿命予測方法および触媒寿命解析装置 |
| JP2015208702A (ja) * | 2014-04-24 | 2015-11-24 | 東ソー株式会社 | 触媒性能の経時変化解析方法および触媒性能の経時変化解析装置 |
| JP2015208704A (ja) * | 2014-04-24 | 2015-11-24 | 東ソー株式会社 | 触媒寿命解析装置を備えた触媒反応製造システム |
| JP2016170117A (ja) * | 2015-03-13 | 2016-09-23 | 東ソー株式会社 | 触媒寿命予測方法と触媒寿命解析装置 |
| JP2016170116A (ja) * | 2015-03-13 | 2016-09-23 | 東ソー株式会社 | 触媒性能の経時変化解析方法および触媒性能の経時変化解析装置 |
| US10915837B2 (en) | 2017-05-25 | 2021-02-09 | Cosmo Oil Co., Ltd. | Method, server, computer-readable command, and recording medium for providing recommended operation condition for plant |
| WO2018216746A1 (ja) * | 2017-05-25 | 2018-11-29 | コスモ石油株式会社 | プラントの推奨運転条件を提供する方法、サーバ、コンピュータ可読命令および記録媒体 |
| WO2019187962A1 (ja) * | 2018-03-28 | 2019-10-03 | コスモ石油株式会社 | Ｒｆ装置の推奨触媒を提案する装置、方法、プログラム、非一時的コンピュータ可読記録媒体 |
| JP2019175135A (ja) * | 2018-03-28 | 2019-10-10 | コスモ石油株式会社 | Ｒｆ装置の推奨触媒を提案する装置、方法、プログラム、非一時的コンピュータ可読記録媒体 |
| JP7169755B2 (ja) | 2018-03-28 | 2022-11-11 | コスモ石油株式会社 | Ｒｆ装置の推奨触媒を提案する装置、方法、プログラム、非一時的コンピュータ可読記録媒体 |
| CN110210117A (zh) * | 2019-05-31 | 2019-09-06 | 西安工程大学 | 一种细纱机剩余运行寿命的预测方法及系统 |
| CN110210117B (zh) * | 2019-05-31 | 2023-03-31 | 西安工程大学 | 一种细纱机剩余运行寿命的预测方法及系统 |
| WO2023048284A1 (ja) * | 2021-09-27 | 2023-03-30 | 日本化薬株式会社 | 多管式反応器の運転またはその準備行為をサポートする方法 |
| JP7345072B2 (ja) | 2021-09-27 | 2023-09-14 | 日本化薬株式会社 | 多管式反応器の運転またはその準備行為をサポートする方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3669573B2 (ja) | 2005-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zamba et al. | A multivariate change-point model for statistical process control | |
| JP5293739B2 (ja) | 重回帰分析による予測モデルの作成方法、作成システムおよび作成プログラム | |
| JP2004532476A (ja) | プロセス装置の監視プログラム | |
| JP2022500778A (ja) | 製造工程の生成物である化合物および／または化合物製剤の品質を予測するシステムおよび方法 | |
| JP2002372507A (ja) | 固定床触媒の寿命予測方法及び寿命予測プログラム | |
| US10431341B2 (en) | Detection device, method, and program for assisting network entropy-based detection of precursor to state transition of biological object | |
| US7908300B2 (en) | Guided entry system for individuals for annotating process deviations | |
| JP6024718B2 (ja) | 高炉炉熱予測装置及び高炉炉熱予測方法 | |
| JP2012178136A (ja) | 予測装置、予測方法、及び、プログラム | |
| JP6206368B2 (ja) | 高炉内状態推定装置および高炉内状態推定方法 | |
| WO2018081465A1 (en) | Systems and methods for characterizing nucleic acid in a biological sample | |
| US11934982B2 (en) | Feedstock processing corrosion management | |
| US10733331B2 (en) | Information processing apparatus, information processing method, and recording medium | |
| CN103472732A (zh) | 一种改进的基于马氏距离的多变量控制器性能监控方法 | |
| TWI802294B (zh) | 實驗點推薦裝置、實驗點推薦方法及半導體裝置製造系統 | |
| JP2007094856A (ja) | プラント運転支援装置 | |
| Fonseca-Pérez et al. | A reconsideration on the resolution of phase stability analysis using stochastic global optimization methods: Proposal of a reliable set of benchmark problems | |
| JP2022102726A (ja) | 化学反応における活性化エネルギーの学習装置及び予測装置 | |
| US20220383160A1 (en) | Sample analysis device, sample analysis method, pharmaceutical analysis device and pharmaceutical analysis method | |
| TWI321657B (ja) | ||
| JP6492814B2 (ja) | 触媒寿命解析装置を備えた触媒反応製造システム | |
| TWI787971B (zh) | 腐蝕率預測裝置、方法及其電腦程式產品 | |
| JP2008128736A (ja) | 波形表示装置及び波形表示方法 | |
| JP7453895B2 (ja) | 探索条件提示装置、探索条件提示方法、及び探索条件提示プログラム | |
| Kim et al. | Flexible Bayesian human fecundity models |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040906 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20041124 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050119 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050406 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050407 |
|
| R150 | Certificate of patent (=grant) or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090422 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090422 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100422 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110422 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120422 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130422 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |