JP7006859B2 - 空調制御装置、空調システム、空調制御方法、空調制御プログラム - Google Patents

Description

本発明は、機械学習モデルを用いて空調機の起動時刻を推定する空調制御装置、空調システム、空調制御方法、及び空調制御プログラムに関する。

消費電力を抑えつつ、ユーザーにとって快適な空調を実行するように空調機を制御する技術の例として、例えば、特許文献１に記載の空調制御装置が挙げられる。当該空調制御装置は、室温変化の履歴を示す室温履歴情報を、空調機の動作履歴情報に対応付け、これらの情報に基づいて、空調機が温度を調節しない場合の室温をオフ時予測室温として予測し、当該オフ時予測室温に基づいて、目標時刻において室温を目標温度に到達させるための制御パラメータを決定する。

上記のオフ時予測室温の予測について、より詳細には、特許文献１に記載の空調制御装置は、室温履歴情報と動作履歴情報とに基づいて、空調機が温度を調節しない場合の居室の将来の室温を予測するための機械学習モデルであるオフ時室温予測モデルを作成し、当該オフ時室温予測モデルを用いて、空調機の起動時刻を決定する。

特開２０１７－６７４２７号公報

しかし、特許文献１に記載の空調制御装置における機械学習モデルの作成は、十分に蓄えられた学習データ、例えば室温履歴情報及び動作履歴情報のデータが存在することを前提としている。そして、十分なデータが蓄えられていない状態で学習を行った機械学習モデルでは、不適切な予測をする可能性があるという問題があった。

本発明は、上記の課題を解決するためになされたものであり、十分なデータが蓄えられていない状態で学習を行った機械学習モデルを用いた場合において、空調機の起動時刻について推定精度を向上した空調制御装置を得ることを目的とする。

密領域と密領域より学習データが少ない疎領域とを含む学習データ群を用いて、空調機の制御に関する情報である空調データから空調機の起動時刻を推定するための学習を行った機械学習モデルに基づいて、空調機の起動時刻を推定する空調制御装置において、空調データを取得する空調データ取得部と、空調データが、疎領域に存在するか判定する判定部と、判定部が、空調データが疎領域に存在すると判定した場合、機械学習モデルに対して、空調データと起動時刻とを関連付ける非学習モデルによる外挿を行い、非学習モデルを用いて空調データから空調機の起動時刻を推定する推定部と、推定部が推定した起動時刻に空調機を起動する制御を行う制御部と、を備えた。

本発明に係る空調制御装置は、判定部が、空調データが疎領域に存在すると判定した場合、密領域と密領域より学習データが少ない疎領域とを含む学習データ群を用いて、空調機の制御に関する情報である空調データから空調機の起動時刻を推定するための学習を行った機械学習モデルに対して、非学習モデルによる外挿を行い、非学習モデルを用いて空調データから空調機の起動時刻を推定する推定部を備えたので、学習データが少ない領域では非学習モデルによる外挿を行い、空調機の起動時刻を推定することにより、推定精度を向上することができる。

実施の形態１に係る空調システム１０００の構成を示す構成図である。 実施の形態１に係る空調制御装置１００の構成を示す構成図である。 乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。 乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。 乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。 乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。 乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。 推定部１６が学習済みの機械学習モデルに対して非学習モデルによる外挿を行い、当該非学習モデルを用いて空調データから空調機の起動時刻を推定する方法の具体例を示す説明図である。 空調制御装置１００を実現するコンピュータのハードウェア構成の例を示す構成図である。 空調制御装置１００が機械学習モデルの学習を行う動作を示すフローチャートである。 空調制御装置１００が空調機ＡＣの起動時刻を推定する動作を示すフローチャートである。 乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。

実施の形態１．
図１は、実施の形態１に係る空調システム１０００の構成を示す構成図である。空調システム１０００は、空調制御装置１００、複数の室外機ＯＵ（ＯＵ１，ＯＵ２，・・・ＯＵｎ）及び複数の室内機ＩＵ（ＩＵ１１，ＩＵ１２，・・・ＩＵ１ｐ，ＩＵ２１，ＩＵ２２，・・・ＩＵ２ｑ，ＩＵｎ１，ＩＵｎ２，・・・ＩＵｎｒ）を備える。ここで、ｎ、ｐ、ｑ、及びｒは任意の正の整数である。以下では、特定の室外機を指す場合には、ＯＵの後に追加の符号を付し、任意の室外機または全室外機を指す場合にはＯＵのみで表記することとする。室内機についても同様に、特定の室内機を指す場合には、ＩＵの後に追加の符号を付し、任意の室内機または全室内機を指す場合にはＩＵのみで表記することとする。符号の付し方については、後述する空調機ＡＣ、室外環境センサーＯＵＳ、室内環境センサーＩＵＳ、受付部ＩＵＲについても同様である。
各室外機ＯＵには複数の室内機ＩＵが接続されており、例えば室外機ＯＵ１には室内機ＩＵ１（ＩＵ１１，ＩＵ１２，・・・ＩＵ１ｐ）が接続される。空調制御装置１００と室内機ＩＵ及び室外機ＯＵは電磁的方法により接続され、接続される方法は有線か無線かを問わない。また、空調制御装置１００は室内機ＩＵ及び室外機ＯＵが設置された建築物と同じ建築物に設置される構成でも、遠隔地に備えられ、ネットワーク経由で接続される構成であってもよい。

各室外機ＯＵは、接続された室内機ＩＵとともに室内の空調を行う空調機ＡＣ（ＡＣ１，ＡＣ２，・・・ＡＣｎ）を構成する。また、各室外機ＯＵは、当該室外機ＯＵが設置された室外の環境を検知し、室外環境情報を取得する室外環境センサーＯＵＳ（ＯＵＳ１，ＯＵＳ２，・・・ＯＵＳｎ）を備える。各室外環境センサーＯＵＳは、取得した室外環境情報を空調データとして空調制御装置１００に出力する。ここで、室外環境情報とは、室外の環境を示すものであり、例えば、各時刻における室外温度や室外湿度等である。

各室内機ＩＵは、当該室内機ＩＵが設置された室内の環境を検知し、室内環境情報を取得する室内環境センサーＩＵＳ（ＩＵＳ１１，ＩＵＳ１２，・・・ＩＵＳ１ｐ，ＩＵＳ２１，ＩＵＳ２２，・・・ＩＵＳ２ｑ，ＩＵＳｎ１，ＩＵＳｎ２，・・・ＩＵＳｎｒ）を備える。また、各室内機ＩＵは、ユーザーから設定情報を受け付ける受付部ＩＵＲ（ＩＵＲ１１，ＩＵＲ１２，・・・ＩＵＲ１ｐ，ＩＵＲ２１，ＩＵＲ２２，・・・ＩＵＲ２ｑ，ＩＵＲｎ１，ＩＵＲｎ２，・・・ＩＵＲｎｒ）を備える。各室内機ＩＵは、取得した室内環境情報、受け付けた設定情報、及び空調機ＡＣの運転状態を示す運転情報を、空調データとして室外機ＯＵを介して空調制御装置１００に出力する。ここで、室内環境情報とは、室内の環境を示すものであり、例えば、各時刻における室内温度及び室内湿度等である。設定情報は、受付部を介してユーザーにより設定される情報であり、例えば、室内機ＩＵが設置された部屋の環境値を目標値に到達させる目標時刻、目標温度、及び目標湿度等である。ここで、部屋の環境値とは、室内環境情報が示す値であり、目標値とは、設定情報が示す部屋の環境値の目標となる値である。
運転情報は、空調機ＡＣの運転状態を示すものであり、例えば、空調機ＡＣの運転能力、空調機ＡＣが起動状態か停止状態かを示す情報、並びに、空調機ＡＣが冷房モード、暖房モード及び除湿モードを含む運転モードのいずれのモードにあるかを示す情報等である。また、室内機ＩＵだけでなく、室外機ＯＵも運転情報を出力する構成としてもよい。ここで、運転能力とは、空調機ＡＣの空調能力の強さを示すものであり、例えば、室外機ＯＵの圧縮機周波数、冷媒蒸発温度、凝縮温度、及び稼働する室内機ＩＵの台数等である。

