JP6589227B1

JP6589227B1 - 制御装置、空調制御システム、制御情報の算出方法及びプログラム

Info

Publication number: JP6589227B1
Application number: JP2018192773A
Authority: JP
Inventors: 純次森川; 成治近藤; 稲葉　隆; 隆稲葉; 忠三蜷川
Original assignee: Mitsubishi Heavy Industries Thermal Systems Ltd
Current assignee: Mitsubishi Heavy Industries Thermal Systems Ltd
Priority date: 2018-10-11
Filing date: 2018-10-11
Publication date: 2019-10-16
Anticipated expiration: 2038-10-11
Also published as: JP2020060340A

Abstract

【課題】空調機の運転を最適化する制御情報を実用時間内に算出することができる制御装置を提供する。【解決手段】制御装置は、空調機の制御情報を算出する制御情報算出部と、前記制御情報に基づいて前記空調機を運転した場合の評価値を算出する評価部と、前記評価値が、所定の基準を満たすかどうかを判定する判定部と、を備え、前記判定部は、前記評価値の分布が正規分布に従う場合に、判定対象の前記評価値が、前記正規分布における所定の発生確率に相当する範囲に含まれるか否かによって前記評価値の判定を行う。【選択図】図２

Description

本発明は、制御装置、空調制御システム、制御情報の算出方法及びプログラムに関する。

電力の需要と供給のバランスを、例えば電力料金を上下させることによって調整し、ピーク時の電力供給量を抑制するＡＤＲ（Automated Demand Response）と呼ばれる技術が存在する。例えば、米国では、ＯｐｅｎＡＤＲによる電力料金のリアルタイム化が一般市場において展開されており、リアルタイムで安い電力供給会社を選択する等の運用が可能になっている。我が国においても政府主導のもとＡＤＲ規格として、米国のＯｐｅｎＡＤＲ２．０の国内版を展開する動きがある。

このような動きに対し、電力料金がリアルタイムに変化する条件下で、空調の快適さを損なわずに電力料金を抑制する空調システムの制御方法（ＲＴＰ（Real Time Pricing）制御）が提案されている（例えば、特許文献１）。

特開２０１８−０４０５１０号公報

ＲＴＰ制御では、空調機が使用する電力の上限値の条件（電力制限指令）を様々に変更しながらそれぞれの条件下における空調機の動作を予測し、当該動作に必要な電力料金および当該動作によって達成される空調の快適性に関する評価値を計算し、評価値が最適となるような電力制限指令を決定する処理を行う。対象とする空調機の台数等によっては、計算量が増加し、ＲＴＰ制御を実環境で利用しようとしても、実用時間内に最適な評価値となる電力制限指令を決定することができない可能性がある。

そこでこの発明は、上述の課題を解決することのできる制御装置、空調制御システム、制御情報の算出方法及びプログラムを提供することを目的としている。

本発明の一態様によれば、制御装置は、空調機の制御情報を算出する制御情報算出部と、前記制御情報に基づいて前記空調機を運転した場合の評価値を算出する評価部と、前記評価値が、所定の基準を満たすかどうかを判定する判定部と、サンプリング後の前記評価値の分布が正規分布と近似するように設定された所定の時間をかけてランダムに生成された前記制御情報に基づいて予めサンプリングされた複数の前記評価値と、を備え、前記判定部は、前記評価部が算出した判定対象の前記評価値が、前記予めサンプリングされた複数の前記評価値の分布を近似して得られる前記正規分布に基づいて設定される所定の範囲の値か否かによって、前記評価値の判定を行う。

本発明の一態様によれば、前記判定部は、判定対象の前記評価値の前記正規分布における下側確率が、所定の値以下となると当該評価値が前記範囲の値であると判定する。

本発明の一態様によれば、前記制御情報算出部は、ＳＡ（Simulated Annealing）アルゴリズムによって、前記評価値が前記範囲の値となるような前記制御情報を算出する。

本発明の一態様によれば、前記判定部は、前記評価値の算出を開始してからの経過時間が所定時間を上回るか、または、前記評価値の算出回数が所定回数を上回ると、前記評価値が所定の基準を満たすと判定する。

本発明の一態様によれば、前記制御装置は、前記制御情報を変更したときの前記評価値を複数サンプリングして、前記評価値の標本平均と標本標準偏差を算出し、前記基準の設定を受け付ける基準設定部、をさらに備える。

本発明の一態様によれば、前記基準設定部は、サンプリングした前記評価値が正規分布に従うように計算時間を定めて前記評価値をサンプリングする。

本発明の一態様によれば、前記制御装置は、時間毎の電力料金単価の情報を取得する料金情報取得部と、室外機と１台以上の室内機を備える空調機について、１台以上の前記空調機の運転状態に関する運転状態情報を取得する運転情報取得部と、前記運転状態情報と所定の予測モデルとに基づいて、未来の所定期間における前記空調機の運転状態を予測する予測部と、をさらに備え、前記評価部は、前記予測部が予測した運転状態に基づく前記所定期間における前記室内機が空調対象とするエリア毎の快適性の評価値と、前記時間毎の電力料金単価の情報に基づく前記所定期間における前記空調機の動作に必要な電力料金の評価値と、に基づく前記空調機の全てを対象とする運転に関する評価値を算出し、前記制御情報算出部は、前記評価値が所定の基準を満たすような前記制御情報を算出する。

本発明の一態様によれば、空調制御システムは、前記制御装置と、前記制御装置が算出する制御情報に基づいて運転する１台以上の空調機とを備える。

本発明の一態様によれば、制御情報の算出方法は、空調機の制御情報を算出するステップと、前記制御情報に基づいて前記空調機を運転した場合の評価値を算出するステップと、サンプリング後の前記評価値の分布が正規分布と近似するように設定された所定の時間をかけて前記制御情報をランダムに生成し、該制御情報に基づいて複数の前記評価値をサンプリングするステップと、前記評価値を算出するステップにおいて算出された判定対象の前記評価値が、前記サンプリングされた複数の前記評価値の分布を近似して得られる正規分布に基づいて設定される所定の範囲の値か否かによって、前記評価値が所定の基準を満たすかどうかを判定するステップと、を有する。

