特許 5794034 障害予測システム及びプログラム

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

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

(11)【特許番号】5794034

(24)【登録日】2015年8月21日

(45)【発行日】2015年10月14日

(54)【発明の名称】障害予測システム及びプログラム

(51)【国際特許分類】

   G03G 21/00 20060101AFI20150928BHJP

   G05B 23/02 20060101ALI20150928BHJP

【ＦＩ】

   G03G21/00 386

   G05B23/02 R

   G05B23/02 T

   G05B23/02 302Y

【請求項の数】6

【全頁数】14

(21)【出願番号】特願2011-178901(P2011-178901)

(22)【出願日】2011年8月18日

(65)【公開番号】特開2013-41173(P2013-41173A)

(43)【公開日】2013年2月28日

【審査請求日】2014年7月11日

(73)【特許権者】

【識別番号】000005496

【氏名又は名称】富士ゼロックス株式会社

(74)【代理人】

【識別番号】100116687

【弁理士】

【氏名又は名称】田村 爾

(74)【代理人】

【識別番号】100098132

【弁理士】

【氏名又は名称】守山 辰雄

(72)【発明者】

【氏名】安川 薫

【審査官】 齋藤 卓司

(56)【参考文献】

【文献】 国際公開第２０１０／１１００２１（ＷＯ，Ａ１）

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

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

【文献】 特開平０８−３２０７２６（ＪＰ，Ａ）

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

Ｇ０３Ｇ ２１／００

Ｇ０５Ｂ ２３／０２

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

【請求項1】

被監視装置の内部状態を示す複数のパラメータについて、値の変化に因果関係がある一対のパラメータから成る組を複数設定する設定手段と、
前記設定手段により設定された組毎に得られる、各組を成すパラメータ間における値の時系列変化の相関度の推移が最も類似する２つの組について、共通するパラメータが存在する場合には当該共通するパラメータを特定し、共通するパラメータが存在しない場合には、当該２つの組について、一方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度、及び、他方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度を調べ、その中で相関度の推移が最も類似する２つの組に共通するパラメータを特定する特定手段と、
前記特定手段により特定されたパラメータに予め対応付けられた障害が被監視装置に発生する旨の通知を出力する出力手段と、
を備えたことを特徴とする障害予測システム。

【請求項2】

前記設定手段は、前記被監視装置の内部状態を示す複数のパラメータを少なくとも帯電電圧、トナー濃度、現像及びレーザパワーの各類に分類したパラメータ群毎に、そのパラメータ群に分類された各パラメータについて、正常状態の前記被監視装置において、前記値の変化に因果関係がある一対のパラメータからなる組を複数設定する、
ことを特徴とする請求項１に記載の障害予測システム。

【請求項3】

前記特定手段は、前記各組を成すパラメータ間の相関係数を対象範囲の期間をスライドさせた複数の期間についてそれぞれ算出し、当該算出した各期間の相関係数を時系列順に並べて相関データ列とし、当該相関データ列間の距離が近いほど相関度の推移が類似すると判断する、
ことを特徴とする請求項１又は請求項２に記載の障害予測システム。

【請求項4】

コンピュータに、
被監視装置の内部状態を示す複数のパラメータについて、値の変化に因果関係がある一対のパラメータから成る組を複数設定した設定機能と、
前記設定機能により設定された組毎に得られる、各組を成すパラメータ間における値の時系列変化の相関度の推移が最も類似する２つの組について、共通するパラメータが存在する場合には当該共通するパラメータを特定し、共通するパラメータが存在しない場合には、当該２つの組について、一方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度、及び、他方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度を調べ、その中で相関度の推移が最も類似する２つの組に共通するパラメータを特定する特定機能と、
前記特定機能により特定されたパラメータに予め対応付けられた障害が被監視装置に発生する旨の通知を出力する出力機能と、
を実現するためのプログラム。

【請求項5】

前記設定機能は、前記被監視装置の内部状態を示す複数のパラメータを少なくとも帯電電圧、トナー濃度、現像及びレーザパワーの各類に分類したパラメータ群毎に、そのパラメータ群に分類された各パラメータについて、正常状態の前記被監視装置において、前記値の変化に因果関係がある一対のパラメータからなる組を複数設定する、
ことを特徴とする請求項４に記載のプログラム。

【請求項6】

