(57)【要約】  状態予測装置１００は、観測センサ３００で得た計測データを、所定の時間間隔で取得し、観測対象の動作モード毎に区切り、学習データを生成する学習データ生成部２１２と、所定数の前記学習データを予め定めたカテゴリに分類してラベルを付し、教師有データを生成するクラスタリング部２１３と、教師有データを、学習アルゴリズムに入力して学習し、学習済みモデル２５２を生成するとともに、各カテゴリの代表データを生成する学習部２１４と、予め定めたモデル生成単位期間内に学習データに付されたラベル２４４を、当該学習データの生成順に時系列に並べ、時系列モデル２７３を生成する時系列モデル生成部２１５と、新たに取得した計測データから、時系列モデル２７３を用いて予測値を算出する予測部２２０と、を備える。
2019163141.tif 000002
【特許請求の範囲】
【請求項１】
観測対象の状態を観測する第一の観測センサで得た第一のセンサデータを、所定の時間間隔で取得し、前記観測対象の動作モード毎に区切り、学習データを生成する学習データ生成部と、
生成した所定数の前記学習データそれぞれを予め定めたカテゴリに分類し、前記カテゴリ毎に予め用意されたラベルを付し、教師有データを生成するクラスタリング部と、
前記所定数の教師有データを、学習アルゴリズムに入力して学習し、当該学習アルゴリズムの学習済みモデルを生成するとともに、各前記カテゴリの代表データを生成する学習部と、
予め定めたモデル生成単位期間内に前記クラスタリング部で前記学習データに付された前記ラベルを、当該学習データの生成順に時系列に並べ、時系列モデルを生成する時系列モデル生成部と、
新たに取得した第一のセンサデータから、前記時系列モデルを用いて未来時刻の前記観測対象の状態として前記第一のセンサデータの予測値を算出する予測部と、を備える状態予測装置。
【請求項２】
請求項１の状態予測装置であって、
前記第一のセンサデータを用いて前記動作モード毎の時間幅を切出時間幅として算出する時間幅算出部をさらに備え、
前記学習データ生成部は、前記時間幅算出部が算出した切出時間幅で前記第一のセンサデータを区切り、
前記時間幅算出部は、前記第一のセンサデータを取得する毎に当該第一のセンサデータの値のヒストグラムを生成し、当該ヒストグラムが山状の分布を示す時間幅を、前記切出時間幅とする状態予測装置。
【請求項３】
請求項２記載の状態予測装置であって、
前記時間幅算出部は、前記ヒストグラムを生成する毎に、当該ヒストグラムの分散の逆数を算出し、前回の切出時間直後から当該分散の逆数が極大値を取る時刻までの期間を前記切出時間幅と算出する状態予測装置。
【請求項４】
請求項１記載の状態予測装置であって、
各カテゴリの前記代表データと、当該カテゴリ毎の前記ラベルとを対応づけて記憶する代表データ記憶部と、
前記モデル生成単位期間毎に生成された前記時系列モデルであって、前記ラベルの並び順の異なる前記時系列モデルを、それぞれ時系列モデル識別情報に対応づけて記憶する時系列モデル記憶部と、をさらに備え、
前記予測部は、
所定の時間間隔で取得した前記第一のセンサデータを、前記動作モード毎に区切り、予測用データを生成する予測用データ生成部と、
生成した前記予測用データを、前記学習済みモデルに入力して前記カテゴリのいずれかに分類し、当該カテゴリに設定される前記ラベルを付与する分類決定部と、
前記ラベルが付与される毎に、当該ラベルを時系列に並べ、前記時系列モデル記憶部に記憶される前記時系列モデルの中で、最も並び順の類似度の高い時系列モデルを特定し、特定された当該時系列モデルに従って、未来時刻のラベルを予測ラベルとして特定する予測ラベル特定部と、
特定された前記予測ラベルに対応づけて前記代表データ記憶部に記憶される前記代表データを前記第一のセンサデータの前記予測値として算出する予測値算出部と、を備える状態予測装置。
【請求項５】
請求項２記載の状態予測装置であって、
前記学習データ生成部は、前記観測対象の、前記第一の観測センサとは異なる状態を観測する第二の観測センサで得た第二のセンサデータを、前記第一のセンサデータを用いて算出した前記切出時間幅で区切り、前記学習データをさらに生成し、
前記クラスタリング部は、前記第二のセンサデータから生成した前記学習データを、当該第二のセンサデータのために予め定められたカテゴリに分類し、前記カテゴリ毎に予め用意され第二のラベルを付し、当該第二のセンサデータの前記教師有データを生成し、
前記学習部は、所定数の前記第一のセンサデータの教師有データと前記第二のセンサデータの教師有データとを前記学習アルゴリズムに入力し、前記学習済みモデルを生成し、
前記時系列モデル生成部は、前記第二のラベルを時系列に並べ、第二の時系列モデルをさらに生成し、
前記予測部は、新たに取得した前記第二のセンサデータから、前記第二の時系列モデルを用いて、前記未来時刻の前記観測対象の状態として、前記第二のセンサデータの予測値をさらに算出する状態予測装置。
【請求項６】
請求項１記載の状態予測装置であって、
前記第一のセンサデータを用いて前記動作モード毎の時間幅を切出時間幅として算出する時間幅算出部を備え、
前記学習データ生成部は、前記時間幅算出部が算出した切出時間幅で前記第一のセンサデータを区切り、
前記時間幅算出部は、前記第一のセンサデータを取得する毎に、当該第一のセンサデータの取得開始時刻からの値の平均値を算出し、前記平均値が極値を取る時刻間を前記切出時間幅と算出する状態予測装置。
【請求項７】
請求項１から６のいずれか１項記載の状態予測装置であって、
前記予測値と前記新たに取得した第一のセンサデータとを比較し、評価するマッチング評価部をさらに備える状態予測装置。
【請求項８】
観測対象の状態を観測する第一の観測センサで得た第一のセンサデータを、所定の時間間隔で取得し、前記観測対象の動作モード毎に区切り、学習データを生成することを繰り返し、所定数の前記学習データを得る学習データ生成ステップと、
所定数の前記学習データそれぞれを予め定めたカテゴリに分類し、前記カテゴリ毎に予め定めたラベルを付し、教師有データを生成するクラスタリングステップと、
前記所定数の教師有データを、学習アルゴリズムに入力して学習し、当該学習アルゴリズムの学習済みモデルを生成するとともに、前記カテゴリそれぞれの代表データを生成する学習ステップと、
予め定めたモデル生成単位期間内に前記学習データに付された前記ラベルを、当該学習データの生成順に時系列に並べ、時系列モデルを生成する時系列モデル生成ステップと、
新たに取得した前記第一のセンサデータから、前記時系列モデルを用いて未来時刻の前記観測対象の状態として前記第一のセンサデータの予測値を算出する予測ステップと、を備える状態予測制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、状態予測技術に係り、特に、観測センサから継時的に出力されるセンサデータ（計測データ）に基づいて、未来の計測データを予測し、観測対象の状態を予測する技術に関する。
【背景技術】
【０００２】
観測センサから出力される計測データに基づいて観測対象の未来の状態を予測するためには、計測データの時系列モデルの作成が必要である。この時系列モデルの作成には、計測データの時間的な区切りを検出する必要がある。例えば、特許文献１には、「時系列データを、均等なサイズでオーバーラップするデータのセグメントに分割し、各々のセグメントに関し、セグメントのデータを表現する画像を生成し、時系列データを利用して、各々の画像に関連するトレンドを決定し、生成された各画像及び関連するトレンドを前記データセットとして保存する工程を含む（要約抜粋）。」技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１７−１５７２１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１に開示の技術では、計測データを時間的に均等なサイズのセグメントに分割している。ところが、観測対象が複数の動作モードを有し、出力される計測データが複雑に変化する場合、均等なサイズでのセグメント化は、必ずしも適切ではない。従って、状態予測の精度を高めることが難しい。【０００５】
本発明は、上記事情に鑑みてなされたもので、観測センサから出力される計測データ（観測対象に関連する時系列観測値）に基づいて、複数の動作モードを有する観測対象の状態を高い精度で予測する技術を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明は請求の範囲に記載の構成を有する。その一態様をあげるならば、本発明に係る状態予測制御装置は、観測対象の状態を観測する第一の観測センサで得た第一のセンサデータを、所定の時間間隔で取得し、前記観測対象の動作モード毎に区切り、学習データを生成する学習データ生成部と、生成した所定数の前記学習データそれぞれを予め定めたカテゴリに分類し、前記カテゴリ毎に予め用意されたラベルを付し、教師有データを生成するクラスタリング部と、前記所定数の教師有データを、学習アルゴリズムに入力して学習し、当該学習アルゴリズムの学習済みモデルを生成するとともに、前記カテゴリそれぞれの代表データを生成する学習部と、予め定めたモデル生成単位期間内に前記クラスタリング部で前記学習データに付された前記ラベルを、当該学習データの生成順に時系列に並べ、時系列モデルを生成する時系列モデル生成部と、新たに取得した第一のセンサデータから、前記時系列モデルを用いて未来時刻の前記観測対象の状態として前記第一のセンサデータの予測値を算出する予測部と、を備える。
【発明の効果】
【０００７】
本発明によれば、観測センサから出力される計測データに基づいて、複数の動作モードを有する観測対象の状態を高い精度で予測できる。なお、上記した以外の課題、解決手段は以下の実施形態により明らかにされる。
【図面の簡単な説明】
【０００８】
【図１】第一の実施形態の状態予測装置のハードウェア構成図である。
【図２】第一の実施形態の状態予測装置の機能ブロック図である。
【図３】（ａ）は、第一の実施形態のラベルデータ例を、（ｂ）および（ｃ）は、第一の実施形態の教師有データ例を、（ｄ）および（ｅ）は、第一の実施形態の代表データベース例を、それぞれ説明するための説明図である。
【図４】（ａ）および（ｂ）は、第一の実施形態の学習の概要を説明するための説明図である。
【図５】第一の実施形態のディープラーニングの概要を説明するための説明図である。
【図６】第一の実施形態の時系列モデルデータ例を説明するための説明図である。
【図７】第一の実施形態の状態予測処理のフローチャートである。
【図８】第一の実施形態の前処理と学習処理のフローチャートである。
【図９】第一の実施形態の時間幅算出処理のフローチャートである。
【図１０】（ａ）は、第一の実施形態の計測データの波形データ例の、（ｂ）は、第一の実施形態の計測データのヒストグラム例の、それぞれグラフである。
【図１１】（ａ）および（ｂ）は、第一の実施形態の時間幅算出処理を説明するためのグラフである。
【図１２】（ａ）および（ｂ）は、第一の実施形態の時間幅算出処理を説明するためのグラフである。
【図１３】（ａ）および（ｂ）は、第一の実施形態の時間幅算出処理を説明するためのグラフである。
【図１４】第一の実施形態の予測処理のフローチャートである。
【図１５】第一の実施形態の変形例を説明するためのグラフである。
【図１６】（ａ）および（ｂ）は、第一の実施形態の変形例を説明するためのグラフである。
【図１７】（ａ）および（ｂ）は、第一の実施形態の変形例を説明するためのグラフである。
【図１８】第二の実施形態の各計測データの切出時間幅を説明するための説明図である。
【図１９】第二の実施形態のディープラーニングの概要を説明するための説明図である。
【図２０】第二の実施形態の学習処理のフローチャートである。
【図２１】第二の実施形態の予測処理のフローチャートである。
【発明を実施するための形態】
【０００９】
＜＜第一の実施形態＞＞
以下、本発明の好適な実施形態を、図面を参照しながら詳細に説明する。以下では、プラント機器を観測対象とする場合を例にあげて説明する。
【００１０】
本実施形態では、プラント機器の様々な部位に取り付けられた各種の観測センサから継時的に得られるセンサデータ（計測データ）を用い、計測データの変化の時系列モデルを生成する。そして、生成した時系列モデルを用いて、未来時刻の計測データを予測し、観測対象であるプラント機器の状態予測を行う。また、最終的には、予測した計測データに基づき、プラント機器の故障等の有無を判別する。
【００１１】
プラント機器は、制御により動作モードが変化し、自動的にフィードバック制御がなされる。例えば、プラント機器が備える生産容器には、液体や気体等の原材料がポンプ等により投入される。この場合の流入量を外部のプログラムから制御すると、制御に応じてポンプ用モータの電流、電圧が変化する。また、生産容器の温度を制御すると、温度およびヒーターや冷却用ポンプモータの電流、電圧が変化する。
【００１２】
このようなプラント機器の各部位の状態を観測するため、プラント機器の各部位には、電流、電圧や温度等を検出する各種センサ（観測センサ）が取り付けられる。
【００１３】
これらの観測センサから時系列に得られる計測データは、時間的に変化する波形形状を示す波形データである。プラント機器は、上述のように複数の動作モード間でフィードバック制御がなされるため、電流や電圧の変化は、非定常である。従って、この波形データにおいては、複数の動作モードが非定常に観測され、継時的に得られる計測データも非定常となる。
【００１４】
プラント機器が生産するモノに対応して外部のプログラムから動作モードの組み合わせや各動作モードの時間幅を自動的なフィードバックを含めて制御する場合に、このような事象は一般的に発生する。
【００１５】
なお、継時的に得られる計測データ（時系列データ）が非定常であるとは、計測データの平均や分散といった統計的なパラメータが一定ではないことや、不連続に変化することを意味する。
【００１６】
非定常な時系列データは、ＡＲＩＭＡ（Ａｕｔｏ  Ｒｅｇｒｅｓｓｉｖｅ  Ｉｎｔｅｇｒａｔｅｄ  Ｍｏｖｉｎｇ  Ａｖｅｒａｇｅ）等の統計的な手法によりモデル化し、予測することが困難である。本実施形態では、このような動作を行うプラント機器に取り付けられた観測センサから得られる計測データから動作モードの時系列モデルを作成し、対象物の予測値を算出することにより、状態を予測する。
【００１７】
ここで、「対象物」とは、観測センサのそれぞれが出力したセンサ情報（センサデータ；計測データ）から直接観測されるものを意味する。例えば、観測センサが温度センサの場合は、対象物は温度である。電流センサの場合は、電流、電圧センサの場合は、電圧、カメラ（ＣＭＯＳセンサ）である場合は、撮影対象のプラント機器の部位そのものもしくはプラント機器の生産物が対象物である。
【００１８】
［ハードウェア構成］
図１に本実施形態に係る状態予測装置のハードウェア構成を示す。図１に示すように本実施形態に係る状態予測装置１００は、ＣＰＵ１１２と、ＲＯＭ１１３と、ＲＡＭ１１４と、ファイルユニット１１５と、外部インタフェース（Ｉ／Ｆ）１１６と、モニタ１１７と、バス１１８と、を備える。
【００１９】
ＣＰＵ１１２は、演算処理を行う。ＲＯＭ１１３およびファイルユニット１１５は、データ記憶装置として機能する。ＲＯＭ１１３には、処理ルーチン等のプログラムが予め記憶される。なお、ファイルユニット１１５は、例えば、ＨＤＤにより構成されてもよい。モニタ１１７は、処理結果を表示する。ＣＰＵ１１２とＲＯＭ１１３とＲＡＭ１１４とファイルユニット１１５と、外部インタフェース１１６と、モニタ１１７とは、バス１１８によりコマンドやデータを授受可能に接続され、コンピュータとして機能する。
【００２０】
なお、本実施形態に係る状態予測装置１００には、以上の構成の他に、電源装置等、多数の電気系の構成要素が含まれているが、周知または一般的なものであるため詳細な説明を省略する。
【００２１】
また、外部Ｉ／Ｆ１１６には、状態予測装置１００が観測対象とするプラント機器の状態情報を取得するための観測センサ３００が接続される。観測センサ３００には、例えば、電流センサ３１１、電圧センサ３１２、温度センサ３１３、圧力センサ３１４、加速度センサ３１５、カメラ３１６等がある。接続される観測センサの種類は、これらに限定されない。
【００２２】
上述のように、本実施形態の状態予測装置１００は、複数の動作モードが非定常に観測される計測データの時系列モデルを作成し、対象物の予測値を算出する。
【００２３】
時系列モデルの作成にあたっては、観測センサ３００から出力される計測データを継続的に取得し、動作モード毎に区切り、動作モード毎に区切られた計測データの時系列モデルを生成する。また、予測値の算出は、時系列モデルの生成時に用いた計測データとは異なる、新たに取得した計測データを用いて行う。以下、これを実現する本実施形態の状態予測装置１００の機能構成を説明する。
【００２４】
［機能ブロック］
本実施形態の状態予測装置１００の機能ブロックを、図２に示す。上記機能を実現するため、本実施形態の状態予測装置１００は、前処理学習部２１０と、予測部２２０と、を備える。
【００２５】
前処理学習部２１０は、所定の時間間隔で継時的に取り込んだ計測データから学習データを生成し、学習アルゴリズムに入力して学習し、学習済みモデル、分類毎の代表データおよび時系列モデルを生成する。学習データは、計測データを動作モード毎に区切り、セグメント化したものである。本実施形態では、動作モード毎に、それぞれ最適な時間幅でセグメント化し、学習データとする。
【００２６】
これを実現するため、前処理学習部２１０は、時間幅算出部２１１と、学習データ生成部２１２と、クラスタリング部２１３と、学習部２１４と、時系列モデル生成部２１５と、を備える。
【００２７】
また、予測部２２０は、前処理学習部２１０と同様に、観測センサ３００から得た計測データを、動作モード毎に最適な時間幅で区切り、前処理学習部２１０が生成した学習済みモデル、分類毎の代表データおよび時系列モデルを用いて未来時刻の計測データを予測、評価する。
【００２８】
これを実現するため、予測部２２０は、時間幅算出部２２１と、予測用データ生成部２２２と、分類決定部２２３と、予測値算出部２２４と、マッチング評価部２２５と、を備える。
【００２９】
上記各機能は、ＲＯＭ１１３やファイルユニット１１５に予め記憶されたプログラムをＣＰＵ１１２がＲＡＭ１１４にロードして実行することにより実現される。
【００３０】
また、ファイルユニット１１５またはＲＯＭ１１３には、上記処理を実現するために必要なプログラムおよびデータが記憶される。本実施形態では、例えば、ラベル記憶部２３０と、学習データセット記憶部２４０と、学習済みモデル記憶部２５０と、代表データ記憶部２６０と、時系列モデル記憶部２７０と、の各記憶部が設けられる。以下、これらがファイルユニット１１５に設けられる場合を例にあげて説明する。
【００３１】
［前処理学習部］
まず、前処理学習部２１０および予測部２２０の各部について、詳細を説明する。
【００３２】
前処理学習部２１０の時間幅算出部２１１は、計測データを用いて、観測対象の各動作モードの時間幅を算出する。算出される時間幅は、計測データを動作モード毎にセグメント化する際の時間幅である。このため、以下、算出される時間幅を切出時間幅と呼ぶ。
【００３３】
本実施形態では、計測データを取得する毎に、計測値のヒストグラムを作成し、その分布が、例えば、正規分布等の一塊の山状の分布を示す時間幅を、切出時間幅として算出する。算出手法の詳細は、後述する。なお、本実施形態では、後述のクラスタリング部２１３の処理のため、所定数のデータを必要とする。このため、これを十分満たす期間の計測データを取得後、順次、切出時間幅を算出してもよい。
【００３４】
学習データ生成部２１２は、計測データを動作モード毎に切り取り、学習アルゴリズムに入力する学習データを生成する。本実施形態では、時間幅算出部２１１が算出した切出時間幅で計測データを区切り、学習データを生成する。本実施形態では、学習データとして、計測データの時間的な変化を示す波形データを生成する。なお、さらに、切出時間幅の計測データの、各値の発生頻度の分布を示すヒストグラム（ヒストグラム）を用いてもよい。波形データは、各計測データの値が現れる順序を示し、ヒストグラムは、各計測データの値の割合を示す。
【００３５】
クラスタリング部２１３は、所定数の学習データが収集されると、動作モード毎に切出した計測データから生成した各学習データから、教師有データを生成する。本実施形態では、予め定めた手法に従って、各学習データをクラスタリング（カテゴライズ、分類）し、学習データと分類結果とを対応づけた教師有データを生成する。生成された教師有データは、学習データセット記憶部２４０に記憶する。
【００３６】
クラスタリングを実行する際に収集する学習データの数（上述の所定数）は、観測するデータ種で異なる。少なくとも、一通りの動作モードを検出可能な数の学習データを収集する。一通りの動作モードが短時間で現れるデータ種の場合は、例えば、数十秒分の計測データから生成された学習データを用い、一通りの動作モードが現れるのに長時間かかるデータ種の場合は、数十日分の計測データから生成された学習データを用いる。収集数は、データ種に応じて、予め定めておく。
【００３７】
クラスタリングには、例えば、ｋ−ｍｅａｎｓ法、ｋ−ｍｅａｎｓ法を改良した方法等を用いる。分類（カテゴリ）数、各カテゴリの識別情報であるラベルは、予め用意し、ラベル記憶部２３０に記憶しておく。
【００３８】
ラベル記憶部２３０に記憶されるラベルデータ２３１例を図３（ａ）に示す。学習データが波形データとヒストグラムとの場合、それぞれについて、ラベルが用意される。ここでは、一例として、８個のラベルが用意される場合、すなわち、８つのカテゴリにクラスタリングする場合のラベル例を示す。
【００３９】
ここで、ｋ−ｍｅａｎｓ法を用い、８つのカテゴリにクラスタリングする場合を例にあげて説明する。なお、動作モードが既知でない場合には、本例のように予め初期値として定めた８つのカテゴリへのクラスタリングが適切ではなく、その結果、後に示す学習部２１４の学習アルゴリズム（ディープラーニングネット）２５１の学習が収束しなくなったり、マッチング評価部２２５の評価結果が所定の範囲を超えてしまい、再学習が必要となったりすることがある。その場合は、ここでの処理（クラスタリング処理）に戻り、ラベル数（カテゴリ数）を増減し、クラスタリングし直す。例えば、カテゴリ数を１０に増やし、１０のカテゴリにクラスタリングをし直したり、ラベル数を６に減らし、６のカテゴリにクラスタリングし直したりする。
【００４０】
まず、クラスタリング部２１３は、各学習データに、８個のラベルを適当に割り振る。同じラベルを付された学習データは、１つのカテゴリに分類されたということである。
【００４１】
次に、カテゴリ毎に、各学習データに対応するヒストグラム特徴量（度数、平均値、分散値、最頻値の度数、中央値等）の重心値を計算する（第一のステップ）。得られた８つの重心値に、それぞれ、当該カテゴリのラベルを対応づける。
【００４２】
その後、各学習データのヒストグラム特徴量について、それぞれ、８個の重心値との距離を計算し、最も近い重心値に対応づけられているラベルを付け直す（第二のステップ）。
【００４３】
第一のステップと第二のステップとを、重心値やラベル付が変化しなくなるまで、あるいは、重心値の変化が所定範囲以下になるまで、繰り返す。最終的に各学習データに付されたラベルに対応するカテゴリが、当該学習データの分類結果（クラスタリング結果）である。すなわち、各学習データと、当該学習データに付されたラベルの組が、教師有データであり、学習データセット記憶部２４０に記憶される。
【００４４】
このとき、学習データセット記憶部２４０に格納される教師有データ群の構成例を、図３（ｂ）および図３（ｃ）に示す。本実施形態では、教師有データ２４１として、図３（ｂ）に示す波形データの教師有データ２４１ｗと、図３（ｃ）に示す、ヒストグラムの教師有データ２４１ｈとがそれぞれ記憶される。
【００４５】
教師有データ２４１ｗは、取得時刻２４２ｗと、学習データ（波形データ）２４３ｗと、ラベル２４４ｗと、が対応づけて記憶される。取得時刻２４２ｗは、例えば、各学習データ２４３ｗの切出時間幅の最初の時刻等とする。また、教師有データ２４１ｈは、それぞれ、取得時刻２４２ｈと、学習データ（ヒストグラムデータ）２４３ｈと、ラベル２４４ｈと、が対応づけて記憶される。
【００４６】
なお、以下、教師有データ２４１ｗと教師有データ２４１ｈとを、区別する必要がない場合は、教師有データ２４１で、また、各項目も、取得時刻２４２、学習データ２４３、ラベル２４４、で代表する。他のデータについても、同様とする。
【００４７】
学習部２１４は、学習用データセットを、予め定めた学習アルゴリズムに入力し、学習済みモデルを生成する。本実施形態では、予め定めた学習アルゴリズムとして、例えば、ディープラーニングを用いる。すなわち、学習部２１４は、教師有データ群である学習用データセットによりディープラーニングネットを学習する。
【００４８】
本実施形態では、図４（ａ）に示すように、前処理学習部２１０において、学習データ２４３群で構成される学習データセット２４３ａを、学習アルゴリズム２５１に入力することにより、学習済みのディープラーニングネット２５２を生成する。以下、学習済みのディープラーニングネットを学習済みモデル２５２とも呼ぶ。そして、図４（ｂ）に示すように、予測部２２０において、その学習済みモデル２５２を用いて、予測用データ２４３ｂを分類２４３ｃする。すなわち、各予測用データ２４３ｂに、ラベル２４４を付す。
【００４９】
ここで、本実施形態の学習部２１４が生成する学習済みのディープラーニングネット２５２の一例を説明する。図５は、学習済みのディープラーニングネット２５２の一例を説明するための図である。
【００５０】
本図に示すように、学習部２１４は、例えば、切出時間幅で切出した計測データから得た学習データ２４３群で構成される学習データセット２４３ａを、予め設定した時間幅に揃え、学習アルゴリズム２５１のデータ入力ウィンドウに入力する。なお、データの時間幅は、時間幅前後のデータを付け加えたり、データを繰り返したり等により、長さを揃える。
【００５１】
図５に示すように、データ入力ウィンドウに入力されたデータは、ディープラーニングネットの畳み込み層、Ｍａｘｏｕｔ層、Ｐｏｏｌｉｎｇ層の各層で特徴抽出され、Ｆｕｌｌｙｃｏｎｎｅｃｔｉｏｎ層に入力される。そして、Ｆｕｌｌｙｃｏｎｎｅｃｔｉｏｎ層で特徴の組合せから分類結果として出力される。なお、畳み込み層、Ｍａｘｏｕｔ層、Ｐｏｏｌｉｎｇ層は、その３つ組の階層を何回か重ねる構成としてもよい。
【００５２】
本実施形態では、最終的に学習済みモデル２５２として、学習済みのディープラーニングネット２５２の階層数、ノード数、各ノード間のつながり、重み等が記憶される。
【００５３】
また、学習部２１４の学習後、または、クラスタリング部２１３によるクラスタリング後、ラベル２４４毎の代表データを生成し、ラベル２４４に対応づけて、代表データ記憶部２６０に記憶する。代表データは、ラベル２４４毎に、当該ラベルが付与された波形データおよび／またはヒストグラムの平均を算出することにより得る。平均は、例えば、波形データの場合、同じラベルが付与された波形データを重ね合わせて、波形数で割ることにより算出される。
【００５４】
代表データ記憶部２６０に記憶される代表データベース２６１（２６１ｗ、２６１ｈ）の例を図３（ｄ）および図３（ｅ）に示す。代表データベース２６１ｗでは、各ラベル２４４ｗに対応づけて、波形データの代表データである代表波形データ２６３ｗが記憶される。また、代表データベース２６１ｈでは、同様に、ラベル２４４ｈに対応づけて、代表ヒストグラム２６３ｈが記憶される。
【００５５】
時系列モデル生成部２１５は、各ラベル２４４を、時系列に並べ、時系列モデルを生成する。本実施形態では、所定期間（モデル生成単位期間）毎に、各学習データに付与されたラベル２４４を時系列に並べ、時系列モデルとして生成する。本実施形態では、学習データセット記憶部２４０に記憶される、各学習データのラベル２４４を、当該学習データの生成順に並べることにより、時系列モデルを生成する。また、生成した時系列モデルは、時系列モデル記憶部２７０に記憶する。
【００５６】
なお、モデル生成単位期間は、例えば、１時間、１日等、計測データの種類に応じて予め定められる。
【００５７】
時系列モデル記憶部２７０に格納される時系列モデルデータ２７１の例を図６に示す。時系列モデルデータ２７１は、時系列モデル２７３毎に、時系列モデル２７３を識別する時系列モデル識別情報（モデルＩＤ）２７２が対応づけて記憶される。時系列モデル２７３は、時系列に並べられたラベル２４４の並び順が異なると異なるものとされ、それぞれ、モデルＩＤ２７２が付与される。
【００５８】
なお、図６に示す例では、１つでもラベル２４４の並び順が異なる時系列モデル２７３は、異なる時系列モデルとして異なるモデルＩＤ２７２を付与しているが、これに限定されない。複数の類似のラベル２４４の並び順の時系列モデル２７３を、１つのグループとして、同じモデルＩＤ２７２を付与するよう構成してもよい。
【００５９】
また、学習データとして、波形データとヒストグラムとの両方を用いる場合、両データの並び順を併せて、１つのラベルの並び順として扱う。例えば、図５の例では、モデルＩＤ２７２のＴＳ＿１とＴＳ＿２とは、ヒストグラムのラベルの並び順は同じであるが、波形データのラベル並び順が異なるため、これらは、異なるラベル並び順として、異なるモデルＩＤ２７２が付与される。
【００６０】
［予測部］
次に、予測部２２０の各構成について説明する。予測部２２０の時間幅算出部２２１と、予測用データ生成部２２２とは、それぞれ、前処理学習部２１０の時間幅算出部２１１と、学習データ生成部２１２と同様の機能を有する。ただし、予測用データ生成部２２２が生成するデータを、予測用データと呼ぶ。
【００６１】
予測部２２０では、前処理学習部２１０が、生成する学習データと同種の学習データを生成する。すなわち、波形データとヒストグラムの２種の学習データが生成される場合、予測部２２０は、波形データとヒストグラムとの２種の予測用データを生成する。一方、波形データのみの学習データが生成される場合は、波形データのみの予測用データを生成する。
【００６２】
分類決定部２２３は、生成された予測用データを、学習アルゴリズムに入力し、当該予測用データの分類を決定し、当該分類に設定されたラベル２４４を付与する。本実施形態では、このとき、学習済みモデル２５２を用いる。
【００６３】
予測値算出部２２４は、ラベル２４４を付与された予測用データが所定数、収集されると、それらを用いて、未来時刻の計測データを予測する。本実施形態では、まず、予測ラベル特定部としての処理を行う。すなわち、収集された予測用データのラベル２４４を、時系列に並べ、時系列モデル記憶部２７０に記憶される時系列モデル２７３と比較し、類似の時系列モデル２７３を特定する。そして、特定された時系列モデル２７３に従って、未来時刻の予測用データのラベル２４４を特定する。
【００６４】
そして、予測値算出部２２４は、特定されたラベル２４４に対応づけて、代表データ記憶部２６０に記憶される代表データ（波形データおよびヒストグラム）２６３を、予測値（期待値）として算出する。
【００６５】
マッチング評価部２２５は、算出された予測値と、当該未来時刻の実測された計測データ（実測データ）とを比較し、評価（マッチング）する。評価は、例えば、偏差σがわかっている場合、予測値からの３σの範囲に入っているか否かで評価を行ってもよい。
【００６６】
評価は、例えば、前処理学習部２１０が生成した学習済みモデルの妥当性について行う。評価においては、予測値と実測データとの差異を算出し、差異が所定の範囲内に収まっていれば、妥当と判断し、所定の範囲を超えている場合、再学習の必要あり、と判断する。
【００６７】
なお、学習済みモデルの妥当性ではなく、故障または異常の発生の有無を判断してもよい。例えば、予測値と実測データとの差異を算出し、差異が所定の範囲内を超えている場合、故障または異常発生の可能性有り、と判断する。
【００６８】
［状態予測処理］
次に、本実施形態の状態予測装置１００による状態予測処理の流れを説明する。図７は、本実施形態の状態予測処理の全体の流れのフローである。
【００６９】
本図に示すように、まず、前処理学習部２１０が、前処理と学習処理を行う（ステップＳ１１００）。ここでは、前処理学習部２１０は、前処理と学習処理として、クラスタリング、学習、時系列モデル生成処理を行い、学習済みモデル、時系列モデル、各カテゴリの代表データの生成を行う。その後、予測部２２０が、最新の計測データと学習済みモデル、代表データおよび時系列モデルを用いて予測処理を行う（ステップＳ１２００）。なお、前処理と学習処理と、予測処理とは、独立して行われる。
【００７０】
以下、前処理と学習処理と、予測処理と、について、それぞれ、詳細を説明する。ここでは、波形データとヒストグラムとの学習データを生成する場合を例にあげて説明する。
【００７１】
［前処理と学習処理］
図８は、本実施形態の前処理学習部２１０による前処理と学習処理との処理フローである。本実施形態では、前処理および学習処理において、学習済みモデルと、各分類の代表データと、時系列モデルとを生成する。この前処理および学習処理は、ユーザからの開始の指示や予め設定されたタイムスケジュールや計測の区切り信号などを契機に実行される。また、観測対象のプラント機器が正常な状態の時の計測データを用いて実行される。前処理学習部２１０は、所定の時間間隔で、観測センサ３００から計測データを受信するものとする。ここで、ｍは、カウンタであり、切出時間幅ＴＷの個数をカウントする。
【００７２】
まず、前処理学習部２１０は、カウンタｍを初期化（ｍ＝１）する（ステップＳ１１００）。
【００７３】
時間幅算出部２１１は、計測データを取り込み（ステップＳ１１０１）、取り込むごとに、ｍ番目の切出時間幅ＴＷ（ｍ）の算出処理を行う（ステップＳ１１０２）。切出時間幅ＴＷ（ｍ）が決定されるまで（ステップＳ１１０３）、時間幅算出部２１１は、計測データの取り込みおよび切出時間幅ＴＷ（ｍ）の算出処理を繰り返す。切出時間幅決定処理の詳細は、後述する。
【００７４】
ｍ番目の切出時間幅ＴＷ（ｍ）が決定したら、学習データ生成部２１２は、学習データを生成する（ステップＳ１１０４）。ここでは、まず、取り込んだ計測データから切出時間幅ＴＷ（ｍ）分の波形データとヒストグラムとを切り出す。そして、切出時間幅ＴＷ（ｍ）で切り出した波形データおよびヒストグラムに対し、所定のデータ値の範囲（ダイナミックレンジ）に収める正規化等の処理を施し、学習データとして生成する。なお、このとき、周期的な変化、季節変動等の要因を除去してもよい。なお、学習データを生成する毎に、直前までに取り込んだ計測データは破棄する。
【００７５】
前処理学習部２１０は、以上の処理を繰り返し、所定数の学習データを収集する。所定数の学習データが生成されたら(ステップＳ１１０５)、クラスタリング部２１３は、収集した所定数の学習データ（波形データおよびヒストグラム）を、クラスタリングし(ステップＳ１１０６)、それぞれを、教師有データとする。これにより、教師有データのセット、すなわち、学習データセット２４３ａが生成される。
【００７６】
なお、所定数の学習データが未収集の場合、ｍを１インクリメントし(ステップＳ１１１０)、ステップＳ１１０１へ戻り、処理を繰り返す。
【００７７】
学習部２１４は、生成された学習データセット２４３ａを、学習アルゴリズム２５１に入力して学習させ(ステップＳ１１０７)、学習済みモデル２５２を生成する。このとき、各分類の代表データ２６３も生成する。
【００７８】
ここで、時系列モデル生成部２１５は、予め定めたモデル生成単位期間ＴＤＴ分の計測データのクラスタリングを終えたか否かを判別する（ステップＳ１１０８）。ここでは、処理開始からモデル生成単位期間ＴＤが経過したか否かを判別する。
【００７９】
モデル生成単位期間ＴＤを経過していない場合（ステップＳ１１０８；Ｎｏ）は、ステップＳ１１００へ戻り、処理を繰り返す。
【００８０】
一方、モデル生成単位期間ＴＤを経過した場合（ステップＳ１１０８；Ｙｅｓ）、時系列モデル生成部２１５は、時系列モデルを生成し（ステップＳ１１０９）、処理を終了する。
【００８１】
時系列モデル生成部２１５は、例えば、モデル生成単位期間ＴＤのラベル２４４の並び順を時系列に抽出する。そして、既に生成された時系列モデル２７３を参照し、同じ並び順の時系列モデル２７３の有無を判別する。同じ並び順の時系列モデル２７３が有る場合は、そのまま処理を終了する。一方、同じ並び順の時系列モデル２７３が無い場合は、新たに、モデルＩＤ２７２を設定し、ラベル２４４の並び順を時系列モデル２７３として時系列モデル記憶部２７０に登録し、処理を終了する。
【００８２】
なお、所定の類似度のラベル２４４の並び順に、同じモデルＩＤ２７２を付与するグループ管理の場合は、以下の通りに行う。まず、予め定めた類似度の範囲内の時系列モデル２７３が登録されているか否かを判別する。そして、登録されている場合は、そのまま処理を終了する。一方、未登録の場合は、新たにモデルＩＤ２７２を設定し、得られた時系列モデル２７３を登録する。
【００８３】
［時間幅算出処理］
次に、上記Ｓ１１０２の時間幅算出処理の流れを、図９に従って説明する。ここでは、ステップＳ１１０１の計測データの取り込み処理も含めて説明する。本実施形態では、計測データを得る毎に、当該計測データ値毎の頻度を示すヒストグラムを作成する。そして、そのヒストグラムが、一塊の山状の分布に近い分布を示す期間を１つの動作モード期間とし、波形データおよびヒストグラムを切り取る切出時間幅ＴＷ（ｍ）とする。具体的には、計測データの収集開始からヒストグラムの分散値の逆数が極大値を取る時刻までを、切出時間幅ＴＷ（ｍ）とする。
【００８４】
ここでは、計測データは、Δｔ毎に収集されるものとする。また、前回（ｍ−１回目）の切出時間幅ＴＷ（ｍ−１）を決定した際の、最終時刻からΔｔを経過した時刻を、計測データの収集開始時刻ｔ０とする。ｎは、カウンタである。また、Ｎ個の計測データを収集しても、切出時間幅ＴＷ（ｍ）が決定されない場合、予め定めた時間幅を切出時間幅ＴＷ（ｍ）とする。
【００８５】
まず、初期処理を行う（ステップＳ１３０１）。ここでは、時刻ｔ０の計測データＤｔ（０）と、時刻ｔ１（ｔ１＝ｔ０＋Δｔ）の計測データＤｔ（１）とを取得し、例えば、ＲＡＭ１１４等に、記憶する。そして、カウンタｎを２とする。
【００８６】
次に、時刻ｔｎ（ｔｎ＝ｔ０＋Δｔ・ｎ）の計測データＤｔ（ｎ）を取得する（ステップＳ１３０２）。
【００８７】
Ｄｔ（０）からＤｔ（ｎ）までの間に取得した計測データを用いてヒストグラムを作成する（ステップＳ１３０３）。
【００８８】
そして、作成したヒストグラムの分散の逆数Ｖ（ｎ）を算出し、記憶する（ステップＳ１３０４）。なお、分散は、算出対象の計測データの平均値を求め、その平均値と各計測データの値との差の二乗和を全データ数で除算することにより得られる。
【００８９】
そして、Ｖ’（ｎ）とＶ”（ｎ）とを算出する（ステップＳ１３０５）。Ｖ’（ｎ）は、１回前に取得した計測データＤｔ（ｎ−１）までのヒストグラムの分散の逆数Ｖ（ｎ−１）からＶ（ｎ）を減算した値である。また、Ｖ”（ｎ）は、Ｖ’（ｎ−１）からＶ’（ｎ）を減算した値である。
【００９０】
そして、Ｖ’（ｎ）の絶対値が予め定めた閾値Ｔｈより小さく、かつ、Ｖ”（ｎ）が負であるという条件を満たすか否かを判別する（ステップＳ１３０６）。
【００９１】
ステップＳ１３０６の判別結果が条件を満たさない場合、所定数の計測データを処理したか否を判別する(ステップＳ１３１０)。所定数の計測データの処理を終えていない場合は、ｎを１増加させ（ステップＳ１３１１）、ステップＳ１３０２へ戻り、処理を繰り返す。
【００９２】
ステップＳ１３０６の判別結果が条件を満たす場合、および、ステップＳ１３１０において、所定数の計測データの処理を終えている場合、その時点の処理対象の計測データＤｔ（ｎ）を取得した時刻ｔｎを時間区切りとする（ステップＳ１３０７）。そして、ｔ０から時刻ｔｎまでを切出時間幅Ｔ（ｍ）とする。
【００９３】
そして、次の時間幅設定処理のために、ｔ０を設定し（ステップＳ１３０８）、処理を終了する。ここでは、ｔ０は、時刻ｔｎに、Δｔを加算した時刻とする。
【００９４】
以上の時間幅算出部２１１による時間幅算出処理を、具体例で説明する。例えば、計測データ３２０が、図１０（ａ）に示すように時間的に変化するものとする。なお、計測データ３２０は、上述のようにΔｔ毎に取得される離散値である。また、図１０（ｂ）は、計測データ３２０の各値の頻度をグラフ化したヒストグラム３３０である。
【００９５】
時間幅算出部２１１は、図１１（ａ）に示すように、所定の時間間隔Δｔで、計測データ３２０を取得する毎に、図１１（ｂ）に示すように、ヒストグラム３３０を作成する。そして、ヒストグラム３３０を作成する毎に、その分散の逆数Ｖ（ｎ）が極大値を取るか否かを判別する。具体的には、Ｖ’（ｎ）、Ｖ”（ｎ）を算出し、ステップＳ１３０６の条件を判断する。
【００９６】
ステップＳ１３０６の条件を満たす、すなわち、ヒストグラムの分散の逆数Ｖ（ｎ）が極大値を取った時点ｔ１までを、１番目の切出時間幅ＴＷ（１）とする。
【００９７】
そして、ｔ１にΔｔを加えた時刻、すなわち、１番目の切出時間幅ＴＷ（１）の最後の計測データ３２０の次の計測データ３２０を取得する時刻を、次の処理ループのデータの収集開始時刻ｔ０とする。そして、新たに、計測データ３２０の収集およびヒストグラム３３０の作成を開始する。これは、２番目の切出時間幅ＴＷ（２）を決定するための処理である。
【００９８】
時間幅算出部２１１は、２番目の切出時間幅ＴＷ（２）の判定を開始する。この場合の計測データ３２０の変化の様子を図１２（ａ）に示す。また、ヒストグラムの変化の様子を図１２（ｂ）に示す。
【００９９】
この場合も、ｔ１＋Δｔをｔ０とし、図１２（ａ）に示すように、ｔ０からΔｔ毎に計測データ３２０を取得する。そして、図１２（ｂ）に示すように、取得する毎に、ヒストグラム３３０を作成し、その都度、分散の逆数Ｖ（ｎ）が極大値を取るか否かを判別する。そして、極大値をとった時点ｔ２までを、２番目の切出時間幅ＴＷ（２）とする。
【０１００】
そして、ｔ２にΔｔを加えた時刻、すなわち、２番目の切出時間幅ＴＷ（２）の最後の計測データ３２０の次の計測データ３２０を取得する時刻をｔ０とし、新たに、計測データ３２０の収集およびヒストグラムの作成を開始する。これは、３番目の切出時間幅ＴＷ（３）を決定するための処理である。
【０１０１】
時間幅算出部２１１は、３番目の切出時間幅ＴＷ（３）の判定を開始する。この場合、の計測データ３２０の変化の様子を図１３（ａ）に示す。また、ヒストグラムの変化の様子を図１３（ｂ）に示す。
【０１０２】
この場合も、ｔ２＋Δｔをｔ０とし、図１３（ａ）に示すように、ｔ０からΔｔ毎に計測データ３２０を取得する。そして、図１３（ｂ）に示すように、取得する毎に、ヒストグラム３３０を作成し、その都度、分散の逆数Ｖ（ｎ）が極大値を取るか否かを判別する。そして、極大値を取る時点ｔ３までを、３番目の切出時間幅ＴＷ（３）と決定する。
【０１０３】
時間幅算出部２１１は、以上の処理を、繰り返し、Ｍ番目の切出時間幅ＴＷ（Ｍ）まで作成する。
【０１０４】
［予測処理］
次に、予測部２２０による、予測処理の流れを説明する。図１４は、本実施形態の予測処理の処理フローである。予測処理は、ユーザからの開始指示や予め設定されたタイムスケジュールや計測の区切り信号などを契機に開始される。
【０１０５】
予測部２２０の時間幅算出部２２１は、学習処理同様、計測データを取り込むごとに（ステップＳ１２０１）、切出時間幅の算出を試みる（ステップＳ１２０２）。切出時間幅を決定するまで、ステップＳ１２０１およびＳ１２０２を繰り返す。切出時間幅の算出手法は、前処理学習部２１０の時間幅算出部２１１による処理と同様である。
【０１０６】
時間幅算出部２２１が切出時間幅を決定したら（ステップＳ１２０３；Ｙｅｓ）、予測用データ生成部２２２は、予測用データを生成する（ステップＳ１２０４）。予測用データの生成手法は、基本的に前処理学習部２１０の学習データ生成部２１２による学習データ生成手法と同様である。例えば、所定の値の幅（ダイナミックレンジ）に値が収まるよう、正規化処理を行う。また、周期的な変化、季節変動等の要因を除去してもよい。
【０１０７】
次に、分類決定部２２３は、生成された予測用データ（波形データおよびヒストグラム）を、学習アルゴリズムと学習済みモデルとを用いて分類する(ステップＳ１２０５)。ここでは、分類されたカテゴリに設定されているラベル２４４を付与する。
【０１０８】
その後、予測値算出部２２４は、それまでに分類した各切出時間幅の予測用データのラベル２４４を、時刻順に並べる。そして、記憶されている時系列モデル２７３と比較し、最も類似度の高い並び順を含む時系列モデル（グループ）２７３を、最類似時系列モデルとして特定する（ステップＳ１２０６）。
【０１０９】
例えば、波形データの予測用データから得た、ラベル２４４の並び順がＡｗ、Ｂｗ、Ｃｗであり、ヒストグラムの予測用データから得たラベル２４４の並び順がＥｈ、Ｂｈ、Ｃｈである場合、図６に示す時系列モデルデータ２７１では、モデルＩＤ２７２がＴＳ＿４の時系列モデル２７３が最も類似度が高い。従って、モデルＩＤ２７２がＴＳ＿４の時系列モデル２７３が、最類似時系列モデルとして特定される。
【０１１０】
また、波形データの予測用データから得た、ラベル２４４の並び順がＢｗ、Ｄｗ、Ａｗであり、ヒストグラムの予測用データから得たラベル２４４の並び順がＡｈ、Ａｈ，Ａｈである場合、図６に示す時系列モデルデータ２７１では、モデルＩＤ２７２がＴＳ＿３の時系列モデル２７３が、この並び順を含み、最も類似度が高い。従って、モデルＩＤ２７２がＴＳ＿３の時系列モデル２７３が、最類似時系列モデルとして特定される。
【０１１１】
なお、学習データおよび予測用データとして、波形データのみを用いる場合は、波形データの時系列モデルのみで最類似時系列モデルを決定する。
【０１１２】
そして、予測値算出部２２４は、決定した最類似時系列モデルを用い、以降の波形データおよびヒストグラムを予測する（ステップＳ１２０７）。ここでは、決定した最類似時系列モデルに従って、以降のラベル（予測ラベル）を特定する。そして、その予測ラベルに対応づけて記憶される代表データ２６３（代表波形データ２６３ｗおよび／または代表ヒストグラム２６３ｈ）を、予測値とする。
【０１１３】
例えば、最類似時系列モデルが、図６に示すモデルＩＤＴＳ＿４の時系列モデル２７３である上記の例の場合、波形データに関し、次のラベル２４４は、Ｄｗである。従って、ラベルＤｗに対応づけて記憶される代表波形データ２６３ｗが予測値として算出される。
【０１１４】
その後、マッチング評価部２２５は、以降の観測センサ３００から出力される実際の計測データと、ステップＳ１２０７で得た予測値とを比較し、評価する(ステップＳ１２０８)。
【０１１５】
以上説明したように、本実施形態の状態予測装置１００は、観測対象の状態を観測する第一の観測センサ３００で得た第一のセンサデータを、所定の時間間隔で取得し、観測対象の動作モード毎に区切り、学習データを生成する学習データ生成部２１２と、生成した所定数の前記学習データそれぞれを予め定めたカテゴリに分類し、カテゴリ毎に予め用意されたラベルを付し、教師有データを生成するクラスタリング部２１３と、所定数の教師有データを、学習アルゴリズムに入力して学習し、当該学習アルゴリズムの学習済みモデル２５２を生成するとともに、前記カテゴリそれぞれの代表データを生成する学習部２１４と、予め定めたモデル生成単位期間内にクラスタリング部２１３で学習データに付されたラベル２４４を、当該学習データの生成順に時系列に並べ、時系列モデル２７３を生成する時系列モデル生成部２１５と、新たに取得した第一のセンサデータから、時系列モデル２７３を用いて未来時刻の観測対象の状態として第一のセンサデータの予測値を算出する予測部２２０と、を備える。
【０１１６】
このため、本実施形態によれば、観測対象が複数の動作モードを有し、計測データが複雑に変化する場合であっても、動作モード毎に区切って学習データを生成するため、適切に状態を反映して区切られたデータから状態予測のための学習データを生成できる。適切な学習データで学習された結果から予測のための学習済みモデル、代表データおよび時系列モデルが生成されるため、精度の高い予測値を得ることができる。
【０１１７】
また、本実施形態の状態予測装置１００は、第一のセンサデータを用いて動作モード毎の時間幅を切出時間幅として算出する時間幅算出部をさらに備えてもよい。時間幅算出部を備えることにより、本実施形態の状態予測装置１００は、実際の計測データから動作モードの区切りを検出できる。このため、より高い精度で、動作モードを検出でき、検出結果に従って、計測データを区切ることができる。従って、より精度の良い予測値を得ることができる。
【０１１８】
このように、本実施形態によれば、観測対象物に関連する物理量（計測データ）の時系列観測値に基づいて、観測対象物の状態を精度よく予測できる。
【０１１９】
また、本実施形態では、ヒストグラムについても、併せてクラスタリングし、学習させて分類してもよい。計測データには、その性質により、所定の時間幅内の、値が現れる順序を示す波形（波形データ）より、値の割合を示す分布（ヒストグラム）に意味があるタイプの計測データがある。ヒストグラムについても、波形データ同様の処理を行うことで、このようなタイプの計測データについても、その特徴に応じて、的確に分類できる。すなわち、本実施形態によれば、多面的に計測データを解析できる。
【０１２０】
＜変形例＞
なお、上記実施形態では、時間幅算出処理において、計測データの収集開始から、ヒストグラムの分散値の逆数が極大値を取る時刻までを、切出時間幅ＴＷ（ｍ）に決定している。しかしながら、切出時間幅の決定手法は、これに限定されない。例えば、計測データを取得する毎に、取得開始時刻からの計測データ値（計測値）の平均値を算出し、平均値の変化が極値を取る時刻間を、切出時間幅ＴＷ（ｍ）としてもよい。
【０１２１】
この場合、具体的には、計測データ値の平均値の変化量（微分）を算出し、変化量が予め定めた閾値未満となる時刻までを、切出時間幅ＴＷ（ｍ）とする。なお、閾値は、０に近い値とする。
【０１２２】
例えば、計測データを、１０００回取得した場合（１０００度数）の、計測データ値の平均値の推移を図１５に示す。本グラフにおいて、横軸は、計測データの度数、縦軸は、計測データの平均値である。
【０１２３】
本図に示すように、１から５０度数の間、３５１〜４００度数の間、４５１〜５００の間、６５１〜７００度数の間、８００〜８５１度数の間等に、平均値の変化は、極値を取る。これらの各期間のうち、１から５０度数の間の平均値の推移、５１〜４００度数の間の平均値の推移、３５１〜５００度数の間の平均値の推移、４８１〜７００度数の間の平均値の推移の拡大図を、それぞれ、図１６（ａ）〜図１７（ｂ）に示す。各図において、縦軸は、計測データの平均値であり、横軸は、グラフの観測期間における観測数の度数（累計）である。すなわち、横軸の最左値（０）は、図１６（ａ）においては、１、図１６（ｂ）においては、５１、図１７（ａ）においては、３５１、図１７（ｂ）においては、４８１である。
【０１２４】
本変形例では、計測データを得る毎に、時刻ｔ０からその時点までの各計測データの平均値を算出し平均値の推移とする。１回前に算出した平均値と差を算出する。その差がこれらの図に示すように、略０に近い場合、極値を取るものとして、時間区切りとする。
【０１２５】
＜＜第二の実施形態＞＞
次に、本発明の第二の実施形態を説明する。本実施形態では、同時に取得される複数の観測センサ３００からの計測データそれぞれを予測する。ただし、切出時間幅は、特定の計測データを用いて決定し、その切出時間幅を、他の全ての計測データに適用する。
【０１２６】
以後、切出時間幅を決定する計測データを主計測データ（第一の観測センサで得た第一のセンサデータ）、その他のデータを、従計測データ（第二の観測センサで得た第二のセンサデータを含む）と呼び、本実施形態を説明する。例えば、第一の実施形態の図１０（ａ）の例では、３２０が主計測データであり、その他の３２１および３２２が従計測データである。
【０１２７】
なお、本実施形態では、複数の計測データ全てについて、それぞれ、当該計測データを主計測データとし、残りの計測データを従計測データとして、以下の処理を行う。以下、特定の計測データを主計測データとし、残りの計測データを従計測データとする場合を例に、本実施形態を説明する。
【０１２８】
本実施形態の状態予測装置１００のハードウェア構成および機能ブロックは、第一の実施形態と同様である。以下、本実施形態について、第一の実施形態と異なる構成に主眼をおいて説明する。
【０１２９】
ただし、本実施形態では、主計測データ３２０以外に、従計測データ３２１、３２２についても、主計測データ３２０の切出時間幅で切出し、学習データを生成する。そして、それぞれ、学習済みモデルを生成し、各カテゴリ（ラベル）毎の代表データおよび時系列モデルを生成する。
【０１３０】
本実施形態の概要を、図１８を用いて説明する。各計測データａ〜ｅについて、それぞれ、第一の実施形態の手法で切出時間幅を決定し、それぞれ、ラベル２４４を付した状態を図１８に示す。
【０１３１】
ここで、計測データａを主計測データとし、残りの計測データｂ〜ｅを従計測データとする。本実施形態では、例えば、主計測データを計測データａとした場合、主計測データａの切出時間幅で、従計測データｂ〜ｅの計測データも切り出し、それぞれ、学習データを生成する。
【０１３２】
例えば、主計測データａのラベルａ２が付された切出時間幅について、従計測データｂ〜ｅは、その計測データ固有の時間区切り（切出時間幅）に関わらず、この切出時間幅で切出され、その後の処理に用いられる。
【０１３３】
本実施形態の学習部２１４による学習の様子を、図１９を用いて説明する。本実施形態では、主計測データから得た学習データと、従計測データから得た学習データとを、全て、データ入力ウィンドウに入力し、第一の実施形態同様、学習を行う。そして、最終的に学習済みモデル２５２を得る。
【０１３４】
図２０に、本実施形態の学習処理の処理フローを示す。本実施形態の学習処理については、ステップＳ１１０３の、切出時間幅の決定までは、第一の実施形態と同様である。ただし、切出時間幅の決定に用いるデータは、主計測データ３２０である。
【０１３５】
切出時間幅を決定した後、学習データ生成部２１２は、主計測データ３２０および従計測データ３２１、３２２それぞれの、学習データを生成する（ステップＳ２１０１）。このとき、切出時間幅ＴＷ（ｍ）には、第一の実施形態同様の手法で、主計測データ３２０で決定したものを用いる。その他の生成手法は、第一の実施形態の学習データ生成手法と同様である。本実施形態においても、計測データそれぞれについて、波形データおよびヒストグラムの学習データを生成してもよい。
【０１３６】
その後、所定量の学習データが蓄積されたら(ステップＳ１１０５)、クラスタリング部２１３は、主計測データ３２０から生成した学習データのクラスタリングと、従計測データ３２１、３２２のクラスタリングと、を行う（ステップＳ２１０２）。主計測データのクラスタリングは、第一の実施形態と同様である。
【０１３７】
一方、従計測データについては、従計測データを、切出時間幅ＴＷ（ｍ）で切り取った学習データを、主計測データ同様の手法で、それぞれ、分類する。これにより、各計測データについて、教師有データを生成する。なお、カテゴリ数、各カテゴリのラベルは、従計測データ種毎に、予め用意する。
【０１３８】
そして、学習部２１４は、図１９に示すように、生成された学習データセット２４３ａ全てを学習アルゴリズム２５１に入力して学習させ(ステップＳ２１０３)、学習済みモデル２５２を生成するとともに、代表データを算出する。
【０１３９】
以後、本実施形態においても、時系列モデル生成部２１５が、モデル生成単位期間ＴＤが経過したか否かを判別する（ステップＳ１１０６）。そして、未経過の場合、切出時間幅を算出するカウンタｍの値を１インクリメントし（ステップＳ１１０８）、ステップＳ１１０１へ戻り、処理を繰り返す。
【０１４０】
一方、モデル生成単位期間ＴＤが経過した場合、時系列モデル生成部２１５は、主計測データおよび従計測データそれぞれの波形データおよびヒストグラムについて、時系列モデル２７３を生成し（ステップＳ１１０７）、処理を終了する。
【０１４１】
また、予測処理においても、学習処理同様、主計測データを用いて切出時間幅を算出する。その後、主計測データおよび従計測データそれぞれについて、波形データおよびヒストグラムを決定した切出時間幅で切り取り、予測用データを生成する。そして、それぞれの予測用データを用いて、分類し、予測を行う。
【０１４２】
本実施形態の予測処理の流れを、図２１を用いて説明する。本実施形態の予測処理も、ステップＳ１２０３の切出時間幅の決定までは、第一の実施形態と同様である。ただし、切出時間幅の決定に用いるデータは、主計測データ３２０である。
【０１４３】
切出時間幅を決定した後、予測用データ生成部２２２は、主計測データ３２０および従計測データ３２１、３２２それぞれの、予測用データを生成する（ステップＳ２２０４）。予測用データの生成手法は、第一の実施形態と同様である。
【０１４４】
その後、分類決定部２２３は、第一の実施形態と同様に、全ての予測用データを学習済みモデルのディープラーニングネットに入力し、各予測用データの分類を決定し、付与するラベル２４４を特定する（ステップＳ２２０５）。ここでは、主計測データ３２０および従計測データ３２１、３２２それぞれから生成した予測用データについて、分類の決定およびラベル２４４の特定を行う。
【０１４５】
その後、第一の実施形態同様、予測値算出部２２４は、主計測データ３２０、従計測データ３２１、３２２それぞれについて、類似時系列モデルを特定する（ステップＳ２２０６）。そして、未来時刻の波形データおよびヒストグラムを予測し（ステップＳ２２０７）、予測値とする。
【０１４６】
そして、マッチング評価部２２５は、以降の観測センサ３００から出力される実際の計測データと、ステップＳ２２０７で得た予測値とを比較し、第一の実施形態と同様の評価を行う(ステップＳ１２０８)。
【０１４７】
上述したように、本実施形態では、全ての種類の計測データについて、それぞれを主計測データ、残りの計測データを従計測データとして、上記処理を行う。その結果に応じ、以降、主計測データを１つ、または、特定のグループの計測データに固定してもよい。主計測データは、複数の種類の計測データの中で、例えば、計測データ値の変化が、最も動作モードの変化に追従するデータを選択することが望ましい。
【０１４８】
以上説明したように、本実施形態の状態予測装置１００は、第一の実施形態の構成に加え、さらに、従計測データについても、主計測データと同様にクラスタリングし、学習済みモデル２５２、代表データ２６３および時系列モデル２７３を生成する。
【０１４９】
このように、本実施形態によれば、第一の実施形態の効果に加え、簡易な構成で多種多様な計測データを、高精度に予測をすることができる。
【０１５０】
本発明は上記の実施形態に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施形態は、本発明を分かりやすく説明するためのものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施形態の構成の一部を他の実施形態の構成に置き換えることが可能であり、また、ある実施形態の構成に他の実施形態の構成を加えることも可能である。また、各実施形態の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
【０１５１】
また、上記の各構成、機能、処理部、処理手段等は、それらの一部または全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。
【０１５２】
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。
【符号の説明】
【０１５３】
１００：状態予測装置、１１２：ＣＰＵ、１１３：ＲＯＭ、１１４：ＲＡＭ、１１５：ファイルユニット、１１６：外部インタフェース、１１７：モニタ、１１８：バス、
２１０：前処理学習部、２１１：時間幅算出部、２１２：学習データ生成部、２１３：クラスタリング部、２１４：学習部、２１５：時系列モデル生成部、
２２０：予測部、２２１：時間幅算出部、２２２：予測用データ生成部、２２３：分類決定部、２２４：予測値算出部、２２５：マッチング評価部、
２３０：ラベル記憶部、２３１：ラベルデータ、２４０：学習データセット記憶部、２４１：教師有データ、２４１ｈ：教師有データ、２４１ｗ：教師有データ、２４２：取得時刻、２４２ｈ：取得時刻、２４２ｗ：取得時刻、２４３：学習データ、２４３ａ：学習データセット、２４３ｂ：予測用データ、２４３ｃ：分類、２４３ｗ：波形データの学習データ、２４３ｈ：ヒストグラムの学習データ、２４４：ラベル、２４４ｈ：ラベル、２４４ｗ：ラベル、２５０：学習済みモデル記憶部、２５１：学習アルゴリズム、２５２：学習済みモデル、２６０：代表データ記憶部、２６１：代表データベース、２６１ｗ：代表データベース、２６１ｈ：代表データベース、２６３：代表データ、２６３ｈ：代表ヒストグラム、２６３ｗ：代表波形データ、２７０：時系列モデル記憶部、２７１：時系列モデルデータ、２７２：モデルＩＤ、２７３：時系列モデル、
３００：観測センサ、３１１：電流センサ、３１２：電圧センサ、３１３：温度センサ、３１４：圧力センサ、３１５：加速度センサ、３１６：カメラ、３２０：計測データ（主計測データ）、３２１：従計測データ、３２２：従計測データ、３３０：ヒストグラム
【図１】
2019163141.tif 000003
【図２】
2019163141.tif 000004
【図３】
2019163141.tif 000005
【図４】
2019163141.tif 000006
【図５】
2019163141.tif 000007
【図６】
2019163141.tif 000008
【図７】
2019163141.tif 000009
【図８】
2019163141.tif 000010
【図９】
2019163141.tif 000011
【図１０】
2019163141.tif 000012
【図１１】
2019163141.tif 000013
【図１２】
2019163141.tif 000014
【図１３】
2019163141.tif 000015
【図１４】
2019163141.tif 000016
【図１５】
2019163141.tif 000017
【図１６】
2019163141.tif 000018
【図１７】
2019163141.tif 000019
【図１８】
2019163141.tif 000020
【図１９】
2019163141.tif 000021
【図２０】
2019163141.tif 000022
【図２１】
2019163141.tif 000023
【手続補正書】
【提出日】平成31年1月28日(2019.1.28)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】０００６
【補正方法】変更
【補正の内容】
【０００６】
本発明は請求の範囲に記載の構成を有する。その一態様をあげるならば、本発明に係る状態予測装置は、観測対象の状態を観測する第一の観測センサで得た第一のセンサデータを、所定の時間間隔で取得し、前記観測対象の動作モード毎に区切り、学習データを生成する学習データ生成部と、生成した所定数の前記学習データそれぞれを予め定めたカテゴリに分類し、前記カテゴリ毎に予め用意されたラベルを付し、教師有データを生成するクラスタリング部と、前記所定数の教師有データを、学習アルゴリズムに入力して学習し、当該学習アルゴリズムの学習済みモデルを生成するとともに、前記カテゴリそれぞれの代表データを生成する学習部と、予め定めたモデル生成単位期間内に前記クラスタリング部で前記学習データに付された前記ラベルを、当該学習データの生成順に時系列に並べ、時系列モデルを生成する時系列モデル生成部と、新たに取得した第一のセンサデータから、前記時系列モデルを用いて未来時刻の前記観測対象の状態として前記第一のセンサデータの予測値を算出する予測部と、を備える。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００７８
【補正方法】変更
【補正の内容】
【００７８】
ここで、時系列モデル生成部２１５は、予め定めたモデル生成単位期間ＴＤ分の計測データのクラスタリングを終えたか否かを判別する（ステップＳ１１０８）。ここでは、処理開始からモデル生成単位期間ＴＤが経過したか否かを判別する。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】０１２４
【補正方法】変更
【補正の内容】
【０１２４】
本変形例では、計測データを得る毎に、時刻ｔ０からその時点までの各計測データの平均値を算出し平均値の推移とする。１回前に算出した平均値との差を算出する。その差がこれらの図に示すように、略０に近い場合、極値を取るものとして、時間区切りとする。
【手続補正４】
【補正対象書類名】図面
【補正対象項目名】図１
【補正方法】変更
【補正の内容】
【図１】
2019163141.tif 000024
【手続補正５】
【補正対象書類名】図面
【補正対象項目名】図２０
【補正方法】変更
【補正の内容】
【図２０】
2019163141.tif 000025
【国際調査報告】
2019163141.tif 000026
2019163141.tif 000027
2019163141.tif 000028
2019163141.tif 000029