また、上記において、室外機ＯＵ及び室内機ＩＵが、それぞれ室外環境センサーＯＵＳ及び室内環境センサーＩＵＳを備える構成としたが、各種センサーは室外機ＯＵ及び室内機ＩＵと独立して設ける構成としても良い。あるいは、室外環境情報については、室外環境センサーＯＵＳからではなく、外部のサーバー等から取得する構成としても良い。

図２は、実施の形態１に係る空調制御装置１００の構成を示す構成図である。空調制御装置１００は、密領域と密領域より学習データが少ない疎領域とを含む学習データ群を用いて、空調機の制御に関する情報である空調データから空調機の起動時刻を推定するための学習を行った機械学習モデルに基づいて、空調機の起動時刻を推定するものであり、送受信部１０、制御部１１、記憶部１２、空調データ取得部１３、機械学習部１４、判定部１５、及び推定部１６を備える。
ここで、学習データが多い領域を密領域、学習データが少ない領域を疎領域と呼ぶこととする。密領域には学習を行うのに十分な量の学習データが存在し、推定部１６は当該学習データを用いて学習を行った機械学習モデルから適切な起動時刻を推定することができる。

機械学習モデルの学習方法は周知技術を用いればよく、以下では、具体例として教師あり学習の場合について説明する。他の具体例として、例えば、強化学習を用いても良く、その場合、環境値が目標値に到達するまでに要した何らかの指標値、例えば実際に要した所要時間と推定値との差分や、消費電力量、単位時間あたりの温度変化など、もしくはそれらを複合した式を報酬とし、強化学習を実施するようにすれば良い。

送受信部１０は、室内機ＩＵ及び室外機ＯＵを含む外部の装置と各種情報の送受信を行うものである。特に、実施の形態１において、送受信部１０は、室内機及び室外機から空調データを受信し、室内機及び室外機に対して制御信号を送信する。
ここで、空調データとは、空調機ＡＣの制御に関する情報であり、実施の形態１において、室内機ＩＵから送信される室内環境情報、設定情報、運転情報、及び室外機ＯＵから送信される室外環境情報を含む。判定部１５に入力される空調データは、現在時刻の空調データであることが望ましいが、いつ計測された空調データか識別できるものであれば、空調機ＡＣの前回電源ＯＦＦ時から現在時刻までの過去の空調データであってもよい。また、送受信部１０は、環境値が目標値に到達するまでに実際に要した所要時間を室内機ＩＵから取得する。空調データに、当該所要時間を加えたデータを学習データと呼ぶ。推定部１６が推定に用いた空調データに上記の所要時間が加えられた学習データは記憶部１２に蓄えられ、機械学習部１４が行う機械学習モデルの学習に用いられる。

制御部１１は、複数の室内機ＩＵと複数の室外機ＯＵを集中管理するものである。制御部１１は、室外機ＯＵ及び室内機ＩＵから空調データを取得する。制御部１１は、複数の室外機ＯＵ及び複数の室内機ＩＵから取得した空調データを集約し、集約した空調データを機械学習部１４及び判定部１５に送信する。また、実施の形態１において、制御部１１は、後述する推定部１６の推定結果に基づいて空調機ＡＣの制御を行う。より具体的には、制御部１１は、推定部１６が推定した起動時刻に空調機ＡＣを起動する制御を行う。

記憶部１２は、各種情報を記憶するものであり、実施の形態１においては、後述する機械学習部１４で学習を行った機械学習モデルや、当該機械学習モデルの学習を行うための学習データ等を記憶する。

空調データ取得部１３は、送受信部１０を介して空調機ＡＣから空調データを取得するものである。実施の形態１においては、空調データ取得部１３は、現在時刻の空調データを取得し、取得した空調データを起動時刻の推定に用いるための空調データとして乖離度算出部１５１に送信する。また、空調データ取得部１３は、推定部１６の推定結果に基づき制御部１１が空調機ＡＣの制御を行った後、推定部１６の推定に用いた空調データに、環境値が目標値に到達するまでに要した所要時間を加え、機械学習部１４が機械学習モデルの学習を行うための学習データとして記憶部１２に送信する。

機械学習部１４は、空調データから空調機の起動時刻を推定するための機械学習モデルの学習を行うものであり、学習データ取得部１４１と学習部１４２とを備える。実施の形態１において、機械学習部１４が学習を行う機械学習モデルは、室内温度と目標温度との温度差、室内温度と室外温度との温度差、及び空調機の運転能力が入力され、空調機ＡＣを起動させてから目標時刻に部屋の室内温度が目標温度に到達するまでに要する所要時間を出力するものである。
また、学習部１４２が用いる機械学習モデルは、非線形関数を表現できるものであれば何でもよいが、実施の形態１においては、ニューラルネットワークを用いるものとする。

機械学習モデルへの入力は上記のものに限らず、差をとらずに、室内温度や目標温度及び室外温度を直接入力したり、室内湿度や室外湿度、室外機ＯＵや室内機ＩＵの台数を入力したりするようにしてもよい。また、機械学習モデルの出力も同様に、上記のものに限らず、その他の環境値、例えば、室内湿度等が目標値に到達するまでに要する所要時間を出力するようにしてもよいし、複数の環境値が目標値に到達するまでに要する所要時間を出力するようにしてもよい。

上記の機械学習モデルは、目標時刻や現在時刻も入力に用い、起動時刻を直接出力する構成としてもよい。実施の形態１において、後述する推定部１６は、学習済みの機械学習モデルを用いて環境値が目標値に到達するまでの所要時間を取得し、当該所要時間と目標時刻とから空調機ＡＣの起動時刻を推定するものであるが、空調機ＡＣの起動時刻は上述のように所要時間と目標時刻とから算出されるものであるため、目標時刻が設定された状態において、機械学習モデルが起動時刻を直接出力する場合も、機械学習モデルは所要時間を出力し、後処理を行うことにより起動時刻を算出する場合も、推定部１６としては、起動時刻を推定することとなる。