前記特定機能は、前記各組を成すパラメータ間の相関係数を対象範囲の期間をスライドさせた複数の期間についてそれぞれ算出し、当該算出した各期間の相関係数を時系列順に並べて相関データ列とし、当該相関データ列間の距離が近いほど相関度の推移が類似すると判断する、
ことを特徴とする請求項４又は請求項５に記載のプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、障害予測システム及びプログラムに関する。

【背景技術】

【0002】

従来、プラントなどのように連続的に運転されるシステムについて、観測により得られたデータが閾値を超えるか否かで異常を特定する方法や、２つのパラメータ間の相互相関係数を利用して異常を特定する方法などが考案されている。

【0003】

例えば、木構造的に因果関係の関連する多数のプロセス値を有するプロセスにおいて、結果側から原因側に１つの経路に沿って順次プロセス値を辿って、隣接する２つのプロセス値毎に夫々の経時的測定データの相互相関係数を求めると共に、上記経路に沿った複数の相互相関係数の相乗平均を求め、全ての可能な経路につき求めた上記相乗平均の大きさにより、結果側のプロセス値の変動に対する原因側のプロセス値の影響の程度を分析することを特徴とするプロセスにおける変動の分析方法の発明が提案されている（特許文献１参照）。

【0004】

例えば、各々が複数のセンサを備える複数の設備を遠隔監視する遠隔監視システムであって、前記設備から送信される、前記複数のセンサで検知された実測センサ値を、随時取得するセンサ値取得手段と、前記設備が正常運転をしている際に、前記複数のセンサで検知され前記センサ値取得手段で取得された実測センサ値に基づいて、前記複数のセンサで検知される複数の実測センサ値間の相関関係を求め、当該相関関係を予測モデルとして格納する予測モデル構築手段と、前記センサ値取得手段が前記複数のセンサの実測センサ値を取得した際に、当該取得した実測センサ値と前記予測モデル構築手段が格納する予測モデルとから当該複数のセンサの予測センサ値を求め、当該実測センサ値と予測センサ値との差異に基づいて、当該設備の故障兆候を検知する故障兆候検知手段とを有することを特徴とする遠隔監視システムの発明が提案されている（特許文献２参照）。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特開平０８−３２８６４４号公報

【特許文献2】特開２００５−１４９１３７号公報

【発明の概要】

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

【0006】

本発明は、被監視装置の内部状態を示す複数のパラメータの時系列データに基づいて障害予測を行う技術に関し、或るパラメータのデータ異常を原因にして他のパラメータのデータにも異常が生じる場合であっても、原因となるパラメータを効果的に特定して障害予測を精度向上に役立てることを可能にする技術を提案することを目的とする。

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

【0007】

請求項１に係る本発明は、被監視装置の内部状態を示す複数のパラメータについて、値の変化に因果関係がある一対のパラメータから成る組を複数設定した設定手段と、前記設定手段により設定された組毎に得られる、各組を成すパラメータ間における値の時系列変化の相関度の推移が最も類似する２つの組について、共通するパラメータが存在する場合には当該共通するパラメータを特定し、共通するパラメータが存在しない場合には、当該２つの組について、一方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度、及び、他方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度を調べ、その中で相関度の推移が最も類似する２つの組に共通するパラメータを特定する特定手段と、前記特定手段により特定されたパラメータに予め対応付けられた障害が被監視装置に発生する旨の通知を出力する出力手段と、を備えたことを特徴とする障害予測システムである。

【0008】

請求項２に係る本発明は、請求項１に係る本発明において、前記設定手段は、前記被監視装置の内部状態を示す複数のパラメータを少なくとも帯電電圧、トナー濃度、現像及びレーザパワーの各類に分類したパラメータ群毎に、そのパラメータ群に分類された各パラメータについて、正常状態の前記被監視装置において、前記値の変化に因果関係がある一対のパラメータからなる組を複数設定する、ことを特徴とする障害予測システムである。

【0009】

請求項３に係る本発明は、請求項１、２に係る本発明において、前記特定手段は、前記各組を成すパラメータ間の相関係数を対象範囲の期間をスライドさせた複数の期間についてそれぞれ算出し、当該算出した各期間の相関係数を時系列順に並べて相関データ列とし、当該相関データ列間の距離が近いほど相関度の推移が類似すると判断する、ことを特徴とする障害予測システムである。

【0011】