本発明の一態様によれば、プログラムは、制御装置のコンピュータを、空調機の制御情報を算出する手段、前記制御情報に基づいて前記空調機を運転した場合の評価値を算出する手段、サンプリング後の前記評価値の分布が正規分布と近似するように設定された時間をかけてランダムに生成された前記制御情報に基づいて複数の前記評価値をサンプリングする手段、前記評価値を算出する手段が算出した判定対象の前記評価値が、前記サンプリングされた複数の前記評価値の分布を近似して得られる正規分布に基づいて設定される所定の範囲の値か否かによって、前記評価値が所定の基準を満たすかどうかを判定する手段、として機能させる。

本発明によれば、空調機の運転を最適化する制御情報を短時間で算出することができる。

本発明の一実施形態における空調制御システムの一例を示す図である。本発明の一実施形態における制御装置の一例を示す機能ブロック図である。本発明の一実施形態における電力料金スケジュール情報の一例を示す図である。本発明の一実施形態における運転制限スケジュール情報の一例を示す図である。本発明の一実施形態における電力制限運転を説明する第一の図である。本発明の一実施形態における電力制限運転を説明する第二の図である。本発明の一実施形態における探索打ち切り処理を説明する第一の図である。本発明の一実施形態における探索打ち切り処理を説明する第二の図である。本発明の一実施形態における探索打ち切り処理の効果を説明する図である。本発明の一実施形態における制御装置の処理の一例を示す第一のフローチャートである。本発明の一実施形態における制御装置の処理の一例を示す第二のフローチャートである。本発明の一実施形態における制御装置のハードウェア構成の一例を示す図である。

＜実施形態＞
以下、本発明の一実施形態による空調制御システムおよび最適な制御情報の算出方法について図１〜図１２を参照して説明する。
図１は、本発明の第一実施形態における空調制御システムの一例を示す図である。空調制御システム１００は、ビル４に備えられた１台以上空調機３１〜３ｍを、ユーザの快適性を保ちつつ低コストで運転するための制御システムである。さらに本実施形態の空調制御システム１００は、エリア毎にそのエリアに対する空調の優先度を設定したうえで、空調機３１〜３ｍ全体での経済性やエリア毎の快適性を所望の状態に制御することを特徴とする。ここでエリアとは、空調対象の単位となる空間のことである。エリアは、ビル４内の１つの部屋、１つのフロアであってもよいし、複数のビルを跨る複数の部屋や複数のフロアを１つのエリアとして設定してもよい。また、１つのエリアには、異なる冷媒系統（ブロック）に属する複数の室内機が設置されていてもよい。

空調制御システム１００は、ＤＲＡＳサーバ１と、ＤＲＡＳクライアント２と、ゲートウェイ３と、空調機３１〜３ｍと、を含んで構成される。ＤＲＡＳサーバ１は、例えば電力会社が運用するサーバ端末装置である。電力会社は、電力料金単価を例えば１０分毎に変更する。例えば、電力会社は、電力の需要による電力不足が見込まれる場合には、電力料金単価を上昇させて需要家に節電を促し、電力需要が比較的少ないと見込まれる時間には電力料金単価を下げて電力の消費を促す。ＤＲＡＳサーバ１とＤＲＡＳクライアント２とは通信可能に接続されており、一般に公開されたプロトコル（例えばＯｐｅｎＡＤＲ２．０等）に基づく通信を行う。例えば、ＤＲＡＳサーバ１は、未来における所定期間分（例えば、３０分間）の電力料金のスケジュール情報を、ＤＲＡＳクライアント２に送信する。

ＤＲＡＳクライアント２は、例えば空調機３１等の運用を行う企業が運用するサーバ端末装置である。ＤＲＡＳクライアント２は、ＤＲＡＳサーバ１から取得した電力料金のスケジュール情報に基づいて、どのような運転を行えば、ユーザの設定温度を実現しつつ安い電力料金で空調機３１〜３ｍを運転できるかを示す運転制限スケジュール情報を算出する。ＤＲＡＳクライアント２は、インターネット等のネットワークを介して、ゲートウェイ３と通信可能に接続されている。ＤＲＡＳクライアント２は、算出した運転制限スケジュール情報をゲートウェイ３に送信する。

ゲートウェイ３は、ビル４に設けられた通信装置である。ビル４には、複数の室内機を有するマルチ型の空調機３１〜３ｍが設けられている。空調機３１は、室外機３０１と、室内機３１１〜３１ｎを備える。空調機３２は、室外機３０２と、室内機３２１〜３２ｎを備える。空調機３３，３４，３ｍについても同様である。例えば、室外機３０１は、圧縮機、熱交換器等を有しており、冷媒を室内機３１１，３１ｎ等へ送出する。室内機３１１等は、ビル４の部屋に設けられており、その部屋の温度を、ユーザ所望の温度となるように空調を行う。空調機３１〜３ｍの各々は、１つの冷媒系統である。つまり、室外機３０１と、室内機３１１〜３１ｎ等は冷媒配管で接続された１つの制御単位である。本実施形態では、この制御単位をブロックと呼ぶ。例えば、空調機３１をブロック１、空調機３２をブロック２、空調機３ｍをブロックｍと呼ぶ。空調制御システム１００では、空調機３１等の制御はブロック単位で行うが、温度制御はエリア単位で行う。例えば、エリアａ１は１つの部屋であって、エリアａ１の温度は複数の室内機３１１、３１２、３２１、３２２を用いて制御する。室内機３１１、３１２はブロック１、室内機３２１、３２２はブロック２に属している。同様にエリアａ２は１つのフロア全体であって、エリアａ２の温度を制御する室内機３１ｎ〜３ｍｎは、それぞれブロック１〜ｍに属している。エリアが複数のブロックに跨ることは必須ではない。例えば、エリアａ３は、室内機３４２〜３４３によって制御されるが、室内機３４２〜３４３は何れもブロック４に属している。また、エリアは、エリアａ４のように１つの室内機３ｍ１で制御される空間であってもよい。本実施形態では、エリア単位で室温達成の優先度を設定することができる。例えば、エリアａ１の温度を重視する場合、エリアａ１の優先度を高く設定することで、他のエリアの温度を多少犠牲にしても、エリアａ１の温度をより高い精度で所望の温度に制御することができる。また、室温の達成をあまり重視しない部屋については、各部屋をそれぞれ１つずつのエリアとして設定してもよいし、それらの部屋をまとめて一つのエリアとして設定してもよい。複数ブロック全体での電力料金を抑えつつ、優先度の高いエリアの快適性を重視して空調快適性を実現する制御を、エリア複合ＲＴＰ（Real Time Pricing）制御と呼ぶ。