学習データ取得部１４１は、記憶部１２から上記機械学習モデルの学習を行うための学習データと、機械学習モデルとを取得するものである。また、学習データ取得部１４１は、取得した学習データに対して、機械学習モデルの入力に適した形式に変換する前処理を行う。実施の形態１において、機械学習モデルの入力は、室内温度と目標温度との温度差、室内温度と室外温度との温度差、及び空調機の運転能力であるので、学習データ取得部１４１は、空調データに含まれる室内環境情報が示す室内温度と設定情報が示す目標温度との温度差を算出し、同様に、室内温度と室外環境情報が示す室外温度との温度差を算出する。ここで、空調機の運転能力は運転情報に含まれるものとしたため、上記のような前処理は必要ない。また、以下では、前処理後のデータも学習データと呼ぶこととする。

学習部１４２は、学習データ取得部１４１が取得した学習データを用いて、機械学習モデルの学習を行うものである。学習は、誤差逆伝搬法等の既存の学習方法を用いて行えばよい。ここで、機械学習モデルの学習は初期学習だけでなく、再学習も含む。学習部１４２は、学習後の機械学習モデルを記憶部１２に送信し、記憶部１２に記憶させる。

判定部１５は、空調データ取得部１３が取得した空調データが疎領域に存在するか判定するものである。判定部１５は、予め設定された基準に基づいて、空調データが疎領域に存在するか判定するものであり、実施の形態１において、乖離度算出部１５１と、乖離度判定部１５２とを備える。
また、空調データが疎領域に存在するとは、空調データが学習データに対する外れ値である場合だけでなく、空調データの付近に学習データが存在するものの、その数が少ない場合も含む。

実施の形態１において、記憶部１２は、空調機ＡＣから受信した空調データに実際の所要時間を加えたものを学習データとして記憶し、機械学習部１４は、記憶部１２に追加された学習データを用いて機械学習モデルを更新する。そのため、学習データの疎領域は最初に決められてから変化しない静的なものではなく、学習データの増加に伴い動的に変化するものである。また、実施の形態１において、学習データの追加及び更新は一日毎に行うものとする。ただし、既に学習データが十分にあり、あまり頻繁な更新が必要でないような場合には、数日毎や一週間毎の更新としてもよい。

乖離度算出部１５１は、学習データに対する空調データの乖離度を算出するものである。ここで、学習データに対する空調データの乖離度とは、空調データと学習データに含まれる空調データの乖離具合を示すものであり、当該学習データを用いて学習した機械学習モデルに当該空調データを入力した場合に適切な推定結果を得ることができるかの指標となる量である。空調データと学習データの乖離具合は、空調データに含まれるパラメータで定義されるパラメータ空間で表される。ここで、空調データに含まれるパラメータとは、空調データが含む各種情報のパラメータであり、例えば、室内環境情報が示す室内温度や室外環境情報が示す室外温度等である。乖離度の具体例としては、パラメータ空間における空調データと学習データとの距離、学習データをデータ数に基づいて複数の領域に分類したときに、空調データがどの領域に存在するかを示す値、及び空調データの周囲にどれだけの学習データが存在するかを示す値等である。ここで、空調データと学習データとの距離とは、空調データと、学習データに含まれる各データとの距離だけでなく、学習データが多く集まる密領域との距離、及び学習データを確率密度関数で近似したときの確率密度関数の最頻値との距離等も含む。

乖離度判定部１５２は、乖離度算出部１５１が算出した乖離度に基づいて、空調データが疎領域に存在するか判定するものである。より具体的には、乖離度判定部１５２は、乖離度が所定の閾値以上か判定し、乖離度が所定の閾値以上である場合、空調データは疎領域に存在すると判定する。

以下で、乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例について説明する。

一つ目の方法は、起動時刻推定用の空調データと、学習データに含まれる空調データに対する最近傍点との距離を乖離度として算出する方法である。この方法について、図３を参照しながら説明する。図３は、乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。

図３において、黒丸で示される点が現在時刻の空調データＤ３１であり、白丸で示される点が学習データである。白丸で示される点は、正確には、学習データを空調データと同じパラメータ空間、すなわち所要時間を取り除いた空間に射影したデータであるが、ここでは、学習データを射影したデータも学習データと呼ぶこととする。
横軸と縦軸はともに、空調データに含まれるパラメータであり、例えば、室内温度や室外温度等である。空調データに含まれるパラメータは、一般にはより高次元であるが、図３においては説明のために２次元で表している。乖離度算出部１５１は、空調データＤ３１と、各学習データとの距離を算出する。ここで、学習データのうち、空調データＤ３１と最も近い点は、学習データＤ３２であり、乖離度算出部１５１は、空調データＤ３１と学習データＤ３２との距離Ｌ３１を乖離度として算出する。

乖離度判定部１５２は、乖離度算出部１５１が算出した乖離度、すなわち、空調データＤ３１と学習データＤ３２との距離Ｌ３１が所定の第一閾値Ｔ１以上か判定し、第一閾値Ｔ１以上である場合、空調データＤ３１は疎領域に存在すると判定する。ここで、第一閾値Ｔ１はゼロ以上の実数である。

二つ目の方法は、学習データの密領域を設定し、空調データと密領域との距離を算出する方法であり、この方法について、図４を参照しながら説明する。図４は、乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。

乖離度算出部１５１は、学習データが多数存在する領域である密領域Ａ４１を設定する。ここで、密領域Ａ４１の設定方法としては、例えば、適当な大きさの矩形領域を取り、その矩形領域の内側に存在する学習データの個数をカウントし、何割の学習データが矩形領域に含まれるか算出する。そして、矩形領域に含まれる学習データが所定の割合未満であれば、矩形領域を広げ、所定の割合以上であれば矩形領域を狭くする処理を繰り返す。この処理において、前回の処理では矩形領域を広げ、今回の処理では矩形領域を狭くすることになった場合、前回の処理後の矩形領域を密領域Ａ４１として設定すればよい。前回の処理で矩形領域を狭くし、今回の処理で矩形領域を広げた場合も同様である。
あるいは、空調データが外れ値かどうかを判定したい場合には、学習データに含まれる各データのパラメータの最大値及び最小値を検出し、当該最大値及び最小値により形成される領域を密領域Ａ４１として設定するようにしてもよい。図４においては、後者の方法を採用し、パラメータ１の最大値は学習データＤ４２、パラメータ１の最小値は学習データＤ４３、パラメータ２の最大値は学習データＤ４４、パラメータ２の最小値は学習データＤ４５により設定される。図４において、乖離度算出部１５１は、空調データＤ４１と密領域Ａ４１との距離Ｌ４１を乖離度として算出する。

乖離度判定部１５２は、乖離度算出部１５１が算出した乖離度、すなわち、空調データＤ４１と密領域Ａ４１との距離Ｌ４１が所定の第二閾値Ｔ２以上か判定し、第二閾値Ｔ２以上である場合、空調データＤ４１は疎領域に存在すると判定する。第二閾値Ｔ２はゼロ以上の実数である。

三つ目の方法は、学習データを確率密度関数により近似し、当該確率密度関数に基づいて、乖離度を算出するものである。この方法について、図５を参照しながら説明する。図５は、乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。