請求項４に係る本発明は、コンピュータに、被監視装置の内部状態を示す複数のパラメータについて、値の変化に因果関係がある一対のパラメータから成る組を複数設定した設定機能と、前記設定機能により設定された組毎に得られる、各組を成すパラメータ間における値の時系列変化の相関度の推移が最も類似する２つの組について、共通するパラメータが存在する場合には当該共通するパラメータを特定し、共通するパラメータが存在しない場合には、当該２つの組について、一方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度、及び、他方の組における相関度の推移と当該組を成すパラメータについて設定された他の組における相関度の推移との類似度を調べ、その中で相関度の推移が最も類似する２つの組に共通するパラメータを特定する特定機能と、前記特定機能により特定されたパラメータに予め対応付けられた障害が被監視装置に発生する旨の通知を出力する出力機能と、を実現するためのプログラムである。

【0012】

請求項５に係る本発明は、請求項４に係る本発明において、前記設定機能は、前記被監視装置の内部状態を示す複数のパラメータを少なくとも帯電電圧、トナー濃度、現像及びレーザパワーの各類に分類したパラメータ群毎に、そのパラメータ群に分類された各パラメータについて、正常状態の前記被監視装置において、前記値の変化に因果関係がある一対のパラメータからなる組を複数設定する、ことを特徴とするプログラムである。

【0013】

請求項６に係る本発明は、請求項４，５に係る本発明において、前記特定機能は、前記各組を成すパラメータ間の相関係数を対象範囲の期間をスライドさせた複数の期間についてそれぞれ算出し、当該算出した各期間の相関係数を時系列順に並べて相関データ列とし、当該相関データ列間の距離が近いほど相関度の推移が類似すると判断する、ことを特徴とするプログラムである。

【発明の効果】

【0015】

請求項１，４に係る本発明によれば、或るパラメータのデータ異常を原因にして他のパラメータのデータにも異常が生じる場合であっても、原因となるパラメータを効果的に特定して障害予測を精度向上に役立てることが可能になる。また、相関度の時系列変化が最も類似する２つの組に共通するパラメータが存在しない場合であっても、原因となるパラメータの特定を行うことができる。

【0016】

請求項２，５に係る本発明によれば、原因となるパラメータの特定に係る計算量を削減することができる。また、被監視装置における異常の発生を高い感度で検出できるようになる。

【0017】

請求項３，６に係る本発明によれば、元データの僅かな変化を敏感に捉えることができる。

【図面の簡単な説明】

【0019】

【図1】本発明の一実施形態に係る障害予測システムの構成例を示す図である。

【図2】パラメータペアの一例を示す図である。

【図3】原因パラメータ特定部による処理を説明する図である。

【図4】原因パラメータ特定部による処理を説明する図である。

【図5】原因パラメータとその異常に関連する障害とを対応付けた対応表の一例を示す図である。

【図6】本例の障害予測システムによる処理フローの一例を示す図である。

【図7】原因パラメータの特定処理に係る処理フローの一例を示す図である。

【図8】２つのパラメータの各時系列データとその相関データ列との関係の一例を示す図である。

【図9】３つのパラメータから原因パラメータを特定する例を説明する図である。

【図10】３つのパラメータから原因パラメータを特定する例を説明する図である。

【図11】サーバ装置のハードウェア構成を例示する図である。

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

【0020】

本発明の一実施形態に係る障害予測システムについて、図面を参照して説明する。
まず、本例に係る障害予測システムの説明に先立ち、その背景について説明しておく。
例えば、複写機では、機械メンテナンスのために、プロセス制御パラメータなどの機械内部の状態を示す多くのデータが稼働中に収集されている。これらの時系列データは、複雑なフィードバックループを持つシステムパラメータのデータであることが多く、或るパラメータのデータ異常を原因にして他のパラメータのデータにも異常が生じることがある。このため、時系列データの相互相関係数の推移から障害予測する際に、時系列データに異常があるパラメータが複数検出された場合には、どのパラメータが原因であったかを特定できなければ、障害予測に支障が生ずる。
そこで、本例の障害予測システムでは、以下のような構成により、原因パラメータを特定することで、障害予測の精度向上に役立てるようにしている。

【0021】

