(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024111329
(43)【公開日】2024-08-16
(54)【発明の名称】要因分析装置及び要因分析方法
(51)【国際特許分類】
G05B 23/02 20060101AFI20240808BHJP
【FI】
G05B23/02 T
【審査請求】有
【請求項の数】5
【出願形態】OL
(21)【出願番号】P 2024102724
(22)【出願日】2024-06-26
(62)【分割の表示】P 2024023566の分割
【原出願日】2020-05-18
(31)【優先権主張番号】P 2019185739
(32)【優先日】2019-10-09
(33)【優先権主張国・地域又は機関】JP
(71)【出願人】
【識別番号】000005234
【氏名又は名称】富士電機株式会社
(74)【代理人】
【識別番号】100107766
【弁理士】
【氏名又は名称】伊東 忠重
(74)【代理人】
【識別番号】100070150
【弁理士】
【氏名又は名称】伊東 忠彦
(72)【発明者】
【氏名】石橋 直人
(72)【発明者】
【氏名】村上 賢哉
(57)【要約】
【課題】性能変化の要因を適切に抽出すること。
【解決手段】一実施形態に係る要因分析装置は、対象を複数のセンサにより計測した計測データを収集する収集手段と、前記収集手段により収集した計測データを用いて、前記対象の性能値を推定する性能推定手段と、前記性能推定手段により推定された性能値の分析期間を設定し、前記分析期間において性能値が変化した計測データを検知する変化検知手段と、前記変化検知手段により設定された前記分析期間において、性能変化の要因を前記計測データから分析する分析手段と、を有することを特徴とする。
【選択図】
図1
【特許請求の範囲】
【請求項1】
対象を複数のセンサにより計測した計測データを収集する収集手段と、
前記収集手段により収集した計測データを用いて、前記対象の性能値を推定する性能推定手段と、
前記性能推定手段により推定された性能値の分析期間を設定し、前記分析期間において性能値が変化した計測データを検知する変化検知手段と、
前記変化検知手段により設定された前記分析期間において、性能変化の要因を前記計測データから分析する分析手段と、
を有することを特徴とする要因分析装置。
【請求項2】
前記分析手段は、
前記分析期間における前記計測データから分類二進木を構築し、該分類二進木を分析することで、前記性能変化の要因を抽出する、ことを特徴とする請求項1に記載の要因分析装置。
【請求項3】
前記変化検知手段は、
前記分析期間において、所定のクラスタリング手法により前記計測データから前記性能値を複数のクラスタに分割し、
各クラスタにそれぞれ属する計測データに対してラベルを付与することで、前記性能値が変化した計測データを検知する、ことを特徴とする請求項1又は2に記載の要因分析装置。
【請求項4】
前記変化検知手段は、
前記分析期間において、所定の異常検知手法により前記計測データから前記性能値を複数の正常又は異常の種別に分割し、
各種別にそれぞれ分割された計測データに対してラベルを付与することで、前記性能値が変化した計測データを検知する、ことを特徴とする請求項1又は2に記載の要因分析装置。
【請求項5】
対象を複数のセンサにより計測した計測データを収集する収集手順と、
前記収集手順で収集した計測データを用いて、前記対象の性能値を推定する性能推定手順と、
前記性能推定手順で推定された性能値の分析期間を設定し、前記分析期間において性能値が変化した計測データを検知する変化検知手順と、
前記変化検知手順で設定された前記分析期間において、性能変化の要因を前記計測データから分析する分析手順と、
をコンピュータが実行することを特徴とする要因分析方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、要因分析装置及び要因分析方法に関する。
【背景技術】
【0002】
近年、センサの発達や通信費の低下等に伴ってセンサデータの収集が容易となってきており、これらのセンサデータを活用して様々な対象(例えば、プラント、各種設備、機器、装置等)の効率的な運転を実現するサービスが期待されている。このような対象の効率的な運転を実現するために、当該対象の性能の変化要因を分析する技術が提案されている。
【0003】
例えば、特許文献1には、外気温や海水温度等の外的要因を除去し、実測値に基づいて発電プラントの熱効率、発電プラントを構成する機器の性能値、及び当該性能値の変化が発電プラントの熱効率に与える影響を解析することができる技術が開示されている。また、例えば、特許文献2には、プラント固有の知見や経験によるノウハウで補間することなく、精度の高い効率診断とその対応を支援するための技術が開示されている。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特許第5618322号公報
【特許文献2】特許第6472339号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかしながら、例えば、特許文献1に開示されている技術では発電プラントの熱効率に影響を与える機器を特定することは可能であるが、当該機器をセンシングするどのセンサが熱効率に影響を与えるのかまでは特定することができなかった。また、例えば、特許文献2に開示されている技術では分析木により性能低下の要因を分析しているが、性能低下の判断と分析木による要因抽出とが異なるロジックであるため、性能の低下が判断された際に分析木にてその要因を必ずしも抽出できない可能性がある。
【0006】
本発明の一実施形態は、上記の点に鑑みてなされたもので、性能変化の要因を適切に抽出することを目的とする。
【課題を解決するための手段】
【0007】
上記目的を達成するため、一実施形態に係る要因分析装置は、対象を複数のセンサにより計測した計測データを収集する収集手段と、前記収集手段により収集した計測データを用いて、前記対象の性能値を推定する性能推定手段と、前記性能推定手段により推定された性能値の分析期間を設定し、前記分析期間において性能値が変化した計測データを検知する変化検知手段と、前記変化検知手段により設定された前記分析期間において、性能変化の要因を前記計測データから分析する分析手段と、を有することを特徴とする。
【発明の効果】
【0008】
性能変化の要因を適切に抽出することができる。
【図面の簡単な説明】
【0009】
【
図1】第一の実施形態に係る要因分析装置の全体構成の一例を示す図である。
【
図2】第一の実施形態に係る要因分析装置のハードウェア構成の一例を示す図である。
【
図3】第一の実施形態に係る要因分析処理の流れの一例を示すフローチャートである。
【
図4】性能変化期間の検知結果の一例を示す図である。
【
図5】第二の実施形態に係る要因分析装置の全体構成の一例を示す図である。
【
図6】第二の実施形態に係る要因分析処理の流れの一例を示すフローチャートである。
【
図7】第三の実施形態に係る要因分析装置の全体構成の一例を示す図である。
【
図8】第三の実施形態に係る要因分析処理の流れの一例を示すフローチャートである。
【発明を実施するための形態】
【0010】
以下、本発明の各実施形態について説明する。本発明の各実施形態では、性能変化の要因を適切に抽出することが可能な要因分析装置10について説明する。
【0011】
[第一の実施形態]
以降では、本発明の第一の実施形態について説明する。第一の実施形態では、性能を効率化する対象(例えば、プラント、各種設備、機器、装置等)をセンシングすることで得られた計測データを用いて、性能の変化を検知すると共にその変化の要因を抽出するための決定木として分類二進木を構築し、この分類二進木を分析することによって性能変化の要因を抽出する場合について説明する。
【0012】
なお、以降では、性能を効率化する対象を「効率化対象」とも表す。効率化対象には、上述したように、例えば、プラント、各種設備、機器、装置等が挙げられる。また、各効率化対象は、1以上のセンサによってその状態や動作等がセンシングされているものとする。
【0013】
<全体構成>
まず、本実施形態に係る要因分析装置10の全体構成について、
図1を参照しながら説明する。
図1は、第一の実施形態に係る要因分析装置10の全体構成の一例を示す図である。
【0014】
図1に示すように、本実施形態に係る要因分析装置10は、収集部101と、性能推定部102と、変化検知部103と、要因分析部104と、記憶部105とを有する。
【0015】
収集部101は、効率化対象をセンシングすることで得られた計測データを収集し、記憶部105に格納する。ここで、計測データには、例えば、効率化対象を識別する識別情報(例えば、プラントIDや設備ID、機器ID、装置ID等)と、当該効率化対象の状態や動作等がセンシング(計測)された日時と、これらの計測を行ったセンサの識別情報(例えば、センサID等)と、当該センサの計測値(センサ値)とが含まれる。なお、例えば、効率化対象の状態や動作をN個のセンサで計測している場合、計測データには、これらN個のセンサのセンサIDと、これらN個のセンサのセンサ値とが含まれる。
【0016】
なお、収集部101は、任意の方法で計測データを収集すればよい。例えば、収集部101は、LAN(Local Area Network)等の通信ネットワークを介して計測機器やセンサ等から計測データを受信することで収集してもよいし、キーボード等の入力装置によって計測データが入力されることで収集してもよい。
【0017】
また、収集部101は、計測データ以外にも、例えば、効率化対象の運転データを収集してもよい。運転データには、例えば、効率化対象の計測値の正常範囲や異常範囲、効率化対象の各種イベント(例えば、操業日や設備点検日等)を示す情報等が含まれる。運転データが収集された場合には、これらの運転データも記憶部105に格納される。
【0018】
性能推定部102は、記憶部105に格納されている計測データを用いて、効率化対象の性能値を推定する。これにより、各計測データのそれぞれに対して性能値が推定される。
【0019】
性能値は、効率化対象の入出力関係から、物理的に定義された定量的な性能換算式によって推定される。例えば、性能値が効率である場合は、一般に、プラントや設備毎に出力値を入力値で割った換算式で推定される。
【0020】
具体的には、例えば、効率化対象がガスタービン発電所のプラントである場合、その性能値は、以下の式(1)により推定することができる。
【0021】
【数1】
ここで、Performanceは性能値、GeneratorOutputは発電出力、GasFlowは燃料ガス流量を表す。発電出力及び燃料ガス流量は、計測データに含まれる計測値である。
【0022】
ただし、上記の式(1)は一例であって、性能値の換算式は、効率化対象の種類等に応じて異なる。例えば、効率化対象がタービン室であれば、その性能値(タービン室効率)は、タービン室出力をタービン室入力で割ることで推定することができる。タービン室出力及びタービン室入力は、計測データに含まれる計測値である。なお、本実施形態は特定の種類の性能値に限定されるものではなく、任意の種類の性能値に対して適用可能である。
【0023】
変化検知部103は、効率化対象毎に、性能推定部102によって推定された性能値を複数の期間に分割した上で、各期間で所定の統計量を算出し、この統計量によって各期間をラベル付けすることで、性能が変化(例えば、性能が低下)した前後の期間を性能変化期間として検知する。
【0024】
例えば、期間P1~期間Pnのn個の期間に性能値が分割され、性能値の高いことを示す第1のラベル又は性能値が低いことを示す第2のラベルのいずれかが各期間にラベル付けされる場合、変化検知部103は、第1のラベルが付与されている期間Pmと第2のラベルが付与されている期間Pm+1とが存在するときに、これらの期間Pm及びPm+1を性能変化期間として検知する。
【0025】
要因分析部104は、変化検知部103によって検知された性能変化期間において、計測データから性能変化の要因を抽出する。このとき、要因分析部104は、当該性能変化期間における計測データから分類二進木を構築し、この分類二進木を分析することによって性能変化の要因を抽出する。
【0026】
記憶部105は、収集部101によって収集された計測データを記憶する。記憶部105には、例えば、効率化対象毎に、計測データが時系列データとして記憶されている。なお、記憶部105には、例えば、効率化対象毎に、運転データが記憶されていてもよい。また、これら以外にも、記憶部105には、種々のデータ(例えば、性能変化の要因を抽出するための決定木、この決定木を分析することによって性能変化の要因を抽出する際における途中の計算結果等)が記憶されてもよい。
【0027】
なお、
図1に示す要因分析装置10の構成は一例であって、他の構成であってもよい。例えば、要因分析装置10が複数の装置で構成されており、これら複数の装置に各部(収集部101、性能推定部102、変化検知部103、要因分析部104及び記憶部105)が分散されていてもよい。
【0028】
<ハードウェア構成>
次に、第一の実施形態に係る要因分析装置10のハードウェア構成について、
図2を参照しながら説明する。
図2は、第一の実施形態に係る要因分析装置10のハードウェア構成の一例を示す図である。
【0029】
図2に示すように、第一の実施形態に係る要因分析装置10は、入力装置201と、表示装置202と、外部I/F203と、通信I/F204と、プロセッサ205と、メモリ装置206とを有する。これら各ハードウェアは、バス207により相互に通信可能に接続されている。
【0030】
入力装置201は、例えば各種ボタンやタッチパネル、キーボード、マウス等であり、ユーザが要因分析装置10に各種操作を入力するのに用いられる。表示装置202は、例えばディスプレイ等であり、要因分析装置10の処理結果等を表示する。なお、要因分析装置10は、入力装置201及び表示装置202のうちの少なくとも一方を有していなくてもよい。
【0031】
外部I/F203は、外部装置とのインタフェースである。外部装置には、記録媒体203a等がある。要因分析装置10は、外部I/F203を介して、記録媒体203aの読み取りや書き込み等を行うことができる。要因分析装置10が有する各機能部(例えば、収集部101、性能推定部102、変化検知部103及び要因分析部104等)を実現する1以上のプログラムは、記録媒体203aに記録されていてもよい。
【0032】
記録媒体203aには、例えば、SDメモリカード(SD memory card)やUSBメモリ、CD(Compact Disk)、DVD(Digital Versatile Disk)等がある。
【0033】
通信I/F204は、要因分析装置10が他の装置や機器等とデータ通信を行うためのインタフェースである。要因分析装置10が有する各機能部を実現する1以上のプログラムは、通信I/F204を介して、所定のサーバ装置等から取得(ダウンロード)されてもよい。
【0034】
プロセッサ205は、例えばCPU(Central Processing Unit)等であり、メモリ装置206等からプログラムやデータを読み出して各種処理を実行する演算装置である。要因分析装置10が有する各機能部は、メモリ装置206等に格納されている1以上のプログラムがプロセッサ205に実行させる処理により実現される。
【0035】
メモリ装置206は、例えばHDD(Hard Disk Drive)やSSD(Solid State Drive)、RAM(Random Access Memory)、ROM(Read Only Memory)等であり、プログラムやデータを格納している記憶装置である。要因分析装置10が有する記憶部105は、例えば、メモリ装置206等を用いて実現可能である。ただし、当該記憶部105は、例えば、要因分析装置10と通信ネットワークを介して接続される記憶装置等を用いて実現されていてもよい。
【0036】
第一の実施形態に係る要因分析装置10は、
図2に示すハードウェア構成を有することにより、後述する要因分析処理を実現することができる。なお、
図2に示すハードウェア構成は一例であって、要因分析装置10は、他のハードウェア構成を有していてもよい。例えば、要因分析装置10は、複数のプロセッサ205を有していてもよいし、複数のメモリ装置206を有していてもよい。
【0037】
<要因分析処理の流れ>
次に、本実施形態に係る要因分析装置10によって性能変化の要因を抽出する要因分析処理の流れについて、
図3を参照しながら説明する。
図3は、第一の実施形態に係る要因分析処理の流れの一例を示すフローチャートである。なお、記憶部105には、収集部101によって収集された計測データが効率化対象毎に記憶されているものとする。
【0038】
以降では、或る1つの効率化対象の性能変化の要因を抽出する場合について説明する。分析対象とする効率化対象は、例えば、ユーザ等によって指定されてもよいし、全ての効率化対象であってもよい。
【0039】
まず、性能推定部102は、記憶部105に格納されている計測データを用いて、効率化対象の性能値を推定する(ステップS101)。これにより、各計測データのそれぞれに対して性能値が推定される。
【0040】
なお、1つの計測データに対して複数種類の性能値が推定されてもよいが、以降では、1つの計測データに対して1種類の性能値が推定されたものとして説明する。1つの計測データに対して複数種類の性能値が推定された場合には、以降のステップS102~ステップS103を性能値の種類毎に繰り返し実行すればよい。
【0041】
次に、変化検知部103は、以下のStep11~Step15により性能変化期間を検知する(ステップS102)。
【0042】
Step11:変化検知部103は、上記のステップS101で推定された性能値を複数の期間に分割する。ここで、期間の単位としては任意の単位を設定することができるが、例えば、年、季節、月、週、又は日を単位とすることが挙げられる。又は、例えば、記憶部105に運転データが格納されている場合は、この運転データを参照して、設備点検日を境として性能値を複数の期間に分割してもよい。
【0043】
Step12:次に、変化検知部103は、上記のStep11で分割された各期間で、性能値に影響を及ぼす外的要因を除去する補正を行う。変化検知部103は、例えば、上記の特許文献1に開示されている方法と同様に、外的要因と性能値との回帰式の傾きを用いて性能値を補正すればよい。なお、このStep12の処理の実行は任意であり、必ずしも実行されなくてもよい。以降では、上記の補正が行われた性能値を「補正後性能値」とも表す。
【0044】
Step13:次に、変化検知部103は、上記のStpe11で分割された各期間で、補正後性能値(又は、上記のStep12の処理が実行されなかった場合は性能値)の統計量を算出する。統計量としては任意の統計量を算出すればよいが、例えば、当該期間における補正後性能値(又は、性能値)の平均、標準偏差、最大、最小等が挙げられる。又は、例えば、異なる期間の間(例えば、期間Pmと期間Pm+1との間)における補正後性能値(又は、性能値)の相関係数であってもよい。
【0045】
Step14:次に、変化検知部103は、各期間で算出された統計量を用いて、これらの各期間に対して、当該期間における性能を表すラベルを付与する。このとき、変化検知部103は、例えば、統計検知や閾値判定等のアルゴリズムによりラベル付けを行えばよい。これにより、各期間に対して、例えば、当該期間における性能が高いか又は低いかを表すラベルが付与される。ただし、ラベルが2値であることは一例であって、各期間を任意の個数のクラスに分類するためラベルが付与されてもよい。各期間に対してラベルが付与されることにより、当該期間における計測データ及び補正後性能値(又は性能値)に対しても同様のラベルが付与されることになる。
【0046】
なお、このStep14では、人手によりラベル付けが行われてもよい。この場合、上記のStep13で算出された期間毎の統計量や期間毎の補正後性能値(又は性能値)をユーザに提示することで、当該ユーザによってラベル付けが行われる。
【0047】
以降では、簡単のため、補正後性能値と、補正を行わなかった場合の性能値とを区別せずに、単に「性能値」と表す。
【0048】
Step15:そして、変化検知部103は、上記のStep14で各期間に付与されたラベルを用いて、性能が変化(例えば、性能が低下)した前後の期間を性能変化期間として検知する。
【0049】
上記のStep11~Step15で検知された性能変化期間の一例を
図4に示す。
図4に示す例では、上記のStep11で年を単位とした期間が設定された場合を示している。
図4に示す例では、期間P
mの統計量がS
m、期間P
m+1の統計量がS
m+1で、これらの統計量S
m及びS
m+1によって期間P
mに対してラベル「高」が、期間P
m+1に対してラベル「低」が付与されている。このため、期間P
mと期間P
m+1とで性能が変化しているため、期間P
m及びP
m+1が性能変化期間として検知される。
【0050】
なお、上記のステップS102で性能変化期間が検知されなかった場合、以降のステップS103の処理は実行されない。
【0051】
最後に、要因分析部104は、以下のStep21~Step23により上記のステップS102で検知された性能変化期間において性能変化の要因を抽出する(ステップS103)。
【0052】
Step21:まず、要因分析部104は、分析対象とする因子を任意に設定する。例えば、要因分析部104は、分析対象とする因子として、任意のセンサIDを設定する。なお、要因分析部104は、任意の個数のセンサIDを因子として設定する。このとき、どのセンサIDを分析対象の因子として設定するかは任意に決めることが可能であるが、例えば、ユーザ等によって指定されたセンサIDを分析対象の因子として設定すればよい。
【0053】
Step22:次に、要因分析部104は、上記のステップS102で検知された性能変化期間でラベル付けされた性能値と、上記のStep21で設定した因子とを用いて、各因子を分岐ノード、ラベルの値をリーフノードとして、計測データの性能値を分類する分類二進木を構築する。すなわち、要因分析部104は、各因子を分岐ノード(ルートノードも含む)、ラベルの値をリーフノードとして、計測データをルートノードに入力した場合に、当該計測データの性能値を正しいラベルの値(つまり、上記のステップS102で当該性能値に対して付与されたラベルの値)に分類する分類二進木を構築する。
【0054】
Step23:そして、要因分析部104は、上記のStep22で構築した分類二進木を用いて、性能変化の要因を抽出する。ここで、要因分析部104は、分岐ノードを説明変数、リーフノードを目的変数として以下の式(2)により変数重要度VI(x)を算出し、この変数重要度VI(x)が最も高い(又は、所定の閾値以上)の説明変数に対応する要因を性能変化の要因として抽出すればよい。なお、変数重要度VI(x)をユーザに提示し、このユーザの選択によって性能変化の要因が抽出されてもよい。
【0055】
【数2】
ここで、xは説明変数、tは分岐ノード番号、Xは説明変数の集合、N
sは分岐ノードの集合、Δi(x,t)は説明変数xのΔi(t)である。また、Δi(t)は分岐ノードtの改善度である。
【0056】
以上のように、本実施形態に係る要因分析装置10では、計測データを複数の期間に分割した上で、これらの期間を性能に応じてラベル付けすることで性能変化期間を検知する。そして、本実施形態に係る要因分析装置10では、ラベル付けされた性能変化期間から分類二進木を構築し、その分析を行うことで、性能変化の要因(例えば、センサID)を抽出する。このため、本実施形態に係る要因分析装置10では、例えば、或る効率化対象を計測する複数のセンサの中で、性能変化の要因となっているセンサを特定することが可能となる。
【0057】
[第二の実施形態]
以降では、本発明の第二の実施形態について説明する。第二の実施形態では、計測データを用いて、性能変化の要因を分析するための期間を設定すると共にその期間での性能変化の要因を抽出するための決定木として回帰二進木を構築し、この回帰二進木を分析することによって性能変化の要因を抽出する場合について説明する。
【0058】
なお、第二の実施形態では、主に、第一の実施形態との相違点について説明し、第一の実施形態と同様の構成要素についてはその説明を省略するものとする。
【0059】
<全体構成>
まず、本実施形態に係る要因分析装置10の全体構成について、
図5を参照しながら説明する。
図5は、第二の実施形態に係る要因分析装置10の全体構成の一例を示す図である。
【0060】
図5に示すように、本実施形態に係る要因分析装置10は、第一の実施形態と異なり、変化検知部103を有しない。また、本実施形態に係る要因分析装置10は、分析期間設定部106を有する。
【0061】
分析期間設定部106は、性能変化の要因を分析する期間(以降、「分析期間」とも表す。)を設定する。すなわち、第二の実施形態では、各期間に対してラベル付けを行って性能変化期間を検知することは行わずに、性能変化の要因を分析する分析期間のみを設定する。
【0062】
また、本実施形態に係る要因分析部104は、分析期間設定部106によって設定された分析期間において、計測データから性能変化の要因を抽出する。このとき、要因分析部104は、当該分析期間における計測データから回帰二進木を構築し、この回帰二進木を分析することによって性能変化の要因を抽出する。
【0063】
<要因分析処理の流れ>
次に、本実施形態に係る要因分析装置10によって性能変化の要因を抽出する要因分析処理の流れについて、
図6を参照しながら説明する。
図6は、第二の実施形態に係る要因分析処理の流れの一例を示すフローチャートである。なお、記憶部105には、収集部101によって収集された計測データが効率化対象毎に記憶されているものとする。また、以降では、或る1つの効率化対象の性能変化の要因を抽出する場合について説明する。
【0064】
まず、性能推定部102は、
図3のステップS101と同様に、記憶部105に格納されている計測データを用いて、効率化対象の性能値を推定する(ステップS201)。
【0065】
次に、分析期間設定部106は、上記のステップS201で推定された性能値に対して分析期間を設定する(ステップS202)。なお、このとき、当該分析期間における性能値に対して、上記のStep12と同様の補正を行ってもよい。
【0066】
なお、本実施形態では、上記のステップS201で性能値を推定した後に、上記のステップS202で分析期間を設定したが、これらの順序は逆であってもよい。すなわち、例えば、計測データに対して分析期間を設定した後に、当該分析期間における各計測データからそれぞれ性能値を推定してもよい。また、以降では、簡単のため、補正を行った場合の補正後性能値と、補正を行わなかった場合の性能値とを区別せずに、単に「性能値」と表す。
【0067】
最後に、要因分析部104は、以下のStep31~Step33により上記のステップS202で設定された分析期間において性能変化の要因を抽出する(ステップS203)。
【0068】
Step31:まず、要因分析部104は、上記のStep21と同様に、分析対象とする因子(例えば、センサID)を任意に設定する。
【0069】
Step32:次に、要因分析部104は、上記のステップS202で設定された分析期間における性能値と、上記のStep31で設定した因子とを用いて、各因子を分岐ノード、性能値をリーフノードとして、計測データと性能値との関係を表す回帰二進木を構築する。
【0070】
Step33:そして、要因分析部104は、上記のStep32で構築した回帰二進木を用いて、性能変化の要因を抽出する。なお、要因分析部104は、上記のStep23と同様に、分岐ノードを説明変数、リーフノードを目的変数として上記の式(2)により変数重要度VI(x)を算出し、この変数重要度VI(x)が最も高い(又は、所定の閾値以上)の説明変数に対応する要因を性能変化の要因として抽出すればよい。
【0071】
以上のように、本実施形態に係る要因分析装置10では、計測データに対して分析期間を設定した上で、この分析期間から回帰二進木を構築し、その分析を行うことで、性能変化の要因(例えば、センサID)を抽出する。このため、本実施形態に係る要因分析装置10では、第一の実施形態と同様に、例えば、或る効率化対象を計測する複数のセンサの中で、性能変化の要因となっているセンサを特定することが可能となる。
【0072】
[第三の実施形態]
以降では、本発明の第三の実施形態について説明する。第三の実施形態では、計測データを用いて、分析期間において性能変化のサンプルを検知すると共にその分析期間での性能変化の要因を抽出するための決定木として分類二進木を構築し、この分類二進木を分析することによって性能変化の要因を抽出する場合について説明する。
【0073】
なお、第三の実施形態では、主に、第一の実施形態との相違点について説明し、第一の実施形態と同様の構成要素についてはその説明を省略するものとする。
【0074】
<全体構成>
まず、本実施形態に係る要因分析装置10の全体構成について、
図7を参照しながら説明する。
図7は、第三の実施形態に係る要因分析装置10の全体構成の一例を示す図である。
【0075】
図7に示すように、本実施形態に係る要因分析装置10は、第一の実施形態と異なり、変化検知部103を有しない。また、本実施形態に係る要因分析装置10は、変化サンプル検知部107を有する。
【0076】
変化サンプル検知部107は、分析期間を設定すると共に性能値が変化したサンプル(以降、「性能変化サンプル」とも表す。)を検知する。すなわち、第三の実施形態では、各期間に対してラベル付けを行って性能変化期間を検知することは行わずに、各サンプルに対してラベル付けを行うことで性能変化が変化したサンプルを検知する。なお、サンプルとは、記憶部105に記憶されている計測データのことである。
【0077】
また、本実施形態に係る要因分析部104は、変化サンプル検知部107によって設定された分析期間において、サンプル(計測データ)から性能変化の要因を抽出する。このとき、要因分析部104は、第一の実施形態と同様に、当該分析期間における計測データから分類二進木を構築し、この分類二進木を分析することによって性能変化の要因を抽出する。
【0078】
<要因分析処理の流れ>
次に、本実施形態に係る要因分析装置10によって性能変化の要因を抽出する要因分析処理の流れについて、
図8を参照しながら説明する。
図8は、第三の実施形態に係る要因分析処理の流れの一例を示すフローチャートである。なお、記憶部105には、収集部101によって収集された計測データが効率化対象毎に記憶されているものとする。また、以降では、或る1つの効率化対象の性能変化の要因を抽出する場合について説明する。
【0079】
まず、性能推定部102は、
図3のステップS101と同様に、記憶部105に格納されている計測データを用いて、効率化対象の性能値を推定する(ステップS301)。
【0080】
次に、変化サンプル検知部107は、以下のStep41~Step45により性能変化サンプルを検知する(ステップS302)。
【0081】
Step41:変化サンプル検知部107は、上記のステップS301で推定された性能値に対して分析期間を設定する。
【0082】
Step42:次に、変化サンプル検知部107は、上記のStep12と同様に、上記のStep41で設定された分析期間で、性能値に影響を及ぼす外的要因を除去する補正を行う。なお、このStep42の処理の実行は任意であり、必ずしも実行されなくてもよい。以降では、簡単のため、補正を行った場合の性能値と補正を行わなかった場合の性能値とを区別せずに、単に「性能値」と表す。
【0083】
Step43:次に、変化サンプル検知部107は、上記のStep41で設定された分析期間において、計測データから性能値を複数のクラスタ又は種別に分割する。複数のクラスタに分割する際は、例えば、k-means法やGaussian Mixture Model等のクラスタリング手法を用いて、性能値の特徴に基づきサンプル毎に複数のクラスタに分割すればよい。一方で、複数の種別に分割する際は、例えば、Isolation ForestやRandom Cut Forest等の異常検知手法を用いて、性能値の異常度合いに基づきサンプル毎に複数の正常、異常の種別に分割すればよい。
【0084】
Step44:次に、変化サンプル検知部107は、上記のステップStep43で分割されたクラスタ又は種別に対して性能を表すラベルを付与する。これにより、各クラスタ又は各種別に属するサンプル(計測データ)に対してラベルが付与される。なお、ラベルとしては、例えば、正常又は異常のいずれであるかを表す2個のラベルを付与すればよい。ただし、ラベルが2個であることは一例であって、任意の個数のラベルが付与されてもよい。
【0085】
Step45:そして、変化サンプル検知部107は、上記のStep44で各サンプルに付与されたラベルを用いて、性能が変化したサンプル(例えば、ラベルの値が正常を示す値から異常を示す値に変化したサンプル)を性能変化サンプルとして検知する。
【0086】
最後に、要因分析部104は、
図3のステップS103と同様に、上記のStep41で設定された分析期間において性能変化の要因を抽出する(ステップS303)。すなわち、要因分析部104は、分析対象とする因子を任意に設定した上で、当該分析期間でラベル付けされた性能値と当該因子とを用いて分類二進木を構築し、この分類二進木を用いて、性能変化の要因を抽出する。
【0087】
以上のように、本実施形態に係る要因分析装置10では、計測データに対して分析期間を設定した上で、この分析期間から分類二進木を構築し、その分析を行うことで、性能変化の要因(例えば、センサID)を抽出する。このため、本実施形態に係る要因分析装置10では、第一の実施形態と同様に、例えば、或る効率化対象を計測する複数のセンサの中で、性能変化の要因となっているセンサを特定することが可能となる。
【0088】
本発明は、具体的に開示された上記の各実施形態に限定されるものではなく、特許請求の範囲から逸脱することなく、種々の変形や変更、組み合わせ等が可能である。
【符号の説明】
【0089】
10 要因分析装置
101 収集部
102 性能推定部
103 変化検知部
104 要因分析部
105 記憶部
106 分析期間設定部
107 変化サンプル検知部