(57)【要約】【課題】空調制御システム及び空調制御方法において、電子機器を冷却するファンの電力を削減すること。【解決手段】発熱部品１４ａを備えた電子機器１４と、電子機器１４に冷却風Cを供給するファン１２ａと、ファン１２ａの回転数を制御する制御部３０とを有し、制御部３０は、電子機器１４の消費電力の上昇率vp_maxが第１の閾値θ1よりも小さい第１の場合に、発熱部品１４ａの将来の温度を予測することにより、発熱部品１４ａの温度が許容範囲に収まるようにファン１２ａの回転数を制御し、上昇率vp_maxが第１の閾値θ1以上となる第２の場合に、第１の場合におけるよりも大きな回転数u0となるようにファン１２ａの回転数の制御を切り替えることを特徴とする空調制御システムによる。【選択図】図４
2015162098.tif 000002
【特許請求の範囲】
【請求項１】
発熱部品を備えた電子機器と、
前記電子機器に冷却風を供給するファンと、
前記ファンの回転数を制御する制御部とを有し、
前記制御部は、
前記電子機器の消費電力の上昇率が第１の閾値よりも小さい第１の場合に、前記発熱部品の将来の温度を予測することにより、前記発熱部品の前記温度が許容範囲に収まるように前記ファンの回転数を制御し、
前記上昇率が前記第１の閾値以上となる第２の場合に、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数の制御を切り替えることを特徴とする空調制御システム。
【請求項２】
前記制御部は、
前記上昇率が前記第１の閾値以上となってから所定の時間が経過するまで、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数を制御し、
前記所定の時間が経過したとき、前記発熱部品の将来の温度を予測することにより、前記発熱部品の温度が前記許容範囲に収まるように前記ファンの回転数の制御を切り替えることを特徴とする請求項１に記載の空調制御システム。
【請求項３】
前記所定の時間は、前記上昇率が前記第１の閾値以上となるときの該上昇率に応じて長く設定されることを特徴とする請求項２に記載の空調制御システム。
【請求項４】
前記電子機器に電力を供給する電源線と、
前記電源線に接続され、前記消費電力を測定する電力センサとを更に有し、
前記制御部は、前記電力センサで測定された前記消費電力から前記上昇率を算出することを特徴とする請求項１乃至請求項３のいずれか１項に記載の空調制御システム。
【請求項５】
電子機器に冷却風を供給するファンの回転数を制御する制御部が、前記電子機器の消費電力の上昇率が第１の閾値以上か否かを判断する処理と、
前記上昇率が前記第１の閾値よりも小さいと判断された第１の場合に、前記制御部が、前記電子機器が備える発熱部品の将来の温度を予測することにより、前記発熱部品の前記温度が許容範囲に収まるように前記ファンの回転数を制御する処理と、
前記上昇率が前記第１の閾値以上と判断された第２の場合に、前記制御部が、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数の制御を切り替える処理と、
を有することを特徴とする空調制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、空調制御システム及び空調制御方法に関する。
【背景技術】
【０００２】
データセンタにおいては、サーバ等の複数の電子機器にジョブが分配され、各電子機器においてジョブが実行される。電子機器にはCPU(Central Processing Unit)等の発熱部品が設けられているが、大量のジョブを処理するとCPU温度も上昇し、CPUが熱暴走するおそれがある。
【０００３】
このようなCPU温度の上昇を防止するために、データセンタには電子機器を冷却するための空調機が設けられる。空調機にはいくつかのタイプがある。例えば、データセンタの外気をそのままデータセンタに取り込んで冷却風を生成し、その冷却風で電子機器を冷却する空調機がある。また、データセンタ内で空気を循環させつつ、熱交換器でその空気を冷却して冷却風を冷却するパッケージエアコンも空調機として使用されることがある。
【０００４】
いずれのタイプの空調機においても、冷却風を生成するためのファンが設けられる。データセンタの消費電力を抑制するには、そのファンの消費電力をいかにして削減するかが鍵となる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００１−２８４８６８号公報
【特許文献２】特開平５−９５０６３号公報
【特許文献３】特開２０１１−６５４４４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
空調制御システム及び空調制御方法において、電子機器を冷却するファンの電力を削減することを目的とする。
【課題を解決するための手段】
【０００７】
以下の開示の一観点によれば、発熱部品を備えた電子機器と、前記電子機器に冷却風を供給するファンと、前記ファンの回転数を制御する制御部とを有し、前記制御部は、前記電子機器の消費電力の上昇率が第１の閾値よりも小さい第１の場合に、前記発熱部品の将来の温度を予測することにより、前記発熱部品の前記温度が許容範囲に収まるように前記ファンの回転数を制御し、前記上昇率が前記第１の閾値以上となる第２の場合に、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数の制御を切り替える空調制御システムが提供される。
【０００８】
また、その開示の別の観点によれば、電子機器に冷却風を供給するファンの回転数を制御する制御部が、前記電子機器の消費電力の上昇率が第１の閾値以上か否かを判断する処理と、前記上昇率が前記第１の閾値よりも小さいと判断された第１の場合に、前記制御部が、前記電子機器が備える発熱部品の将来の温度を予測することにより、前記発熱部品の前記温度が許容範囲に収まるように前記ファンの回転数を制御する処理と、前記上昇率が前記第１の閾値以上と判断された第２の場合に、前記制御部が、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数の制御を切り替える処理とを有する空調制御方法が提供される。
【発明の効果】
【０００９】
以下の開示によれば、電子機器の消費電力の上昇率が第１の閾値以上となった場合、発熱部品の温度上昇を見越してファンの回転数を大きくするので、発熱部品の温度の上昇にファンの回転数が追従でき、発熱部品の冷却不足を抑制できる。
【００１０】
その結果、発熱部品の冷却不足を防止する目的でファンの回転数を常に高めておく必要がなくなり、ファンの消費電力を低減することが可能となる。
【図面の簡単な説明】
【００１１】
【図１】図１は、本実施形態に係るデータセンタの模式上面図である。
【図２】図２は、本実施形態に係るデータセンタの模式側面図である。
【図３】図３は、本実施形態に係る空調制御システムのハードウェア構成図である。
【図４】図４は、本実施形態に係る空調制御システムの制御部の機能ブロック図である。
【図５】図５は、本実施形態に係る空調制御システムのモデル予測部の機能ブロック図である。
【図６】図６は、本実施形態に係る空調制御方法について示すフローチャートである。
【図７】図７（ａ）〜（ｃ）は、比較例に係る調査結果を示すグラフである。
【図８】図８（ａ）〜（ｃ）は、比較例に係る調査結果を示すグラフである。
【図９】図９は、比較例と本実施形態の各々についてファンユニットの電力量を調査して得られたグラフである。
【発明を実施するための最良の形態】
【００１２】
本実施形態の説明に先立ち、本願発明者が検討した事項について説明する。
【００１３】
サーバ等の電子機器に搭載されるCPUには、CPUが熱暴走を起こさないための上限温度が設定されている。理想的には、CPUを冷却するためのファンの回転数をなるべく低下させてファンの冷却能力を落とすことで、CPU温度が前述の上限温度よりもわずかに低い状態を維持し、ファンの消費電力を抑制するのが好ましい。
【００１４】
しかしながら、CPU温度はサーバの稼働率によって常に変動しており、CPU温度が急激に上昇する場合がある。このような場合に前述のようにファンの回転数が低下していると、CPU温度がその上限温度を超えてしまう。
【００１５】
これを防ぐためには、例えば、前述の上限温度よりも十分に低い温度をCPU温度の目標温度として設定し、CPU温度がその目標温度になるようにファンの回転数を増大させればよい。但し、これではファンの消費電力も増加するため、データセンタの省エネルギ化を実現することはできない。
【００１６】
以下に、ファンの消費電力を抑制することが可能な本実施形態について説明する。
【００１７】
（本実施形態）
図１は、本実施形態に係るデータセンタの模式上面図である。
【００１８】
なお、以下では、外気を利用してサーバ等の電子機器を冷却するモジュール型データセンタを例にして説明する。
【００１９】
図１に示すように、このデータセンタ１は、直方体形状のコンテナ１０と、そのコンテナ１０内に配置されたファンユニット１２と、複数のラック１３とを有する。ラック１３の各々には、サーバ等の電子機器１４が収容されている。
【００２０】
また、コンテナ１０の相互に対向する二つの面のうちの一方には吸気口１０ａが設けられており、他方には排気口１０ｂが設けられている。
【００２１】
ファンユニット１２は複数のファン１２ａを備えており、各ファン１２ａが回転することにより、吸気口１０ａからコンテナ１０内に外気が取り込まれ、その外気から冷却風Cが生成される。
【００２２】
冷却風Cは、各電子機器１４を冷却した後、前述の排気口１０ｂから排出される。
【００２３】
図２は、このデータセンタ１の模式側面図である。
【００２４】
なお、図２において、図１で説明したのと同じ要素には図１におけるのと同じ符号を付し、以下ではその説明を省略する。
【００２５】
図２に示すように、冷却ユニット１２とラック１３との間の空間はコールドアイル２２として供せられ、ラック１３と排気口１０ｂとの間の空間はホットアイル２３として供せられる。
【００２６】
前述のコールドアイル２２の上には仕切り板１５が設けられている。そして、この仕切り板１５と、ラック１３の上面と、コンテナ１０の天井面とにより流路２４が画定される。
【００２７】
その流路２４の端部には開閉自在なダンパ１７が設けられる。冬季等のように外気の温度が低い場合には、そのダンパ１７が開くことにより、各電子機器１４から排出された暖かな冷却風Cが流路２４を通ってファンユニット１２の上流側に導かれ、各電子機器１４の過剰冷却を防止できる。
【００２８】
図３は、このデータセンタ１を冷却するための空調制御システムのハードウェア構成図である。
【００２９】
図３に示すように、空調制御システム１００には、電源線１６、制御部３０、目標温度設定部３１、発熱温度測定部３２、エアー温度測定部３３、回転数設定部３５、及び制御パラメータ設定部３６が設けられる。
【００３０】
電源線１６は、複数のコンセント１９ａを備えたテーブルタップ１９に接続されており、そのテーブルタップ１９から各電子機器１４の電源線１８に電力が分配される。
【００３１】
テーブルタップ１９には、複数の電源線１８ごとに電力センサ３４が設けられている。各電力センサ３４は、電源線１８を流れる電流等に基づいて、電源線１８が接続されている電子機器１４の各々の消費電力vPを測定し、その測定結果を制御部３０に伝送する。
【００３２】
制御部３０は、ファンユニット１２の各ファン１２ａの回転数を制御するものであり、例えば、マイコン、FPGA(Field Programmable Gate Array)、及びPLC(Programmable Logic Controller)のいずれかである。これらに代えて、制御部３０が行うプログラムを汎用コンピュータにロードすることにより、ソフトウェア的に制御部３０を実現してもよい。
【００３３】
各電子機器１４には、CPU等の発熱部品１４ａが設けられる。発熱温度測定部３２は、発熱部品１４ａと同一チップ内に搭載されており、電子機器１４内に設けられた不図示の通信機を介して発熱部品１４ａの温度を制御部３０に伝送する。その伝送のプロトコルとしては例えばUDP(User Datagram Protocol)があるが、これ以外のプロトコルを用いて電子機器１４と制御部３０との間で通信を行ってもよい。
【００３４】
なお、発熱部品１４ａと発熱温度測定部３２としてそれぞれ別々の半導体パッケージを用い、これらを密着させることで発熱部品１４ａの温度を測定するようにしてもよい。
【００３５】
また、この例では前述のようにテーブルタップ１９の電力センサ３４で発熱部品１４ａの凡その消費電力が推定できる。そのため、発熱部品１４ａの消費電力を測定する目的で発熱部品１４ａにアクセスする必要がなく、発熱部品１４ａがジョブを実行するのを妨げることがない。
【００３６】
各発熱部品１４ａの温度y0は、ファンユニット１２で生成される冷却風Cの風量によって変動する。以下では温度y0をCPU温度y0とも呼ぶ。
【００３７】
この例では、目標温度設定部３１においてCPU温度y0の目標温度rを定め、後述のようにCPU温度y0がその目標温度rに近づくように、ファン１２ａの回転数を設定する。
【００３８】
また、エアー温度測定部３３は、電子機器１４に供給される直前の冷却風Cの温度vTを測定し、その測定結果を制御部３０に伝送する。
【００３９】
なお、エアー温度測定部３３は電子機器１４ごとに設けられており、前述の温度vTは電子機器１４ごとに測定される。
【００４０】
そして、回転数設定部３５は各ファン１２ａの最大の回転数u0を記憶する。その回転数u0については後述する。
【００４１】
また、制御パラメータ設定部３６は、各ファン１２ａの回転数を制御するときに用いる各種の制御パラメータを記憶する。
【００４２】
図４は、制御部３０の機能ブロック図である。
【００４３】
なお、図４において、図３で説明したのと同じ要素には図３におけるのと同じ符号を付し、以下ではその説明を省略する。
【００４４】
制御部３０は、温度収集部４１、最高温度検出部４２、平均温度演算部４３、電力収集部４４、最大電力検出部４５、モデル予測部４６、及び回転数固定部４７を有する。そのほかに、制御部３０は、電力上昇率演算部４８、継続時間演算部４９、切替判定部５０、切替操作部５１、及び操作量記憶部５２を有する。
【００４５】
温度収集部４１は、発熱温度測定部３２で測定した各発熱部品１４ａのCPU温度y0を収集し、それらを最高温度検出部４２に出力する。
【００４６】
そして、最高温度検出部４２は、複数のCPU温度y0のうちで最も高い温度である最高温度yrealを検出する。
【００４７】
一方、平均温度演算部４３は、エアー温度測定部３３で測定された電子機器１４ごとの冷却風Cの温度vTを受け、それらの温度vTの平均温度v1を算出する。
【００４８】
電力収集部４４は、電力センサ３４で測定した各電子機器１４の消費電力vPを収集し、それらを最大電力検出部４５と電力上昇率演算部４８とに出力する。
【００４９】
そして、最大電力検出部４５は、複数の消費電力vPのうちで最も大きな最大電力v2を算出する。
【００５０】
また、モデル予測部４６は、最高温度yreal、平均温度v1、及び最大電力v2に基づいて将来のCPU温度を予測することにより、CPU温度が許容範囲に収まるようなファン１２ａの回転数を算出する。その算出方法については後述する。
【００５１】
一方、回転数固定部４７は、回転数設定部３５に予め記憶されていた各ファン１２ａの最大の回転数u0に各ファン１２ａの回転数を固定する。
【００５２】
なお、本実施形態では、ある時刻においてはモデル予測部４６と回転数固定部４７のいずれか一方のみで各ファン１２ａの回転数の制御を行う。モデル予測部４６と回転数固定部４７のどちらを選択するかについては後述する。
【００５３】
そして、上記の最大の回転数u0としては、モデル予測部４６の制御下にある各ファン１２ａの回転数よりも大きい回転数が採用される。
【００５４】
また、回転数固定部４７によって各ファン１２ａの回転数の制御が行われている場合、回転数固定部４７による制御の残り時間を継続時間tconと呼ぶ。
【００５５】
更に、電力上昇率演算部４８は、複数の電子機器１４の各々の消費電力vPの上昇率を計算し、そのうちで最も大きな最大電力上昇率vp_maxを特定する。なお、消費電力vPの上昇率は、消費電力vPの時間微分で定義される。その最大電力上昇率vp_maxは、後段の継続時間演算部４９と切替判定部５０に出力される。
【００５６】
継続時間演算部４９は、上記の最大電力上昇率vp_maxに基づいて、前述の継続時間tconを算出する。継続時間tconの具体的な算出方法については後述する。
【００５７】
切替判定部５０は、モデル予測部４６と回転数固定部４７のどちらに基づいて各ファン１２ａの回転数を制御するかを判定する。その判定には、前述の最大電力上昇率vp_maxや継続時間tconが用いられる。
【００５８】
切替操作部５１は、切替判定部５０の判定結果に基づいて、各冷却ファン１２ａの回転数の制御主体をモデル予測部４６と回転数固定部４７のいずれか一方に切り替える。
【００５９】
なお、各冷却ファン１２ａの回転数の制御主体がモデル予測部４６の場合、切替操作部５１は、各冷却ファン１２ａの現在の回転数uを操作量記憶部５２に出力する。
【００６０】
図５は、前述のモデル予測部４６の機能ブロック図である。
【００６１】
図５に示すように、モデル予測部４６は、予測モデル６０、補正部６１、コスト関数６２、及び最適化部６３を有する。
【００６２】
これらのうち、予測モデル６０は、前述の最高温度yreal、平均温度v1、最大電力v2、及び各ファン１２ａの回転数uに基づいて、発熱部品１４ａの将来の温度の予測値
【００６３】
【数１】
2015162098.tif 000003
を算出する。
【００６４】
また、補正部６１は、この予測値
【００６５】
【数２】
2015162098.tif 000004
を補正することにより、この予測値が発熱部品１４ａの実際の温度に近づくようにする。
【００６６】
そして、コスト関数６２は、予測値
【００６７】
【数３】
2015162098.tif 000005
と目標温度rとの差分を重み付けしてなる関数であり、その形については後述する。
【００６８】
更に、最適化部６３は、現在から将来にわたる所定の時間範囲について、後述の制約条件を満たし、かつコスト関数６２の値Jを最小とする冷却ファン１２ａの操作量Δuを算出する。
【００６９】
次に、本実施形態に係る空調制御方法について説明する。
【００７０】
図６は、本実施形態に係る空調制御方法について示すフローチャートである。
【００７１】
このフローチャートは、制御部３０により所定の制御周期Δtで実行される。制御周期Δtは、このフローチャートを実行する周期を表す整数である。
【００７２】
まず、ステップS11において、複数の発熱部品１４ａの各々のCPU温度y0と、電子機器１４ごとに測定した冷却風Cの温度vTと、電子機器１４の各々の消費電力vPとを制御部３０が取得する。
【００７３】
次いで、ステップS12に移り、最高温度検出部４２が、複数のCPU温度y0のうちで最も高い温度である最高温度yrealを検出する。
【００７４】
また、これと共に、平均温度演算部４３が前述の温度vTの平均温度v1を算出する。
【００７５】
更に、最大電力検出部４５が、複数の消費電力vPのうちで最も大きな消費電力である最大電力v2を算出する。
【００７６】
次に、ステップS13に移る。
【００７７】
本ステップにおいては、まず、目標温度設定部３１において予め定められているCPU温度y0の目標温度rをモデル予測制御部４６が取得する。
【００７８】
また、回転数固定部４７が、回転数設定部３５に予め記憶されている各ファン１２ａの最大の回転数u0を取得する。
【００７９】
更に、制御部３０が、制御パラメータ設定部３６に記憶されている各種の制御パラメータを取得する。なお、その制御パラメータの具体的な内容については後述する。
【００８０】
次いで、ステップS14において、電力上昇率演算部４８が、複数の電子機器１４の各々の消費電力vPの上昇率を計算し、そのうちで最も大きな最大電力上昇率vp_maxを監視する。
【００８１】
次に、ステップS15に移り、継続時間演算部４９が以下のように継続時間tcon(k)を算出する。
【００８２】
前述のように、継続時間tcon(k)は、回転数固定部４７によって各ファン１２ａの回転数の制御が行われている場合における回転数固定部４７による制御の残り時間であって、時点kに依存する。なお、時点kは、制御部３０が図６のフローチャートを実行するのが何回目であるのかを示す整数である。
【００８３】
まず、次の式（１）に従い、継続時間tcon(k-1)をΔtだけ減ずる。
【００８４】
【数４】
2015162098.tif 000006
Δtは、前述の制御周期であって、この例では１秒とする。
【００８５】
式（１）は、制御部３０が図６のフローチャートを一回実行するたびに、回転数固定部４７によって制御が行われる残りの時間がΔtだけ少なくなることを示す。
【００８６】
更に、このtcon(k-1)とは別に、次の式（２）に従ってtcon(k)を算出する。
【００８７】
【数５】
2015162098.tif 000007
式（２）に示すように、本実施形態では、最大電力上昇率vp_maxが第１の閾値θ1以上の場合に継続時間tcon(k)を最大電力上昇率vp_maxに比例させる。第１の閾値θ1の意義については後述する。また、βrは正の比例定数である。これは、最大電力上昇率vp_maxが大きくなるほど発熱部品１４ａの温度も高くなるので、回転数固定部４７によって各ファン１２ａの回転数を最大にしておく継続時間を長くすべきとの考えによる。
【００８８】
なお、最大電力上昇率vp_maxが第１の閾値θ1よりも小さい場合は継続時間tcon(k)を０とする。
【００８９】
そして、次の式（３）に従い、式（１）のtcon(k-1)と式（２）のtcon(k)のうちで大きい方をk時点の継続時間tcon(k)として採用する。
【００９０】
【数６】
2015162098.tif 000008
このように二つの値tcon(k-1)、tcon(k)のうちで大きい方を採用するのは、安全性を考慮してのことであり、二つのうちで大きな方を採用しておけば電子部品１４ａの冷却不足が防止できるためである。
【００９１】
次に、ステップS16に移る。
【００９２】
ステップS16以降では、モデル予測部４６と回転数固定部４７のいずれか一方を選択して各ファン１２ａの制御をする。
【００９３】
前述のように、回転数固定部４７は、各ファン１２ａを最大の回転数u0で回転させる。よって、各ファン１２ａの制御主体として回転数固定部４７を選択すると、モデル予測部４６を選択した場合よりも積極的に発熱部品１４ａを冷却できる。
【００９４】
そこで、この例では発熱部品１４ａを積極的に冷却する必要がある場合には回転数固定部４７を選択し、そうでない場合にはモデル予測部４６を選択する。
【００９５】
また、どちらを選択するかの判断材料として、前述の最大電力上昇率vp_maxと継続時間tcon(k)のいずれかを採用する。
【００９６】
最大電力上昇率vp_maxが急激に上昇すると将来的に発熱部品１４ａの温度が上昇する蓋然性が高いので、最大電力上昇率vp_maxは発熱部品１４ａを積極的に冷却すべきかどうかの判断材料となり得る。
【００９７】
また、継続時間tcon(k)は、式（２）のように最大電力上昇率vp_maxに比例するので、最大電力上昇率vp_max自身と同様に発熱部品１４ａを積極的に冷却すべきかどうかの判断材料となり得る。
【００９８】
そこで、ステップS16においては、切替判定部５０が、最大電力上昇率vp_maxが第１の閾値θ1以上（vp_max≧θ1）か否かについてと、継続時間tcon(k)が第２の閾値θ2よりも大きい（tcon(k)＞θ2）か否かについての判断をする。
【００９９】
各閾値θ1、θ2は、発熱部品１４ａを積極的に冷却しなければならない程度に最大電力上昇率vp_maxや継続時間tcon(k)が高いかどうかを判断する指標であって、実験やシミュレーションにより予め求めておけばよい。
【０１００】
ここで、上記の二つの条件vp_max≧θ1、tcon(k)＞θ2のいずれか一方でも成立すればステップS18に移り、二つの条件の両方とも成立しない場合にはステップS17に移る。
【０１０１】
まず、ステップS17に移った場合について説明する。
【０１０２】
ステップS17においては、モデル予測部４６が発熱部品１４ａの将来の温度を予測することにより、発熱部品１４ａの温度が許容範囲に収まるように各ファン１２ａの回転数を制御する。この制御は、以下のように予測モデルを用いて行われる。
【０１０３】
その予測モデルの一般式は、次の式（４）で表される。
【０１０４】
【数７】
2015162098.tif 000009
式（４）の左辺の
【０１０５】
【数８】
2015162098.tif 000010
は、k+1時点における発熱部品１４ａの温度の予測値である。
【０１０６】
また、式（４）の右辺のu(k)はk時点での冷却ファン１２ａの回転数であり、v1(k)はk時点における冷却風Cの平均温度である。そして、v2(k)は、k時点における各電子機器１４の消費電力のうちで最も大きい最大電力である。
【０１０７】
なお、回転数u(k)は、操作量記憶部５２に記憶されているものを使用する。
【０１０８】
本実施形態では、式（４）の一般式を次の式（５）、（６）のように具体化する。
【０１０９】
【数９】
2015162098.tif 000011
【０１１０】
【数１０】
2015162098.tif 000012
但し、
【０１１１】
【数１１】
2015162098.tif 000013
である。
【０１１２】
式（５）、（６）におけるx(k)は、k時点における状態変数であり、n次元（nは自然数）のベクトルである。また、Aはn×nの行列、Buはn次元ベクトル、Bvはn×２次元の行列、Cはn次元ベクトルである。
【０１１３】
なお、A、Bu、Bv、Cの各々の成分の値は、将来の発熱部品１４ａの温度の予測値
【０１１４】
【数１２】
2015162098.tif 000014
を最も良く近似するように、実験データによりシステム同定により求め得る。システム同定の方法としては、例えば、予測誤差法や部分空間同定法がある。
【０１１５】
また、発熱部品１４ａの温度の動的特性を表現する物理モデルの微分方程式が導出できる場合は、その微分方程式をテーラー展開して線形化することによってもA、Bu、Bv、Cの各々の成分の値を求め得る。
【０１１６】
更に、nは、モデルの次数ndとむだ時間dtで定まり、n=nd+dtとなることが知られている。この理由については、後述の参考例において説明する。
【０１１７】
なお、むだ時間dtは、ファン１２ａの回転数のむだ時間であり、ファン１２ａに対して回転数を指示してからファン１２ａがその回転数に達するまでの時間として定義される。この例では、その時間を四捨五入することによりむだ時間dtを整数値とし、例えばdt＝１２秒とする。
【０１１８】
また、上記では状態空間モデルを用いたが、モデルの表現方法は、重回帰式モデルであってもよいし、マップ関数のようなデータであってもよい。
【０１１９】
次に、補正部６１が、次の式（８）に基づいて、k+1時点における発熱部品１４ａの温度の予測値
【０１２０】
【数１３】
2015162098.tif 000015
を補正し、補正後の予測値y(k+1|k)を算出する。
【０１２１】
【数１４】
2015162098.tif 000016
式（８）とこれ以降の式では、p時点の変数αをq時点の情報で算出する場合にα(p|q)と書く。
【０１２２】
式（８）の右辺第１項の
【０１２３】
【数１５】
2015162098.tif 000017
は、補正前の発熱部品１４ａの温度のk+1時点での予測値である。
【０１２４】
そして、式（８）の右辺第２項は補正項である。その補正項において、yreal(k)はk時点での発熱部品１４ａの実際の最高温度であって、ステップS12で取得した値である。そして、その補正項におけるy(k|k-1)はk時点での発熱部品１４ａの温度の予測値である。k時点においては実際の値がyreal(k)−y(k|k-1)だけ予測値から乖離していたので、yreal(k)−y(k|k-1)を式（８）の右辺に加えることで、k+１時点において予測値が実際の値から乖離するのを防止できる。
【０１２５】
ここで、将来区間pを導入する。将来区間pとは、現在から発熱部品１４ａの温度を予測する将来までの時間を表す整数である。以下では将来区間pを例えば１００とする。
【０１２６】
そして、ファン１２ａの回転数の変化量Δuを次の式（９）で定義する。
【０１２７】
【数１６】
2015162098.tif 000018
式（９）において、iは、将来区間pをp個に等分するインデックスである。
【０１２８】
式（９）から理解されるように、変化量Δu(k+i|k)は、k+i時点におけるファン１２ａの回転数u(k+i|k)と、k+iよりも一つ前のk+i-１時点におけるファン１２ａの回転数u(k+i-1|k)で定義される。
【０１２９】
また、式（９）における各回転数u(k)としては、操作量記憶部５２に記憶されているものを使用し得る。
【０１３０】
なお、ファン１２ａの回転数は制御部３０が操作するものであるため、以下では変化量Δuのことを操作量Δuと呼ぶ場合もある。
【０１３１】
式（９）のインデックスiを用いると、前述の式（５）、（６）、（８）は、それぞれ次の式（１０）〜（１２）のように表現できる。
【０１３２】
【数１７】
2015162098.tif 000019
【０１３３】
【数１８】
2015162098.tif 000020
【０１３４】
【数１９】
2015162098.tif 000021
次に、各パラメータの許容範囲を以下のように定義する。
【０１３５】
【数２０】
2015162098.tif 000022
【０１３６】
【数２１】
2015162098.tif 000023
【０１３７】
【数２２】
2015162098.tif 000024
式（１３）は、発熱部品１４ａの最高温度yの許容範囲を定めたものであって、yminはその許容範囲の下限温度を示し、ymaxはその上限温度を示す。これらの値は特に限定されないが、上限温度ymaxとしては発熱部品１４ａが熱暴走を起こさない８０℃程度の温度を採用し得る。また、下限温度yminは例えば２０℃である。
【０１３８】
式（１４）は、ファン１２ａの操作量Δuの許容範囲を定めたものである。その許容範囲の最小値Δuminと最大値Δumaxは、一回の操作でファン１２ａの回転数を変化させ得る限界値である。
【０１３９】
そして、式（１５）は、ファン１２ａの回転数uそのものの許容範囲を定めたものであり、uminはその許容範囲の下限値を示し、umaxはその上限値を示す。
【０１４０】
各パラメータy、Δu、uは、それぞれ式（１３）〜（１５）の制約条件を受けることになる。
【０１４１】
また、上記の制約条件の他に、本実施形態では操作量Δuに対して更に次の式（１６）の制約条件を課す。
【０１４２】
【数２３】
2015162098.tif 000025
式（１６）は、k+m時点以降において操作量Δuが０になることを示す。これは、将来区間の終点において初めて操作量Δuが０になるのではなく、終点に向かって徐々に操作量Δuを０に近づけていくべきとの考えによる。
【０１４３】
また、mの値は特に限定されないが、この例ではmは１とする。
【０１４４】
次に、最適化部６３が、次の式（１７）で表されるコスト関数６２を呼び出す。
【０１４５】
【数２４】
2015162098.tif 000026
式（１７）において、Q、RΔu、Ruは重みを示すスカラである。式（１７）の右辺第一項は、制御対象である発熱部品１４ａの温度yを前述の目標温度rに近づける操作であり、Qはその操作の重み、すなわち目標値追従パラメータである。
【０１４６】
式（１７）の右辺第二項は、操作量uの変化量Δuを０に近づける操作であり、RΔuはその操作の重み、すなわち操作量低減パラメータである。RΔuが小さいと操作量Δuは大きく、RΔuが大きいと操作量Δuは小さくなる。
【０１４７】
式（１７）の右辺第三項は、回転数uを目標の回転数utargetに近づける操作である。目標の回転数utargetが小さいほどファン１２ａの消費電力を抑制できるので、本実施形態ではutargetを０とする。Ruは回転数を目標の回転数utargetに近づける操作の重み、すなわち操作量変動幅パラメータである。
【０１４８】
これらの制御パラメータQ、RΔu、Ruは、前述の制御パラメータ設定部３６に格納されており、最高温度検出部４２を介してモデル予測部４６が取得する。
【０１４９】
そして、最適化部６３が、次の式（１８）に基づいて、このコスト関数６３の値Jを最小にする操作量Δuの入力列を算出する。
【０１５０】
【数２５】
2015162098.tif 000027
次に、式（１８）で求まった最適入力列
【０１５１】
【数２６】
2015162098.tif 000028
の先頭要素Δuopt(k|k)を最適化部６３が抽出する。
【０１５２】
更に、最適化部６３は、k時点におけるファン１２ａの回転数u(k)を次の式（１９）から算出する。
【０１５３】
【数２７】
2015162098.tif 000029
コスト関数６２を最小にする最適化ソルバーは、遺伝的アルゴリズム（GA: Genetic Algorithm）や粒子群最適化（PSO： Particle Swarm Optimization：）等の近似解を探索するメタヒューリスティックな数値解法を用いることもできる。但し、この例では２次計画問題を解くための逐次２次計画法（SQP: Sequential Quadratic Programming）を用いた。
【０１５４】
以上でステップS17を終了する。
【０１５５】
この後は、ステップS19に移り、制御部３０がファン１２ａの回転数を制御するための制御信号を生成し、各ファン１２ａの回転数を式（１９）のu(k)にする。
【０１５６】
次に、ステップS16からステップS18に移った場合について説明する。
【０１５７】
ステップS18においては、回転数固定部４７の制御下において、回転数設定部３５に記憶されている最大の回転数u0で各ファン１２ａを回転させる。
【０１５８】
最大の回転数u0は、ステップS17でモデル予測部４６の制御下にある各ファン１２ａの回転数よりも大きければ特に限定されず、シミュレーションや実機を用いた実験で任意に設定し得る。
【０１５９】
その実験では、例えば、予め冷却風Cの温度vTを電子機器１４の仕様で許容できる最高温度にしておく。そして、各ファン１２ａの風量を、各ファン１２ａで送風できる最大風量にしておく。この状態で、CPUの最高温度yをモニタリングしながら、少しずつファン１２ａの回転数を下げていく。そして、最高温度yが前述の上限値ymaxを超えない最大の回転数を見つけ、その回転数を最大の回転数u0とし得る。
【０１６０】
なお、このようにステップS18を実行し続けると継続時間tcon(k)が徐々に少なくなっていくので、図６のフローチャートを繰り返し実行していくうちにステップS16において条件式tcon(k)＞θ2がいずれ不成立となり、再びステップS17が行われる。最大電力上昇率vp_maxが第１の閾値θ1以上となってステップS18が開始された後、このようにステップS17が再び行われるまでの時間を所定の時間Tと呼ぶ。
【０１６１】
これによれば、最大電力上昇率vp_maxが第１の閾値θ1以上となってステップS18が開始された後、所定の時間Tが経過するまではステップS18が継続して実行され、ファン１２ａの回転数が最大の回転数u0となる。
【０１６２】
そして、上記の所定の時間Tが経過したときは、再びステップS17が行われることになる。
【０１６３】
なお、継続時間tcon(k)が長いほど、回転数固定部４７がファン１２ａの回転数を制御している時間が長くなるので、上記の所定の時間Tも長くなる。式（２）のように本実施形態では継続時間tcon(k)を最大電力上昇率vp_maxに比例させたので、結局、所定の時間Tは、最大電力上昇率vp_maxに応じて長く設定されることになる。
【０１６４】
また、ステップS16でvp_max≧θ1が成立するか否かを判定するときに使用する最大電力上昇率vp_maxは、直前のステップS14で取得する。よって、所定の時間Tの設定に使用する最大電力上昇率vp_maxは、当該最大電力上昇率vp_maxが第１の閾値θ1以上となるときの最大電力上昇率vp_maxの値となる。
【０１６５】
この後は、前述のステップS19に移り、制御部３０がファン１２ａの回転数を制御するための制御信号を生成し、各ファン１２ａの回転数を最大の回転数u0に固定する。
【０１６６】
なお、このように回転数を固定せずに、最大の回転数u0よりも大きな回転数で各ファン１２ａを回転させてもよい。
【０１６７】
以上により、本実施形態に係る空調制御方法の基本ステップを終了する。
【０１６８】
上記した本実施形態によれば、ステップS16で最大電力上昇率vp_maxが第１の閾値θ1以上（vp_max≧θ1）と判定されたときに、ステップS18〜S19において最大の回転数u0にファン１２ａの回転数を固定する。
【０１６９】
電子機器１４の消費電力が増大すると将来的に発熱部品１４ａの温度が上昇する可能性が高くなる。よって、このように電力の上昇率を判断基準にすることで、発熱部品１４ａの温度が実際に上昇するのを待たずにファン１２ａの回転数を最大の回転数u0にすることができ、最大の回転数u0に固定するタイミングが発熱部品１４ａの温度上昇に遅れるのを防止できる。
【０１７０】
この効果を確かめるために、本願発明者は次のような調査を行った。
【０１７１】
この調査では、比較例と本実施形態の各々について、複数のCPU温度のうちの最高温度yが時間と共にどのように変化するかが調べられた。
【０１７２】
比較例の調査結果を図７（ａ）に示す。
【０１７３】
この比較例では、図４の回転数固定部４７を設けずに、常にモデル予測部４６でファン１２ａの回転数の制御を行った。なお、CPUの目標温度rは８９℃に設定した。
【０１７４】
また、図７（ｂ）はファン１２ａの回転数uの時間変化を示すグラフであり、図７（ｃ）は最大電力v2の時間変化を示すグラフである。
【０１７５】
図７（ａ）の点線円Aに示すように、比較例においては、発熱部品１４ａの使用率が上昇して最大電力v2が急激に増加したときに、最高温度yが目標温度rを超えるオーバーシュートが発生している。
【０１７６】
オーバーシュートが発生しても最高温度yが目標温度rを超えないようにするには、目標温度rを８９℃よりも低く設定する必要がある。例えば、この例ではオーバーシュート時の温度上昇が１０℃程度なので、目標温度rを８９℃よりも１０℃低い７９℃程度にしなければならない。
【０１７７】
しかし、このように目標温度rを低くすると、最高温度yが目標温度rを超えないようにするためにファン１２ａの回転数を常に高めておく必要があり、ファン１２ａの消費電力が増加してしまう。
【０１７８】
一方、図８（ａ）〜（ｃ）は、これと同じ調査を本実施形態に対して行って得られたグラフである。
【０１７９】
これらのグラフのうち、図８（ａ）は、CPU温度のうちの最高温度yの時間変化を示すグラフである。
【０１８０】
そして、図８（ｂ）はファン１２ａの回転数uの時間変化を示すグラフであり、図８（ｃ）は最大電力v2の時間変化を示すグラフである。
【０１８１】
図８（ａ）の点線円Bに示すように、本実施形態では最大電力v2が急激に増加しても、最高温度yはほぼ目標温度r以下に保たれており、比較例のような顕著なオーバーシュートは発生していない。
【０１８２】
これは、最大電力v2が急激に上昇した１０００秒の時点において、ファン１２ａの制御主体がモデル予測部４６から回転数固定部４７に切り替わり、その回転数固定部４７の制御下でファン１２ａが最大の回転数u0で回転したためと考えられる。
【０１８３】
このようにオーバーシュートが殆ど発生しないため、本実施形態ではオーバーシュートを見越して目標温度rを低く設定する必要がなく、比較例よりもファン１２ａの消費電力を抑制することが可能となる。
【０１８４】
図９は、図７及び図８の比較例と本実施形態の各々について、ファンユニット１２の電力量を調査して得られたグラフである。
【０１８５】
比較例では前述のように目標温度rを７９℃としている。一方、本実施形態の目標温度rは８９℃である。
【０１８６】
図９に示すように、本実施形態においては、ファンユニット１２の電力量が比較例のそれよりも約５２％も削減できる。これは、前述のように本実施形態では目標温度rを低くする必要がないので、各ファン１２ａの回転数を高める必要がないからである。
【０１８７】
図７〜図９の結果により、本実施形態のように最大電力上昇率vp_maxを判断基準にしてファン１２ａの回転数を最大の回転数u0にするか否かを判断することが、ファン１２ａの消費電力を低減するのに有効であることが確認できた。
【０１８８】
以上、本実施形態について詳細に説明したが、本実施形態は上記に限定されない。
【０１８９】
例えば、図１では外気から冷却風Cを生成するタイプのモジュール型データセンタ１を例示したが、そのモジュール型データセンタ１に気化式冷却装置を設け、冷却風Cの温度や湿度をその気化式冷却装置で調節するようにしてもよい。
【０１９０】
更に、モジュール型データセンタに代えて、パッケージエアコンで冷却風Cを生成するタイプのデータセンタに本実施形態を適用することにより、パッケージエアコンのファンの回転数を制御してもよい。
【０１９１】
また、データセンタの空調に代えて、発熱部を多く含む施設の空調に本実施形態を適用してもよい。
【０１９２】
（参考例）
上記した本実施形態のステップS17においては、状態変数x(k)の次元n、モデルの次数nd、及びむだ時間dtが、n=nd+dtなる関係を満たすことを説明した。以下に、その理由について説明する。
【０１９３】
まず、次の式（２０）のような離散時間の1入力1出力、モデル次数１の状態空間モデルを考える。
【０１９４】
【数２８】
2015162098.tif 000030
ここで、入力uのむだ時間が1秒、kの周期が1秒であるとき、次の式（２１）のように入力uの値を状態変数の2つ目の要素に記憶し、次の周期に1行目にシフトするように表現できる。
【０１９５】
【数２９】
2015162098.tif 000031
式（２１）の例では、状態変数の次数は、モデルの次数１とむだ時間を考慮するための１を足し合わせた２となる。
【０１９６】
また、入力uのむだ時間が2秒のときは、上記と同様に、次の式（２２）のように、状態変数の2つ目の要素、3つ目の要素、及び入力uの値をシフトさせて表現できる。
【０１９７】
【数３０】
2015162098.tif 000032
式（２２）の例では、状態変数の次数は、モデルの次数１とむだ時間を考慮するための２を足し合わせた３となる。
【０１９８】
入力uのむだ時間が3秒であるときは、上記と同様に、次の式（２３）のように、状態変数の2つ目の要素、3つ目の要素、4つ目の要素、及び入力uの値をシフトさせて表現できる。
【０１９９】
【数３１】
2015162098.tif 000033
式（２３）の例では、状態変数の次数は、モデルの次数１とむだ時間を考慮するための３を足し合わせた４となる。
【０２００】
次に、以下の式（２４）のような離散時間の1入力1出力、モデル次数２の状態空間モデルを考える。
【０２０１】
【数３２】
2015162098.tif 000034
ここで、入力uのむだ時間が1秒、kの周期が1秒であるとき、次の式（２５）のように入力uの値を状態変数の3つ目の要素に記憶し、次の周期に1行目と２行目にシフトするように表現できる。
【０２０２】
【数３３】
2015162098.tif 000035
よって、状態変数の次数は、モデルの次数２とむだ時間を考慮するための次数１を足し合わせた３となる．
また、入力uのむだ時間が２秒でkの周期が1秒のときは、次の式（２６）のように入力uの値を状態変数の3つ目の要素に記憶し、更に次の周期に4つ目の要素に記憶し、次に1行目と２行目にシフトするように表現できる。
【０２０３】
【数３４】
2015162098.tif 000036
式（２６）の例では、状態変数の次数は、モデルの次数２とむだ時間を考慮するための次数２を足し合わせた４となる。
【０２０４】
以上のように、n=nd+dtなる関係が満たされる。
【０２０５】
以上説明した各実施形態に関し、更に以下の付記を開示する。
【０２０６】
（付記１）  発熱部品を備えた電子機器と、
前記電子機器に冷却風を供給するファンと、
前記ファンの回転数を制御する制御部とを有し、
前記制御部は、
前記電子機器の消費電力の上昇率が第１の閾値よりも小さい第１の場合に、前記発熱部品の将来の温度を予測することにより、前記発熱部品の前記温度が許容範囲に収まるように前記ファンの回転数を制御し、
前記上昇率が前記第１の閾値以上となる第２の場合に、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数の制御を切り替えることを特徴とする空調制御システム。
【０２０７】
（付記２）  前記制御部は、
前記上昇率が前記第１の閾値以上となってから所定の時間が経過するまで、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数を制御し、
前記所定の時間が経過したとき、前記発熱部品の将来の温度を予測することにより、前記発熱部品の温度が前記許容範囲に収まるように前記ファンの回転数の制御を切り替えることを特徴とする付記１に記載の空調制御システム。
【０２０８】
（付記３）  前記所定の時間は、前記上昇率が前記第１の閾値以上となるときの該上昇率に応じて長く設定されることを特徴とする付記２に記載の空調制御システム。
【０２０９】
（付記４）  前記電子機器に電力を供給する電源線と、
前記電源線に接続され、前記消費電力を測定する電力センサとを更に有し、
前記制御部は、前記電力センサで測定された前記消費電力から前記上昇率を算出することを特徴とする付記１乃至付記３のいずれかに記載の空調制御システム。
【０２１０】
（付記５）  前記電子機器が複数設けられ、
前記制御部は、前記電子機器の前記消費電力として、複数の前記電子機器の消費電力のうちで最も大きな最大電力を採用することを特徴とする付記１乃至付記４のいずれかに記載の空調制御システム。
【０２１１】
（付記６）  前記発熱部品の温度と、前記電子機器に入る直前の前記冷却風の温度と、前記電子機器の電力とに基づいて、前記発熱部品の前記将来の温度を予測する予測部を更に有することを特徴とする付記１乃至付記５のいずれかに記載の空調制御システム。
【０２１２】
（付記７）  前記第１の場合におけるよりも大きな前記回転数に前記ファンの前記回転数を固定する回転数固定部を更に有することを特徴とする付記１乃至付記６のいずれかに記載の空調制御システム。
【０２１３】
（付記８）  電子機器に冷却風を供給するファンの回転数を制御する制御部が、前記電子機器の消費電力の上昇率が第１の閾値以上か否かを判断する処理と、
前記上昇率が前記第１の閾値よりも小さいと判断された第１の場合に、前記制御部が、前記電子機器が備える発熱部品の将来の温度を予測することにより、前記発熱部品の前記温度が許容範囲に収まるように前記ファンの回転数を制御する処理と、
前記上昇率が前記第１の閾値以上と判断された第２の場合に、前記制御部が、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数の制御を切り替える処理と、
を有することを特徴とする空調制御方法。
【０２１４】
（付記９）前記制御部が、前記上昇率が前記第１の閾値以上となってから所定の時間が経過するまで、前記第１の場合におけるよりも大きな回転数となるように前記ファンの回転数を制御する処理と、
前記所定の時間が経過したとき、前記制御部が、前記発熱部品の将来の温度を予測することにより、前記発熱部品の温度が前記許容範囲に収まるように前記ファンの回転数の制御を切り替える処理とを更に有することを特徴とする付記８に記載の空調制御方法。
【０２１５】
（付記１０）  前記消費電力の前記上昇率が第１の閾値よりも小さいか否かを判断する処理において、
前記電子機器に電力を供給する電源線に接続された電力センサを用いて前記消費電力を取得し、該消費電力から前記上昇率を算出することを特徴とする付記８又は付記９に記載の空調制御方法。
【符号の説明】
【０２１６】
１…データセンタ、１０…コンテナ、１０ａ…吸気口、１０ｂ…排気口、１２…ファンユニット、１２ａ…ファン、１３…ラック、１４…電子機器、１４ａ…発熱部品、１５…仕切り板、１６、１８…電源線、１７…ダンパ、１９…テーブルタップ、１９ａ…コンセント、２２…コールドアイル、２３…ホットアイル、２４…流路、３０…制御部、３１…目標温度設定部、３２…発熱温度測定部、３３…エアー温度測定部、３４…電力センサ、３５…回転数設定部、３６…制御パラメータ設定部、４１…温度収集部、４２…最高温度検出部、４３…平均温度演算部、４４…電力収集部、４５…最大電力検出部、４６…モデル予測部、４７…回転数固定部、４８…電力上昇率演算部、４９…継続時間演算部、５０…切替判定部、５１…切替操作部、５２…操作量記憶部、６０…予測モデル、６１…補正部、６２…コスト関数、６３…最適化部。
【図１】
2015162098.tif 000037
【図２】
2015162098.tif 000038
【図３】
2015162098.tif 000039
【図４】
2015162098.tif 000040
【図５】
2015162098.tif 000041
【図６】
2015162098.tif 000042
【図７】
2015162098.tif 000043
【図８】
2015162098.tif 000044
【図９】
2015162098.tif 000045