図１には、本例の障害予測システムの構成を例示してある。
本例の障害予測システムは、障害予測の対象となる複数の複写機と、各複写機から収集したデータに基づいて複写機に障害が発生することを予測するサーバ装置とを、通信網を介して有線又は無線により通信可能に接続したシステムであり、サーバ装置側に、データ収集部１、パラメータ分類部２、グループ別パラメータペア抽出部３、グループ別パラメータペア格納部４、相互相関係数算出部５、閾値比較部６、類似度算出部７、原因パラメータ特定部８、障害予測部９、を設けてある。本例では、１台のサーバ装置に上記の各機能部１〜９を設けた構成であるが、複数台のサーバ装置に分散して設けるようにしてもよい。

【0022】

データ収集部１は、各々の複写機から、その内部状態を示す複数のパラメータのデータを収集して蓄積する。収集データには、収集元の複写機を識別する装置ＩＤ、収集日時、パラメータ種別及びその値（設定値や測定値）が含まれている。この収集データに基づいて、同一の装置ＩＤ及び同一のパラメータ種別についてパラメータの値を収集日時順に並べることで、パラメータの時系列データが得られる。

【0023】

パラメータ分類部２は、複数のパラメータを、共通性を有するパラメータ同士をまとめたパラメータ群（グループ）に分類（グルーピング）する。分類は、パラメータが収集されている場所（部位）や内容などのパラメータ属性に基づいて行う。例えば、帯電電圧系、トナー濃度系、現像系、レーザパワー系等の系毎に関連するパラメータを分類する。パラメータによっては、２つ以上のパラメータ群に分類されることもある。
このように、複数のパラメータを予め分類しておくことで、後述する原因パラメータの特定に係る計算量を大幅に削減することができる。

【0024】

グループ別パラメータペア抽出部３は、パラメータ分類部２により分類されたパラメータ群別（グループ別）に、正常状態（障害及びその予兆が検出されていない状態）の複写機において時系列データの相関度が強いパラメータペア（２つのパラメータの組み合わせ）を抽出し、そのペアの情報をグループ別パラメータペア格納部４に格納する。すなわち、値の変化に因果関係がある一対のパラメータから成る組（ペア）を抽出して予め設定しておく。なお、パラメータペアの抽出は、１台の複写機（例えば、試験用の複写機）から収集された時系列データに基づいて行ってもよく、複数の複写機（例えば、運用中の複写機）から収集された時系列データに基づいて行ってもよい。

【0025】

本例では、データ収集部１により収集された各パラメータの時系列データに基づいて、正常状態における２つのパラメータの値の時系列変化の相関度として、その時系列データ間の相互相関係数（以下、相関係数）を算出し、正常状態の相関係数が閾値（本例では、０．９）以上となるパラメータペアを抽出する。
ここで、相関係数とは、２つのデータ間の相関（類似性の度合）を示す統計学的指標であり、２つのデータが全く同じ（或いは、変化の仕方が同じ）場合には、相関係数は１となり、逆に全く正反対に変化する場合には、−１となる。また、２つのデータに全く関連性が無い場合には、ゼロとなる。

【0026】

例えば、データ列（ｘ_１，ｘ_２，・・・，ｘ_ｎ）とデータ列（ｙ_１，ｙ_２，・・・，ｙ_ｎ）の相関係数は、（式１）により求めることができる。

【数1】

ここで、ｘ￣，ｙ￣は、それぞれ、ｘ＝｛ｘ_ｉ｝，ｙ＝｛ｙ_ｉ｝の相加平均である。
これは、各データの平均からのずれを表すベクトル（ｘ−ｘ￣）＝（ｘ_１−ｘ￣，・・・，ｘ_ｎ−ｘ￣），（ｙ−ｙ￣）＝（ｙ_１−ｙ￣，・・・，ｙ_ｎ−ｙ￣）のなす角の余弦である。また、この式は、共分散をそれぞれの標準偏差で割ったものに等しい。

【0027】

相関係数の計算は、注目している期間（例えば、６ヶ月）について、その期間全体に亘って一度に算出（すなわち、６ヶ月全体のデータ推移について相関係数を算出）してもよいが、その期間内で算出期間を或る程度の幅（例えば、１週間）で順次シフトしながら算出（すなわち、１週間単位のデータ推移について相関係数を算出）する方法を採用してもよく、注目している期間で連続して強い相関を示すパラメータペアを検出できればよい。なお、正常状態でどのパラメータとも強い相関が得られないパラメータは利用しない。
このように、正常状態において相関が強いパラメータをペアにしておくことで、異常発生時に高い感度で相関係数に変化が現れる効果がある。

【0028】