乖離度算出部１５１は、図５で示すように、学習データの出現頻度を正規分布ＤＦ５１で近似し、正規分布ＤＦ５１の最頻値からの距離を乖離度として算出する。ここで、正規分布ＤＦ５１の最頻値をμとする。
乖離度判定部１５２は、乖離度算出部１５１が算出した乖離度が所定の第三閾値Ｔ３以上か判定し、第三閾値Ｔ３以上である場合、空調データＤ５１は疎領域に存在すると判定する。第三閾値Ｔ３はゼロ以上の実数であり、図５においては、σで示されている。

また、図５においては、パラメータが一次元の場合について示したが、パラメータが二次元以上の場合も同様である。
また、単に最頻値からの距離を乖離度として算出するのではなく、確率密度関数に基づき、パラメータ空間をいくつかの領域に分類し、どの領域に存在するかを乖離度として算出するようにしてもよい。この方法について、図６を参照しながら説明する。
まず、乖離度算出部１５１は、学習データを確率密度関数で近似する。ここでも、確率密度関数は正規分布を採用する。そして、当該確率密度関数に基づき、例えば、５０％領域Ａ６１、７５％領域Ａ６２、９５％領域Ａ６３、９５％領域外Ａ６４を設定する。ここで、５０％Ａ６１領域とは、その領域の内部に学習データの５０％が存在する領域であり、７５％領域Ａ６２、９５％領域Ａ６３についても同様である。そして９５％領域外Ａ６４は、９５％領域Ａ６３の外側の領域であり、この領域を疎領域として設定する。乖離度算出部１５１は、空調データの位置に基づき、例えば、空調データが５０％領域Ａ６１に存在するときは５０、７５％領域Ａ６２に存在するときは７５、９５％領域Ａ６３に存在するときは９５、９５％領域外に存在する場合には１００を乖離度として算出する。そして、９６以上１００未満の値を所定の閾値として設定することにより、乖離度判定部１５２は空調データが疎領域に存在するか判定することができる。

四つ目の方法は、空調データの周囲にどれだけの学習データが存在するかを示す値を乖離度として算出する方法である。この方法について図７を参照しながら説明する。図７は、乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。

乖離度算出部１５１は、空調データＤ７１の周囲に領域Ａ７１を設定する。図７において、領域Ａ７１は空調データＤ７１を中心とした半径Ｒ７１以下の領域として設定された円盤である。ここで、半径Ｒ７１の値は、空調制御装置１００の設計者、あるいは空調機ＡＣのユーザーにより予め設定されているものとする。乖離度算出部１５１は、まず、この領域Ａ７１の内部に存在する学習データの数をカウントし、学習データの数、あるいは領域Ａ７１の内部に存在する学習データ数の割合を算出する。上記により算出された値は、その値が小さいほど疎領域に存在する可能性が高くなるようになっているため、その値に負号を付けるか、逆数を取ったものを乖離度として算出する。この操作を行うことにより、四つ目の方法においても、乖離度として算出された値は、その値が大きいほど疎領域に存在する可能性が高くなるものとして扱うことができる。乖離度判定部１５２は、乖離度算出部１５１が算出した乖離度が所定の第四閾値Ｔ４以上か判定し、第四閾値Ｔ４以上である場合、空調データＤ７１は疎領域に存在すると判定する。例えば、図７においては、空調データＤ７１の周囲に設定された領域Ａ７１には学習データが存在しないため、空調データＤ７１は疎領域に存在すると判定され、空調データＤ７２の周囲に設定された領域Ａ７２には学習データが多数存在するため、空調データＤ７２は疎領域に存在しない、あるいは密領域に存在すると判定される。ここで領域Ａ７２は、空調データＤ７２を中心とした半径Ｒ７１以下の領域として設定された円盤である。

上記の四つの方法において、空調データのパラメータ空間上の位置に応じて乖離度を算出する構成としたが、空調データのパラメータ空間上の位置座標、すなわちパラメータそのものを直接乖離度として用いる構成としてもよい。この場合において、乖離度判定部１５２は、予め疎領域を設定しておき、空調データの各パラメータが疎領域の範囲内に含まれるか否かを判定することにより、空調データが疎領域に存在するかを判定することができる。例えば、密領域が矩形領域で設定され、疎領域が当該密領域の外側の領域として設定されている場合には、各パラメータが所定の閾値以上かをそれぞれ独立に判定すればよい。また、乖離度判定部１５２が設定する疎領域は、最初に設定してから固定する必要はなく、学習データの増加に伴い動的に設定するようにしてもよい。
また、上記の方法において用いる位置座標は、空調データに含まれるパラメータそのものである必要はなく、適当な座標変換を行った後の座標成分を用いるようにしてもよい。

推定部１６は、機械学習部１４が学習を行った学習済みの機械学習モデルを用いて、空調データ取得部１３が取得した空調データから空調機の起動時刻を推定するものである。実施の形態１において、機械学習モデルは、環境値が目標値に到達するまでに要する所要時間を出力するものであるので、推定部１６は、目標時刻から所要時間前の時刻を算出することにより、空調機ＡＣの起動時刻を推定する。
また、推定部１６は、判定部１５が、空調データが疎領域に存在すると判定した場合、学習データを用いて空調データから空調機の起動時刻を推定するための学習を行った機械学習モデルに対して、非学習モデルによる外挿を行い、非学習モデルを用いて空調データから空調機の起動時刻を推定する。
ここで、非学習モデルによる外挿とは密領域が一つの区間で表されるときに、その区間の外側の領域において、空調データと起動時刻との関係を非学習モデルで表現するだけでなく、密領域が複数の区間として離散的に存在する場合に、その複数の区間の間の領域において、空調データと起動時刻との関係を非学習モデルで表現することも含む。
また、実施の形態１において、推定部１６は、非学習モデルとして線形モデルを用いて機械学習モデルに対する外挿を行う。ここでは、線形モデルを用いて外挿を行うものとしたが、ある領域において、経験則的に、または一般的に知られている関係式から空調データと所要時間との関係がどのようなものとなるのか既知の場合には、線形モデル以外のモデルを用いて外挿を行うようにしてもよい。

推定部１６は、学習データ取得部１４１と同様に、取得した空調データに対して、機械学習モデルの入力に適した形式に変換する前処理を行う。実施の形態１において、機械学習モデルの入力は、室内温度と目標温度との温度差、室内温度と室外温度との温度差、及び空調機の運転能力であるので、推定部１６は、空調データに含まれる室内環境情報が示す室内温度と設定情報が示す目標温度との温度差、及び室内温度と室外環境情報が示す室外温度との温度差を算出する。ここで、空調機の運転能力は運転情報に含まれるものとしたため、上記のような前処理は必要ない。また、学習データと同様に、前処理後の空調データも空調データと呼ぶこととする。

推定部１６が、機械学習モデルに対して非学習モデルによる外挿を行う方法の具体例について図８を参照しながら説明する。図８は、推定部１６が学習済みの機械学習モデルに対して非学習モデルによる外挿を行い、当該非学習モデルを用いて空調データから空調機の起動時刻を推定する方法の具体例を示す説明図である。

