(57)【要約】【課題】設備に関する未知のパラメータの推定精度を向上できるようにすること。【解決手段】設備３に関するパラメータを推定するパラメータ推定装置１は、設備に関する物理量を測定するセンサ２から観測値を取得し、観測値を対象パラメータに基づいて算出する設備モデル式を取得し（４１）、設備モデル式に対象パラメータの推定値を代入して得られる推定観測値の観測値に対する尤度を第１尤度関数により算出し、対象パラメータの推定値があらかじめ用意された所定範囲内（４２）に属する確率を求める第２尤度関数と第１尤度関数とに基づいて、対象パラメータの推定値の尤度を算出し、算出された尤度に基づいて対象パラメータを推定し（４３）、推定された対象パラメータを所定の評価基準にしたがって評価する（４４）。【選択図】図１
2020095451.tif 000002
【特許請求の範囲】
【請求項１】
設備に関するパラメータを推定する装置であって、
前記設備に関する物理量を測定するセンサから観測値を取得し、
前記観測値を対象パラメータに基づいて算出する設備モデル式を取得し、
前記設備モデル式に前記対象パラメータの推定値を代入して得られる推定観測値の前記観測値に対する尤度を第１尤度関数により算出し、
前記対象パラメータの推定値があらかじめ用意された所定範囲内に属する確率を求める第２尤度関数と前記第１尤度関数とに基づいて、前記対象パラメータの推定値の尤度を算出し、
前記算出された尤度に基づいて前記対象パラメータを推定し、
前記推定された対象パラメータを所定の評価基準にしたがって評価する、
パラメータ推定装置。
【請求項２】
前記対象パラメータの評価結果を出力する、
パラメータ推定装置。
【請求項３】
前記対象パラメータの推定値の前記所定範囲における分布密度が前記所定範囲外での分布密度よりも大きくなるように前記対象パラメータの尤度が設定されている、
請求項１に記載のパラメータ推定装置。
【請求項４】
前記設備モデル式は、設備の種類に応じてあらかじめ用意されている、
請求項１に記載のパラメータ推定装置。
【請求項５】
前記センサから得られる観測値の数よりも前記対象パラメータの数が大きく設定されている、
請求項１に記載のパラメータ推定装置。
【請求項６】
前記設備モデル式には、異なる複数の対象パラメータが含まれている、
請求項１に記載のパラメータ推定装置。
【請求項７】
前記所定範囲は、一つの連続した範囲である、
請求項１に記載のパラメータ推定装置。
【請求項８】
前記所定範囲は、複数の連続した範囲からなる、
請求項１に記載のパラメータ推定装置。
【請求項９】
前記所定範囲は、あらかじめ設定される単位時間で変化する、
請求項１に記載のパラメータ推定装置。
【請求項１０】
前記所定範囲は、複数の対象パラメータそれぞれに関する前記所定範囲の組合せにしたがって前記第２尤度関数と前記第１尤度関数とに基づいて、前記複数の対象パラメータの推定値の尤度をそれぞれ算出する、
請求項１に記載のパラメータ推定装置。
【請求項１１】
前記推定された複数の対象パラメータを所定の評価基準にしたがって評価する場合、前記所定範囲の組合せに含まれていない所定範囲を持つ対象パラメータに対しても当該対象パラメータに関する前記所定範囲を設定して評価する、
請求項１０に記載のパラメータ推定装置。
【請求項１２】
設備に関するパラメータを計算機により推定する方法であって、
前記計算機は、
前記設備に関する物理量を測定するセンサから観測値を取得し、
前記観測値を対象パラメータに基づいて算出する設備モデル式を取得し、
前記設備モデル式に前記対象パラメータの推定値を代入して得られる推定観測値の前記観測値に対する尤度を第１尤度関数により算出し、
前記対象パラメータの推定値があらかじめ用意された所定範囲内に属する確率を求める第２尤度関数と前記第１尤度関数とに基づいて、前記対象パラメータの推定値の尤度を算出し、
前記算出された尤度に基づいて前記対象パラメータを推定し、
前記推定された対象パラメータを所定の評価基準にしたがって評価する、
パラメータ推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、パラメータ推定装置および方法に関する。
【背景技術】
【０００２】
エネルギを効率的に使用するため、ビルや商店施設にＢＥＭＳ（Ｂｕｉｌｄｉｎｇ  Ｅｎｅｒｇｙ  Ｍａｎａｇｅｍｅｎｔ  Ｓｙｓｔｅｍ）が普及しつつある。ＢＥＭＳは、建物の熱源設備、空調設備、電力設備等のユーティリティ設備を管理する。ＢＥＭＳは、ユーティリティ設備にセンサおよび通信装置を付けて計測し、見える化し、管理することで、設備の保全および設備の最適化を実行する。
【０００３】
さらに近年では、ＣＥＭＳ（Ｃｏｍｍｕｎｉｔｙ  Ｅｎｅｒｇｙ  Ｍａｎａｇｅｍｅｎｔ  Ｓｙｓｔｅｍ）の重要性が高まっている。ＣＥＭＳは、複数のＢＥＭＳを管理し、地域全体のエネルギを管理する。ＢＥＭＳおよびＣＥＭＳの主要な機能は、設備状態の見える化である。設備状態とは、消費電力だけではなく、効率、流量、温度、圧力など、ＢＥＭＳおよびＣＥＭＳが管理するユーティリティ設備の動作状態を表す値である。しかし、設備の持つすべての状態を知るために、多種多様なセンサをユーティリティ設備に設けることは、高コストであり、現実的ではない。
【０００４】
そこで、一つの状態を検知するセンサデータに基づいて、他の状態の値を推定することが考えられる。すなわち、ＢＥＭＳまたはＣＥＭＳの管理下にある計測データを用いて、未知の状態を推定することにより、検知対象の状態毎にセンサを取り付けずに、設備状態を管理することができる。
【０００５】
ところで、ユーティリティ設備の状態を推定するためには、パラメータを推定する必要がある。しかし、未知のパラメータを推定する場合、物理的な値ではない値に推定値が収束する場合がある。未知のパラメータの推定精度を向上する技術として、特許文献１の技術が提案されている
【０００６】
特許文献１では、「分電盤接続されている複数の電気機器による消費電力の合計値を用いて、前記複数の電気機器の稼動状態を確率モデルによりモデル化したときのモデルパラメータを推定するパラメータ推定ステップを含み、前記パラメータ推定ステップでは、前記複数の電気機器それぞれの動作傾向から事前知識として予め定めることのできる電力データの特徴に基づいて、尤度関数によって算出される尤度が最大になるようなモデルパラメータを推定し、前記確率モデルは、Ｆａｃｔｏｒｉａｌ  ＨＭＭ（Ｈｉｄｄｅｎ  Ｍａｒｋｏｖ  Ｍｏｄｅｌ）であり、前記尤度は、前記Ｆａｃｔｏｒｉａｌ  ＨＭＭによってモデル化された前記消費電力の合計値のパターンの、実測した前記消費電力の合計値に対する確からしさを示す値である。」という記載がある。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１７−４９２２１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかし、特許文献１では、Ｆａｃｔｏｒｉａｌ  ＨＭＭ以外のモデルを用いる場合に限定されてしまい、他のモデルを用いて未知パラメータを推定することは難しい。
【０００９】
ところで、ユーティリティ設備モデルの多数のパラメータを容易に推定するためには、観測値の種類の数が多い必要がある。観測値とは、センサなどの計測装置により実測された値である。しかし多種類の観測値を得るためには、上述のように、センサなどの計測装置をユーティリティ設備に取り付ける必要があり、コスト上の問題がある。
【００１０】
さらに、観測値の種類の数が少ない場合は、いわゆる劣問題となるため、推定値が収束しない可能性がある。したがって、限られた種類の観測値を用いて、真値に近い推定値を取得する手法が必要である。
【００１１】
本発明は、上記課題に鑑みてなされたもので、その目的は、設備に関する未知のパラメータの推定精度を向上できるようにしたパラメータ推定装置および方法を提供することにある。
【課題を解決するための手段】
【００１２】
上記課題を解決すべく、本発明に従うパラメータ推定装置は、設備に関するパラメータを推定する装置であって、設備に関する物理量を測定するセンサから観測値を取得し、観測値を対象パラメータに基づいて算出する設備モデル式を取得し、設備モデル式に対象パラメータの推定値を代入して得られる推定観測値の観測値に対する尤度を第１尤度関数により算出し、対象パラメータの推定値があらかじめ用意された所定範囲内に属する確率を求める第２尤度関数と第１尤度関数とに基づいて、対象パラメータの推定値の尤度を算出し、算出された尤度に基づいて対象パラメータを推定し、推定された対象パラメータを所定の評価基準にしたがって評価する。
【発明の効果】
【００１３】
本発明によれば、第１尤度関数と第２尤度関数とに基づいて算出される対象パラメータの推定値の尤度から、対象パラメータを推定し、評価することができる。
【図面の簡単な説明】
【００１４】
【図１】パラメータ推定装置を含むシステム全体の構成図である。
【図２】パラメータ推定装置の構成図である。
【図３】パラメータ推定処理を示すフローチャートである。
【図４】設備モデル式を決定する処理を示すフローチャートである。
【図５】対象パラメータの適正範囲を設定する処理を示すフローチャートである。
【図６】対象パラメータの適正範囲の設定例を示す説明図である。
【図７】対象パラメータの適正範囲の他の設定例を示す説明図である。
【図８】対象パラメータの適正範囲のさらに別の設定例を示す説明図である。
【図９】第２実施例に係り、適正範囲の設定処理を示すフローチャートである。
【発明を実施するための形態】
【００１５】
以下、図面に基づいて、本発明の実施の形態を説明する。本実施形態に係るパラメータ推定装置では、対象パラメータ（以下、パラメータまたは未知パラメータと呼ぶ場合がある）の推定値が所定範囲での分布密度が高くなるように、尤度を設定して学習する。
【００１６】
本実施形態では、未知パラメータの推定値を確率的に自動算出し、評価する。パラメータの推定値が、所定範囲において尤度が高くなる所定の分布に従うと仮定することにより、対象パラメータの良否を判定する。
【００１７】
本実施形態によれば、パラメータの適正範囲に基づく尤度関数を設定し、観測値の推定値の尤度関数と適正範囲に基づく尤度関数とを用いてパラメータを推定し、パラメータの良否を判断することができる。したがって、本実施形態では、ユーティリティ設備における観測値の種類の数よりも多い数のパラメータを推定する場合に、真値に近い推定値を算出することができる。
【００１８】
本実施形態に係るパラメータ推定装置およびパラメータ推定方法は、例えば、以下のように表現することもできる。
【００１９】
「表現１」ユーティリティ設備のパラメータ推定装置において、パラメータの適正範囲に基づき、尤度関数を設定し、パラメータを推定し、パラメータを判定するパラメータ推定装置。
【００２０】
「表現２」ユーティリティ設備のパラメータ推定方法において、パラメータの適正範囲に基づき、尤度関数を設定し、パラメータを推定し、パラメータを判定するパラメータ推定方法。
【００２１】
「表現３」表現２記載のパラメータ推定方法において、一つの連続した適正範囲に基づく尤度設定をするステップを含むパラメータ推定方法。
【００２２】
「表現４」表現２記載のパラメータ推定方法において、複数の連続区間を有する適正範囲に基づく尤度設定をするステップを含むパラメータ推定方法。
【００２３】
「表現５」表現２記載のパラメータ推定方法において、推定パラメータに対して、一定時間の変化量の適正範囲に基づく尤度設定をするステップを含むパラメータ推定方法。
【００２４】
「表現６」表現２記載のパラメータ推定方法において、設備モデル式を記憶部から検索し、推定するパラメータを決定するステップを含むパラメータ推定方法。
【００２５】
「表現７」ユーティリティ設備のパラメータを推定する方法において、複数のパラメータの適正範囲に基づき、複数の尤度関数の組合せ生成して設定し、パラメータを推定し、パラメータを判定するパラメータ推定方法。
【００２６】
「表現８」表現２記載のパラメータ推定方法において、各尤度設定の組合せで、適正範囲を用いてないパラメータに対しても、パラメータ判定ステップで、適正範囲設定ステップで設定した適正範囲を適用し、良否を判定するパラメータ推定方法。
【実施例１】
【００２７】
図１〜図８を用いて第１実施例を説明する。本実施例では、ビルディング内のユーティリティ設備についての未知パラメータを推定する場合を例に挙げて説明するが、これに限らず、例えば、工場の生産設備、商業施設の冷蔵設備、遊戯設備、エレベータ、エスカレータ、食品配送設備等の他の施設の他の設備にも適用することができる。
【００２８】
＜全体の構成＞
【００２９】
図１は、パラメータ推定装置１を含むシステムの全体構成を示す。パラメータ推定装置１は、例えば、ＢＥＭＳ４とセンサ２とに通信ネットワークＣＮを介して接続される。センサ２は、ユーティリティ設備（以下、設備）３と接続されており、計測情報（センサデータ、観測値とも呼ぶ）を、ＢＥＭＳ４またはパラメータ推定装置１へ送信する。
【００３０】
設備３は、例えば、ボイラ設備、空調装置、ガスタービン発電装置、電力関連機器等である。センサ２は、設備３の各種物理量を計測し、その計測結果である計測情報を出力する装置である。センサ２は、例えば、消費電力、効率、流量、温度、圧力、液位、質量、音、光量などを計測する。色彩、味、匂いといった心理物理量を計測するセンサであってもよい。
【００３１】
センサ２で計測されたデータ（計測情報）は、通信ネットワークＣＮを介してパラメータ推定装置１へ直接的に送信されてもよいし、ＢＥＭＳ４などの装置と通信ネットワークＣＮとを介してパラメータ推定装置１へ送信されてもよい。さらに、センサ２は、定期的にデータを送信してもよいし、パラメータ推定装置１またはＢＥＭＳ４からの要求に応じてデータを送信してもよい。あるいは、センサ２の計測したデータを記録媒体に格納し、記録媒体からパラメータ推定装置１へ読み込ませてもよい。
【００３２】
なお、本実施例では、ＢＥＭＳ４を例に挙げるが、これに限らず、ＣＥＭＳなどのユーティリティ設備３を管理するシステムであればよい。
【００３３】
パラメータ推定装置１は、例えばＣＰＵ（マイクロプロセッサ）、メモリ、通信機能、表示機能、入力機能（いずれも図示せず）などを備えた装置である。パラメータ推定装置１は、パーソナルコンピュータを利用して構成することもできるし、専用コンピュータとして構成することもできる。
【００３４】
パラメータ推定装置１の機能構成の詳細は、図２と共に後述する。ここでは、パラメータ推定装置１の概要を先に説明する。パラメータ推定装置１は、設備モデル式決定部４１と、適正範囲設定部４２と、パラメータ推定部４３と、パラメータ判定部４４といった機能を有している。
【００３５】
図２は、パラメータ推定装置１の機能構成を示す。パラメータ推定装置１は、例えば、取得部１０と、インタフェース部２０と、記憶部３０と、演算部４０とを含む。これら各機能１０〜４０は、メモリに格納された所定のコンピュータプログラム（不図示）をマイクロプロセッサが読み込んで実行することにより、実現される。機能１０〜４０のいずれか少なくとも一部をハードウェア回路として構成してもよい。
【００３６】
インタフェース部２０は、例えば、入力部２１と出力部２２とを有する。演算部４０は、例えば、設備モデル式決定部４１と、適正範囲設定部４２と、パラメータ推定部４３と、パラメータ判定部４４とを備える。
【００３７】
取得部１０は、設備３のセンサ２、ＢＥＭＳ４（あるいはＣＥＭＳ）、インターネット上の情報源（不図示）からデータ（情報とも呼ぶ）を取得する。取得されたデータは、記憶部３０に格納される。
【００３８】
インタフェース部２０の入力部２１は、パラメータ推定装置１に対してユーザが情報を入力するための機能である。ユーザは、例えばキーボードまたはマウスなどの入力装置（不図示）を用いて、設備３などの情報を入力部２１へ入力する。キーボードなどに代えて、あるいはキーボードなどと共に、音声認識装置などの入力装置を用いてもよい。
【００３９】
入力部２１に入力された設備３などの情報は、記憶部３０に転送されて記憶される。入力部２１に入力された情報の一部または全部を、出力部２２により表示出力させることもできる。さらに、入力部２１に入力された情報は、記憶部３０を介して演算部４０に送られる。
【００４０】
入力される情報としては、例えば、設備３の仕様、設備３の設置場所、設備３と他の設備との接続関係、設備３について推定するパラメータがある。入力情報の例は、これらに限定されない。
【００４１】
出力部２２は、ユーザに対して情報を提供する機能である。出力部２２は、例えば、ディスプレイ、プリンタ、音声合成装置などの出力装置（不図示）を用いて、パラメータの推定結果などを表示する。出力部２２は、入力部２１へ入力された情報を表示したり、演算部４０から記憶部３０に格納された情報を読み出して表示したり、記憶部３０に記憶されている情報を読み出して表示したりすることができる。
【００４２】
記憶部３０は、取得部１０と、インタフェース部２０と、演算部４０との間で、情報を送受信する。記憶部３０は、入力された情報を格納すると共に、要求された情報を取り出して要求元へ送信する。記憶部３０には、設備３の仕様情報などのほかに、後述する設備モデル式と適正範囲とパラメータ初期値などが格納される。
【００４３】
設備モデル式とは、取得部１０から取得するデータとインタフェース部２０から取得する推定パラメータとの関係を示す式である。設備モデル式は、取得部１０から取得されるデータと推定パラメータとの関係を表現可能な複数の式として構成することができ、物理式を含んでもよい。「所定範囲」としての適正範囲とは、設備種類、設備仕様、設置環境といったタグ（種類情報）が設定された、各パラメータの適正範囲である。
【００４４】
演算部４０は、設備モデル式決定部４１と、適正範囲設定部４２と、パラメータ推定部４３と、パラメータ判定部４４とを含む。演算部４０は、インタフェース部２０と記憶部３０とから演算に必要な情報を取得して対象パラメータを推定し、インタフェース部２０に推定結果を送信する。対象パラメータの推定結果は、出力部２２に表示される。
【００４５】
パラメータの推定結果を含む情報を、設備導入のシミュレーションシステムなどの外部装置（不図示）へ提供してもよい。パラメータ推定装置１から外部装置へ有線または無線により情報を転送してもよいし、記録媒体に格納して外部装置へ提供してもよい。
【００４６】
＜パラメータ推定装置１の全体の処理の概要＞
【００４７】
図３は、パラメータ推定装置１の全体処理を示すフローチャートである。ステップＳ１２〜Ｓ１５の詳細は後述する。
【００４８】
パラメータ推定装置１は、インタフェース部２０と記憶部３０とからデータを取得する（Ｓ１１）。パラメータ推定装置１は、ステップＳ１１で取得されたデータに基づいて、設備モデル式を決定する（Ｓ１２）。パラメータ推定装置１は、ステップＳ１１で取得されたデータに基づいて、適正範囲を設定する（Ｓ１３）。
【００４９】
ここで、ステップＳ１２とステップＳ１３とは、並列で処理されてもいい。ステップＳ１２を先に処理し、ステップＳ１３を後で処理してもよい。ステップＳ１２とステップＳ１３の処理の順番は問わない。
【００５０】
パラメータ推定装置１は、ステップＳ１２で決定された設備モデル式とステップＳ１３で設定された適正範囲とに基づいて、対象パラメータを推定する（Ｓ１４）。
【００５１】
パラメータ推定装置１は、ステップＳ１４で算出されたパラメータの推定値を評価する（Ｓ１５）。評価とは、例えば、算出されたパラメータの推定値の良否を判定することである。
【００５２】
パラメータ推定装置１は、ステップＳ１５で「否」と判定されると、ステップＳ１４へ戻ってパラメータ推定を再度行う。一方、パラメータ推定装置１は、ステップＳ１５で「良」と判定された場合、ステップＳ１４で算出されたパラメータの推定値を出力する（Ｓ１６）。
【００５３】
＜設備モデル式決定部４１の処理＞
【００５４】
図４は、図３中にステップＳ１２で示された、設備モデル式決定部４１での処理の詳細を示すフローチャートである。
【００５５】
設備モデル式決定部４１は、図３のステップＳ１１で取得されたデータに基づいて、記憶部３０に格納されている設備モデル式の中から演算に使用する設備モデル式を一つ以上検索する（Ｓ１２１）。
【００５６】
記憶部３０に格納されている設備モデル式には、例えば、「設備種類」、「式の左辺のパラメータ」、「式の右辺のパラメータ」といったタグが付いている。そこで、設備モデル式決定部４１は、タグに基づいて設備モデル式を検索する。設備モデル式の検索には、クラスタリングを用いればいいが、この限りではない。
【００５７】
設備モデル式決定部４１は、ステップＳ１２１で検索された設備モデル式のパラメータを確認し、推定対象のパラメータを決める（Ｓ１２２）。以下、ステップＳ１２１とステップＳ１２２の一例を、蒸気ボイラのガス消費量を計算する式を用いて説明する。
【００５８】
Ｇ＝（Ｈ＋Ｈｂ）／Ｈｒｏｕｔ×Ｇｒｏｕｔ・・・式（１）
【００５９】
Ｈ＝ｆ（Ｓ，Ｅｓ，Ｅｗ）        ・・・式（２）
【００６０】
Ｈｂ＝ｆ（Ｓ，Ｅｋ，Ｅｗ，Ｂｒ）・・・式（３）
【００６１】
Ｅｓ＝ｆ（Ｐ）                  ・・・式（４）
【００６２】
Ｅｗ＝ｆ（Ｔｗ）                ・・・式（５）
【００６３】
Ｅｋ＝ｆ（Ｔｋ）                ・・・式（６）
【００６４】
Ｔｋ＝ｇ（Ｐ）                  ・・・式（７）
【００６５】
ここで、「Ｇ」はガス消費量、「Ｈ」は熱量、「Ｈｂ」はブロー損失熱量、「Ｈｒｏｕｔ」は定格熱出力、「Ｇｒｏｕｔ」は定格ガス消費量、「Ｓ」は蒸気量、「Ｅｓ」は蒸気エンタルピ、「Ｅｗ」は給水エンタルピ、「Ｅｋ」は缶水エンタルピ、「Ｂｒ」はブロー係数、「Ｐ」は使用圧力、「Ｔｗ」は給水温度、「Ｔｋ」は缶水温度である。
【００６６】
各式は、物理モデル式または近似式である。例えば、給水エンタルピＥｗの式（５）は給水温度の二次近似式であり、ガス消費量Ｇの式（１）は物理モデル式である。
【００６７】
ステップＳ１１において、インタフェース部２０から、推定パラメータとしてブロー係数Ｂｒを、設備種類として蒸気ボイラをそれぞれ取得するものとする。さらに、ステップＳ１１では、取得部１０と入力部２１とから、蒸気量Ｓ、使用圧力Ｐ、給水温度Ｔｗ、定格熱出力Ｈｒｏｕｔ、定格ガス消費量Ｇｒｏｕｔを取得するものとする。この場合、ステップＳ１２１では、記憶部３０の設備モデル式を検索することにより、式（１）、式（２）、式（３）、式（４）、式（５）、式（６）、式（７）を取得する。
【００６８】
設備モデル式決定部４１は、ステップＳ１１において推定パラメータとして取得されたブロー係数Ｂｒ以外に推定が必要なパラメータを決める（Ｓ１２２）。ここでは、ステップＳ１１において、ブロー係数Ｂｒ以外のすべてのパラメータを取得できたため、推定パラメータはブロー係数Ｂｒと決定する。これに対し、もし蒸気量Ｓが取得部１０から取得できない場合は、蒸気量Ｓも推定パラメータであると決定される。
【００６９】
＜適正範囲設定部４２の処理＞
【００７０】
図５は、図３中にステップＳ１３で示された、適正範囲設定部４２での処理の詳細を示すフローチャートである。
【００７１】
適正範囲設定部４２は、設備モデル式決定部４１により決定された推定パラメータの適正範囲を記憶部３０から検索し、検索された適正範囲を推定対象パラメータに設定する（Ｓ１３１）。すなわち、適正範囲設定部４２は、設備種類と推定対象パラメータの名称とをタグとして、記憶部３０が格納している適正範囲の中から検索し、各推定対象パラメータに対する適正範囲を設定する（Ｓ１３１）。
【００７２】
適正範囲の検索には、例えばクラスタリングを用いることができるが、検索手法は問わない。適正範囲を見つけることのできない推定対象パラメータについては、適正範囲を設定しない。
【００７３】
図６〜図８に示すように、適正範囲には、様々なパターンがある。図６は、例えば、給水温度の適正範囲が「２０℃〜４０℃」と一つの連続した範囲の場合である。図７は、給水温度の適正範囲が「１０℃〜２０℃」と「３０℃〜４０℃」との、複数の連続した区間を有する場合である。
【００７４】
図８は、ある一定時間（ｔ１−ｔ２）において温度変化をするパラメータの適正範囲を表す。点Ａは、時刻ｔ１における温度Ｔ１である。点Ｂ１は、時刻ｔ２において取り得る最低温度Ｔ２である。点Ｂ２は、時刻ｔ２において取り得る最高温度Ｔ３である。
【００７５】
図８の例において、ある一定時間（ｔ１とｔ２の差分）による温度変化範囲は、ΔＴ１からΔＴ２の間である。このように、時間差分（ｔ１−ｔ２）によるパラメータの変化量ΔＴを適正範囲として設定することも可能である。
【００７６】
図５に戻る。適正範囲設定部４２は、ステップＳ１３１で設定された適性範囲を用いて、尤度関数を設定する（Ｓ１３２）。尤度関数の設定方法を、蒸気ボイラを例にあげて詳細に説明する。
【００７７】
設備モデル式決定部４１により決定された設備モデル式が、式（１）、式（２）、式（３）、式（４）、式（５）、式（６）、式（７）であるとし、さらに、ガス消費量Ｇの計測値が記憶部３０に格納されていると仮定する。
【００７８】
式（１）、式（２）、式（３）、式（４）、式（５）、式（６）、式（７）の各パラメータに値を入力することにより、ガス消費量Ｇの推定値（推定観測値）を算出することができる。算出されたガス流量Ｇの推定値の尤度は、尤度関数により計算できる。例えば、計測値と推定値との誤差である観測誤差が正規分布であるとき、尤度関数は以下の式（８）で求められる。式（８）において、「Ｘ」は確率関数、「μ」は平均、「σ２」は分散である。ステップＳ１３２の説明として、図５にも式（８）を示す。
【００７９】
ｆ（Ｘ：μ，σ２）＝（１／（２πσ２）１／２）ｅｘｐ（−（Ｘ−μ）２／２σ２）
・・・式（８）
【００８０】
適正範囲設定部４２は、ステップＳ１３１で設定されたパラメータの適正範囲を用いて、尤度関数を更新する。
【００８１】
例えば、給水温度が推定対象パラメータであり、ステップＳ１３１で適正範囲を「２０℃〜４０℃」と設定した場合、パラメータの尤度関数として一様分布を適用する。ここでは一様分布を用いるが、一様分布以外の分布も適用できる。給水温度の推定値が２０℃から４０℃の範囲内の場合、その尤度は「１／（４０−２０）」であり、それ以外の場合は尤度が「０」である。
【００８２】
適正範囲設定部４２は、パラメータの尤度関数を式（８）に示す尤度関数に乗算することにより、尤度関数を設定する。ここでは図６の適正範囲に基づく尤度関数の設定を説明したが、図７または図８に示した適正範囲に基づく尤度関数設定も、図６の場合の尤度関数設定を例として応用することが可能である。推定対象パラメータが複数であり、それぞれの適正範囲がある場合は、式（８）の尤度関数にそれぞれのパラメータの尤度関数を乗算することにより、尤度関数を設定すればよい。
【００８３】
＜パラメータ推定部４３の処理＞
【００８４】
パラメータ推定部４３は、設備モデル式決定部４１により決定された設備モデル式と適正範囲設定部４２により設定された推定対象パラメータの適正範囲とに基づく尤度関数を用いることにより、パラメータを推定する。パラメータ推定手法としては、例えば粒子フィルタ法や、マルコフ連鎖モンテカルロ法を用いることができる。
【００８５】
＜パラメータ判定部４４の処理＞
【００８６】
パラメータ判定部４４は、パラメータ推定部４３により算出されたパラメータの推定値の良否を判定する。良否の判定基準としては、例えば、推定値の収束度、適正範囲、観測誤差などがある。
【００８７】
推定値の収束度には、推定ステップごとの推定値の分散を用いればよい。これに代えて、標準偏差などを用いてもよい。推定値の収束度に基づく良否判定では、予め設定した閾値を基準に判定される。適正範囲に基づく良否判定では、推定値が適正範囲内に収束した場合に「良」と判定し、収束しない場合に「否」と判定する。観測誤差に基づく良否判定では、予め設定された閾値を基準に判定される。推定値の収束度、適正範囲、観測誤差の複数の判定基準（評価基準）に基づく良否判定において、いずれか一つでも「否」と判定された場合には、パラメータ推定部４３において再度推定処理を行う（Ｓ１４）。
【００８８】
＜出力部２２の処理＞
【００８９】
出力部２２は、入力部２１、記憶部３０、および演算部４０からの情報を受信して、出力装置に表示する。出力部２２は、演算部４０のパラメータ判定部４４による判定基準毎の判定結果と、推定値と、適正範囲とをそれぞれ出力装置に表示する。
【００９０】
このように構成される本実施例によれば、推定対象パラメータを効率よく正確に推定することができる。本実施例によれば、観測値の種類の数よりも多いパラメータを推定する場合であっても、適正範囲内の推定値の尤度が高くなるように設定された尤度関数を用いることにより、各推定対象パラメータの推定値を適正範囲内に収束させることができ、真値に近い推定値を効率的に取得することができる。
【００９１】
なお、本実施例では、適正範囲設定部４２において様々なパターンの適正範囲を用いることにより、推定対象パラメータについてより多くの情報を得ることができる。したがって、この場合は、より一層、真値に近い推定値を効率的に取得することができる。
【実施例２】
【００９２】
図９を用いて第２実施例を説明する。本実施例は、第１実施例の変形例であるため、第１実施例との相違を中心に説明する。
【００９３】
第１実施例では、適正範囲設定部４２において尤度関数を設定する場合、すべての推定対象パラメータの尤度関数を乗算することにより、全体としての尤度関数を設定する場合を説明した。しかし、もしパラメータ判定部４４での判定結果が「否」であった場合、パラメータ推定部４３で再度推定することになるが、再推定してもその判定結果が「否」となってしまい、推定処理（Ｓ１４）と判定処理（Ｓ１５）とが繰り返される可能性が考えられる。
【００９４】
図９は、本実施例における適正範囲設定部４２の処理Ｓ１３Ａを示す。図９に示すフローチャートと図５で述べたフローチャートとを比較すると、本実施例ではステップＳ１３２Ａが第１実施例のステップＳ１３２と異なる。そこで、この異なるステップＳ１３２Ａについて説明する。
【００９５】
適正範囲設定部４２は、ステップＳ１３１で設定された適正範囲に基づいて尤度関数を設定する場合に、各推定対象パラメータの適正範囲の組合せを計算して、尤度関数を設定する（Ｓ１３２Ａ）。
【００９６】
例えば、推定対象のパラメータがパラメータＰａとパラメータＰｂとの２つある場合、尤度設定の組合せは、パラメータＰａの適正範囲による尤度設定、パラメータＰｂの適正範囲による尤度設定、パラメータＰａの適正範囲とパラメータＰｂの適正範囲とによる尤度設定の、合計３つある。
【００９７】
そして、本実施例のパラメータ推定部４３は、それぞれの組合せの尤度設定を用いて、パラメータ推定処理を行う。本実施例のパラメータ判定部４４は、それぞれの組合せの尤度設定によるパラメータ推定値の良否を判定する。
【００９８】
パラメータ判定部４４は、第１実施例で述べた判定基準を用いることができる。その判定基準の中で、適正範囲の設定されていない推定対象パラメータに対しても、ステップＳ１３１で設定された適正範囲を適用して、良否を判定する。
【００９９】
例えば、パラメータＰａの適正範囲だけに基づいて尤度関数を設定する組合せの場合、パラメータ判定部４４は、ステップＳ１３１で設定されたパラメータＰｂの適正範囲に基づいて、パラメータＰｂの推定値の良否を判定する。
【０１００】
このように構成される本実施例も第１実施例と同様の作用効果を奏する。さらに、本実施例では、複数の尤度設定の組合せを用いて推定するため、パラメータ判定部４４で「良」と判定されるパラメータ推定値を取得する確率を高めることができる。
【０１０１】
なお、本発明は、上述した実施形態に限定されない。当業者であれば、本発明の範囲内で、種々の追加や変更等を行うことができる。上述の実施形態において、添付図面に図示した構成例に限定されない。本発明の目的を達成する範囲内で、実施形態の構成や処理方法は適宜変更することが可能である。
【０１０２】
また、本発明の各構成要素は、任意に取捨選択することができ、取捨選択した構成を具備する発明も本発明に含まれる。さらに特許請求の範囲に記載された構成は、特許請求の範囲で明示している組合せ以外にも組み合わせることができる。
【符号の説明】
【０１０３】
１：パラメータ推定装置、２：センサ、３：設備、４１：設備モデル式決定部、４２：適正範囲設定部、４１：パラメータ推定部、４２：パラメータ判定部
【図１】
2020095451.tif 000003
【図２】
2020095451.tif 000004
【図３】
2020095451.tif 000005
【図４】
2020095451.tif 000006
【図５】
2020095451.tif 000007
【図６】
2020095451.tif 000008
【図７】
2020095451.tif 000009
【図８】
2020095451.tif 000010
【図９】
2020095451.tif 000011