図２には、パラメータペアの一例を示してある。
図２（ａ）は帯電電圧系のパラメータ群におけるパラメータペアの例であり、ＰＡ＿ａ、ＰＡ＿ｂ、ＰＡ＿ｃ、・・・等のパラメータにより構成されている。図２（ｂ）はトナー濃度系のパラメータ群におけるパラメータペアの例であり、ＰＢ＿ａ、ＰＢ＿ｂ、ＰＢ＿ｃ、・・・等のパラメータにより構成されている。図２（ｃ）は現像系のパラメータ群におけるパラメータペアの例であり、ＰＣ＿ａ、ＰＣ＿ｂ、ＰＣ＿ｃ、・・・等のパラメータにより構成される。各パラメータ群は、最小で３種類以上のパラメータを含むように構成される。なお、本例では、各パラメータは、Ｙ（Ｙｅｌｌｏｗ），Ｍ（Ｍａｇｅｎｔａ），Ｃ（Ｃｙａｎ），Ｋ（Ｋｅｙ Ｐｌａｔｅ）といった色要素別に設けられており、パラメータ名の末尾に付した数値（１）〜（４）で色要素を区別している。

【0029】

相互相関係数算出部５は、グループ別パラメータペア格納部４に格納されているパラメータ群別（グループ別）のパラメータペアの情報に基づいて、対象の複写機について直近の期間（本例では、直近過去５日間）におけるパラメータペアの時系列データ間の相関度（本例では、相関係数）を算出する。なお、相関係数を算出する対象範囲の期間（算出期間）は、相関係数の変化の急峻さ（日数が少ないほど急峻）と関係が有り、当該算出期間は厳密に決められるものではなく、異常検出との関係で決めればよい。
本例の相互相関係数算出部５では、データ収集部１により新たなデータが収集されるに従い（例えば、１日毎）、パラメータ群別のパラメータペア毎に、算出期間をスライドさせながら相関係数の算出を順次行う。これによりパラメータペア毎に得られる相関係数を時系列順に並べたデータ列を、そのパラメータペアに係る相関データ列（相関度の推移を示すデータ列）とする。すなわち、本例の相互相関係数算出部５によれば、対象の複写機について、パラメータペア毎に相関データ列を得ることができる。

【0030】

ここで、相関データ列について図８を参照して説明しておく。
図８（ａ）には、パラメータＡ，Ｂの各々の時系列データ（各パラメータの値の時系列変化）を示すグラフを例示してあり、図８（ｂ）には、パラメータＡ，Ｂのペアに係る相関データ列（相関係数の推移）を示すグラフを例示してある。
図８（ａ）によれば、パラメータＡの時系列データとパラメータＢの時系列データは同様な変化をしているように見える。しかしながら、図８（ｂ）に示すように、これらの時系列データに基づいて作成される相関データ列では、元データの僅かな変化の相違が大きな変動として現れている。例えば、図８（ａ）の点線で囲った部分の相違は僅かだが、図８（ｂ）ではその相違が大きな変動で表されている。このように、相関データ列を用いることで、元データの僅かな変化を敏感に捉えることができる。

【0031】

閾値比較部６は、相互相関係数算出部５により得られるパラメータペア毎に相関データ列を常時監視して、予め定められた閾値（本例では、０．４）と比較する。そして、閾値を超えた（或いは当該閾値以上の）相関係数を相関データ列に含むパラメータペアが検出された場合に異常状態であると判定し、当該パラメータペアが属するパラメータ群（グループ）を特定する。すなわち、データ異常の原因となった原因パラメータの候補が属するパラメータ群を特定する。なお、原因パラメータが複数のパラメータ群に属する場合には、複数のパラメータ群が特定され得る。

【0032】

類似度算出部７は、閾値比較部６により特定された原因パラメータを含むパラメータ群（グループ）について、そのパラメータ群における全てのパラメータペアの相関データ列に基づいて、各相関データ列間の類似度を算出する。
本例では、相関データ列間の類似度として、予め定められた期間（例えば、５日間）における相関データ列間のユークリッド距離を用いる。すなわち、相関データ列をベクトルと見做すことで、相関データ列間（ベクトル間）の距離を、各々の相関データ列間の波形の変化の類似性を表す指標として用いることができる。この場合、相関データ列間の距離が小さいほど、両データ列の波形の変化が似ていることを意味し、距離がゼロであれば、両データ列は等しいことを意味する。

【0033】