図８に実線で示されるのが機械学習モデルＭ８１であり、非線形関数で表される。図８において、横軸は室内温度と目標温度との温度差、縦軸は室内温度が目標温度に到達するまでの所要時間であるが、一般には、横軸は機械学習モデルＭ８１の入力パラメータ、縦軸は機械学習モデルＭ８１の推定値である。また、図８は説明のために二次元で示しているが、入力パラメータ及び出力パラメータは、一般にはより高次元である。例えば、実施の形態１においては、機械学習モデルへの入力は、室内温度と目標温度との温度差、室内温度と室外温度との温度差、空調機ＡＣの運転能力の３つであるので、入力パラメータは三次元であり、出力は室内温度が目標温度に到達するまでの所要時間であるので、出力パラメータは一次元である。

ここで、入力された空調データに推定値を加えたデータを推定データと呼ぶこととする。図８において、横軸上に黒四角で示されているのが空調データであり、実線の白丸が学習データ、点線の白丸が機械学習モデルＭ８１を用いて得られた推定値を含む推定データ、黒丸が入力された空調データから非学習モデルＭ８２を用いて得られた推定値を含む推定データである。
また、図８において、室内温度と目標温度との温度差がＰ８２以上の領域が密領域、Ｐ８２未満の領域が疎領域として設定されている。

判定部１５で空調データが密領域に存在すると判定された場合は、単に当該空調データのパラメータ値に対応する推定値を機械学習モデルＭ８１から算出するだけである。以下では、判定部１５で空調データが疎領域に存在すると判定された場合について説明する。

まず、推定部１６は、密領域に存在する空調データを２点選択する。ここで、密領域に存在する空調データを選択するとは、密領域と設定されているパラメータ空間の領域において、適当なパラメータを選択することを意味し、記憶部１２に記憶された学習データに含まれる空調データを２点選択するだけでなく、密領域において適当なパラメータを選択し、当該パラメータを有する空調データを擬似的に生成することを含む。また、ここでのパラメータ空間は、機械学習モデルの入力パラメータに対応するものであり、判定部１５における空調データのパラメータ空間とは一般には異なるが、機械学習モデルの入力パラメータは空調データのパラメータから生成されるものであり、機械学習モデルのパラメータ空間の場合においても、疎領域及び密領域の設定は空調データのときと同様に行えばよい。

ここで、選択される２点は任意であるが、疎領域の境界付近に存在することが望ましい。図８においては、一点は密領域と疎領域の境界に存在する点（空調データＤ８２）を選択している。もう一点は、空調データＤ８２の近くで選択することが望ましく、図８においては、空調データＤ８２が示す温度差Ｐ８２より所定の値ＰＤ８１だけ大きい値Ｐ８３を示す空調データＤ８３を選択している。

次に、推定部１６は、選択した２点の空調データに対応する推定データを機械学習モデルＭ８１を用いて生成する。図８において、機械学習モデルＭ８１を用いると、空調データＤ８２からは推定値ＥＶ８２が得られ、推定データＥ８２が生成される。また、空調データＤ８３からは推定値ＥＶ８３が得られ、推定データＥ８３が生成される。
そして、推定部１６は生成された推定データＥ８２と推定データＥ８３を通る直線を非学習モデルＭ８２として生成する。図８において、非学習モデルＭ８２は一点鎖線で示される。

最後に、推定部１６は、入力された推定対象である空調データＤ８１の推定値ＥＶ１を推定する。ここで、空調データＤ８１に推定値ＥＶ１を追加した推定データＥ１を記憶部１２に記憶するようにしてもよい。

また、上記において、推定部１６は推定データを生成したが、機械学習モデルあるいは非学習モデルから空調データに対応する推定値を算出する動作のみ行い、推定データは生成しない構成であってもよい。

図８において、非学習モデルＭ８２は直線で示されているが、一般には非学習モデルは超平面である。上記においては、非学習モデルＭ８２を生成するために、空調データを２点選択したが、より高次元な場合には、より多く、具体的には超平面の次元より１だけ多く、空調データを選択する必要がある。
また、上記においては、機械学習モデルＭ８１上の２点を用いて、非学習モデルＭ８２を生成したが、非学習モデルの生成方法はこれに限らず、例えば、密領域において、疎領域との境界付近に存在する学習データを複数個選択し、これらの学習データから回帰分析により得られる超平面を非学習モデルとしてもよい。

非学習モデルによる外挿を行った場合と、外挿を行わなかった場合の違いについて説明する。例えば、図８に示すように、疎領域においては、パラメータの増加に従い、単調増加すべき状況であっても、学習データが少ないと適切に学習を行うことができず、減少後増加するような機械学習モデルが生成されてしまう可能性がある。このような場合において、非学習モデルによる外挿を行わない場合、空調データＤ８１から得られる推定値は推定値ＥＶ２（得られる推定データは推定データＥ２）となり、本来予想される推定値より大きな値を推定してしまう。しかしながら、実施の形態１に係る空調制御装置１００では、密領域の機械学習モデル、すなわち密領域の学習データに基づき疎領域に非学習モデルを生成し、密領域の機械学習モデルと接続することにより、非学習モデルを外挿するので、より適切な値、図８においては、推定値ＥＶ１であるが、外挿を行わない場合よりも小さな値を推定することができる。

次に、実施の形態１における空調制御装置１００のハードウェア構成について説明する。空調制御装置１００の各機能は、コンピュータにより実現される。図９は、空調制御装置１００を実現するコンピュータのハードウェア構成の例を示す構成図である。
図９に示したハードウェアには、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）等の処理装置１００００と、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）やハードディスク等の記憶装置１０００１が備えられる。
図２に示す、送受信部１０、制御部１１、空調データ取得部１３、機械学習部１４、判定部１５、及び推定部１６は、記憶装置１０００１に記憶されたプログラムが処理装置１００００で実行されることにより実現され、記憶部１２は記憶装置１０００１により実現される。ここで、上記の構成は、単数の処理装置１００００及び記憶装置１０００１により実現する構成に限らず、複数の処理装置１００００及び記憶装置１０００１により実現する構成であってもよい。
また、空調制御装置１００の各機能を実現する方法は、上記したハードウェアとプログラムの組み合わせに限らず、処理装置にプログラムをインプリメントしたＬＳＩ（Ｌａｒｇｅ Ｓｃａｌｅ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）のような、ハードウェア単体で実現するようにしてもよいし、一部の機能を専用のハードウェアで実現し、一部を処理装置とプログラムの組み合わせで実現するようにしてもよい。

以上のように、空調制御装置１００は構成される。

次に、空調制御装置１００の動作について学習動作と推定動作に分けて説明する。ここで、空調制御装置１００の動作が空調制御方法であり、学習動作に対応する空調制御方法を学習方法、推定動作に対応する方法を推定方法と呼ぶこととする。また、空調制御方法をコンピュータに実行させるプログラムが空調制御プログラムである。
まず、空調制御装置１００の学習動作について、図１０を参照しながら説明する。図１０は、空調制御装置１００が機械学習モデルの学習を行う動作を示すフローチャートである。

