(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-10-16
(45)【発行日】2023-10-24
(54)【発明の名称】情報処理装置、評価方法、および評価プログラム
(51)【国際特許分類】
G06N 20/00 20190101AFI20231017BHJP
G05B 23/02 20060101ALI20231017BHJP
【ＦＩ】
G06N20/00
G05B23/02 T
【請求項の数】
7
(21)【出願番号】P 2022522488
(86)(22)【出願日】2020-05-15
(86)【国際出願番号】 JP2020019549
(87)【国際公開番号】W WO2021229815
(87)【国際公開日】2021-11-18
【審査請求日】2022-09-28
(73)【特許権者】
【識別番号】000005223
【氏名又は名称】富士通株式会社
(74)【代理人】
【識別番号】100121083
【弁理士】
【氏名又は名称】青木 宏義
(74)【代理人】
【識別番号】100138391
【弁理士】
【氏名又は名称】天田 昌行
(74)【代理人】
【識別番号】100074099
【弁理士】
【氏名又は名称】大菅 義之
(72)【発明者】
【氏名】雨宮 智
(72)【発明者】
【氏名】阿部 秀城
【審査官】山田 辰美
(56)【参考文献】
【文献】特開２０１９－１５３０１８（ＪＰ，Ａ）
【文献】特開２０１９－０９６０１４（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０６Ｎ ２０／００
Ｇ０５Ｂ ２３／０２
(57)【特許請求の範囲】
【請求項1】
評価対象物の状態を計測する複数のセンサから得られた時系列の計測データのうちで、基準ウィンドウ期間に計測されたデータから生成したスパース構造モデルに対する、前記基準ウィンドウ期間よりも後のウィンドウ期間に計測されたデータから生成したスパース構造モデルのセンサ毎の変化度を成分とする変化スコアを、学習期間に含まれる複数のウィンドウ期間のそれぞれで求めて複数の変化スコアを取得する取得部と、前記複数の変化スコアをクラスタリングしてクラスタの位置を特定する特定部と、
評価対象期間に含まれる少なくとも１つのウィンドウ期間で求めた少なくとも１つの変化スコアを、前記クラスタに追加することで生じる前記クラスタの位置の移動に基づいて、前記評価対象物の劣化を評価する評価部と、
を含む、情報処理装置。
【請求項2】
前記評価部は、前記クラスタに追加することで生じる前記クラスタの位置の移動距離、移動速度、および移動加速度の少なくとも１つに基づいて、前記評価対象物の劣化を評価する、請求項１に記載の情報処理装置。【請求項3】
前記特定部は、前記クラスタリングにより得られた複数のクラスタの位置を特定し、前記評価部は、前記複数のクラスタのうちの少なくとも１つのクラスタの位置の移動に基づいて、前記評価対象物の劣化を評価する、
請求項１または２に記載の情報処理装置。
【請求項4】
前記評価部は、前記クラスタの位置の移動が所定の条件を満たす場合に、前記評価対象物に対する保守作業を促す情報を出力する、請求項１から３のいずれか１項に記載の情報処理装置。【請求項5】
前記評価部は、前記クラスタの位置の移動が第１の条件を満たす場合に、前記評価対象物に対する第１の種別の劣化の保守作業を促す情報を出力し、
前記クラスタの位置の移動が第２の条件を満たす場合に、前記評価対象物に対する第２の種別の劣化の保守作業を促す情報を出力する、
請求項１から３のいずれか１項に記載の情報処理装置。
【請求項6】
評価対象物の状態を計測する複数のセンサから得られた時系列の計測データのうちで、基準ウィンドウ期間に計測されたデータから生成したスパース構造モデルに対する、前記基準ウィンドウ期間よりも後のウィンドウ期間に計測されたデータから生成したスパース構造モデルのセンサ毎の変化度を成分とする変化スコアを、学習期間に含まれる複数のウィンドウ期間のそれぞれで求めて複数の変化スコアを取得し、前記複数の変化スコアをクラスタリングしてクラスタの位置を特定し、
評価対象期間に含まれる少なくとも１つのウィンドウ期間で求めた少なくとも１つの変化スコアを、前記クラスタに追加することで生じる前記クラスタの位置の移動に基づいて、前記評価対象物の劣化を評価する、
ことを含む、情報処理装置が実行する評価方法。
【請求項7】
評価対象物の状態を計測する複数のセンサから得られた時系列の計測データのうちで、基準ウィンドウ期間に計測されたデータから生成したスパース構造モデルに対する、前記基準ウィンドウ期間よりも後のウィンドウ期間に計測されたデータから生成したスパース構造モデルのセンサ毎の変化度を成分とする変化スコアを、学習期間に含まれる複数のウィンドウ期間のそれぞれで求めて複数の変化スコアを取得し、前記複数の変化スコアをクラスタリングしてクラスタの位置を特定し、
評価対象期間に含まれる少なくとも１つのウィンドウ期間で求めた少なくとも１つの変化スコアを、前記クラスタに追加することで生じる前記クラスタの位置の移動に基づいて、前記評価対象物の劣化を評価する、
処理を情報処理装置に実行させる評価プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、情報処理装置、評価方法、および評価プログラムに関する。
【背景技術】
【0002】
例えば、産業機器、産業プラント、情報処理システム、並びにその他の機器および設備などの異常や劣化を検出するための技術が開発されている。例えば、機器や設備などにセンサを設置し、センサから得られる時系列のデータに基づいて、機器や設備などの異常を検出したり、劣化の度合いを推定したりすることが考えられる。
【0003】
これに関し、例えば、機器および設備などの異常や劣化を監視したり診断したりするための技術が知られている（例えば、特許文献１から特許文献３）。
【先行技術文献】
【特許文献】
【0004】
【文献】特開２０１６－９５７５１号公報
【文献】特開２００４－８６８９７号公報
【文献】特開２０１６－１０５６８７号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかしながら、センサで計測した計測データには機器の異常や劣化以外にも、例えば、季節の移り変わりに起因する計測値の変動といった機器の異常や劣化とは相関の薄い事象に基づく変動が含まれている。こうした要因は無数に存在し得るため、そのすべてを特定して計測データから排除することは難しい。そのため、多様な要因の影響を受けて変動し得る計測データから、機器や設備などの劣化を評価することのできる技術の提供が望まれている。
【0006】
１つの側面では、本発明は、多様な要因の影響を受けて変動し得る計測データから、機器や設備などの劣化を評価することを目的とする。
【課題を解決するための手段】
【0007】
本発明の一つの態様の情報処理装置は、評価対象物の状態を計測する複数のセンサから得られた時系列の計測データのうちで、基準ウィンドウ期間に計測されたデータから生成したスパース構造モデルに対する、基準ウィンドウ期間よりも後のウィンドウ期間に計測されたデータから生成したスパース構造モデルのセンサ毎の変化度を成分とする変化スコアを、学習期間に含まれる複数のウィンドウ期間のそれぞれで求めて複数の変化スコアを取得する取得部と、複数の変化スコアをクラスタリングしてクラスタの位置を特定する特定部と、評価対象期間に含まれる少なくとも１つのウィンドウ期間で求めた少なくとも１つの変化スコアを、クラスタに追加することで生じるクラスタの位置の移動に基づいて、評価対象物の劣化を評価する評価部と、を含む。
【発明の効果】
【0008】
多様な要因の影響を受けて変動し得る計測データから、機器や設備などの劣化を評価することができる。
【図面の簡単な説明】
【0009】
【図1】ボイラーチューブの複数の種類の劣化を例示する図である。
【図2】実施形態に係る情報処理装置のブロック構成を例示する図である。
【図3】実施形態に係る劣化の評価処理の流れを例示する図である。
【図4】実施形態に係る基準となるウィンドウ期間のスパース構造モデルを生成する処理の動作フローを例示する図である。
【図5】実施形態に係る劣化の評価処理の動作フローを例示する図である。
【図6】実施形態に係る劣化の評価処理の流れを例示する図である。
【図7】複数のクラスタの移動を例示する図である。
【図8】実施形態に係る判定情報を例示する図である。
【図9】実施形態に係る情報処理装置を実現するためのコンピュータのハードウェア構成を例示する図である。
【発明を実施するための形態】
【0010】
以下、図面を参照しながら、本発明のいくつかの実施形態について詳細に説明する。なお、複数の図面において対応する要素には同一の符号を付す。
【0011】
上述のように、例えば、機器や設備などの状態を計測するセンサを設置し、センサから得られる時系列の計測データに基づいて、機器や設備などの劣化の状態を推定することが考えられる。なお、以下では、機器や設備などの劣化を評価する対象物を、評価対象物と呼ぶことがある。
【0012】
一例として、火力発電所のボイラーチューブ群の劣化および余寿命を管理したいとする。この場合に、例えば、光ファイバ温度センサなどの温度センサや振動センサなどのセンサを、ボイラーチューブ群の温度や振動などの状態を計測できる配置で設置する。そして、センサで計測したボイラーチューブ群の各所の温度や振動などの計測データの変化から火力発電所のボイラーチューブ群の劣化を評価して、余寿命などを推定することが考えられる。そして、例えば、劣化を評価したり、余寿命を推定したりするがことができれば、点検、交換、および修理といった保守作業の実行タイミングを最適化することができる。また、例えば、ボイラーチューブのリークに起因する計画外の稼働停止といった状況を回避することができる。
【0013】
しかしながら、センサの計測データの変化には、評価対象物の劣化以外の要因に起因する情報も含まれることがある。一例として、センサで１年間の時系列の計測データを取得したとする。この場合、計測データは、例えば、昼夜の差、季節および天気などの差に起因して変動し得る。また、評価対象物が、高速稼働、低速稼働、および停止というように稼働状態に複数の稼働モードを含む場合、それらの稼働モードの変更によってもセンサの計測データは変動し得る。更には、評価対象物の周辺に配置されているエアーコンディショナーの稼働および停止というように、周辺設備の状態の変化によってもセンサの計測データは変動し得る。
【0014】
こうした評価対象物の劣化とは相関の弱い要因に基づく計測データの変動は、評価対象物の劣化を評価する上で排除することが望ましい。しかしながら、こうした要因のすべてを特定することは困難であり、評価対象物の劣化とは相関の弱い要因に起因する変動を計測データから排除することは現実的には難しい。また、排除することができたとしても、評価対象物の劣化とは相関の弱い要因に応じた変動を計測データから取り除くには大変な手間がかかる。
【0015】
そのため、多様な要因の影響を受けて変動し得るセンサの計測データから、評価対象物の劣化を評価することのできる技術の提供が望まれている。
【0016】
また更に、例えば、評価対象物の劣化の形態に、複数の種類が存在することもある。図１は、ボイラーチューブの複数の種類の劣化を例示する図である。図１では、劣化の要因と、劣化の分類と、劣化の内容とが対応づけられている。ボイラーチューブの劣化要因には、クリープ、疲労、腐食・摩耗などがある。また、それぞれの劣化要因について、その発生の仕方などに応じた劣化の分類が示されている。例えば、クリープは、高温下での長時間にわたる応力作用のもとで低速度の変形が起こる長時間クリープと、スケールによる蒸気の流れの障害などによって過熱されたときに急激な変形が生じる短時間クリープとに分類することができる。この様に、ボイラーチューブに発生する劣化といっても多岐にわたる。そして、劣化の種類のよって、センサの計測データへの劣化の表れ方も異なり得る。
【0017】
そのため、計測データから、評価対象物の複数の劣化の種別に応じて劣化を評価することのできる技術の提供も望まれている。以下、実施形態に係る評価対象物の劣化の評価について説明する。
【0018】
図２は、実施形態に係る情報処理装置２００のブロック構成を例示する図である。情報処理装置２００は、例えば、サーバーコンピュータ、パーソナルコンピュータ（ＰＣ）、モバイルＰＣ、タブレット端末などのコンピュータであってよい。情報処理装置２００は、例えば、制御部２０１、記憶部２０２、および通信部２０３を含む。制御部２０１は、例えば取得部２１１、特定部２１２、および評価部２１３などを含み、また、その他の機能部を含んでもよい。情報処理装置２００の記憶部２０２は、例えば、後述する判定情報８００などの情報を記憶している。通信部２０３は、例えば、制御部２０１の指示に従ってセンサなどの他の装置と通信してよい。これらの各部の詳細および記憶部２０２に格納されている情報の詳細については後述する。
【0019】
図３は、実施形態に係る劣化の評価処理の流れを例示する図である。図３（ａ）には、評価対象物の状態を計測する複数のセンサから取得された計測データが示されている。情報処理装置２００の制御部２０１は、例えば、複数のセンサから計測データを直接収集してもよいし、別の通信装置から複数のセンサで計測された計測データを取得してもよい。
【0020】
続いて、制御部２０１は、例えば、計測データに前処理を実行してよい。前処理では、例えば、ノイズとなり得る望ましくないデータの除去などが行われる。例えば、評価対象物を清掃している時、点検している時、起動したばかりでまだ動作が安定していない時の計測データのように通常の安定した稼働状態とは異なる状態でのデータを含む期間の計測データを劣化の評価に持ち込むと、ノイズとなり得る。そのため、制御部２０１は、これらの非通常稼働時の計測データを予め除去する前処理を実行してもよい。なお、別の実施形態では、制御部２０１は、以下に述べるウィンドウ期間を非通常稼働時の計測データを含まないように配置してもよい。
【0021】
続いて、制御部２０１は、計測データに複数のウィンドウ期間を設定する。ウィンドウ期間は、例えば、所定の長さの期間であってよく、一実施形態では制御部２０１は、ウィンドウ期間の単位で計測データを処理してよい。
【0022】
図３（ａ）では、時系列にウィンドウ期間１～ウィンドウ期間ｎのｎ個のウィンドウ期間が設定されている。ここで、時系列の最初のウィンドウ期間１は、基準となるウィンドウ期間として用いてよい。基準となるウィンドウ期間は、例えば、評価対象物の劣化が少なく、評価対象物が正常に稼働している期間に設定されてよい。一例では、基準となるウィンドウ期間は、評価対象物を設置または交換し、稼働を開始したばかりの初期の時期に設定されていてよい。また、複数のウィンドウ期間の他のウィンドウ期間は、例えば、基準となるウィンドウ期間の後に設定されていてよい。
【0023】
そして、制御部２０１は、基準としたウィンドウ期間１の計測データ１と、他の複数のウィンドウ期間の計測データとを取得する。図３（ｂ）では、ウィンドウ期間２の計測データ２が取得された場合を例示している。
【0024】
続いて、制御部２０１は、取得した複数のウィンドウ期間の計測データのそれぞれについて、ウィンドウ期間におけるセンサ間の関係をスパース構造モデルとしてモデリングする。なお、スパース構造モデルは、例えば、多変量データの中から変数同士の本質的な繋がりを抽出して得られる構造を持つモデルであってよい。スパース構造モデルは、例えば、全ての変数同士が何らかの相関関係で繋がっている「密な構造」ではない構造であってよく、弱い相関関係を省略しつつ強い相関関係のみを残した結果として得られる「疎な構造」であってよい。スパース構造モデルは、例えば、行列表現された相関関係（例えば、多変量データの分散共分散行列に対する精度行列）において、ほとんどの成分が０の疎行列（スパース行列）であるような構造を持ったモデルで表すことができる。例えば、上述のウィンドウ期間における計測データに基づくセンサ間の関係をモデル化したスパース構造モデルは、ガウシアングラフィカルモデル（ＧＧＭ：Gaussian Graphical Model）、およびグラフィカルラッソ（Graphical lasso）などのアルゴリズムを用いて生成することができる。図３（ｃ）には、スパース構造モデルがグラフ構造で例示されている。なお、スパース構造モデルは、イジングモデル（Ising Model）、ポッツモデル（Potts Model）などのその他の表現で示されてもよい。
【0025】
そして、制御部２０１は、基準としたウィンドウ期間１の計測データ１から得たスパース構造モデルと、他のウィンドウ期間の計測データから得たスパース構造モデルとの間で生じた関係の崩れを変化スコアとして検出する。例えば、制御部２０１は、基準としたウィンドウ期間１と、ウィンドウ期間２のスパース構造モデルからセンサ毎の変化度を成分とする変化スコアを算出してよい（図３（ｃ））。なお、変化スコアは、基準としたウィンドウ期間１と、他のウィンドウ期間のスパース構造モデルから求めることができるため、ウィンドウ期間がｎ個であれば、ｎ－１個の変化スコアを得ることができる。また、変化スコアは、例えば、ＫＬダイバージェンス（Kullback-Leibler divergence）、ＪＳダイバージェンス（Jensen-Shannon divergence）などを用いて分布間の差を計算することで求めることができる。
【0026】
続いて、制御部２０１は、得られた複数の変化スコアをクラスタリングする。例えば、変化スコアをＫＬダイバージェンスで求める場合、変化スコアは複数のセンサの各センサごとの変化度を成分とするベクトルで表される。そのため、変化スコアは、変化スコアの成分のそれぞれを軸とする多次元空間上に点としてプロットすることができる。図３（ｄ）では、変化スコアの点を２次元でプロットした例が示されているが、実際にはセンサ数分の次元を有する多次元空間に変化スコアの点がプロットされていてよい。そして、制御部２０１は、多次元空間上にプロットされた複数の変化スコアと対応する点をクラスタリングしてよい。例えば、制御部２０１は、混合ガウスモデルを用いて学習を行い、変化スコアのクラスタリングを実行してよい。或いは、制御部２０１は、例えば、ＤＢＳＣＡＮ（Density-based spatial clustering of applications with noise）、Ｋ平均法（k-means）などのその他のアルゴリズムを用いてクラスタリングを実行してもよい。図３（ｄ）では、クラスタリングの結果得られた２つのクラスタが示されている。
【0027】
続いて、制御部２０１は、クラスタリングに用いた複数のウィンドウ期間よりも後の計測データに、更に複数のウィンドウ期間を設け、各ウィンドウ期間と対応する変化スコアの点をクラスタに追加してゆく。なお、追加する変化スコアの点がクラスタリングされるクラスタは、例えば、追加する変化スコアの点と最も距離の近いクラスタであってよい。或いは、別の実施形態では、追加する変化スコアの点がクラスタリングされるクラスタは、再度クラスタリングを実行して決定するなどその他の手法で決定されてもよい。また、制御部２０１は、クラスタに点を追加する際に、例えば、そのクラスタに属する変化スコアの点のうちウィンドウ期間が最も古い点を除去するなど、クラスタへの点の追加に応じて、クラスタに属する古い点を除去してもよい。
【0028】
そして、この様にクラスタに新たな変化スコアの点を追加してゆくと、図３（ｅ）に示すように、クラスタの位置が移動する。実施形態では制御部２０１は、このクラスタの位置の移動に基づいて、評価対象物の劣化を評価する。例えば、制御部２０１は、クラスタの位置の移動距離、移動速度、および移動加速度などに基づいて、評価対象物の劣化を評価してよい。クラスタの移動距離は、例えば、稼働初期の計測データに基づいて特定したクラスタと、その後の計測データに基づいて特定した移動後のクラスタとの間の距離であってよい。クラスタ間の距離は、一例では、クラスタの重心間の距離で表すことができる。また、図３（ｅ）に示すように、複数のクラスタで移動が起こる場合には、制御部２０１は、複数のクラスタのそれぞれでの移動距離、移動速度、および移動加速度などの値の平均および分散などの統計値に基づいて、評価対象物の劣化を評価してもよい。
【0029】
上述のように、センサで計測される計測データには、評価対象物の劣化以外の要因に起因する変動も含まれ得る。しかしながら、劣化以外の要因による計測値の変動には、元の状態に戻り得る可逆的な変化を示す要因が多く含まれている。例えば、昼夜の差、天気の差、周辺機器が稼働しているか否かなどに起因するセンサの計測値の変動は、可逆的な変化であり、長期的に見れば元の値に戻り得る変動である。そのため、これらの可逆的な変化を示す要因による計測値の変動は、クラスタ分けには寄与し得るが、長期的にみればクラスタの位置をあまり変化させない成分であることが推定される。
【0030】
一方で、評価対象物に生じる劣化は多くの場合、不可逆な現象であり、劣化に起因して一度発生したセンサの計測値の変動は、もとに戻る可能性が低い。そのため、上述のように、変化スコアをクラスタリングし、例えば、数カ月、数年、および数十年などの長期的な期間にわたってクラスタの位置の移動を追うことで、クラスタの移動に基づいて評価対象物の劣化を評価することが可能である。そして、このようにクラスタの移動に基づいて評価対象物の劣化を評価すると、劣化以外の要因に起因する計測データの変動の影響を抑えて評価対象物の劣化を評価することができる。
【0031】
以下、実施形態に係る劣化の評価処理について更に詳細に説明する。
【0032】
なお、以下では、評価対象物の状態を計測する複数のセンサの或るウィンドウ期間ｊにおける計測データを含むデータ行列Ｘｊを式１で表すものとする。
Ｘｊ＝｛ｘ１，ｘ２，…，ｘｋ，…，ｘＮ｝ ・・・式１
Ｘｊ＝｛ｘ１，ｘ２，…，ｘｋ，…，ｘＮ｝ ・・・式１
【0033】
ここで、ウィンドウ期間は、所定の長さの期間であってよく、ｘ１，ｘ２，…，ｘｋ，…，ｘＮ（ここで、ｋ＝１～Ｎ）は、例えば、ウィンドウ期間の各時点における複数のセンサの計測値を要素とするベクトルである。例えば、ｘｋは、時点ｋにおける複数のセンサの各センサの計測値を要素とするベクトルであり、以下の式で表すことができる。
ｘｋ＝（ｘｋ１，ｘｋ２，…，ｘｋｉ，…，ｘｋＭ）Ｔ ・・・式２
ｘｋ＝（ｘｋ１，ｘｋ２，…，ｘｋｉ，…，ｘｋＭ）Ｔ ・・・式２
【0034】
ここで、Ｍはセンサの総数である。また、ｉは、ｉ＝１～Ｍであり、センサを識別するための番号を表す。従って、ｘｋ１，ｘｋ２，…，ｘｋｉ，…，ｘｋＭは、ウィンドウ期間内の時点ｋにおける複数のセンサの各センサｉの計測値である。
【0035】
図４は、実施形態に係る基準となるウィンドウ期間のスパース構造モデルを生成する処理の動作フローを例示する図である。制御部２０１は、例えば、図４の基準となるウィンドウ期間のスパース構造モデルを生成する動作フローの実行指示が入力されると、図４の処理を開始してよい。
【0036】
Ｓ４０１において制御部２０１は、評価対象物の情報を計測する複数のセンサにより計測された計測データから、基準となるウィンドウ期間に計測された計測値を含むデータ行列Ｘrefを取得する。なお、以下では基準となるウィンドウ期間を、例えば、基準ウィンドウ期間と呼ぶことがある。基準ウィンドウ期間は、例えば、評価対象物の劣化が少なく、評価対象物が正常に稼働している期間に設定されてよい。例えば、制御部２０１は、記憶部２０２に記憶されている計測データからデータ行列Ｘrefを読み出してよい。
【0037】
Ｓ４０２において制御部２０１は、基準ウィンドウ期間におけるセンサの計測データを含むデータ行列Ｘrefからスパース構造モデルを生成し、本動作フローは終了する。例えば、制御部２０１は、基準ウィンドウ期間の複数のセンサの各センサの計測データに基づき分散共分散行列を生成し、生成した分散共分散行列をスパース化してスパース構造モデルを生成してよい。一例では、スパース構造モデルは、ＧＧＭにより導出することができる。以下では、制御部２０１が、データ行列Ｘrefからスパース構造モデルとしてＧＧＭ：ｐ（ｘ｜Ｘref）を導出する場合の例を述べる。しかしながら、実施形態はこれに限定されるものではなく、その他の手法でスパース構造モデルが生成されてもよい。
【0038】
以上において、基準ウィンドウ期間の計測データから得られたスパース構造モデル（例えば、ＧＧＭ：ｐ（ｘ｜Ｘref））を、評価対象物が劣化していない状態の基準となるモデルとして用い、変化スコアを計算する。
【0039】
変化スコアは、例えば、上述の基準ウィンドウ期間に計測された複数のセンサの計測データに対して、基準ウィンドウ期間よりも後のウィンドウ期間に計測された複数のセンサの計測データがどれくらい変化しているかを評価したスコアであってよい。変化スコアは、基準ウィンドウ期間のスパース構造モデルと、評価対象のウィンドウ期間のスパース構造モデルとの間で、分布間の差を計算することで求めることができる。分布間の差は、例えば、ＫＬダイバージェンス、ＪＳダイバージェンスなどを用いて計算することができる。以下では、ＫＬダイバージェンスを用いて変化スコアを求める場合を例に実施形態を説明する。
【0040】
【0041】
Ｓ５０１において制御部２０１は、クラスタリングを行うための学習期間における複数のウィンドウ期間のデータ行列Ｘｊ（ｊ＝１…Ｌ）を取得する。なお、例えば、ｊはウィンドウ期間を識別する番号であってよく、Ｌは学習に用いるウィンドウ期間の数であってよい。例えば、制御部２０１は、記憶部２０２に予め記憶されている計測データからデータ行列Ｘｊを読み出してよい。
【0042】
Ｓ５０２において制御部２０１は、各ウィンドウ期間ｊごとにスパース構造モデルを生成する。例えば、制御部２０１は、データ行列ＸｊからＧＧＭ：ｐ（ｘ｜Ｘｊ）を導出してよい。
【0043】
Ｓ５０３において制御部２０１は、各ウィンドウ期間ｊ（ｊ＝１…Ｌ）のそれぞれにおいて、複数のセンサの各センサｉごとに変化度ａｉを算出する。制御部２０１は、例えば、以下の式３に示すＫＬダイバージェンスの式を用いて、基準ウィンドウ期間のＧＧＭ：ｐ（ｘ｜Ｘref）と、評価対象のウィンドウ期間ｊのＧＧＭ：ｐ（ｘ｜Ｘｊ）との間でセンサｉごとの変化度ａｉを算出してよい。そして、制御部２０１は、ウィンドウ期間ｊに対して、各センサｉごとの変化度ａｉを成分とする変化スコアＡｊ＝｛ａ１，ａ２，…，ａｉ，…，ａＭ｝を取得する。
【数1】
【0044】
Ｓ５０４において制御部２０１は、得られた変化スコアＡｊをクラスタリングする。例えば、制御部２０１は、変化スコアＡｊの成分のそれぞれを軸とする多次元空間上に、各ウィンドウ期間ｊごとの変化スコアＡｊと対応する点がプロットしてよい。図６は、実施形態に係る劣化の評価処理の流れを例示する図であり、図６（ａ）に示すように、制御部２０１は、各ウィンドウ期間ｊのそれぞれについて得られた変化スコアＡｊを、多次元空間上にプロットしてよい。なお、図６では、変化スコアの２つの成分を軸として平面に点がプロットされている例が示されているが、実際には変化スコアＡｊの成分のそれぞれを軸とする多次元空間上に変化スコアＡｊと対応する点がプロットされてよい。
【0045】
そして、制御部２０１は、プロットされた複数の変化スコアの点をクラスタリングしてよい。一例では、制御部２０１は、混合ガウスモデルを用いて学習を行い、複数の変化スコアの点に対してクラスタリングを実行してよい。
【0046】
なお、制御部２０１は、混合ガウスモデルにおいてクラスタリングに用いるクラスタの数を、例えば、赤池情報量基準（ＡＩＣ）およびベイズ情報量基準（ＢＩＣ）などの指標を用いて決定してよい。また、実施形態は、混合ガウスモデルによるクラスタリングに限定されるものではなく、制御部２０１は、ＤＢＳＣＡＮ、Ｋ平均法などのその他のアルゴリズムでクラスタリングを実行してもよい。
【0047】
図６（ｂ）は、変化スコアのクラスタリングを例示しており、変化スコアが２つのクラスタに分類されている。
【0048】
Ｓ５０５において制御部２０１は、クラスタリングモデルから、各クラスタの位置を特定する。例えば、制御部２０１は、クラスタの位置として、各クラスタを代表する代表点μｒ（ｒ＝１…ｓ，ｓはクラスタ数）の位置を算出してよい。クラスタを代表する代表点には、例えば、クラスタの重心を用いることができる。図６（ｂ）では、クラスタの代表点として、重心が示されている。なお、クラスタを代表する代表点は、重心以外のその他のクラスタ内の点を代表する点であってもよい。
【0049】
以上のＳ５０５までの処理により、制御部２０１は、例えば、クラスタの初期の位置を特定することができる。そして、続く、Ｓ５０６以降の処理において、制御部２０１は、例えば、評価対象物の劣化の状態を評価する対象の期間（以下、評価対象期間と呼ぶことがある）におけるセンサの計測データから評価対象物の劣化を評価する。
【0050】
Ｓ５０６において制御部２０１は、評価対象期間の計測データがあるか否かを判定する。評価対象期間は、例えば、評価対象物の劣化を評価する対象とする期間であってよく、クラスタリングを行うための学習期間の後に計測された計測データを含む期間であってよい。評価対象期間の計測データが無い場合（Ｓ５０６がＮＯ）、フローはＳ５０６の処理を繰り返す。一方、評価対象期間の計測データがある場合（Ｓ５０６がＹＥＳ）、フローはＳ５０７に進む。
【0051】
続く、Ｓ５０７において制御部２０１は、評価対象期間の計測データからデータ行列Ｘｏを取得する。ｏは、ｏ＝１…Ｑであってよく、評価対象期間の計測データに含まれるウィンドウ期間を識別する番号であってよい。また、Ｑは、評価対象期間の計測データに含まれるウィンドウ期間の数であってよい。例えば、制御部２０１は、記憶部２０２に予め記憶されている計測データからデータ行列Ｘｏを読み出してよい。
【0052】
Ｓ５０８において制御部２０１は、各ウィンドウ期間ｏごとのデータ行列Ｘｏからスパース構造モデルを生成する。例えば、制御部２０１は、データ行列ＸｏからＧＧＭ：ｐ（ｘ｜Ｘｏ）を導出してよい。
【0053】
Ｓ５０９において制御部２０１は、例えば、評価対象期間の各ウィンドウ期間ｏのそれぞれにおいて、複数のセンサの各センサｉごとに変化度を算出する。変化度は、例えば、上述の式３に示すＫＬダイバージェンスの式においてＸｊの代わりにＸｏを用いることで算出することができ、制御部２０１は、それによりウィンドウ期間ｏごとに変化スコアＡｏを取得する。
【0054】
Ｓ５１０において制御部２０１は、得られた変化スコアＡｏをクラスタに追加する。一例では、制御部２０１は、多次元空間上に変化スコアＡｏをプロットした点から最も近いクラスタに、変化スコアＡｏをクラスタリングしてよい。或いは、別の実施形態では、制御部２０１は、変化スコアＡｏを加えて再度、混合ガウスモデルなどによるクラスタリングを実行し、変化スコアＡｏがクラスタリングされるクラスタを決定してもよい。そして、例えば、変化スコアＡｏをクラスタに追加することにより、クラスタが移動し得る。図６（ｃ）には、変化スコアＡｏをクラスタに追加した例が示されており、クラスタへの変化スコアＡｏの追加により、クラスタが移動している。
【0055】
Ｓ５１１において制御部２０１は、クラスタの位置を再度特定する。例えば、制御部２０１は、クラスタの位置として、各クラスタを代表する代表点μｔ（ｔ＝１…ｓ，ｓはクラスタ数）を特定してよい。代表点には、例えば、クラスタの重心など、Ｓ５０５で用いたのと同様の点を用いることができる。なお、クラスタには評価対象期間の計測データから得られた変化スコアＡｏの点が含まれるため、クラスタの代表点の位置には、Ｓ５０５と比較して変化が生じ得る。
【0056】
Ｓ５１２において制御部２０１は、クラスタの位置の移動に基づいて、評価対象物の劣化を評価し、評価結果を出力し、本動作フローは終了する。評価結果は、例えば、評価対象物の劣化の度合いを表す劣化度であってよい。例えば、Ｓ５１２において制御部２０１は、Ｓ５０５で特定した学習期間におけるクラスタの位置から、Ｓ５１１で特定した評価対象期間の変化スコアの点を含むクラスタの位置までのクラスタの移動を表す値を、評価対象物の劣化度として出力してよい。移動を表す値は、例えば、クラスタの移動距離、移動速度、および移動加速度などであってよい。例えば、図６（ｄ）では、２つのクラスタの重心間の距離を劣化度として用いる例が示されており、ユーザは劣化度として出力された重心間の距離の大きさから、評価対象物の劣化の状態を推定することができる。
【0057】
なお、図６（ｄ）では、１つのクラスタの移動を表す値を劣化度として用いる例が示されているが、クラスタリングの結果得られた複数のクラスタが移動することも起こり得る。この場合、クラスタの移動を表す値は、例えば、複数のクラスタのそれぞれのクラスタの移動を表す値の平均値、中央値、最頻値、および分散などの統計値で表されてよい。図７は、複数のクラスタの移動を例示する図である。例えば、図７（ａ）に示すように、３つのクラスタＡ、クラスタＢ、およびクラスタＣがあり、それぞれのクラスタが移動しているとする。この場合に、制御部２０１は、３つのクラスタのそれぞれで取得した移動を表す値の統計値（例えば、平均値）を、劣化度として出力してよい。
【0058】
また、別の実施形態では、制御部２０１は、クラスタの移動を表す値が所定値以下で移動の小さいクラスタについては、劣化の評価において利用しなくてもよい。例えば、図７（ｂ）では、クラスタＢの移動は非常に小さい。このように、移動の小さいクラスタは、評価対象物の劣化に基づく変化の成分をあまり含んでいないことが推定される。そのため、制御部２０１は、例えば、移動を表す値が所定値以下のクラスタについては除いて、残りのクラスタの移動を表す値から、劣化度を取得してもよい。例えば、図７（ｂ）の例では、制御部２０１は、移動を表す値が所定値以下のクラスタＢを除いて、クラスタＡおよびクラスタＣの移動を表す値の統計値（例えば、平均値）を、劣化度として出力してよい。
【0059】
また、Ｓ５１２において制御部２０１が出力する評価対象物の劣化の評価結果の別な例として、制御部２０１は、クラスタの移動に基づいて、評価対象物に対する保守の要否を判定して出力してもよい。例えば、制御部２０１は、判定情報８００を参照し、評価対象物に対する保守の要否を判定してよい。
【0060】
図８は、実施形態に係る判定情報８００を例示する図である。判定情報８００には、例えば、評価対象物に対する保守作業の要否を判定するための情報が登録されていてよい。例えば、図８（ａ）の判定情報８００には閾値が登録されている。閾値には、例えば、評価対象物の劣化が少ない状態でのクラスタの位置から、評価対象物に保守作業を行うことが望ましい状態でのクラスタの位置までのクラスタの移動を表す値が予め調べられ、登録されている。それにより、例えば、制御部２０１は、Ｓ５１２において特定したクラスタの移動を表す値が閾値よりも大きい場合に、評価対象物に対する保守作業を促す情報を出力してもよい。保守作業を促す情報の出力は、例えば、情報処理装置２００と接続されたディスプレイなどの表示装置に、評価対象物の保守作業を実行するよう促すメッセージなどの情報を表示させることであってよい。また、別な例では、保守作業を促す情報の出力は、メールなどで評価対象物の管理者に保守作業を促す情報を通知することであってよい。
【0061】
また、判定情報８００には、例えば、評価対象物の劣化の種別ごとに保守作業の要否を判定するための情報が登録されてもよい。例えば、図８（ｂ）の判定情報８００では、評価対象物の劣化の種別と対応づけてクラスタの移動距離の閾値に加えて、移動方向の情報が登録されている。また、図８（ｃ）の判定情報８００では、評価対象物の劣化の種別と対応づけてクラスタの移動距離の閾値と、移動距離の特定に用いるクラスタを指定する情報とが対応づけて登録されている。図８（ｄ）の判定情報８００では、評価対象物の劣化の種別と対応づけてクラスタの移動距離の閾値と、移動速度の範囲と、移動の加速度の範囲とが登録されている。
【0062】
図１を参照して述べたように、評価対象物の劣化には様々な種類があることがあり、劣化の過程も劣化の種類に応じて異なることがある。そして、劣化の種類に応じてクラスタの移動の傾向も異なり得る。この場合にも、例えば、図８（ｂ）に示すクラスタの移動方向や、図８（ｃ）に示す移動距離の特定に用いるクラスタを指定する情報などのクラスタの移動を特徴付ける別の情報を用いることで、劣化の種別を特定することが可能である。また、同様に、例えば、図８（ｄ）に示すように、クラスタの移動を表す値としてクラスタの移動距離の閾値、移動速度の範囲、移動の加速度の範囲といった複数の情報を用いることで、劣化の種別を特定することが可能である。そして、Ｓ５１２の処理では、制御部２０１は、例えば、判定情報８００を参照し、劣化の種別に応じて評価対象物に対する保守の要否を判定し、判定結果を出力してもよい。例えば、制御部２０１は、クラスタの移動が、判定情報８００に登録されている長時間クリープの劣化の条件を満たす場合、長時間クリープの劣化の保守を促す情報を出力してよい。一方、制御部２０１は、クラスタの移動が、判定情報８００に登録されている腐食・摩耗の劣化の条件を満たす場合、腐食・摩耗の劣化の保守を促す情報を出力してよい。
【0063】
以上で述べたように、図５の動作フローによれば、制御部２０１は、変化スコアをクラスタリングして得られたクラスタの移動に基づいて、評価対象物の劣化を評価することができる。
【0064】
例えば、定期的に保守作業を実行して評価対象物の状態を維持する場合、たまたま劣化が早く進んでしまうことも起こり得るため、保守作業の実行間隔は、評価対象物の平均的な寿命よりも短めに設定され得る。しかしながら、実施形態によれば、クラスタの移動に基づいて、評価対象物の劣化状態を高い精度で推定することができるため、保守作業の実行間隔を長くしたりすることが可能である。その結果、保守作業にかかるコストや、作業負担を軽減することができる。
【0065】
また、たまたま評価対象物の劣化が早く進んでしまった場合にも、実施形態によればクラスタの移動に基づいて、評価対象物の劣化を高い精度で検出することができるため、ユーザは評価対象物が壊れる前に保守作業を実施することができる。そのため、評価対象物の計画外の停止などを回避することができる。なお、上記においては評価対象物の例としてボイラーチューブ群の劣化を述べたが、実施形態はこれに限定されるものではなく、その他の部品、機器、および設備の劣化にも実施形態は適用することができる。
【0066】
以上において、実施形態を例示したが、実施形態はこれに限定されるものではない。例えば、上述の動作フローは例示であり、実施形態はこれに限定されるものではない。可能な場合には、動作フローは、処理の順番を変更して実行されてもよく、別に更なる処理を含んでもよく、または、一部の処理が省略されてもよい。
【0067】
例えば、図５では、Ｓ５０１からＳ５０５までの処理と、Ｓ５０６からＳ５１２までの処理とを続けて実行する例が示されているが、実施形態はこれに限定されるものではない。別の実施形態では、Ｓ５０１からＳ５０５までの処理と、Ｓ５０６からＳ５１２までの処理とは別の動作フローとして実行されてもよい。一例として、Ｓ５０１からＳ５０５までの処理は、評価対象物の稼働初期の動作が安定している時期に１度だけ実行されればよい。また、例えば、Ｓ５０６からＳ５１２までの処理は、月に１度、半年に１度などの所定期間の計測データが蓄積すると、その都度、その所定期間の計測データを評価対象期間のデータとして実行されてよい。
【0068】
また、上述の実施形態では、例えば、クラスタの移動距離を、クラスタの重心間の距離を用いて特定する例を述べているが実施形態はこれに限定されるものではなく、その他の値をクラスタの移動距離として用いることもできる。例えば、別の実施形態では、移動前のクラスタ内の複数の変化スコアの点と、移動後のクラスタ内の複数の変化スコアの点との間で、最も短い距離または最も長い距離をクラスタの移動距離として用いてもよい。或いは、別の実施形態では、移動前と移動後の２つのクラスタに含まれる変化スコアの点の全ての組み合わせでの距離の平均を、クラスタの移動距離として用いてもよい。
【0069】
また、上述の実施形態において、評価対象期間のウィンドウ期間と対応する変化スコアの点のクラスタへの追加は、一例では、計測データの計測の時系列に従う順序で実行されてよい。クラスタへの変化スコアの点の追加を、時系列に従う順序で実行することで、クラスタの位置の移動の時系列な変化を捉えることができる。しかしながら、実施形態はこれに限定されるものではない。例えば、別の実施形態では、制御部２０１は、所定期間ごとにまとめて、所定期間に含まれるウィンドウ期間と対応する変化スコアの点を時系列に従ってクラスタに追加してもよい。この場合、所定期間内のウィンドウ期間と対応する変化スコアの点は任意の順序で追加されてもよい。
【0070】
なお、上述の実施形態においてＳ５０３の処理では制御部２０１は、例えば、取得部２１１として動作する。また、Ｓ５０５の処理では制御部２０１は、例えば、特定部２１２として動作する。Ｓ５１２の処理では制御部２０１は、例えば、評価部２１３として動作する。
【0071】
図９は、実施形態に係る情報処理装置２００を実現するためのコンピュータ９００のハードウェア構成を例示する図である。図９の情報処理装置２００を実現するためのハードウェア構成は、例えば、プロセッサ９０１、メモリ９０２、記憶装置９０３、読取装置９０４、通信インタフェース９０６、および入出力インタフェース９０７を備える。なお、プロセッサ９０１、メモリ９０２、記憶装置９０３、読取装置９０４、通信インタフェース９０６、入出力インタフェース９０７は、例えば、バス９０８を介して互いに接続されている。
【0072】
プロセッサ９０１は、例えば、シングルプロセッサであっても、マルチプロセッサやマルチコアであってもよい。プロセッサ９０１は、メモリ９０２を利用して例えば上述の動作フローの手順を記述したプログラムを実行することにより、上述した制御部２０１の一部または全部の機能を提供する。例えば、情報処理装置２００のプロセッサ９０１は、記憶装置９０３に格納されているプログラムを読み出して実行することで、取得部２１１、特定部２１２、および評価部２１３として動作する。
【0073】
メモリ９０２は、例えば半導体メモリであり、ＲＡＭ領域およびＲＯＭ領域を含んでよい。記憶装置９０３は、例えばハードディスク、フラッシュメモリ等の半導体メモリ、または外部記憶装置である。なお、ＲＡＭは、Random Access Memoryの略称である。また、ＲＯＭは、Read Only Memoryの略称である。
【0074】
読取装置９０４は、プロセッサ９０１の指示に従って着脱可能記憶媒体９０５にアクセスする。着脱可能記憶媒体９０５は、例えば、半導体デバイス、磁気的作用により情報が入出力される媒体、光学的作用により情報が入出力される媒体などにより実現される。なお、半導体デバイスは、例えば、ＵＳＢ（Universal Serial Bus）メモリである。また、磁気的作用により情報が入出力される媒体は、例えば、磁気ディスクである。光学的作用により情報が入出力される媒体は、例えば、ＣＤ－ＲＯＭ、ＤＶＤ、Blu-ray Disc等（Blu-rayは登録商標）である。ＣＤは、Compact Discの略称である。ＤＶＤは、Digital Versatile Diskの略称である。
【0075】
記憶部２０２は、例えばメモリ９０２、記憶装置９０３、および着脱可能記憶媒体９０５を含んでよい。例えば、情報処理装置２００の記憶装置９０３には、例えば、判定情報８００が格納されている。
【0076】
通信インタフェース９０６は、プロセッサ９０１の指示に従って、他の装置と通信する。一例では、通信インタフェース９０６は、評価対象物の状態を計測するセンサ９５０と有線または無線通信で通信し、センサ９５０から計測データを収集してよい。また、別の実施形態では、通信インタフェース９０６は、センサ９５０で計測された計測データを記憶する別の装置と有線または無線通信で通信し、計測データを取得してもよい。通信インタフェース９０６は、上述の通信部２０３の一例である。
【0077】
入出力インタフェース９０７は、例えば、入力装置および出力装置との間のインタフェースである。入力装置は、例えばユーザからの指示を受け付けるキーボード、マウス、タッチパネルなどのデバイスである。出力装置は、例えばディスプレイなどの表示装置、およびスピーカなどの音声装置である。
【0078】
実施形態に係る各プログラムは、例えば、下記の形態で情報処理装置２００に提供される。
（１）記憶装置９０３に予めインストールされている。
（２）着脱可能記憶媒体９０５により提供される。
（３）プログラムサーバなどのサーバから提供される。
（１）記憶装置９０３に予めインストールされている。
（２）着脱可能記憶媒体９０５により提供される。
（３）プログラムサーバなどのサーバから提供される。
【0079】
なお、図９を参照して述べた情報処理装置２００を実現するためのコンピュータ９００のハードウェア構成は、例示であり、実施形態はこれに限定されるものではない。例えば、上述の構成の一部が削除されてもよく、また、新たな構成が追加されてもよい。また、別の実施形態では、例えば、上述の制御部２０１の一部または全部の機能がＦＰＧＡ、ＳｏＣ、ＡＳＩＣ、およびＰＬＤなどによるハードウェアとして実装されてもよい。なお、ＦＰＧＡは、Field Programmable Gate Arrayの略称である。ＳｏＣは、System-on-a-chipの略称である。ＡＳＩＣは、Application Specific Integrated Circuitの略称である。ＰＬＤは、Programmable Logic Deviceの略称である。
【0080】
以上において、いくつかの実施形態が説明される。しかしながら、実施形態は上記の実施形態に限定されるものではなく、上述の実施形態の各種変形形態および代替形態を包含するものとして理解されるべきである。例えば、各種実施形態は、その趣旨および範囲を逸脱しない範囲で構成要素を変形して具体化できることが理解されよう。また、前述した実施形態に開示されている複数の構成要素を適宜組み合わせることにより、種々の実施形態が実施され得ることが理解されよう。更には、実施形態に示される全構成要素からいくつかの構成要素を削除して、または実施形態に示される構成要素にいくつかの構成要素を追加して種々の実施形態が実施され得ることが当業者には理解されよう。
【符号の説明】
【0081】
２００ 情報処理装置
２０１ 制御部
２０２ 記憶部
２０３ 通信部
２１１ 取得部
２１２ 特定部
２１３ 評価部
８００ 判定情報
９００ コンピュータ
９０１ プロセッサ
９０２ メモリ
９０３ 記憶装置
９０４ 読取装置
９０５ 着脱可能記憶媒体
９０６ 通信インタフェース
９０７ 入出力インタフェース
９０８ バス
９５０ センサ
２０１ 制御部
２０２ 記憶部
２０３ 通信部
２１１ 取得部
２１２ 特定部
２１３ 評価部
８００ 判定情報
９００ コンピュータ
９０１ プロセッサ
９０２ メモリ
９０３ 記憶装置
９０４ 読取装置
９０５ 着脱可能記憶媒体
９０６ 通信インタフェース
９０７ 入出力インタフェース
９０８ バス
９５０ センサ
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】