(57)【特許請求の範囲】
【請求項１】
外気冷房を含む複数の熱源機器を備えた空調熱源システムの運転を制御する熱源制御装置であって、
所定の計測周期で計測された前記空調熱源システムの空調負荷および外気条件の計測値を時系列データとして格納する計測値データベースと、
前記計測値の時系列データに基づいて、前記計測周期の所定回数先までの所定期間における前記空調負荷および前記外気条件の予測値を時系列データとして求める時系列予測部と、
前記予測値の時系列データに基づいて、前記外気条件と対応付けられた前記複数の熱源機器のモデルパラメータを推定するパラメータ推定部と、
前記モデルパラメータに基づいて、所定の制御周期ごとの前記複数の熱源機器の運転状態を状態変数とし、少なくとも前記空調熱源システムの前記所定期間における運用コストの最小化を目的関数として定式化された混合整数計画問題を解き、少なくとも次の制御周期における前記状態変数を求める演算部と、
前記演算部の演算結果に基づいて前記複数の熱源機器の運転を制御する制御部と、
を有することを特徴とする熱源制御装置。
【請求項２】
前記演算部は、直前の制御周期における演算結果を利用して前記混合整数計画問題を解くことを特徴とする請求項１に記載の熱源制御装置。
【請求項３】
前記状態変数は、前記複数の熱源機器のそれぞれの停止または起動を表す０−１変数、および前記複数の熱源機器のそれぞれの出力を表す連続変数を含み、
前記演算部は、分枝限定法を用いて前記混合整数計画問題を解き、
前記分枝限定法においては、前記演算部は、前記０−１変数を０または１に固定する優先順位を決定するための指標値を求め、当該求めた指標値を、前記混合整数計画問題を解くのに利用するとともに、次の制御周期における前記指標値の初期値として利用することを特徴とする請求項２に記載の熱源制御装置。
【請求項４】
前記外気条件は、外気温度および外気湿度を含み、
前記パラメータ推定部は、前記モデルパラメータが前記外気温度および前記外気湿度と対応付けられたパラメータテーブルを参照して前記モデルパラメータを求めることを特徴とする請求項１ないし請求項３の何れかに記載の熱源制御装置。
【請求項５】
請求項４に記載の熱源制御装置と、
前記空調熱源システムと、
前記外気温度および前記外気湿度を計測する温湿度計と、
を備えることを特徴とする空調システム。
【請求項６】
前記複数の熱源機器は、冷凍機および冷却塔をさらにを含むことを特徴とする請求項５に記載の空調システム。
【請求項７】
外気冷房を含む複数の熱源機器を備えた空調熱源システムの運転を制御する制御部と、
所定の計測周期で計測された前記空調熱源システムの空調負荷および外気条件の計測値を時系列データとして格納する記憶部と、
を備えるコンピュータに、
前記計測値の時系列データに基づいて、前記計測周期の所定回数先までの所定期間における前記空調負荷および前記外気条件の予測値を時系列データとして求める時系列予測処理と、
前記予測値の時系列データに基づいて、前記外気条件と対応付けられた前記複数の熱源機器のモデルパラメータを推定するパラメータ推定処理と、
前記モデルパラメータに基づいて、所定の制御周期ごとの前記複数の熱源機器の運転状態を状態変数とし、少なくとも前記空調熱源システムの前記所定期間における運用コストの最小化を目的関数として定式化された混合整数計画問題を解き、少なくとも次の制御周期における前記状態変数を求め、前記制御部に入力する演算処理と、
を実行させることを特徴とする熱源制御プログラム。
【請求項８】
前記外気条件は、外気温度および外気湿度を含み、
前記記憶部は、前記モデルパラメータが前記外気温度および前記外気湿度と対応付けられたパラメータテーブルを予め格納しており、
前記パラメータ推定処理は、前記パラメータテーブルを参照して前記モデルパラメータを求めることを特徴とする請求項７に記載の熱源制御プログラム。
【請求項９】
外気冷房を含む複数の熱源機器を備えた空調熱源システムの運転を制御する熱源制御方法であって、
所定の計測周期で計測された前記空調熱源システムの空調負荷および外気条件の計測値の時系列データに基づいて、前記計測周期の所定回数先までの所定期間における前記空調負荷および前記外気条件の予測値を時系列データとして求め、
前記予測値の時系列データに基づいて、前記外気条件と対応付けられた前記複数の熱源機器のモデルパラメータを推定し、
前記モデルパラメータに基づいて、所定の制御周期ごとの前記複数の熱源機器の運転状態を状態変数とし、少なくとも前記空調熱源システムの前記所定期間における運用コストの最小化を目的関数として定式化された混合整数計画問題を解き、少なくとも次の制御周期における前記状態変数を求め、
求めた前記状態変数に基づいて前記複数の熱源機器の運転を制御することを特徴とする熱源制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、熱源制御装置、空調システム、熱源制御プログラム、および熱源制御方法に関する。
【背景技術】
【０００２】
空調熱源システムを運用する場合、省エネルギーや温室効果ガス削減のため、運用コストを最小化するように作成された運転スケジュールに基づいて、熱源機器の運転制御が行われている。例えば、特許文献１では、混合整数（線型）計画問題を定式化して解くことによって、蓄熱槽を含む熱源システムの最適運転スケジュールを決定する熱源運転支援制御方法が開示されている。
【０００３】
このようにして、混合整数計画法などの数理計画法を用いて運転スケジュールを作成し、空調熱源システムの運用コストを最小化するように熱源機器の運転制御を行うことができる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００４−３１７０４９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
近年、さらなる省エネルギー化を図るため、熱源機器として、自然エネルギーを活用する外気冷房を備えた空調熱源システムも知られている。
【０００６】
しかしながら、このような空調熱源システムでは、外気冷房の冷却能力は、温度や湿度などの外気条件から直接影響を受け、大きく変動してしまう。さらに、外気条件の変動に追従するため、熱源機器の運転制御を行う制御周期をより短くする必要がある。そのため、特許文献１の熱源運転支援制御方法を、外気冷房を備えた空調熱源システムに適用することは困難である。
【課題を解決するための手段】
【０００７】
前述した課題を解決する主たる本発明は、外気冷房を含む複数の熱源機器を備えた空調熱源システムの運転を制御する熱源制御装置であって、所定の計測周期で計測された前記空調熱源システムの空調負荷および外気条件の計測値を時系列データとして格納する計測値データベースと、前記計測値の時系列データに基づいて、前記計測周期の所定回数先までの所定期間における前記空調負荷および前記外気条件の予測値を時系列データとして求める時系列予測部と、前記予測値の時系列データに基づいて、前記外気条件と対応付けられた前記複数の熱源機器のモデルパラメータを推定するパラメータ推定部と、前記モデルパラメータに基づいて、所定の制御周期ごとの前記複数の熱源機器の運転状態を状態変数とし、少なくとも前記空調熱源システムの前記所定期間における運用コストの最小化を目的関数として定式化された混合整数計画問題を解き、少なくとも次の制御周期における前記状態変数を求める演算部と、前記演算部の演算結果に基づいて前記複数の熱源機器の運転を制御する制御部と、を有することを特徴とする熱源制御装置である。
【０００８】
本発明の他の特徴については、添付図面及び本明細書の記載により明らかとなる。
【発明の効果】
【０００９】
本発明によれば、外気条件を考慮して混合整数計画法を用いることにより、外気冷房を備えた空調熱源システムの運転制御を的確に行うことができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の一実施形態における空調システム全体の構成の概略および熱源制御装置の構成を示すブロック図である。
【図２】時系列予測部１３、パラメータ推定部１４、および混合整数計画演算部１５の機能をコンピュータに実現させるためのプログラムの動作の概略を説明するフローチャートである。
【図３】混合整数計画演算部１５の動作（混合整数計画演算処理）を説明するフローチャートである。
【発明を実施するための形態】
【００１１】
本明細書および添付図面の記載により、少なくとも以下の事項が明らかとなる。
【００１２】
＝＝＝空調システム全体の構成の概略および熱源制御装置の構成＝＝＝
以下、図１を参照して、本発明の一実施形態における空調システム全体の構成の概略および熱源制御装置の構成について説明する。
【００１３】
図１に示されている空調システムは、空調熱源システム２の運転制御を行って、対象となる施設内や室内の温度・湿度などの環境を調整するシステムであり、例えば、消費電力の一層の削減が求められているデータセンタなどで用いられる。また、当該空調システムは、空調熱源システム２以外に、熱源制御装置１および温湿度計３を含んで構成されている。
【００１４】
空調熱源システム２は、一例として、ＮＲ台の冷凍機Ｒ１ないしＲＮＲ、ＮＣ台の冷却塔Ｃ１ないしＣＮＣ、およびＮＡ台の外気冷房Ａ１ないしＡＮＡを備えている。さらに、空調熱源システム２は、搬送機器として、各冷凍機および各冷却塔の共有ポンプＰＭを備えている。そして、各外気冷房および共有ポンプＰＭは、冷媒配管などを介して冷却対象（負荷）Ｌと接続されている。
【００１５】
熱源制御装置１は、計測値データベース１１、パラメータテーブル１２、時系列予測部１３、パラメータ推定部１４、混合整数計画演算部１５、および制御部１６を含んで構成されている。
【００１６】
計測値データベース１１には、温湿度計３から、時刻ｋ＝ｋ０（現時刻）における外気温度（計測値）ｔｋ０および外気湿度（計測値）ｍｋ０が入力されている。また、計測値データベース１１には、現時刻ｋ０における空調熱源システム２の空調負荷（計測値）ｙｋ０も入力されている。なお、計測値データベース１１には、所定の計測周期で計測されたこれらの計測値が時系列データとして格納され、当該計測値の時系列データは、時系列予測部１３に適宜読み出されている。
【００１７】
時系列予測部１３からは、所定期間における外気温度、外気湿度、および空調負荷の予測値が時系列データとして出力されている。ここで、所定期間は、計測周期を１ステップとして、現時刻の次のステップ（ｋ＝ｋ０＋１）から所定のステップ数先のステップまでの期間であり、当該所定のステップ数をＫとすると、例えば、外気温度の予測値の時系列データは、ｔｋ０＋１ないしｔｋ０＋Ｋと表される。そして、これらの予測値の時系列データは、パラメータ推定部１４に入力されている。
【００１８】
パラメータテーブル１２には、外気温度および外気湿度と対応付けられた空調熱源システム２のモデルパラメータが予め格納されており、当該パラメータテーブル１２は、パラメータ推定部１４に適宜参照されている。また、後述するように、パラメータ推定部１４からは、パラメータ推定部１４によって推定された空調熱源システム２のモデルパラメータが反映されたベクトルｆｋ０、ｈｋ０、および行列Ｇｋ０が出力されている。そして、これらの行列およびベクトルは、混合整数計画演算部１５に入力されている。
【００１９】
混合整数計画演算部１５からは、混合整数計画演算部１５によって求められた、所定の制御周期ごとの熱源機器の運転状態を表す状態変数が出力されている。後述するように、本実施形態では、制御周期を計測周期と同一周期とし、混合整数計画演算部１５からは、現時刻の次のステップにおける状態変数（ベクトル）ｕｋ０＋１およびδｋ０＋１が出力される。また、これらの状態変数は、制御部１６に入力され、制御部１６からは、空調熱源システム２の熱源機器（冷凍機、冷却塔、外気冷房）および搬送機器（共有ポンプＰＭ）の制御信号ＣＮＴが出力されている。
【００２０】
＝＝＝混合整数計画問題の定式化＝＝＝
本実施形態の熱源制御装置１は、予め空調熱源システム２の運用コストの最小化問題を混合整数計画問題として定式化しておき、制御周期ごとにこの問題を解くことによって、各熱源機器および搬送機器の運転制御を行うための状態変数を求めている。そこで、まず、一例として図１に示した空調熱源システム２に対して、このような混合整数計画問題を定式化する方法について説明する。
【００２１】
目的は空調熱源システム２（熱源機器および搬送機器）の運用コストの最小化であり、目的関数を時刻ｋ＝１，２，…，Ｋにおける空調熱源システム２の運用コストＪとし、以下の式（１）のように定義する。  
0005874297.tif 000002
このとき、時刻ｋにおける運用コスト特性を、以下の式（２）のように、時変パラメータを有する線形離散事象モデルで表す。  
0005874297.tif 000003
ここで、冷凍機、冷却塔、および外気冷房の出力（冷却能力）を、それぞれ
0005874297.tif 000004
の連続変数で表し、時刻ｋにおける制御入力（ベクトル）ｕｋを、以下の式（４）のように定義する。  
0005874297.tif 000005
また、冷凍機、冷却塔、外気冷房、および共有ポンプＰＭの起動／停止（１／０）を、それぞれ
0005874297.tif 000006
の０−１変数で表し、時刻ｋにおける起動／停止（ベクトル）δｋを、以下の式（６）のように定義する。  
0005874297.tif 000007
なお、図１に示した空調熱源システム２においては、共有ポンプの台数ＮＰ＝１となっている。
【００２２】
次に、制約条件として、時刻ｋにおける空調負荷ｙｋ∈Ｒｎと、冷却能力ｕｋおよび起動／停止δｋとは、各時刻ｋにおいて需給バランスを満たす必要があるため、以下の制約条件式（７）を満たす。  
0005874297.tif 000008
また、各熱源機器の冷却能力は以下の上下限制約式（８）を満たす。  
0005874297.tif 000009
さらに、熱源機器の急激な起動／停止の切り替えを防ぐため、以下の連続起動時間制約式（９）
0005874297.tif 000010
および以下の連続停止時間制約式（１０）
0005874297.tif 000011
を課す。
【００２３】
一方、搬送機器である共有ポンプＰＭは、それを利用する熱源機器の何れかが起動する際、同時に起動される必要があるため、以下の運用制約式（１１）を満たす。  
0005874297.tif 000012
ここで、共有ポンプＰＭを利用する熱源機器の冷却能力の和の最大値を
0005874297.tif 000013
と定義している。ただし、太字の“ｅｋ”は太字の“ｕｋ”と同じサイズの単位ベクトルである。
【００２４】
以上を整理すると、空調熱源システム２の最適制御問題、すなわち、運用コストの最小化問題は、以下の式（１３）ないし式（１５）のように、混合（０−１）整数計画問題として定式化される。  
0005874297.tif 000014
ここで、
0005874297.tif 000015
とする。また、
0005874297.tif 000016
とすれば、
0005874297.tif 000017
である。なお、式（１８）において、“diag”はdiagonal matrixの略号であり、対角成分のみ保持する対角行列を意味する。また、式（１６）と式（１８）における、斜体の“0”は零行列であり、該当する箇所（行と列成分）が全て0により構成される。
【００２５】
＝＝＝熱源制御装置の動作の概略＝＝＝
以下、図２を適宜参照して、本実施形態における熱源制御装置の動作の概略について説明する。
【００２６】
計測値データベース１１には、温湿度計３から現時刻ｋ０における外気温度ｔｋ０および外気湿度ｍｋ０が入力され、空調熱源システム２から現時刻ｋ０における空調負荷ｙｋ０が入力される。そして、計測値データベース１１は、これらの計測値を所定の計測周期（サンプリング周期）の時系列データとして格納する。
【００２７】
前述したように、パラメータテーブル１２には、空調熱源システム２のモデルパラメータが予め格納されている。より具体的には、式（１３）ないし式（１５）のように定式化された混合整数計画問題では、パラメータテーブル１２は、時変パラメータａ，ｂ，ｃ，ｄを外気温度ｔおよび外気湿度ｍと対応付けた外気温度・湿度テーブルとなっている。
【００２８】
熱源制御装置１のうち、時系列予測部１３、パラメータ推定部１４、および混合整数計画演算部１５の機能は、例えば、計測値データベース１１およびパラメータテーブル１２を格納する記憶部と、制御部１６を備えるコンピュータによって実現することができる。図２は、時系列予測部１３、パラメータ推定部１４、および混合整数計画演算部１５に相当する機能をコンピュータに実現させるためのプログラムの動作の概略を示している。
【００２９】
プログラムの処理が開始されると（Ｓ１）、まず、記憶部に格納された計測値データベース１１を参照して、前述した所定期間（時刻ｋ０＋１から時刻ｋ０＋Ｋまで）における外気温度、外気湿度、および空調負荷を予測し、当該Ｋステップ分の予測値を時系列データとして求める（Ｓ２）。したがって、当該Ｓ２の処理は、時系列予測部１３に相当する機能（時系列予測処理）を実現する。
【００３０】
時系列予測処理においては、例えば、公知の記憶ベース推論を用いることができる。当該記憶ベース推論に基づく予測方法では、例えば、現時刻ｋ０からＫ’ステップ前までの計測値列（直近計測値列）を計測値データベース１１から抽出し、さらに、直近計測値列と類似する同じ時間帯の計測値列（類似計測値列）を、最小二乗法などを用いて抽出する。そして、抽出した類似計測値列の次のＫステップ分の計測値列を直近計測値列と連続するように結合することによって、予測計測値列（予測値の時系列データ）を求めることができる。なお、Ｋ’＝Ｋとしてもよい。
【００３１】
次に、記憶部に格納されたパラメータテーブル１２を参照して、空調熱源システム２のモデルパラメータを推定する（Ｓ３）。したがって、当該Ｓ３の処理は、パラメータ推定部１４に相当する機能（パラメータ推定処理）を実現する。
【００３２】
より具体的には、時系列予測処理によって予測された時刻ｋ＝ｋ０＋１，…，ｋ０＋Ｋにおける外気温度（予測値）ｔｋおよび外気湿度（予測値）ｍｋに対応する外気温度・湿度テーブルを選択し、混合整数計画問題の時変パラメータａ，ｂ，ｃ，ｄを決定する。そして、これらの時変パラメータを式（１６）および式（１８）に当てはめることによって、推定されたモデルパラメータが反映された、時変パラメータベクトルｆｋ０および行列Ｇｋ０を求めることができる。一方、時変パラメータベクトルｈｋ０は、時刻ｋにおける空調負荷（予測値）ｙｋを当てはめることによって求めることができる。
【００３３】
最後に、パラメータ推定処理によって時変パラメータが決定された混合整数計画問題を解き、少なくとも、現時刻の次のステップにおける状態変数（ベクトル）ｕｋ０＋１およびδｋ０＋１を決定し（Ｓ４）、処理を終了する（Ｓ５）。したがって、当該Ｓ４の処理は、混合整数計画演算部１５に相当する機能（混合整数計画演算処理）を実現する。なお、混合整数計画演算処理についての詳細な説明は後述する。
【００３４】
制御部１６は、混合整数計画演算処理によって求められた状態変数に基づいて制御信号ＣＮＴを出力し、空調熱源システム２の各熱源機器および搬送機器の運転制御を行う。そして、以上の処理を制御周期ごとに繰り返すことによって、熱源制御装置１は、空調熱源システム２の運用コストを最小化するように各機器の運転制御を行うことができる。
【００３５】
＝＝＝混合整数計画演算部の動作＝＝＝
以下、図３を参照して、混合整数計画演算部１５の動作（混合整数計画演算処理）について説明する。本実施形態において、混合整数計画演算部１５は、図３に示すような分枝限定法を用いて、定式化された混合整数計画問題を解いている。
【００３６】
まず、前回の演算結果を利用できない最初の制御周期における動作について説明する。なお、この場合の処理手順は、公知の分枝限定法と同様である。
【００３７】
パラメータ推定処理によって混合整数計画問題の時変パラメータが決定されると、混合整数計画演算部１５は、当該問題Ｐをリストに追加するとともに、最適値（目的関数値），最適解の暫定値をそれぞれｚ＝∞，ｘ＝φ（空集合）として、混合整数計画演算処理を開始する（Ｓ４０１）。また、前回の演算結果がなく（Ｓ４０２：ＮＯ）、リストに問題Ｐがあれば（Ｓ４０４：ＮＯ）、リストから問題ＰをＰｉとして選択する（Ｓ４０５）。そして、選択された問題Ｐｉの緩和問題Ｐｉ’を解き、最適解（緩和解）ｘｉ’、およびその場合の最適値ｚ’を求める（Ｓ４０６）。
【００３８】
緩和問題Ｐｉ’を解いた結果、Ｓ４０７ないしＳ４０９の判定結果がいずれもＮＯとなった場合には、状態変数のうち１つの０−１変数を０または１に固定した部分問題Ｐｉ１，Ｐｉ２を生成し（分枝操作）、リストに追加する。この際、後述の擬コストを算出する（Ｓ４１０）。そして、再びＳ４０４の判定を行い、リストが空でなければ（Ｓ４０４：ＮＯ）、リストの順序に従って新たに問題Ｐｉを選択し（Ｓ４０５）、その緩和問題Ｐｉ’を解く（Ｓ４０６）。
【００３９】
一方、緩和問題Ｐｉ’を解いた結果、Ｓ４０７ないしＳ４０９の何れかの判定結果がＹＥＳとなった場合には、問題Ｐｉについてこれ以上部分問題を生成して解の探索を行う必要がないため、問題Ｐｉをリストから除去し、探索範囲を限定する（限定操作）。
【００４０】
緩和問題Ｐｉ’が許容解（実行可能解）を持たず、実行不可能である場合（Ｓ４０７：ＹＥＳ）には、問題Ｐｉも許容解を持たないため、問題Ｐｉをリストから除去し（Ｓ４１２）、再びＳ４０４の判定を行う。
【００４１】
また、緩和問題Ｐｉ’の最適値ｚ’は、問題Ｐｉ、および問題Ｐｉについてのすべての部分問題の最適値の下界値となるため、下界値テストと呼ばれるＳ４０８の判定を行う。すなわち、ｚ’≧ｚが成立する場合（Ｓ４０８：ＹＥＳ）には、問題Ｐｉやその部分問題からは、既に得られている最適値の暫定値ｚ≠∞よりも良い（小さい）最適値は得られないため、問題Ｐｉをリストから除去し（Ｓ４１２）、再びＳ４０４の判定を行う。
【００４２】
また、緩和問題Ｐｉ’の最適解ｘｉ’が整数（０−１）条件を満たし、問題Ｐｉの許容解である場合（Ｓ４０９：ＹＥＳ）には、問題Ｐｉの最適解でもあるため、暫定値をそれぞれｚ＝ｚ’，ｘ＝ｘｉ’として更新する（Ｓ４１１）。そして、問題Ｐｉをリストから除去し（Ｓ４１２）、再びＳ４０４の判定を行う。
【００４３】
以上の処理を繰り返し、リストが空になると（Ｓ４０４：ＹＥＳ）、その時点の暫定値ｘを最適解として求め、混合整数計画演算処理を終了する（Ｓ４１３）。そして、制御部１６は、当該最適解ｘに含まれる状態変数（ベクトル）ｕｋ０＋１およびδｋ０＋１に基づいて、空調熱源システム２の運転制御を行う。
【００４４】
ところで、Ｓ４１０の分枝操作において部分問題を生成する際、０または１に固定する０−１変数（分枝変数）を何れとするかは、混合整数計画演算処理の効率に大きな影響を与える。本実施形態では、一例として、０−１変数を０または１に固定する優先順位を決定する際に、擬コストと呼ばれる指標値を用いることとする。
【００４５】
ここで、ｊ番目の０−１変数の緩和解ｘｊ’を０，１に固定した場合の擬コストをそれぞれｈｄ（ｊ），ｈｕ（ｊ）とすると、目的関数値（最適値）の増分は、それぞれ
0005874297.tif 000018
と推定される。そして、例えば、それぞれの差の絶対値を
0005874297.tif 000019
と求め、ｄｊが最大となる０−１変数を優先的に固定して部分問題を生成する。当該部分問題のうち、推定される目的関数値が小さい方の部分問題について先に解の探索を行うことによって、新たな暫定値ｚ＝ｚ’が得られた場合（Ｓ４０９：ＹＥＳ、Ｓ４１１）には、他方の部分問題については、下界値テストによってただちに探索が終了する（Ｓ４０８：ＹＥＳ、Ｓ４１２）可能性が高くなる。
【００４６】
次に、２回目以降の制御周期における動作について説明する。  
前述したように、計測値データベース１１には、所定の計測周期で計測された外気温度・湿度などの計測値が入力され、混合整数計画演算部１５は、当該計測値を考慮に入れて定式化された混合整数計画問題を、分枝限定法を用いて制御周期ごとに解いている。このように、オンラインで繰り返し分枝限定法を用いて混合整数計画問題を解く場合には、時間的に隣接する制御周期における問題の類似性を活用することができる。
【００４７】
本実施形態では、２回目以降の制御周期で、前回（直前の制御周期）の演算結果がある場合（Ｓ４０２：ＹＥＳ）には、直前の制御周期において求めた擬コストを現制御周期における擬コストの初期値として設定し（Ｓ４０３）、解の探索を開始する（Ｓ４０４〜）。したがって、何れの０−１変数を０または１に固定して部分問題を生成すべきかの優先順位情報を初期段階から利用することができるため、解の探索を効率よく行うことができる。
【００４８】
また、前回の演算結果を利用する他の方法として、外気温度・湿度や空調負荷の変化が微小である場合には、混合整数計画問題の最適解の変化も微小であると想定されるため、新たに混合整数計画問題を解くことなく、前回の演算結果を利用する構成としてもよい。例えば、外気温度・湿度や空調負荷の変化量が閾値以下の場合に、制御部１６は、直前の制御周期において求めた最適解ｘに含まれる、時刻ｋ０＋２における状態変数（ベクトル）ｕｋ０＋２およびδｋ０＋２に基づいて、空調熱源システム２の運転制御を行う構成とすることもできる。
【００４９】
このように、前回の演算結果を積極的に利用することによって、混合整数計画演算部１５の計算量を低減することができる。
【００５０】
前述したように、熱源制御装置１において、予め空調熱源システム２の運用コストの最小化問題を混合整数計画問題として定式化しておき、計測値の時系列データが格納された計測値データベース１１を参照して所定期間における外気温度・湿度および空調負荷を予測し、これらの予測値に基づいて推定されたモデルパラメータが反映された混合整数計画問題を制御周期ごとに解き、少なくとも、現時刻の次のステップにおける熱源機器の運転状態を表す状態変数を求めることによって、外気冷房を含む熱源機器の外気温度・湿度への依存性を考慮した混合整数計画問題を定式化して解くことができる。そのため、混合整数計画法を用いて、外気冷房を備えた空調熱源システム２の運用コストを最小化するように各機器の運転制御を的確に行うことができる。
【００５１】
また、混合整数計画問題を解く際に、直前の制御周期における演算結果を利用することによって、計算量を低減することができる。
【００５２】
また、各熱源機器の起動／停止を０−１変数で表した混合（０−１）整数計画問題を分枝限定法を用いて解く際に、擬コストと呼ばれる指標値を求めて、０−１変数を０または１に固定する優先順位を決定するとともに、求めた擬コストを次の制御周期における擬コストの初期値として利用することによって、何れの０−１変数を０または１に固定して部分問題を生成すべきかの優先順位情報を初期段階から利用することができるため、解の探索を効率よく行うことができる。そのため、空調熱源システム２の運用コストの最小化問題について、最適性が保証された厳密な最適解を実用的な計算時間で得ることが可能となる。
【００５３】
また、時変パラメータを外気温度・湿度と対応付けたパラメータテーブル１２を参照することによって、外気温度・湿度の予測値に対応する外気温度・湿度テーブルを選択し、混合整数計画問題の時変パラメータを決定することができる。
【００５４】
また、図１に示した空調システムにおいて、空調熱源システム２および温湿度計３から計測値データベース１１に、それぞれ現時刻ｋ０における空調負荷および外気温度・湿度の計測値を入力することによって、熱源制御装置１は、当該計測値を考慮に入れて定式化された混合整数計画問題を解き、各熱源機器の運転制御を行うための状態変数を求めることができる。
【００５５】
また、熱源機器として、自然エネルギーを活用する外気冷房のほか、冷凍機および冷却塔を組み合わせた空調熱源システム２を構築し、外気温度・湿度に応じて各熱源機器の運転制御を行うことによって、外気条件の変動に追従して需給バランスを満たしつつ省エネルギー化を図ることができる。
【００５６】
また、熱源制御装置１の時系列予測部１３、パラメータ推定部１４、および混合整数計画演算部１５に相当する機能をコンピュータに実現させるためのプログラムにおいて、予め空調熱源システム２の運用コストの最小化問題を混合整数計画問題として定式化しておき、記憶部に格納された計測値の時系列データに基づいて所定期間における外気温度・湿度および空調負荷を予測し、これらの予測値に基づいて推定されたモデルパラメータが反映された混合整数計画問題を制御周期ごとに解き、少なくとも、現時刻の次のステップにおける熱源機器の運転状態を表す状態変数を求めることによって、外気冷房を含む熱源機器の外気温度・湿度への依存性を考慮した混合整数計画問題を定式化して解くことができる。そして、制御部１６は、求めた状態変数に基づいて、外気冷房を備えた空調熱源システム２の運用コストを最小化するように各機器の運転制御を的確に行うことができる。
【００５７】
また、時変パラメータが外気温度・湿度と対応付けられたパラメータテーブル１２を記憶部に予め格納しておくことによって、外気温度・湿度の予測値に対応する外気温度・湿度テーブルを選択し、混合整数計画問題の時変パラメータを決定することができる。
【００５８】
また、予め空調熱源システム２の運用コストの最小化問題を混合整数計画問題として定式化しておき、計測値の時系列データを参照して所定期間における外気温度・湿度および空調負荷を予測し、これらの予測値に基づいて推定されたモデルパラメータが反映された混合整数計画問題を制御周期ごとに解き、少なくとも、現時刻の次のステップにおける熱源機器の運転状態を表す状態変数を求めることによって、外気冷房を含む熱源機器の外気温度・湿度への依存性を考慮した混合整数計画問題を定式化して解き、求めた状態変数に基づいて、外気冷房を備えた空調熱源システム２の運用コストを最小化するように各機器の運転制御を的確に行うことができる。
【００５９】
なお、上記実施形態は、本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明は、その趣旨を逸脱することなく、変更、改良され得るとともに、本発明にはその等価物も含まれる。
【符号の説明】
【００６０】
１          熱源制御装置
２          空調熱源システム
３          温湿度計
１１        計測値データベース
１２        パラメータテーブル
１３        時系列予測部
１４        パラメータ推定部
１５        混合整数計画演算部
１６        制御部
Ｌ          冷却対象（負荷）
ＰＭ        共有ポンプ
Ｒ１〜ＲＮＲ  冷凍機
Ｃ１〜ＣＮＣ  冷却塔
Ａ１〜ＡＮＡ  外気冷房
【図１】
0005874297.tif 000020
【図２】
0005874297.tif 000021
【図３】
0005874297.tif 000022