ゲートウェイ３は、ＤＲＡＳクライアント２から取得した空調機毎の運転制限スケジュール情報を、対応する空調機３１〜３ｍへ出力する。空調機３１〜３ｍの各々は取得した運転制限スケジュール情報に基づいて、自装置の運転を制御する。

図２は、本発明の一実施形態における制御装置の一例を示す機能ブロック図である。
図２を用いてＤＲＡＳクライアント２の構成について説明する。ＤＲＡＳクライアント２は、運転情報取得部２１と、料金情報取得部２２と、制御情報算出部２３と、判定部２４と、運転状態予測部２５と、評価部２６と、基準設定部２７と、通信部２８と、記憶部２９とを備える。
運転情報取得部２１は、空調機３１〜３ｍの運転情報を取得する。運転情報とは、例えば、空調機３１の運転による消費電力量、運転時の室内温度、設定温度および室外温度等である。運転情報取得部２１は、例えば所定の時間毎に空調機３１〜３ｍの運転情報を取得し、記憶部２９に取得した運転情報を記録する。
料金情報取得部２２は、ＤＲＡＳサーバ１が送信した電力料金のスケジュール情報を取得する。電力料金のスケジュール情報については、後に図３を用いて説明する。

制御情報算出部２３は、所定の最適化手法を用いて、ユーザの快適性を保ちつつ、電力料金を抑制する空調機３１〜３ｍの運転を実現する運転制限スケジュール情報を算出する。所定の最適化手法とは、例えば、ＳＡ（Simulated Annealing）アルゴリズム（焼きなまし法）、遺伝的アルゴリズム、粒子群最適化などの手法である。また、ユーザの快適性を保ちつつ、電力料金を抑制する運転制限スケジュール情報とは、後述する評価部２６が算出する評価値Ｊｃの値が最も小さな値となる場合の運転制限スケジュールのことをいう。本実施形態の場合、運転制限スケジュール情報とは、所定時間毎の空調機３１〜３ｍのそれぞれが消費する電力の上限値である。

判定部２４は、評価部２６が算出する評価値Ｊｃが、所定の基準を満たすかどうかを判定する。判定部２４は、評価値の分布が正規分布に従う場合に評価対象となる評価値が、正規分布に基づいて設定される所定の範囲の値か否かによって、基準を満たすかどうかの判定を行う。より具体的には判定部２４は、評価対象となる評価値の発生確率が、正規分布における所定の下側確率以下となると基準を満たすと判定する。

運転状態予測部２５は、運転制限スケジュール情報が示す条件下での空調機３１，３２，３ｍ等の運転状態を予測する。例えば、運転状態予測部２５は、運転情報取得部２１が取得した空調機３１の現在の運転状態に基づいて数分先の空調機３１の運転状態を予測する。また、例えば、運転状態予測部２５は、現在の運転状態に基づいて予測した数分先の空調機３１の運転状態に基づいて、さらに数分先の空調機３１の運転状態を予測する。運転状態の予測には、例えば、ニューラルネットワークと呼ばれる機械学習アルゴリズムによる予測モデルを用いる。

評価部２６は、運転状態予測部２５が予測した所定時間先の空調機３１〜３ｍの運転状態を後述する評価関数によって評価する。より具体的には、評価部２６は、電力料金評価部２６１と快適性評価部２６２とを備えており、電力料金（コスト）とユーザにとっての快適性の両面から運転状態予測部２５が予測した未来における空調機３１〜３ｍの運転に対する評価を行う。

基準設定部２７は、判定部２４が判定に用いる基準を設定する。基準設定部２７は、例えば、ユーザが入力した評価値Ｊｃの正規分布における下側確率を基準として設定する。

通信部２８は、他装置との通信を行う。例えば、通信部２８は、ＤＲＡＳサーバ１とデータ通信を行い、電力料金のスケジュール情報を取得する。あるいは、通信部２８は、ゲートウェイ３とデータ通信を行い、運転制限スケジュール情報をゲートウェイ３へ送信する。

記憶部２９は、運転制限スケジュール情報の算出に必要な種々のデータを記憶する。例えば、記憶部２９は、運転状態予測部２５が運転状態の予測に用いる予測モデルを記憶している。この予測モデルは、ある運転制限スケジュール情報を与えた場合の過去の運転状態の実績情報に基づいてニューラルネットワーク等の手法で生成した予測モデルである。この予測モデルは、電力制限スケジュール情報と、実際の消費電力量、室内温度、室外温度を入力すると、所定時間先の未来における室内温度とそのときまでに消費される消費電力量の予測値を出力する。本実施形態の場合、予測モデルは、１分毎の室内温度と消費電力量を５分先まで予測する。予測モデルは、例えば、空調機３１〜３ｍ毎に用意されている。

図３は、本発明の一実施形態における電力料金スケジュール情報の一例を示す図である。電力料金のスケジュール情報には、所定時間（例えば１０分）毎の電力料金単価が含まれている。図３の例の場合、本日の１４：００〜１４：１０では１ｋＷｈあたりの電力料金が３０円、１４：１０〜１４：２０では１０円、１４：２０〜１４：３０では１００円となっている。ＤＲＡＳサーバ１は、例えば、時間が１０分経過する毎に、３０分先までの１０分毎の電力料金単価が含まれたスケジュール情報をＤＲＡＳクライアント２へ送信する。ＤＲＡＳクライアント２では、料金情報取得部２２が電力料金単価のスケジュール情報を、通信部２８を介して取得し、記憶部２９へ記録する。