まず、ステップＳ１において、学習データ取得部１４１は、記憶部１２から学習データ及び機械学習モデルを取得する。学習データ取得部１４１は、取得した学習データと機械学習モデルとを学習部１４２に送信する。ここで、初期学習の場合は、全ての学習データを取得し、再学習の場合には、追加された学習データのみを取得すればよい。
また、学習データ取得部１４１は、取得した学習データに対して、機械学習モデルの入力に適した形式に変換する前処理を行う。実施の形態１において、機械学習モデルの入力は、室内温度と目標温度との温度差、室内温度と室外温度との温度差、及び空調機の運転能力であるので、学習データ取得部１４１は、空調データに含まれる室内環境情報が示す室内温度と設定情報が示す目標温度との温度差を算出し、同様に、室内温度と室外環境情報が示す室外温度との温度差を算出する。

次に、ステップＳ２において、学習部１４２は、受信した学習データを用いて機械学習モデルの学習を行う。実施の形態１においては、機械学習モデルはニューラルネットワークであり、学習は誤差逆伝搬法等の既存の学習方法を用いて行うことができる。

最後に、ステップＳ３において、学習部１４２は、学習済みの機械学習モデルを記憶部１２に保存する。
以上のステップを行い、空調制御装置１００は学習動作を終了する。

次に、空調制御装置１００の推定動作について、図１０を参照しながら説明する。
図１０は、空調制御装置１００が空調機ＡＣの起動時刻を推定する動作を示すフローチャートである。

まず、ステップＳ１０１の空調データ取得工程において、空調データ取得部１３は送受信部１０を介して空調機ＡＣから空調データを取得する。より具体的には、空調データ取得部１３は、室外機ＯＵから室外環境を取得し、室内機ＩＵから室内環境情報、設定情報及び運転情報を取得する。

次に、ステップＳ１０２の判定工程において、判定部１５は、ステップＳ１０１で取得した空調データが疎領域に存在するか判定する。
より具体的には、まず、乖離度算出部１５１が、当該空調データの乖離度を算出する。ここでは、上記の一つ目の方法において、空調データと学習データとの距離を乖離度として算出したものとする。
そして、乖離度判定部１５２が、当該乖離度が所定の閾値以上か判定し、閾値以上であると判定した場合、当該空調データは疎領域に存在すると判定する。

ステップＳ１０２の判定工程において、空調データが疎領域に存在すると判定された場合には、空調制御装置１００の動作はステップＳ１０３に進み、空調データが疎領域に存在しないと判定された場合には、空調制御装置１００の動作はステップＳ１０４に進む。

以下のステップＳ１０３とステップＳ１０４とを合わせて推定工程と呼ぶこととする。
ステップＳ１０３において、推定部１６は、空調機の起動時刻の推定を行う。ステップＳ１０３は、空調データが疎領域に存在すると判定された場合なので、機械学習モデルに対して非学習モデルで外挿を行い、当該外挿モデルを用いて起動時刻を推定する。
ステップＳ１０４においても、推定部１６は、空調機の起動時刻の推定を行うが、ステップＳ１０４は、空調データが疎領域に存在しないと判定された場合なので、非学習モデルによる外挿は行わず、機械学習モデルを直接用いて起動時刻を推定する。

次に、ステップＳ１０５の制御工程において、制御部１１は、推定部１６が推定した起動時刻に空調機ＡＣが起動するように制御信号を送信する。制御信号は送受信部を介して、空調機ＡＣに送信され、空調機ＡＣは当該制御信号を受信すると、当該起動時刻において起動するように設定される。

最後に、ステップＳ１０５の学習データ保存工程において、送受信部１０は、室内機ＩＵから実際に室内温度が目標温度に到達するまでにかかった所要時間を取得し、当該所要時間をステップＳ１０１で取得した空調データと合わせて学習データとして記憶部１２に保存する。

以上のような空調制御装置１００の推定動作により、入力された空調データが学習データから乖離した値であっても、学習済みの機械学習モデルに非学習モデルを外挿することにより、外挿を行わずに機械学習モデルをそのまま用いるよりも、より適切に起動時刻を推定することができる。
すなわち、実施の形態１に係る空調制御装置１００は、学習データが少ない領域では非学習モデルによる外挿を行い、空調機の起動時刻を推定することにより、不適切な推定をする可能性を低減することができる。

実施の形態１における空調制御装置１００の変形例について、以下で説明する。

上記において、制御部１１は、空調制御装置１００が備える他の構成と同一のハードウェアにより実現される構成としたが、他の構成と別のハードウェア（処理装置及び記憶装置）により実現される構成としてもよい。ここで、制御部１１を別のハードウェアで実現する構成において、制御部１１は、同一のハードウェアで実現されるときと同様、複数の室内機ＩＵ及び室外機ＯＵを集中管理する構成でも、個々の室内機ＩＵ及び室外機ＯＵにそれぞれ備えられ、それぞれに備えられた制御部１１が室内機ＩＵ及び室外機ＯＵを制御する構成としてもよい。
また、制御部１１だけでなく、空調制御装置１００は、図１に示すように、室内機ＩＵ及び室外機ＯＵの外部に備える構成だけでなく、室内機ＩＵ及び室外機ＯＵの内部に備える構成であっても良い。

上記において、空調制御装置１００は、複数の室外機ＯＵに接続され、室外機ＯＵを介して室内機ＩＵと接続される構成としたが、空調制御装置１００と室内機ＩＵ及び室外機ＯＵとがそれぞれ直接接続される構成や、空調制御装置１００と室内機ＩＵが接続され、室内機ＩＵを介して空調制御装置１００と室外機ＯＵが接続される構成であってもよい。また、複数の室外機ＯＵ及び複数の室内機ＩＵを備える構成としたが、室外機ＯＵ及び室内機ＩＵは単数であっても良い。

推定部１６及び学習データ取得部１４１で空調データを機械学習モデルの入力に適した形に変換する前処理を行う構成としたが、推定部１６ではなく判定部１５で上記の前処理を行う構成としても良い。
あるいは、推定部１６及び学習データ取得部１４１で前処理を行うのではなく、空調データ取得部１３が前処理を行い、推定に用いる空調データだけでなく、学習データについても前処理を行ったデータを記憶部１２に保存する構成としてもよい。
これらの場合において、判定部１５における空調データが疎領域に存在するかの判定は、前処理後の空調データで行うようにしてもよい。すなわち、判定部１５は空調データに含まれるパラメータではなく、機械学習モデルの入力パラメータにより張られるパラメータ空間において、空調データが疎領域に存在するか判定するようにしてもよい。

実施の形態１において、空調制御装置１００の機械学習モデルの学習に用いる学習データ群は、密領域と疎領域を含むことを前提としており、学習初期の段階において、学習データ群全体のデータ数が十分な量存在せず、密領域と疎領域を明確に区別できないようなときには、十分に効果を発揮できない場合が考えられる。そのため、判定部１５は、空調データの乖離度を算出する前に、学習データ群に含まれる学習データの数をカウントし、学習データの数が所定の第五閾値以上の場合に、空調データが疎領域に存在するか判定を行い、推定部１６は、非学習モデルの外挿を行うようにしてもよい。このようにすることで、学習データ群全体のデータ数が少ないときには、学習データを蓄積することに専念し、学習データが十分蓄積できた時点から、非学習モデルによる外挿を行うようにすることができる。

