JP2013543569A - Hvacシステムに対するエネルギー最適制御決定 - Google Patents

Hvacシステムに対するエネルギー最適制御決定 Download PDF

Info

Publication number
JP2013543569A
JP2013543569A JP2013525010A JP2013525010A JP2013543569A JP 2013543569 A JP2013543569 A JP 2013543569A JP 2013525010 A JP2013525010 A JP 2013525010A JP 2013525010 A JP2013525010 A JP 2013525010A JP 2013543569 A JP2013543569 A JP 2013543569A
Authority
JP
Japan
Prior art keywords
value
actuator
change
sensor
new
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2013525010A
Other languages
English (en)
Inventor
クリフォード, シー. フェデルスピール,
ジェリー チン,
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Vigilent Corp
Original Assignee
Vigilent Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Vigilent Corp filed Critical Vigilent Corp
Publication of JP2013543569A publication Critical patent/JP2013543569A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/30Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
    • F24F11/46Improving electric energy efficiency or saving
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/50Control or safety arrangements characterised by user interfaces or communication
    • F24F11/56Remote control
    • F24F11/58Remote control using Internet communication
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/62Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
    • F24F11/63Electronic processing
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/70Control systems characterised by their outputs; Constructional details thereof
    • F24F11/72Control systems characterised by their outputs; Constructional details thereof for controlling the supply of treated air, e.g. its pressure
    • F24F11/74Control systems characterised by their outputs; Constructional details thereof for controlling the supply of treated air, e.g. its pressure for controlling air flow rate or air velocity
    • F24F11/77Control systems characterised by their outputs; Constructional details thereof for controlling the supply of treated air, e.g. its pressure for controlling air flow rate or air velocity by controlling the speed of ventilators
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0205Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
    • G05B13/026Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system using a predictor
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D23/00Control of temperature
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D23/00Control of temperature
    • G05D23/19Control of temperature characterised by the use of electric means
    • G05D23/1927Control of temperature characterised by the use of electric means using a plurality of sensors
    • G05D23/193Control of temperature characterised by the use of electric means using a plurality of sensors sensing the temperaure in different places in thermal relationship with one or more spaces
    • G05D23/1932Control of temperature characterised by the use of electric means using a plurality of sensors sensing the temperaure in different places in thermal relationship with one or more spaces to control the temperature of a plurality of spaces
    • G05D23/1934Control of temperature characterised by the use of electric means using a plurality of sensors sensing the temperaure in different places in thermal relationship with one or more spaces to control the temperature of a plurality of spaces each space being provided with one sensor acting on one or more control means
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05KPRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
    • H05K7/00Constructional details common to different types of electric apparatus
    • H05K7/14Mounting supporting structure in casing or on frame or rack
    • H05K7/1485Servers; Data center rooms, e.g. 19-inch computer racks
    • H05K7/1498Resource management, Optimisation arrangements, e.g. configuration, identification, tracking, physical location
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05KPRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
    • H05K7/00Constructional details common to different types of electric apparatus
    • H05K7/20Modifications to facilitate cooling, ventilating, or heating
    • H05K7/20709Modifications to facilitate cooling, ventilating, or heating for server racks or cabinets; for data centers, e.g. 19-inch computer racks
    • H05K7/20836Thermal management, e.g. server temperature control
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/30Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/62Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Signal Processing (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Human Computer Interaction (AREA)
  • Fluid Mechanics (AREA)
  • Fuzzy Systems (AREA)
  • Remote Sensing (AREA)
  • Thermal Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Air Conditioning Control Device (AREA)
  • Feedback Control In General (AREA)
  • Cooling Or The Like Of Electrical Apparatus (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

複数のセンサ及び複数のアクチュエータを含む環境維持システムを制御する方法、システム及び装置が提供される。アクチュエータの動作レベルは、ペナルティ関数を最適化することで決定可能である。ペナルティ関数の一部として、センサ値は基準値と比較可能である。動作レベルの最適化値は、種々の動作レベルでのアクチュエータのエネルギー使用及び種々の動作レベルでの基準値に対するセンサ値の予測差を考慮できる。予測差は、伝達モデルを使用して決定可能である。伝達モデルの精度は、予測値を測定値と比較することで決定可能である。この精度は、伝達モデルの出力から新しい動作レベルを決定する(例えば、精度に基づいて伝達モデルの出力を減衰する)ために使用可能である。

Description

本発明は、一般に、例えばサーバコンピュータを含むデータセンタに存在できるような共通の空間の温度及び湿度を調節するために使用可能な冷暖房空調設備(HVAC)システム等の環境制御システムに関する。本発明は、特に、所望の環境条件を提供しつつ、最適な方法でモジュール(例えば、エネルギーコスト及び/又は資金コスト)を効率的に動作させることに関する。
近年のデータセンタは、室内の温度、湿度及び他の変動要素を調節するためにHVACシステムを使用する。データセンタにわたり配置された多くのHVACユニットを有することは一般的である。HVACユニットは、床設置型ユニットであることが多いが、壁取付け型ユニット、ラック取付け型ユニット又は天井取付け型ユニットであってもよい。更にHVACユニットは、揚床プレナム、空気ダクトのネットワーク又はデータセンタの戸外のいずれかに冷却された空気を提供することが多い。データセンタ自体又は大きなデータセンタの大部分は、一般にオープンプラン構成を有し、すなわちデータセンタの一方の部分の空気を他方の部分の空気と分離する永続的な隔壁を有さない。従って、多くの場合、共通の空間を有するこれらのデータセンタは、複数のHVACユニットにより温度調節及び湿度調節される。
データセンタ用HVACユニットは、一般に分散制御及び独立制御で動作する。各ユニットがデータセンタからユニットに入る空気の温度及び湿度を調節しようとして動作するのは一般的である。例えばHVACユニットは、ユニットに入る空気の温度及び湿度を決定するセンサを含んでもよい。このセンサの測定値に基づいて、そのHVACを制御することにより、ユニットに入る空気の温度及び湿度をそのユニットに対する設定値と整合するように変更しようとしてユニットの動作を変更する。
信頼性を確保するため、殆どのデータセンタは過剰な数のHVACユニットを有するように設計される。オープンプラン構成により空気がデータセンタにわたり自由に流れるため、一方のユニットの動作は他方のユニットの動作に結合されてもよい。過剰なユニット及びそれらがほぼ重なり合う領域に空気を分配するという事実は、冗長性を提供する。これにより、単一のユニットが故障する場合、データセンタの機器(サーバ、ルータ等)は依然として十分な冷却能力を有することを保証する。
上述したように、HVACユニットの標準的な操作手順は、各HVACユニットへの戻り空気温度を調節することである。そのような操作は、データセンタのサーバ及び他の演算機器に入る空気の温度を認識しない。HVACユニット制御又はデータセンタのオペレータがこの情報を使用できないため、一般に余分なHVACユニットは、サーバが冷却され続けていることを保証するために常に動作する。しかし、そのような予防措置はエネルギーを無駄にする。
冗長なHVACユニットを手動でoffにしてエネルギーを節約することが可能であるが、そうすることによりその機器を過熱させる危険性が生じる。更にいくつかのデータセンタは可変負荷を有し、その場合、1日中いくつかのユニットを停止及び開始する必要があるだろう。また、データセンタの密に結合された複雑な空気流パターンにより、オペレータは、必要に応じて再始動するのに最も効率的であるHVACユニットを認識するのが困難になる。
更に現在の方法は、HVACユニットの制御方法を決定する際にユニットのコスト(例えば、エネルギーの使用量又は維持)を無視する。より複雑なアルゴリズムが試行される場合、現在の方法は更に不安定になりうる。
従って、HVACシステムを動作するコストを説明し且つ冗長なHVACユニットを説明する際に安定した方法及びシステムを提供することが望ましい。
本発明の実施形態は、センサを使用して環境維持モジュール(例えば、HVACユニット)を制御するシステム及び方法を提供できる。種々の実施形態において、システムは、環境を加熱及び/又は冷却できる。センサは、環境内の種々の場所において温度又は他の状態変数(湿度、気圧、光等)を測定でき、この温度フィードバックは、最小のエネルギーを使用しつつ温度を所望の範囲内又はそれに近い範囲に維持するようにモジュールの動作レベルを変更するためにエネルギー消費情報と共に使用可能である。これらのフィードバック制御動作は発見的であってもよく、例えば、温度をオペレータが指定した設定内又はそれ以下に維持するように、環境維持モジュールにより生成された空気の温度を変更するか、あるいは送風機の速度を減速する。フィードバック制御動作は、多くのアクチュエータが同一の制御コマンドを受信するようにこれらの制御動作に影響を与えるために従来の比例積分偏差(PID)フィードバック制御装置を利用してもよい。
一態様において、予測制御決定は、ペナルティ関数に取り入れられる状態変化への制御動作の影響(例えば、アクチュエータの動作レベルの変化)を予測するために伝達モデル(例えば、伝達行列等の非線形関数又は線形関数)を使用して行われる。コンピュータシステムは、ペナルティ関数を最小限にする制御動作変化を算出できる。一般に、ペナルティ関数を最小限にする制御動作は、全てのアクチュエータを同時に変更することを含む。変更されるアクチュエータの数及び種類を制限することが可能である。例えば状況標識は、アクチュエータが変更不可能であるようにリモートシステムにより無効にされることを示してもよい。
発見的フィードバック制御動作は、予測制御動作と組み合わせられる。これは、システムをマルチレートシステムとして構成することで達成される。マルチレートシステムにおいて、フィードバック制御動作はより高い頻度で実行され、予測制御動作はより低い頻度で実行される。一実現例において、モデル予測制御動作は、多くのアクチュエータに送出された共通のフィードバック制御動作を変調するために使用される。
別の態様において、実際に実行される制御動作変化(例えば、ペナルティ関数等による発見的フィードバック又は予測からの)は、測定された過去の予測の精度に基づいて減衰される。測定された最近の過去の予測誤差の規模が小さい場合、最小減衰量が適用される。測定された最近の過去の予測誤差の規模が増大するにつれ、減衰量は増加する。測定された最近の過去の予測誤差の規模が十分に大きい場合、制御動作は場合によっては全く変化しない。種々の実現例において、減衰は、算出された変化から伝達モデルを用いて決定されたような動作レベルに値を加算又は減算することで実行される。別の実施形態において、減衰は、小数又は0〜1の他の値(終点を含む)を掛けることで得られるものであってよい。
ペナルティ関数を減少するが必ずしもそれを最小限にしない制御動作が算出されることもある。次善の予測制御決定を使用するという決定は、ランダム化されてもよく且つ/あるいはランダム化レベルの部分最適性を有してもよい。これらの次善の予測制御動作は、学習処理が伝達モデル予測の精度を向上できる「持続励振性」データを有することを保証するために使用される。
予測制御動作は、ゆっくりと(例えば、15分おきに1回)実行可能である。動作は、次の予測制御決定が実行される前にシステムが準定常状態に達するのに十分なほどゆっくりと実行されなければならない。予測制御動作は、次の予測制御動作が実行可能になる前に十分な時間が経過した時を決定するために定常状態検出器を非同期的に使用して更に実行されてもよい。
一実施形態によると、環境維持システムを制御する方法が提供される。システムは、複数のアクチュエータ及び複数のセンサを含む。各センサは、環境の物理的状態の値を測定する。センサに対応するセンサ値Sが受信される。アクチュエータの動作レベルUへの潜在的な変化dUが決定される。コンピュータシステムは、予測されたセンサ値S’の変化dSを決定するために伝達モデルをdUに適用し、予測された変化dSから新しいセンサ値SNewを決定し、且つペナルティ関数への第1の寄与を決定するためにSNewを使用することにより、ペナルティ関数への第1の寄与を算出する。伝達モデルは、アクチュエータの動作レベルの変化と結果として得られるセンサ値の変化との関係を提供する。第1の寄与は、1つ以上のそれぞれの設定値に対する新しいセンサ値の各々の関係に基づく。コンピュータシステムは、動作レベルUへの潜在的な変化dUに対応する新しい動作レベルUnewを決定し、対応する新しい動作レベルで各アクチュエータを動作するコストを決定し、且つコストを集約して第2の寄与を取得することにより、ペナルティ関数への第2の寄与を更に算出する。ペナルティ関数の第1の値は、潜在的な変化dUに対する第1の寄与及び第2の寄与に基づいて決定される。アクチュエータの動作レベルUへの最適な変化は、動作レベルUへの複数の潜在的な変化に対するペナルティ関数の値に基づいて決定される。アクチュエータの動作レベルは、最適な変化に基づいて設定される。
別の実施形態によると、センサに対応するセンサ値Sが受信される。センサ値Sは、アクチュエータの動作レベルU1に対応する。アクチュエータの新しい動作レベルU2が受信され、この場合U2−U1=dUである。センサ値Sの変化dSが予測される。変化dSは変化dUに対応し、この場合、予測は伝達モデルに従って実行される。次に、アクチュエータは新しい動作レベルU2で動作される。新しいセンサ値Snewが後続して測定される。受信、予測、動作及び後続測定は、複数回繰り返される。伝達モデルの精度の計測値は、予測センサ値を測定センサ値と比較することで決定される。伝達モデルに基づく提案されたアクチュエータの動作レベルへの変更が受信される。提案されたアクチュエータの動作レベルへの変更は、精度計測値に基づいて調整される。アクチュエータの動作レベルは、提案された変更に対する調整に基づいて設定される。
本発明の他の実施形態は、本明細書において説明する方法と関連付けられたシステム、デバイス及びコンピュータ可読媒体に関する。
本明細書において使用されるように、環境維持モジュールは、HVACユニット、コンピュータルーム空調装置(CRAC)ユニット、あるいは制御信号を受信し且つ特定の領域において環境条件を達成又は維持するように設定される出力を提供するあらゆるデバイスであってもよい。本明細書における使用値、動作レベルは、冷却の加熱等の動作を実行するためにモジュールにより費やされている労力を測定する動作パラメータ(入力又は出力)の値であってもよい。モジュールの動作レベルは、停止した全電力で動作する入力電力レベルの割合、モジュール内の特定のデバイスの動作に対する割合又は値(例えば、送風機速度、温度設定値、湿度設定値又はバルブ位置)及び例えば基準レベルと比較されるような出力加熱、出力空気流の冷却電力の割合又は値を含んでもよい。
以下の詳細な説明及び添付の図面を参照して、本発明の本質及び利点がより適切に理解されるだろう。
図1は、本発明の一実施形態に係るデータセンタを示す平面図である。 図2は、本発明の一実施形態に係る単一の屋上ユニットにより調整されたオープンプランの建物を示す平面図である。 図3は、本発明の一実施形態に係るコンピュータルーム空気処理ユニット300を示す概略図である。 図4は、本発明の一実施形態に従って建物内の環境条件を維持する制御システム400を示すブロック図である。 図5は、本発明の一実施形態に従って複数のモジュール(例えば、CRACユニット)及びセンサを含む環境維持システムを初期化する方法500を示すフローチャートである。 図6は、本発明の一実施形態に従って高い効率で所望の範囲内にセンサ値を維持するように環境維持システムを制御する方法600を示すフローチャートである。 図7は、本発明の一実施形態に従ってアクチュエータの動作を軽減することにより環境維持システムを制御する方法700を示すフローチャートである。 図8は、本発明の一実施形態に従ってモジュールを停止することにより複数のモジュール及びセンサを含む環境維持システムを制御する方法800を示すフローチャートである。 図9は、本発明の一実施形態に従って複数のモジュール及びセンサを含むモジュールを開始し且つ停止することにより環境維持システムを制御する方法900を示すフローチャートである。 図10は、予測誤差の大きさに基づいて共分散行列の追跡を更新する式を示すグラフである。 図11は、本発明の実施形態に従ってペナルティ関数を使用して環境維持システムを制御する方法1100を示すフローチャートである。 図12は、本発明の実施形態に従って環境維持システムを制御する方法1200を示すフローチャートである。 図13は、本発明の実施形態に従って環境維持システムを制御するために使用される最適なdUを決定するシステム1300を示すブロック図である。 図14は、本発明の実施形態に従って伝達モデルの精度尺度を使用して環境維持システムを制御する方法1400を示すフローチャートである。 図15は、本発明の実施形態に従って伝達モデルの精度を使用して環境維持システムを制御する方法1500を示すフローチャートである。 図16は、本発明の実施形態に従って中点及び極値制御を実現するシステム1600を示す図である。 図17は、本発明の実施形態に従って動作の中点及び極値レベルを決定するために予測制御を使用する方法1700を示すフローチャートである。 図18は、本発明の実施形態に係るシステム及び方法と共に使用可能な例示的なコンピュータ装置を示すブロック図である。
環境(例えば、データセンタ)が十分に低温であること又は温暖であることを保証するため、標準的な動作手順は、最低限度必要とされるものを超える余分なHVACユニット(又は他の環境維持モジュール)を動作させることである。しかし、そのような予防措置はエネルギーを無駄にする。種々の実施形態により、環境にわたりセンサ(例えば、モジュールにわたるセンサ又はモジュールの外側の場所におけるセンサ)を解析し、特定のモジュールの動作レベルを安全に低下できるか及びセンサを範囲に戻すのに最適なモジュールを含む更なる動作が必要か(例えば、範囲外の測定値のために)を決定できる。そのような実施形態により、エネルギー消費を減少しつつ安定した環境を提供できる。伝達モデル(例えば、行列又は行列の集合)は、上述の決定を実行するために使用可能である。一態様において、伝達モデルは、システムのセンサ毎にモジュールの動作レベルの変化と結果として得られる温度差との関係を提供できる。伝達モデルは、センサ値を範囲内に維持しつつ、最適な動作レベルを決定するようにペナルティ関数(動作のコストを含む)を最適化するために使用可能である。
予測制御動作のみを使用することに関する1つの問題は、予測が準定常挙動に基づく場合に予測制御決定を迅速に行えないことである。これは、制御動作間に妨害がある場合、システムが十分迅速に反応しなくてもよく且つ制御性能が受け入れ不可能であってもよいことを意味する。この問題は、例えばマルチレートシステムを使用して発見的フィードバック制御動作をモデル予測制御動作と組み合わせることで修正される。マルチレートシステムにおいて、発見的フィードバック制御動作はより高い頻度で実行され、モデル予測制御動作はより低い頻度で実行される。
予測制御動作のみを使用することに関する別の問題は、予測が常に正確ではない可能性があることである。動作点が変化するにつれ、予測の精度は低下するだろう。予測子を適応可能にすることで精度を向上できるが、予測がより不正確になる動作点が変化した後の期間が存在する可能性がある。不正確な予測に基づく制御動作により、受け入れ不可能な制御性能を招く可能性がある。この問題は、次善であるが受け入れ可能な制御性能を有するのもとして既知である発見的フィードバック制御動作と予測制御動作を組み合わせることで修正される。最近の過去の予測精度尺度を使用することにより、予測制御決定に対する重要度は、システムが、予測精度が低い場合には主に発見的フィードバック制御動作に依存するが、予測精度が高い場合には最適な制御性能を達成するためにモデル予測を使用するように変調可能である。この解決方法は、システムを故障許容にするという更なる利点を有する。モデル予測が完全に失敗する場合、システムは、次善であるが受け入れ可能な制御性能の発見的フィードバック制御動作のみを使用することに戻ってよい。
特に多くの冗長なHVACユニットが空気を床下プレナムに分配する際にそれらを実行することに関する別の問題は、より少ないHVACユニットが使用された場合にHVACユニットからの排出空気温度がより高くなる可能性があることである。より多くのHVACユニットが使用されている状態で依然としてサーバから同一の量の熱を抽出している間に温度がより高くなる可能性があるため、排出空気温度はより高くなる可能性がある。換言すると、サーバから熱を抽出するより大きな空気流(より多くのHVACユニットを含む)は、HVACユニットからの排出温度がより高い可能性があることを意味する。その結果、コンクリートスラブフロア又は揚床はそれほど低温ではなく、それにより災害復旧機能を減らす。例えば、スラブ及び床が蓄冷媒体(ヒートシンク)として動作するため、低温の床を有することにより、冷却障害(例えば、建物全体に対して電力が遮断される場合)から復旧するのに使用可能な時間を増加させる。例えば電力が遮断されているためにHVACユニットが冷気を押し出していない場合、床のより低い温度によりサーバをより低温に維持できる。
I.システム概要
図1は、本発明の一実施形態に係るデータセンタの間取り図を示す。周囲の壁1は、データセンタの周囲の壁又は小売店等の商業用建物の周囲の壁であってもよい。データセンタは、複数のHVACユニット2(又は他の環境維持モジュール)、複数の環境センサ3及び複数のサーバラック4を含む。示されるように、HVACユニットは、コンピュータルーム空調装置(CRAC)ユニットである。
一実施形態において、HVACユニット2は、空気流をデータセンタに提供し、サーバラック4のサーバを冷却する単一の機器である。一態様において、HVACユニット2は、それらを通過する空気を冷却、加熱、加湿又は除湿できる。環境センサ3は、温度又は湿度等の環境パラメータを測定するデバイスである。環境センサ3は、有線通信手段又は無線通信手段(例えば、Modbus、BACnet、Wi−Fi、WiMAX、ZigBee又は他のあらゆる適当なプロトコル)等のあらゆる手段により測定値(読み取り値とも呼ばれる)を送信できる。HVACユニット2(及び本明細書において示される他のモジュール)は、データセンタ内で特定の状態を提供するように1つ以上のプロセッサを含むコンピュータシステムにより制御されてもよい。
図2は、本発明の一実施形態に係る単一の屋上ユニットにより調整されたオープンプランの建物を示す平面図である。この例において、HVACユニットは屋上ユニット(RTU)2である。周囲の壁5は、小売店等の商業用建物の外側又は内側の壁、あるいはそのような建物又は店内の空間である。示されるように、有線通信はその特定のRTUに近接するRTU2とセンサ3との間で行われるが、無線通信が更に使用されてもよい。商品棚6及び店のチェックアウトカウンター7が更に示される。
図3は、本発明の一実施形態に係るコンピュータルーム空気処理ユニット300を示す概略図である。コンピュータルーム空気処理ユニット300は、環境維持モジュールの一例である。示されるように、コンピュータルーム空気処理ユニット300は、冷水バルブ320により変調された冷水を含んでもよい冷却コイル310を有する。コンピュータルーム空気処理ユニット300は、再熱コイル330(例えば、電気コイル)及び加湿器340(例えば、赤外線加湿器)を更に有する。
一実施形態において、送風機350は、A/C誘導電動機により駆動された遠心送風機である。誘導モータは、その速度を変化させる可変速(周波数)ドライブVSD355を有してもよい。無線センサ360は戻り空気温度を測定し、無線センサ370は排出空気温度を測定し、無線制御380はVSD355を制御する。排出空気温度センサ370及び戻り空気温度センサ360は、別個の無線センサではなく無線制御380につながれたプローブであってもよい。
動作の一実施形態において、無線センサ360、370は、例えば図4の監視制御器15を含む制御コンピュータに信号を渡す無線ネットワークゲートウェイに無線装置を介して読み取り値を送出する。監視制御15は、コンピュータシステム自体であってもよい。制御コンピュータは、作動コマンドを無線制御380に中継し、可変速ドライブ355の速度を変化させる無線ゲートウェイにコマンドを送出できる。
図4は、本発明の一実施形態に従って建物内の環境条件を維持する制御システム400を示すブロック図である。この例において、制御システム400は、HVACユニット2(例えば、ユニット300)と、複数の環境センサ3と、計算を行う1つ以上のプロセッサ10を含む監視制御器15とを備える。HVACユニット2は、空間の環境を維持するのに使用されてもよい送風機、バルブ又は温度要素等に対する最終制御要素(アクチュエータとも呼ばれる)を含む。本明細書において示されるように、アクチュエータの入力及び出力はモジュールの動作レベルに対応してもよい。一態様において、監視制御器15は、減少したエネルギー量又は最小のエネルギー量を使用して安定した環境条件を提供する動作レベル(on及びoff、並びに中間の変動を含む)を有するように最終制御要素を制御できる。
モジュール(HVACユニット)
いくつかの実施形態において、監視制御器15は、各HVACユニット2の入力12へのコマンドを算出することにより複数のHVACユニット2の動作を調整できる。コマンドは、センサ3からの環境センサ読み取り値に基づいて算出される。入力12は、種々の異なるHVACユニット2及び/あるいはHVACユニット2内のデバイス又は回路に対応してもよい。
一実施形態において、HVACユニット1の入力1は、1つのアクチュエータの動作パラメータ(例えば、送風機、温度設定値、湿度設定値又はバルブ位置)に対応してもよく、HVACユニット1の入力2は、同一のHVACユニット1の異なるアクチュエータに対応してもよい。動作パラメータは種々の動作値(レベル)を有してもよく、その結果、各々が種々のエネルギー量を消費する。別の実施形態において、HVACユニット2のうちのいくつかは、動作レベルの制御に対して1つの入力しか有さない。
他の実施形態において、HVACユニット2の温度に対する設定値は、監視制御器15から更に提供されてもよい。例えば設定値は、HVACユニット2により排出された空気の所望の温度又はユニットに戻る空気の所望の温度であってもよい。他の入力は、湿度(又は加湿器コマンド)、あるいは可変周波数ドライブ(VFD)へのコマンドに対する設定値であってもよい。
一実施形態において、各HVACユニットは、各々がそのHVACユニットの1つのアクチュエータに対応する同数の入力を有する。別の実施形態において、種々のHVACユニットは異なる数のアクチュエータを有してもよい。そのような一実施形態において、センサの数は、アクチュエータの総数に関係なく同一であってもよい。部分的に、センサの数が同一のままであってもよい理由は、各センサが各アクチュエータに影響を及ぼしてもよく、且つ各アクチュエータが各センサに影響を及ぼしてもよいためである。例えば温度アクチュエータ(冷却バルブ等)は、環境が十分に低温である場合に凝固物が冷却コイル上に形成される際に発生する可能性があるような湿度に影響を及ぼしてもよい。同様に、湿度アクチュエータ(例えば、赤外線加湿器及び蒸発冷却バルブ)は、赤外線加湿器又は蒸発冷却器が湿度を上げる際に発生する可能性があるような温度に影響を及ぼす。
センサ
環境センサ3は、温度、湿度及び気圧等の環境の物理的状態の値を測定できる。環境センサ3は、例えば有線通信手段又は無線通信手段(例えば、Modbus、BACnet、Wi−Fi、WiMAX、ZigBee又は他のあらゆる適当なプロトコル)により読み取り値を監視制御器15に返送できる。センサの例には、温度センサ、湿度センサ及び気圧センサが含まれる。単一のセンサは、複数の環境条件、例えば上記の条件のうちの3つ全てを測定できてもよい。環境センサ3は、規則的なパターンに従って又はランダムに位置付けられてもよい。更に環境センサ3は、センサのクラスタにより編成されてもよく、あるいは個々に配置されてもよい。
いくつかの実施形態において、監視制御器15は、例えば関連付けられた基準値8の集合により特定されるように、温度センサ読み取り値Fが温度範囲R内にあるようにする。範囲は、単にある特定の温度未満(例えば、華氏78度未満)であってもよい。範囲は、2つの温度により更に特定されてもよい。そのような温度範囲は、要望に応じて狭くても広くてもよい。そのような範囲は加熱に更に適用されてもよい。ある特定の実施形態は、最小のエネルギー量を使用しつつ温度毎に特定の温度範囲(その全ては、温度センサ毎に異なってもよくあるいは同一であってもよい)を維持しようとしてもよい。
一実施形態において、監視制御器15は、環境センサ毎の所望の基準値8の集合をフラッシュメモリ、キャッシュ又は他の適切なメモリ等に内部に格納する。他の実施形態において、基準値8は、ディスクドライブ又は光学ドライブ等に外部に格納されてもよい。動作中、監視制御器15は、最小のエネルギー量を使用して特定の範囲で環境センサ3からの値を維持するようにHVACユニット2の動作レベルを調整する(例えば、温度範囲を超えずに動作する最小限のモジュールを有することにより)。
HVACへの入力
一実施形態において、監視制御器15は、入力12に提供され且つHVACユニット2の最終制御要素(例えば、アクチュエータ)に対して直接使用されるコマンドを算出する。入力12に送出されるこれらのコマンドは、例えば有線通信手段又は無線通信手段により提供されてもよい。これらのコマンドは、HVACユニット2のあらゆる数の動作レベルを開始、停止又は変化させてもよい。
別の実施形態において、監視制御器15は、HVACユニット2においてローカルデジタル制御器(例えば、マイクロプロセッサベースの制御を有する)により使用される入力12へのコマンドを算出する。一態様において、ユニットのローカルデジタル制御器への各入力は、ユニットのアクチュエータに対応する。次にローカルデジタル制御器は、最終制御要素に送出された最終コマンドを決定してもよい。例えばローカルデジタル制御器は、アクチュエータに対してデジタル信号をアナログ信号に変換してもよく、あるいは信号のプロトコルをアクチュエータにより使用可能になるように変換してもよい。ローカルデジタル制御器は、ローカル制御ループを介して特定の設定でアクチュエータを維持するように更に動作してもよい。従って、監視制御器15は、最終制御要素に直接命令するのではなく、ローカルデジタル制御器においてローカル制御ループの設定値に命令してもよい。
状況標識
一実施形態において、監視制御器15は、環境センサ3及び/又はHVACユニット2から状況標識13を受信する手段を有する。一態様において、状況標識13は、HVACユニット2又はセンサ3が現在動作可能であるかに関する情報を提供してもよい。別の態様において、状況標識13は、戻り空気温度、排出温度、使用されているユニットの容量の部分(例えば、割合)(これは動作レベルの一例である)及び冷水バルブ(例えば、320)がどのくらい開いているか等のHVACユニットの設定を提供してもよい。状況標識13は、説明のためにHVACユニット2及びセンサ3とは別個に示され、HVACユニット2又はセンサ3自体から実際に受信されてもよい。
一実施形態において、HVACユニット2に対する状況標識13は、HVACユニット2のローカルデジタル制御器から取得されてもよい。これらのローカルデジタル制御器は、監視制御器15によりクエリされ、ローカルデジタル制御器又はHVACユニット2が「on」又は「off」であるかを決定してもよい。ユニットが「off」である場合、そのユニットのアクチュエータに対する状況標識13は、ゼロ等のある特定の値であってもよい。
別の例において、環境センサ3は、明確に規定され且つ容易に検出されるいくつかの故障モードを有する。一態様において、一方の故障モードは「到達不可能」であり、監視制御器15のネットワークインタフェース等のゲートウェイがセンサと通信できないことを意味する。他方の故障モードは範囲外電圧(0ボルト又は1.5ボルトのいずれか)であり、0ボルトはセンサプローブが短絡回路を有することを示し、1.5ボルトはセンサプローブが開回路を有するかあるいは欠落していることを示す。これらの故障のうちのいずれかの結果、そのセンサに対してゼロの状況標識が得られるだろう。
行列等の伝達関数モデル(TM)は、環境センサに対する環境維持モジュールを増加(及び場合によっては同等に減少)する影響への尺度である。行列は、システムの全てのセンサ又はセンサの一部分のみに対して影響を及ぼしてよい。一態様において、TMの行Jの数は、環境センサ(CRAC冷却サーバラックを使用する有する実施形態に対してコールドアイルセンサ又は吸気センサとも呼ばれる)の数に等しくてよく、列の数は環境維持モジュールの数に等しくてよい。従って、一実施形態において、モジュール毎に1つの列しかない。そのような一実施形態において、モジュールのエネルギー消費の1つの尺度、すなわち動作レベルが決定される1つのパラメータしかない。別の実施形態において、モジュールに対して2つ以上の行があってもよいため、各々がモジュールの動作レベルの測定値を提供するパラメータが2つ以上ありうる。尚、行及び列は切り替わってもよい。更に「行列」という用語は、値のあらゆる索引付け可能な配列であってもよい。
II.伝達関数行列及び負荷行列の初期化
伝達関数行列(TM)は、環境センサに対する環境維持モジュールを増加(及び場合によっては同等に減少)する影響の基準である。行列は、システムの全てのセンサ又はセンサの一部分のみに対して影響を及ぼしてもよい。一態様において、TMの行Jの数は、環境センサ(CRACユニットを使用する実施形態に対してコールドアイルセンサとも呼ばれる)の数に等しくてもよく、列の数は環境維持モジュールの数に等しくてもよい。従って、一実施形態において、モジュール毎に1つの列しかない。そのような一実施形態において、モジュールのエネルギー消費の1つの基準、すなわち動作レベルが決定される1つのパラメータしかない。別の実施形態において、モジュールに対して2つ以上の行があってもよいため、各々がモジュールの動作レベルの測定値を提供する2つ以上のパラメータがあってもよい。尚、行及び列は切り替わってもよい。更に「行列」という用語は、値のあらゆる索引付け可能なアレイであってもよい。
本明細書において説明されるように、動作レベルは入力値又は出力値であってもよい。例えば、無電力の入力コマンド(例えば、電圧又はデジタル値)は0の動作レベルであってもよく、全電力の入力は100%の動作レベル又は何らかの最大値であってもよい。更に動作レベルは、送風機速度、温度設定値、湿度設定値又はバルブ位置等の特定のアクチュエータに対する入力値、あるいはそのようなアクチュエータの位置の出力測定値であってもよい。別の実施形態において、更に動作レベルは、提供される冷却又は加熱のレベル等の出力レベルであってもよい。この出力レベルは、上回ることにより100%より高い割合を提供できる設計値に対する実際のレベルの割合であってもよい。パラメータが出力値である場合、出力パラメータを変更するために使用される1つ以上の入力コマンド変数があってもよい。
図5は、本発明の一実施形態に従って複数のモジュール(例えば、CRACユニット)及びセンサを含む環境維持システムを初期化する方法500を示すフローチャートである。方法を説明する際、図4を参照する。初期化は伝達行列TMを作成することを含む。一実施形態において、TMの列は、モジュールの動作レベルを上下(例えば、開始及び停止)させ、且つ動作レベルの変化前と変化後の温度の差を利用することにより初期化される。一実施形態において、モジュールは、逐次(すなわち、一度に1つ)増減する。他の実施形態において、TMの列は、一度に複数のモジュールの動作レベルを変化させることにより初期化される。
ステップ510において、第1のセンサ値Y1(例えば、温度、湿度又は気圧)は、各センサ(例えば、センサ3からの)毎に受信される(例えば、制御器15において)。第1のセンサ値Y1は、測定コマンドの結果コンピュータ(例えば、制御器15)により能動的に測定されてもよく、あるいは送信されたメッセージを受信できるポートを介して受動的に取得されてもよい。特定のセンサJに対する第1のセンサ温度Y1はY1として書き込まれてもよい。第1のセンサ値Y1は、1回又は複数回取得されてもよい。例えば第1の値Y1は、あらゆるモジュールの動作レベルが変化する前に取得されてもよく、あるいは特定のモジュールの動作値が変化する前に毎回取得されてもよい。
第1のセンサ値Y1に加え、モジュールの1つ以上の動作パラメータの動作レベル等の他の値が記録されてもよい。戻り温度(例えば、図3のセンサ360からの)、排出温度(例えば、センサ370からの)及び流量(例えば、モジュールがVFDを有する場合に流量×VFDパーセントコマンドを設計する)等の動作パラメータ入力設定値及び測定された出力値の例は、測定され且つ制御器15と通信可能に結合されたメモリ等に格納されてもよい。
ステップ520において、1つ以上のモジュールの1つ以上の動作レベルが変化する。一実施形態において、モジュールの全ての動作レベルが変化し、その変化は1レベルだけであってもよい。別の実施形態において、モジュールの動作レベルのうちの一部のみが変化する。変化しない動作レベルは、別の時に変化してもよく、あるいは伝達行列の計算に含まれなくてもよい。種々の実施形態において、複数のモジュールからの動作レベルは同時に変化してもよい。
動作レベルの変化量は、何らかの方法で異なってもよくあるいは同等であってもよい。例えば、変更される各動作パラメータは、同一の割合(例えば、100%、50%、25%等)だけ変更されてもよい。一実現例において、100%は、パラメータに対する最大設計値に対して測定される(例えば、最大設計空気流の100%)。設計空気流より大きい空気流を実現することが可能であるため、割合は100%を上回るだろう。
特定のパラメータの動作レベルに対して実現される厳密な変化量を決定するため、変化コマンドが出力された後にパラメータが測定されてもよい。動作パラメータが出力値である場合、新しい動作レベルは、変化コマンドから直接認識されなくてもよい。例えば、変化コマンドが空気流を大きくするためのものであってもよく、且つ一般にコマンドに対応する空気流を認識するためのいくつかの校正された設定値があってもよいが、測定により空気流に対する実際の値をより正確に取得できる。更にいくつかのパラメータは、出力値の測定によることを除いて全く認識されなくてもよい。他の実施形態において、動作レベルの変化は自動的に認識されてもよい(例えば、変化がoffにするためのものである場合、特にon及びoffのみがある場合)。
ステップ530において、第2のセンサ値Y2はセンサ毎に受信される。一実施形態において、モジュールの動作レベルが変化した後にタイマ(例えば、ウェブが設定可能な期間で)を開始する。デフォルトとして、期間は15分であってもよい。一実施形態において、第2の値Y2はタイマの終了後に測定される。別の実施形態において、値は、変化コマンドが与えられた後に継続的に測定され、第2の値Y2は、測定値が準定常ステート状態になった後に格納される。例えばセンサの変更出力は、摂動後にある特定の変化率を有してもよい。変化率が閾値を下回って低下すると、準定常ステート状態が決定されてもよい。センサの出力における変化の閾値に対する絶対値は別の例である
ステップ540において、各々が対応する第1のセンサ温帯Y2と第2のセンサ値Y2との差が算出される(例えば、プロセッサ10により)。従って、使用されているN個のセンサがある場合、変更された動作パラメータ毎にY2−Y1のN個の各値があり、Jは1〜Nにわたる。この差は正又は負であってもよい。一般に冷却するために、変化が動作レベルの上昇である場合、温度(冷却のための関心センサ値)は低下し、Y2−Y1の各々は負である。更に一般に冷却するために、変化が動作レベルの低下である場合、温度は上昇し、Y2−Y1の各々は正である。しかし、これらの関係は常に当てはまらない。例えば、モジュールの戻り温度が所望の排出温度以下である場合、モジュールは、冷却能力をoffにしてもよい。従って、空気が依然として送風機により送り出されることによりいくらかの熱が空気に伝えられるため、モジュールは、冷却しないが実際には加熱する。従って、いくつかの伝達行列要素は他の反対符号を有してもよく、それは非直感的である。
ステップ550において、TM行列要素は、対応する差を使用してセンサ毎に算出される。例えば、プロセッサ10はN×K行列要素を算出してもよく、Nはセンサの数であり、Kは変更された動作パラメータの数である。一実施形態において、1つの動作パラメータが一度に変更される場合、単一の列は、行列要素毎に式を使用して一度に更新されてもよい。別の実施形態において、2つ以上の動作パラメータが一度に変更される場合、複数の列が一度に更新され、組み合わされた式(例えば、再帰最小二乗)は行列要素を更新するために使用される。
ステップ560において、それ以上の動作パラメータがまだ変更されていないかが決定される。まだ変更されていない動作パラメータがある場合、方法500は繰り返されてもよい。一実施形態において、方法500の第1の反復により1つの動作パラメータPだけ減少したと仮定する。次の反復において、ステップ520は、パラメータPのレベルを上昇させ、最後の反復の前のパラメータPの動作レベルを有することを含んでもよく、パラメータPのレベルを低下させることを含んでもよい。従って、1つの動作レベルは第1の反復に基づいて変化し、第2の動作レベルは第2の反復に基づいて変化する。他の実施形態は、全ての反復において上下する複数の動作レベルを有してもよい。以下に説明するように、そのような実施形態は、行列要素を決定するために再帰最小二乗を使用してもよい。
一実施形態において、それ以上の動作レベルが変化する必要があるかをこのように決定することは、TMのそれ以上の列が初期化される必要があるかと同等である。モジュールに対して1つの動作レベルがあり且つ動作レベルオプションがon又はoffである一実施形態において、モジュールを開始又は停止するかを決定する。そのような一実施形態において、停止する必要があるより多くのモジュールがある場合、停止したモジュールは再開されてもよく、他のモジュールを停止してTM行列の他の要素を決定してもよい。モジュールが再開される場合、スタート/ストップタイマが再開されてもよく、スタート/ストップタイマの期限が切れた後、次のモジュールに対してこのような初期化が実行される。
ステップ570において、TMの全ての行列要素が算出された後、伝達行列TMは、環境維持システムのメモリに格納されてもよい。伝達行列は、アクチュエータの動作レベルを変化させるべきかを決定するために検索されてもよい。そのような決定は、例えば方法600〜900で実行されてもよい。次に、伝達行列及びLOAD行列の行列要素を算出するある特定の実施形態を説明する。
行列要素の算出
1つの動作パラメータが一度に変更される一実施形態において、変更された動作パラメータに対応するTMの1つの列の行列要素は、第2のセンサ値が受信された後に決定される。動作パラメータがUindx(本実施形態において、モジュールと一対一の対応関係を有する)であり且つセンサインデックスがSindxである場合、行列要素TM(Sindx,Uindx)は、
Figure 2013543569
・・・式(1)
のように算出されてもよく、式中、Y2は動作パラメータUindxが変更された後のSindxに対応するセンサ値であり、Y1は動作パラメータUindxが変更される前のSindxに対応するセンサ値であり、ΔlevelはパラメータUindxの動作レベルの変化である。例えば変化がアクチュエータ毎に同一でなかった場合、正規化係数が更に使用されてもよい。以下に説明するように、正規化係数であると考えられてもよいエネルギー係数が含まれてもよい。
センサ値が温度であり且つモジュール機能が冷却するように機能する場合、例えばモジュールを遮断すること(又は他の減少)によりY2がY1より大きくなるはずであり且つ動作レベル(Δlevel)の変化は負であるため、一般にTMの値は負になる。温度差は負であるがΔlevelが正であるようなモジュール(又は他の増加)を開始することにより、同様の結果が得られる。しかし、ステップ540に対して上述したように、逆の結果が生じる可能性があり、それは非直感的である。
TMは、全ての行列要素がΔlevelの同一のユニットに対応してもよいように正規化されてもよい。例えば、特定の行列要素を算出するために使用されたレベルの実際の変化(例えば、10%、10rpm)に関係なく、行列要素が係数で乗算されてもよいため、全ての行列要素は分母に同一の値を有する。従って、後続のステップにおいて、動作レベルの変化は、行列要素毎に種々のユニットにある動作レベルの変化とは対照的に予測された温度の変化を決定するために一様に使用されてもよい。
一実施形態において、Δlevelは動作レベルの変化の割合であり、例えば動作レベルの最大容量にonにする場合は100%及び最大動作レベルからoffにする場合は−100%である。別の実施形態において、Δlevelは、伝達行列TMを作成するために使用された最小増分に対する単位である。例えば、増分が5V、100rpm又は他の値(百分率を含む)である場合、伝達行列TMが最小増分の単位であった場合、10Vの変化は2の値であり、300rpmは3の値である。別の実施形態において、伝達行列が全てのモジュール(例えば、最大電力又は送風機設定)にわたる動作の最大レベルに基づく単位で作成された場合、Δlevelは小数値であってもよい。各モジュールは異なる範囲の動作レベルを有してもよく、例えば一方の送風機は2000rpmの最大速度を有してもよく、他方の送風機は1000rpmの最大速度を有してもよい。一態様において、Δlevelは自身で正規化してもよい。例えば、変更が特定のパラメータに対して常に同一である場合、正規化係数は必要なくてもよい。
一実施形態において、動作パラメータは、モジュールが操作している加熱/冷却流の百分率容量%Capであってもよい。本実施形態において、変化によりモジュールが閉じている場合、Δlevelは%Capに対応してもよい。一態様において、全てのユニットが同一の容量で動作していない場合、%Capの値は正規化された測定値を考慮してもよい。%Capは現在の動作レベルの一例である。
種々の実施形態において、%Capは、ユニットのクエリから戻されるか、あるいは、
Figure 2013543569
・・・式(2)
のように算出され、式中、Fは停止する前の停止したモジュールの流量であり、TRは停止する前の各モジュールPの戻り温度であり、TDは動作レベルが変化する前の各モジュールPの排出温度であり、Fはユニットの設計流量であり、TR−TDは各モジュールPの設計値ΔTである。
一実施形態において、設計値ΔTは、冷却バルブが完全に開いているかあるいは好適な設定値で開いている場合の温度差であってもよい。そのような値は、更なる要素として含まれてもよい使用される冷却物質(例えば、水)の温度に依存してもよい。別の実施形態において、固定の流れの送風機を含むモジュールは、F及びFの同一の値を有する。そのような実施形態において、モジュールにより戻された%Capは、冷却バルブ(例えば、バルブ320)の設定に対応してもよい。
可変送風機を含むモジュールにおいて、Fは、100%の容量又は他の何らかの割合、あるいはユニットの好適な(例えば、最適な)動作が行われる値であってもよい。一態様において、Fの値は、停止前の送風機の速度に対する設定に対応する。他の実施形態において、加熱(又は冷却)流量FΔT(設計又は停止前)は、他の手段により決定されてもよく、あるいは単にこの値を測定するモジュールから受信されるだけであってもよい。
LOAD行列の算出
一実施形態において、LOAD行列は、システムの容量に対するモジュールの動作レベルを低下させる影響の基準を提供する。例えばLOAD行列は、他のCRACユニットのうちの1つをoffにする場合にコンピュータルームのサーバの熱負荷に対処するように増加するCRACユニットの容量の割合に厳密に関連してもよい。一実施形態において、LOADの行及び列の番号は環境維持モジュールの数に等しい。
負荷行列は伝達行列TMと同時に算出されてもよい。従って、TMの列はLOAD行列の列と同時に算出されてもよい。一実施形態において、負荷伝達関数行列(LOAD)は、
Figure 2013543569
・・・式(3)
のように算出され、式中、Δ%CapCindx,KindxはモジュールKindxを停止(又は減少)することにより発生したモジュールCindxの百分率容量の変化であり、%CapKindxは停止前のモジュールKindxの百分率容量である。LOAD(Cindx,Cindx)の値は定義上−100%に等しい。一実施形態において、LOAD(Cindx,Cindx)の値は、一般にこの値が使用されないため、算出されないかあるいは0(又は他のデフォルト値)に更に設定されてもよい。LOAD行列は、環境維持システムのメモリに格納されてもよく、動作レベルが上昇又は低下したモジュールを決定するために後で使用されてもよい。
III.エネルギー使用量を削減するためのTMの使用
伝達行列は、センサを許容可能な範囲内に維持するために使用されてもよい。伝達行列は、削減されたエネルギー量を使用しつつセンサを許容可能な範囲内に維持する動作レベルを決定するために更に使用されてもよい。エネルギーに対する影響力を決定するため、動作レベルの変化は、アクチュエータ毎に固定値であると仮定される(しかし、全て同一ではない)。このように、各アクチュエータは、互いに比較されてもよく、センサに最も影響を及ぼすアクチュエータを識別してもよい。従って、アクチュエータに伝えられたエネルギー量が最も効率的であるため、このセンサが範囲外である場合、最も影響力を有するアクチュエータは、最もエネルギー効率のよいアクチュエータとして解釈されてもよい。
いくつかの実施形態において、各アクチュエータの増加に対して費やされたエネルギー量は同一である。例えば、アクチュエータの各動作パラメータがモジュールの電力レベルであり、且つ各モジュールが同一のエネルギー効率を有する場合、エネルギー効率の差はない。差が小さい場合、エネルギー効率は同等であると更に仮定されてもよい。
他の実施形態において、アクチュエータは異なるエネルギー消費を有する。そのような状況において、センサ値(Y2−Y1)の変化は、エネルギー係数により変調されてもよい。変調された結果を比較することにより、変化するアクチュエータを決定する際のエネルギー消費を説明してもよい。例えば、所定の変化に対してより少ないエネルギーを使用するアクチュエータがセンサ値の差に対して増加する変調された結果を有することにより、そのアクチュエータを変更する優先度が上がる。一方、エネルギー効率がより低いアクチュエータは、センサ値の差を相対的に減らす。エネルギー消費に対するそのような変調は、伝達行列事態に含まれてもよい。変調の値は、同一の割合だけ各アクチュエータを変更し且つエネルギー使用量を記録することにより決定されてもよい。種々の実施形態において、変調は、乗算、除算、加算、並びに加算演算、乗算演算及び更なる機能演算を含んでもよい座標系の更なる点を含んでもよい。いくつかの実施形態は、変調ではなく、変更するアクチュエータを決定するのに更なる要素としてエネルギー消費値を使用してもよい。例えば、最も少ないエネルギー消費を有しつつ適切なセンサ値の変化を提供するアクチュエータが選択されてもよい。
一実施形態において、加熱/冷却流のモジュールの効率は、エネルギー消費要素の一例である。効率は、モジュールのアクチュエータのある特定の入力設定(例えば、冷水バルブ設定)又は式2の百分率容量として解釈されてもよい。別の実施形態において、モジュールCindxの効率ηは、
Figure 2013543569
・・・式(5)
のように算出され、式中FanSpdは、モジュールが現在使用中の送風機最高速度の割合である。一実施形態において、送風機速度は、可変送風機を有さないモジュールに対して100%である。可変送風機を含むモジュールに対して、送風機は、最大設定を下回る設定、すなわち100%を下回って動作していてもよい。
一態様において、分母に送風機速度を使用する場合に送風機速度が100%となるため、その場合に可変送風機速度を有さないモジュールを停止することが優先されてもよく、それにより効率が落ちる。示されるように、冷却能力の量が所定の送風機速度に対して少ないため、あまり効率的でないモジュールはより低い効率を有する。他の効率が入力のレベルで除算されたあらゆる百分率出力を含んでもよく、それにより効率を測定する。例えば、特定の冷水バルブ設定に対する冷却出力は、バルブのアクチュエータに対する効率を算出するために使用されてもよい。
IV.センサ値を範囲に維持
図6は、本発明の一実施形態に従って高い効率で所望の範囲内にセンサ値を維持するように環境維持システムを制御する方法600を示すフローチャートである。方法600は、特定のセンサ値が所望の範囲外のセンサ値を有する場合に動作レベルを変化(上下)させるのに最適な1つ以上のアクチュエータを決定する。例えば、センサ(例えば、高温すぎるコールドアイルセンサ)の温度が閾値を上回る場合、方法600は、可能な限り迅速に又はエネルギー効率のよい方法でセンサを範囲に戻すように動作レベルを上昇(例えば、開始)させるべきモジュールを決定してもよい。一実施形態において、方法600の全て又は一部は、制御器15、あるいは本明細書において説明される別のコンピュータ又はプロセッサにより実行される。
ステップ610において、センサ値が所望の範囲外のセンサ(例えば、センサ3のうちの1つ)が識別される(例えば、温度が閾値を上回る)。一実施形態において、アプリケーション(例えば、システムのプロセッサ上で実行するソフトウェア)は、範囲外のセンサ値があるかを周期的にチェックしてもよい。部屋を冷却しているシステムの場合、センサ値が所望されるより高温なため、そのようなセンサはホットスポットと呼ばれてもよい。例えば高温すぎるセンサは、センサを監視し且つセンサが閾値より高くなる等、範囲外になる場合に警告信号を送出させることにより識別されてもよい。警告信号は、制御器15の内部又は制御システムの別の部分に生成されてもよい。
いくつかの実施形態において、所望の範囲は、部屋に対する目標値により規定されてもよい。範囲の例には、目標値からの+又は−のある特定の値、目標値を下回るあらゆる値、あるいは目標値を上回るあらゆる値が含まれる。所望の範囲は、デフォルト値(例えば、データセンタの温度に対して華氏83度)でウェブが設定可能であってもよい。別の実施形態において、例えば環境が特定の温度を上回るように要求される場合、センサは低温すぎるものとして識別されてもよい。そのような一実施形態において、モジュールは、冷却するのではなく加熱するだろう。
ステップ620において、動作レベルが上昇してもよい1つ以上のアクチュエータは、増加可能なアクチュエータとして識別される。動作レベルは、モジュールのアクチュエータに対してあらゆる動作パラメータを示してもよい。動作レベルが上昇しなくてもよい場合の例には、アクチュエータが最大動作レベルであるか、あるいはほぼ最大動作レベルである場合が含まれる。一実施形態において、この決定に対する基準は、動作レベルが最大レベルから所定の量(例えば、割合)内であるかであってもよい。いくつかの例において、特定のモジュールの動作パラメータのうちの一部のみが増加してもよく、他の例において、全ての動作パラメータが増加してもよい。いくつかの実施形態において、増加が停止されるモジュールに限定されてもよいため、増加は開始コマンドになる。
いくつかの実施形態において、アクチュエータを減少することにより、実際にはセンサを範囲内又は少なくとも範囲内により近接して移動させてもよい。そのような例は、方法500のステップ540に対して、例えば伝達行列要素が反対符号を有する場合に上述される。そのような実施形態において、ステップ620は、動作レベルの低下が識別されたセンサ値を有益に変化させてもよいアクチュエータを更に識別するように変更されてもよい。以下において、動作レベルの上昇を参照するが、動作レベルの低下が更に実行されてもよい。
ステップ630において、システムは、増加可能なアクチュエータ(すなわち、ステップ620で識別された)が範囲外のセンサ値を変化させる量(変化値)を予測する。この予測は、増加可能なアクチュエータ毎に実行されてもよい。予測された変化値は、アクチュエータの動作レベルが上昇する量に対して所定の値(例えば、50%又は100%)を使用することにより推定されてもよい。一実施形態において、予測された変化値は、それぞれ停止したモジュールを開始することによりセンサ温度が変化することが決定される程度を予測する。従って、一実施形態において、方法600は、各々が停止したCRACユニットを開始することによりホットスポットを消去することに対する衝撃を推定してもよい。
いくつかの実施形態において、アクチュエータ毎の所定の動作レベルの上昇は、予測された変化値を決定するために同一である。他の実施形態において、所定の上昇は、増加可能なアクチュエータ間で異なってもよい。例えば上昇は、動作レベルの上昇のある特定の割合(例えば、10%、30%又は100%)であってもよい。実際に実現される上昇の実際の量は、予測された変化値を決定するために使用された所定の上昇と異なってもよい。一実施形態において、所定の上昇は、アクチュエータ(モジュールと同等であってもよい)の動作レベルが上昇できる全ての量である。例えば、種々のモジュールが変化前に種々のレベルで動作していてもよいため、各モジュールは異なる増加を有してもよい。ステップ620は、所定の上昇及び/又は実際の上昇に対応し、識別されたアクチュエータが実際には適当な量だけ動作レベルを上昇できることを保証してもよい。
一実施形態において、予測された変化値ΔYは、例えば上述したように決定された伝達行列を使用して取得される。一実施形態は、
post−Ypre=TM(Sindx,UindxΔlevel(Uindx
・・・式(6)
を使用し、式中、Ypostは変化後の推定されたセンサ値であり、Ypreは範囲外の現在のセンサ値であり、Sindxは範囲外のセンサに対応し、Uindxは動作レベルを上昇するために考慮されているアクチュエータに対応する。所定の上昇は、伝達行列TMを作成するために使用された値とは異なってもよいΔlevelである。更に上述したように、Δlevelは、アクチュエータが増加する実際の量とは異なってもよい。
従って、いくつかの実施形態において、Ypostは、CRAC Uindxを開始した後の推定されたホットスポット温度であってもよく、Ypreはホットスポット温度であってもよく、Sindxはホットスポット温度を有するセンサに対応してもよく、Uindxは開始するために考慮されているCRACに対応してもよい。一実施形態において、Δlevel(Uindx)は、モジュールを開始する結果得られる推定された容量の変化である。
図6を再度参照すると、ステップ640において、アクチュエータは、予測された変化値に基づいて増加するために選択される。一実施形態において、Ypostの値は、予測された変化値を考慮する。別の実施形態において、Ypost−Ypreの値は予測された変化値である。TMを使用する他の予測された変化値が更に使用されてもよい。変化値は、センサ値を範囲内にすることを保証するために使用されてもよい。一態様において、Ypostは、特定の量により範囲を規定する最大/最小値より低く/高くなるように選択されてもよい。
例えば、Ypost又はYpost−Ypreは、センサ値を範囲内に変化させられるアクチュエータを決定するために使用されてもよい。変化値は、センサSindxに対して最も大きな値の変化を有するアクチュエータを更に決定してもよい。一実施形態において、最大の予測された変化値を含むモジュールは、おそらく最も少ない量の動作変化、すなわち最も少ないエネルギー量を含む範囲外状態を直すために使用される。別の実施形態において、最大の予測された変化値を含むモジュールは、センサ値を最速に変化させ、且つ範囲状態の別の出力がセンサSindxに対してより可能性が低いように値を最も自然に変化させると更に仮定されてもよい。例えば、最も低いYpostを提供するCRACユニットが使用されてもよい。変化値は、正の値又は負の値であってもよい。従って、最大という用語は、値が負である場合により小さな数を示してもよい。別の実施形態において、変化閾値より大きい予測された変化値を有するユニットのうちのいずれかが使用されてもよい。
上述したように、伝達行列は、例えば式4に対して上述したように、各アクチュエータに対してエネルギー消費要素を含んでもよい。これらのエネルギー消費要素は、例えば要素が既に伝達行列にない場合に変化値を変更するために更に使用されてもよい。例えばモジュールは、依然としてYpostを範囲内にしつつ、送風機消費電力又は他の電力の最小化に基づいて選択されてもよい。一実施形態において、これらのモジュールが最大容量を下回って動作するため、可変送風機速度を有するモジュールを開始するための優先度がある。
ステップ650において、選択されたアクチュエータの動作レベルは上昇するか、あるいは場合によっては上述したように低下する。上述したように、一実施形態において、アプリケーションは、ホットスポットを消去する(例えば、ホットスポットにおいて推定された最も低い温度を生成するCRACユニットを開始する)ことに対する最も大きな衝撃を有するアクチュエータの動作レベルを上昇させるものとする。実際の変化Δlevel(Uindx)を含む推定されたYpostが目標値を下回る場合、2つ以上のアクチュエータ又はモジュールは増加してもよい。一実施形態において、Ypostがどれも目標値を超えない(すなわち、−デッドバンド)場合、2つのアクチュエータは同時に増加してもよい。デッドバンドは、所望の範囲の境界を超える十分な変化が達成されることを保証できる。
上述したように、レベルが上昇している動作パラメータは、上述の百分率容量等の出力パラメータであってもよい。出力レベルを使用する問題は、特に出力レベルが他のモジュール及びアクチュエータの動作に依存してもよいと仮定すると、出力レベルが所定の入力に対するものを厳密に予測できないという問題であってもよい。従って、達成される動作レベルの厳密な変化を予測するのは困難である。いくつかの実施形態において、動作レベルの変化は既存の動作レベルから推定される。
いくつかの実施形態は、その低温の空気が増加後に再配布される方法を推定できる。均一な再配布の結果、動作モジュール時間の平均出力レベルのデフォルト値N/(N+1)が得られてもよく、式中Nは、実際の出力レベルに対して新しいモジュールを増加する前の動作ユニットの数である。平均は、各モジュールから取得されてもよく、あるいは全てのモジュールに対して同一であると仮定されてもよいため、1つの出力レベルのみが決定されてもよい。別の実施形態において、出力レベルは、上述(式3)のLOAD行列又はモジュールを増加し且つ変化した容量を測定することにより作成される同様のLOAD行列から決定されてもよい。実際の容量は実際の容量とは異なってもよいが、この推定された値は、2つ以上のモジュールを増加する必要があるかを決定するのを助長してもよい。
ステップ660において、他のアクチュエータは、例えばある特定の期間又は1又は複数の基準(例えば、方法500に対して上述されたような準定常ステート状態)に基づいて増加することを防止される。一実施形態において、選択された1つ以上のアクチュエータが増加した後にタイマを開始する。一態様において、タイマは、タイマの期限が切れるまで別の増加を防止するものとするウェブが設定可能なタイマ(デフォルト15分)である。別の実施形態において、モジュールは、アクチュエータを増加することに応答してセンサ値が平衡状態に達するまで増加することを防止されてもよい。一態様において、これにより、非常に多くのアクチュエータが同様の状態の時に増加し必要以上に多くのユニットを使用して非常に多くのエネルギーを使用することを防止する。
ステップ670において、伝達関数行列TMが更新される。例えば、増加したアクチュエータ(例えば、開始したモジュール)に対するTMの値は、式1を使用して更新されてもよい。このように、初期化がどれほど正確であったかを決定できるか、あるいは温度が維持されている負荷の変化に適応でき、TMはそれに応答して改善されてもよい。従って、このような更新は、伝達行列の変化する関係に対処できる。このように、システムは、システムの変数が変化するのに伴い変化するため、正確な推定を提供し続けることができる。
いくつかの実施形態において、更新の第1のステップは、いずれかの動作レベル及び選択されたアクチュエータの動作レベル(例えば、百分率容量)が上昇する前にセンサ値を記録することであってもよい(offにされるアクチュエータを識別することを含んでもよい)。更新の次のステップは、設定期間後又は基準を満たした後のセンサ値及び増加したアクチュエータの動作レベルを記録することであってもよい。一実施形態において、式1のTMの新しい行列要素は、変更されるアクチュエータに対応する列(j)毎に算出される。新しいTMは、
TM(i,j)new=gTM(i,j)new+(1−g)TM(i,j)old
・・・式(7)
のように更新され、式中gは、0〜1(例えば、デフォルト0.3)である。一実施形態において、gの値は、時間で変動してもよくあるいは別の変数に依存してもよい。他の実施形態において、他のより古いTM行列(すなわち、TMolder等の前の更新からの行列要素)は、組み合わせてTMnewを取得するために使用されてもよい。
一実施形態において、各アクチュエータの動作レベルが測定され、各列が更新される。従って、各行はセンサ値において同一の変化を有するが、動作レベルの変化は異なる。動作レベルがたとえ唯一のアクチュエータへの入力が変更されても変化してもよい出力レベルである場合、これらの他の行列要素が算出されてもよい。
他の実施形態において、低温すぎるアイルセンサ(又は他のセンサ)は、コールドスポットを識別するために更に使用されてもよい。上述の方法は、例えばシステムが上述したように周囲温度に対して環境を冷却するのではなく加熱しようとしている場合、加熱するHVACユニットをonにするために使用されてもよい。システムが冷却される実施形態において、コールドスポットを感知することにより、HVACユニットのうちの1つの停止をトリガしてもよい。低い温度が未熟なディスクドライブの故障と関連付けられることが示されているため、コールドスポットをそのように防止することにより助長できる。モジュールの停止は、他の理由でも発生する可能性がある。
方法600は繰り返されてもよい。例えばチェックは、いずれかのセンサが範囲外であるかを決定するために周期的に実行されてもよい。一実施形態において、増加量は少なくとも所定の量(100%未満であってもよい)であってもよい。増加が最大値未満である場合、次回では同一の1つ以上のアクチュエータが選択されてもよい。
いくつかの実施形態において、単に中間値だけ増減させることとは対照的に、アクチュエータ(場合によってはモジュール全体)を開始又は停止することが優先されてもよい。一実施形態において、アクチュエータをonにすることでセンサ値が範囲に戻らない場合、現在実行中の他のアクチュエータが解析されてもよい。別の実施形態において、on又はoffにされるだけでよいものの組合せがある場合にのみ優先され、他のアクチュエータは他の中間設定を有する。これらの実施形態は、本明細書において記載され且つ説明される他の実施形態に適用されてもよい。2つの異なる伝達行列が使用されてもよく、一方は開始又は停止すべきかを決定するためのものであり、他方は100%未満であってもよい変数増加又は変数減少すべきかを決定するためのものである。
ランダム性を使用した変更すべきパラメータの選択
最大変化値を提供する動作パラメータ(例えば、モジュールへの電力)が常に選択される場合、ある特定のパラメータは常に変更されてもよい。例えば変更される同一のアクチュエータは、多くの場合、その伝達行列がそのアクチュエータを間違って重み付けするようにしてもよいため、将来その同一のアクチュエータをより多く選択するようにしてもよい。偏りを除去するため、いくつかの実施形態はランダムに別のアクチュエータを選択する。選択は、依然としてセンサを範囲内にするアクチュエータに依然として限定されてもよい。ランダム性により、アクチュエータの増加の何らかの特定のパターンを防止でき、システムの自然関数ではなく偏りを反映する。
一実施形態において、選択により、第2の次善のアクチュエータをランダムに選択する。該次善のアクチュエータが選択される全体的な割合(例えば、時間の20%)が選択されてもよいが、該次善の選択に対する厳密な時間はランダムであってもよい。別の実施形態において、センサ値を範囲に戻すことが予測される各アクチュエータは、ランダムな選択中に選定されることを許可される。
エネルギー節約対信頼性
上述したように、アクチュエータを減少することにより、範囲外のセンサ値を所望に変化させることができる場合もある。変化は、センサ値を範囲に戻すようにするほど十分でない場合もあるが、依然として所望の方向である。例えばデータセンタを冷却する場合、モジュールを停止することにより、実際にはセンサ温度が低下する可能性がある。
いくつかの実施形態において、センサ値を可能な限り迅速に範囲に戻すよりもエネルギーを節約するのが望ましい。従って、アクチュエータを増加するよりも減少することに優先度があってもよい。一実施形態において、そのような優先度は、センサを範囲に戻す全てのアクチュエータの変更の間で行われてもよい。別の実施形態において、変化よりセンサが範囲に戻らない場合でも変化が所望の範囲に対するものである限り、減少することがそのように優先されてもよい。例えば、範囲が83°未満であり且つホットスポット温度が87°である場合、温度を85°に変化させる予測されたアクチュエータの変更は、アクチュエータを減少することを選択するのに許容可能だろう。一実施形態において、減少は常にアクチュエータを停止させている。
V.動作レベルを低下させることによるエネルギーの節約
全てのセンサが範囲内にある場合、アクチュエータの動作はエネルギーを節約するために軽減されてもよい。そのように軽減された動作は、センサ値に対して範囲外状態を防止するように制御されてもよい。例えば予測された変化値は、センサ値を範囲に戻すアクチュエータを決定するために先に使用されたように、範囲外状態を発生させないアクチュエータを決定するために使用されてもよい。
図7は、本発明の一実施形態に従ってアクチュエータの動作を軽減することにより環境維持システムを制御する方法700を示すフローチャートである。方法700は、減少するのに許容可能なアクチュエータを決定する。一実施形態において、方法700及び本明細書において示される他の方法の全て又は一部は、制御器15又は本明細書において説明される別のプロセッサにより実行される。
ステップ710において、現在エネルギーを使用中(すなわち、動作中)のアクチュエータが識別される。種々の実施形態において、エネルギー使用は、電気的(例えば、送風機への電力)、例えば加熱/冷却(例えば、冷水を冷却要素に供給する)等熱的又は可燃性(例えば、ガスを燃やして熱を提供する)であってもよい。一実施形態において、識別されたアクチュエータは、減少するのに使用可能なアクチュエータに更に限定される。例えば、モジュールへの電力が減少する動作パラメータのみである場合、モジュールは単一のアクチュエータとしてとらえられる。
ステップ720において、動作中のモジュールの動作レベルが少なくとも所定の量(例えば、50%又は100%)だけ低下した場合、システムは結果として得られるセンサ値を推定する。結果として得られるセンサ値は、動作中のモジュール毎に決定されてもよい。動作レベルは、モジュール自体であると考えられてもよいモジュールのアクチュエータに対応する。モジュールに対して2つ以上のアクチュエータがあってもよいため、結果として得られるセンサ値は、アクチュエータ毎の所定の変化量に対して決定されてもよい。一実施形態において、以下の式
indx:Ypost=Ypre+TM(Sindx,PindxΔlevel(Uindx
・・・式(8)
は、パラメータに対する所定の変更Δlevelに対してセンサ値の変化を決定するために使用される。
ステップ730において、どれも第1の範囲外にない推定されたセンサ値を有する現在実行中の複数のモジュールのうちの1つ以上の第1の集合が決定される。例えば、一実施形態において、モジュールをoffにすることにより範囲外状態を発生させない場合、モジュールは第1の集合の一部になる。別の実施形態において、モジュールの1つ以上のアクチュエータを減少することにより範囲外状態を発生させない場合、ハットモジュールは第1の集合の一部となる。本明細書において説明されるように、各センサは異なる範囲を有してもよく、第1の範囲という用語はこれらの種々の範囲を含む。範囲は、本明細書において説明されるような閾値及びデッドバンド値を含んでもよい。
ステップ740において、第1の集合の少なくとも1つのモジュールの動作レベルは、現在実行中の複数のモジュールのうちの1つ以上がどれも第1の範囲外にない推定されたセンサ値を有するとの決定に応答して低下する。一実施形態において、動作レベルは、モジュール全体に対する動作レベルであってもよく、例えばモジュールをオフにすることであってもよい。別の実施形態において、アクチュエータの2つ以上の動作レベルは、モジュールに対して低下してもよい。更に別の実施形態において、複数のモジュールの種々のアクチュエータが減少してもよい。一実施形態において、複数のアクチュエータを減少することの相乗効果は、線形性を仮定し且つ単に減少する結果得られる変化を合計することにより達成されてもよい。
いくつかの実施形態において、動作レベルが低下した可能性のあるモジュールの第1の集合は、閾値(例えば、40%)未満の効率(例えば、式5から算出されるような)を有するモジュールに更に限定される。一実施形態において、最も効率の低いアクチュエータが減少する。別の実施形態において、減少後のシステムに対する負荷が閾値を上回る場合、動作レベルは低下しない。一実現例において、負荷は、式3のLOAD行列を使用して算出されてもよい。
実行中の各モジュールのモジュールCindxの推定された百分率容量
Figure 2013543569
は、
Figure 2013543569
・・・式(9)
のように算出されてもよく、負荷行列は、他のアクチュエータの使用の容量に対する動作レベルを低下させる影響の基準を提供する。負荷行列要素LOAD(Cindx,Kindx)は、モジュールCindxの容量に対するモジュールKindxを停止する影響の基準を提供する。推定された百分率容量は、合計されるかあるいは組み合わされ、合計負荷を決定してもよい。%Cap(Sindx)は停止されるCRACユニットであり、%Cap(Cindx)は効率が決定されているCRACユニットに対するものである。
Figure 2013543569
に対するオーババーは、値が推定値であることを示す。
一実施形態において、別のモジュールの動作レベルを低下させることは、第1の集合の少なくとも1つのモジュールを減少した後に所定の時間の長さの期限が切れるまで防止される。別の実施形態において、伝達行列要素は、上述したのと同様の方法で更新されてもよい。例えば減少前後のセンサ値は、1つ以上のアクチュエータの動作レベルの変化と共に使用されてもよい。
VI.モジュールの停止
上述したように、動作パラメータの動作レベルを低下させるモジュールを識別する処理の結果、モジュールを完全に停止させてもよい。図8は、本発明の一実施形態に従って複数のモジュール及びセンサを含むモジュールを停止することにより環境維持システムを制御する方法800を示すフローチャートである。方法800は、停止させるのに最適なモジュールを決定できる。特に方法800は、システムが冷却でき且つセンサが温度センサである一例を説明する。
ステップ810において、効率が閾値未満であるモジュールが識別される。一実施形態において、閾値はウェブが設定可能である(例えば、デフォルトが40%)。効率は、式5又は本明細書において説明される他の方法を使用して測定されてもよい。一実施形態において、方法800は、第1の反復で効率の最も低い値を含むモジュールを最初に見つけることができる。このように、最も少ない量の冷却(実施形態に依存して加熱)を実行しているユニットは、場合によっては遮断するために目標にされてもよい。他の実施形態において、2つ以上のモジュールは反復毎に識別されてもよい。
ステップ820において、識別されたモジュールが遮断された場合に依然として存在し続ける各モジュールの百分率容量が推定される。一実施形態において、アプリケーションは、目標ユニットが式9を使用して遮断された場合に依然として存在し続ける各CRACユニットの百分率容量を推定するものとする。
ステップ830において、エネルギー負荷は百分率容量から決定される。一実施形態において、エネルギー負荷はシステムの平均百分率容量として示される。別の実施形態において、熱負荷は、例えばユニットの各容量で乗算された百分率容量の項を合計することにより決定されるようなシステムの合計エネルギー負荷である。別の実施形態において、エネルギー負荷はモジュールの最大百分率容量に対応する。
ステップ840において、推定されたエネルギー負荷が負荷閾値を上回るかが決定される。推定されたエネルギー負荷が負荷閾値を上回らない場合、方法は進む。推定されたエネルギー負荷が負荷閾値を上回る場合、ステップ810に戻り、場合によっては停止するために別のモジュールを識別する。この決定は、モジュールを停止することによりシステムに負荷をかけ過ぎるのを防止することを助長する。
一例において、目標モジュールを遮断することにより、ユニットが遮断された後も依然として存在し続ける容量の目標の割合(デフォルト=90%)をエネルギー負荷が超えない場合、そのモジュールは遮断される。一実施形態において、最も効率の低いモジュールを遮断することにより、依然として存在し続けるモジュールに負荷をかけ過ぎるかあるいはそのモジュールを範囲外状態(例えば、ホットスポット)にする場合、次に最も効率の低いモジュールは、停止するために評価されるものとする。
エネルギー負荷がモジュールの最大百分率容量である一実施形態において、効率閾値は、合計容量(例えば、95%)に対して使用された効率閾値より大きくてもよい。双方の種類の閾値が更に使用されてもよい。
ステップ850において、識別されたモジュールを停止した場合のセンサ温度が推定される。一実施形態において、目標CRACユニットが式8を使用して遮断された場合、アプリケーションはコールドアイル温度を推定する。
ステップ860において、推定されたセンサ値のうちのいずれかが閾値を上回るか、すなわち範囲外にあるかが決定される。このステップは、ホットスポットの作成を防止するのを助長する。推定されたセンサ値が閾値を上回らない場合、方法は進む。推定されたセンサ値が閾値を上回る場合、ステップ810に戻り、別のモジュールを識別する。
一例において、目標モジュールを遮断することによりいかなるセンサ値も閾値−デッドバンド(例えば、83F°−2F°)を上回らない場合、目標モジュールを停止してもよい。デッドバンドを使用することにより、センサが閾値に近接しないことを保証できる。
一実施形態において、ホットスポットを決定するために使用された目標温度(又は他のセンサ値)は、センサ温度が目標閾値−デッドバンドを上回るかを決定するために使用された同一の目標温度である。
ステップ870において、識別されたモジュールを停止する。一実施形態において、識別されたモジュールを停止した後、例えばウェブが設定可能な制限(デフォルト15分)を含むスタート/ストップタイマを開始し、且つスタート/ストップタイマの期限が切れるまで別のCRACユニットを停止しないものとする。
また、スタート/ストップタイマの期限が切れた後、目標CRACユニットに対応するTMの値は上述したように更新されてもよい。一態様において、更新の第1のステップは、スタート/ストップタイマが再開される前に停止されているCRACユニットを含むonにされるCRACユニットのコールドアイル温度及び百分率容量を記録するステップである。更新の第2のステップは、スタート/ストップタイマの期限が切れた後のコールドアイル温度及び百分率容量を記録するステップである。式1及び式2におけるTMの値は、例えば式5に従って算出され且つ更新される。
VII.開始及び停止
いくつかの実施形態において、制御システムは、変更(増減)すべきアクチュエータを追跡し、センサ値を所望の範囲内に効率的に維持し、且つ種々の時間において範囲外状態を発生させずにエネルギーを節約するために減少されてもよいアクチュエータを追跡できる。一実施形態において、センサ値を範囲内に維持するための論理は、エネルギー使用を減少するための論理より優先される。例えば、効率の低いCRACがあり且つ同時にホットスポットがある場合、CRACを停止するのではなく開始する。
図9は、発明の一実施形態に従ってセンサ値を所望の範囲内に維持し且つエネルギー使用を減少するようにモジュールの動作レベルを上下させることにより環境維持システムを制御する方法900を示すフローチャートである。方法900は、種々の状態に基づいて増減するのに最適なアクチュエータを決定できる。
ステップ910において、伝達行列TM(式1を参照)及び場合によってはLOAD行列(式3)が初期化される。一実施形態において、伝達行列TMは方法500を使用して初期化される。他の実施形態において、デフォルト値が仮定される場合又は別の制御システムにより初期化が実行される場合、伝達行列TMは、制御システムにより受信されてもよい。
ステップ920において、アクチュエータの動作レベルの変化が防止されるかが決定される。種々の実施形態に対して本明細書において示されるように、防止、例えば準定常ステート状態がセンサ値で達成されたかは、タイマ又は基準により決定されてもよい。一例として、タイマの期限が切れていない場合、「yes」分岐を取り且つモジュールを開始又は停止するのが許容可能になるまで待つか、あるいはモジュールのアクチュエータの動作レベルの他の変更を選択する。タイマの期限が切れた場合、「no」分岐に沿って継続する。変化の防止は、相対的に高い頻度で、例えば10秒毎又は30秒毎に周期的に決定されてもよい。
ステップ930において、センサ値が監視され、範囲外状態が存在するかを決定する。一実施形態において、温度(例えば、サーバに近接するコールドアイル温度)が監視され、ホットスポット(又は閾値の他の侵害)が発生しているかを決定する。範囲外状態が存在する場合、ステップ940に進む。範囲外状態が存在しない場合、ステップ960に進む。
一実施形態において、動作レベルの変化が防止されてもセンサ値は監視されてもよい。このように、防止が解除されると、防止が実施された際に発生していた範囲外状態に対処して前に進んでもよい。
ステップ940において、制御システムは、センサ値を範囲内にするために変更すべき動作パラメータを決定してもよい。そのような決定は、方法600の実施形態のうちのいずれかにより実行されてもよい。一実施形態において、動作パラメータは、モジュールを開始又は停止するかである。別の実施形態において、動作パラメータは、1つ以上のアクチュエータの動作レベルであってもよい。
ステップ950において、タイマを開始し、本明細書において示される他の方法の他のステップに対して更に実行されてもよいように、タイマの期限が切れるまで及び/又は準定常ステート状態が実現されるまで他の動作レベルが変化するのを防止する。タイマの期限が切れる際、TM行列及び/又はLOAD行列は更新されてもよい。
範囲外状態が存在しない場合に分岐に戻ると、ステップ960において、アクチュエータ(場合によってはモジュール全体を含む)のいずれかが効率閾値を下回る効率を有するかが決定される。効率は、本明細書において示される方法のうちのいずれかにより決定されてもよい。一実施形態において、モジュールがエネルギー閾値を下回って動作していない場合、ステップ920に戻る。別の実施形態において、モジュールが閾値を下回る効率で動作していない場合でも処理は継続してもよい。
ステップ970において、アクチュエータが範囲外状態を発生させずに動作レベルを低下(例えば、停止)させることができるかが決定される。範囲外状態が発生するかの予測は、本明細書において示されるように実行されてもよい。例えばそのような決定は、方法700の実施形態のうちのいずれかにより実行されてもよい。一実施形態において、依然として範囲外状態を発生させない最も高い能力を含むモジュールを減少(例えば、停止)する。
ステップ980において、タイマを開始し、タイマの期限が切れるまで及び/又は準定常ステート状態が達成されるまで他のユニットが変化するのを防止する。タイマの期限が切れる際、TM行列及び/又はLOAD行列は更新されてもよい。
VIII. 共分散行列
伝達行列に加え、一実施形態は共分散行列を使用できる。一実施形態において、共分散行列は、制御動作の数の大きさの正方行列である(すなわち、共分散行列は、uの要素数と同数の行及び列を有する)。初期の共分散行列は、大きな数(例えば、1e6)で乗算された恒等行列であってもよい。一態様において、共分散行列は、センサに対応する伝達行列の推定された値の不確実性を説明してもよい。例えば共分散行列は、特定のセンサに対応するパラメータの伝達行列ベクトルの不確実性を提供してもよい。一実施形態において、他方のアクチュエータを超える一方のアクチュエータの動作レベルを変化させることにより、そのパラメータがセンサに影響を及ぼす方法に関するより多くの情報があるため、そのアクチュエータに対応する共分散行列の要素を減少できる。
一実施形態において、伝達行列TMの初期化は以下の方法に従って実行されてもよい。以下のステップのうちのいくつかはオプションであってもよい。
ステップ1において、共分散行列COVは、対角エントリが10である対角行列に設定される。ステップ2において、伝達行列値TMはゼロに設定される。ステップ3において、全てのアクチュエータ(モジュール全体と同等であってもよい)は100%又は他の何らかの一般的な値になる。ステップ4において、「前の」温度が記録される。ステップ5において、最初のアクチュエータを下げるか又はoffにする。
ステップ6において、設定可能なタイマが時間切れになる(例えば、15分のデフォルトで)か、準定常ステート状態に到達するか、あるいはコールドアイル温度が限度を超える(例えば、87F°のデフォルトで)まで初期化は待つ。ステップ7において、「後の」温度が記録される。ステップ8において、温度変化dYが算出される。
ステップ9において、学習更新が実行される。一実施形態において、学習ベクトルLは、
L=COV*du/(λ+du*COV*du)
・・・式(11)
のように算出され、共分散行列は、
COV=(I−L*du)*COV/λ
・・・式(12)
を使用して更新され、式中、duは入力制御コマンドの変化(又は出力レベルの変化)のベクトルであり、Iは恒等行列であり、λは忘却係数である。一態様において、λは初期化中「1」であってもよい。行列は、例えば以下の対角要素を上記の対角要素に設定することにより対称になってもよい。
伝達行列は、予測誤差を算出してTMの行を更新することにより更に算出されてもよく、これは、全ての列又は単に変更されたアクチュエータに対応する列に対して実行されてもよい。一実施形態において、センサSindxに対する誤差はe=dY(Sindx)−TM(Sindx,:)duであり、新しいTM行列要素は、全ての列又は単に変更されたアクチュエータに対する列にわたりTM(sindx,:)=TM(sindx,:)+(Le)である。
ステップ10において、第1のアクチュエータは前のレベルに戻り、第2のアクチュエータはoffにされるかあるいは下げられる。ステップ11において、例えばステップ6のように時間を待つ。ステップ12において、「後の」温度が記録される。ステップ13において、温度変化dYが算出される。ステップ14において、例えばステップ9で上述されたように学習更新が実行される。ステップ15において、次のアクチュエータは、前のレベルに戻り、下げられるかあるいはoffにされる。
ステップ16において、後続ユニットが全て上/下又はon/offに循環されるまで、それらに対するステップ11に進んでもよい。第1のアクチュエータが最初に自身で下/offに循環されたように、最後のアクチュエータは自身で上/onに循環されてもよい。ステップ17において、最後の共分散行列の追跡はTrfinalとして格納されてもよい。いくつかの実施形態において、単に2つ以上のアクチュエータが一度に上に循環され、2つ以上のアクチュエータが一度に下に循環される。
システムが環境を制御している間、行列要素を更新するために同様の処理又は同一の処理が実現されてもよい。そのような処理は、例えば方法500〜900に対して上述されたように、あらゆる変更がアクチュエータに対して実行される際に実行されてもよい。
ステップ101において、変更すべきアクチュエータが決定される。ステップ102において、「前の」温度が記録される。ステップ103において、アクチュエータの動作レベルが変化する。ステップ104において、例えば上述の方法のステップ6のように時間を待つ。ステップ105において、「後の」温度が記録される。ステップ106において、温度変化(dY)が算出される。
ステップ107において、動作レベルの変化(du)が算出される。ステップ108において、例えばステップ9で上述されたように学習更新が実行され、λは1.0であってもよい。ステップ109において、予測誤差のベクトルは、E=dY−TM*duのように算出される。ステップ110において、平均絶対予測誤差|E|及び最大絶対予測誤差|E|\が算出される。ステップ111において、加重平均|E|及び|E|\である予測誤差の大きさは、
||E||=w*|E|+(1−w)*|E|
のように算出される。
ステップ112において、共分散行列の追跡に対する再設定値は、予測誤差の大きさを追跡の再設定値にリンクする式から決定される。図10はそのような式の一例を提供する。ステップ113において、追跡は、追跡が再設定値と同等であるように同一の乗数で既存の値の各々を変倍することにより更新されてもよい。一実施形態において、共分散行列は、その追跡が
COV=COVTrreset/trace(P)
のように追跡再設定値と同等であるように再設定される。ステップ114において、ステップ101を繰り返す。
IX.ペナルティ関数
上述したように、特定のアクチュエータに対する動作パラメータは、種々の動作値(レベル、設定等)を有してもよく、その結果、各々が種々のエネルギー量を消費する。このようにエネルギー(例えば、電気、蒸気、天然ガス等)の使用が異なることにより、システムを動作するコストは増加しうる。また、アクチュエータ毎のエネルギー使用は種々のアクチュエータ間で異なってよく、特定の種類のアクチュエータ内でも異なってよい。アクチュエータを起動するコストは維持コストに更に依存してよく、維持コストはアクチュエータの設定に依存して変動してよい。従って、アクチュエータを起動するコスト(例えば、エネルギー及び/又は維持)に対処することは、特に所望のセンサ範囲を更に所望する場合に非常に複雑になりうる。
いくつかの実施形態において、ペナルティ関数は、アクチュエータ状態U(動作レベル)への最適な変化dUを決定するために使用される。ペナルティ関数は、センサ読み取り値が範囲外である場合にペナルティを増加させるセンサ寄与、特定の設定で動作するためのエネルギー使用量コスト及び/又は維持コストに基づいて増加するコスト寄与及びアクチュエータの大きな変化を阻止する安定性寄与等の種々の寄与のうちの1つ以上を含むことができる。ペナルティ関数の最適化から算出されたような最適な変化dUは、アクチュエータを直接変更するために使用可能であるか、あるいは更なる関数において動作レベルUへの最終的な変化を提供するために使用可能である。動作パラメータの変化を決定するそのような方法は、本明細書で説明する他のあらゆる方法と組み合わせられる。
従って、ペナルティ関数Pは、範囲外のセンサ読み取り値に対するペナルティを含むことができる。範囲外のセンサ読み取り値に対するペナルティは、センサがより範囲から外れるにつれて増加しうる。ペナルティ関数は、センサが故障しているかあるいは信頼できないという状況標識を有するセンサからの読み取り値にペナルティを課さないように構成されてもよい。ペナルティ関数は、殆どが範囲外である設定可能な数のセンサ読み取り値をペナルティ関数に含まないように更に構成されてもよい。このように殆どが範囲外のセンサを含まないことは、故障しているが故障しているものとして診断されていないか、あるいは不適切に配置される可能性のあるセンサからの読み取り値にペナルティを課さないという目的にかないうる。
ペナルティ関数Pは、各動作コマンドの消費電力の推定値を更に含むことができる。例えば送風機の消費電力は、送風機の三乗に対して理論上の関係を有する。一般に、送風機速度と、流れと、気圧と、電力との関係を送風機の相似則と呼ぶ。送風機速度が動作コマンドであり、且つ最高速度の送風機の消費電力が既知である場合、送風機の相似則は、送風機速度が変更された後の送風機電力を推定するために使用可能である。例えば、動作レベルUが送風機速度に対する設定である場合、コストは(Uに比例しうる。各アクチュエータは、そのコスト寄与に対して異なる関数を有してよい。従って、動作レベルUは、コストc=C(U)を取得するための関数Cへの入力であってよい。
ペナルティ関数は、非エネルギーコストを有する制御動作に対するペナルティを更に有してよい。例えば、CRACユニットの開始又は停止は、それと関連付けられた維持コストを有してもよい。例えば、開始はペナルティがかかってもよい。また、送風機又は他のアクチュエータを起動することにより、ある特定の部分の寿命が縮まる可能性があり、そのようにその部分の寿命が縮まることは、アクチュエータに対する動作パラメータに依存してよい。従って、コストは全体的な動作のコストでありうる。
一実施形態において、動作レベルに対する最適化値は制限される。例えば最適化動作レベルは、全てがoffであるようなゼロになることを許可されなくてもよい。そのような制約は、ハード制約(例えば、そのような値は許可されない)であってもよい。別の実現例において、制約は、ペナルティ関数に含まれ且つより高い値を生じさせることにより、ある特定の点において最小値を提供しない(例えば、ゼロにおいて最小値ではない)。そのような制約をソフト制約と考えることができる。これらの他のペナルティ関数の項(一部又は全てがペナルティ関数の項と呼ばれてもよい)は、センサ値を予測するために使用される伝達行列TMから作成されてもよい。センサ値が基準値から離れている場合、項はより大きくなる。この項の増加は、非線形、例えば指数関数的であってよい。指数関数的な増加は、差分が十分に大きくなる(例えば、10度より大きい)まで発生しない可能性があるが、センサ値がより近接する(例えば、10度以内)時に線形であってもよい。
上述したように、ペナルティ関数は、システムの動作パラメータUの関数であってよい。種々の例において、動作パラメータは、システムの各アクチュエータへの入力、アクチュエータの一部に総合的に提供される入力又は全てのアクチュエータへの入力であってもよい。例えば、種々のモジュールの各送風機は、アクチュエータの一部(集合)に対応する入力となる同一の動作パラメータにより起動される。別の動作パラメータは、各モジュールの各冷水バルブを制御できる。あるいは、モジュール毎の送風機又はバルブは個別に制御可能である。
一実施形態において、最適な変化dUはペナルティ関数Pの最小化から決定される。ペナルティ関数の値pを推定するために、伝達行列TMは、結果として得られるセンサ値Sを推定するために使用可能である。このような推定は、特定の変化dUに対するペナルティ値pを決定するために使用可能である。従って、最適なdU(少なくともいくつかの停止基準に係る)は、現在のセンサ値の測定から決定可能である。
図11は、本発明の実施形態に従ってペナルティ関数を使用して環境維持システムを制御する方法1100を示すフローチャートである。方法1100は、図4の監視制御器15を用いて実現される。ステップ1110において、センサ値のフィードバックを受信できる。ステップ1120において、ペナルティ関数Pの初期値が決定される。
ステップ1130において、動作レベルの最適な変化dUは、センサ値S、アクチュエータの現在の動作状態U及び潜在的な変化dUの関数であるペナルティ関数Pを最適化(例えば、最小化)することで決定される。Pを最適化する例の更なる詳細を以下に提供する。一態様において、最適な変化dUは、単にペナルティ関数Pの値を減少するために最適である。ペナルティ関数Pに対して最も小さい値を有するグループのdUを決定する最適化法において、グループの異なるdUが試行されてもよい。
最適化に関して、Uがパラメータの動作レベル(値)を示すベクトルであると仮定すると、ペナルティ関数Pは単一の値Cを出力できる。従って、一実施形態において、ペナルティ関数は、所定の動作レベルの集合に対してP(U)として示される。本明細書で説明するように、伝達モデル(例えば、伝達行列)がUをセンサSに関連付ける(双方がPの決定で使用される場合)ために使用可能であるため、Uはシステムを規定できる。関数P(U)は、Cに対する最小値に対応するUを見つけるように最小化される。同等に、ペナルティ関数Pは、最大値(最小値の極値ではなく)を使用できる場合に作成される。極値は局所的極値であってよい。あらゆる適切な最適化法が使用されてもよい。例えば、極値を見つけるために、傾斜法又はより単純な誤差差異方法が使用可能である。一実施形態において、推計法(例えば、焼きなまし法及びモンテカルロ法)、決定論的方法、発見的方法又は他の適切な方法等の大域的最適化法が使用可能である。
ステップ1140において、更新された動作状態Uupdは、最適なdUoptから決定される。例えば、更新された動作状態Uupdは、U+dUoptとして解釈される。別の例として、Uupd=λU+(1−λ)(U+dUopt)であり、式中、λは混合係数である。ステップ11150において、システムはUの次の更新まで実行される。
上述したように、Pは種々の寄与を含むことができる。例えば、Pは、P+P+P+他の項に等しくてよく、式中、Pはエネルギー寄与であり、Pはセンサ値の差分に対するペナルティであり、Pは変更(例えば、開始及び停止)するための余分な維持コストに対するペナルティである。一実施形態において、維持コストは、動作レベルの少なくとも2つの範囲と関連付けられた種々のコストを有する。特定の種類のあらゆる寄与は、それぞれのセンサ値又はそれぞれのU値に対する部分からの集約(例えば、和、加重和又は他の関数)によって構成される。
上述したように、Pは、センサ値Sを基準値Rと比較することで決定可能である。各センサ値は1つ以上の基準値と比較可能である。例えば、特定のセンサに対するSはRと比較可能であり(例えば、センサ値が範囲外であるかを確認するために)、その差分は、他のセンサからの部分寄与と共に集約可能である(例えば、和)。集約の前に差分に比例因子(例えば、利得)を掛けて部分寄与を取得する。例えば、設定値との温度の差分が乗算され、その特定のセンサに対するペナルティの寄与を取得する。以下は、Pを算出する一例に対する同一の擬似コードであり(僅かに異なる表記が使用される)、その一般的な構造は他の寄与に対して使用可能である。
S_Penalty=0;
foreach sensor i
if((Snew(i)−R(i))>0)
S_Penalty=S_Penalty+(Snew(i)−R(i))S_Gain;//設定値にわたり全てのセンサを合計する。この場合、S_Gainは単に一定の乗数加重である
end if
end for
より一般的な意味において、特定のセンサに対するPへの部分寄与P(S)は、Rに基づく関数であってよい。例えば関数P(S)は、設定値を下回るあらゆる値に対するものであってよく、Rを上回るSの値に対してあらゆる所望の方法で増加してよい。多数の基準値と比較する代わりに、関数は、特定の範囲(例えば、所望の温度範囲)を下回る値に対して正であり、その範囲内でゼロであり、且つその範囲を上回る値に対して正であってよい。その範囲の両側の正数は、場合によっては種々の関数形式を有する。各寄与P(S)は、例えば合計することにより、Pに対する全体的な値を決定するために使用されてもよい。
従って、Sの測定値は、Pを決定するために使用可能である。一実施形態において、最適化法中にdUの種々の値毎にPを決定するために、伝達行列TMが使用されてもよい。例えば、Snew=S+TMdUであり、Snewは、duに対応する新しいPを決定するために使用されてもよい。
同様に、Pは、各アクチュエータの設定Uから決定可能である。Uは、動作レベルをエネルギーコスト(例えば、電気のコスト)に相関する関数に入力される。上述したように、関数は、Uの三乗の形式の項、例えばPE=aU +他の項を含むことができる。合計値Pは、部分PEiの集約から取得可能である。一実現例において、Pは、電気(又は燃料等の使用されているエネルギーは何でも)の現在の価格を考慮してよい。例えば、キロワット毎にPE=(0.05+0.95(U^3))priceである。一態様において、0.05は、アクチュエータをプラグインするための最小のエネルギー使用であると考えられる。尚、動作レベルとエネルギーコストとの関係は、あらゆる関数形式(例えば、指数関数的、二次式等)であってよい。
一実施形態において、最適化法中にdUの種々の値毎にPを決定するために、コンピュータシステムは新しいUを算出できる。例えば、Unew=U+dUである。Unewは、特定の変化dUに対応するPの値を決定するために使用可能である。Pは、Pと同様に決定可能である。P及びPの双方がUの状態に依存するため、状態Uに依存する寄与PはP+Pに等しくてよい。
ペナルティ関数への他の寄与は、変化dUの量に対するペナルティPdUを含むことができる。ペナルティは、Uの大きな変化にペナルティを課してよい。このように、方法はより安定される。一実施形態において、dUの値は、制限du_Limitと比較される。別の寄与は、下限又は上限を上回るUに対するペナルティであってよい。Uに対する制限の一例は、温度設定値を制御する(動作レベルの一例)場合にこの温度設定値が所定の下限を下回るべきではないか、あるいは冷却(及び凝縮の浄化)のために更なるエネルギーを発生させる可能性のある他の凝縮が生じうることである。従って、関数は下限を下回る正であり、上限と下限との間のゼロであり、且つその下限を上回る正であってよい。正の寄与は、単に利得項を掛けられた共役差積であってよく、あるいはより複雑な関数形式を有してよい。
場合によっては、Pへの寄与のうちのいずれかは、アクチュエータ又はセンサ毎の個々の項から構成されてもよい。例えばPE,Jは、エネルギー消費への動作レベルJのマッピングを提供してよい。マッピングは非線形であってもよい。マッピングは、非線形関数の線形係数を提供してもよい行列であってよい。一実現例において、Pにおけるマッピング行列(又はPに対する他の関数部分)は、いくつかのデータ点をサンプリングし且つ関数をパラメータ毎のパラメータ−エネルギー消費の関係に適合させる係数を決定することで更新可能である。
項は、エネルギー消費を減少するのを支援できる。しかし、これが長い間行われているため、それらをこれ以上下げさせない境界が発生しうる。1つの境界は、設定値を上回るSコストである。この境界領域にある間、Uに対する最適な設定があってよい(本明細書のモデルが正しいと仮定する)。他の寄与は、アクチュエータの境界極限又はあらゆる所定の時間でそれらが変化を許可される程度を示しうる。
図12は、本発明の実施形態に従って環境維持システムを制御する方法1200を示すフローチャートである。ステップ1210において、システムに関するデータが受信される。例えば、現在のセンサ値S(環境測定の現在の集合等、すなわち例えば、サーバ吸気温度)、現在の設定値R(測定毎の基準値の現在の集合等、すなわち例えば、コールドアイル設定値)、現在の動作値U(アクチュエータ毎の現在の動作状態等、すなわち例えば、送風機速度)及び現在のモデル(例えば、伝達行列TM)を受信できる。他のデータ、例えばアクチュエータ毎の受け入れ可能な動作変化の制限(アクチュエータ毎に異なってもよい)、アクチュエータ毎の動作の上限/下限(アクチュエータ毎に同一であってよい)、すなわち例えば、0%/100%の送風機速度、特定のペナルティ寄与を決定するために使用される、すなわち例えば、送風機速度制約、吸気温度制約に対するコスト乗数及び所定の動作値(例えば、PEに対して0.05値以上)において最も低いエネルギーコストを算出するために必要な各アクチュエータパラメータの集合を更に受信できる。
ステップ1220において、ペナルティ関数に対する初期値C0。dUがこの初期値に対して0であると解釈されるという点で、この初期値はベースライン値であると考えることができる。初期のペナルティ値は、C0=PenaltyFunction(0、U、S、R、TM)として書き込み可能である。上述したように、ペナルティへの1つの寄与は、各センサに対する基準値を上回るためのものであってよい。他の寄与は、それぞれのアクチュエータに対する動作の下限/上限を上回る動作値及び全てのアクチュエータに対するエネルギーの動作値コストのためのものであってよい。
ステップ1230において、最適化プログラムルーチンは、次のdUを見つけるように実行される。この次のdUは、ペナルティ関数から決定されたような最適なdUを見つけようとするアクチュエータへの可能な変化である。理想的には、最適なdUは、モデルTMを使用して現在のU、S及びRに基づいてペナルティ関数を最小化する。dUは、現在の状態Uからアクチュエータ毎の動作変化を規定してよい。一実施形態において、dUの値は、可変空間を検索してペナルティ関数を最小限にするUを見つけるベクトルとして理解可能である。例えば傾斜法(共役傾斜等)は、dUに対する最良な方向を決定するために使用可能である。一実現例において、適切な長さを決定するために、近似(例えば、二次近似)は、dUの方向のスカラに対して使用可能である。別の例において、dUは、Nelder−Mead法等のシンプレックス法を使用して決定される。
一実施形態において、dUの値は、大きすぎないように制約される。例えば、全てのアクチュエータに対する入力変化制限を上回るペナルティがあってよい。そのような制約は、変化を相対的に小さく維持することでシステムを安定させてよい。一実施形態において、アクチュエータ毎のあらゆる変化が動作変化制限を上回る場合、差分にコスト乗数を掛ける。アクチュエータに対するこの全てのペナルティの和は、コストにより算出可能であり、このペナルティ寄与を取得する。
ステップ1240において、現在のdUに対するペナルティが算出される。例えば、汎関数表記C=PenaltyFunction(du、U、S、R、TM)が使用可能である。上述したように、ペナルティへの1つの寄与は、基準値Rを上回るセンサ値に対するものであってよい。上述したように、動作値Uの変化に基づいて温度値(又は変化)を予測する伝達モデル(例えば、モデルが線形である場合の伝達行列)を取得できる。Snewを決定するために、TMは、動作変化dUを乗算して環境計測において投影変化(Snew)を見つけてよい。一実施形態において、基準値とは異なるセンサに対するペナルティ関数は、P(TM(U))として書き込まれてよく、TM(U)はセンサ値をベクトルSとして出力できる。ベクトルSは基準値Rと比較可能である。従って、Pは、R−TM(U+dU)の関数であってよい。例えば、いずれかの測定が基準値を上回る場合、コンピュータは差分にコスト乗数を掛ける。適用する全ての測定に対する和は、このペナルティ寄与により算出可能である。
アクチュエータに対する現在の投影動作値におけるエネルギー使用Pに対するペナルティの場合、消費されたエネルギーは、ある特定の種類(又はモデル)のアクチュエータに固有であってもよい関数に基づいて決定される。これらのアクチュエータの状態Uは、全体的なエネルギーコストを算出するために使用可能である(例えば、送風機速度に対する相似則等のアクチュエータ毎のk Wレーティングを使用して)。アクチュエータに固有の関数は、例えば動作レベル、電気のコスト、エネルギー効率パラメータ等のエネルギーコストを算出するために2つ以上のパラメータを有してもよい。新しいアクチュエータ状態Unewを決定するために、方法は、動作変化を現在のアクチュエータ状態に追加して投影動作値を取得できる。他の寄与と同様に、エネルギーコストはアクチュエータにわたり合計される。一例において、いずれかの新しい動作レベルがそれぞれの下限を上回る場合、方法は、差分にコスト乗数を掛けることができる。上述したように、上限は更に使用可能である。上述されていない他の寄与も同様に実現可能である。
ペナルティ関数は、最適なdUに対してハード制約又はソフト制約を有してよい。ハード制約は、dUがある特定の値を有することを許可されない場合の制約である。ソフト制約は、dUがある特定の値を有することを許可されるが、その値が絶えず増加し続けるペナルティ値を発生させることにより、ペナルティが高いためにそのようなdU値が調査される場合の制約である。
ステップ1250において、最適化法が収束されるかが決定される。例えば、ペナルティ値Cが1つの反復から次の反復に大幅に変化していないことが決定されてもよい。別の例として、2つの連続したdUは十分に類似する(すなわち、適応性のあるステップサイズが縮小するように算出される)ことが分かる。当業者は、使用可能な種々の収束基準を理解するだろう。収束基準を満たす場合、現在のdUは最適なdUとして解釈される。
ステップ1260において、収束基準を満たす場合、最適なdUは更に変更されてもよい。例えば最適なdUは、正規化、ランダム化(例えば、少量だけ)及び/又は変倍される。一実現例において、dU=rhodUであり、式中、rhoは倍率である。従って、実質的にUnew=U+rhodUである。正規化のために、正規化最適化プログラムは、最大許容動作変化制限に対してdUを算出してよい。ランダム化のために、ランダム化は可能なランダムスケジュール下で発生してよい。例えば、ペナルティ関数を最小限にしないがベースラインに対してペナルティ関数を減少するdUが使用可能である。一実現例において、そのような次善の決定は、確率pで各予測決定サイクルにおいて行われる。別の実現例において、次善の決定は、R予測決定サイクルの各々において行われる。変倍のために、デフォルトは、設定可能な定数だけ縮小可能である。一態様において、rhoは、例えばrhoとして使用可能な回帰テーブル及びピアソン相関係数を使用して伝達モデルの精度に基づいて決定されてもよい。
ステップ1270において、収束基準を満たさない場合、オプションのステップは、最大数の反復に達しているかをチェックするステップである。最大数に達していない場合、方法1200はステップ1230に戻り、新しいdUを選択する。最大数に達している場合、ある特定の基準に依存して現在のdUが使用されてもよい。一実施形態において、ランダム化を適用するために、最適化反復の数は、次善の決定が行われていなかった場合に取られる数を下回るようにランダム化される。
ステップ1280において、現在のペナルティ値CがベースラインC0を下回るかが決定可能である。C<C0の場合、現在のdUは全く変化がないよりは少なくともよい。その後処理は、ステップ1260に進み、あらゆる後処理を実行してもよい。CがC0を下回らない場合、ステップ1290でUは変化しなくてもよい。
図13は、本発明の実施形態に従って環境維持システムを制御するために使用される最適なdUを決定するシステム1300を示すブロック図である。システム1300は、HVACユニット、アクチュエータ、状況標識及びセンサに関するシステム400に類似する。例えば上述したように、更新学習モジュール1310は、状況及びセンサの更新を受信し、且つ伝達モデルTMを更新するためにそれらを使用できる。更新された伝達モデルは、最適なdUを決定するために最適化モジュール1320により使用可能である。
最適化モジュールは、更新学習モジュールを介して又は不図示の他の接続を通して上述したように入力(例えば、U及びS)を受信できる。最適化モジュール1320は、最適なdUを決定するために本明細書で説明したペナルティ関数のうちのいずれかを使用できる。上述したように、変調器1330は、最適なdUに対して後処理を実行するために使用されてもよい。これらのモジュールの出力は、アクチュエータに対して新しい設定Uを提供する。これらのモジュールのうちのいずれかは、1つの回路上で他のモジュールと組み合わせられた自身の回路であってもよく、あるいは汎用プロセッサで実現された異なるソフトウェア機能であってもよい。
一実現例において、現在の電気(又は他のエネルギー源)のコストは、電気供給者又はネットワークへのネットワーク接続を介して取得される。最適化モジュール1320は、動作レベルに基づいてコストを決定する関数において電気のコストを使用できる。
更新学習モジュール1310は、変調器1330を介して実現されたような後処理ステップにおいて使用可能な伝達モデルの精度を更に決定できる。以下の説明は、伝達モデルの精度を決定すること及び新しいUを決定するためにその精度を使用することに関する更に詳細で広範な概念を提供する。
一実施形態において、2つの反復ループが存在してよい。より高速な周波数ループは、測定されたセンサ値の差分に基づいて動作レベルを決定するフィードバックループであり、より低速な周波数ループは、最適化値を決定するフィードバックループである。例えば、最適化モジュール1320により決定されたような予測動作レベルは、第1の頻度で(例えば、15分おきに)且つそのような入力の間に入力され、センサからのフィードバックは、より単純なモデル(例えば、Federspiel 1及び図4において説明したようなTM及び/又は種々の比例積分偏差(PID)オブジェクトの逆を使用するモデル)に入力され、第2の頻度で(例えば、15秒又は30秒おきに)アクチュエータへの入力を決定できる。
X.伝達モデルの精度の使用
上述したように、実施形態は、予測誤差のサイズに基づいて個々の制御コマンド(すなわち、動作レベルを設定する)の変化dU(例えば、最適な変化)又は制御コマンドの累積サイズを変倍できる。米国特許出願第12/396,944号公報におけるフィードバック方法を含むdUの変化を決定する方法に関係なく、そのような実施形態を使用できる。例えばモジュール1330は、予測誤差が増加するにつれ予測制御動作があまり使用されなくなるように、最適な制御コマンドの累積サイズを変倍できる。一態様において、ピアソン相関係数は、伝達モデルにおいて予測誤差を決定するために使用される。
図14は、本発明の実施形態に従って伝達モデルの精度尺度を使用して環境維持システムを制御する方法1400を示すフローチャートである。ステップ1410において、センサ値Sが受信される(例えば、図4又は図13のセンサ3から)。センサ値Sは、本明細書で説明したあらゆる処理デバイスにより受信されてもよい。アクチュエータの第1の動作レベルU1に対応するセンサ値S。
ステップ1420において、アクチュエータの新しい動作レベルU2が受信される。新しい動作レベルU2及び第1の動作レベルU1は、変化dUだけ異なる。新しい動作レベルU2は、種々の方法で算出されてもよく、U1に追加される変化dUとして決定されてもよい。一例として、最適化法1100又は1200、あるいは米国特許出願第12/396,944号公報のフィードバック方法が使用されてもよい。
ステップ1430において、センサ値Sの変化dSが予測される。変化dSは変化dUに対応する。変化dSを単に新しいセンサ値として算出できることにより、変化dSを黙示的に算出する。別の例として、dSはdUから直接算出可能である。予測は、伝達モデルに従って実行される。S2−S1=TM(U2−U1)として書き込み可能なdS=TMdUを使用できる。
ステップ1440において、アクチュエータは、例えば一定の時間の長さの間、新しい動作レベルU2で動作される。新しい動作レベルU2が一般に第1の動作レベルU1とは異なるため、センサ値は変化すべきである。ステップ1450において、新しいセンサ値Snewが後続して測定される。尚、伝達モデルが完全に正確ではないため、新しいセンサ値は予測値と同一でなくてもよい。ステップ1410〜ステップ1450は繰り返され、複数の異なる時間において予測センサ値及び対応する測定値の多くの集合を取得できる。このデータは、今後使用するためにメモリに格納される。尚、このデータは、dS(測定差)及びdS’(予測差)として同等に格納される。そのようなデータ点のテーブルは回帰テーブルと呼べる。
ステップ1460において、伝達モデルの精度の計測値を決定できる(例えば、予測センサ値と測定センサ値との相関のレベルとして決定される)。計測値は、伝達モデルがステップ1410〜ステップ1450で取得されたデータ点にわたりどれほど正確であったかの尺度を提供できる。計測値は種々の方法で決定可能である。例えば、線形関係(例えば、S_predicted=αS_measured、式中、αは勾配である)及び線形関係に対するデータ点の誤差(S_predicted、S_measured)は、最小二乗を使用する等して測定される。従って、予測値と測定値との差分を集約して合計誤差値を取得する。この誤差の規模は、相関のレベルを決定するために使用可能である。誤差は正規化されて相関を得ることができる。一実施形態において、ピアソン相関係数が使用される。
上述したように、Sは、ベクトルであるために多次元データ点に対応する。従って、S_predicted−S_measuredも多次元データ点である。一実施形態において、差分の絶対値は、相関の決定の一環として合計される。別の実施形態において、センサ毎のデータ点(S_predicted、S_measured)は、仮定された関数関係(例えば、オフセットを有してもよい線形)と比較可能である。各センサは、異なる勾配及びオフセット等の異なる機能関係を有してよく、あるいはそれらは同一であってよい。
ステップ1470において、アクチュエータの動作レベルに対して提案された変更(例えば、最適なdU)が受信される。提案された変更は伝達モデルに基づく。例えば提案された変更は、本明細書で説明したように、ペナルティ関数を最小限にするために最適化法を使用して決定可能である。別の例として、伝達モデルは、センサを設定値内に維持する目的のみに基づいてアクチュエータが変更されるべき方法を決定するために、反転可能な伝達行列として使用可能である。
ステップ1480において、アクチュエータの動作レベルに対して提案された変更は、相関のレベルに基づいて調整される。例えば、最適なdUはピアソン係数により変倍可能である。このように、伝達モデルが正確な場合、最適なdUが直接使用される(Unew=U+dU)か、あるいは少なくともその殆どが使用される(例えば.8〜1.0の変倍係数rhoを含むUnew=U+rhodu)。しかし、伝達モデルが不正確になるように提供している場合、変倍係数rhoは小さい(例えば、0.1を下回る)か、あるいはゼロであってよい。
一実施形態において、調整は、ペナルティ関数においてパラメータを変更することで実現可能である。例えば、アクチュエータの動作レベルに対する下限及び上限は変更されてもよい。そのような変化は、この限界を超えるdUが引き起こすペナルティを減少するために提案されたdUを変更させる特定のペナルティ寄与に影響を及ぼす可能性がある。例えば、伝達モデルが正確になるように示されると仮定すると、モデルが正確であるためにセンサを設定値に準拠したままにしつつ、エネルギーを節約するように上限を減少できる。そのような制約により、モデルの誤差のためにセンサ値が設定値を上回る可能性があるため、モデルが正確でない場合、ペナルティ関数は動作レベルUを制限したくないだろう。
ステップ1490において、アクチュエータの動作レベルは、提案された変更への調整に基づいて設定される。このステップにおいて、Unewは、例えば式Unew=U+rhodUを介して先に提供されたように算出される。式中、rhodUは調整された変化である。他の例において、rhodUは、Uに追加する前に再度変倍、オフセット又は変更されてもよい。そのような変倍は、特定のアクチュエータ又はアクチュエータのグループに対する極値又は中点の設定に基づいて行われてもよい。これを以下に更に詳細に説明する。更に別の例において、rhodU及びUは、Unewを決定する関数に入力可能である。
図15は、本発明の実施形態に従って伝達モデルの精度を使用して環境維持システムを制御する方法1500を示すフローチャートである。ステップ1505において、伝達モデルTM、TMを更新する共分散関数(例えば、共分散行列)及び回帰テーブルは、校正試験中等に初期化される。回帰テーブルは、推定(又は予測)センサ値に対して1つの列及び実際の測定に対して別の列を有してよい。別の例において、測定における差分は1つの列にあってよく、予測差が格納される。テーブルの各行は、新しい更新事象に対応してよい。
ステップ1510において、最適なdUを算出する。一実施形態において、dUは、本明細書で説明したペナルティ関数のうちのいずれかのあらゆる最適化を使用して算出可能である。他の実施形態において、PIDオブジェクトを使用できる。一実施形態において、PIDオブジェクト11は、3つの独立パラメータ、すなわち比例値、整数値及び導関数値を算出する。比例値は、現在の誤差に対する反応を決定し、整数は、過去の誤差の合計に基づいて反応を決定し、導関数は、誤差が変化している速度に対する反応を決定する。これらの3つの動作の加重和は、制御バルブの位置、発熱体の電源又はモータドライブの速度等の最終制御要素又はデータ点を介して処理16を調整するか、あるいはHVACユニットの組込み(ローカル)デジタル制御装置の設定値を再設定するために使用される。更なる詳細を出願第12/396,944号公報において提供する。
ステップ1515において、予測精度計測値rho(例えば、相関係数)が算出される。一実施形態において、rhoは、回帰テーブルに格納された過去のデータセットのスライディングウィンドウを使用する予測vs測定センサデータのピアソン相関係数(rho)である。ウィンドウは、テーブルが大きくなりすぎない(例えば、最後の10,000個のデータ点のみが格納及び使用されてもよい)ように使用されてもよい。ユーザは、ウィンドウを規定できてもよい。一実現例において、rho<0の場合、rho=0と設定する。必要に応じて、相関を決定するのに十分なデータがない場合、ゼロ又はrhoの代わりに設定可能な係数を使用できる。
一実施形態において、ピアソン相関係数(又は他のあらゆる予測精度計測値)は、推定vs測定データの精度及び厳密度を推定できる。2次元グラフにおいて、一方の軸は推定され、他方の軸は測定される。グラフが完璧に線形である場合、係数は1となる。グラフが完全に分散する場合、係数は0又は負となる傾向がある。負の場合、値を0に設定する。
ステップ1520において、新しいduはrhoを使用して算出され、例えばdU=rhodUである。他の実施形態において、相関rhoは、dUを掛けられるスカラを決定するために使用可能である。あるいは、双方の種類の変倍を実行できる。スカラ乗数に加えて、他の種類の関数が更に考えられる。
ステップ1525において、Uが変更される前にセンサの値Sbeforeが決定される。ステップ1530において、新しいUが動作される(すなわち、アクチュエータの動作レベルが変更される)。新しいUは、Unew=Uold+dUとして算出可能である。ステップ1535において、アクチュエータは、システムが準平衡状態に達すること等の条件を満たすまで(例えば、タイマの期限が切れるまで)実行される。タイマは、5分、10分又は15分であってよい。ステップ1540において、新しいセンサ値Safterが記録される。
ステップ1545において、データ点は、メモリに格納される回帰テーブルに格納される。一実施形態において、dS=Safter−Sbeforeが算出され、dS’=TMdUが算出される。dS及びdS’のデータセットに対する値は、回帰テーブルに格納される。ステップ1550において、伝達モデルは、例えば上述したような共分散行列を使用して更新される。その後、方法はステップ1510から再度開始してよい。
XI.中点及び極値の変更
ペナルティ関数から取得された最適なdU及びdUを決定するために使用された伝達モデルの精度に関する情報の各々は、種々の方法で使用可能である。一実施形態において、最適なdUは、アクチュエータに対するフィードバック設定を乗算するスカラとして使用されてもよい。以下に説明するように、このフィードバック設定は、伝達モデルに基づいていなくてもよく、例えばセンサ分散の中点又は極値の値で動作する単純な発見的フィードバックループ(例えば、平均温度又はほぼ最高温度で動作するPIDループ)に基づいてよい。同様に、伝達モデルの精度は、伝達モデルを使用する中点又は極値の方法に対して使用可能である。
従って、センサ読み取り値からのフィードバックに基づく制御動作は、予測制御動作と組み合わせられてもよい。その2つが組み合わせられる場合、予測制御動作は、フィードバック制御動作を低減するオフセット及び/又は変倍係数であってよい。単一のフィードバックループは、極値制御に対して設定され(例えば、N番目に大きいセンサ読み取り値に応答する)、場合によっては第2のフィードバックループは、中点制御(例えば、平均センサ読み取り値)に対して設定される。予測制御動作は、制御動作(全てのアクチュエータに対して変更されないまま)とアクチュエータとの間のModブロックに影響を及ぼす。ペナルティ関数が殆どが範囲外であるN個の読み取り値を無視する場合、フィードバック制御動作は、これらのN個の範囲外の読み取り値にも応答すべきでない。
図16は、本発明の実施形態に従って中点及び極値制御を実現するシステム1600を示す図である。システム1600は、システム400及び1300と全体的に同様のアーキテクチャを有する。センサ値は、2つのフィードバック機構を介して極値及び中点比較モジュール1610及び予測決定モジュール1620に提供される。これらの2つのモジュールの各々を実現した実施形態に依存して種々の関数を実行できる。本明細書で説明したように、これらのモジュールのうちのいずれかは、1つの回路上で他のモジュールと組み合わせられた自身の回路であってもよく、あるいは汎用プロセッサで実現された異なるソフトウェア機能であってもよい。
一実施形態において、予測決定モジュール1620は、本明細書で説明したペナルティ関数のうちのいずれかを使用して最適化を実行する。例えば方法1100又は方法1200の実施形態により決定されたように、モジュール1620は、最適なdUの出力を提供できる。モジュール1620は、アクチュエータ又は設定をアクチュエータに送出するデバイスから現在のアクチュエータの状態を受信できる。他のフィードバックループは、センサ値を極値及び中点比較モジュール1610に送出する。極値及び中点比較モジュール1610は、センサ値の中点(例えば、平均値又は中央値)、極値(極大/極小)、並びに/あるいはN番目に大きい値又は小さい値(Nの値に依存して極値又は中点であると考えられる)を算出できる。
次に中点値は、基準中点値と比較可能である(例えば、割合又は差分)。比較の結果は、変化値(新しい設定であってよい)を決定するために使用可能である。例えば、単純な線形(又は他の関数)モデルは、中点及び基準値の差分を全てのアクチュエータの動作のレベルを提供する0%〜100%の出力値(又は同等に動作レベルの変化)に関連付けることができる。一実現例において、中点PID1615は、出力値を受信し、且つ前の出力値に基づいて出力値を変更できる。それにもかかわらず、モジュール1610及び中点PID(又は中点値を使用する他のあらゆる機構)は、動作レベル(変更値であってもよい)を中点変調器1617に提供できる。一態様において、各変調器は、センサ値が総合的に解析されている結果得ることのできる同一の値を受信できる。
一実施形態において、極値及び中点の制御の双方が使用される場合、予測決定モジュール1620は、上(及び/又は下)限の極値を上回ること、並びに所望のレベルからセンサ分散の中点が逸れることの双方に対してペナルティを使用できる。そのように実行する際、モジュール1620は、上(及び/又は下)限を上回る極値センサ読み取り値、並びに所望のレベルからのセンサ分散の中点の偏差を緩和しつつ、ペナルティ関数を最小限にする個々のコマンドを実現するように、(一般的に共通の)極値制御に対する個々のPIDコマンド及び(一般的に共通の)中点制御に対する個々のPIDコマンドを調整できる。
中点変調器1617は、全て又はほんの一部のセンサの中点を介して制御されるある特定のアクチュエータに出力を提供できる。中点変調器1617は、自身が制御するそれぞれのアクチュエータに対して異なる方法で一般に受信した設定を調整(例えば、オフセット又は変倍)できる。例えば中点変調器1は、制御出力をHVACユニット1のアクチュエータ2に提供する。変調器1617は、対応する変調器に固有のdU値に基づいて設定を調整できる。例えば、対応するアクチュエータに対するdU値が他のアクチュエータに対する他のdU値と比較して相対的に小さい場合、設定の変化は、中点PID1615から受信した値と比較して低減されてもよい。これは、中点PID1615の出力が全てのアクチュエータに対して同一であってもよいためであるが、モジュール1620は、特定のアクチュエータが高いセンサ値に対して小さな影響を与えるか、動作レベルを上昇するためにより高いコストを有するか、あるいは別のアクチュエータのdUを更に大きくできるようにしつつ中点を所望のレベルに近接させ続けることを認識する。この場合、この他のアクチュエータは、高いセンサ値に対してより大きな影響を与える。従って、特定のアクチュエータに対する中点PID1615からの変化値は、他のアクチュエータに対するものと同じほど高いべきではない。この場合、変化値は、より簡単にあるいはより安価に中点を設定値に低下(又は上昇)させる。
他のアクチュエータは、極値制御を介して同様に制御可能である。例えばモジュール1610は、それぞれの設定値からの差分から決定されてもよい極大及び/又は極小の(あるいはN番目に大きい/小さい)センサ値を識別できる。中点解析に対して上述したように、この極値センサ値は基準極値と比較可能であり、極値PID1625は動作設定を提供できる。最適なdUは、極値方法を介して制御されているそれぞれのアクチュエータに対して極値PID1625の出力を個々に調整するために使用可能である。
別の実施形態において、モジュール1610は、例えば本明細書及び米国非特許出願第12/396,944号公報において説明したように、伝達モデルに基づいてアクチュエータに対して提案された変更を決定できる。伝達モデルは、センサ値の中点又は極値を基準値内にするようにアクチュエータが伝達モデルに基づいて変更されるべき方法を決定するために使用可能である。中点又は極値が使用されているかに依存して、異なる伝達モデルが使用可能である。変化の量は、あらゆる特定のアクチュエータに対して異なってよく、あるいは全てのアクチュエータに対して同一であってよい。決定は、それぞれPID1615及びPID1625に提供される。しかし、本実施形態において、変調器1617及び1627は、伝達モデルの精度に基づいて設定を調整できる。
一実現例において、伝達モデルが予測値と測定値との比較に基づいて(例えば、ピアソン相関係数を使用して)正確であると決定される場合、モジュール1620(PIDが使用されなかった場合)の出力又はPIDからの出力は、殆ど変更されずに、例えば小さな混合係数で使用可能である。しかし、伝達モデルが正確でない場合、伝達モデルから提案された変更は殆ど又は全て無視される。
別の実現例において、変調器の変倍係数は、伝達モデルの精度に基づいて決定されてもよい。例えば変倍係数は、エネルギーコストを削減するために80%に設定可能である。伝達モデルが正確な場合、システムは、依然としてセンサ値を受け入れ可能な範囲内に維持できるべきである。従って、この動作最大値は、伝達モデルのより高いレベルの精度で低下しうる。しかし、伝達モデルの精度が低い場合、変倍係数はより高いべきであり、さもなければセンサ値は範囲外になる可能性が高い。同様に、全ての動作値の中点に対する変倍係数は、伝達モデルの精度に基づいて設定可能である。根本的な理由は、モデルの精度が良好な場合、方法が可能な限り新しい調整を維持できるか、さもなければシステムが、現在の値からそれぞれの最大値(又は最小値)にそれらを線形に再調整することである。
一実施形態において、可変速度送風機は、極値センサ値を変更し、且つ空気温度又は排出空気温度のセンサ設定値を戻して平均値、中央値又は他の統計値を調整するために使用される。別の実施形態において、アクチュエータ毎にPIDがあってよく、各PID出力の減衰は潜在的に異なる。従って、アクチュエータ毎の変調器。また、極値PIDは、各モジュール(例えば、各送風機)の第1のアクチュエータに対して使用可能であるが、第1のアクチュエータの集合毎の変調は別個に変調可能である。
冷却ユニットからの排出温度に影響を及ぼすアクチュエータに対する一例として、システムは、排出されている空気の温度に影響を及ぼしうる(例えば、排出温度設定値を調整することで)。調整は、ラック温度の平均(中点)に基づいて行われる。中点温度が高すぎる場合、システムはdU又は伝達モデルの精度に基づいて排出空気温度設定値を低下させる。基本的な一例において、全ての排出空気温度設定値は、中点モデルに基づいて同一となり、この共通の設定は、相対dU値に基づいて又は伝達モデルの同量の精度だけ(しかし、アクチュエータ毎の精度は更に隔離可能である)変更される。別の例として、多くの冷却ユニットの送風機の速度を調整することに関して、中点(又は極値)は、設定値を上回るN番目に高いラック温度として更に決定可能である(あるいは、誤差が設定値−フィードバックである場合、ラックセンサのN番目に小さい高温誤差はゼロを下回る)。
図17は、本発明の実施形態に従って動作の中点及び極値レベルを決定するために予測制御を使用する方法1700を示すフローチャートである。方法1700は、最適化法と予測精度方法との組合せを提供し、全ては極値/中点のフレームワーク内にある。
ステップ1705において、伝達モデルTM、TMを更新する共分散関数(例えば、共分散行列)及び回帰テーブルは、例えば関数テスト中に初期化可能である。このステップは、ステップ1505のように実行されてもよい。ステップ1710において、最適なdUは、ペナルティ関数の最適化等の伝達モデルTMを使用して算出される。ステップ1715において、例えば本明細書で説明したように、伝達モデルの予測精度の計測値(相関係数等)が決定される。測定値及び予測値が同一の平均値及び/又は同一の分散量を有する確率等の予測精度を評価する他の計測値は、相関係数の代わりに1715で使用可能である。ステップ1720において、dUはrhoにより変更されてもよい。このステップは、オプションのステップの一例である。
ステップ1725において、dUは、必要に応じてPID1615又はPID1625の出力等のフィードバック内側ループの出力と組み合わせられてもよい(例えば、結果がPID出力+dUとなるように、変調器を使用してPID出力を変倍することにより)。フィードバックループは、より多く、すなわち図17に示したループの反復の間に複数回繰り返すため、「内側」と呼ばれる。尚、方法1500には内側ループはなかった。他の実施形態において、dUは、アクチュエータの下/上限を調整するために(調整は、どちらの限界も上下させてもよい)、例えば動作レベルがどちらの限界も超えないように使用可能である。尚、下/上限(極値、すなわち最大値及び最小値とも呼ばれる)の変化は、アクチュエータの動作レベルの中点に対する全体的な値にも影響を及ぼしうる。同様に、中点レベルの変化は、アクチュエータの極値動作レベルに対する全体的な値にも影響を及ぼしうる。dUに対する値は、アクチュエータ毎に提案された動作レベルを調整するために使用可能なアクチュエータの集合に対する適切な中点を決定するために更に使用可能である。
ステップ1730において、アクチュエータ変調器の最大値及び最小値は、rhoを使用して調整可能である。変調器に対する最大値及び最小値は、特定のアクチュエータの動作の範囲を規定できるため、動作レベルに対して提案された値を変倍するために使用可能である。各アクチュエータ変調器は、動作に対する上限及び下限を有してよい。一実施形態において、これらの限界は変調器1627において規定される。上限は、伝達モデルが正確な場合はエネルギーを節約するように下げられ、且つ/あるいは伝達モデルが不正確な場合には上げられる。一実現例において、以下の式、すなわちMAXNew=MAXDefault+rho(MAX−MAXDefault)を使用できる。そのような式を使用することにより、上限は、MAXDefaultが>=MAXの場合には決して大きくならない。MAXDefault<=MAXの場合、新しい極値は決して小さくならない。どちらの構成も変調器の下限又は上限に対して使用されてもよい。しかし、伝達モデルが正確であるかあるいは不正確であるかに依存して、他の式により、変調器において上又は下の移動が可能になってもよい。アクチュエータの集合に対する中点制御のために使用された変調器に対する調整が同様に行われてよい。変調器に対するこれらの調整は、極値PIDオブジェクト及び中点PIDオブジェクトの双方に接続された変調器に適用できる。
ステップ1735において、センサに対する現在の値はSbeforeとして記録される。また、Uの新しい値は、提案された動作レベルU0、例えば最適化法の出力(すなわち、最適なdU)又は最適なdUを受信してもよいPIDの出力から算出される。一実施形態において、新しいU(Ubeforeと呼ばれる)は、極値内側PIDループにより制御される変調器に対するMAXを使用して算出される。一実現例において、Ubefore=MIN+PIDOutput(MAX−MIN)である。このUbeforeは、それぞれのアクチュエータに対する設定として使用可能である。極値PIDオブジェクトにより制御されている他のアクチュエータは、同様に算出されたUbefore値で更に制御可能である。従って、式は、ベクトル値を含むことができ、一般的にアクチュエータの集合全体に対する単一の式としてとらえられる。
ステップ1740において、アクチュエータは、タイマの期限が切れることにより測定されたような準定常状態条件等の条件に達するまで実行される。一実施形態において、Uの値は、この時間の間依然として更新可能である。例えば、センサ値の中点又は極値は基準値と比較可能であるか、あるいはセンサ値自体は基準値と比較可能であり、差分は、アクチュエータを増加又は減少するかを決定するために使用可能である。一実施形態において、この更新は、例えばFederspiel 1で説明したような伝達モデルの逆を使用する伝達モデルを使用できる。尚、新しいMAXは、内側ループのみが更新しているこのような期間中同一のままである。そのような二段階処理により、より精細なレベルでシステムを制御し、且つタイマの期限が切れるまで又は他の停止条件までシステムが一定のままでなければならなかった場合よりも適切な応答で負荷の変化に適応できる。
ステップ1745において、条件に達した後、センサの現在の値はSafterとして決定される。Uの現在の値(Uafterと呼ばれる)は、PID1615又はPID1625等の内側制御ループから出力された最後に提案された動作レベルU1から算出可能である。I番目のアクチュエータに対するUの値は、変化していないMAX値のベクトルのI番目の値を使用して算出可能である。一実現例において、Uafter_i=MAX_i+U1_i(MAX_i−MIN_i)である。
ステップ1750において、伝達モデルの更新及び伝達モデルの精度の決定に関連した種々の値が決定される。例えば、タイマの前後の動作レベルの変化を算出でき、すなわちduLearning=Uafter−Ubeforeである。時間の前後のセンサ値を算出でき、すなわちdS=Safter−Sbeforeである。伝達モデルがセンサ値の変化に対して予測したであろうものに対する値は、dS’=TMduLearningから算出可能である。dS及びdS’のデータセットは、回帰テーブル等に格納される。例えばdS及びduLearningを使用して本明細書で説明するように、伝達モデルTMは更に更新可能である。伝達モデルが更新されず、且つ現在のデータ点(dS、dS’)が例えば回帰テーブルに格納されていないことで相関を決定するために使用されないように、dUlearning=0の場合に例外を設けることができる。
その後、方法はステップ1710を繰り返すことができ、Uafterを最適なdUを決定するための開始点として使用できる。一実施形態において、前のMAX_I/MIN_I及びPID出力は、Uの大きな変化、あるいは上限又は下限の超過に対するペナルティ関数への寄与に対するパラメータを決定するために使用可能である。
環境維持システムのうちのいずれかは、HVACシステムに加えて多くの他の種類のシステムに適用可能である。例えば、制御システムを照明制御システムに適用できる。照明制御システムにおいて、アクチュエータは調光可能なバラストであってよく、センサは照明領域にわたり分散されたフォトセルであってよい。従って、制御されている環境のうちの1つは領域の照明であってよい。発見的制御は、共通の制御動作をPIDオブジェクトから全ての調光可能なバラストに送出することにより、中点(例えば、平均)照明レベル又は極値(例えば、最小)照明レベルを調整するためにPIDフィードバックループを使用できる。照明電力レベル、並びに超過及び中点設定からの偏差に対する潜在的なペナルティを含むペナルティ関数は、各バラストに対するPIDフィードバックコマンドを変調するためにモデル予測と共に使用される。これらの最適なPIDコマンドの予測変調は、最近の過去の予測精度が高くなかった場合に減衰され、予測精度が悪化するにつれ減衰量は増加する。
XII.コンピュータシステム
本明細書で説明したコンピュータシステムのうちのいずれかは、あらゆる適切な数のサブシステムを利用してもよい。そのようなサブシステムの例をコンピュータ装置1800の図18に示す。いくつかの実施形態において、コンピュータシステムは単一のコンピュータ装置を含み、サブシステムはコンピュータ装置の構成要素であってよい。他の実施形態において、コンピュータシステムは、各々が内部構成要素を含むサブシステムである多数のコンピュータ装置を含むことができる。
図18に示したサブシステムは、システムバス1875を介して相互接続される。更なるサブシステム、例えばプリンタ1874、キーボード1878、固定ディスク1879、ディスプレイアダプタ1882に結合されるモニタ1876及びその他が示される。入力/出力(I/O)制御器1871に結合する周辺デバイス及びI/Oデバイスは、シリアルポート1877等の当技術分野において既知であるあらゆる数の手段によりコンピュータシステムに接続される。例えば、シリアルポート1877又は外部インタフェース1881(例えば、イーサネット、Wi−Fi等)は、コンピュータシステム1800をインターネット等のワイドエリアネットワーク、マウス入力デバイス又はスキャナに接続するために使用可能である。システムバス1875を介した相互接続により、中央処理装置1873は、各サブシステムと通信し且つシステムメモリ1872又は固定ディスク11879からの命令の実行及びサブシステム間の情報のやり取りを制御できる。システムメモリ1872及び/又は固定ディスク1879は、コンピュータ可読媒体を組み込んでもよい。本明細書で説明した値のうちのいずれかは、1つの構成要素から別の構成要素及びユーザに出力される。
コンピュータシステムは、例えば外部インタフェース1881又は内部インタフェースにより共に接続された複数の同一の構成要素又はサブシステムを含むことができる。いくつかの実施形態において、コンピュータシステム、サブシステム又は装置は、ネットワークを介して通信できる。そのような例において、1つのコンピュータはクライアントであり、且つ別のコンピュータはサーバであると考えられる。この場合、各々は同一のコンピュータシステムの一部でありうる。クライアント及びサーバの各々は、多数のシステム、サブシステム又は構成要素を含むことができる。
本発明の実施形態のうちのいずれかは、モジュラーにおいて又は統合された方法でハードウェア(例えば、特定用途向け集積回路又はフィールドプログラマブルゲートアレイ)、並びに/あるいは一般にプログラム可能なプロセッサを含むコンピュータソフトウェアを使用して制御論理の形式で実現可能であることが理解されるべきである。本明細書において提供された開示内容及び教示に基づいて、当業者は、ハードウェア、並びにハードウェアとソフトウェアとの組合せを使用して本発明の実施形態を実現する他のやり方及び/又は方法を認識し且つ理解するだろう。
本出願において説明されたソフトウェアコンポーネント又はソフトウェア機能のうちのいずれかは、例えば従来の技術又はオブジェクト指向の技術等を使用するJava、C++又はPerl等の何らかの適切なコンピュータ言語を使用するプロセッサにより実行されるソフトウェアコードとして実現されてもよい。ソフトウェアコードは、格納及び/又は送信のために一連の命令又はコマンドとしてコンピュータ可読媒体に格納されてもよく、適切な媒体は、ランダムアクセスメモリ(RAM)、読み出し専用メモリ(ROM)、ハードドライブ又はフロッピディスク等の磁気媒体、あるいはコンパクトディスク(CD)又はDVD(デジタル汎用ディスク)等の光学記憶媒体及びフラッシュメモリ等を含む。コンピュータ可読媒体は、そのような記憶装置又は送信装置のあらゆる組合せであってもよい。
更にそのようなプログラムは、インターネットを含む種々のプロトコルに適合する有線ネットワーク、光ネットワーク及び/又は無線ネットワークを介して送信するように構成された搬送波信号を使用して符号化され且つ送信されてもよい。従って、本発明の一実施形態に係るコンピュータ可読媒体は、そのようなプログラムにより符号化されたデータ信号を使用して作成されてもよい。プログラムコードにより符号化されたコンピュータ可読媒体は、互換デバイスとパッケージ化されてもよく、あるいは他のデバイスとは別個に提供されてもよい(例えば、インターネットダウンロードを介して)。そのようなあらゆるコンピュータ可読媒体は、単一のコンピュータプログラム(例えば、ハードドライブ又はコンピュータシステム全体)上又は単一のコンピュータプログラム内に常駐してもよく、システム又はネットワーク内の種々のコンピュータプログラム上又は種々のコンピュータプログラム内に存在してもよい。コンピュータシステムは、本明細書で説明した結果のうちのいずれかをユーザに提供するモニタ、プリンタ又は他の適切なディスプレイを含んでもよい。
本明細書で説明した方法のうちのいずれかは、ステップを実行するように構成可能なプロセッサを含むコンピュータシステムを用いて全体的に又は部分的に実行されてもよい。従って、実施形態は、本明細書で説明した方法のうちのいずれかのステップを実行するように構成されたコンピュータシステムに関し、場合によっては種々の構成要素はそれぞれのステップ又はそれぞれのステップのグループを実行する。本明細書の方法のステップは、番号付きステップとして提示されるが、同時に又は異なる順序で実行されてよい。更にこれらのステップの一部は、他の方法の他のステップの一部と共に使用されてもよい。また、全て又は一部のステップはオプションであってもよい。また、方法のうちのいずれかのステップのうちのいずれかは、これらのステップを実行するモジュール、回路又は他の手段により実行される。
特定の実施形態の特定の詳細は、本発明の実施形態の趣旨及び範囲から逸脱せずに何らかの適切な方法で組み合わせられてもよい。しかし、本発明の他の実施形態は、各々の個々の態様又はこれらの個々の態様の特定の組み合わせに関連する特定の実施形態に関してもよい。
本発明の例示的な実施形態の上記の説明は、例示及び説明のために提示された。本発明の例示的な実施形態の上記の説明は、本発明を網羅すること又は説明された厳密な形式に限定することを意図せず、多くの変更及び変形は上記教示に鑑みて可能である。本発明の原理及びその実際的な適用範囲を最適に説明することにより、他の当業者が考えられる特定の用途に適するような種々の実施形態における本発明及び種々の変更を含む本発明を最適に利用できるように、実施形態は選択され且つ説明された。
特に指示のない限り、単数形の記述は「1つ以上」を意味することを意図する。
上述の全ての特許、特許出願、公開及び明細書は、全ての目的のために全ての内容がこの引用により本明細書に組み込まれる。従来技術であると認められるものはない。

Claims (30)

  1. 複数のアクチュエータ及び各々が環境の物理的状態の値を測定する複数のセンサを含む環境維持システムを制御する方法であって、
    前記センサに対応するセンサ値Sを受信することと、
    前記アクチュエータの動作レベルUへの潜在的な変化dUを決定することと、
    予測された前記センサ値Sの変化dSを決定するために、アクチュエータの動作レベルの変化と結果として得られるセンサ値の変化との関係を提供する伝達モデルをdUに適用し、
    前記予測された変化dSから新しいセンサ値SNewを決定し、
    1つ以上のそれぞれの設定値に対する新しいセンサ値の各々の関係に基づくペナルティ関数への第1の寄与を決定するためにSNewを使用することにより、コンピュータシステムを用いて前記ペナルティ関数への前記第1の寄与を算出することと、
    前記動作レベルUへの前記潜在的な変化dUに対応する新しい動作レベルUnewを決定し、
    新しい動作レベル毎に前記新しい動作レベルで前記対応するアクチュエータを動作するコストを決定し、
    前記コストを集約して第2の寄与を取得することにより、前記コンピュータシステムを用いて前記ペナルティ関数への前記第2の寄与を算出することと、
    前記潜在的な変化dUに対する前記第1の寄与及び前記第2の寄与に基づいて前記ペナルティ関数の第1の値を決定することと、
    前記動作レベルUへの複数の潜在的な変化に対する前記ペナルティ関数の値に基づいて前記アクチュエータの前記動作レベルUへの最適な変化を決定することと、
    前記最適な変化に基づいて前記アクチュエータの前記動作レベルを設定することと、
    を備えることを特徴とする方法。
  2. 前記最適な変化に基づいて前記アクチュエータの前記動作レベルを設定することは、
    前記最適な変化を調整することと、
    前記調整された最適な変化をUに追加することとを含むことを特徴とする請求項1記載の方法。
  3. 前記最適な変化を調整することは、正規化、ランダム化及び変倍のうちの少なくとも1つを含むことを特徴とする請求項2記載の方法。
  4. 予測センサ値を測定センサ値と比較することで前記伝達モデルの精度を決定することを更に備え、前記最適な変化は、前記伝達モデルの精度に基づいて調整されることを特徴とする請求項2又は3記載の方法。
  5. 前記最適な変化に基づいて前記アクチュエータの前記動作レベルを設定することは、
    前記アクチュエータに対する新しい設定を受信することと、
    前記最適な変化に基づいて前記新しい設定を調整することとを含むことを特徴とする請求項1記載の方法。
  6. 前記センサ値の中点と中点基準との比較又は
    前記センサ値の極値と極値基準との比較に基づいて新しい設定を提供するPIDオブジェクトを使用して前記新しい設定を算出することを更に備えることを特徴とする請求項5記載の方法。
  7. 最大数の前記動作レベルUへの潜在的な変化は、前記最適な変化を決定するように解析され、前記対応する潜在的な変化のうちの少なくとも1つがUへの変化を含まない前記ペナルティ関数を下回る前記ペナルティ関数に対する値を有する場合、前記最適な変化は前記複数の潜在的な変化から選択されることを特徴とする請求項1から6のいずれか1項に記載の方法。
  8. 新しい動作レベルでアクチュエータを動作する前記コストは、エネルギーコスト及び維持コストのうちの少なくとも1つを含むことを特徴とする請求項1から7のいずれか1項に記載の方法。
  9. 前記アクチュエータを動作する前記コストは、
    エネルギー使用率毎のコストを受信することと、
    エネルギー使用率毎の前記コスト及び前記新しい動作レベルを前記対応するアクチュエータを動作する前記コストを出力する関数に入力することとを含むことを特徴とする請求項8記載の方法。
  10. 前記第1の寄与を決定するためにSNewを使用することは、
    センサ毎に、
    Newをパラメータとしてセンサの少なくとも1つの設定値を有するそれぞれの関数に入力し、部分的な第1の寄与の出力を取得することと、
    前記センサから前記部分的な第1の寄与を合計することとを含むことを特徴とする請求項1から9のいずれか1項に記載の方法。
  11. それぞれの関数は、
    前記対応する新しいセンサ値SNewと設定値との差分を利用することと、
    前記差分に利得係数を掛けて前記センサに対する前記部分的な寄与を取得することとを含むことを特徴とする請求項10記載の方法。
  12. 前記アクチュエータの動作レベルUへの潜在的な変化dUを決定することは、前記新しい動作レベルUnewを受信することで達成され、この場合Unew=U+dUであり、予測されたセンサ値Sの変化dSを決定するために伝達モデルをdUに適用することは、前記伝達モデルを前記潜在的なUに対する新しい値に適用することで達成されることを特徴とする請求項1から11のいずれか1項に記載の方法。
  13. 前記アクチュエータの動作レベルUへの潜在的な変化dUを決定することは、前記現在の潜在的な変化を決定するために過去の潜在的な変化を取り入れる最適化アルゴリズムを使用することを特徴とする請求項1から12のいずれか1項に記載の方法。
  14. 前記設定値は種々のセンサ値を含むことを特徴とする請求項1から13のいずれか1項に記載の方法。
  15. 前記伝達モデルは数値の行列であることを特徴とする請求項1から14のいずれか1項に記載の方法。
  16. 実行される時に複数のアクチュエータ及び各々が環境の物理的状態の値を測定する複数のセンサを含む環境維持システムを制御するようにコンピュータシステムを制御する請求項1から15のいずれか1項に記載のステップを含む複数の命令を格納する有形のコンピュータ可読媒体を備えることを特徴とするコンピュータプログラム。
  17. 複数のアクチュエータと、
    各々が環境の物理的状態の値を測定する複数のセンサと、
    請求項1から15のいずれか1項に記載のステップを実行するように構成された少なくとも1つのプロセッサと、
    を備えることを特徴とする環境維持システム。
  18. 複数のアクチュエータ及び各々が環境の物理的状態の値を測定する複数のセンサを含む環境維持システムを制御する方法であって、
    前記アクチュエータの動作レベルUに対応する前記センサに対応するセンサ値Sを受信することと、
    前記アクチュエータの新しい動作レベルUを受信することであり、この場合U−U=dUであることと、
    伝達モデルに従って実行される前記変化dUに対応する前記センサ値Sの変化dSを予測することと、
    前記新しい動作レベルUで前記アクチュエータを動作することと、
    後続して新しいセンサ値Snewを測定することと、
    前記受信、予測、動作及び後続測定を複数回繰り返すことと、
    前記予測センサ値を前記測定センサ値と比較することで前記伝達モデルの精度の計測値を決定することと、
    前記伝達モデルに基づく提案された前記アクチュエータの前記動作レベルへの変更を受信することと、
    前記精度計測値に基づいて、前記提案された前記アクチュエータの前記動作レベルへの変更を調整することと、
    前記提案された変更に対する前記調整に基づいて前記アクチュエータの前記動作レベルを設定することと、
    を備えることを特徴とする方法。
  19. 前記精度計測値を決定することは、
    前記予測センサ値及び前記測定センサ値を格納することと、
    センサ毎に、
    前記予測センサ値と前記測定センサ値との差分を算出することと、
    前記差分を集約することと、
    前記集約された差分の規模に基づいて前記計測値を決定することとを含むことを特徴とする請求項18記載の方法。
  20. 前記予測センサ値は予測センサ値の差分として格納され、前記測定センサ値は測定センサ値の差分として格納されることを特徴とする請求項19記載の方法。
  21. 予測センサ値と測定センサ値との差分を算出することは、前記差分を算出する前に前記値のうちの少なくとも1つを変倍及び/又はオフセットすることを含むことを特徴とする請求項19又は20記載の方法。
  22. 前記精度計測値は、ピアソン相関係数を使用して決定されることを特徴とする請求項19から21のいずれか1項に記載の方法。
  23. 前記提案された変更は、
    予測センサ値と設定値との差分を含む第1の寄与と、
    前記提案された変更の結果得られる動作レベルと関連付けられたコストを含む第2の寄与とを含むペナルティ関数の最適化から決定されることを特徴とする請求項18から22のいずれか1項に記載の方法。
  24. 新しいセンサ値Snewを測定することは、前記新しい動作レベルUで前記アクチュエータを動作し始めた後の所定の時間において実行されることを特徴とする請求項18から23のいずれか1項に記載の方法。
  25. 前記センサ値Sの変化dSを予測することは、dUに前記伝達モデルを掛けることで達成されることを特徴とする請求項18から24のいずれか1項に記載の方法。
  26. 前記アクチュエータの新しい動作レベルUを受信することは、dU及びUを受信することで達成されることを特徴とする請求項18から25のいずれか1項に記載の方法。
  27. 前記アクチュエータの前記動作レベルへの前記提案された変更を調整することは、
    前記アクチュエータの動作の範囲を指定することと、
    前記範囲に基づいて前記提案された変更を調整することとを含むことを特徴とする請求項18から26のいずれか1項に記載の方法。
  28. 前記範囲は前記動作を完全な動作性より狭い範囲に限定し、前記提案された変更は動作の百分率であり、前記調整することは前記限定された範囲に基づいて前記百分率を変更することを特徴とする請求項27記載の方法。
  29. 前記伝達モデルは数値の行列であることを特徴とする請求項18から28のいずれか1項に記載の方法。
  30. 前記提案された変更は、前記伝達行列を反転し、且つ前記反転された伝達行列にセンサ値と設定値との差分を掛けることで決定されることを特徴とする請求項29記載の方法。
JP2013525010A 2010-08-20 2011-08-22 Hvacシステムに対するエネルギー最適制御決定 Pending JP2013543569A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US37577810P 2010-08-20 2010-08-20
US61/375,778 2010-08-20
PCT/US2011/048677 WO2012024692A2 (en) 2010-08-20 2011-08-22 Energy-optimal control decisions for hvac systems

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2017002098A Division JP6514242B2 (ja) 2010-08-20 2017-01-10 環境維持システムの制御方法、環境維持システム及びプログラム

Publications (1)

Publication Number Publication Date
JP2013543569A true JP2013543569A (ja) 2013-12-05

Family

ID=45605713

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2013525010A Pending JP2013543569A (ja) 2010-08-20 2011-08-22 Hvacシステムに対するエネルギー最適制御決定
JP2017002098A Active JP6514242B2 (ja) 2010-08-20 2017-01-10 環境維持システムの制御方法、環境維持システム及びプログラム

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2017002098A Active JP6514242B2 (ja) 2010-08-20 2017-01-10 環境維持システムの制御方法、環境維持システム及びプログラム

Country Status (6)

Country Link
US (2) US8924026B2 (ja)
EP (1) EP2606406B1 (ja)
JP (2) JP2013543569A (ja)
KR (1) KR20130098346A (ja)
CA (1) CA2808503A1 (ja)
WO (1) WO2012024692A2 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017506734A (ja) * 2014-03-28 2017-03-09 三菱電機株式会社 蒸気圧縮システムを制御するための極値探索コントローラー及び方法
JP2017110904A (ja) * 2010-08-20 2017-06-22 ヴィジレント コーポレイションVigilent Corporation 環境維持システムの制御
WO2020121551A1 (ja) * 2018-12-12 2020-06-18 日本電信電話株式会社 マルチデバイス連携制御装置、マルチデバイス連携制御方法およびマルチデバイス連携制御プログラム、並びに、学習装置、学習方法および学習プログラム
JPWO2021059412A1 (ja) * 2019-09-25 2021-04-01
CN115978722A (zh) * 2023-03-17 2023-04-18 四川港通医疗设备集团股份有限公司 一种基于人工智能的净化机组监测管理方法及系统
WO2023176151A1 (ja) * 2022-03-17 2023-09-21 トリニティ工業株式会社 空調制御装置及びそのためのプログラム

Families Citing this family (89)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009117741A1 (en) 2008-03-21 2009-09-24 The Trustees Of Columbia University In The City Of New York Decision support control centers
US9709965B2 (en) 2008-12-04 2017-07-18 Baselayer Technology, Llc Data center intelligent control and optimization
WO2010096783A1 (en) 2009-02-20 2010-08-26 The Trustees Of Columbia University In The City Of New York Dynamic contingency avoidance and mitigation system
EP2916622B1 (en) 2009-10-08 2019-09-11 Delos Living, LLC Led lighting system
FR2964204B1 (fr) * 2010-08-25 2012-08-17 Schneider Electric Ind Sas Procede de determination de parametres de regulation d'un systeme hvac
US8997847B2 (en) * 2010-09-10 2015-04-07 Ford Global Technologies, Llc Cooling in a liquid-to-air heat exchanger
JP6313217B2 (ja) * 2011-12-12 2018-04-18 ヴィジレント コーポレイションVigilent Corporation Hvacユニットの気温制御
US9317026B2 (en) * 2013-05-31 2016-04-19 Patrick Andrew Shiel Method for determining the unique natural thermal lag (NTL) of a building
WO2013130956A1 (en) 2012-03-02 2013-09-06 Vigilent Corporation Multi-dimensional optimization for controlling environmental maintenance modules
US9582009B2 (en) * 2012-04-30 2017-02-28 SmrtEn, LLC System and method for optimizing and reducing the energy usage of an automatically controlled HVAC system
WO2014022593A1 (en) 2012-08-03 2014-02-06 Synapsense Corporation Apparatus and method for controlling computer room air conditioning units (cracs) in data centers
JP2015534701A (ja) 2012-08-28 2015-12-03 デロス リビング エルエルシーDelos Living Llc 居住環境に関連するウェルネスを増進するためのシステム、方法、及び物品
US9898917B2 (en) * 2012-10-15 2018-02-20 Vigilent Corporation Method and apparatus for providing environmental management using smart alarms
WO2014075108A2 (en) * 2012-11-09 2014-05-15 The Trustees Of Columbia University In The City Of New York Forecasting system using machine learning and ensemble methods
WO2014075280A1 (en) * 2012-11-16 2014-05-22 Honeywell International Inc. Fuse multiple drawings into an equipment (bim) model
EP2746888B1 (en) * 2012-12-20 2019-02-20 ABB Research Ltd. Method and system for fluid flow control in a fluid network system
CN109582109A (zh) * 2013-03-07 2019-04-05 贝斯莱尔科技有限公司 数据中心控制方法和系统
US9810442B2 (en) 2013-03-15 2017-11-07 Google Inc. Controlling an HVAC system in association with a demand-response event with an intelligent network-connected thermostat
US9595070B2 (en) 2013-03-15 2017-03-14 Google Inc. Systems, apparatus and methods for managing demand-response programs and events
US9807099B2 (en) 2013-03-15 2017-10-31 Google Inc. Utility portals for managing demand-response events
EP3961342A1 (en) * 2013-04-19 2022-03-02 Google LLC Controlling an hvac system in association with a demand-response event
JP6247746B2 (ja) * 2013-05-08 2017-12-13 ヴィジレント コーポレイションVigilent Corporation 環境に管理されるシステムにおける影響の学習
US9070274B2 (en) * 2013-05-09 2015-06-30 International Business Machines Corporation Method for early detection of cooling-loss events
JP6179196B2 (ja) * 2013-05-31 2017-08-16 富士通株式会社 データセンター
US9395712B2 (en) * 2013-05-31 2016-07-19 Patrick Andrew Shiel Building energy usage reduction by automation of optimized plant operation times and sub-hourly building energy forecasting to determine plant faults
US9506666B2 (en) * 2013-06-13 2016-11-29 Trane International Inc. System and method for monitoring HVAC system operation
US9416987B2 (en) 2013-07-26 2016-08-16 Honeywell International Inc. HVAC controller having economy and comfort operating modes
US9618224B2 (en) 2013-07-26 2017-04-11 Honeywell International Inc. Air quality based ventilation control for HVAC systems
EP3033239A1 (en) * 2013-08-18 2016-06-22 Sensibo Ltd. Power consumption assesment of an hvac system
US10114721B2 (en) 2013-08-18 2018-10-30 Sensibo Ltd. Power consumption assesment of an HVAC system
CN106105169B (zh) * 2013-12-17 2020-03-31 贝利莫控股公司 移动通信装置和用于管理多个致动器的操作的方法
US9883009B2 (en) * 2013-12-27 2018-01-30 International Business Machines Corporation Automatic computer room air conditioning control method
WO2015130786A1 (en) * 2014-02-28 2015-09-03 Delos Living Llc Systems, methods and articles for enhancing wellness associated with habitable environments
WO2015134655A2 (en) * 2014-03-05 2015-09-11 Adeptdc Co. Systems and methods for intelligent controls for optimal resource allocation for data center operations
WO2015171624A1 (en) 2014-05-05 2015-11-12 Vigilent Corporation Point-based risk score for managing environmental systems
US10126009B2 (en) 2014-06-20 2018-11-13 Honeywell International Inc. HVAC zoning devices, systems, and methods
US20160048399A1 (en) * 2014-08-15 2016-02-18 At&T Intellectual Property I, L.P. Orchestrated sensor set
FR3031401B1 (fr) * 2015-01-06 2017-07-14 Ubiant Sa Systeme de gestion de la consommation energetique d'un batiment
WO2016115230A1 (en) 2015-01-13 2016-07-21 Delos Living Llc Systems, methods and articles for monitoring and enhancing human wellness
EP3256913A1 (en) 2015-02-11 2017-12-20 NEC Europe Ltd. A method for operating a thermal system and a thermal system
US10078610B2 (en) * 2015-05-04 2018-09-18 Dell Products, L.P. System and method for optimized thermal control for management controller offline
US11153183B2 (en) * 2015-06-11 2021-10-19 Instana, Inc. Compacted messaging for application performance management system
KR102053841B1 (ko) * 2015-07-20 2019-12-09 주식회사 엘지화학 Hvac 제어 시스템 및 이를 이용한 hvac 제어 방법
US10174957B2 (en) * 2015-07-27 2019-01-08 Mitsubishi Electric Research Laboratories, Inc. System and method for controlling multi-zone vapor compression systems
US10186889B2 (en) * 2015-10-08 2019-01-22 Taurus Des, Llc Electrical energy storage system with variable state-of-charge frequency response optimization
US11210617B2 (en) 2015-10-08 2021-12-28 Johnson Controls Technology Company Building management system with electrical energy storage optimization based on benefits and costs of participating in PDBR and IBDR programs
US9568923B1 (en) 2015-10-27 2017-02-14 International Business Machines Corporation Determining a time for corrective action in a data center
US10101050B2 (en) 2015-12-09 2018-10-16 Google Llc Dispatch engine for optimizing demand-response thermostat events
US10502470B2 (en) 2016-03-22 2019-12-10 Vertiv Corporation System and method to maintain evaporator superheat during pumped refrigerant economizer operation
US10215665B2 (en) 2016-05-03 2019-02-26 General Electric Company System and method to model power output of an engine
WO2018039433A1 (en) 2016-08-24 2018-03-01 Delos Living Llc Systems, methods and articles for enhancing wellness associated with habitable environments
EP3523701A1 (en) 2016-10-07 2019-08-14 Vertiv Corporation Adaptive pid control for chilled water crac units
US10337753B2 (en) * 2016-12-23 2019-07-02 Abb Ag Adaptive modeling method and system for MPC-based building energy control
US11076509B2 (en) 2017-01-24 2021-07-27 The Research Foundation for the State University Control systems and prediction methods for it cooling performance in containment
US10458669B2 (en) 2017-03-29 2019-10-29 Johnson Controls Technology Company Thermostat with interactive installation features
CN107036237B (zh) * 2017-03-30 2020-04-10 青岛海尔空调器有限总公司 空调器制热控制方法和控制装置
US10712038B2 (en) 2017-04-14 2020-07-14 Johnson Controls Technology Company Multi-function thermostat with air quality display
US10731885B2 (en) 2017-04-14 2020-08-04 Johnson Controls Technology Company Thermostat with occupancy detection via proxy measurements of a proxy sensor
WO2018191688A2 (en) 2017-04-14 2018-10-18 Johnson Controls Techology Company Thermostat with exhaust fan control for air quality and humidity control
WO2018191699A1 (en) 2017-04-14 2018-10-18 Johnson Controls Technology Company Multi-function thermostat with intelligent supply fan control for maximizing air quality and optimizing energy usage
US10866003B2 (en) 2017-04-14 2020-12-15 Johnson Controls Technology Company Thermostat with preemptive heating, cooling, and ventilation in response to elevated occupancy detection via proxy
US10837665B2 (en) 2017-04-14 2020-11-17 Johnson Controls Technology Company Multi-function thermostat with intelligent ventilator control for frost/mold protection and air quality control
US10921010B2 (en) * 2017-06-21 2021-02-16 Johnson Controls Technology Company Building management system with simulation and user action reinforcement machine learning
US11668481B2 (en) 2017-08-30 2023-06-06 Delos Living Llc Systems, methods and articles for assessing and/or improving health and well-being
US10660241B2 (en) 2017-09-06 2020-05-19 Vertiv Corporation Cooling unit energy optimization via smart supply air temperature setpoint control
CN107859646B (zh) * 2017-10-27 2019-11-29 英业达科技有限公司 温度控制装置及其控制方法
TWI648610B (zh) * 2017-11-09 2019-01-21 向暘科技股份有限公司 小型送風機溫度調節的控制方法
CN107701483A (zh) * 2017-11-17 2018-02-16 英业达科技有限公司 温度控制装置及其方法
WO2019102404A1 (en) * 2017-11-23 2019-05-31 L&T Technology Services Limited Multi-zone flexi-positioning air-conditioning system
US11131474B2 (en) 2018-03-09 2021-09-28 Johnson Controls Tyco IP Holdings LLP Thermostat with user interface features
CN110469925B (zh) 2018-05-10 2022-11-29 开利公司 Hvac系统及其控制方法
EP3850458A4 (en) 2018-09-14 2022-06-08 Delos Living, LLC AIR CLEANING SYSTEMS AND PROCEDURES
US10876752B2 (en) * 2018-12-27 2020-12-29 Trane International Inc. Discharge air reset for a constant volume air changeover bypass
US10876756B2 (en) * 2018-12-27 2020-12-29 Trane International Inc. Discharge air reset for a constant volume air changeover bypass
WO2020176503A1 (en) 2019-02-26 2020-09-03 Delos Living Llc Method and apparatus for lighting in an office environment
WO2020198183A1 (en) 2019-03-25 2020-10-01 Delos Living Llc Systems and methods for acoustic monitoring
US11209808B2 (en) * 2019-05-21 2021-12-28 At&T Intellectual Property I, L.P. Systems and method for management and allocation of network assets
US11639804B2 (en) 2019-12-13 2023-05-02 Trane International Inc. Automated testing of HVAC devices
US12004327B2 (en) * 2019-12-27 2024-06-04 Solaredge Technologies Korea Limited Company Containerized HVAC control
US11751360B2 (en) * 2020-03-17 2023-09-05 International Business Machines Corporation Intelligently deployed cooling fins
US11732916B2 (en) 2020-06-08 2023-08-22 Emerson Climate Technologies, Inc. Refrigeration leak detection
US11359846B2 (en) 2020-07-06 2022-06-14 Emerson Climate Technologies, Inc. Refrigeration system leak detection
US11885516B2 (en) 2020-08-07 2024-01-30 Copeland Lp Refrigeration leak detection
US11754324B2 (en) 2020-09-14 2023-09-12 Copeland Lp Refrigerant isolation using a reversing valve
US11609032B2 (en) 2020-10-22 2023-03-21 Emerson Climate Technologies, Inc. Refrigerant leak sensor measurement adjustment systems and methods
EP4012421A1 (en) * 2020-12-11 2022-06-15 Nxp B.V. System and method for monitoring feedback circuits
US11940188B2 (en) 2021-03-23 2024-03-26 Copeland Lp Hybrid heat-pump system
US11761663B2 (en) * 2021-11-19 2023-09-19 Johnson Controls Tyco IP Holdings LLP HVAC system for reducing intra-space variation of controlled environmental conditions
CN114857747A (zh) * 2022-03-16 2022-08-05 湖北合合能源科技发展有限公司 一种基于LoRa技术的中央空调智慧能源管理系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5416081A (en) * 1977-07-06 1979-02-06 Toshiba Corp Control of equipment
JP2001227797A (ja) * 2000-02-16 2001-08-24 Daikin Ind Ltd 蓄熱空調システムにおける負荷予測方法およびその装置
JP2004293844A (ja) * 2003-03-26 2004-10-21 Hitachi Plant Eng & Constr Co Ltd 空調設備
JP2007206060A (ja) * 2005-12-21 2007-08-16 General Electric Co <Ge> 制約問題に対する提案される解のロバストネスを評価し、制約問題の解を展開する際にロバストネスを考慮する方法および装置
US7676280B1 (en) * 2007-01-29 2010-03-09 Hewlett-Packard Development Company, L.P. Dynamic environmental management
JP2010523167A (ja) * 2007-01-31 2010-07-15 メドトロニック・ミニメッド・インコーポレーテッド インスリン注入を制御および監視するモデル予測の方法およびシステム
JP2010181082A (ja) * 2009-02-05 2010-08-19 Yamatake Corp モデル関数処理装置および方法
JP2013511694A (ja) * 2009-11-18 2013-04-04 ペースコントロールズ,エルエルシー デューティサイクル化されたhvac&r機器の自動制御及び最適化のためのコントローラー、並びに該コントローラーを用いるシステム及び方法
JP5602878B2 (ja) * 2009-12-31 2014-10-08 アーベーベー・リサーチ・リミテッド 発電プラントに対する負荷スケジューリングを最適化する方法及び制約解析モジュール

Family Cites Families (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4873649A (en) 1988-06-10 1989-10-10 Honeywell Inc. Method for operating variable speed heat pumps and air conditioners
US5233540A (en) * 1990-08-30 1993-08-03 The Boeing Company Method and apparatus for actively reducing repetitive vibrations
JPH04270854A (ja) 1991-02-26 1992-09-28 Hitachi Ltd 空気調和機の制御方法
US5170935A (en) 1991-11-27 1992-12-15 Massachusetts Institute Of Technology Adaptable control of HVAC systems
US5303080A (en) * 1993-04-01 1994-04-12 Eastman Kodak Company Beam scanning system including actively-controlled optical head
JPH06323595A (ja) 1993-05-12 1994-11-25 Daikin Ind Ltd 空気調和装置の運転制御装置
WO1995001592A1 (en) * 1993-06-30 1995-01-12 Ford Motor Company Method and system for controlling an automotive hvac system
US5464369A (en) 1994-02-25 1995-11-07 Johnson Service Company Method and apparatus for estimating the rate at which a gas is generated within an enclosed space
US5610843A (en) * 1995-03-01 1997-03-11 Sri International Methods and apparatuses for multi input/multi output control systems
US5875109A (en) 1995-05-24 1999-02-23 Johnson Service Company Adaptive flow controller for use with a flow control system
US5768121A (en) 1995-05-24 1998-06-16 Johnson Service Company Adaptive flow controller for use with a flow control system
US5920478A (en) 1997-06-27 1999-07-06 Oakleaf Engineering, Inc. Multi-input multi-output generic non-interacting controller
US6101459A (en) 1997-08-15 2000-08-08 Compaq Computer Corporation System and associated method for cooling components within a computer system
US5862982A (en) 1997-09-24 1999-01-26 Johnson Service Company Optimal ventilation control strategy
US6470230B1 (en) 2000-01-04 2002-10-22 Advanced Micro Devices, Inc. Supervisory method for determining optimal process targets based on product performance in microelectronic fabrication
US6557574B2 (en) 2000-04-04 2003-05-06 Clifford C. Federspiel Pressure based flow rate measurement device integrated with blades of a damper
US7097111B2 (en) 2000-07-21 2006-08-29 Gun Valley Temperature Controls Llc Environmental control system and method for storage buildings
US7797062B2 (en) * 2001-08-10 2010-09-14 Rockwell Automation Technologies, Inc. System and method for dynamic multi-objective optimization of machine selection, integration and utilization
US9729639B2 (en) * 2001-08-10 2017-08-08 Rockwell Automation Technologies, Inc. System and method for dynamic multi-objective optimization of machine selection, integration and utilization
US6719625B2 (en) 2001-09-26 2004-04-13 Clifford Conrad Federspiel Method and apparatus for controlling variable air volume supply fans in heating, ventilating, and air-conditioning systems
US6574104B2 (en) 2001-10-05 2003-06-03 Hewlett-Packard Development Company L.P. Smart cooling of data centers
US6402043B1 (en) 2001-10-18 2002-06-11 John F. Cockerill Method for controlling HVAC units
US6718277B2 (en) 2002-04-17 2004-04-06 Hewlett-Packard Development Company, L.P. Atmospheric control within a building
US7089087B2 (en) 2002-05-17 2006-08-08 Carrier Corporation Limited access comfort control
US6865449B2 (en) 2002-05-17 2005-03-08 Carrier Corporation Location adjusted HVAC control
US6775997B2 (en) 2002-10-03 2004-08-17 Hewlett-Packard Development Company, L.P. Cooling of data centers
US7035763B2 (en) * 2003-09-03 2006-04-25 Siemens Westinghouse Power Corporation Systems and methods for selecting training data and generating fault models for use in use sensor-based monitoring
US7664573B2 (en) 2003-09-26 2010-02-16 Siemens Industry, Inc. Integrated building environment data system
US7194337B2 (en) 2003-10-31 2007-03-20 Hewlett-Packard Development Company, Lp Agent-based operation of a robotic device
DE602005026178D1 (de) 2004-11-09 2011-03-10 Truveon Corp Gebäudeklimatisierungssteuerungsverfahren und -system
US7058477B1 (en) 2004-11-23 2006-06-06 Howard Rosen Thermostat system with remote data averaging
US20060116067A1 (en) 2004-12-01 2006-06-01 Federspiel Clifford C Method and apparatus for determining critical pressure of variable air volume heating, ventilating, and air-conditioning systems
JP2006207929A (ja) * 2005-01-28 2006-08-10 Daikin Ind Ltd 空調システムの最適運転制御システムおよび最適運転制御方法
EP1856453B1 (en) 2005-03-10 2016-07-13 Aircuity Incorporated Dynamic control of dilution ventilation in one-pass, critical environments
US7117129B1 (en) 2005-03-11 2006-10-03 Hewlett-Packard Development Company, L.P. Commissioning of sensors
US7894943B2 (en) * 2005-06-30 2011-02-22 Sloup Charles J Real-time global optimization of building setpoints and sequence of operation
US7726144B2 (en) 2005-10-25 2010-06-01 Hewlett-Packard Development Company, L.P. Thermal management using stored field replaceable unit thermal information
US7363094B2 (en) 2006-01-09 2008-04-22 General Electric Company Multivariable controller design method for multiple input/outputs systems with multiple input/output constraints
JP4661640B2 (ja) * 2006-03-09 2011-03-30 株式会社日立製作所 空調制御システム
US7689296B2 (en) * 2006-04-28 2010-03-30 Honeywell Asca Inc. Apparatus and method for controlling a paper machine or other machine using measurement predictions based on asynchronous sensor information
US7890215B2 (en) * 2006-12-22 2011-02-15 Duncan Scot M Optimized control system for cooling systems
US7847681B2 (en) 2007-03-23 2010-12-07 Johnson Controls Technology Company Building automation systems and methods
WO2009111489A2 (en) 2008-03-03 2009-09-11 Federspiel Corporation Methods and systems for coordinating the control of hvac units
US7966152B2 (en) 2008-04-23 2011-06-21 Honeywell International Inc. System, method and algorithm for data-driven equipment performance monitoring
US8374731B1 (en) 2008-12-24 2013-02-12 Emc Corporation Cooling system
US8255085B2 (en) * 2009-02-05 2012-08-28 Johnson Controls Technology Company Asymmetrical control system and method for energy savings in buildings
WO2011022696A1 (en) 2009-08-21 2011-02-24 Federspiel Corporation Method and apparatus for efficiently coordinating data center cooling units
US8639651B2 (en) * 2009-10-30 2014-01-28 Hewlett-Packard Development Company, L. P. Manipulating environmental conditions in an infrastructure
US20110161059A1 (en) 2009-12-30 2011-06-30 Ankur Jain Method for Constructing a Gray-Box Model of a System Using Subspace System Identification
JP2013543569A (ja) * 2010-08-20 2013-12-05 ヴィジレント コーポレイション Hvacシステムに対するエネルギー最適制御決定

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5416081A (en) * 1977-07-06 1979-02-06 Toshiba Corp Control of equipment
JP2001227797A (ja) * 2000-02-16 2001-08-24 Daikin Ind Ltd 蓄熱空調システムにおける負荷予測方法およびその装置
JP2004293844A (ja) * 2003-03-26 2004-10-21 Hitachi Plant Eng & Constr Co Ltd 空調設備
JP2007206060A (ja) * 2005-12-21 2007-08-16 General Electric Co <Ge> 制約問題に対する提案される解のロバストネスを評価し、制約問題の解を展開する際にロバストネスを考慮する方法および装置
US7676280B1 (en) * 2007-01-29 2010-03-09 Hewlett-Packard Development Company, L.P. Dynamic environmental management
JP2010523167A (ja) * 2007-01-31 2010-07-15 メドトロニック・ミニメッド・インコーポレーテッド インスリン注入を制御および監視するモデル予測の方法およびシステム
JP2010181082A (ja) * 2009-02-05 2010-08-19 Yamatake Corp モデル関数処理装置および方法
JP2013511694A (ja) * 2009-11-18 2013-04-04 ペースコントロールズ,エルエルシー デューティサイクル化されたhvac&r機器の自動制御及び最適化のためのコントローラー、並びに該コントローラーを用いるシステム及び方法
JP5602878B2 (ja) * 2009-12-31 2014-10-08 アーベーベー・リサーチ・リミテッド 発電プラントに対する負荷スケジューリングを最適化する方法及び制約解析モジュール

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017110904A (ja) * 2010-08-20 2017-06-22 ヴィジレント コーポレイションVigilent Corporation 環境維持システムの制御
JP2017506734A (ja) * 2014-03-28 2017-03-09 三菱電機株式会社 蒸気圧縮システムを制御するための極値探索コントローラー及び方法
WO2020121551A1 (ja) * 2018-12-12 2020-06-18 日本電信電話株式会社 マルチデバイス連携制御装置、マルチデバイス連携制御方法およびマルチデバイス連携制御プログラム、並びに、学習装置、学習方法および学習プログラム
JPWO2020121551A1 (ja) * 2018-12-12 2021-09-02 日本電信電話株式会社 マルチデバイス連携制御装置、マルチデバイス連携制御方法およびマルチデバイス連携制御プログラム、並びに、学習装置、学習方法および学習プログラム
JP7088313B2 (ja) 2018-12-12 2022-06-21 日本電信電話株式会社 マルチデバイス連携制御装置、マルチデバイス連携制御方法およびマルチデバイス連携制御プログラム
US11874634B2 (en) 2018-12-12 2024-01-16 Nippon Telegraph And Telephone Corporation Multi-device coordination control device, multi-device coordinaton control method, and multi-device coordination control program, and learning device, learning method, and learning program
JPWO2021059412A1 (ja) * 2019-09-25 2021-04-01
WO2021059412A1 (ja) * 2019-09-25 2021-04-01 日本電信電話株式会社 空調システム及び空調制御方法
JP7347519B2 (ja) 2019-09-25 2023-09-20 日本電信電話株式会社 空調システム及び空調制御方法
WO2023176151A1 (ja) * 2022-03-17 2023-09-21 トリニティ工業株式会社 空調制御装置及びそのためのプログラム
CN115978722A (zh) * 2023-03-17 2023-04-18 四川港通医疗设备集团股份有限公司 一种基于人工智能的净化机组监测管理方法及系统

Also Published As

Publication number Publication date
US9291358B2 (en) 2016-03-22
EP2606406A2 (en) 2013-06-26
JP2017110904A (ja) 2017-06-22
JP6514242B2 (ja) 2019-05-15
EP2606406A4 (en) 2017-10-25
US20150100165A1 (en) 2015-04-09
WO2012024692A2 (en) 2012-02-23
CA2808503A1 (en) 2012-02-23
KR20130098346A (ko) 2013-09-04
US8924026B2 (en) 2014-12-30
EP2606406B1 (en) 2020-04-01
US20120101648A1 (en) 2012-04-26
WO2012024692A3 (en) 2012-05-18

Similar Documents

Publication Publication Date Title
JP6514242B2 (ja) 環境維持システムの制御方法、環境維持システム及びプログラム
JP5820375B2 (ja) データセンタ冷却ユニットを効率的に調整する方法及び装置
JP7184797B2 (ja) ビルディング機器用のモデル予測的メンテナンスシステム
JP6254953B2 (ja) 環境維持モジュールを制御するための多次元最適化
KR101973652B1 (ko) Hvac 시스템에 대한 가변 공기 용적 모델링
US10108154B2 (en) Influence learning for managing physical conditions of an environmentally controlled space by utilizing a calibration override which constrains an actuator to a trajectory
CN108476473B (zh) 用于延长建筑物控制系统中的无线传感器的电池寿命的系统和方法
JP6313217B2 (ja) Hvacユニットの気温制御
US20150178421A1 (en) Systems for and methods of modeling, step-testing, and adaptively controlling in-situ building components
JP2018181331A (ja) 制約ハンドリングを伴う極値探索制御システム
US11408626B2 (en) Central plant control system with dynamic computation reduction
Kumar et al. Data center air handling unit fan speed optimization using machine learning techniques

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140814

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20151005

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20151215

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160401

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20160909