図４は、本発明の一実施形態における運転制限スケジュール情報の一例を示す図である。運転制限スケジュール情報には、所定時間（例えば５分）毎の電力制限値が含まれている。
この電力制限値は、図３で例示した電力料金単価の場合に、ユーザが設定した室内の設定温度をなるべく達成しつつ、最安の電力料金で空調機３１等を運転するための空調機３１等が消費する電力の上限値である。
図４の例の場合、単価が３０円の時間帯（１４：００〜１４：１０）では、前半が１４ｋＷ、後半が７ｋｗとなっている。続いて、単価が１０円の時間帯（１４：１０〜１４：２０）では前半が３ｋＷ、後半が７ｋｗとなっている。最後に単価が１００円の時間帯（１４：２０〜１４：３０）では前半、後半ともに０ｋｗとなっている。
ＤＲＡＳクライアント２では、制御情報算出部２３が算出した運転制限スケジュールが示す条件下での３０分後までの空調機３１〜３ｍの運転状態を運転状態予測部２５が予測し、評価部２６がその３０分間の運転による電力料金、エリア毎の快適性を評価する。そして、制御情報算出部２３は、評価部２６による評価結果に基づいて、より適切な運転制限スケジュール情報を再設定する。例えば、快適性が良好でも電力料金が高額になる場合には、制御情報算出部２３は運転制限スケジュールの電力制限値を抑え、より低電力で運転を行う条件を再設定する。また、例えば、電力料金は非常に安くなるが快適性も低い場合、制御情報算出部２３は運転制限スケジュールの電力制限値を上昇させ、電力料金が相対的に上昇してもユーザが設定した設定温度を実現することにより、快適性を向上させる条件を再設定する。図４に例示する電力制限運転のスケジュール情報は、このようにして最適化された所定時間毎の電力制限値である。

次に評価部２６が使用する評価関数について説明する。以下の式（１）〜（４）が本実施形態の評価関数の一例である。

ここで、αは電力料金評価に対する重み係数、ｂはブロックの識別番号（ｂ＝１,２,・・・，ｂ_ｍａｘ）、ｍは単位時間（例えば５分間）の番号、Ｒ（ｍ）はｍ単位時間時の電力単価、Ｗ_５ ^ｂ（ｍ）はｍ単位時間時におけるブロックｂでの空調機３ｂによる使用電力量［Ｗｈ］、ａはエリアの識別番号（ａ＝１,２,・・・，ａ_ｍａｘ）、Ｚ^ａ（ｍ）はｍ単位時間時におけるエリアａの快適性評価値である。β^ａはエリアａに対する空調の優先度（−１．０≦β^ａ≦１．０）であって式（３）が成立する。Ｔ^ａ _ＳＡＡ（ｍ）はｍ単位時間時におけるエリアａの平均室温偏差（℃）である。ｉは室内機の識別番号（ｉ＝１，２・・・Ｉｍａｘ）、Ｔ^ｂ _Ａｉ（ｍ）はｍ単位時間時における対象となるエリアａに属する室内機ｉのセンサが計測した室内温度（吸込み温度）（℃）、Ｔ^ｂ _Ｓｉ（ｍ）はｍ単位時間時における対象となるエリアａについて設定された室内機ｉのユーザ設定温度（℃）、Ｃ^b _Ｐｉはブロックｂの室内機ｉの定格冷暖房能力［ｋＷ］、Ｓ^b _ＲＳｉはブロックｂの室内機ｉの運転状態（０が停止、１が運転中）である。

電力料金評価部２６１は、式（１）右辺の第１項の式により電力料金の評価値を算出する。快適性評価部２６２は、式（１）右辺の第２項の式と式（２）〜式（４）により快適性の評価値を算出する。式（１）右辺の第１項はブロック別の単位時間あたりの電力料金を所定期間にわたって合計し、合計した値を全ブロックについて集計する式である。電力料金が高い程、第１項の値は大きくなる。第２項は、単位時間あたりの室内機毎のユーザによる設定温度と実際の室内温度（空調温度）の温度差に基づく値に対して、エリア別の優先度係数β^ａを乗じて２乗した値を所定期間にわたって合計し、合計した値を全エリアについて集計する式である。第２項は、ユーザが設定した設定温度と実際の室内温度の温度差が大きい程（快適性が低い程）、大きな値となる。また、優先度係数β^ａが大きい、つまり優先されるエリアａでの温度差が大きいほど、第２項の値は大きな値となる。従って、電力料金が安く、各エリア（特に優先度の高いエリア）での快適性が高い程、式（１）の値は小さな値となる。制御情報算出部２３は、評価部２６が算出した式（１）による評価値Ｊｃがなるべく小さくなるように電力制限運転のスケジュール情報を調整する。

図５は本発明の一実施形態における電力制限運転を説明する第一の図である。
図６は本発明の一実施形態における電力制限運転を説明する第二の図である。
図５、図６において縦軸は消費電力量（ｋＷ）、横軸は時間（分）を示している。図６の場合、縦軸はさらに電力制限値（ｋＷ）を示している。
図５に示すグラフ５１は、平常運転時の制御（本実施形態の電力制限値による運転制限を設けない制御）によって空調機３１を運転したときの消費電力量の時間推移を示している。
図６に示すグラフ６１は、本実施形態のエリア複合ＲＴＰ制御によって空調機３１〜３ｍを運転したときの空調機３１についての消費電力量の時間推移を示している。また、グラフ６２は、制御情報算出部２３が算出した空調機３１に対する運転制限スケジュール（電力制限値の時間推移）を示している。

まず、運転状態予測部２５が、運転情報取得部２１が取得した最新の空調機３１の運転状態情報（電力制限値、実際の電力、室内温度、設定温度、室外温度）と予測モデルによって５分先までの１分毎の空調機３１の運転状態を予測する。次に運転状態予測部２５が、自らが予測した５分後の空調機３１の運転状態情報と予測モデルによってさらに５分先（つまり時間５分から１０分）までの１分毎の空調機３１の運転状態を予測する。運転状態予測部２５は、同様の５分単位の予測を計６回繰り返し、時間０を基準として３０分後までの運転状態を模擬する。なお、これら６回の模擬において電力制限値については任意の値を与える。
次に評価部２６が式（１）に、運転状態予測部２５が模擬した３０分間の運転状態情報を代入し評価値を算出する。なお、α、β^ａ、Ｃ^b _Ｐｉ等の定数については予め設定され記憶部２９に記録されている。また時間毎の電力料金単価Ｒ（ｍ）には料金情報取得部２２が取得した電力料金スケジュール情報を用いる。各エリアの目標温度Ｔ^ｂ _Ｓｉ（ｍ）については、各室内機ｉを通じて設定された値を用いる。

