特許 7238476 設備管理支援装置、設備管理支援方法、プログラム及びコンピュータ読み取り可能な記録媒体

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-03-06

(45)【発行日】2023-03-14

(54)【発明の名称】設備管理支援装置、設備管理支援方法、プログラム及びコンピュータ読み取り可能な記録媒体

(51)【国際特許分類】

   G05B 19/418 20060101AFI20230307BHJP

   G06Q 10/06 20230101ALI20230307BHJP

   C10B 29/06 20060101ALN20230307BHJP

【ＦＩ】

G05B19/418 Z

G06Q10/06

C10B29/06

【請求項の数】 9

(21)【出願番号】P 2019038324

(22)【出願日】2019-03-04

(65)【公開番号】P2020024657

(43)【公開日】2020-02-13

【審査請求日】2021-11-04

(31)【優先権主張番号】P 2018139162

(32)【優先日】2018-07-25

(33)【優先権主張国・地域又は機関】JP

(73)【特許権者】

【識別番号】000006655

【氏名又は名称】日本製鉄株式会社

(74)【代理人】

【識別番号】110002044

【氏名又は名称】弁理士法人ブライタス

(72)【発明者】

【氏名】阿部 快洋

(72)【発明者】

【氏名】枝野 寛弘

(72)【発明者】

【氏名】塩谷 政典

【審査官】石田 宏之

(56)【参考文献】

【文献】米国特許出願公開第２００９／００９１４４３（ＵＳ，Ａ１）

【文献】特開２０１３－２４６８１８（ＪＰ，Ａ）

【文献】特開２０１７－１１７４３８（ＪＰ，Ａ）

【文献】米国特許出願公開第２００２／０１６１７６３（ＵＳ，Ａ１）

【文献】特開２０１２－０５８８９０（ＪＰ，Ａ）

【文献】国際公開第２０１７／１３４９０８（ＷＯ，Ａ１）

【文献】特許第５８３８８９６（ＪＰ，Ｂ２）

【文献】特許第５７７２５４６（ＪＰ，Ｂ２）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０５Ｂ １９／４１８

Ｇ０６Ｑ １０／２０

Ｃ１０Ｂ ２９／０６

(57)【特許請求の範囲】

【請求項1】

コークス炉を構成する複数の窯のうち１つの窯を評価対象として、複数の異なる期間に取得された前記評価対象の窯の押出位置毎の押出負荷を含む押出負荷データを、前記コークス炉における押出位置及び押出負荷を変数とした多次元のヒストグラムに変換するヒストグラム変換処理部と、
複数の前記ヒストグラムの差異を表す評価指標を算出する評価指標算出部と、
押出負荷の最大値または押出負荷のピーク位置に基づいて決定される前記評価対象の窯の状態が悪化したのか良化したのかを表す符号を、前記評価指標に対して付与した値を、異常度として算出する異常度算出部と、
複数の前記窯のそれぞれに対して算出された前記異常度の高い順に、点検の必要な前記窯の優先度を決定する優先度決定部と、
前記優先度決定部により決定された優先度の高い窯を作業者に通知するための通知情報を生成する情報生成部と、
を備える、設備管理支援装置。

【請求項2】

前記評価指標算出部は、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅを用いて前記評価指標を算出する、請求項１に記載の設備管理支援装置。

【請求項3】

前記評価指標算出部は、χ２距離を用いて前記評価指標を算出する、請求項１に記載の設備管理支援装置。

【請求項4】

前記評価指標算出部は、ある期間の前記押出位置に対する前記押出負荷の分布を基づき算出される重みを考慮した重みつきχ２距離を用いて前記評価指標を算出する、請求項１に記載の設備管理支援装置。

【請求項5】

前記評価指標算出部は、ある期間の前記押出位置に対する前記押出負荷の分布を基づき算出されるＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅと、χ２距離または重みを考慮した重みつきχ２距離のうちいずれか一方とを組み合わせて、前記評価指標を算出する、請求項１に記載の設備管理支援装置。

【請求項6】

前記情報生成部は、前記窯のトラブル発生を検知するために予め設定された閾値に基づき、前記評価指標を時系列に並べた情報あるいは前記異常度を時系列に並べた情報に対して閾値判定を行い、トラブル発生を予測した情報を生成する、請求項１に記載の設備管理支援装置。

【請求項7】

コークス炉を構成する複数の窯のうち１つの窯を評価対象として、複数の異なる期間に取得された前記評価対象の窯の押出位置毎の押出負荷を含む押出負荷データを、前記コークス炉における押出位置及び押出負荷を変数とした多次元のヒストグラムに変換するヒストグラム変換処理ステップと、
複数の前記ヒストグラムの差異を表す評価指標を算出する評価指標算出ステップと、
押出負荷の最大値または押出負荷のピーク位置に基づいて決定される前記評価対象の窯の状態が悪化したのか良化したのかを表す符号を、前記評価指標に対して付与した値を、異常度として算出する異常度算出ステップと、
複数の前記窯のそれぞれに対して算出された前記異常度の高い順に、点検の必要な前記窯の優先度を決定する優先度決定ステップと、
前記優先度決定ステップにて決定された優先度の高い窯を作業者に通知するための通知情報を生成する情報生成ステップと、
を含む、設備管理支援方法。

【請求項8】

コンピュータを、
コークス炉を構成する複数の窯のうち１つの窯を評価対象として、複数の異なる期間に取得された前記評価対象の窯の押出位置毎の押出負荷を含む押出負荷データを、前記コークス炉における押出位置及び押出負荷を変数とした多次元のヒストグラムに変換するヒストグラム変換処理部と、
複数の前記ヒストグラムの差異を表す評価指標を算出する評価指標算出部と、
押出負荷の最大値または押出負荷のピーク位置に基づいて決定される前記評価対象の窯の状態が悪化したのか良化したのかを表す符号を、前記評価指標に対して付与した値を、異常度として算出する異常度算出部と、
複数の前記窯のそれぞれに対して算出された前記異常度の高い順に、点検の必要な前記窯の優先度を決定する優先度決定部と、
前記優先度決定部により決定された優先度の高い窯を作業者に通知するための通知情報を生成する情報生成部と、
を備える、設備管理支援装置として機能させるプログラム。

【請求項9】

コンピュータに、
コークス炉を構成する複数の窯のうち１つの窯を評価対象として、複数の異なる期間に取得された前記評価対象の窯の押出位置毎の押出負荷を含む押出負荷データを、前記コークス炉における押出位置及び押出負荷を変数とした多次元のヒストグラムに変換するヒストグラム変換処理部と、
複数の前記ヒストグラムの差異を表す評価指標を算出する評価指標算出部と、
押出負荷の最大値または押出負荷のピーク位置に基づいて決定される前記評価対象の窯の状態が悪化したのか良化したのかを表す符号を、前記評価指標に対して付与した値を、異常度として算出する異常度算出部と、
複数の前記窯のそれぞれに対して算出された前記異常度の高い順に、点検の必要な前記窯の優先度を決定する優先度決定部と、
前記優先度決定部により決定された優先度の高い窯を作業者に通知するための通知情報を生成する情報生成部と、
を備える、設備管理支援装置として機能させるプログラムを記録したコンピュータ読み取り可能な記録媒体。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、製造業で用いられる設備を構成する複数のユニットの状態管理を支援する設備管理支援装置、設備管理支援方法、プログラム及びコンピュータ読み取り可能な記録媒体に関する。

【背景技術】

【0002】

製造業で用いられる設備では、安定した操業を行うことができるように定期的にメンテナンスが実施される。定期的なメンテナンス以外にも設備に異常な傾向がみられるときには、追加のメンテナンスが実施される。

【0003】

