(57)【要約】【課題】所定区域に設置された対象設備それぞれの消費電力を、より簡易な構成でより高精度に推定できる消費電力推定装置を提供する。【解決手段】運転状態が監視可能な１以上の対象設備１００それぞれの消費電力を推定する消費電力推定装置１０は、合計消費電力の時系列データを取得する合計電力ベクトル生成部３２と、１以上の対象設備１００の運転パラメータの時系列データを取得する状態行列生成部３４と、１以上の参照信号を生成する参照信号生成部３６と、合計消費電力を目的変数、1以上の対象設備１００の運転パラメータおよび１以上の参照信号の成分値を説明変数として重回帰分析することで、１以上の対象設備１００それぞれの寄与度を算出する寄与度推定部３８と、対象設備１００の寄与度と運転パラメータとを乗算することで、対象設備１００の消費電力を算出する内訳算出部４２と、を備える。【選択図】図２
2020004041.tif 000002
【特許請求の範囲】
【請求項１】
所定区域に設置されるとともにその運転状態が監視可能な１以上の対象設備それぞれの消費電力を推定する消費電力推定装置であって、
前記所定区域全体での消費電力である合計消費電力の時系列データを取得する合計消費電力取得部と、
前記１以上の対象設備の運転状態を数値化した運転パラメータの時系列データを取得する運転状態取得部と、
それぞれが所定の基底関数で表される１以上の参照信号を生成する参照信号生成部と、
取得された前記合計消費電力を目的変数、取得された前記1以上の対象設備の運転パラメータおよび生成された前記１以上の参照信号の成分値を説明変数として重回帰分析することで、前記１以上の対象設備それぞれの前記合計消費電力に対する寄与度を算出する寄与度推定部と、
前記対象設備の前記寄与度と前記運転パラメータとを乗算することで、前記対象設備の消費電力を算出する消費電力算出部と、
を備えることを特徴とする消費電力推定装置。
【請求項２】
請求項１に記載の消費電力推定装置であって、
前記消費電力算出部は、前記１以上の対象設備それぞれの消費電力を各時間ごとに算出することで、前記合計消費電力の各時間ごとの内訳を生成する、ことを特徴とする消費電力推定装置。
【請求項３】
請求項１または２に記載の消費電力推定装置であって、
前記基底関数は、略山形または略矩形の立ち上がりが１回だけ発生する関数である、ことを特徴とする消費電力推定装置。
【請求項４】
請求項３に記載の消費電力推定装置であって、
各参照信号は、前記基底関数で表される信号を、単位時間ごとに繰り返す、ことを特徴とする消費電力推定装置。
【請求項５】
請求項３または４に記載の消費電力推定装置であって、
前記１以上の参照信号は、１以上の信号グループを含み、
各信号グループは、互いに同一の形状、かつ、互いに位相が異なる複数の参照信号を含む、
ことを特徴とする消費電力推定装置。
【請求項６】
請求項４に記載の消費電力推定装置であって、
同一の信号グループに属する前記複数の参照信号は、幅Ｗの山を持つ略山形の関数で表されており、Ｗ／２ずつ位相がずれている、
ことを特徴とする消費電力推定装置。
【請求項７】
請求項５に記載の消費電力推定装置であって、
前記１以上の参照信号は、２以上の信号グループを含んでおり、
一つの信号グループに属する参照信号の山の幅は、他の信号グループに属する参照信号の山の幅と異なる、
ことを特徴とする消費電力推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本明細書は、所定区域に設置されるとともにその運転状態が監視可能な１以上の対象設備それぞれの消費電力を推定する消費電力推定装置を開示する。
【背景技術】
【０００２】
ビルや工場等の建物において、省エネを目的として電力の管理を行う場合には、その内部に設置されている複数の電気設備それぞれの消費電力を把握する必要がある。しかし、このような消費電力を計測するためには、各電気設備に電力量メータを設ける必要があり、コストがかかっていた。
【０００３】
そこで、一部では、所定区域で消費される合計消費電力の時系列データと、所定区域に設置される設備の駆動状況を示す時系列データと、に基づいて回帰分析することで、各設備の消費電力を推定する技術が提案されている。
【０００４】
例えば、特許文献１では、所定区域に設置された複数の空調機それぞれの消費電力を算出する装置が開示されている。この特許文献１では、所定区域で消費される合計消費電力と、空調機以外で消費される電力であるベース電力と、から複数の空調機で消費する全空調消費電力を求め、この全空調消費電力の時系列データと、複数の空調機それぞれの駆動状態を示す時系列データと、に基づいて各空調機の消費電力を求めている。かかる技術によれば、設備ごと電力計測機器を設ける必要がないため、低コストで、複数の設備それぞれの消費電力を取得できる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００４−１３８２６５号公報
【特許文献２】特開２０１７−０５０９７１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１の技術では、空調機以外で消費される電力であるベース電力の推定精度が低いため、結果として、設備ごとの消費電力の推定精度も低いという問題があった。すなわち、特許文献１では、特定の日時における合計消費電力（例えば、特定の条件を満たすある一日のうち、空調機が駆動していない時間帯の合計消費電力）を、空調機以外で消費される電力、すなわち、ベース電力と推定しており、このベース電力の値は、日、月が変わっても一定であるとみなしている。
【０００７】
しかし、空調機以外で消費される電力としては、照明や、パソコン、オフィス機器、調理機器など様々あり、これらの使用量（消費電力）は、刻々と変化するのが普通である。したがって、特許文献１のように、ベース電力を一定とみなした場合には、現実との乖離が大きくなり、空調機ごとの消費電力の推定精度が悪化する。
【０００８】
そこで、一部では、消費電力推定装置で、その駆動状況が監視できない設備、例えば、オフィス機器などで消費される電力量を、所定区域内に存在する人員の数から推定することも提案されている（例えば特許文献２など）。しかし、この場合、人員の数を管理するためのシステムが別途必要であり、コストアップを招く。
【０００９】
そこで、本明細書では、所定区域に設置された対象設備それぞれの消費電力を、より簡易な構成でより高精度に推定できる消費電力推定装置を開示する。
【課題を解決するための手段】
【００１０】
本明細書で開示する消費電力推定装置は、所定区域に設置されるとともにその運転状態が監視可能な１以上の対象設備それぞれの消費電力を推定する消費電力推定装置であって、前記所定区域全体での消費電力である合計消費電力の時系列データを取得する合計消費電力取得部と、前記１以上の対象設備の運転状態を数値化した運転パラメータの時系列データを取得する運転状態取得部と、それぞれが所定の基底関数で表される１以上の参照信号を生成する参照信号生成部と、取得された前記合計消費電力を目的変数、取得された前記1以上の対象設備の運転パラメータおよび生成された前記１以上の参照信号の成分値を説明変数として重回帰分析することで、前記１以上の対象設備それぞれの前記合計消費電力に対する寄与度を算出する寄与度推定部と、前記対象設備の前記寄与度と前記運転パラメータとを乗算することで、前記対象設備の消費電力を算出する消費電力算出部と、を備えることを特徴とする。
【００１１】
この場合、前記消費電力算出部は、前記１以上の対象設備それぞれの消費電力を各時間ごとに算出することで、前記合計消費電力の各時間ごとの内訳を生成してもよい。
【００１２】
また、前記基底関数は、略山形または略矩形の立ち上がりが１回だけ発生する関数であってもよい。この場合、各参照信号は、前記基底関数で表される信号を、単位時間ごとに繰り返してもよい。さらに、前記１以上の参照信号は、１以上の信号グループを含み、各信号グループは、互いに同一の形状、かつ、互いに位相が異なる複数の参照信号を含んでもよい。この場合、同一の信号グループに属する前記複数の参照信号は、幅Ｗの山を持つ略山形の関数で表されており、Ｗ／２ずつ位相がずれていてもよい。さらに、前記１以上の参照信号は、２以上の信号グループを含んでおり、一つの信号グループに属する参照信号の山の幅は、他の信号グループに属する参照信号の山の幅と異なっていてもよい。
【発明の効果】
【００１３】
本明細書に開示の消費電力推定装置によれば、複数の参照信号で、対象設備以外で消費される電力を模擬するため、所定区域に設置された対象設備それぞれの消費電力を、より簡易な構成でより高精度に推定できる。
【図面の簡単な説明】
【００１４】
【図１】消費電力推定装置の物理構成を示すブロック図である。
【図２】消費電力推定装置の機能構成を示すブロック図である。
【図３】参照信号の一例を示す図である。
【図４】複数の参照信号を重ね合わせた図である。
【図５】回帰モデルを説明する図である。
【図６】寄与度の推定処理を示すフローチャートである。
【図７】消費電力の内訳の推定処理を示すフローチャートである。
【図８】本例による消費電力の実測値と推定値とを示す図である。
【図９】他の参照信号の一例を示す図である。
【図１０】他の参照信号の一例を示す図である。
【図１１】従来技術による消費電力の実測値と推定値とを示す図である。
【発明を実施するための形態】
【００１５】
以下、図面を参照して消費電力推定装置１０の構成について説明する。図１は、消費電力推定装置１０の物理的構成を、図２は、機能的構成を示すブロック図である。この消費電力推定装置１０は、所定区域に設置された１以上の対象設備それぞれの消費電力を回帰分析により、推定する装置である。すなわち、一般に、ビルや工場等の建物には、その建物全体、あるいは、フロア全体での消費電力を計測する電力量メータが設けられている。この電力量メータを用いれば、所定区域（例えば建物全体またはフロア全体）の合計消費電力は、把握できる。
【００１６】
しかし、省エネを目的として電力を管理する場合には、合計消費電力だけでなく、所定区域内に設置された１以上の対象設備それぞれの消費電力（以下「個別消費電力」という）も把握することが望まれる。対象設備とは、建物管理システムで運転管理されている設備であり、例えば、空調設備である。こうした対象設備の個別消費電力を正確に把握するためには、対象設備ごとに電力量メータを設ける必要がある。しかし、対象設備ごとに電力量メータを設けることは、コスト増加を招く。そこで、消費電力推定装置１０は、こうした対象設備の個別消費電力を、合計消費電力、および、対象設備の運転状態に基づいて推定する。以下では、一つのビル全体を「所定区域」とし、このビル全体に設置された複数の空調設備を「対象設備」とする消費電力推定装置１０を例に挙げて説明する。
【００１７】
図１に示すように、消費電力推定装置１０は、特定のプログラムに従って駆動する１以上のコンピュータで構築される。すなわち、消費電力推定装置１０は、図１に示すように、ＣＰＵ１２と、記憶装置１４と、入力装置１６と、出力装置１８と、通信インターフェース（以下「通信Ｉ／Ｆ」という）２０と、を備えており、これらが、バス２２により接続されている。
【００１８】
ＣＰＵ１２は、各種演算を行なう。具体的には、ＣＰＵ１２は、必要に応じて、記憶装置１４に記憶される消費電力推定プログラムを読み込んで、消費電力の推定に必要な各種演算を行なう。このＣＰＵ１２による処理の具体的内容については、後述する。
【００１９】
記憶装置１４は、各種データ（プログラム含む）を記憶するもので、例えば、ＲＡＭやＲＯＭといった半導体メモリ、ハードディスクドライブ、ソリッドステートドライブ等を１以上組み合わせて構成される。この記憶装置１４には、後述する消費電力の推定を行なうための消費電力推定プログラムが記憶されている。
【００２０】
入力装置１６は、オペレータからの操作指示およびデータ入力を受け付けるもので、例えば、キーボード、マウス、タッチパネル、マイク等を１以上、組み合わせて構成される。出力装置１８は、オペレータに各種情報を出力するもので、例えば、ディスプレイ、スピーカー等を１以上、組み合わせて構成される。
【００２１】
通信Ｉ／Ｆ２０は、他の電子機器および情報端末とデータを遣り取りする。例えば、建物管理システムで検知された合計消費電力や、対象設備（空調設備）の運転状態は、この通信Ｉ／Ｆ２０を介して、消費電力推定装置１０に入力される。また、消費電力推定装置１０による診断結果は、この通信Ｉ／Ｆ２０を介して、他の情報端末に送られてもよい。なお、図１では、消費電力推定装置１０を、１台のコンピュータとして図示しているが、消費電力推定装置１０は、互いに通信可能な複数のコンピュータで構成されてもよい。また、消費電力推定装置１０は、建物管理システムの一部であってもよい。
【００２２】
次に、消費電力推定装置１０の機能構成について図２を参照して説明する。ビルには、複数の対象設備１００と、複数の非監視設備１０２と、が設置されている。対象設備１００は、個別消費電力の推定対象となる設備であり、建物管理システム１１０で、その運転が管理されている設備である。本例では、ビルに設置された空調設備、例えば、空調室内機、空調室外機等を対象設備１００としている。
【００２３】
非監視設備１０２は、電力を消費する設備のうち、建物管理システム１１０で管理されていない設備である。したがって、非監視設備１０２としては、例えば、照明や、コンセント等が含まれる。対象設備１００および非監視設備１０２は、いずれも、電力供給ライン１０３を介して電源１０６に接続されており、当該電源１０６から電力が供給される。対象設備１００および非監視設備１０２に供給された電力の合計値、すなわち、合計消費電力は、電力量メータ１０４で計測される。
【００２４】
建物管理システム１１０は、予め設定された管理プログラムに従って、対象設備１００の運転を管理する。この建物管理システム１１０の構成は、適宜、変更可能であるが、建物管理システム１１０は、少なくとも、運転状態記憶部１１４と、合計消費電力記憶部１１２と、を有している。
【００２５】
運転状態記憶部１１４は、複数の対象設備１００それぞれの運転状態の時系列データを記憶する。より具体的には、建物管理システム１１０は、各対象設備１００の運転状態を定期的に取得して、日時と対応づけて運転状態記憶部１１４に記憶する。このとき、建物管理システム１１０は、運転状態を数値で表す。例えば、運転中は「１」、運転停止中は「０」という数値で表す。以下では、この運転状態を示す数値を、運転パラメータと呼ぶ。この運転パラメータの値は、上述したように、運転中なら「１」、運転停止中なら「０」という２値であってもよい。また、別の形態として、運転パラメータの値は、暖房運転中なら「１」、冷房運転中なら「０．８」、換気運転中なら「０．３」、運転停止中なら「０」というように、運転種類に応じて多段階に設定されてもよい。さらに、別の形態として、運転パラメータの値は、空調の目標温度と現在の室温との差分値や、コンプレッサーの回転数などに応じて設定されてもよい。いずれにしても、運転状態記憶部１１４には、複数の対象設備１００それぞれの運転状態を示す数値（運転パラメータ値）の時系列データが記憶されている。
【００２６】
合計消費電力記憶部１１２は、ビル全体の消費電力の合計値（合計消費電力）の時系列データを記憶する。具体的には、建物管理システム１１０は、電力量メータ１０４で計測された単位時間当たりの合計消費電力を定期的に取得して、日時と対応付けて合計消費電力記憶部１１２に記憶する。この合計消費電力のサンプリングタイミングと、対象設備１００の運転状態のサンプリングタイミングは、一致させることが望ましい。また、サンプリング周期は、特に限定されないが、３０秒以上、１時間以内の値に設定されることが望ましく、１分以上、１０分以内の値に設定されることがより望ましい。
【００２７】
また、本例では、運転パラメータ値および合計消費電力の時系列データを、建物管理システム１１０が記憶するとしているが、これらのデータは、消費電力推定装置１０に設けられた記憶装置１４に記憶されてもよい。
【００２８】
消費電力推定装置１０は、合計消費電力を目的変数、運転パラメータ値および後述する参照信号の成分を説明変数として重回帰分析することで、複数の対象設備１００それぞれの個別消費電力を推定する。この消費電力推定装置１０の機能構成を説明する前に、個別消費電力の推定の原理について説明する。
【００２９】
あるサンプリングタイミングｎにおける合計消費電力ｙｎは、複数の対象設備１００の個別消費電力の合計値Ｘｎと、複数の非監視設備１０２の個別消費電力の合計値Ｒｎと、の和となる。すなわち、ｙｎ＝Ｘｎ＋Ｒｎとなる。以下では、Ｘｎを「対象消費電力」、Ｒｎを「非監視消費電力」と呼ぶ。
【００３０】
第ｍ対象設備１００の個別消費電力は、その運転状態（運転パラメータｘｍｎ）に依存すると考えると、第ｍ対象設備１００の個別消費電力は、（ｘｍｎ・ａｍ）と表すことができる。ここで、ａｍは、対象設備１００それぞれに設定される寄与度である。そして、第１〜第Ｍ対象設備１００の個別消費電力の合計値である対象消費電力Ｘｎは、式１で表すことができる。
【数１】
2020004041.tif 000003
【００３１】
一方、非監視設備１０２は、運転状態を把握できないため、対象設備１００のようなモデルを構築することはできない。そこで、本例では、非監視設備１０２の個別消費電力の合計である非監視消費電力Ｒｎを、複数の参照信号ｒｉ（ｉ＝１，２，・・・，Ｉ）を用いて模擬する。図３、図４は、参照信号の一例を示す図である。図３は、第１〜第５参照信号ｒ１〜ｒ５を示す図であり、図４は、第１〜第５参照信号ｒ１〜ｒ５を重ね合わせた図である。
【００３２】
参照信号は、所定の基底関数で表される信号を単位時間ごと（例えば２４時間ごと）に繰り返した信号である。この基底関数の種類は、特に限定されないが、略山形または略矩形の立ち上がりが１回だけ発生する関数であることが望ましい。本例では、基底関数として、式２に示すガウス関数を用いている。
【数２】
2020004041.tif 000004
【００３３】
ガウス関数は、山形（より正確には釣鐘形）の関数であるが、その山の幅Ｗは、σに依存し、山の中心は、μに依存する。第ｉ参照信号ｒｉは、式２で表される信号を、２４時間ごとに繰り返した信号である。図３、図４の例では、各参照信号ｒ１〜ｒ５の山の幅Ｗは、約１２時間＝７２０分となっている。また、複数の参照信号の位相は、山の幅Ｗの約１／２である６時間＝３６０分ずつずれている。なお、ここで、山の幅Ｗとは、ｒｉの値が、ピーク値の１％を超えてから１％未満になるまでの期間とする。
【００３４】
本例では、非監視消費電力Ｒｎを、第１〜第Ｉ参照信号ｒ１〜ｒＩを用いて、式３のように表す。なお、ｂｉは、参照信号ｒｉそれぞれに設定される寄与度である。また、ｃは、定数項であり、時間に依存することなく、常時、発生する消費電力を意味している。
【数３】
2020004041.tif 000005
【００３５】
ここで、これまでの説明において、対象設備１００の運転パラメータｘｍｎと、参照信号ｒｉｎと、合計消費電力ｙｎは、既知であるが、寄与度ａｍおよびｂｉと、定数項ｃが未知である。この寄与度ａｍおよびｂｉと定数項ｃを解くことができれば、対象設備１００の個別消費電力ｘｍｎ・ａｍを求めることができる。そこで、本例では、消費電力を式４でモデル化し、重回帰分析により、寄与度ａｍおよびｂｉと定数項ｃとを求める。この式４の各項目についての説明を図５に示す。
【数４】
2020004041.tif 000006
【００３６】
この式４は、Ｎ＞（Ｍ＋Ｉ＋１）であれば解くことができる。したがって、（Ｍ＋Ｉ＋１）回以上、合計消費電力ｙｎおよび運転パラメータｘｍｎが収集できれば、寄与度ａｍ，ｂｉと定数項ｃを算出できる。そして、各対象設備１００の運転パラメータｘｍｎに、算出された寄与度ａｍを乗算することで、サンプリングタイミングｎにおける第ｍ対象設備１００の消費電力を算出できる。なお、以下では、式４の左辺を「合計電力ベクトル」と呼ぶ。また、右辺にある一つ目の行列を「状態行列」と呼ぶ。
【００３７】
消費電力推定装置１０は、こうした演算が行なえるような機能構成となっている。そこで、再び、図２を参照して消費電力推定装置１０の機能構成について説明する。合計電力ベクトル生成部３２は、合計消費電力記憶部１１２に記憶されている合計消費電力ｙｎの時系列データから、合計電力ベクトル、すなわち、式４における左辺のベクトルを生成する。したがって、合計電力ベクトル生成部３２は、合計消費電力ｙｎの時系列データを取得する合計消費電力取得部として機能する。なお、このとき、合計電力ベクトルを構成するサンプル数Ｎは、対象設備１００の個数Ｍと、参照信号の個数Ｉと、１と、を加算した値（Ｍ＋Ｉ＋１）以上であれば、特に限定されない。ただし、サンプル数Ｎは、２４時間分以上の個数であることが望ましい。
【００３８】
状態行列生成部３４は、運転状態記憶部１１４に記憶されている運転パラメータ値ｘｍｎの時系列データ、および、参照信号生成部３６で生成された参照信号ｒｉｎに基づいて、状態行列を生成する。なお、状態行列生成部３４は、運転パラメータ値ｘｍｎの時系列データを取得する運転状態取得部として機能する。参照信号生成部３６は、状態行列生成部３４からの要求に応じて、参照信号ｒｉｎを生成する。
【００３９】
生成された合計電力ベクトルおよび状態行列は、寄与度推定部３８に入力される。寄与度推定部３８は、合計電力ベクトルおよび状態行列ベクトルを、式４に当てはめたうえで、寄与度ａｍおよびｂｉと、定数項ｃを解く。この算出には、公知の重回帰分析技術を用いることができるため、詳細な説明は省略する。
【００４０】
寄与度推定部３８で算出された寄与度ａｍ，ｂｉと定数項ｃは、寄与度記憶部４０に記憶される。内訳算出部４２は、対象設備１００の運転パラメータｘｍｎに、寄与度記憶部４０に記憶されているａｍを乗算して、対象設備１００の個別消費電力を算出するもので、消費電力推定部として機能する。また、式３に、参照信号ｒｉｎ、および、寄与度記憶部４０に記憶されているｂｉ，ｃを当てはめて、複数の非監視設備１０２全ての消費電力、すなわち非監視消費電力Ｒｎを算出する。そして、これにより、各時間の消費電力の内訳が算出される。算出された内訳（各対象設備１００の個別消費電力と非監視消費電力Ｒｎ）は、時刻と対応付けられて消費電力記憶部４４に記憶される。
【００４１】
次に、この消費電力推定装置１０における消費電力推定処理の流れについて図６、図７を参照して説明する。消費電力推定処理は、寄与度の推定処理と、消費電力の内訳の推定処理と、に大別される。図６は、寄与度の推定処理の流れを示している。この寄与度の推定処理は、少なくとも、１回、実行されればよい。ただし、推定精度を向上するためには、この寄与度の推定処理を、ある程度の時間を空けて、繰り返し実行し、寄与度記憶部４０に記憶されている寄与度および定数項の値を更新するのが望ましい。すなわち、複数の対象設備１００それぞれの寄与度ａｍは、単位時間（上記の例ではサンプリング周期）当たりの個別消費電力を表している。この単位時間当たりの個別消費電力は、周辺環境や、対象設備１００の劣化等によって変化する。寄与度を定期的に推定し、更新することで、こうした対象設備１００の劣化等も考慮した推定結果を得ることができる。
【００４２】
寄与度を推定する場合には、まず、合計電力ベクトル生成部３２が、合計消費電力ｙｎの時系列データに基づいて、合計電力ベクトルを生成する（Ｓ１０）。また、消費電力推定装置１０が、運転パラメータｘｍｎの時系列データと、参照信号ｒｉｎと、に基づいて、状態行列を生成する（Ｓ１２）。
【００４３】
合計電力ベクトルおよび状態行列ベクトルが生成されれば、寄与度推定部３８は、式４の回帰モデルに従って、寄与度ａｍ，ｂｉ、定数項ｃを重回帰分析により求める（Ｓ１４）。そして、算出された寄与度ａｍ，ｂｉ、定数項ｃが、寄与度記憶部４０に記憶されれば、寄与度の推定処理は、終了となる（Ｓ１６）。
【００４４】
図７は、消費電力の内訳の推定処理の流れを示すフローチャートである。この消費電力の内訳の推定処理は、寄与度が算出された後であれば、どのようなタイミングで何回実行されてもよい。したがって、内訳の推定処理は、定期的に行われてもよいし、ユーザから指示されたタイミングで行われてもよい。いずれにしても、内訳を推定する場合、内訳算出部４２は、内訳を推定したい期間中に取得された複数の対象設備１００それぞれの運転パラメータｘｍｎの時系列データを読み込む（Ｓ２０）。また、内訳算出部４２は、寄与度記憶部４０に記憶されている複数の寄与度ａｍ，ｂｉおよび定数項ｃも読み込む（Ｓ２２）。続いて、内訳算出部４２は、対象設備１００の運転パラメータｘｍｎに、寄与度ａｍを乗算した値を、対象設備１００の個別消費電力として算出する（Ｓ２４）。また、式３に、参照信号ｒｉｎ、および、寄与度記憶部４０に記憶されているｂｉ，ｃを当てはめて、非監視消費電力Ｒｎを算出する（Ｓ２６）。そして、これらの消費電力を消費電力記憶部４４に記憶すれば、内訳の推定処理は終了となる（Ｓ２８）。
【００４５】
以上の説明から明らかなとおり、本例では、非監視設備１０２の消費電力の合計値である非監視消費電力Ｒｎを、複数の参照信号ｒｉｎで模擬している。そのため、非監視消費電力Ｒｎを一定とみなす特許文献１に比べて、誤差を低減することができ、対象設備１００の個別消費電力の推定精度をより向上できる。また、消費電力を、人員の数で模擬する特許文献２と異なり、本例では、人員の数を消費電力推定に利用していない。そのため、人員の数をモニタリングするシステム（例えば入館システムなど）が不要となり、対象設備１００の個別消費電力をより簡易に推定できる。
【００４６】
そこで、本例と従来技術との推定精度の比較について、図８、図１１を参照して説明する。図８、図１１は、それぞれ、あるビルの２４時間分の消費電力の実測値Ｒ＃，Ｘ＃と推定値Ｒ＊，Ｘ＊とを示すグラフである。より詳しく説明すると、図８、図１１では、複数の非監視設備１０２の消費電力の合計値（非監視消費電力Ｒ）および複数の対象設備１００の個別消費電力の合計値（対象消費電力Ｘ）の積み上げグラフである。図８、図１１において、太破線は、非監視消費電力の実測値Ｒ＃を、太実線は、当該実測値Ｒ＃の上に積み上げて描画した対象消費電力の実測値Ｘ＃を示している。同様に、図８、図１１において、細破線は、非監視消費電力の推定値Ｒ＊を、細実線は、当該推定値Ｒ＊の上に積み上げて描画した対象消費電力の推定値Ｘ＊を示している。
【００４７】
ここで、図１１では、特許文献１に開示の技術に基づいて、推定値Ｒ＊，Ｘ＊を推定している。すなわち、図１１では、推定値Ｒ＊が常に一定であるとみなしたうえで、合計消費電力の実測値と、複数の対象設備１００の運転状態と、に基づいて、複数の対象設備１００それぞれの個別消費電力を重回帰分析で求めている。すなわち、図１１は、式５に示す回帰モデルに基づいて、各寄与度ａｍを求めている。この式５から明らかなとおり、従来技術において、非監視消費電力の推定値Ｒ＊は、回帰モデルでは既知の定数項として扱われる。
【数５】
2020004041.tif 000007
【００４８】
図１１から明らかなとおり、従来技術では、非監視消費電力の推定値Ｒ＊を一定とみなすため、この非監視消費電力の推定値Ｒ＊と実測値Ｒ＃との誤差が大きくなる。回帰モデルで用いられる定数項であるＲ＊の誤差が大きいと、当然ながら、回帰モデルを用いた推定精度も悪化し、対象設備１００の個別消費電力の推定精度が悪化する。これは、図１１において、対象設備１００の個別消費電力の合計値である対象消費電力の推定値Ｘ＊と実測値Ｘ＃との誤差が大きいことからも明らかである。
【００４９】
一方、本例では、複数の参照信号ｒｉを用いて、非監視消費電力を模擬している。その結果、図８に示す通り、本例では、非監視消費電力の推定値Ｒ＊と実測値Ｒ＃との誤差が小さい。そして、非監視消費電力を精度よく推定できるため、本例では、対象設備１００の個別消費電力の推定精度も向上している。これは、図８において、対象消費電力の推定値Ｘ＊と実測値Ｘ＃との誤差が小さいことからも明らかである。
【００５０】
なお、これまで説明した構成は一例であり、少なくとも、合計消費電力の時系列データを目的変数、複数の対象設備の運転パラメータの時系列データおよび複数の参照信号の成分を説明変数として、重回帰分析して、複数の対象設備それぞれの消費電力を求めるのであれば、その他の構成は、適宜、変更されてもよい。例えば、複数の参照信号の形態は、適宜、変更されてもよい。例えば、上述の説明では、幅Ｗが１２時間の山形の信号を参照信号として用いているが、参照信号の山の幅Ｗは、図９に示すように、適宜、変更されてもよい。図９において、参照信号ｒｉの幅Ｗは、上段から順に、幅Ｗが、２４時間、１６時間、１２時間、８時間となっている。また、いずれの場合でも、複数の参照信号ｒｉは、Ｗ／２ずつ位相がずれている。
【００５１】
また、複数の参照信号ｒｉは、２以上の信号グループを含んでおり、一つの信号グループに属する参照信号ｒｉの山の幅Ｗは、他の信号グループに属する参照信号ｒｉの山の幅Ｗと異なってもよい。例えば、複数の参照信号ｒｉは、図９に示した四種類の参照信号ｒｉを組み合わせてもよい。図１０に示すような複数の参照信号ｒｉで、非監視消費電力を模擬してもよい。
【００５２】
また、これまでは、参照信号として、ガウス関数で表される信号を用いてきたが、参照信号の形は、適宜、変更されてもよい。ただし、参照信号を表す基底関数は、山形または矩形の立ち上がりが１回だけ生じる関数であることが望ましい。
【符号の説明】
【００５３】
１０  消費電力推定装置、１２  ＣＰＵ、１４  記憶装置、１６  入力装置、１８  出力装置、２０  通信Ｉ／Ｆ、２２  バス、３２  合計電力ベクトル生成部（合計消費電力取得部）、３４  状態行列生成部（運転状態取得部）、３６  参照信号生成部、３８  寄与度推定部、４０  寄与度記憶部、４２  内訳算出部（消費電力推定部）、４４  消費電力記憶部、１００  対象設備、１０２  非監視設備、１０３  電力供給ライン、１０４  電力量メータ、１０６  電源、１１０  建物管理システム、１１２  合計消費電力記憶部、１１４  運転状態記憶部。
【図１】
2020004041.tif 000008
【図２】
2020004041.tif 000009
【図３】
2020004041.tif 000010
【図４】
2020004041.tif 000011
【図５】
2020004041.tif 000012
【図６】
2020004041.tif 000013
【図７】
2020004041.tif 000014
【図８】
2020004041.tif 000015
【図９】
2020004041.tif 000016
【図１０】
2020004041.tif 000017
【図１１】
2020004041.tif 000018