次に制御情報算出部２３が、評価部２６が算出した評価値に基づいて、５分毎の電力制限値の最適化を行う。例えば、開始から５分間（１回目の模擬区間）は目標の設定温度を達成するために高めの電力制限値を設定し、空調機３１の運転を促進する。また、４回目の模擬区間（１５〜２０分）については、次に電力料金単価が１００円と高額になる区間なので、高額区間（５〜６回目の模擬区間）での運転を抑制するために予め強めの運転を行って部屋を涼しくしておくなど、電力料金単価の変動と設定温度と室内温度との乖離具合（快適性）とに応じた電力制限値を設定する。

また、制御情報算出部２３は、目標の設定温度を達成すべく、各エリアａに付された優先度係数β^ａを参照して、優先度係数β^ａが大きな室内機ｉの室温が設定温度を達成できるような制御情報を算出する。例えば、制御情報算出部２３は、優先度係数β^ａに大きな値が設定された室内機を多く含むブロックの空調機３１等への電力制限値は、確実に目標温度が達成できるように十分に大きな値を設定し、優先度係数β^ａに小さな値が設定された室内機を多く含む空調機の電力制限値には、相対的に小さな値を設定する。

次に運転状態予測部２５は、制御情報算出部２３が設定した６つの区間（３０分間）に対する新たな運転制限スケジュール情報（時間毎の電力制限値、または時間毎の電力制限値および圧縮機の回転数）を用いて、再度、未来の３０分間における空調機３１の運転状態を模擬する。また、評価部２６はその模擬結果を式（１）の評価関数によって評価する。また、制御情報算出部２３が、ＳＡアルゴリズムなどの最適化手法を用いて評価値Ｊｃを最小にする運転制限スケジュール情報を算出する。

このような運転制限スケジュール情報の最適化を繰り返し、評価値Ｊｃの値を最小にする空調機３１〜３ｍごとの運転制限スケジュールが決定されると、制御情報算出部２３が最適化された運転制限スケジュール情報を、通信部２８を介してゲートウェイ３に送信する。ゲートウェイ３は最適化された運転制限スケジュール情報を空調機３１〜３ｍに送信し、空調機３１〜３ｍは、運転制限スケジュール情報が示す５分毎の電力制限値を超えないように運転を行う。これにより、低コスト且つ高い快適性を提供する空調制御が可能になる。

ここで、制御情報算出部２３がｍ＝Ｎｈまで先の電力制限スケジュール情報を算出すると仮定すると、ブロックｂに対する電力制限指令の順列は（Ｐ^ｂ _L1，Ｐ^ｂ _L2，・・・，Ｐ^ｂ _LNh）となり、ビル全館に配置された全ての空調機を探索対象とする電力指令値の順列は以下のようになる。

すると探索アルゴリズムは、Ｎ_ＡＬＬ＝Ｎ_Ｌ ^{Ｎｂ×Ｎｈ}の状態を有する空間を、最適解を求めて探索することになる。Ｎ_Ｌは電力制限値の段数であり、例えば、電力制限値を５，１０，１５，２０，２５（ＫＷ）の５段階の中から選択する場合はＮ_Ｌ＝５である。Ｎｂはブロック数（室外機の数＝ｂ_ｍａｘ）である。Ｎｈは１度に何分先まで予測するかであり、例えば、１単位時間が５分で１０分先まで予測するならばＮｈ＝２である。例えば、Ｎ_Ｌ＝５、Ｎｂ＝１０、Ｎｈ＝２という典型的な例の場合、探索空間は、Ｎ_ＡＬＬ＝５²⁰となる。この大きさの探索空間を全探索すると、数百台以上のスーパーコンピュータを用いても何時間もかかってしまい実用的ではない。データ取得から電力制限指令値を空調機に与えて、空調機がその電力制限指令値に従って動作を始めるまでに要する時間を考慮すると、実際の空調機の制御においては、最適解（最小の評価値）の算出に割り当てることができる処理時間は１０秒程度である。そこで、本実施形態では、ＳＡアルゴリズムなどの最適化手法による探索に加え、実用時間内に十分な精度を有する解を得られた時点で探索打ち切り処理を行う。

図７は、本発明の一実施形態における探索打ち切り処理を説明する第一の図である。
図７にマルチコアのＣＰＵによって、式（１）の評価値Ｊｃが最小となる最適解を探索したときの様子を示す。縦軸は評価値Ｊｃの値、横軸は探索開始からの経過時間を示している。各コアにＳＡアルゴリズムを実行させて並列に解の探索を実行させたことにより、複数のグラフが表示されている。０．５秒付近で一旦１つの値に収束しているのは、それまでに得られた中で最も良質な解を全てのコアで実行中のアルゴリズムに与え同期させる処理を行ったことによる。図７に示すように１秒の手前から以降は評価値Ｊｃの値に大きな変化は無く、その後わずかに低下して１０秒に至っている。探索開始から１秒を経過する手前の評価値の値は１５００を下回り、その後、１０秒経過時には評価値Ｊｃの値は１３９７を記録した。もし、評価値Ｊｃの値が、例えば１４５０程度であっても実用上問題が無い程度に低コスト且つ高快適性を実現することができれば、上限となる１０秒間をかけて探索処理を行う必要は無い。図７の例では、解の探索は１秒弱で十分と考えられる。そこで、本実施形態では、所定の基準を設け、探索中の解がこの基準に達したら探索を打ち切り、基準を満たしたと判定した時点の準最適解を最終的な解として採用する。