判定部１５は、空調データが疎領域に存在するかを判定したが、疎領域に存在するかだけでなく、どのパラメータが学習データから乖離した値、すなわち異常な値を示しているかを判定するようにしてもよい。例えば、空調データと学習データの各パラメータを比較したり、空調データと予め設定された疎領域のパラメータを比較したりすることにより実現できる。
例えば、図３においては、空調データＤ３１のパラメータ１は、学習データＤ３３のパラメータ１以下であり、異常な値は取っておらず、一方、空調データのパラメータ２は、どの学習データよりも大きな値を取っており、学習データから乖離している。乖離度判定部１５２は、空調データと学習データの各パラメータを比較することにより、どのパラメータが異常な値を取っているか判定することができる。
また、パラメータの判定は、空調データが疎領域に存在すると判定された場合において追加で判定するようにしてもよいし、疎領域の判定と同時に行う構成であってもよい。例えば、異常なパラメータが見つかった場合に、当該パラメータが異常であり空調データは疎領域に存在すると判定するようにしてもよい。

実施の形態１に係る空調制御装置１００は、空調データが疎領域に存在する場合、機械学習モデルに外挿を行い、起動時刻を推定するものであるが、推定された起動時刻に起動しても外乱等により目標時刻に目標温度に到達できない場合が想定される。また、会社や学校等において、長期休暇明けのような場合には、壁の潜熱等の影響により、室内温度の目標温度到達時刻は、通常より遅くなる場合がある。上記のような課題に対しては、空調機ＡＣのユーザーにとって、目標時刻より遅い時刻に目標温度に到達するよりも、目標時刻より若干早い時刻に目標温度に到達した方が、快適性が高いと考えられるため、推定部１６は機械学習モデルにより得られた推定値、あるいは非学習モデルにより得られた推定値に、より早い起動時刻を推定するための補正を行うようにしてもよい。
例えば、推定部１６は、推定値に安全率を加える補正を行うことにより、空調機ＡＣの推定される起動時刻が安全率を加えない場合よりも早くなるように推定を行う。より具体的には、所要時間の推定値に安全率を加えることにより、より長い所要時間を推定値として算出したり、所要時間の推定値に１より大きな値を積算し、より長い所要時間を推定値として算出したりすることにより、補正を行わない場合よりも、より早い時間を起動時刻として推定し、結果として補正を行わない場合よりもより早い時刻に目標時刻に到達するように空調機を制御することができる。

実施の形態２．
実施の形態１において、機械学習モデルは、室内温度と目標温度との温度差、室内温度と室外温度との温度差、及び空調機の運転能力を入力とするものであった。実施の形態２における空調制御装置１００は、機械学習モデルの入力として、さらに空調機ＡＣの能力セーブ値を用いるものである。ここで、能力セーブ値は運転情報に含まれるものとする。能力セーブとは、空調機の運転能力を抑えることにより最大電力や消費電力を抑えるものであり、例えば、圧縮機周波数を低くしたり、ファンの回転数を低下したりすることにより達せられる。また、能力セーブ値とは、能力セーブの尺度であり、実施の形態２においては、空調機ＡＣの最大運転能力に対して実際に運転させる運転能力の比率を能力セーブ値として用いる。実施の形態２において、能力セーブ値は、６０％、８０％、及び１００％に設定可能であるとする。ここで、６０％は６０％の能力で空調機ＡＣを運転している状態、８０％は８０％で空調機ＡＣを運転している状態、１００％は１００％で空調機ＡＣを運転している状態を示す。
また、実施の形態２においては、空調機ＡＣの最大運転能力に対して実際に運転させる運転能力の比率を能力セーブ値として用いるものとしたが、圧縮機周波数やファンの回転数を直接能力セーブ値として用いるようにしてもよい。
実施の形態１と異なる部分について、以下で説明する。

空調データ取得部１３は、取得した空調データに能力セーブ値を加えることにより、空調機ＡＣの能力セーブ値を含む運転情報を取得する。より具体的には、取得した空調データに能力セーブ値が６０％であるという情報を追加したデータ、取得した空調データに能力セーブ値が８０％であるという情報を追加したデータ、及び取得した空調データに能力セーブ値が１００％であるという情報を追加したデータを生成する。ここで、空調データに能力セーブ値に関する情報が追加されたデータも空調データと呼ぶこととする。空調データ取得部１３は、生成した各空調データを判定部１５に送信する。

機械学習部１４は、実施の形態１と同様に、機械学習モデルの学習を行う。ここで、実施の形態２における機械学習モデルは、上記したように、入力パラメータの一つに能力セーブ値を有する。このときの学習データは、後述する推定部１６により選択された能力セーブ値を含む空調データと、実際に目標温度に到達するまでにかかった所要時間を加えたものである。

判定部１５は、取得した空調データに能力セーブ値を加えたデータが疎領域に存在するか判定する。空調データが疎領域に存在するか判定する方法は、実施の形態１に記載された方法をそのまま用いても良いが、実施の形態２における能力セーブ値のようにパラメータが離散的な場合に用いることができる方法について、図１２を参照しながら説明する。

図１２は、乖離度算出部１５１が乖離度を算出する方法、及び乖離度判定部１５２が、空調データが疎領域に存在するか判定する方法の具体例を示す説明図である。
図１２において、空調データＤ１２１は黒丸、学習データは白丸で示されている。
実施の形態２において、能力セーブ値は離散値としているので、能力セーブ値の値により各グループを形成することができる。乖離度算出部１５１は、入力された空調データが属するグループに含まれる学習データのデータ数をカウントし、そのカウント値に負号を付したものを乖離度として算出する。図１１において、空調データＤ１２１は、能力セーブ値が６０％のグループＧ１２１に属しており、そのグループ内の学習データは１つなので、乖離度を－１と算出する。ここで、乖離度算出部１５１は、データ数のカウント値に負号を付したものを乖離度として算出したが、負号を付すのではなく、逆数を取ったものを乖離度として算出するようにしてもよい。

乖離度判定部１５２は、算出された乖離度が所定の第六閾値以上か判定し、第六閾値以上の場合、空調データは疎領域に存在すると判定する。図１２において、例えば、第六閾値が－３と設定されている場合には、空調データＤ１２１は疎領域に存在すると判定される。上記の判定において、疎領域に存在しないと判定された場合、実施の形態１に記載した方法を用いて、他のパラメータに関して疎領域に存在するか判定する。

推定部１６は、能力セーブ値が加えられた各空調データを学習済みの機械学習モデルに入力し、起動時刻を推定する。すなわち、推定部１６は、能力セーブ値が６０％、８０％、及び１００％のときそれぞれに対する起動時刻を得る。

また、推定部１６は、どの能力セーブ値で空調機ＡＣを起動させるか、すなわち、どの能力セーブ値に対応する起動時刻で空調機ＡＣを起動させるかを選択する。ここで、ユーザーにより予め、空調機ＡＣの動作時における能力セーブ値が設定されている場合には、推定部１６は、単に設定された能力セーブ値に対応する起動時刻で起動させることを選択すればよい。あるいは、設定された能力セーブ値に対応する起動時刻が、現在時刻よりも前の時刻の場合には、対応する起動時刻が現在時刻よりも後の時刻で、かつ、最も大きな値を有する能力セーブ値を選択するようにしてもよい。