例えば、製造業で用いられる設備の一例であるコークス炉では、炉体保全のため、メンテナンス予定に従い、炭化室診断補修装置あるいは窯幅計測等を用いて、炉内点検が実施される。例外的に、コークスの押出し時に異常な押出負荷波形が発生した場合には、炉壁の凹凸や肌荒れ、窯口狭小化等といった窯が異常な状態となっている虞があるため、予定以外の炉内点検やメンテナンスが追加で実施される。予定外の炉内点検やメンテナンスを実施するかどうかの判断は、例えばたかだか３サイクル分程度の範囲での押出負荷波形の最大値（最大押出負荷）に着目して行われる。

【0004】

ここで、特許文献１には、炭化室の異常を判定する方法が開示されている。特許文献１に記載の方法では、まず、コークス炉の炭化室毎に、生成したコークスを押出機の押出ラムで押し出す際、押出ラムの炭化室内での押出位置に対する押出ラムを駆動するモータの負荷電流を計測した負荷電流トレンドから、コークス炉の操業状態を考慮して理想負荷電流トレンドを予め求めておく。そして、測定対象となる炭化室におけるコークス押出完了後の実績負荷電流トレンドと理想負荷電流トレンドとの押出位置毎の負荷電流の偏差を求め、求めた偏差から当該炭化室の異常を判定する。

【先行技術文献】

【特許文献】

【0005】

【文献】特開２０１７－１７１７９１号公報

【非特許文献】

【0006】

【文献】Elizaveta Levina; Peter Bickel (2001). "The EarthMover's istance is the Mallows Distance: Some Insights from Statistics". Proceedings of ICCV 2001. Vancouver, Canada: 251‐256.

【発明の概要】

【発明が解決しようとする課題】

【0007】

炭化室の補修直後は、カーボンを取り除いた状態で負荷が増大する場合もあれば、炉壁平滑化がされて負荷が減少する場合もあり、押出負荷が落ち着くまではある程度の期間が必要である。しかし、その期間は窯によって様々であり、上記特許文献１に記載の方法のように、理想負荷トレンドを設定するためのデータ選定は容易ではない。また、カーボン付着状態や炉壁狭小化等によって押出性も変化していくため、上記特許文献１に記載の方法のように、理想負荷トレンドを１回設定するだけでは、炭化室の異常を高精度に検知することができない。

【0008】

そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、設備を構成するユニットの状態が不安定に変化する場合であっても、ユニットの異常の発生を容易に検知することの可能な情報を提供できる、新規かつ改良された設備管理支援装置、設備管理支援方法、プログラム及びコンピュータ読み取り可能な記録媒体を提供することにある。

【課題を解決するための手段】

【0009】

上記課題を解決するために、本発明のある観点によれば、製造業で用いられる設備を構成する複数のユニットのうち１つのユニットを評価対象として、複数の異なる期間に取得された当該評価対象のユニットの観測値データ群を、少なくとも当該設備における時空間データ及び時空間データに関係して変化する観測値を変数とした多次元のヒストグラムにそれぞれ変換するヒストグラム変換処理部と、複数のヒストグラムの差異を表す評価指標を算出する評価指標算出部と、評価指標に基づいて、評価対象のユニットの状態の変化を表す状態情報を生成する情報生成部と、を備える、設備管理支援装置が提供される。

【0010】

評価指標算出部は、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅを用いて評価指標を算出してもよい。

【0011】

また、評価指標算出部は、χ２距離を用いて評価指標を算出してもよい。

【0012】

あるいは、評価指標算出部は、ある期間の時空間データに対する観測値の分布を基づき算出される重みを考慮した重みつきχ２距離を用いて評価指標を算出してもよい。

【0013】

また、評価指標算出部は、ある期間の時空間データに対する観測値の分布を基づき算出されるＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅと、χ２距離または重みを考慮した重みつきχ２距離のうちいずれか一方とを組み合わせて、評価指標を算出してもよい。

【0014】

設備管理支援装置は、評価指標の変化特徴を表す符号を評価指標に付与した値を、ユニットの異常度合いを表す異常度として算出する異常度算出部をさらに備えてもよい。

【0015】

また、設備管理支援装置は、複数のユニットそれぞれに対して算出された異常度の高い順に、点検の必要なユニットの優先度を決定する優先度決定部をさらに備え、情報生成部は、状態情報の代わりに、あるいは状態情報とともに、優先度決定部により決定された優先度の高いユニットを作業者に通知するための通知情報を生成してもよい。

【0016】

情報生成部は、ユニットのトラブル発生を検知するために予め設定された閾値に基づき、評価指標を時系列に並べた情報あるいは異常度を時系列に並べた情報に対して閾値判定を行い、トラブル発生を予測した情報を生成してもよい。

【0017】

また、上記課題を解決するために、本発明の別の観点によれば、製造業で用いられる設備を構成する複数のユニットのうち１つのユニットを評価対象として、複数の異なる期間に取得された当該評価対象のユニットの観測値データ群を、少なくとも当該設備における時空間データ及び時空間データに関係して変化する観測値を変数とした多次元のヒストグラムにそれぞれ変換するヒストグラム変換処理ステップと、複数のヒストグラムの差異を表す評価指標を算出する評価指標算出ステップと、評価指標に基づいて、評価対象のユニットの状態の変化を表す状態情報を生成する情報生成ステップと、を含む、設備管理支援方法が提供される。

【0018】

さらに、上記課題を解決するために、本発明の別の観点によれば、コンピュータを、製造業で用いられる設備を構成する複数のユニットのうち１つのユニットを評価対象として、複数の異なる期間に取得された当該評価対象のユニットの観測値データ群を、少なくとも当該設備における時空間データ及び時空間データに関係して変化する観測値を変数とした多次元のヒストグラムにそれぞれ変換するヒストグラム変換処理部と、複数のヒストグラムの差異を表す評価指標を算出する評価指標算出部と、評価指標に基づいて、評価対象のユニットの状態の変化を表す状態情報を生成する情報生成部と、を備える、設備管理支援装置として機能させるプログラムが提供される。

【0019】

また、上記課題を解決するために、本発明の別の観点によれば、コンピュータに、製造業で用いられる設備を構成する複数のユニットのうち１つのユニットを評価対象として、複数の異なる期間に取得された当該評価対象のユニットの観測値データ群を、少なくとも当該設備における時空間データ及び時空間データに関係して変化する観測値を変数とした多次元のヒストグラムにそれぞれ変換するヒストグラム変換処理部と、複数のヒストグラムの差異を表す評価指標を算出する評価指標算出部と、評価指標に基づいて、評価対象のユニットの状態の変化を表す状態情報を生成する情報生成部と、を備える、設備管理支援装置として機能させるプログラムを記録したコンピュータ読み取り可能な記録媒体が提供される。

【発明の効果】

【0020】

以上説明したように本発明によれば、設備を構成するユニットの状態が不安定に変化する場合であっても、ユニットの異常の発生を容易に検知することの可能な情報を提供できる。

【図面の簡単な説明】

【0021】

【図1】本発明の第１の実施形態に係る設備管理支援装置の機能構成を示す機能ブロック図である。

【図2】同実施形態に係る設備管理支援方法を示すフローチャートである。

【図3】１ヶ月前から２ヶ月前までの押出負荷の傾向を示すデータとそのヒストグラムとを示す説明図である。

【図4】直近から１ヶ月前までの押出負荷の傾向を示すデータとそのヒストグラムとを示す説明図である。

【図5】Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）を説明する説明図である。

【図6】異常度算出部により算出される異常度を説明する説明図である。

【図7】評価指標への符号の与え方において、押出負荷の最大値に着目する場合と押出負荷のピーク位置に着目する場合とを説明するための図である。

【図8】χ２距離を説明するための説明図である。

【図9】異常度算出部により算出される異常度を説明する説明図である。