次に本実施形態の探索打ち切り処理について説明する。
図８は、本発明の一実施形態における探索打ち切り処理を説明する第二の図である。
図８のドットで示す領域Ｋ１は、一例として、Ｎ_Ｌ＝５、Ｎ_ｂ＝２０、Ｎ_ｈ＝６（Ｎ_ＡＬＬ＝５¹²⁰）の探索空間から１００００点をサンプリングしたときの評価値Ｊｃのヒストグラムである。横軸に評価値Ｊｃの値、縦軸に出現率を示す。サンプリングした値の標本平均Ｘ_ＪＣは２０７２、標本標準偏差Ｓ_ＪＣは１２８．３である。図８のグラフＫ２は、標本平均Ｘ_ＪＣと標本標準偏差Ｓ_ＪＣを用いて算出した正規分布のグラフである。図示するようにサンプリングした評価値Ｊｃの分布は正規分布に近似できる。他のサンプリング例でも評価値Ｊｃの分布が概ね正規分布となることが実験的に確認できている。そこで、本実施形態では、探索空間を構成する評価値Ｊｃが正規分布に従うと仮定し、以下の処理手順により解が正規分布の下位何％に到達しているのかを計算し、所定の基準より下位に達していれば（小さな値となっていれば）探索を打ち切る処理を行う。

［１］．電力制限指令値の順列をランダムに複数生成して、各順列で空調機を運転したときの評価値Ｊｃの評価値をサンプリングする。［２］．サンプリングした評価値の標本平均Ｘ_ＪＣと標準偏差Ｓ_ＪＣを算出する。［３］．下側確率Ｐ_ＪＣ ^ＥＮＤを設定する。［４］．解の探索を開始する。［５］．探索中の評価値Ｊ_Ｃの標準化変数Ｚ_ＪＣ＝（評価値Ｊ_Ｃ−Ｘ_ＪＣ）／Ｓ_ＪＣを算出する。［６］．累積分布関数（式（６））により、Ｚ_ＪＣの下側確率Ｐ_ＪＣ（％）を求める。

上記の式（７）は確率密度関数である。［７］．下側確率Ｐ_ＪＣ（％）が上で設定した下側確率Ｐ_ＪＣ ^ＥＮＤに到達したら（下側確率Ｐ_ＪＣ≦下側確率Ｐ_ＪＣ ^ＥＮＤとなったら）探索を打ち切る。例えば、図８に示すＰｘが、下側確率がＰ_ＪＣ ^ＥＮＤとなる点であるとすると、評価値Ｊ_ＣがＰｘ以下の値となっていれば、探索を打ち切る。

図９は、本発明の一実施形態における探索打ち切り処理の効果を説明する図である。
図９に実際に上記手順で探索を打ち切った実行例を示す。図中Ｊ´_ｃｍｉｎは、予め全探索により発見した厳密最小解Ｊ´_ｃｍｉｎ＝１１１である。一方、上記の［１］〜［７］の手順で解の探索及び探索の打ち切りを実行したところ、探索の開始から０．０２秒が経過した時に厳密最小解Ｊ´_ｃｍｉｎと同じ値の解Ｊ_ｃｍｉｎ＝１１１を得ることができた。

図１０は、本発明の一実施形態における制御装置の処理の一例を示す第１のフローチャートである。
図１０を用いて、解の探索および探索打ち切り処理に先立ち必要となる処理の流れを説明する。まず、基準設定部２７が、これから行う空調と同様の条件（例えば、各エリアの設定温度や何分先までの評価値を算出するか等）を設定して、ランダムに生成した複数の電力制限スケジュール情報に基づいて評価値Ｊｃのサンプリングを行う（ステップＳ１）。運転状態の予測には、運転状態予測部２５を用いることができる。また、１つの評価値Ｊｃのサンプルを得るための計算時間を定め、その計算時間で得られた評価値Ｊｃを複数（例えば、１万個）算出する。サンプルを得るための計算時間には、サンプリングした複数の評価値Ｊｃの分布が正規分布に従うような適切な時間を設定することができる。例えば、図７に例示した例では計算時間を１０秒として設定すると、多数の解が最適解付近に収束してしまい、サンプリングした評価値Ｊｃが正規分布とならない。従って、図７に示す例では、計算時間を、例えば０．８〜１秒程度に制限してサンプリングを行う。次に基準設定部２７は、サンプリングした複数の評価値Ｊｃに基づいて、標本平均Ｘ_ＪＣと標本標準偏差Ｓ_ＪＣを算出する（ステップＳ２）。次に基準設定部２７は、基準となる所定の下側確率を設定する（ステップＳ３）。下側確率Ｐ_ＪＣ ^ＥＮＤ（例えば、０．０１％）は、ユーザが入力してもよい。次に基準設定部２７は、標本平均Ｘ_ＪＣと標本標準偏差Ｓ_ＪＣと下側確率Ｐ_ＪＣ ^ＥＮＤを記憶部２９に記録する（ステップＳ４）。

図１１は、本発明の一実施形態における制御装置の処理の一例を示す第２のフローチャートである。
図１１を用いて運転制限スケジュール情報の生成処理の流れについて説明する。記憶部２９には、式（１）〜（４）に用いられる定数（α，β^ａ等）が記録されている。また、予め図１０で説明した処理により、記憶部２９には標本平均Ｘ_ＪＣと標本標準偏差Ｓ_ＪＣと下側確率Ｐ_ＪＣ ^ＥＮＤが記録されている。空調制御システム１００の構成は、図１の構成を前提とする。

まず、料金情報取得部２２が、ＤＲＡＳサーバ１が生成した電力料金単価スケジュール情報を取得し（ステップＳ１１）、記憶部２９に記録する。電力料金スケジュール情報には、例えば、３０分先までの１０分毎の電力料金単価が含まれている。
次に運転情報取得部２１が、複数の空調機３１〜３ｍそれぞれの運転状態情報を取得する（ステップＳ１２）。運転状態情報には、ブロック単位での単位時間毎の消費電力量、外気温度、室内機単位での室内温度、ユーザによる設定温度等の情報が含まれている。運転情報取得部２１は、取得した運転状態情報を記憶部２９に記録する。