また、能力セーブ値及び能力セーブ値に対応する起動時刻の別の決定方法として、能力セーブ値と対応する起動時刻の組に基づいて、どの能力セーブ値及び起動時刻で空調機ＡＣを起動させるか決定するようにしてもよい。
例えば、あらかじめ能力セーブ値に対する単位消費電力等の指標地のテーブルや関数などを保持しておき、指標値と起動時刻との積が最小かつ起動時刻が現在時刻より後の時刻となる能力セーブ値と起動時刻の組を選択するようにしてもよい。

また、起動後の一定時間経過後も起動時刻推定を行い、現在の能力セーブ値から推定される起動時刻が現在時刻よりも早い時間であれば能力セーブ値を小さくするように制御し、逆に現在の能力セーブ値から推定される起動時刻が現在時刻よりも遅い時間であれば能力セーブ値を大きくするように制御してもよい。上記の制御を行うことにより、例えば、冬期に能力セーブ値を下げた状態で空調機ＡＣを運転させることによる不暖（居室が暖まらないこと）の発生を防ぐことができる。また、一定時間経過後の推定は一度だけでなく、複数回に渡って行うようにしてもよい。

推定部１６における非学習モデルの外挿についても、実施の形態１と同様の方法を用いても良いが、能力セーブ値のようにパラメータが離散的な場合に用いることができる外挿方法について、以下で説明する。
パラメータが離散的な場合、実施の形態１のように密領域の空調データを２点選択する方法は行わず、単に入力された空調データのパラメータと密領域の空調データのパラメータとの比を、離散的なパラメータを密領域のパラメータに置き換えて機械学習モデルから得られた推定値に、乗算した値を非学習モデルの外挿により得られた推定値としてもよい。例えば、図１２の場合には、密領域として、能力セーブ値が１００％のグループＧ１２３を用いると、空調データＤ１１１のパラメータは６０％であり、密領域に属する空調データのパラメータは１００％であるので、空調データＤ１１１の能力セーブ値を１００％に置き換えて機械学習モデルに入力し得られた推定値に、１００（％）／６０（％）を乗算した値を非学習モデルの外挿により得られた推定値としてもよい。
ここで、能力セーブ値と所要時間との関係は、能力セーブ値が大きいほど所要時間は短くなるため、密領域に属する空調データの能力セーブ値／入力された空調データの能力セーブ値を乗算したが、各パラメータと推定値がどのような比例関係にあるかに応じて、密領域に属する空調データと入力された空調データのどちらを分母、あるいは分子にするかを設定する必要がある。

推定部１６は、起動時刻の推定を行った後、制御部１１に推定結果と、どの空調データが選択されたかを通知し、制御部１１は、通知された起動時刻において、選択された能力セーブ値で空調機ＡＣを動作させるように制御を行う。
また、推定部１６は、空調データ取得部１３にどの空調データが選択されたかを通知し、空調データ取得部１３は選択された空調データと、目標温度に到達するまでの実際の所要時間とを合わせて学習データとして、記憶部１２に保存する。

実施の形態２に係る空調制御装置１００は以上のように構成され、適切な能力セーブ値で空調機ＡＣの運転を制御することにより、消費電力を低下することができる。

実施の形態１に係る空調制御装置１００の変形例については、実施の形態２に係る空調制御装置１００にも適用可能である。

本発明に係る空調制御装置は、目標時刻に目標温度に達するように空調を管理する空調システムに用いるのに適している。

１００ 空調制御装置、１０００ 空調システム、１０ 送受信部、１１ 制御部、１２ 記憶部、１３ 空調データ取得部、１４ 機械学習部、１４１ 学習データ取得部、１４２ 学習部、１５ 判定部、１５１ 乖離度算出部、１５２ 乖離度判定部、１６ 推定部

Claims (9)

1. 密領域と前記密領域より学習データが少ない疎領域とを含む学習データ群を用いて、空調機の制御に関する情報である空調データから前記空調機の起動時刻を推定するための学習を行った機械学習モデルに基づいて、前記空調機の起動時刻を推定する空調制御装置において、
   前記空調データを取得する空調データ取得部と、
   前記空調データが、前記疎領域に存在するか判定する判定部と、
   前記判定部が、前記空調データが前記疎領域に存在すると判定した場合、前記機械学習モデルに対して、前記空調データと前記起動時刻とを関連付ける非学習モデルによる外挿を行い、前記非学習モデルを用いて前記空調データから前記空調機の起動時刻を推定する推定部と、
   前記推定部が推定した起動時刻に前記空調機を起動する制御を行う制御部と、
   を備えた空調制御装置。
2. 前記空調データ取得部は、前記空調データとして、室内の環境を示す室内環境情報と、室外の環境を示す室外環境情報と、前記空調機の運転状態を示す運転情報と、を取得する
   ことを特徴とする請求項１に記載の空調制御装置。
3. 前記空調データ取得部は、前記空調機の能力セーブ値を含む前記運転情報を取得する
   ことを特徴とする請求項２に記載の空調制御装置。
4. 前記判定部は、
   前記学習データに対する前記空調データの乖離度を算出する乖離度算出部と、
   前記乖離度に基づいて、前記空調データが前記疎領域に存在するか判定する乖離度判定部と、
   を備えたことを特徴とする請求項１から３のいずれか一項に記載の空調制御装置。
5. 前記乖離度算出部は、前記空調データに含まれるパラメータで定義される空間における前記空調データと前記学習データとの距離を前記乖離度として算出し、
   前記乖離度判定部は、前記前記乖離度が所定の閾値以上か判定し、前記乖離度が前記閾値以上である場合、前記空調データは前記疎領域に存在すると判定する
   ことを特徴とする請求項４に記載の空調制御装置。
6. 前記推定部は、前記機械学習モデルの出力に対して、より早い起動時刻を推定するための補正を行う
   ことを特徴とする請求項１から５のいずれか一項に記載の空調制御装置。
7. 室内機と、
   前記室内機に接続され、前記空調機を構成する室外機と、
   前記室内機と前記室外機とを制御する、請求項１から６のいずれかを１項に記載の空調制御装置と、
   を備えた空調システム。
8. 空調機の制御に関する情報である空調データを取得する空調データ取得工程と、
   前記空調データが学習データの少ない疎領域に存在するか判定する判定工程と、
   前記判定工程で、前記空調データが前記疎領域に存在すると判定した場合、前記学習データを用いて、前記空調データから前記空調機の起動時刻を推定するための学習を行った機械学習モデルに対して、非学習モデルによる外挿を行い、前記非学習モデルを用いて前記空調データから前記空調機の起動時刻を推定する推定工程と、
   前記推定工程で推定した起動時刻に前記空調機を起動する制御を行う制御工程と、
   を含む空調制御方法。
9. 請求項８に記載の全工程をコンピュータに実行させる空調制御プログラム。

Images