【図10】ヒストグラムの一例を示す説明図である。

【図11】同実施形態に係る設備管理支援装置のハードウェア構成例を示すブロック図である。

【図12】実施例として、ＥＭＤを時系列に並べた情報を示すグラフである。

【図13】実施例として、異常度を時系列に並べた情報を示すグラフである。

【図14】実施例として、第１の実施形態に係る設備管理支援方法により算出されるＥＭＤと、上記式（４）に基づき算出されたχ２距離とについての指標値の時間変化を示すグラフである。

【図15】実施例として、図１４に示したχ２距離と重みつきχ２距離とについての指標値の時間変化を示すグラフである。

【発明を実施するための形態】

【0022】

以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。

【0023】

また、以下では、一例として、コークス炉の炉体保全のため、コークス炉を構成する複数の窯の状態を管理するために本発明の設備管理支援装置を適用する場合について説明する。なお、本発明の設備管理支援装置は、例えば、以下の操業、品質、設備状態等を管理する設備に対しても適用可能である。

【0024】

（操業管理）
・コークス炉
：コークスケーキの表面温度、燃焼室フリュー温度や炭化室炉壁温度チャートを用いた燃焼状態悪化の検出
・連続鋳造設備
：連続鋳造のモールドにおける熱電対温度チャートを用いた凝固不良状態悪化の検出
（品質管理）
・熱間圧延設備、冷間圧延設備
：コイル片の疵発生数データ情報を用いた品質悪化傾向の検出
・棒線設備
：連続鋳造後の鋳片における位置毎の表面疵発生数データ情報を用いた品質悪化傾向の検出
・めっき設備
：連続溶融亜鉛めっきライン（ＣＧＬ）におけるコイル位置毎の合金化度のデータチャートを用いた品質悪化傾向の検出
（設備状態管理）
・原料：鉄鋼石石炭ヤードにおけるコンベヤのモータ電流を用いた設備異常検知

【0025】

＜１．第１の実施形態＞
［１－１．設備管理支援装置］
まず、図１を参照して、本発明の第１の実施形態に係る設備管理支援装置１００の機能構成について説明する。図１は、本実施形態に係る設備管理支援装置１００の機能構成を示す機能ブロック図である。

【0026】

本実施形態に係る設備管理支援装置１００は、図１に示すように、ヒストグラム変換処理部１１０と、評価指標算出部１２０と、異常度算出部１３０と、優先度決定部１４０と、情報生成部１５０と、出力部１６０とを備える。

【0027】

ヒストグラム変換処理部１１０は、設備を構成する複数のユニットのうち１つのユニットを評価対象として、複数の異なる期間に取得された当該評価対象のユニットの観測値データ群を、少なくとも当該設備における時空間データ及び時空間データに関係して変化する観測値を変数とした多次元のヒストグラムにそれぞれ変換する。時空間データは、位置及び時間のうち少なくともいずれか一方を示すデータであり、１次元情報（直線位置情報）、２次元情報（平面位置情報）であってもよい。また、時空間データに関係して変化する観測値は、１つの観測値に限らず、複数の観測値であってもよい。

【0028】

例えば、ヒストグラム変換処理部１１０は、コークス炉を構成する各窯について、コークスの押出時に取得された押出位置毎の押出負荷を含む押出負荷データを、少なくとも押出位置と押出負荷とを含む複数軸の座標系で表される多次元のヒストグラムにそれぞれ変換する。押出位置は時空間データであり、押出負荷は時空間データに関係して変化する観測値の１つである。押出負荷データは、操業実績データとして操業実績データ記憶部２００に蓄積されている。ヒストグラム変換処理部１１０は、所定期間に取得された複数の押出負荷データを用いて、押出位置毎の押出負荷の波形データをヒストグラムに変換して特徴量化する。このようなヒストグラムは、複数の異なる期間毎にそれぞれ作成される。ヒストグラム変換処理部１１０は、作成したヒストグラムを評価指標算出部１２０へ出力する。

【0029】

評価指標算出部１２０は、ヒストグラム変換処理部１１０により作成された複数のヒストグラムから、２つの異なる期間のヒストグラムの差異を表す評価指標を算出する。この評価指標により、評価対象とする窯における各押出位置での押出負荷の分布の変化を捉えることが可能となる。２つの異なる期間のヒストグラムの差異を評価する評価指標としては、本実施形態では、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）を用いる。評価指標算出部１２０は、算出した２つの異なる期間のヒストグラムの差異を表す評価指標を、異常度算出部１３０または情報生成部１５０のうち少なくともいずれか一方に出力する。

【0030】

異常度算出部１３０は、設備を構成するユニットの異常度合いを表す異常度を算出する。例えば、情報生成部１５０は、コークス炉の各窯について、窯の異常度合いを表す異常度を算出する。本実施形態では、異常度を、評価指標算出部１２０により算出された評価指標に対して、ユニット（例えば、窯）の状態が悪化したのか良化したのかを表す符号を付与することで表す。これは、評価指標の波形特徴変化だけでは窯状態が悪化したのか良化したのか判断できない場合があることによる。そこで、押出負荷の最大値（すなわち、コークスケーキ移動開始の評価）あるいは押出負荷が最大値を取るときの押出位置（すなわち、ピーク位置の評価）の変化を考慮した指標を、評価指標算出部１２０により算出された評価指標に対して付与する符号として用いて、異常度を算出する。異常度算出部１３０は、算出した異常度を、優先度決定部１４０または情報生成部１５０のうち少なくともいずれか一方に出力する。

【0031】

優先度決定部１４０は、複数のユニットそれぞれに対して算出された異常度に基づいて、点検の必要なユニットの優先度を決定する。例えば、優先度決定部１４０は、コークス炉の複数の窯それぞれに対して算出された異常度に基づいて、点検の必要な窯の優先度を決定する。優先度は、例えば１つの炉団を構成する複数の窯内で設定してもよい。優先度決定部１４０は、異常度が高いほどトラブル発生の可能性が高いとして、異常度の高い順に窯の点検の優先度を決定する。優先度決定部１４０は、決定した優先度を情報生成部１５０に出力する。

【0032】

情報生成部１５０は、評価対象のユニットの状態の変化を表す状態情報を生成する。状態情報は、例えば、設備の状態、設備における操業の状態、または、設備を用いて製造される製品の品質の状態のうち少なくともいずれか１つの変化を表す。

【0033】

例えば、情報生成部１５０は、状態情報として、評価指標に基づいて窯の状態を表す窯状態情報を生成する。窯状態情報としては、例えば、評価指標を時系列に並べた情報、異常度を時系列に並べた情報等がある。これらの情報は、時間の経過とともに窯の状態がどのように変化しているかを表している。また、情報生成部１５０は、破孔あるいは押詰りといったトラブルが発生する可能性が高くなったことを検知するために予め設定された閾値に基づき、評価指標を時系列に並べた情報あるいは異常度を時系列に並べた情報に対して閾値判定を行い、トラブル発生を予測した情報を生成してもよい。

【0034】

さらに情報生成部１５０は、優先度決定部１４０により決定されたユニット（例えば窯）の点検の優先度に基づいて、優先度の高いユニットを作業者に通知するための通知情報を生成してもよい。情報生成部１５０は、例えば、優先度が高い順から所定数の窯（例えば上位５つの窯）を通知する通知情報を生成してもよい。通知情報は、少なくとも窯を特定可能な情報（例えば窯番号）を含み、さらに通知する窯の窯状態情報も含んでもよい。情報生成部１５０により生成された情報は、出力部１６０に出力される。

【0035】

出力部１６０は、情報生成部１５０により生成された情報を作業者に通知するための機能部である。出力部１６０は、例えば情報を視覚的に提示するディスプレイであり、音声を出力するスピーカを備えていてもよい。作業者は、出力部１６０を介して提示された情報に基づいて、ユニットの状態を把握したり、点検あるいはメンテナンスの必要なユニットを認識したりすることができる。