次に制御情報算出部２３が、空調機３１〜３ｍのそれぞれについて、未来の所定期間（例えば３０分間）における所定時間毎（例えば５分毎）の運転制限スケジュール情報を算出する（ステップＳ１３）。例えば、制御情報算出部２３は、初回では任意の初期値を３０分間にわたる５分毎の電力制限値として設定してもよい。

次に運転状態予測部２５が、空調機３１〜３ｍのそれぞれについて、未来の所定期間における所定時間毎の運転状態を予測する（ステップＳ１４）。より具体的には、運転状態予測部２５は、ステップＳ１２で取得した運転状態情報と、ステップＳ１３で設定した運転制限スケジュール情報とを予測モデルに入力し、未来の５分間における１分毎の空調機単位の消費電力量と、室内機単位の室内温度を予測する。次に運転状態予測部２５は、ステップＳ１２で取得した運転状態情報とステップＳ１３で設定した運転制限スケジュール情報とのうち、ステップＳ１２で取得した実際の消費電力量と室内温度とを、予測モデルによって予測した５分後の予測値で置き換えた情報を、予測モデルに入力し、さらに５分先（開始時間から５分後〜１０分後までの５分間）までの１分毎の空調機単位の消費電力量と、室内機単位の室内温度を予測する。以下、同様にして５分毎の予測を計６回繰り返し、３０分間分の消費電力量と、室内温度の予測を行う。

次に評価部２６が、運転状態予測部２５によって予測された運転状態に対する評価値を算出する（ステップＳ１５）。具体的には、ステップＳ１４の予測により、１分毎の消費電力量、室内機毎の室内温度の予測値が算出されているので、評価部２６は、式（１）〜式（４）を用いてビル４全体の評価値Ｊｃを算出する。評価部２６は、評価値Ｊｃを判定部２４へ出力する。

次に判定部２４が、評価値Ｊｃが基準に達したか否かを判定する（ステップＳ１６）。具体的には、判定部２４が、評価対象の評価値Ｊｃと、記憶部２９が記憶する標本平均Ｘ_ＪＣおよび標本標準偏差Ｓ_ＪＣに基づいて、標準化変数Ｚ_ＪＣを算出する。続いて判定部２４は、上記の式（６）、（７）を用いて、標準化変数Ｚ_ＪＣの下側確率Ｐ_ＪＣを算出する。次に判定部２４は、下側確率Ｐ_ＪＣと記憶部２９が記憶する下側確率Ｐ_ＪＣ ^ＥＮＤを比する。下側確率Ｐ_ＪＣ≦下側確率Ｐ_ＪＣ ^ＥＮＤが成立すると、判定部２４は評価値Ｊｃが基準に達したと判定する。下側確率Ｐ_ＪＣ＞下側確率Ｐ_ＪＣ ^ＥＮＤの場合、判定部２４は評価値Ｊｃが基準に達していないと判定する。

予め基準設定部２７によってサンプリングされた評価値Ｊｃによっては、しばらく計算を繰り返しても、下側確率Ｐ_ＪＣ≦下側確率Ｐ_ＪＣ ^ＥＮＤが成立しない可能性がある。そのような場合、判定部２４は、解の探索から所定の上限時間（例えば、１０秒）が経過すると評価値Ｊｃが基準に達したと判定する。又は、判定部２４は、評価値の算出回数が規定回数に到達すると、評価値Ｊｃが基準に達したと判定してもよい。

評価値が基準に到達した場合（ステップＳ１６；Ｙｅｓ）、制御情報算出部２３は、評価値Ｊｃが基準を満たす場合の空調機３１〜３ｍ毎の運転制限スケジュール情報を、ゲートウェイ３へ出力する（ステップＳ１７）。ゲートウェイ３は、運転制限スケジュール情報を対応する空調機３１〜３ｍへ出力し、各空調機３１〜３ｍは、運転制限スケジュール情報に従って、次の３０分間の運転を行う。これにより、ユーザ所望の設定温度を達成しつつ、空調機の電力料金を低減することができる。

評価値が基準に到達しない場合（ステップＳ１６；Ｎｏ）、評価値が基準に到達するまで、ステップＳ１３からの処理を繰り返す。具体的には、制御情報算出部２３がＳＡアルゴリズムなどの最適化手法を用いて、式（１）による評価値を最小化する新たな運転制限スケジュール情報を算出する（ステップＳ１３）。

本実施形態の空調制御システム１００によれば、電力料金単価が変動する環境下でも、優先度係数β^ａに応じたユーザの快適性を実現しつつ、建物全体での電力料金を抑えた複数の空調機の運転が可能になる。また、その際に最適な運転制限スケジュール情報を算出する必要があるが、本実施形態の探索打ち切り処理によれば、最適化手法によって最適な運転制限スケジュール情報を算出するにあたり、所定の基準により解の探索を継続するか否かを判定する。当該基準に基づいて解の探索を継続するか否かを判定することにより、厳密な最適解の発見に必要な時間よりも大幅に短縮した時間内で、十分実用に耐える準最適解（低コスト且つ高快適性を十分に実現できる運転制限スケジュール情報）を算出することができる。これにより、制御対象となる空調機（ブロック）が多数存在する場合や電力制限値の段数が多い場合などでも効果的に計算量を削減し、処理の高速化、計算負荷の低減を図り、遅延なく空調機に電力制限指令を行うことができる。

なお、上記例では、運転制限スケジュール情報の算出を行う機能を、ＤＲＡＳクライアント２に実装する場合を例に説明を行ったが、これに限定されない。例えば、これらの機能をゲートウェイ３に実装してもよい。
また、上記例では最適化手法として、ＳＡアルゴリズムを用いたが、探索中の解の分布が正規分布に従うならば、他の方法であっても本実施形態の探索打ち切り処理を適用することができる。

図１２は、本発明の一実施形態における制御装置のハードウェア構成の一例を示す図である。
コンピュータ９００は、ＣＰＵ９０１、主記憶装置９０２、補助記憶装置９０３、入出力インタフェース９０４、通信インタフェース９０５を備える。
上述のＤＲＡＳクライアント２（制御装置）は、コンピュータ９００に実装される。そして、上述した各機能は、プログラムの形式で補助記憶装置９０３に記憶されている。ＣＰＵ９０１は、プログラムを補助記憶装置９０３から読み出して主記憶装置９０２に展開し、当該プログラムに従って上記処理を実行する。また、ＣＰＵ９０１は、プログラムに従って、記憶領域を主記憶装置９０２に確保する。また、ＣＰＵ９０１は、プログラムに従って、処理中のデータを記憶する記憶領域を補助記憶装置９０３に確保する。