例えば、パラメータＡ，Ｂのペアに係る相関データ列ＡＢ（ｘ_１，ｘ_２，・・・，ｘ_ｎ）と、パラメータＣ，Ｄのペアに係る相関データ列ＣＤ（ｙ_１，ｙ_２，・・・，ｙ_ｎ）との距離Ｄ_{ＡＢ，ＣＤ}は、（式２）により求めることができる。

【数2】

本例では、ｎ＝３０で計算して距離を算出している。

【0034】

原因パラメータ特定部８は、閾値比較部６により特定された原因パラメータを含むパラメータ群（グループ）について、類似度計算部７により算出された各パラメータペアの相関データ列間の類似度（本例では、距離）に基づいて、データ異常の原因となった原因パラメータを特定（推定）する。

【0035】

原因パラメータ特定部８による原因パラメータの特定処理について、図３、図４を参照して説明する。
例えば、図３に例示するように、パラメータＡ，Ｂのペアに係る相関データ列ＡＢと、パラメータＥ，Ａのペアに係る相関データ列ＥＡとの距離Ｄ_{ＡＢ，ＥＡ}が、全ての相関データ列間の距離の中で最も短いとする。このとき、その距離Ｄ_{ＡＢ，ＥＡ}に係る各パラメータペアでは、パラメータＡが共通する。この場合、パラメータＡが原因となって他のパラメータに影響を与えていると推定される。このように、最短距離となる２つの相関データ列に共通するパラメータが存在する場合には、そのパラメータがデータ異常（障害等によって発生した相関係数の低下）の原因パラメータと推定できる。

【0036】

また、例えば、図４に例示するように、パラメータＡ，Ｂのペアに係る相関データ列ＡＢと、パラメータＥ，Ｆのペアに係る相関データ列ＥＦとの距離Ｄ_{ＡＢ，ＥＦ}が、全ての相関データ列間の距離の中で最も短いとする。このとき、その距離Ｄ_{ＡＢ，ＥＡ}に係る各パラメータペアには、共通するパラメータがない。この場合、以下の手順により、条件付きで最短距離を調べ直して原因パラメータを推定する。
（手順１）一方の相関データ列ＡＢと、その対象であるパラメータＡ，Ｂのいずれかが関わる他の相関データ列との距離を求める。すなわち、パラメータＡに関して距離Ｄ_{ＡＢ，ＡＣ}、Ｄ_{ＡＢ，ＡＥ}、Ｄ_{ＡＢ，ＡＦ}を求め、パラメータＢに関して距離Ｄ_{ＡＢ，ＢＣ}、Ｄ_{ＡＢ，ＢＥ}、Ｄ_{ＡＢ，ＢＦ}を求める。
（手順２）他方の相関データ列ＥＦと、その対象であるパラメータＥ，Ｆのいずれかが関わる他の相関データ列との距離を求める。すなわち、パラメータＥに関して距離Ｄ_{ＥＦ，ＥＡ}、Ｄ_{ＥＦ，ＥＢ}、Ｄ_{ＥＦ，ＥＣ}を求め、パラメータＦに関して距離Ｄ_{ＥＦ，ＦＡ}、Ｄ_{ＥＦ，ＦＢ}、Ｄ_{ＥＦ，ＦＣ}を求める。
（手順３）上記の（手順１）及び（手順２）で求めた距離の中から最も短い距離の相関データ列に共通するパラメータを原因パラメータとして推定する。例えば、相関データ列ＡＢと相関データ列ＢＣとの距離Ｄ_{ＡＢ，ＢＣ}が最短であれば、これら２つの相関データ列に共通するパラメータＢを原因パラメータとして推定する。
なお、もともと初めに全ての組み合わせで距離計算を行うので、（手順１）及び（手順２）での距離を改めて求め直す必要は無い。

【0037】

障害予測部９は、原因パラメータ特定部８により特定された原因パラメータについて、当該原因パラメータのデータ異常に起因して発生することが予測される障害を特定（推定）する。本例では、図５に例示するように、原因パラメータとその異常に関連する障害（及び施すべき処置の内容）とを対応付けた対応表をメモリに保持しており、当該対応表を参照することで、対象の複写機に発生することが予測される障害（及び施すべき処置の内容）を推定する。
障害予測部９により推定された障害の情報は、対象の複写機の保守を担当する保守担当者などに通知される。これにより、実際に障害が発生する前に、計画的にメンテナンス等の保守作業を行うことができる。