【0036】

［１－２．設備管理支援方法］
図２～図７に基づいて、本実施形態に係る設備管理支援装置１００による設備管理支援方法を説明する。図２は、本実施形態に係る設備管理支援方法を示すフローチャートである。図３及び図４は、異なる時期における押出負荷の傾向を示すデータとそのヒストグラムとを示す説明図である。図５は、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）を説明する説明図である。図６は、異常度算出部１３０により算出される異常度を説明する説明図である。図７は、評価指標への符号の与え方において、押出負荷の最大値に着目する場合と押出負荷のピーク位置に着目する場合とを説明するための図である。

【0037】

（データ取得：Ｓ１００、ヒストグラム変換処理：Ｓ１１０）
設備管理支援装置１００は、操業実績データ記憶部２００から観測値データ群として所定期間の操業実績データが入力されると（Ｓ１００）、ヒストグラム変換処理部１１０により、操業実績データに含まれる押出負荷データがヒストグラムに変換される（Ｓ１１０）。

【0038】

ヒストグラムに変換される押出負荷データは、設備を構成するユニットである窯毎に、時系列順に配列された状態で、所定の期間毎に区分される。押出負荷データを区分する期間の長さは同一とし、データ数もほぼ同一であるのが望ましい。なお、各押出負荷データ群の期間は重複しないように設定されるが、多少の重複はあってもよい。例えば、直近２か月分の操業実績データの押出負荷データを取り扱う場合、直近から１か月前までの期間の押出負荷データと、１か月前から２か月前までの期間の押出負荷データとの２つの押出負荷データ群を作成する。ヒストグラム変換処理部１１０は、各期間の押出負荷データ群から、押出位置と押出負荷とを変数とした２次元のヒストグラムを生成する。

【0039】

押出負荷データ群と、その押出負荷データ群から変換して生成されたヒストグラムとを、図３及び図４に例示する。図３は１か月前から２か月前までの期間の押出負荷データに関する情報であり、図４は直近から１か月前までの期間の押出負荷データに関する情報である。図３及び図４の上側には、横軸に窯長手方向位置、縦軸に時間をとり、各窯長手方向位置での押出負荷の値を色合いで表現した図を示している。色が淡い部分ほど押出負荷が高いことを表している。また、図３及び図４の下側には、押出負荷データ群を変換して生成されたヒストグラムを示している。ヒストグラムは、押出位置と押出負荷とに応じた発生頻度（データ分布）を表している。

【0040】

図３及び図４より、所定の期間の押出負荷データ群をヒストグラムにより表すことで、各期間における押出位置と押出負荷とがとる値の傾向を捉えやすくなる。例えば、図３下側に示すように、１か月前から２か月前までの期間では、各押出位置における押出負荷の値におおきなばらつきは見られない。しかし、図４下側に示すように、直近から１か月前までの期間では、押出開始後、押出負荷が高まる傾向がみられる。

【0041】

（評価指標算出処理：Ｓ１２０）
ヒストグラム変換処理部１１０により各期間のヒストグラムが生成されると、評価指標算出部１２０は、２つのヒストグラムの差異を評価する（Ｓ１２０）。窯異常が発生し始めている場合には、押出負荷波形の特定位置において押出負荷にばらつきが現れ出すことがある。２つの異なる期間のヒストグラムの差異を見ることで、このような押出負荷のばらつき現象を捉えることができる。

【0042】

本実施形態では、２つの異なる期間のヒストグラムの差異を評価する評価指標としてＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）を用いる（非特許文献１参照）。Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅは、「特徴量」と「重み」の集合で与えられる複数の分布間の距離を表し、この距離が遠いほど分布間の差異が大きいことを意味する。本実施形態では、ステップＳ１１０で生成されたヒストグラムを評価対象の分布とし、各ヒストグラムの変数である押出位置及び押出負荷を「特徴量」として用い、それぞれの発生頻度（度数）を示すビンを「重み」として用いる。Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅは、高次元情報の比較に適した指標であり、外れ値に強く適切にヒストグラムの差異を評価することができる。

【0043】

より少し詳細に説明すると、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅは、最適化問題の１つである「輸送問題」の考え方に基づいている。例えば、図５に示すヒストグラムＰとヒストグラムＱとの差異を評価する場合には、ヒストグラムＰの場所ｉには、重み（すなわちビンの度数）の量だけ荷物が積まれているものとし、ヒストグラムＱとなるように、荷物を積み替えていくことを考える。つまり、図５で、ヒストグラムＱの各場所ｊにおいて、決められた重み（すなわち、比較するときのヒストグラムＱのビンの度数；点線のヒストグラム）の量の荷物をヒストグラムＰから倉庫へ格納していくと考える。このとき、ヒストグラムＰにある全ての荷物をヒストグラムＱに輸送するのに、どの荷物をどこへ輸送すれば最も効率良いのかを求め、その最適解（後述の仕事量を最小にする解）がヒストグラムＰとヒストグラムＱの差異を評価するための指標となる。

【0044】

そこで、まず、その輸送にかかる仕事量Ｗｏｒｋ（Ｐ，Ｑ）を下記式（１）のように定義し、この仕事量Ｗｏｒｋ（Ｐ，Ｑ）を最小化するときの荷物の輸送量（場所ｉから場所ｊへ輸送されるそれぞれの荷物の量）を求める。なお、ヒストグラムＰの場所ｉの荷物（ビン）について、ヒストグラムＱの場所ｊに移動させるときの輸送コストをｄ_ｉｊで表し、その輸送量（ビンの度数）をｆ_ｉｊで表す。

【0045】

【数1】

【0046】

なお、輸送コストｄ_ｉｊについては、多くの場合ではユークリッド距離が使用されており、本実施形態でもユークリッド距離を使用する。ここで、輸送コストｄ_ｉｊは事前に与えられているものとする。そして、下記式（２）を目的関数とする線形計画問題を解き、仕事量Ｗｏｒｋ（Ｐ，Ｑ）が最小のときの輸送量ｆ_＊ｉｊを求める。

【0047】

【数2】

【0048】

そして、求めたｆ_＊ｉｊを用いて、下記式（３）よりＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅの値ＥＭＤ（Ｐ，Ｑ）を求めることができる。

【0049】

【数3】

【0050】

ここで求められるＥＭＤ（Ｐ，Ｑ）は、上述したように、ヒストグラムＰの荷物（ビン）をヒストグラムＱに輸送するときの輸送問題の最適値と考えられる。このＥＭＤ（Ｐ，Ｑ）を２つのヒストグラムＰ、Ｑの類似度とみなし、ＥＭＤ（Ｐ，Ｑ）の値から２つのヒストグラムＰ、Ｑの差異を評価することができる。ＥＭＤ（Ｐ，Ｑ）の値が小さいほど、ヒストグラムＰとヒストグラムＱとは類似した分布を有しており、分布の変化が少ないといえる。

【0051】

評価指標算出部１２０は、上記式（３）より、２つのヒストグラムの差異を表す評価指標（すなわちＥＭＤ（Ｐ，Ｑ））を算出し、異常度算出部１３０または情報生成部１５０のうち少なくともいずれか一方に出力する。

【0052】