なお、ＤＲＡＳクライアント２の全部または一部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各機能部による処理を行ってもよい。ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、ＣＤ、ＤＶＤ、ＵＳＢ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また、このプログラムが通信回線によってコンピュータ９００に配信される場合、配信を受けたコンピュータ９００が当該プログラムを主記憶装置９０２に展開し、上記処理を実行しても良い。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。なお、ＤＲＡＳクライアント２は、複数のコンピュータ９００によって構成されていても良い。

その他、本発明の趣旨を逸脱しない範囲で、上記した実施の形態における構成要素を周知の構成要素に置き換えることは適宜可能である。また、この発明の技術範囲は上記の実施形態に限られるものではなく、本発明の趣旨を逸脱しない範囲において種々の変更を加えることが可能である。なお、ＤＲＡＳクライアント２は制御装置の一例である。

１００・・・空調制御システム
１・・・ＤＲＡＳサーバ
２・・・ＤＲＡＳクライアント
３・・・ゲートウェイ
３１、３２、３３、３４、３ｍ・・・空調機
４・・・ビル
３０１、３０２、３０３、３０４、３０ｍ・・・室外機
３１１、３１２、３１３、３１ｎ、３２１，３２２、３２３、３２ｎ、３３１、３３２、３３３、３３ｎ、３４１、３４２、３４３、３４ｎ、３ｍ１、３ｍ２、３ｍ３、３ｍｎ・・・室内機
２１・・・運転情報取得部
２２・・・料金情報取得部
２３・・・制御情報算出部
２４・・・判定部
２５・・・運転状態予測部
２６・・・評価部
２６１・・・電力料金評価部
２６２・・・快適性評価部
２７・・・基準設定部
２８・・・通信部
２９・・・記憶部

Claims

空調機の制御情報を算出する制御情報算出部と、
前記制御情報に基づいて前記空調機を運転した場合の評価値を算出する評価部と、
前記評価値が、所定の基準を満たすかどうかを判定する判定部と、
サンプリング後の前記評価値の分布が正規分布と近似するように設定された所定の時間をかけてランダムに生成された前記制御情報に基づいて予めサンプリングされた複数の前記評価値と、
を備え、
前記判定部は、前記評価部が算出した判定対象の前記評価値が、前記予めサンプリングされた複数の前記評価値の分布を近似して得られる正規分布に基づいて設定される所定の範囲の値か否かによって、前記評価値の判定を行う、
制御装置。
前記判定部は、判定対象の前記評価値の前記正規分布における下側確率が、所定の値以下となると当該評価値が前記範囲の値であると判定する、
請求項１に記載の制御装置。
前記制御情報算出部は、ＳＡ（Simulated Annealing）アルゴリズムによって、前記評価値が前記範囲の値となるような前記制御情報を算出する、
請求項１または請求項２に記載の制御装置。
前記判定部は、前記評価値の算出を開始してからの経過時間が所定時間を上回るか、または、前記評価値の算出回数が所定回数を上回ると、前記評価値が所定の基準を満たすと判定する、
請求項１から請求項３の何れか１項に記載の制御装置。
前記制御情報を変更したときの前記評価値を複数サンプリングして、前記評価値の標本平均と標本標準偏差を算出し、前記基準の設定を受け付ける基準設定部、
をさらに備える請求項１から請求項４の何れか１項に記載の制御装置。
前記基準設定部は、サンプリングした前記評価値が正規分布に従うように計算時間を定めて前記評価値をサンプリングする、
請求項５に記載の制御装置。
時間毎の電力料金単価の情報を取得する料金情報取得部と、
室外機と１台以上の室内機を備える空調機について、１台以上の前記空調機の運転状態に関する運転状態情報を取得する運転情報取得部と、
前記運転状態情報と所定の予測モデルとに基づいて、未来の所定期間における前記空調機の運転状態を予測する予測部と、をさらに備え、
前記評価部は、前記予測部が予測した運転状態に基づく前記所定期間における前記室内機が空調対象とするエリア毎の快適性の評価値と、前記時間毎の電力料金単価の情報に基づく前記所定期間における前記空調機の動作に必要な電力料金の評価値と、に基づく前記空調機の全てを対象とする運転に関する評価値を算出し、
前記制御情報算出部は、前記評価値が所定の基準を満たすような前記制御情報を算出する、
請求項１から請求項６の何れか１項に記載の制御装置。
請求項７に記載の制御装置と、
前記制御装置が算出する制御情報に基づいて運転する１台以上の空調機と、
を備える空調制御システム。
空調機の制御情報を算出するステップと、
前記制御情報に基づいて前記空調機を運転した場合の評価値を算出するステップと、
サンプリング後の前記評価値の分布が正規分布と近似するように設定された所定の時間をかけて前記制御情報をランダムに生成し、該制御情報に基づいて複数の前記評価値をサンプリングするステップと、
前記評価値を算出するステップにおいて算出された判定対象の前記評価値が、前記サンプリングされた複数の前記評価値の分布を近似して得られる正規分布に基づいて設定される所定の範囲の値か否かによって、前記評価値が所定の基準を満たすかどうかを判定するステップと、
を有する制御情報の算出方法。
制御装置のコンピュータを、
空調機の制御情報を算出する手段、
前記制御情報に基づいて前記空調機を運転した場合の評価値を算出する手段、
サンプリング後の前記評価値の分布が正規分布と近似するように設定された時間をかけてランダムに生成された前記制御情報に基づいて複数の前記評価値をサンプリングする手段、
前記評価値を算出する手段が算出した判定対象の前記評価値が、前記サンプリングされた複数の前記評価値の分布を近似して得られる正規分布に基づいて設定される所定の範囲の値か否かによって、前記評価値が所定の基準を満たすかどうかを判定する手段、
として機能させるためのプログラム。