【0038】

本例の障害予測システムによる処理について、図６に例示する処理フローを参照して説明する。
本例の障害予測システムでは、前処理として、パラメータ分類部２が、サブシステム（帯電電圧系、トナー濃度系、現像系など）単位でパラメータを分類（グルーピング）し（ステップＳ１）、グループ別パラメータペア抽出部３が、パラメータ群（グループ）別に、正常状態の相関係数が強いパラメータペアを抽出してグループ別パラメータペア格納部４に格納しておく（ステップＳ２）。これらの前処理は、パラメータ群の構成の変更や、パラメータの追加、修正、削除などが無い限り実施し直す必要は無い。

【0039】

上記の前処理が実施された後、対象の複写機毎に、以下の処理を繰り返す。
まず、相互相関係数算出部５が、直近の期間におけるパラメータペアの時系列データ間の相関データ列を生成し（ステップＳ３）、閾値比較部６が、相互相関係数算出部５により算出された相関データ列に閾値を超える相関係数が含まれるか否かを判定し（ステップＳ４）、閾値を超える相関係数を含むパラメータペアが検出された場合に異常状態であると判断して、当該パラメータペアを含むパラメータ群（グループ）を特定する（ステップＳ５）。その後、類似度算出部７及び原因パラメータ特定部８が、特定されたパラメータ群（グループ）について、原因パラメータの特定処理を行い（ステップＳ６）、障害予測部９が、特定された原因パラメータのデータ異常に起因して発生することが予測される障害を特定（推定）する（ステップＳ７）。

【0040】

原因パラメータの特定処理について、図７の処理フローを参照して更に説明する。
まず、特定されたパラメータ群（グループ）に属するパラメータの数が２つか否か（すなわち、パラメータペアが１つか否か）を判定する（ステップＳ１１）。
ステップＳ１１において、パラメータの数が２つと判定された場合には、そのいずれか一方が原因パラメータであると特定する（ステップＳ１２）。
一方、ステップＳ１１において、パラメータの数が２つ以上と判定された場合には、パラメータ群における全てのパラメータペアの相関データ列に基づいて、各相関データ列間の距離（類似度）を算出し（ステップＳ１３）、最短距離となる２つの相関データ列に共通するパラメータが存在するか否かを判定する（ステップＳ１４）。
ステップＳ１４において、最短距離となる２つの相関データ列に共通するパラメータが存在すると判定された場合には、その共通パラメータを原因パラメータとして特定する（ステップＳ１６）。
一方、ステップＳ１４において、最短距離となる２つの相関データ列に共通するパラメータが存在しないと判定された場合には、この２つの相関データ列とその対象である４つのパラメータが関わる他の相関データ列との距離を調べ、その中で最短距離となる２つの相関データ列を特定し（ステップＳ１５）、その共通パラメータを原因パラメータとして特定する（ステップＳ１６）。すなわち、条件付きで最短距離となる相関データ列における共通パラメータを原因パラメータとして特定する。

【0041】

ここで、一例として、パラメータＡ，Ｂ，Ｃの３つのパラメータからなるパラメータ群の中から原因パラメータを特定する場合について、図９、図１０を参照して説明する。
まず、相互相関係数算出部５が、図９（ａ）に示すように、パラメータＡ，Ｂのペアに係る相関データ列ＡＢ、パラメータＢ，Ｃのペアに係る相関データ列ＢＣ、パラメータＣ，Ａのペアに係る相関データ列ＣＡをそれぞれ算出する。
その結果、相関データ列ＡＢについて図１０（ａ）のデータ推移が得られ、相関データ列ＢＣについて図１０（ｂ）のデータ推移が得られ、相関データ列ＣＡについて図１０（ｃ）のデータ推移が得られたとする。
この場合、閾値比較部６により、相関データ列ＡＢ及び相関データ列ＢＣに、閾値を超える相関係数が含まれることが検出され、そのパラメータ群について原因パラメータを特定する処理に移行する。

【0042】