（異常度算出処理：Ｓ１３０）
次いで、異常度算出部１３０は、ユニット（例えば、窯）の異常度合いを表す異常度を算出する（Ｓ１３０）。本実施形態では、異常度を、評価指標算出部１２０により算出された評価指標に対して、ユニットの状態が悪化したのか良化したのかを表す符号を付与することで表す。これは、評価指標の波形特徴変化だけではユニットの状態が悪化したのか良化したのか判断できない場合があることによる。そこで、異常度算出部１３０は、観測値の１つである押出負荷について、最大値（すなわち、コークスケーキ移動開始の評価）あるいは押出負荷が最大値を取るときの押出位置（すなわち、ピーク位置の評価）の変化を考慮した指標を、ステップＳ１２０にて算出された評価指標に対して付与する符号として用いて、異常度を算出する。つまり、異常度は、図６に示すように、ステップＳ１２０にて算出されたＥＭＤ（Ｐ，Ｑ）に、ユニットの状態の変化の傾向を表す符号を付与したものといえる。以下では、符号が正（＋）の場合にはユニットの状態が悪化していることを表すものとし、符号が負（－）の場合にはユニットの状態が良化していることを表すものとする。

【0053】

評価指標への符号の与え方としては、例えば押出負荷の最大値に注目する場合には、図７に示すように、１回の押出し毎に取得された押出負荷の最大値（図７のピーク値Ａ）についての中央値を２つの異なる期間で比較し、その結果に基づき評価指標の符号を決定してもよい。例えば、１か月前から２か月前までの期間（図７の期間Ｔ２）の押出負荷データと、直近から１か月前までの期間（図７の期間Ｔ１）の押出負荷データとを考える。この場合、まず、異常度算出部１３０は、直近から１か月前までの期間における押出負荷データ毎に最大値を抽出して、それらの値の中央値を算出するとともに、１か月前から２か月前までの期間における押出負荷データ毎に最大値を抽出して、それらの値の中央値を算出する。そして、異常度算出部１３０は、これらの中央値の大小を判定する。

【0054】

１か月前から２か月前までの期間での中央値に対して直近から１か月前までの期間での中央値が大きい場合には、評価指標に正の符号が付与される。一方、１か月前から２か月前までの期間での中央値に対して直近から１か月前までの期間での中央値が小さい場合には、評価指標に負の符号が付与される。なお、本実施形態では、日ごとに押出負荷に関するデータを取得しており、その日ごとに押出負荷の最大値を求めて蓄積しているものとし、各期間の押出負荷の最大値の中央値はこのように蓄積されたデータを用いて算出できる。

【0055】

また、評価指標への符号の与え方として、押出負荷のピーク位置に注目する場合には、図７に示すように、１回の押出し毎に取得された各ヒストグラムでの押出負荷のピーク位置（図７のピーク位値Ｂ）についての中央値を２つの異なる期間で比較し、その結果に基づき評価指標の符号を決定してもよい。例えば、１か月前から２か月前までの期間（図７の期間Ｔ２）の押出負荷データと、直近から１か月前までの期間（図７の期間Ｔ１）の押出負荷データとを考える。この場合、まず、異常度算出部１３０は、直近から１か月前までの期間における押出負荷データ毎にピーク位置を抽出して、それらの値の中央値を算出するとともに、１か月前から２か月前における押出負荷データ毎にピーク位置を抽出して、それらの値の中央値を算出する。そして、異常度算出部１３０は、これらの中央値の変化を判定する。

【0056】

直近から１か月前までの期間でのピーク位置の中央値が窯の排出口側へ移動している場合（すなわち、１か月前から２か月前までの期間でのピーク位置の中央値よりも大きい場合）には、評価指標に正の符号が付与される。一方、直近から１か月前までの期間でのピーク位置の中央値が押出機側へ移動している場合（すなわち、１か月前から２か月前までの期間でのピーク位置の中央値よりも小さい場合）には、評価指標に負の符号が付与される。なお、本実施形態では、日ごとに押出負荷に関するデータを取得しており、その日ごとに押出負荷のピーク位置を特定して蓄積しているものとし、各期間の押出負荷のピーク位置の中央値はこのように蓄積されたデータを用いて算出できる。

【0057】

このように、評価指標に対してユニットの状態の変化の傾向を示す符号を付与することで、ユニットの異常の発生をより適切に検知できるようになる。

【0058】

（確認処理：Ｓ１４０）
ステップＳ１３０の処理を終えると、設備管理支援装置１００は、評価対象となるすべてのユニットについて、ステップＳ１１０～ステップＳ１３０の処理が行われたか否か判定する（Ｓ１４０）。例えば、コークス炉の例では、評価対象とする窯は、例えば１つの炉団を構成している窯すべてとしてもよい。ステップＳ１４０にて未処理のユニットがあると判定された場合には、評価対象となるすべてのユニットについて処理を終えるまでステップＳ１１０～ステップＳ１３０の処理が繰り返し実行される（ステップＳ１４０；Ｎｏ）。一方、ステップＳ１４０にてすべてのユニットについて処理が行われたと判定された場合には（ステップＳ１４０；Ｙｅｓ）、ステップＳ１５０以降の処理が実行される。

【0059】

（優先度決定処理：Ｓ１５０）
評価対象とするすべてのユニットについてステップＳ１１０～Ｓ１３０の処理を終えると、優先度決定部１４０は、ステップＳ１３０にて複数のユニットそれぞれに対して算出された異常度に基づいて、点検の必要なユニットの優先度を決定する（Ｓ１５０）。例えば、コークス炉の例では、優先度は、例えば１つの炉団を構成する複数の窯に対して設定してもよい。優先度決定部１４０は、異常度が高いほどトラブル発生の可能性が高いとして、異常度の高い順にユニットの点検の優先度を決定する。

【0060】

（出力処理：Ｓ１６０）
ステップＳ１５０にてユニットの点検の優先度が決定されると、情報生成部１５０は、優先度の高いユニットを作業者に通知するための通知情報を生成する。出力部１６０は、情報生成部１５０により生成された通知情報を作業者に通知する。通知情報により、例えば、コークス炉の例では、窯の点検の優先度の高い窯の窯番号が作業者に通知される。作業者は、通知情報に基づき、優先的に炉内点検あるいはメンテナンスすべき窯を把握することができ、早期に必要なアクションをとることが可能となる。

【0061】

以上、本実施形態に係る設備管理支援方法について説明した。本実施形態に係る設備管理支援方法を用いることで、複数の異なる期間の観測値データを評価範囲として、それらの期間を逐次的に変更しながら解析を行うことが可能となる。すなわち、予め設定されたトレンドに対して評価を行うのではなく、評価対象とする所定期間の観測値データを変更して評価することで、設備を構成するユニットが不安定に変化する場合であっても、ユニットの異常を捉えることが可能となる。これにより、コークス炉の例では、例えば押詰りあるいは破孔といった窯異常を高精度に予見することができるようになる。その結果、窯異常に対して適切な対策を早期に実施することが可能となり、非稼働となる窯の発生を防止して操業安定化を図ることが可能となる。

【0062】

＜２．第２の実施形態＞
次に、図８～図１０に基づいて、本発明の第２の実施形態に係る設備管理支援装置について説明する。本実施形態に係る設備管理支援装置は、第１の実施形態に係る設備管理支援装置１００と比較して、２つの異なる期間のヒストグラムの差異を評価する評価指標として、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）の代わりにχ２距離を用いる点で相違する。

【0063】

なお、本実施形態に係る設備管理支援装置の構成は、図１に示した設備管理支援装置１００の機能構成と同一である。このため、本実施形態に係る設備管理支援装置についても、図１に示した第１の実施形態の設備管理支援装置１００と同一の符号を用いて説明し、各機能部の詳細な説明については省略する。また、設備管理支援方法の流れは、図２に示したフローチャートと同一である。以下、図２に沿って本実施形態に係る設備管理支援方法について説明する。

【0064】

