(57)【要約】【課題】状態量の値の確率分布を考慮して、各々のモデルにより算出された状態量の推定値を基に、対象装置の管理に用いる状態量の値を適切に設定する。【解決手段】推定部は、対象装置を説明する複数のモデルを用いて、目的状態量に係る複数の推定値をそれぞれ推定する。確率特定部は、目的状態量に関する状態量の値の確率分布に基づいて、複数の推定値に対応する複数の確率をそれぞれ特定する。管理値特定部は、複数の推定値と複数の確率とに基づいて、対象装置の管理に用いる値を特定する。【選択図】図１
2018169748.tif 000002
【特許請求の範囲】
【請求項１】
対象装置を説明する複数のモデルを用いて、目的状態量に係る複数の推定値をそれぞれ推定する推定部と、
前記目的状態量に関する状態量の値の確率分布に基づいて、前記複数の推定値に対応する複数の確率をそれぞれ特定する確率特定部と、
前記複数の推定値と前記複数の確率とに基づいて、前記対象装置の管理に用いる値を特定する管理値特定部と
を備える情報処理装置。
【請求項２】
前記複数の推定値に基づいて、前記対象装置の管理の評価項目の値である複数の評価値をそれぞれ算出する評価値算出部をさらに備え、
前記確率特定部は、前記評価項目の値の確率分布に基づいて、前記複数の評価値それぞれに対応する前記複数の確率を特定する
請求項１に記載の情報処理装置。
【請求項３】
前記評価値算出部は、値が不明な状態量である不明状態量がとり得る複数の値に基づいて、前記複数の推定値別に複数の評価値をそれぞれ算出し、
前記確率特定部は、前記不明状態量の値を前提条件とした前記評価項目の値の条件付き確率分布に基づいて、前記複数の推定値別の複数の評価値それぞれに対応する前記複数の確率を特定し、
前記管理値特定部は、前記複数の推定値別の複数の評価値それぞれに対応する確率の総和に基づいて、前記対象装置の管理に用いる値を特定する
請求項２に記載の情報処理装置。
【請求項４】
前記評価値算出部は、前記複数の推定値に基づいて、複数種類の前記評価項目それぞれに係る複数の評価値を算出し、
前記管理値特定部は、前記複数の目的状態量の値それぞれについての、前記複数種類の評価項目に対応する各確率に基づいて、前記対象装置の管理に用いる値を特定する
請求項２または請求項３に記載の情報処理装置。
【請求項５】
前記管理値特定部は、最も高い前記確率に関する前記目的状態量の値を前記対象装置の管理に用いる値として特定する
請求項１から請求項４の何れか１項に記載の情報処理装置。
【請求項６】
前記複数のモデルは、少なくとも統計モデルおよび物理モデルのいずれか一方を含む
請求項１から請求項５の何れか１項に記載の情報処理装置。
【請求項７】
過去の前記状態量の値に基づいて前記統計モデルを更新するモデル更新部をさらに備え、
前記確率特定部は、更新された前記統計モデルを用いて推定された推定値に対応する確率を特定する
請求項６に記載の情報処理装置。
【請求項８】
対象装置を説明する複数のモデルを用いて、目的状態量に係る複数の推定値をそれぞれ推定することと、
前記目的状態量に関する値の確率分布に基づいて、前記複数の推定値に対応する複数の確率をそれぞれ特定することと、
前記複数の推定値と前記複数の確率とに基づいて、前記対象装置の管理に用いる値を特定することと
を有する情報処理方法。
【請求項９】
コンピュータに、
対象装置を説明する複数のモデルを用いて、目的状態量に係る複数の推定値をそれぞれ推定することと、
前記目的状態量に関する値の確率分布に基づいて、前記複数の推定値に対応する複数の確率をそれぞれ特定することと、
前記複数の推定値と前記複数の確率とに基づいて、前記対象装置の管理に用いる値を特定することと
を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理装置、情報処理方法およびプログラムに関する。
【背景技術】
【０００２】
発電プラント等のプラントにおいて、プラントを構成する対象装置の運用時の温度および圧力等の対象装置の状態量を監視装置が収集し、収集した状態量を用いて装置の保守および監視等に用いることが検討されている。
監視装置が、装置の運用者が使用し易いように収集した状態量を加工して、装置の保守および監視等を実行することが提案されている。例えば、特許文献１には、コンピュータが、収集された産業用プラントの状態量の欠落部分を検出すると、補完処理を実行し、欠落した状態量の値を算出することが提案されている。
また、補完処理により、欠落した状態量の値の設定を行う際には、物理モデルおよび統計モデル等のモデルが用いられることが知られている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】米国特許出願公開第２０１６／０００４７９４号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１には、情報処理装置が補完処理を実行することにより、欠落した状態量の値を補完することは開示されているものの、欠落した状態量の値の具体的な設定方法については開示されていない。
また、装置の状態量が取る値は、一定の確率分布に従うが、このような確率分布を考慮して、各々のモデルにより推定された状態量の推定値を基に、欠落した状態量の値を設定することは行われていなかった。
本発明は、上記の課題に鑑みてされたものであって、状態量の値の確率分布を考慮して、各々のモデルにより算出された状態量の推定値を基に、対象装置の管理に用いる状態量の値を適切に設定することを目的とする。
【課題を解決するための手段】
【０００５】
本発明の第１の態様によれば、情報処理装置は、対象装置を説明する複数のモデルを用いて、目的状態量に係る複数の推定値をそれぞれ推定する推定部と、前記目的状態量に関する状態量の値の確率分布に基づいて、前記複数の推定値に対応する複数の確率をそれぞれ特定する確率特定部と、前記複数の推定値と前記複数の確率とに基づいて、前記対象装置の管理に用いる値を特定する管理値特定部とを備える。
【０００６】
本発明の第２の態様によれば、第１の態様に係る情報処理装置は、前記複数の推定値に基づいて、前記対象装置の管理の評価項目の値である複数の評価値をそれぞれ算出する評価値算出部をさらに備え、前記確率特定部は、前記評価項目の値の確率分布に基づいて、前記複数の評価値それぞれに対応する前記複数の確率を特定するものであってよい。
【０００７】
本発明の第３の態様によれば、第２の態様に係る情報処理装置は、前記評価値算出部は、値が不明な状態量である不明状態量がとり得る複数の値に基づいて、前記複数の推定値別に複数の評価値をそれぞれ算出し、前記確率特定部は、前記不明状態量の値を前提条件とした前記評価項目の値の条件付き確率分布に基づいて、前記複数の推定値別の複数の評価値それぞれに対応する前記複数の確率を特定し、前記管理値特定部は、前記複数の推定値別の複数の評価値それぞれに対応する確率の総和に基づいて、前記対象装置の管理に用いる値を特定するものであってよい。
【０００８】
本発明の第４の態様によれば、第２または第３の態様に係る情報処理装置は、前記評価値算出部は、前記複数の推定値に基づいて、複数種類の前記評価項目それぞれに係る複数の評価値を算出し、前記管理値特定部は、前記複数の目的状態量の値それぞれについての、前記複数種類の評価項目に対応する各確率に基づいて、前記対象装置の管理に用いる値を特定するものであってよい。
【０００９】
本発明の第５の態様によれば、第１から第４の何れかの態様に係る情報処理装置は、前記管理値特定部は、最も高い前記確率に関する前記目的状態量の値を前記対象装置の管理に用いる値として特定するものであってよい。
【００１０】
本発明の第６の態様によれば、第１から第５の何れかの態様に係る情報処理装置は、前記複数のモデルは、少なくとも統計モデルおよび物理モデルのいずれか一方を含むものであってよい。
【００１１】
本発明の第７の態様によれば、第６の態様に係る情報処理装置は、過去の前記状態量の値に基づいて前記統計モデルを更新するモデル更新部をさらに備え、前記確率特定部は、更新された前記統計モデルを用いて推定された推定値に対応する確率を特定するものであってよい。
【００１２】
本発明の第８の態様によれば、情報処理方法は、対象装置を説明する複数のモデルを用いて、目的状態量に係る複数の推定値をそれぞれ推定することと、前記目的状態量に関する値の確率分布に基づいて、前記複数の推定値に対応する複数の確率をそれぞれ特定することと、前記複数の推定値と前記複数の確率とに基づいて、前記対象装置の管理に用いる値を特定することとを有する。
【００１３】
本発明の第９の態様によれば、プログラムは、コンピュータに、対象装置を説明する複数のモデルを用いて、目的状態量に係る複数の推定値をそれぞれ推定することと、前記目的状態量に関する値の確率分布に基づいて、前記複数の推定値に対応する複数の確率をそれぞれ特定することと、前記複数の推定値と前記複数の確率とに基づいて、前記対象装置の管理に用いる値を特定することとを実行させる。
【発明の効果】
【００１４】
上記態様のうち少なくとも１つの態様に係る情報処理装置は、状態量の値の確率分布を考慮して、各モデルにより算出された目的状態量の推定値を基に、対象装置の管理に用いる状態量の値を適切に設定することができる。
【図面の簡単な説明】
【００１５】
【図１】第１の実施形態に係る管理システムの構成を示す概略ブロック図である。
【図２】第１の実施形態に係る管理装置の構成を示す概略ブロック図である。
【図３】第１の実施形態に係る管理装置の動作を示すフローチャートである。
【図４】第１の実施形態に係る管理値の特定方法の具体例を示す図である。
【図５】第２の実施形態に係る管理装置の構成を示す概略ブロック図である。
【図６】第２の実施形態に係る管理装置の動作を示すフローチャートである。
【図７】第２の実施形態に係る管理値の特定方法の具体例を示す図である。
【図８】第３の実施形態に係る管理装置の動作を示すフローチャートである。
【図９】第３の実施形態に係る管理値の特定方法の具体例を示す図である。
【図１０】第４の実施形態に係る管理装置の動作を示すフローチャートである。
【図１１】第４の実施形態に係る管理値の特定方法の具体例を示す図である。
【図１２】第５の実施形態に係る管理装置の構成を示す概略ブロック図である。
【図１３】第５の実施形態に係る管理装置の動作を示すフローチャートである。
【図１４】少なくとも１つの実施形態に係るコンピュータの構成を示す概略ブロック図である。
【発明を実施するための形態】
【００１６】
〈第１の実施形態〉
以下、図面を参照しながら実施形態について詳しく説明する。
図１は、第１の実施形態に係る管理システムの構成を示す概略ブロック図である。
管理システム１は、対象装置１０と、複数の計測器２０と、通信装置３０と、管理装置４０とを備える。
対象装置１０は、管理装置４０による管理対象の装置である。対象装置１０の例としては、例えばガスタービン、蒸気タービン、ボイラ、石炭ガス化炉などが挙げられる。また、環境プラント、化学プラント、および航空機のような交輸システムであってもよい。
計測器２０は、対象装置１０に設けられ、対象装置１０の状態量を計測する。
通信装置３０は、計測器２０が計測した状態量の計測値をネットワークＮを介して管理装置４０に送信する。
管理装置４０は、通信装置３０から受信した計測値に基づいて対象装置１０を管理する。管理装置４０は、情報処理装置の一例である。
【００１７】
《管理装置の構成》
図２は、第１の実施形態に係る管理装置の構成を示す概略ブロック図である。
管理装置４０は、計測値取得部４０１と、欠落検出部４０２と、モデル記憶部４０３と、推定部４０４と、確率分布記憶部４０５と、確率特定部４０６と、管理値特定部４０７と、管理部４０８ととを備える。
【００１８】
計測値取得部４０１は、通信装置３０から、複数の計測器２０によって計測された状態量の計測値を受信する。
欠落検出部４０２は、計測値取得部４０１が取得した複数の計測値に基づいて、管理すべき状態量のうち値が欠落しているものを検出する。ここで、値の欠落とは、時間的または空間的な欠落をいう。例えば、管理部４０８が時間Δｔごとの状態量を管理する場合において、時刻Ｔの計測値と時刻Ｔ＋２Δｔの計測値が取得された場合、欠落検出部４０２は、時刻Ｔ＋Δｔの計測値の欠落を検出する。また例えば、管理部４０８が距離Δｄごとの状態量を管理する場合において、位置（０，０）、位置（２Δｄ，０）、位置（０，２Δｄ）、位置（２Δｄ，２Δｄ）の計測値が取得された場合、位置（０，Δｄ）、位置（Δｄ，０）、位置（Δｄ，Δｄ）、位置（Δｄ，２Δｄ）、位置（２Δｄ，Δｄ）の計測値の欠落を検出する。
【００１９】
モデル記憶部４０３は、対象装置１０の挙動を説明する複数のモデルを記憶する。モデルとしては、統計モデルや物理モデルを用いることができる。このほか、ルールモデルや知識モデルを用いても良い。統計モデルとは、対象装置１０の過去の運用における状態量の値に基づいて、統計学的に対象装置１０の挙動を再現したモデルである。統計モデルは、蓄積された過去の運用における状態量の値に基づいて更新される。物理モデルとは、対象装置１０の設計情報に基づいて、自然法則に倣った数式（例えば、熱力学方程式）により対象装置１０の挙動を再現したモデルである。
推定部４０４は、計測値取得部４０１が取得した計測値に基づいて、モデル記憶部４０３が記憶する各モデルごとに状態量の値を推定する。以下、推定部４０４による推定の対象となる状態量を目的状態量という。つまり、推定部４０４は、異なるモデルを用いて複数の目的状態量の値を算出する。
【００２０】
確率分布記憶部４０５は、目的状態量の値とその値をとる確率とを関連付けた確率分布テーブルを記憶する。確率分布テーブルは、対象装置１０の設計情報または過去の状態量の統計によって予め求められたものである。なお、確率分布記憶部４０５は、確率分布テーブルに代えて確率分布関数を記憶してもよい。
確率特定部４０６は、確率分布記憶部４０５が記憶する確率分布に基づいて、推定部４０４による目的状態量の推定値ごとに、その値を取る確率を特定する。
管理値特定部４０７は、確率特定部４０６が特定した確率に基づいて推定部４０４が推定した複数の推定値の中から１つを選択し、対象装置１０の管理に用いる値（管理値）とする。
管理部４０８は、計測値取得部４０１が取得した計測値および管理値特定部４０７が特定した値に基づいて対象装置１０を管理する。対象装置１０の管理の例としては、対象装置１０の状態量が運転許容範囲を逸脱していないか監視すること、対象装置１０の評価項目に係る出力が目標を満たしているかを監視すること、および対象装置１０に制御信号を出力することなどが挙げられる。評価項目の例としては、ＮＯｘ排出量、売電収入額、ガス温度などが挙げられる。
【００２１】
《管理装置の動作》
図３は、第１の実施形態に係る管理装置の動作を示すフローチャートである。
管理装置４０が対象装置１０の管理を開始すると、計測値取得部４０１は、通信装置３０から計測器２０による状態量の計測値を取得する（ステップＳ１）。次に、欠落検出部４０２は、計測値取得部４０１が取得した計測値の欠落を検出する（ステップＳ２）。推定部４０４は、計測値取得部４０１が取得した計測値を複数のモデルのそれぞれに適用し、欠落が検出された状態量（目的状態量）の推定値をそれぞれ求める（ステップＳ３）。
【００２２】
次に、確率特定部４０６は、確率分布記憶部４０５が記憶する確率分布に基づいて、各推定値の出現確率を特定する（ステップＳ４）。そして、管理値特定部４０７は、確率特定部４０６が特定した確率のうち最も高いものを特定し、当該確率に係る推定値を選択することで、欠落が検出された状態量の値を特定する（ステップＳ５）。そして管理部４０８は、計測値取得部４０１が取得した計測値および管理値特定部４０７が特定した値に基づいて対象装置１０を管理する（ステップＳ６）。対象装置１０がガスタービンの場合、例えば、ガスタービン出力指令値を変更する、ＩＧＶの開度設定を変更する、燃料流量を変更する、等により、特定された管理値に基づいて、対象装置１０を管理する。
【００２３】
《動作の具体例》
ここで、具体的な例を用いて第１の実施形態に係る管理値の特定方法について説明する。
図４は、第１の実施形態に係る管理値の特定方法の具体例を示す図である。
目的状態量の値の確率分布が図４に示す分布であり、推定部４０４が第１のモデルに基づく推定値ｅ１と、第２のモデルに基づく推定値ｅ２とを出力した場合について説明する。図４に含まれるグラフＧ１は、縦軸に確率密度をとり、横軸に目的状態量の値をとるグラフである。確率特定部４０６は、目的状態量の確率分布に基づいて、推定値ｅ１の出現確率を求める。図４に示す例においては、推定値ｅ１の出現確率の確率密度は０．２である。また確率特定部４０６は、目的状態量の確率分布に基づいて、推定値ｅ２の出現確率を求める。図４に示す例においては、推定値ｅ２の出現確率の確率密度は０．３である。そして、管理値特定部４０７は、特定した各推定値の出現確率のうち、大きい方に係る推定値を、管理値と特定する。図４に示す例においては、推定値ｅ１の出現確率より推定値ｅ２の出現確率の方が大きいため、管理値特定部４０７は、管理値を推定値ｅ２に決定する。
【００２４】
《作用・効果》
このように、第１の実施形態によれば、管理装置４０は、目的状態量の値の確率分布に基づいて、複数の推定値から対象装置１０の管理に用いる値を特定する。つまり、第１の実施形態によれば、管理装置４０は、目的状態量の値の確率分布を考慮して、各々のモデルにより算出された目的状態量の推定値を基に、対象装置１０の管理に用いる目的状態量の値を適切に設定することができる。
【００２５】
〈第２実施形態〉
第１の実施形態に係る管理装置４０は、目標状態量の値の確率分布に基づいて対象装置１０の管理に用いる目的状態量の値を設定する。これに対し、第２の実施形態に係る管理装置４０は、対象装置１０の管理の評価項目の値の確率分布に基づいて、対象装置１０の管理に用いる目的状態量の値を設定する。評価項目の例としては、ＮＯｘ排出量、売電収入額、ガス温度などが挙げられる。
【００２６】
《管理装置の構成》
図５は、第２の実施形態に係る管理装置の構成を示す概略ブロック図である。
第２の実施形態に係る管理装置４０は、第１の実施形態の構成に加えて、評価値算出部４０９をさらに備える。評価値算出部４０９は、推定部４０４が推定した複数の推定値のそれぞれについて、当該推定値と計測値取得部４０１が計測した計測値とを用いて、対象装置１０の評価項目の値を算出する。
【００２７】
第２の実施形態に係る確率分布記憶部４０５は、評価項目の値とその値をとる確率とを関連付けた確率分布テーブルを記憶する。
第２の実施形態に係る確率特定部４０６は、確率分布記憶部４０５が記憶する確率分布に基づいて、評価値算出部４０９が算出した評価値ごとに、その値を取る確率を特定する。
【００２８】
《管理装置の動作》
図６は、第２の実施形態に係る管理装置の動作を示すフローチャートである。
管理装置４０が対象装置１０の管理を開始すると、計測値取得部４０１は、通信装置３０から計測器２０による状態量の計測値を取得する（ステップＳ１０１）。次に、欠落検出部４０２は、計測値取得部４０１が取得した計測値の欠落を検出する（ステップＳ１０２）。推定部４０４は、計測値取得部４０１が取得した計測値を複数のモデルのそれぞれに適用し、欠落が検出された状態量の推定値をそれぞれ求める（ステップＳ１０３）。
【００２９】
次に、評価値算出部４０９は、推定部４０４が推定した複数の推定値のそれぞれに基づいて、複数の評価項目の値を算出する（ステップＳ１０４）。評価項目の値は、複数の状態量の値を説明変数とする関数によって求めることができる。評価値算出部４０９は、当該関数の説明変数に、計測値または推定値を代入することで、評価値を算出する。
【００３０】
次に、確率特定部４０６は、確率分布記憶部４０５が記憶する確率分布に基づいて、各評価値の出現確率を特定する（ステップＳ１０５）。そして、管理値特定部４０７は、確率特定部４０６が特定した確率のうち最も高いものを特定し、当該確率に係る評価値の算出に用いた推定値を選択することで、欠落が検出された状態量の値を特定する（ステップＳ１０６）。そして管理部４０８は、計測値取得部４０１が取得した計測値および管理値特定部４０７が特定した値に基づいて対象装置１０を管理する（ステップＳ１０７）。
【００３１】
《動作の具体例》
ここで、具体的な例を用いて第２の実施形態に係る管理値の特定方法について説明する。
図７は、第２の実施形態に係る管理値の特定方法の具体例を示す図である。
評価対象の値の確率分布が図７に示す分布であり、推定部４０４が第１のモデルに基づく推定値ｅ１と、第２のモデルに基づく推定値ｅ２とを出力した場合について説明する。図６に含まれるグラフＧ２は、縦軸に確率密度をとり、横軸に確率分布の値をとるグラフである。評価値算出部４０９は、推定値ｅ１に基づいて評価項目の値である評価値ｆ（ｅ１）を算出する。また、評価値算出部４０９は、推定値ｅ２に基づいて評価項目の値である評価値ｆ（ｅ２）を算出する。次に、確率特定部４０６は、評価項目の値の確率分布に基づいて、評価値ｆ（ｅ１）の出現確率を求める。図７に示す例においては、評価値ｆ（ｅ１）の出現確率の確率密度は０．３０である。また確率特定部４０６は、評価項目の値の確率分布に基づいて、評価値ｆ（ｅ２）の出現確率を求める。図７に示す例においては、評価値ｆ（ｅ２）の出現確率の確率密度は０．３５である。そして、管理値特定部４０７は、特定した各評価値の出現確率のうち、大きい方の算出に用いられた推定値を、管理値と特定する。図７に示す例においては、評価値ｆ（ｅ１）の出現確率より評価値ｆ（ｅ２）の出現確率の方が大きいため、管理値特定部４０７は、評価値ｆ（ｅ２）の算出に用いた推定値ｅ２を、管理値とする。
【００３２】
《作用・効果》
このように、第２の実施形態によれば、管理装置４０は、目的状態量に基づいて算出される評価対象の値の確率分布に基づいて、複数の推定値から対象装置１０の管理に用いる値を特定する。つまり、第２の実施形態によれば、管理装置４０は、評価対象の値の確率分布を考慮して、各々のモデルにより算出された目的状態量の推定値を基に、対象装置１０の管理に用いる目的状態量の値を適切に設定することができる。
【００３３】
〈第３の実施形態〉
第２の実施形態に係る管理装置４０は、対象装置１０の管理の評価項目の値の確率分布に基づいて、対象装置１０の管理に用いる目的状態量の値を設定する。ここで、評価項目によっては、計測することも予測することもできない状態量によってその値が変動するものがある。例えば、評価項目の１つであるＮＯｘ排出量は、燃焼時の酸素濃度や高音域燃焼ガス滞在時間によって変動するが、これらの値は計測することも予測することもできない。以下、計測することも予測することもできない状態量を不明状態量とよぶ。
第３の実施形態に係る管理装置４０は、不明状態量に鑑みて対象装置１０の管理に用いる目的状態量の値を設定する。なお、管理装置４０の構成は、第２の実施形態と同様である。ただし、第３の実施形態に係る確率分布記憶部４０５は、不明状態量の値ごとに、評価項目の値とその値をとる確率とを関連付けた確率分布テーブルを記憶する。つまり、第３の実施形態に係る確率分布テーブルは、不明状態量の値を前提条件とする条件付き確率分布を示すテーブルである。第３の実施形態においては、確率分布記憶部４０５は、不明状態量の値が第１の範囲内である（値が比較的大きい）場合における確率分布テーブルと、不明状態量の値が第２の範囲内である（値が中程度である）場合における確率分布テーブルと、不明状態量の値が第３の範囲内である（値が比較的小さい）場合における確率分布テーブルとを記憶する。
【００３４】
《管理装置の動作》
図８は、第３の実施形態に係る管理装置の動作を示すフローチャートである。
管理装置４０が対象装置１０の管理を開始すると、計測値取得部４０１は、通信装置３０から計測器２０による状態量の計測値を取得する（ステップＳ２０１）。次に、欠落検出部４０２は、計測値取得部４０１が取得した計測値の欠落を検出する（ステップＳ２０２）。推定部４０４は、計測値取得部４０１が取得した計測値を複数のモデルのそれぞれに適用し、欠落が検出された状態量の推定値をそれぞれ求める（ステップＳ２０３）。
【００３５】
次に、管理装置４０は、不明状態量の値（第１の範囲、第２の範囲、第３の範囲）を１つずつ選択し（ステップＳ２０４）、ステップＳ２０５およびステップＳ２０６の処理を実行する。すなわち、評価値算出部４０９は、推定部４０４が推定した複数の推定値のそれぞれについて、当該推定値、ステップＳ２０４で選択された不明状態量の値、およびステップＳ１０１で取得された計測値に基づき、複数の評価項目の値を算出する（ステップＳ２０５）。確率特定部４０６は、ステップＳ２０４で選択された不明状態量の値に関連付けられた確率分布テーブルに基づいて、各評価値の出現確率を特定する（ステップＳ２０６）。
【００３６】
管理装置４０が不明状態量の値ごとに複数の出現確率を算出すると、管理値特定部４０７は、同一の状態量から算出された評価値ごとに、出現確率の総和を算出する（ステップＳ２０７）。つまり、管理値特定部４０７は、同一の状態量から算出された評価値ごとに、不明状態量の値が第１の範囲にあることを前提条件とする評価値の出現確率と、不明状態量の値が第２の範囲にあることを前提条件とする評価値の出現確率と、不明状態量の値が第３の範囲にあることを前提条件とする評価値の出現確率との和を算出する。
そして、管理値特定部４０７は、確率特定部４０６が算出した確率の総和のうち最も高いものを特定し、当該確率に係る評価値の算出に用いた推定値を選択することで、欠落が検出された状態量の値を特定する（ステップＳ２０８）。そして管理部４０８は、計測値取得部４０１が取得した計測値および管理値特定部４０７が特定した値に基づいて対象装置１０を管理する（ステップＳ２０９）。
【００３７】
《動作の具体例》
ここで、具体的な例を用いて第３の実施形態に係る管理値の特定方法について説明する。
図９は、第３の実施形態に係る管理値の特定方法の具体例を示す図である。
不明状態量の値によって評価対象の値の確率分布が図９に示すように変化し、推定部４０４が第１のモデルに基づく推定値ｅ１と、第２のモデルに基づく推定値ｅ２とを出力した場合について説明する。
図９に含まれるグラフＧ２−１、グラフＧ２−２、グラフＧ２−３は、いずれも縦軸に確率密度をとり、横軸に確率分布の値をとるグラフである。グラフＧ２−１は、不明状態量の値が第１の範囲にあるときにおける評価対象の値の出現確率の分布を表す。グラフＧ２−２は、不明状態量の値が第２の範囲にあるときにおける評価対象の値の出現確率の分布を表す。グラフＧ２−３は、不明状態量の値が第３の範囲にあるときにおける評価対象の値の出現確率の分布を表す。
【００３８】
管理装置４０は、ステップＳ２０４において、不明状態量の値として第１の範囲の値を選択する。評価値算出部４０９は、推定値ｅ１に基づいて不明状態量の値が第１の範囲にある場合における評価項目の値である評価値ｆ１（ｅ１）を算出する。また、評価値算出部４０９は、推定値ｅ２に基づいて不明状態量の値が第１の範囲にある場合における評価項目の値である評価値ｆ１（ｅ２）を算出する。次に、確率特定部４０６は、グラフＧ２−１に示す確率分布に基づいて、評価値ｆ１（ｅ１）の出現確率を求める。図９に示す例においては、評価値ｆ１（ｅ１）の出現確率の確率密度は０．３０である。また確率特定部４０６は、グラフＧ２−１に示す確率分布に基づいて、評価値ｆ１（ｅ２）の出現確率を求める。図９に示す例においては、評価値ｆ１（ｅ２）の出現確率の確率密度は０．３５である。
【００３９】
次に、管理装置４０は、不明状態量の値として第２の範囲の値を選択する。評価値算出部４０９は、推定値ｅ１に基づいて不明状態量の値が第２の範囲にある場合における評価項目の値である評価値ｆ２（ｅ１）を算出する。また、評価値算出部４０９は、推定値ｅ２に基づいて不明状態量の値が第２の範囲にある場合における評価項目の値である評価値ｆ２（ｅ２）を算出する。次に、確率特定部４０６は、グラフＧ２−２に示す確率分布に基づいて、評価値ｆ２（ｅ１）の出現確率を求める。図９に示す例においては、評価値ｆ２（ｅ１）の出現確率の確率密度は０．５０である。また確率特定部４０６は、グラフＧ２−２に示す確率分布に基づいて、評価値ｆ２（ｅ２）の出現確率を求める。図９に示す例においては、評価値ｆ２（ｅ２）の出現確率の確率密度は０．１０である。
【００４０】
次に、管理装置４０は、不明状態量の値として第３の範囲の値を選択する。評価値算出部４０９は、推定値ｅ１に基づいて不明状態量の値が第３の範囲にある場合における評価項目の値である評価値ｆ３（ｅ１）を算出する。また、評価値算出部４０９は、推定値ｅ２に基づいて不明状態量の値が第３の範囲にある場合における評価項目の値である評価値ｆ３（ｅ２）を算出する。次に、確率特定部４０６は、グラフＧ２−３に示す確率分布に基づいて、評価値ｆ３（ｅ１）の出現確率を求める。図９に示す例においては、評価値ｆ３（ｅ１）の出現確率の確率密度は０．１８である。また確率特定部４０６は、グラフＧ２−３に示す確率分布に基づいて、評価値ｆ３（ｅ２）の出現確率を求める。図９に示す例においては、評価値ｆ３（ｅ２）の出現確率の確率密度は０．１５である。
【００４１】
そして、管理値特定部４０７は、推定値ｅ１に基づいて算出した評価値に係る出現確率の総和（０．３０＋０．５０＋０．１８＝０．９８）を算出する。また、管理値特定部４０７は、推定値ｅ２に基づいて算出した評価値に係る出現確率の総和（０．３５＋０．１０＋０．１５＝０．６０）を算出する。そして、管理値特定部４０７は、特定した各出現確率の総和のうち、大きい方の算出に用いられた推定値を、管理値と特定する。図９に示す例においては、推定値ｅ２に基づいて算出した評価値の出現確率より推定値ｅ１に基づいて算出した評価値の出現確率の方が大きいため、管理値特定部４０７は、推定値ｅ１を、管理値とする。
【００４２】
《作用・効果》
このように、第３の実施形態によれば、管理装置４０は、不明状態量がとり得る複数の値に基づいて、評価値の出現確率を特定し、推定値別の複数の評価値それぞれに対応する確率の総和に基づいて、対象装置１０の管理に用いる値を特定する。これにより、管理装置４０は、評価対象の値を算出するにあたって計測することも予測することもできない不明な値がある場合にも、対象装置１０の管理に用いる目的状態量の値を適切に設定することができる。
【００４３】
〈第４の実施形態〉
第２、第３の実施形態によれば、管理装置４０は、ある評価項目の値の出現確率に基づいて目的状態量の値を特定する。これに対し、第４の実施形態に係る管理装置４０は、複数の評価項目の値に基づいて目的状態量の値を特定する。例えば、管理装置４０は、ＮＯｘ排出値の値、売電収入額の値、および排出ガス温度の値に基づいて目的状態量の値を特定する。なお、管理装置４０の構成は、第２の実施形態と同様である。
【００４４】
《管理装置の動作》
図１０は、第４の実施形態に係る管理装置の動作を示すフローチャートである。
管理装置４０が対象装置１０の管理を開始すると、計測値取得部４０１は、通信装置３０から計測器２０による状態量の計測値を取得する（ステップＳ３０１）。次に、欠落検出部４０２は、計測値取得部４０１が取得した計測値の欠落を検出する（ステップＳ３０２）。推定部４０４は、計測値取得部４０１が取得した計測値を複数のモデルのそれぞれに適用し、欠落が検出された状態量の推定値をそれぞれ求める（ステップＳ３０３）。
【００４５】
次に、管理装置４０は、評価項目の種類を１つずつ選択し、以下のステップＳ２０５〜ステップＳ２０７の処理を実行する（ステップＳ２０４）。
まず、評価値算出部４０９は、推定部４０４が推定した複数の推定値のそれぞれに基づいて、ステップＳ２０４で選択した種類に係る評価項目の値を算出する（ステップＳ２０５）。次に、確率特定部４０６は、確率分布記憶部４０５が記憶する確率分布に基づいて、各評価値の出現確率を特定する（ステップＳ２０６）。次に、管理値特定部４０７は、特定した出現確率が所定の閾値（例えば確率密度０．３）以上であるか否かを判定する（ステップＳ２０７）。
【００４６】
管理値特定部４０７は、各種類の評価項目について、出現確率が所定の閾値以上であるか否かを判定すると、出現確率が所定の閾値以上であるとされた項目が最も多い推定値を選択することで、欠落が検出された状態量の値を特定する（ステップＳ２０８）。そして管理部４０８は、計測値取得部４０１が取得した計測値および管理値特定部４０７が特定した値に基づいて対象装置１０を管理する（ステップＳ２０９）。
【００４７】
《動作の具体例》
ここで、具体的な例を用いて第４の実施形態に係る管理値の特定方法について説明する。
図１１は、第４の実施形態に係る管理値の特定方法の具体例を示す図である。
推定部４０４が第１のモデルに基づく推定値ｅ１と、第２のモデルに基づく推定値ｅ２とを出力し、算出すべき評価項目の種類が、ＮＯｘ排出量、売電収入額、排出ガス温度である場合について説明する。
評価値算出部４０９は、推定値ｅ１に基づいてＮＯｘ排出量の評価値、売電収入額の評価値、排出ガス温度の評価値をそれぞれ算出する。また評価値算出部４０９は、推定値ｅ２に基づいてＮＯｘ排出量の評価値、売電収入額の評価値、排出ガス温度の評価値をそれぞれ算出する。次に、確率特定部４０６は、推定値ｅ１から求められたＮＯｘ排出量の評価値、売電収入額の評価値、排出ガス温度の評価値のそれぞれについて、出現確率を求める。同様に、確率特定部４０６は、推定値ｅ２から求められたＮＯｘ排出量の評価値、売電収入額の評価値、排出ガス温度の評価値のそれぞれについて、出現確率を求める。
【００４８】
ここで、管理値特定部４０７は、推定値ｅ１から求められたＮＯｘ排出量の評価値の出現確率、推定値ｅ１から求められた売電収入額の評価値の出現確率、推定値ｅ１から求められた排出ガス温度の評価値の出現確率、推定値ｅ２から求められたＮＯｘ排出量の評価値の出現確率、推定値ｅ２から求められた売電収入額の評価値の出現確率、推定値ｅ２から求められた排出ガス温度の評価値の出現確率、のそれぞれについて、所定の閾値以上であるか否かを判定する。ここでは、図１１に示すように、推定値ｅ１から求められたＮＯｘ排出量の評価値の出現確率、推定値ｅ２から求められたＮＯｘ排出量の評価値の出現確率、および推定値ｅ２から求められた排出ガス温度の評価値の出現確率が閾値以上であり、その他は閾値未満であったものとする。
そして、管理値特定部４０７は、出現確率が閾値以上であるとされた項目が最も多い推定値を管理値と特定する。図１１に示す例においては、推定値ｅ１から求めた評価項目のうち出現確率が閾値以上であるとされた項目の数は１つ、推定値ｅ２から求めた評価項目のうち出現確率が閾値以上であるとされた項目の数は２つであるため、管理値特定部４０７は、推定値ｅ２を、管理値とする。
【００４９】
《作用・効果》
このように、第４の実施形態によれば、管理装置４０は、複数の評価項目に基づいて、対象装置１０の管理に用いる値を特定する。これにより、管理装置４０は、対象装置１０の管理に用いる評価項目が妥当な値となるように、対象装置１０の管理に用いる目的状態量の値を適切に設定することができる。
【００５０】
なお、第４の実施形態では、管理装置４０は、出現確率が閾値以上であるとされた項目の数に基づいて管理値を特定するが、これに限られない。例えば、他の実施形態においては、管理装置４０は、評価項目ごとの出現確率の総和や加重平均に基づいて管理値を特定してもよいし、出現確率が最も大きい項目の数に基づいて管理値を特定してもよい。
【００５１】
〈第５の実施形態〉
第１から第４の実施形態によれば、管理装置４０は、複数のモデルに基づいて状態量の推定値を生成する。第５の実施形態では、複数のモデルの１つが統計モデルであるときの動作について説明する。
【００５２】
《管理装置の構成》
図１２は、第５の実施形態に係る管理装置の構成を示す概略ブロック図である。
第５の実施形態に係る管理装置４０は、第１の実施形態の構成に加え、さらに状態量記憶部４１０とモデル更新部４１１とを備える。モデル更新部４１１は、モデル記憶部４０３が記憶する複数のモデルのうち統計モデルを状態量記憶部４１０が記憶する過去の状態量の値に基づいて更新する。
【００５３】
《管理装置の動作》
図１３は、第５の実施形態に係る管理装置の動作を示すフローチャートである。
管理装置４０が対象装置１０の管理を開始すると、計測値取得部４０１は、通信装置３０から計測器２０による状態量の計測値を取得する（ステップＳ３０１）。次に、欠落検出部４０２は、計測値取得部４０１が取得した計測値の欠落を検出する（ステップＳ３０２）。推定部４０４は、計測値取得部４０１が取得した計測値を統計モデルを含む複数のモデルのそれぞれに適用し、欠落が検出された状態量（目的状態量）の推定値をそれぞれ求める（ステップＳ３０３）。
【００５４】
次に、確率特定部４０６は、確率分布記憶部４０５が記憶する確率分布に基づいて、各推定値の出現確率を特定する（ステップＳ３０４）。そして、管理値特定部４０７は、確率特定部４０６が特定した確率のうち最も高いものを特定し、当該確率に係る推定値を選択することで、欠落が検出された状態量の値を特定する（ステップＳ３０５）。管理部４０８は、計測値取得部４０１が取得した計測値および管理値特定部４０７が特定した値に基づいて対象装置１０を管理する（ステップＳ３０６）。計測値取得部４０１および管理値特定部４０７は、対象装置１０の管理に用いた値を、状態量記憶部４１０に蓄積させる（ステップＳ３０７）。そして、モデル更新部４１１は、モデル記憶部４０３が記憶する統計モデルを、状態量記憶部４１０に蓄積された値に基づいて更新する（ステップＳ３０８）。
【００５５】
《作用・効果》
このように、第５の実施形態によれば、推定部４０４は、各管理のタイミングにおいて、前回の管理のタイミングにおいて更新された統計モデルを用いて、状態量の値を推定することができる。確率特定部４０６は、更新された統計モデルによる推定値の出現確率を特定することとなる。つまり、第５の実施形態によれば、各管理のタイミングにおいて統計データだけでなく、統計モデル自体も更新することで、統計推定値をより精度よく推定することができる。
なお、第５の実施形態に係る管理装置４０は、統計モデルを過去の状態量の値に基づいて更新するが、これに限られない。例えば他の実施形態においては管理装置４０は、状態量記憶部４１０に状態量を蓄積させる一方で統計モデルの更新をしないものであってもよい。この場合においても、過去の状態量の値が蓄積することにより、統計モデルによる推定の精度が向上し得る。例えば、データの蓄積により「大数の法則」で平均値の推定値が真値に近づくこと、分散の範囲を狭く絞り込むことにより、推定精度の向上が期待できる。
【００５６】
〈他の実施形態〉
以上、図面を参照して一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、様々な設計変更等をすることが可能である。
例えば、上述した実施形態に係る管理システム１における管理装置４０は、対象装置１０の管理に用いる値の抽出および特定をする機能を有するが、これに限られない。例えば、他の実施形態に係る管理システム１においては、管理装置４０と別個に対象装置１０の管理に用いる値の抽出および特定をする情報処理装置を備え、管理装置４０は、情報処理装置が特定した値を用いて対象装置１０を管理してもよい。
また例えば、上述した実施形態に係る管理装置４０は、ネットワークＮを介して計測値を取得するが、これに限られない。例えば、他の実施形態に係る管理装置４０は、計測器２０から直接計測値を取得してもよい。この場合、管理システム１は通信装置３０を備えなくてもよい。
【００５７】
また、上述した実施形態によれば、管理装置４０は、複数の推定値の中から１つを選択することで、これを対象装置１０の管理に用いる状態量の値とするが、これに限られない。例えば、他の実施形態においては、管理装置４０は、出現確率に応じた重み係数による推定値の加重平均値を求めることで、これを対象装置１０の管理に用いる状態量の値としてもよい。各推定値の重み係数は、出現確率に対して単調増加する。管理装置４０は、例えば、推定値の出現確率をそのまま重み係数として用いてもよい。
【００５８】
また、上述した実施形態によれば、管理装置４０は、欠落が検出された値を推定により求めるが、これに限られない。例えば、他の実施形態においては、管理装置４０は、欠落の有無に関わらず推定により状態量の値を求め、計測値と推定値のそれぞれについて確率分布に基づいて対象装置１０の管理に用いる値を特定してもよい。
【００５９】
〈コンピュータ構成〉
図１２は、少なくとも１つの実施形態に係るコンピュータの構成を示す概略ブロック図である。
コンピュータ９０は、ＣＰＵ９１、主記憶装置９２、補助記憶装置９３、インタフェース９４を備える。
上述の管理装置４０は、コンピュータ９０に実装される。そして、上述した各処理部の動作は、プログラムの形式で補助記憶装置９３に記憶されている。ＣＰＵ９１は、プログラムを補助記憶装置９３から読み出して主記憶装置９２に展開し、当該プログラムに従って上記処理を実行する。また、ＣＰＵ９１は、プログラムに従って、モデル記憶部４０３および確率分布記憶部４０５に対応する記憶領域を主記憶装置９２に確保する。
【００６０】
補助記憶装置９３の例としては、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ（Compact Disc Read Only Memory）、ＤＶＤ−ＲＯＭ（Digital Versatile Disc Read Only Memory）、半導体メモリ等が挙げられる。補助記憶装置９３は、コンピュータ９０のバスに直接接続された内部メディアであってもよいし、インタフェース９４または通信回線を介してコンピュータ９０に接続される外部メディアであってもよい。また、このプログラムが通信回線によってコンピュータ９０に配信される場合、配信を受けたコンピュータ９０が当該プログラムを主記憶装置９２に展開し、上記処理を実行してもよい。少なくとも１つの実施形態において、補助記憶装置９３は、一時的でない有形の記憶媒体である。
【００６１】
また、当該プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、当該プログラムは、前述した機能を補助記憶装置９３に既に記憶されている他のプログラムとの組み合わせで実現するもの、いわゆる差分ファイル（差分プログラム）であってもよい。
【符号の説明】
【００６２】
１  管理システム
１０  対象装置
２０  計測器
３０  通信装置
４０  管理装置（情報処理装置）
４０１  計測値取得部
４０２  欠落検出部
４０３  モデル記憶部
４０４  推定部
４０５  確率分布記憶部
４０６  確率特定部
４０７  管理値特定部
４０８  管理部
４０９  評価値算出部
４１０  状態量記憶部
４１１  モデル更新部
【図１】
2018169748.tif 000003
【図２】
2018169748.tif 000004
【図３】
2018169748.tif 000005
【図４】
2018169748.tif 000006
【図５】
2018169748.tif 000007
【図６】
2018169748.tif 000008
【図７】
2018169748.tif 000009
【図８】
2018169748.tif 000010
【図９】
2018169748.tif 000011
【図１０】
2018169748.tif 000012
【図１１】
2018169748.tif 000013
【図１２】
2018169748.tif 000014
【図１３】
2018169748.tif 000015
【図１４】
2018169748.tif 000016