原因パラメータの特定では、まず、類似度算出部７が、図９（ｂ）に示すように、相関データ列ＡＢと相関データ列ＢＣとの距離Ｄ_{ＡＢ，ＢＣ}、相関データ列ＢＣと相関データ列ＣＡとの距離Ｄ_{ＢＣ，ＣＡ}、相関データ列ＣＡと相関データ列ＡＢとの距離Ｄ_{ＣＡ，ＡＢ}をそれぞれ算出する。
その結果、図９（ｃ）に示すように、距離Ｄ_{ＡＢ，ＢＣ}＝０．０８５、距離Ｄ_{ＢＣ，ＣＡ}＝０．２５１、距離Ｄ_{ＣＡ，ＡＢ}＝０．１９１が得られたとする。
この場合、各距離の関係は、距離Ｄ_{ＡＢ，ＢＣ}＜Ｄ_{ＣＡ，ＡＢ}＜Ｄ_{ＢＣ，ＣＡ}であるため、原因パラメータ特定部８により、最短の距離Ｄ_{ＡＢ，ＢＣ}について共通するパラメータであるパラメータＢが原因パラメータとして特定される。

【0043】

以上のように、本例の障害予測システムでは、複数のパラメータからなるパラメータ群毎に、相互相関係数算出部５が、各パラメータペアについて相関データ列を生成し、閾値比較部６が、相関データ列に閾値を超える相関係数が含まれるか否かを判定し、閾値を超える相関係数を含むパラメータペアが検出された場合に、類似度算出部７が、そのパラメータが属するパラメータ群について各パラメータペアの相関データ列間の距離を算出し、原因パラメータ特定部８が、距離が最短（若しくは、条件付きで距離が最短）となる２つの相関データ列に共通するパラメータを原因パラメータとして特定し、障害予測部９が、当該原因パラメータのデータ異常に起因する障害の発生を予測するようにした。
このような構成によれば、原因パラメータを効果的に特定することができるため、そのパラメータのデータ異常に起因する障害の発生予測を精度良く行うことが可能になる。

【0044】

ここで、本例の障害予測システムは、複写機について障害予測を行うものであるが、他の装置について障害予測を行うようにしてもよく、例えば、複写機と同様に、用紙等の記録材に画像を形成して出力する画像形成機能を備えた、プリンタ（文書印刷装置）やファクシミリ装置（文書転送装置）、或いはこれらの装置の機能を複合的に備えた複合機など、他の種別の画像形成装置について障害予測を行うようにしてもよい。また、このような画像形成装置への適用に限定するものではなく、複雑なフィードバックループを持つ多数のパラメータを用いて制御される他の装置にも適用することができる。

【0045】

また、本例の障害予測システムでは、被監視装置（本例では、複写機）とは別体のサーバ装置により障害予測を行う構成にしてあるが、各々の被監視装置が自己について障害予測を行う構成にしてもよい。この場合には、被監視装置側に、データ収集部１、相互相関係数算出部５、閾値比較部６、類似度算出部７、原因パラメータ特定部８、障害予測部９などの機能を設ければよい。

【0046】

図１１には、本例の障害予測システムにおけるサーバ装置のハードウェア構成を例示してある。
本例では、各種演算処理を行うＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）２１、ＣＰＵ２１の作業領域となるＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）２２や基本的な制御プログラムを記録したＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）２３等の主記憶装置、本発明の一実施形態に係るプログラムや各種データを記憶するＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）２４等の補助記憶装置、各種情報を表示出力するための表示装置及び操作者により入力操作に用いられる操作ボタンやタッチパネル等の入力機器とのインタフェースである入出力Ｉ／Ｆ２５、他の装置との間で有線又は無線により通信を行うインタフェースである通信Ｉ／Ｆ２６、等のハードウェア資源をサーバ装置のコンピュータが有している。
そして、本発明の一実施形態に係るプログラムを補助記憶装置２４等から読み出してＲＡＭ２２に展開し、これをＣＰＵ２１により実行させることで、上述した各機能部をサーバ装置のコンピュータ上に実現している。

【0047】

なお、本発明の一実施形態に係るプログラムは、例えば、当該プログラムを記憶したＣＤ−ＲＯＭ等の外部記憶媒体から読み込む形式や、通信網等を介して受信する形式などにより、本例に係るサーバ装置のコンピュータに設定される。
また、本例のようなソフトウェア構成により各機能部を実現する態様に限られず、それぞれの機能部を専用のハードウェアモジュールで実現するようにしてもよい。

【符号の説明】

【0048】

１：データ収集部、 ２：パラメータ分類部、 ３：グループ別パラメータペア抽出部、 ４：グループ別パラメータペア格納部、 ５：相互相関係数算出部、 ６：閾値比較部、 ７：類似度算出部、 ８：原因パラメータ特定部、 ９：障害予測部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】