［２－１．設備管理支援方法］
（データ取得：Ｓ１００、ヒストグラム変換処理：Ｓ１１０）
本実施形態に係る設備管理支援方法では、まず、設備管理支援装置１００に操業実績データ記憶部２００から観測値データ群として所定期間の操業実績データが入力されると（Ｓ１００）、ヒストグラム変換処理部１１０により、操業実績データに含まれる押出負荷データがヒストグラムに変換される（Ｓ１１０）。ステップＳ１００及びＳ１１０は、第１の実施形態と同様に行えばよいため、ここでは詳細な説明を省略する。ステップＳ１１０の処理により、図３及び図４に示したようなヒストグラムが取得される。

【0065】

（評価指標算出処理：Ｓ１２０）
ヒストグラム変換処理部１１０により各期間のヒストグラムが生成されると、評価指標算出部１２０は、２つのヒストグラムの差異を評価する（Ｓ１２０）。窯異常が発生し始めている場合には、押出負荷波形の特定位置において押出負荷にばらつきが現れ出すことがある。２つの異なる期間のヒストグラムの差異を見ることで、このような押出負荷のばらつき現象を捉えることができる。

【0066】

本実施形態では、２つの異なる期間のヒストグラムの差異を評価する評価指標として、χ２距離を用いる。χ２距離は独立性の検定等で用いられるχ２統計量に類似する指標であり、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅと同様にヒストグラムの類似比較を行う際に利用される。χ２距離は、下記式（４）により表される。また、図８に、χ２距離を説明するための説明図を示す。なお、ｉは、ヒストグラムＰ及びヒストグラムＱでの場所を表している。

【0067】

【数4】

【0068】

χ２距離では比較する２つのヒストグラムＰ、Ｑについて、互いのビンの度数が０となる場合にそのビンの差分の計算は行わず、図８に示すように、ヒストグラムＰのビンまたはヒストグラムＱのビンのうち少なくともいずれか一方の度数が０ではない場合に、上記式（４）に示すようにビンの差分の平方値を計算し、その平方値をＰ、Ｑのビン数の平均値で割った値の総和をχ２距離とする。評価指標算出部１２０は、上記式（４）より、２つのヒストグラムの差異を表す評価指標（すなわちχ２距離）を算出し、異常度算出部１３０または情報生成部１５０のうち少なくともいずれか一方に出力する。

【0069】

（異常度算出処理：Ｓ１３０）
次いで、異常度算出部１３０は、ユニット（例えば、窯）の異常度合いを表す異常度を算出する（Ｓ１３０）。本実施形態でも第１の実施形態と同様、異常度を、評価指標算出部１２０により算出された評価指標に対して、ユニットの状態が悪化したのか良化したのかを表す符号を付与することで表す。これは、評価指標の波形特徴変化だけではユニットの状態が悪化したのか良化したのか判断できない場合があることによる。そこで、異常度算出部１３０は、観測値の１つである押出負荷について、最大値（すなわち、コークスケーキ移動開始の評価）あるいは押出負荷が最大値を取るときの押出位置（すなわち、ピーク位置の評価）の変化を考慮した指標を、ステップＳ１２０にて算出された評価指標に対して付与する符号として用いて、異常度を算出する。つまり、異常度は、図９に示すように、ステップＳ１２０にて算出された評価指標（χ２距離）に、ユニットの状態の変化の傾向を表す符号を付与したものといえる。以下では、符号が正（＋）の場合にはユニットの状態が悪化していることを表すものとし、符号が負（－）の場合にはユニットの状態が良化していることを表すものとする。

【0070】

評価指標への符号の与え方としては、例えば押出負荷の最大値に注目する場合には、２つの異なる期間のヒストグラムにおいて、各ヒストグラムでの押出負荷の中央値を比較した結果に基づき、評価指標の符号を決定してもよい。あるいは、評価指標への符号の与え方として、押出負荷のピーク位置に注目する場合には、２つの異なる期間のヒストグラムにおいて、各ヒストグラムでの押出負荷のピーク位置の中央値を比較した結果に基づき、評価指標の符号を決定してもよい。かかる評価指標の与え方は、第１の実施形態と同様に行えばよい。このように、評価指標に対してユニットの状態の変化の傾向を示す符号を付与することで、ユニットの異常の発生をより適切に検知できるようになる。

【0071】

（確認処理：Ｓ１４０）
ステップＳ１３０の処理を終えると、設備管理支援装置１００は、評価対象となるすべてのユニットについて、ステップＳ１１０～ステップＳ１３０の処理が行われたか否か判定する（Ｓ１４０）。ステップＳ１４０の処理は、第１の実施形態と同様に行えばよい。例えば、コークス炉の例では、評価対象とする窯は、例えば１つの炉団を構成している窯すべてとしてもよい。ステップＳ１４０にて未処理のユニットがあると判定された場合には、評価対象となるすべてのユニットについて処理を終えるまでステップＳ１１０～ステップＳ１３０の処理が繰り返し実行される（ステップＳ１４０；Ｎｏ）。一方、ステップＳ１４０にてすべてのユニットについて処理が行われたと判定された場合には（ステップＳ１４０；Ｙｅｓ）、ステップＳ１５０以降の処理が実行される。

【0072】

（優先度決定処理：Ｓ１５０）
評価対象とするすべてのユニットについてステップＳ１１０～Ｓ１３０の処理を終えると、優先度決定部１４０は、ステップＳ１３０にて複数のユニットそれぞれに対して算出された異常度に基づいて、点検の必要なユニットの優先度を決定する（Ｓ１５０）。ステップＳ１５０の処理は、第１の実施形態と同様に行えばよい。例えば、コークス炉の例では、優先度は、例えば１つの炉団を構成する複数の窯内で設定してもよい。優先度決定部１４０は、異常度が高いほどトラブル発生の可能性が高いとして、異常度の高い順にユニットの点検の優先度を決定する。

【0073】

ステップＳ１５０にてユニットの点検の優先度が決定されると、情報生成部１５０は、優先度の高いユニットを作業者に通知するための通知情報を生成する。出力部１６０は、情報生成部１５０により生成された通知情報を作業者に通知する。通知情報により、例えば、コークス炉の例では、窯の点検の優先度の高い窯の窯番号が作業者に通知される。作業者は、通知情報に基づき、優先的に炉内点検あるいはメンテナンスすべき窯を把握することができ、早期に必要なアクションをとることが可能となる。

【0074】

［２－２．変形例］
例えば、コークス炉の例において、押出負荷波形には、図１０に示すように、押出負荷の値が取りうる分布に特徴がある。その傾向としては、押出位置が押出機側（ＰＳ；Ｐｕｓｈｅｒ Ｓｉｄｅ）では負荷分布が大きい値の範囲でばらつき、コークス排出側（ＣＳ；Ｃｏｋｅ Ｓｉｄｅ）では負荷分布が小さい値の範囲でばらつくように分布が変化しており、押出位置によって押出負荷のとりうる範囲が異なる。

【0075】

そこで、χ２距離の計算において、押出負荷の値として取りうる頻度が低いビンにおける変化の影響度を大きく評価し、取りうる頻度が高いビンの位置における変化の影響度を小さく評価することで、分布変化をより敏感に捉えられるようにする。これにより、例えば図１０において、コークスケーキの押出しの開始時の挙動の良し悪しを見る上で重要な、押出負荷がピークとなる付近（領域Ｃ）の押出負荷のばらつき傾向を感度よく把握することができる。また、図１０において、コークスケーキの押出し中盤での挙動の良し悪しを見る上で重要な、押出負荷がピークとなるピーク位置がコークス排出側へ移動したときの押出負荷の変化傾向を感度よく把握することもできる。具体的には、図１０の領域Ｄの押出負荷が大きくなる方へばらつくほど、押出し性が悪化していることがわかる。

【0076】

具体的には、χ２距離を算出する上記式（４）にてビン毎の重みｗ_ｉが考慮された、重みつきχ２距離を求める。重みつきχ２距離は、下記式（５）により表される。

【0077】

【数5】

【0078】

ここで、重みｗ_ｉは、下記式（６）に示すように、例えば、ある期間（例えば、補修後等の押出負荷が安定している時期、等）における押出負荷の２次元ヒストグラムを蓄積しておき、各ビンの位置についてビン数Ｘ_ｉの逆数を取ることにより設定してもよい。なお、ビン数Ｘ_ｉが０の場合には１が与えられる。

【0079】

【数6】

【0080】

あるいは、下記式（７）のように、ビン数Ｘ_ｉを全体のヒストグラムの最大ビン数Ｘ_Ｍａｘで割った値を１から引き、重みｗ_ｉとしてもよい。

【0081】

【数7】

【0082】

このように、重みつきχ２距離を用いることで、観測値として取りうる頻度が低いビンにおける変化の影響度を大きく評価し、取りうる頻度が高いビンの位置における変化の影響度を小さく評価することで、観測値の分布変化をより敏感に捉えることができる。

【0083】

以上、本実施形態に係る設備管理支援方法について説明した。本実施形態に係る設備管理支援方法を用いることで、複数の異なる期間の観測値データを評価範囲として、それらの期間を逐次的に変更しながら解析を行うことが可能となる。すなわち、予め設定されたトレンドに対して評価を行うのではなく、評価対象とする所定期間の観測値データを変更して評価することで、設備を構成するユニットが不安定に変化する場合であっても、ユニットの異常を捉えることが可能となる。これにより、コークス炉の例では、例えば押詰りあるいは破孔といった窯異常を高精度に予見することができるようになる。その結果、窯異常に対して適切な対策を早期に実施することが可能となり、非稼働となる窯の発生を防止して操業安定化を図ることが可能となる。

【0084】

また、本実施形態では、２つの異なる期間のヒストグラムの差異を評価する評価指標として、χ２距離を用いることで、感度よく観測値の変化をとらえることが可能となる。

【0085】

なお、２つの異なる期間のヒストグラムの差異を評価する評価指標として、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅと、χ２距離または重みを考慮した重みつきχ２距離のうちいずれか一方との、２つの値を組み合わせて得られる値を用いてもよい。例えば、本実施形態で用いたχ２距離または重みつきχ２距離と、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）との２つの平均値を用いてもよい。あるいは、評価指標として、χ２距離または重みつきχ２距離とＥＭＤとを切り替えながら使用することも考えられる。この場合、例えば、χ２距離を評価指標として用いることをベースとして、ヒストグラムの差異の評価に基づき、引き続き評価指標としてχ２距離を用いるか、あるいは、ＥＭＤを評価指標として用いるかを決定してもよい。具体的には、χ２距離が急峻に増加した場合には、評価指標をＥＭＤの値に切り替えて、ＥＭＤの値に同様の増加傾向があるかを確認して、その分布差の変動の妥当性を評価することにより、ユニットの異常を検知する。

【0086】

このように、分布差の変動を高感度にとらえるχ２距離または重みつきχ２距離と、ロバスト性の高いＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）とを組み合わせて使用することにより、ユニットの状況の変化をより的確にとらえることが可能となる。

【0087】

＜３．ハードウェア構成例＞
以下、図１１を参照しながら、本実施形態に係る設備管理支援装置１００のハードウェア構成について、詳細に説明する。図１１は、本発明の実施形態に係る設備管理支援装置１００のハードウェア構成例を示すブロック図である。

【0088】

設備管理支援装置１００は、主に、ＣＰＵ９０１と、ＲＯＭ９０３と、ＲＡＭ９０５と、を備える。また、設備管理支援装置１００は、更に、バス９０７と、入力装置９０９と、出力装置９１１と、ストレージ装置９１３と、ドライブ９１５と、接続ポート９１７と、通信装置９１９とを備える。

【0089】

ＣＰＵ９０１は、演算処理装置および制御装置として機能し、ＲＯＭ９０３、ＲＡＭ９０５、ストレージ装置９１３、またはリムーバブル記録媒体９２１に記録された各種プログラムに従って、設備管理支援装置１００内の動作全般またはその一部を制御する。ＲＯＭ９０３は、ＣＰＵ９０１が使用するプログラムや演算パラメータ等を記憶する。ＲＡＭ９０５は、ＣＰＵ９０１が使用するプログラムや、プログラムの実行において適宜変化するパラメータ等を一次記憶する。これらはＣＰＵバス等の内部バスにより構成されるバス９０７により相互に接続されている。

【0090】

バス９０７は、ブリッジを介して、ＰＣＩ（Ｐｅｒｉｐｈｅｒａｌ Ｃｏｍｐｏｎｅｎｔ Ｉｎｔｅｒｃｏｎｎｅｃｔ／Ｉｎｔｅｒｆａｃｅ）バスなどの外部バスに接続されている。

【0091】

入力装置９０９は、例えば、マウス、キーボード、タッチパネル、ボタン、スイッチおよびレバーなどユーザが操作する操作手段である。また、入力装置９０９は、例えば、赤外線やその他の電波を利用したリモートコントロール手段（いわゆる、リモコン）であってもよいし、設備管理支援装置１００の操作に対応したＰＤＡ等の外部接続機器９２３であってもよい。さらに、入力装置９０９は、例えば、上記の操作手段を用いてユーザにより入力された情報に基づいて入力信号を生成し、ＣＰＵ９０１に出力する入力制御回路などから構成されている。設備管理支援装置１００のユーザは、この入力装置９０９を操作することにより、設備管理支援装置１００に対して各種のデータを入力したり処理動作を指示したりすることができる。

【0092】

出力装置９１１は、取得した情報をユーザに対して視覚的または聴覚的に通知することが可能な装置で構成される。このような装置として、ＣＲＴディスプレイ装置、液晶ディスプレイ装置、プラズマディスプレイ装置、ＥＬディスプレイ装置およびランプなどの表示装置や、スピーカおよびヘッドホンなどの音声出力装置や、プリンタ装置、携帯電話、ファクシミリなどがある。出力装置９１１は、例えば、設備管理支援装置１００が行った各種処理により得られた結果を出力する。具体的には、表示装置は、設備管理支援装置１００が行った各種処理により得られた結果を、テキストまたはイメージで表示する。他方、音声出力装置は、再生された音声データや音響データ等からなるオーディオ信号をアナログ信号に変換して出力する。

【0093】

ストレージ装置９１３は、設備管理支援装置１００の記憶部の一例として構成されたデータ格納用の装置である。ストレージ装置９１３は、例えば、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）等の磁気記憶部デバイス、半導体記憶デバイス、光記憶デバイス、または光磁気記憶デバイス等により構成される。このストレージ装置９１３は、ＣＰＵ９０１が実行するプログラムや各種データ、および外部から取得した各種のデータなどを格納する。

【0094】

ドライブ９１５は、記録媒体用リーダライタであり、設備管理支援装置１００に内蔵、あるいは外付けされる。ドライブ９１５は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体９２１に記録されている情報を読み出して、ＲＡＭ９０５に出力する。また、ドライブ９１５は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体９２１に記録を書き込むことも可能である。リムーバブル記録媒体９２１は、例えば、ＣＤメディア、ＤＶＤメディア、Ｂｌｕ－ｒａｙ（登録商標）メディア等である。また、リムーバブル記録媒体９２１は、コンパクトフラッシュ（登録商標）（ＣｏｍｐａｃｔＦｌａｓｈ：ＣＦ）、フラッシュメモリ、または、ＳＤメモリカード（Ｓｅｃｕｒｅ Ｄｉｇｉｔａｌ ｍｅｍｏｒｙ ｃａｒｄ）等であってもよい。また、リムーバブル記録媒体９２１は、例えば、非接触型ＩＣチップを搭載したＩＣカード（Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ ｃａｒｄ）または電子機器等であってもよい。

【0095】

接続ポート９１７は、機器を設備管理支援装置１００に直接接続するためのポートである。接続ポート９１７の一例として、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ Ｓｅｒｉａｌ Ｂｕｓ）ポート、ＩＥＥＥ１３９４ポート、ＳＣＳＩ（Ｓｍａｌｌ Ｃｏｍｐｕｔｅｒ Ｓｙｓｔｅｍ Ｉｎｔｅｒｆａｃｅ）ポート、ＲＳ－２３２Ｃポート等がある。この接続ポート９１７に外部接続機器９２３を接続することで、設備管理支援装置１００は、外部接続機器９２３から直接各種のデータを取得したり、外部接続機器９２３に各種のデータを提供したりする。

【0096】

通信装置９１９は、例えば、通信網９２５に接続するための通信デバイス等で構成された通信インタフェースである。通信装置９１９は、例えば、有線または無線ＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）、Ｂｌｕｅｔｏｏｔｈ（登録商標）、またはＷＵＳＢ（Ｗｉｒｅｌｅｓｓ ＵＳＢ）用の通信カード等である。また、通信装置９１９は、光通信用のルータ、ＡＤＳＬ（Ａｓｙｍｍｅｔｒｉｃ Ｄｉｇｉｔａｌ Ｓｕｂｓｃｒｉｂｅｒ Ｌｉｎｅ）用のルータ、または、各種通信用のモデム等であってもよい。この通信装置９１９は、例えば、インターネットや他の通信機器との間で、例えばＴＣＰ／ＩＰ等の所定のプロトコルに則して信号等を送受信することができる。また、通信装置９１９に接続される通信網９２５は、有線または無線によって接続されたネットワーク等により構成され、例えば、インターネット、家庭内ＬＡＮ、赤外線通信、ラジオ波通信または衛星通信等であってもよい。

【0097】

以上、本発明の実施形態に係る設備管理支援装置１００の機能を実現可能なハードウェア構成の一例を示した。上記の各構成要素は、汎用的な部材を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。従って、本実施形態を実施する時々の技術レベルに応じて、適宜、利用するハードウェア構成を変更することが可能である。

【実施例】

【0098】

本発明の第１の実施形態に係る設備管理支援方法により算出されるＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ及び異常値について、コークス炉の窯状態との関係の有無について確認した。図１２は、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅを時系列に並べた情報であり、図１３は、異常度を時系列に並べた情報である。１日１回、直近２か月の操業実績データ（押出負荷データ）を取得し、取得する毎にその日から１か月前までの押出負荷データから作成されたヒストグラムと、１か月前から２か月前までの押出負荷データから作成されたヒストグラムとを比較した。そして、上記式（２）からＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅを算出し、さらにＥａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅに対して符号を付与し、異常値を算出した。符号は、押出負荷の最大値に基づき付与した。

【0099】

図１２の例では、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅが増加した後、破孔が発生した。すなわち、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅの増加傾向をみることで、破孔予測が可能であるといえる。また、図１３の例では、異常値が増加した後、コークスの押詰りが発生した。異常値の場合も同様に、異常値の増加傾向をみることで、押詰り予測が可能であるといえる。これより、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅあるいは異常値の変化を監視し、閾値判定を行うことで、破孔あるいは押詰りといったトラブル発生による窯の非稼働化を防ぐことが可能となる。

【0100】

また、本発明の第２の実施形態に係る設備管理支援方法により算出されるχ２距離または重みつきχ２距離及び異常値について、コークス炉の窯状態との関係の有無について確認した。図１４に、第１の実施形態に係る設備管理支援方法により算出されるＥＭＤと、上記式（４）に基づき算出されたχ２距離とについての指標値の時間変化を示す。図１５に、図１４に示したχ２距離と、上記式（５）に基づき算出された重みつきχ２距離とについての指標値の時間変化を示す。図１４及び図１５の指標値は、各評価指標（すなわち、ＥＭＤ、χ２距離、重みつきχ２距離）の値をそれぞれが取り得る最大値で除した値である。かかる指標値は、日ごとに０．１ほど変化すれば、その変化は急峻であるといえる状況になっていることを表している。なお、重みについては、炉体補修直後の正常な状態と推測される期間データを用いて算出した。

【0101】

図１４をみると、ＥＭＤの指標値はなだらかに増加しているのに対して、χ２距離の指標値はＥＭＤよりも早く指標値に上昇傾向があることがわかる。このように、評価指標としてχ２距離を用いることで、評価指標としてＥＭＤを用いた場合と同様、押詰り前に指標値が上昇し、より早い段階で閾値判定によってトラブル発生を検知することができる。

【0102】

また、図１５をみると、重みつきχ２距離を用いた場合、χ２距離を用いる場合よりもさらに早い段階で指標値が増加していることがわかる。また、重みつきχ２距離の指標値がピークとなった後、少し低下した後に押詰りが発生した。重みつきχ２距離はχ２距離よりも観測値の変化を感度よくとらえることができるため、このような指標値の傾向が現れたものと考えられる。この傾向に基づきトラブルが発生する可能性が高くなったことを検知することができる。

【0103】

図１４及び図１５の結果からも、例えば、χ２距離を評価指標として用いることをベースとして、χ２距離が急峻に増加した場合には、評価指標をロバスト性の高いＥＭＤの値に切り替えて、ＥＭＤの値に同様の増加傾向があるかを確認することで、状況の変化をより的確にとらえられることがわかる。

【0104】

以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。

【0105】

例えば、上記実施形態では、２つの異なる期間のヒストグラムの差異を評価する評価指標として、Ｅａｒｔｈ Ｍｏｖｅｒ’ｓ Ｄｉｓｔａｎｃｅ（ＥＭＤ）あるいはχ２距離を用いたが、本発明はかかる例に限定されない。例えば、ＫＬ情報量、ＪＳ情報量等を評価指標として用いてもよい。

【0106】

また、上記実施形態では、出力処理（Ｓ１６０）においてユニットの異常度に基づき決定された点検の優先度の高いユニットを通知する例を取り上げたが、本発明はかかる例に限定されず、例えば、上記の通知情報に代えて、あるいは通知情報とともに、評価対象のユニットの状態の変化を表す状態情報を作業者に通知してもよい。コークス炉の例では、状態情報として、例えば、評価指標を時系列に並べた情報、異常度を時系列に並べた情報等の窯状態情報が通知される。これらの情報は、時間の経過とともに窯の状態がどのように変化しているかを表しているため、設備管理支援において有益な情報といえる。

【0107】

さらに、上記実施形態では、本実施形態に係る情報生成部１５０は、破孔あるいは押詰りといったトラブルが発生する可能性が高くなったことを検知するために予め設定された閾値に基づき、評価指標を時系列に並べた情報あるいは異常度を時系列に並べた情報に対して閾値判定を行い、トラブル発生を予測した情報を生成することも可能である。したがって、出力処理では、閾値判定の結果からトラブル発生が予測されるユニットを通知するようにしてもよい。閾値判定を行うことで、コークス炉の例であれば破孔あるいは押詰りといったトラブル発生によるユニットの非稼働化を防ぐことが可能となる。

【0108】

また、上記実施形態では、２つのヒストグラムの差異を、評価対象のユニットの状態の変化を表す評価指標としたが、本発明はかかる例に限定されない。例えば、３つ以上のヒストグラムを比較することにより算出されたこれらのヒストグラムの差異を評価指標としてもよい。

【符号の説明】

【0109】

１００ 設備管理支援装置
１１０ ヒストグラム変換処理部
１２０ 評価指標算出部
１３０ 異常度算出部
１４０ 優先度決定部
１５０ 情報生成部
１６０ 出力部
２００ 操業実績データ記憶部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】