(57)【特許請求の範囲】
【請求項１】
少なくとも１つの発電設備と少なくとも１つの排熱利用機器とを含むコジェネレーション・システムを含む熱源機器ネットワークにおいて、該コジェネレーション・システムの該排熱利用機器の機器特性モデルを生成する機器特性モデル生成装置であって、
該少なくとも１つの発電設備の負荷率と、該少なくとも１つの排熱利用機器の負荷率及び複数の離散的な負荷率における入出力特性を表す値とを取得する第１の入出力特性取得部と、
該少なくとも１つの排熱利用機器の該入出力特性を該発電設備の負荷率及び該排熱利用機器の負荷率の関数として表した、該排熱利用機器の機器特性モデルを生成する排熱利用機器特性モデル生成部とを含む、機器特性モデル生成装置。
【請求項２】
請求項１に記載の機器特性モデル生成装置において、
該少なくとも１つの発電設備の負荷率及び複数の離散的な負荷率における入出力特性を表す値を取得する第２の入出力特性取得部と、
該少なくとも１つの発電設備の該入出力特性を該発電設備の該負荷率の関数として表した、該発電設備の機器特性モデルを生成する発電設備機器特性モデル生成部とをさらに含む、機器特性モデル生成装置。
【請求項３】
請求項１に記載の機器特性モデル生成装置において、
該熱源機器ネットワークは、複数の種類の熱源機器をさらに含み、該機器特性モデル生成装置は、
該複数の種類の熱源機器のうちの少なくとも１つについて、該熱源機器の負荷率、冷却水温の予測値、及び複数の離散的な冷却水温値の各値ごとに複数の離散的な負荷率における入出力特性を表す値を取得する第３の入出力特性取得部と、
該少なくとも１つの熱源機器の入出力特性を表す値を冷却水温及び負荷率の関数として表した、該熱源機器の機器特性モデルを生成する熱源機器特性モデル生成部とを含む、機器特性モデル生成装置。
【請求項４】
請求項２に記載の機器特性モデル生成装置において、
該発電設備機器特性モデル生成部は、
該少なくとも１つの発電設備の総排熱量と、該少なくとも１つの発電設備から該少なくとも１つの排熱利用機器に送られる回収可能な排熱量を、該少なくとも１つの発電設備の該負荷率の関数として決定する排熱量決定部を含み、
該排熱利用機器特性モデル生成部は、
該回収可能な排熱量と該少なくとも１つの排熱利用機器の排熱回収容量とに基づいて、該少なくとも１つの発電設備の負荷率と該少なくとも１つの排熱利用機器の排熱利用分の負荷率の関数として該少なくとも１つの排熱利用機器の実際の排熱回収量を決定する排熱回収量決定部と、
該少なくとも１つの排熱利用機器の排熱利用分の負荷率の関数として、該少なくとも１つの排熱利用機器の排熱利用分の機器特性モデルを生成する排熱利用分機器特性モデル生成部と、
該熱源機器ネットワーク全体の熱需要のうち該少なくとも１つの排熱利用機器に割り当てられた熱需要と、該少なくとも１つの排熱利用機器の該排熱利用分の機器特性モデルと、該実際の排熱回収量とに基づいて、該少なくとも１つの排熱利用機器の燃料利用分の負荷率を算出し、該燃料利用分の負荷率の関数として、該少なくとも１つの排熱利用機器の燃料利用分の機器特性モデルを生成する燃料利用分機器特性モデル生成部とを含む、機器特性モデル生成装置。
【請求項５】
請求項３に記載の機器特性モデル生成装置において、
該熱源機器特性モデル生成部は、
該複数の離散的な冷却水温値の各値ごとの該複数の種類の熱源機器の入出力特性を表す値を該第３の入出力特性取得部から受け取り、該複数の離散的な冷却水温値の各値ごとにパラメータの異なるセットを使用して、該複数の種類の熱源機器の入出力特性を表す値を、該複数の離散的な冷却水温値の各値ごとの負荷率の連続的な関数として表す関数推定部と、
該複数の離散的な冷却水温値の各値ごとに異なる該パラメータのセットを該関数推定部から受け取り、該冷却水温値の各値ごとに異なる該パラメータのセットの中の各パラメータ値を冷却水温の連続的な関数として表すパラメータ推定部と、
該入出力特性を表す該冷却水温値の各値ごとの該負荷率の連続的な関数に関する情報を該関数推定部から、及び各パラメータ値を表す該冷却水温の連続的な関数に関する情報を該パラメータ推定部からそれぞれ受け取り、該冷却水温の連続的な関数として表した各パラメータ値を、該入出力特性を表す該冷却水温値の各値ごとの該負荷率の連続的な関数に適用することにより、該入出力特性を冷却水温及び負荷率の連続的な関数として表す関数モデル決定部とを含む、機器特性モデル生成装置。
【請求項６】
請求項４に記載の機器特性モデルによってモデル化されるコジェネレーション・システムにおいて、
該少なくとも１つの発電設備は、
該少なくとも１つの発電設備の発電量と発電機の効率と該総排熱量とに基づいて、該少なくとも１つの発電設備の所要燃料を算出する発電設備所要燃料計算部と、
該算出された所要燃料と該少なくとも１つの発電設備の補機類駆動用電力と電力単価と燃料単価とに基づいて、該少なくとも１つの発電設備のコストを算出する発電設備コスト計算部とを含み、
該少なくとも１つの排熱利用機器は、
該熱源機器ネットワーク全体の熱需要のうち該少なくとも１つの排熱利用機器に割り当てられた熱需要と、該少なくとも１つの排熱利用機器の該排熱利用分の機器特性モデルと、該実際の排熱回収量とに基づいて、該少なくとも１つの排熱利用機器の所要燃料を算出する排熱利用機器所要燃料計算部と、
該算出された所要燃料と該少なくとも１つの排熱利用機器の補機類駆動用電力と燃料単価と電力単価とに基づいて、該少なくとも１つの排熱利用機器のコストを算出する排熱利用機器コスト計算部とを含む、コジェネレーション・システム。
【請求項７】
請求項５に記載の関数モデル決定部によって冷却水温及び負荷率の連続的な関数として表される該入出力特性を用いてモデル化される該熱源機器において、
該冷却水温及び該負荷率の連続的な関数として表された該入出力特性と、該熱源機器全体の熱需要のうち該熱源機器に割り当てられた熱需要である熱量設定値とに基づいて、該熱源機器の所要燃料を算出する熱源機器所要燃料計算部と、
該算出された所要燃料と該熱源機器の補機類駆動用電力と燃料単価と電力単価とに基づいて、該熱源機器のコストを算出する熱源機器コスト計算部とを含む、熱源機器。
【請求項８】
なくとも１つの発電設備と少なくとも１つの排熱利用機器とを含むコジェネレーション・システムを含む熱源機器ネットワークにおいて、該コジェネレーション・システムの該排熱利用機器の機器特性モデルを生成する機器特性モデル生成方法であって、
該少なくとも１つの発電設備の負荷率と、該少なくとも１つの排熱利用機器の負荷率及び複数の離散的な負荷率における入出力特性を表す値とを取得するステップと、
該少なくとも１つの排熱利用機器の該入出力特性を該発電設備の負荷率及び該排熱利用機器の負荷率の関数として表した、該排熱利用機器の機器特性モデルを生成するステップとを含む、機器特性モデル生成方法。
【請求項９】
請求項８に記載の機器特性モデル生成方法において、
該少なくとも１つの発電設備の負荷率及び複数の離散的な負荷率における入出力特性を表す値を取得するステップと、
該少なくとも１つの発電設備の該入出力特性を該発電設備の該負荷率の関数として表した、該発電設備の機器特性モデルを生成するステップとをさらに含む、機器特性モデル生成方法。
【請求項１０】
請求項８に記載の機器特性モデル生成方法において、
該複数の種類の熱源機器のうちの少なくとも１つについて複数の離散的な冷却水温値の各値ごとに複数の離散的な負荷率における入出力特性を表す値を取得するステップと、
該少なくとも１つの熱源機器の入出力特性を表す値を冷却水温及び負荷率の関数として表した、該熱源機器の機器特性モデルを生成するステップとをさらに含む、機器特性モデル生成方法。
【請求項１１】
請求項９に記載の機器特性モデル生成方法において、
該発電設備の機器特性モデルを生成するステップは、
該少なくとも１つの発電設備の総排熱量と、該少なくとも１つの発電設備から該少なくとも１つの排熱利用機器に送られる回収可能な排熱量を、該少なくとも１つの発電設備の該負荷率の関数として決定するステップを含み、
該排熱利用機器の機器特性モデルを生成するステップは、
該回収可能な排熱量と該少なくとも１つの排熱利用機器の排熱回収容量とに基づいて、該少なくとも１つの発電設備の該負荷率と該少なくとも１つの排熱利用機器の排熱利用分の負荷率の関数として該少なくとも１つの排熱利用機器の実際の排熱回収量を決定するステップと、
該少なくとも１つの排熱利用機器の排熱利用分の負荷率の関数として、該少なくとも１つの排熱利用機器の排熱利用分の機器特性モデルを生成するステップと、
該熱源機器ネットワーク全体の熱需要のうち該少なくとも１つの排熱利用機器に割り当てられた熱需要と、該少なくとも１つの排熱利用機器の該排熱利用分の機器特性モデルと、該実際の排熱回収量とに基づいて、該少なくとも１つの排熱利用機器の燃料利用分の負荷率を算出し、該燃料利用分の負荷率の関数として、該少なくとも１つの排熱利用機器の燃料利用分の機器特性モデルを生成するステップとを含む、機器特性モデル生成方法。
【請求項１２】
請求項１０に記載の機器特性モデル生成方法において、
該熱源機器の機器特性モデルを生成するステップは、
該複数の離散的な冷却水温値の各値ごとの該複数の種類の熱源機器の入出力特性を表す値を受け取り、該複数の離散的な冷却水温値の各値ごとにパラメータの異なるセットを使用して、該複数の種類の熱源機器の入出力特性を表す値を、該複数の離散的な冷却水温値の各値ごとの負荷率の連続的な関数として表すステップと、
該複数の離散的な冷却水温値の各値ごとに異なる該パラメータのセットを受け取り、該冷却水温値の各値ごとに異なる該パラメータのセットの中の各パラメータ値を冷却水温の連続的な関数として表すステップと、
該入出力特性を表す該冷却水温値の各値ごとの該負荷率の連続的な関数に関する情報、及び各パラメータ値を表す該冷却水温の連続的な関数に関する情報をそれぞれ受け取り、該冷却水温の連続的な関数として表した各パラメータ値を、該入出力特性を表す該冷却水温値の各値ごとの該負荷率の連続的な関数に適用することにより、該入出力特性を冷却水温及び負荷率の連続的な関数として表すステップとを含む、機器特性モデル生成方法。
【請求項１３】
請求項１１に記載の機器特性モデル生成方法において、
該少なくとも１つの発電設備の発電量と発電機の効率と該総排熱量とに基づいて、該少なくとも１つの発電設備の所要燃料を算出するステップと、
該算出された所要燃料と該少なくとも１つの発電設備の補機類駆動用電力と電力単価と燃料単価とに基づいて、該少なくとも１つの発電設備のコストを算出するステップと、
該熱源機器ネットワーク全体の熱需要のうち該少なくとも１つの排熱利用機器に割り当てられた熱需要と、該少なくとも１つの排熱利用機器の該排熱利用分の機器特性モデルと、該実際の排熱回収量とに基づいて、該少なくとも１つの排熱利用機器の所要燃料を算出するステップと、
該算出された所要燃料と該少なくとも１つの排熱利用機器の補機類駆動用電力と燃料単価と電力単価とに基づいて、該少なくとも１つの排熱利用機器のコストを算出するステップとを含む、機器特性モデル生成方法。
【請求項１４】
請求項１２に記載の機器特性モデル生成方法において、
該冷却水温及び該負荷率の連続的な関数として表された該入出力特性と、該熱源機器全体の熱需要のうち該熱源機器に割り当てられた熱需要である熱量設定値とに基づいて、該熱源機器の所要燃料を算出するステップと、
該算出された所要燃料と該熱源機器の補機類駆動用電力と燃料単価と電力単価とに基づいて、該熱源機器のコストを算出するステップとを含む、機器特性モデル生成方法。
【請求項１５】
少なくとも１つの発電設備と少なくとも１つの排熱利用機器とを含むコジェネレーション・システムを含む熱源機器ネットワークにおいて、該コジェネレーション・システムの該排熱利用機器の機器特性モデルを生成する方法を機器特性モデル生成装置に実行させる機器特性モデル生成プログラムであって、
該機器特性モデル生成装置は、第１の入出力特性取得部と排熱利用機器特性モデル生成部とを含み、該方法は、
該第１の入出力特性取得部が、該少なくとも１つの発電設備の負荷率と、該少なくとも１つの排熱利用機器の負荷率及び複数の離散的な負荷率における入出力特性を表す値とを取得するステップと、
該排熱利用機器特性モデル生成部が、該少なくとも１つの排熱利用機器の該入出力特性を該発電設備の負荷率及び該排熱利用機器の負荷率の関数として表した、該排熱利用機器の機器特性モデルを生成するステップとを含む、機器特性モデル生成プログラム。
【請求項１６】
請求項１５に記載の機器特性モデル生成プログラムにおいて、
該機器特性モデル生成装置は、第２の入出力特性取得部と発電設備機器特性モデル生成部とをさらに含み、該方法は、
該第２の入出力特性取得部が、該少なくとも１つの発電設備の負荷率及び複数の離散的な負荷率における入出力特性を表す値を取得するステップと、
該発電設備機器特性モデル生成部が、該少なくとも１つの発電設備の該入出力特性を該発電設備の該負荷率の関数として表した、該発電設備の機器特性モデルを生成するステップとをさらに含む、機器特性モデル生成プログラム。
【請求項１７】
請求項１５に記載の機器特性モデル生成プログラムにおいて、
該機器特性モデル生成装置は、第３の入出力特性取得部と熱源機器特性モデル生成部とをさらに含み、該方法は、
該第３の入出力特性取得部が、該複数の種類の熱源機器のうちの少なくとも１つについて複数の離散的な冷却水温値の各値ごとに複数の離散的な負荷率における入出力特性を表す値を取得するステップと、
該熱源機器特性モデル生成部が、該少なくとも１つの熱源機器の入出力特性を表す値を冷却水温及び負荷率の関数として表した、該熱源機器の機器特性モデルを生成するステップとを含む、機器特性モデル生成プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、コジェネレーション・システムの機器特性を作成するための方法、システム、及びプログラムに関する。
【背景技術】
【０００２】
近年、環境問題と資源問題の双方の観点から、エンジン、タービン、燃料電池等の発電設備をエネルギーの需要地の近くに配置し、発電した電力と共に、発電に伴って発生する熱エネルギーを有効に利用するコジェネレーション・システム（Ｃｏ−Ｇｅｎｅｒａｔｉｏｎ  Ｓｙｓｔｅｍ（ＣＧＳ））が注目されている。現在、我が国の電力需給システムの主流は、電力会社の大規模集中発電によるものであるが、需要地に隣接して分散配置されるこのような小規模な発電システム全般は、「分散型エネルギーシステム」と呼ばれる。
【０００３】
一般的に、分散型エネルギーシステムの構成の問題は、最適配置問題と最適運用問題が複合した非常に複雑な問題となる。分散型エネルギーシステムの設計する際には、設計条件となるパラメータが多く、また時間ごとの煩雑なエネルギー計算を年間にわたって行う必要がある。特に、コジェネレーション・システムの場合には、発電を行なうためのガスエンジン、ガスタービン、燃料電池等の原動機と、発電時に発生する熱を利用する熱交換器、熱で駆動する吸収式冷温水器等の排熱利用機器から構成され、発電出力と熱の両方を有効に利用しなければ、本来の効率を発揮できないことから、コジェネレーション・システムのシステム構成を決定する問題は重要である。
【０００４】
この運用問題は、互いに関連する複数の機器の時間毎の運転スケジュール等を決める問題であり、熟練したオペレータでなければ、人手で運転計画を策定するのは不可能である。このため、分散型エネルギーシステムの運用問題においては、パーソナル・コンピュータ等による支援ツールを活用するのが望ましい。
【０００５】
例えば、特許文献１は、 電力供給源（ＥＰ）と都市ガス供給減（ＵＧ）からのエネルギーに基づいて、ガスエンジン（ＧＥ）３１とヒートポンプ（ＨＰ）３０と排熱利用吸収式冷凍機（ＡＲ）３２を用いて、電力負荷と冷房負荷と暖房負荷にエネルギーを供給するコジェネレーション・システムを教示している。特許文献１のコジェネレーション・システムは、制御処理プログラムを実行して、ハミルトンアルゴリズムに従って、ガスエンジン（ＧＥ）を含む発電機５１とヒートポンプ（ＨＰ）を含む電気式熱源５２と排熱利用吸収式冷凍機（ＡＲ）を含む排熱利用機５４とからなる建物エネルギーシステムのコストを表す評価関数を最小化する。その際に、ＧＥ、ＨＰ、及びＡＲそれぞれの効率と、ＧＥの排熱回収率をパラメータとして上記の評価関数に組み込んでいる。評価関数の最小化の際に、特許文献１のシステムは、複数の異なる種類の熱源機器の効率及び発電機の排熱回収率を入力パラメータとして使用して、評価関数の最小化又は最大化を行う。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０００‐１４６２５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、これらの入力パラメータとして使用されるＧＥ、ＨＰ、及びＡＲそれぞれの効率及びＧＥの排熱回収率は負荷率に応じて変動するものであるが、特許文献１に記載されたシステムでは、これらのパラメータは最適化手順の最初にキーボードから入力されており、負荷率に応じた変動は何ら考慮されていない。
【０００８】
従って、特許文献１のシステムでは、省エネ機器として普及が進んでいる排熱投入型吸収冷温水機の排熱回収特性が負荷率に依存することへの対応ができていない。
【課題を解決するための手段】
【０００９】
本発明の第１の態様によれば、本発明は、少なくとも１つの発電設備と少なくとも１つの排熱利用機器とを含むコジェネレーション・システムを含む熱源機器ネットワークにおいて、該コジェネレーション・システムの該排熱利用機器の機器特性モデルを生成する機器特性モデル生成装置であって、
該少なくとも１つの発電設備の負荷率と、該少なくとも１つの排熱利用機器の負荷率及び複数の離散的な負荷率における入出力特性を表す値とを取得する第１の入出力特性取得部と、
該少なくとも１つの排熱利用機器の該入出力特性を該発電設備の負荷率及び該排熱利用機器の負荷率の関数として表した、該排熱利用機器の機器特性モデルを生成する排熱利用機器特性モデル生成部とを含む、機器特性モデル生成装置である。
【００１０】
本発明の第２の態様によれば、本発明は、少なくとも１つの発電設備と少なくとも１つの排熱利用機器とを含むコジェネレーション・システムを含む熱源機器ネットワークにおいて、該コジェネレーション・システムの該排熱利用機器の機器特性モデルを生成する機器特性モデル生成方法であって、
該少なくとも１つの発電設備の負荷率と、該少なくとも１つの排熱利用機器の負荷率及び複数の離散的な負荷率における入出力特性を表す値とを取得するステップと、
該少なくとも１つの排熱利用機器の該入出力特性を該発電設備の負荷率及び該排熱利用機器の負荷率の関数として表した、該排熱利用機器の機器特性モデルを生成するステップとを含む、機器特性モデル生成方法である。
【００１１】
本発明の第３の態様によれば、本発明は、少なくとも１つの発電設備と少なくとも１つの排熱利用機器とを含むコジェネレーション・システムを含む熱源機器ネットワークにおいて、該コジェネレーション・システムの該排熱利用機器の機器特性モデルを生成する方法を機器特性モデル生成装置に実行させる機器特性モデル生成プログラムであって、
該機器特性モデル生成装置は、第１の入出力特性取得部と排熱利用機器特性モデル生成部とを含み、該方法は、
該第１の入出力特性取得部が、該少なくとも１つの発電設備の負荷率と、該少なくとも１つの排熱利用機器の負荷率及び複数の離散的な負荷率における入出力特性を表す値とを取得するステップと、
該排熱利用機器特性モデル生成部が、該少なくとも１つの排熱利用機器の該入出力特性を該発電設備の負荷率及び該排熱利用機器の負荷率の関数として表した、該排熱利用機器の機器特性モデルを生成するステップとを含む、機器特性モデル生成プログラムである。
【００１２】
その他、本願が開示する課題、及びその解決方法は、発明を実施するための形態の欄の記載、及び図面の記載等により明らかにされる。
【発明の効果】
【００１３】
排熱投入型吸収冷温水機は発電機からの排熱を利用しているため、発電機の機器特性と排熱投入型吸収冷温水機の機器特性をモデル化する必要があるが、本発明に係るシステム及び方法は、メーカーから提供される特性カーブを、負荷率を変数として、発電機における排熱量と排熱投入型吸収冷温水機における排熱回収量及び成績係数（ＣＯＰ）を回帰して係数を推計するように構成される。従って、メーカーから提供される特性データを用いて特性関数を推計してモデル化することで、任意の負荷率における特性値を容易に取得できるようになる。すなわち、負荷率がわかれば、精度の高い発電機排熱量を得た後、排熱回収量や成績係数を容易に得ることが可能となる。
【図面の簡単な説明】
【００１４】
【図１】本発明に従ったコジェネレーション・システムを含む熱源機器ネットワークを示す図である。
【図２】コジェネレーション・システムの機能ブロック図である。
【図３】電動冷凍機（ＥＲＵ）１０８の機能ブロック図である。
【図４】ガス吸収式冷温水機（ＧＡＢ）１１４の機能ブロック図である。
【図５】本願発明に従った熱源機器ネットワークの運転計画を最適化するためのシステムの概念図である。
【図６】需要予測部のハードウェア構成の機能ブロック図を示している。
【図７】需要予測部のソフトウェア構成を示す図である。
【図８】外気冷房（ナイトパージ）ブロック２０６の動作を示す概念図である。
【図９】比エンタルピ予測ブロック２０８及び外気冷房（ナイトパージ）ブロック２０６によって行われる排熱量予測の概念図である。
【図１０】外気冷房（ナイトパージ）制御装置を示す機能ブロック図である。
【図１１】建物外の気温の予測値、湿度の予測値、空気密度、及び比エンタルピ、建物内の気温の予測値、湿度の予測値、空気密度、比エンタルピ、換気量、及び排熱量の関係を示す図である。
【図１２】各熱源機器の冷却水温予測の概念図である。
【図１３】予測モデルを用いて予測した冷却水温Ｔ（ｉ）の予測値Ｔ＊（ｉ，ｔ）と実績値との比較を示す図である。
【図１４】発電設備（ＥＧＥ）１１０の機器特性モデルの作成及びコスト計算方法を示すブロック図である。
【図１５】排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率―排熱回収量特性及び負荷率―ＣＯＰ特性を示す図である。
【図１６】排熱投入型吸収冷温水機（ＥＸＨ）１１２のコスト計算方法を示すブロック図である。
【図１７】機器特性モデル化ブロック２１２を用いて電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器モデルを求める際の各機能のブロック図である。
【図１８】冷房機器のＣＯＰ特性を示す図である。
【図１９】電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器特性を示す図である。
【図２０】電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器特性を示す図である。
【図２１】電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器特性を示す図である。
【図２２】電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４のコストの計算の手順を示す機能ブロック図である。
【図２３】本発明に従った熱源機器ネットワークの運転計画の最適化方法のフローチャートである。
【図２４】本発明に従った熱源機器ネットワークの運転計画の最適化方法のフローチャートである。
【発明を実施するための形態】
【００１５】
本発明は様々な変更形態及び代替形態が可能であるが、本発明の一例としての実施形態が実施例として図面に示されており、本明細書で詳細に説明される。しかし、一例としての実施形態の本明細書における説明は、開示された特定の形態に本発明を限定することを意図するものではなく、本発明は、特許請求の範囲によって規定された本発明の趣旨及び範囲に入る全ての変更形態、同等形態、及び代替形態を含む。
【００１６】
（１）  熱源機器ネットワークの全体構成
図１は、本発明に従ったコジェネレーション・システムを含む熱源機器ネットワークを示す図である。
【００１７】
熱源機器ネットワーク１００は、受電設備１０２、燃料供給設備１０４、電気負荷１０６、電動冷凍機１０８、コジェネレーション・システム（ＣＧＳ）１２０、ガス吸収式冷温水機１１４、及び空調負荷１１６を含む。コジェネレーション・システム（ＣＧＳ）１２０は、発電設備（ＥＧＥ）１１０及び排熱投入型吸収冷温水機（ＥＸＨ）１１２を含む。
【００１８】
受電設備１０２は、ｔ番目（ｔ＝１，…，ｎ）の単位時間ｔにおいて、総計Ｐｒｅｃ（ｔ）の電力を外部から受け取り、電力Ｐｒｅｃｌ（ｔ）を電気負荷（ＥＬ）１０６に供給し、電力Ｐｒｅｃｒｅ（ｔ）を電動冷凍機（ＥＲＵ）１０８に供給する。燃料供給設備（ＦＳＦ）１０４は、合計ＦＳ（ｔ）の供給燃料のうち、燃料ＦＳＥＧＥ（ｔ）を発電設備（ＥＧＥ）１１０に供給し、燃料ＦＳＥＸＨ（ｔ）を排熱投入型吸収冷温水機（ＥＸＨ）１１２に供給し、燃料ＦＳＧＡＢ（ｔ）をガス吸収式冷温水機（ＧＡＢ）１１４に供給する。発電設備（ＥＧＥ）１１０は、発電電力ＰＥＧ（ｔ）のうち、電力ＰＥＧｌ（ｔ）を電気負荷（ＥＬ）１０６に供給し、電力ＰＥＧＥＲ（ｔ）を電動冷凍機（ＥＲＵ）１０８に供給する。また、発電設備（ＥＧＥ）１１０は、排熱ＨｈＥＧ（ｘＥＧＥ（ｔ））を排熱投入型吸収冷温水機（ＥＸＨ）１１２に供給する。電動冷凍機（ＥＲＵ）１０８は、冷水により冷熱ＨｃＥＲＵ（ｔ）を空調負荷１１６に供給し、ガス吸収式冷温水機（ＧＡＢ）１１４は、冷温水により冷熱ＨｃＧＡＢ（ｔ）を空調負荷１１６に供給する。排熱投入型吸収冷温水機（ＥＸＨ）１１２は、燃料供給設備（ＦＳＦ）１０４から供給された燃料ＦＳＥＸＨ（ｔ）と発電設備（ＥＧＥ）１１０から供給された排熱ＨｈＥＧ（ｔ）を使用して、冷温水により冷熱ＨｃＥＸＨ（ｔ）を空調負荷１１６に供給する。
【００１９】
動作中、発電設備（ＥＧＥ）１１０は、温度ＴＥＧＥ（ｔ）の冷却水で設備を冷却しつつ、指定された負荷率ｘＥＧＥ（ｔ）に応じた電力ＰＥＧ（ｔ）を生成する。また、発電設備（ＥＧＥ）１１０は、負荷率ｘＥＧＥ（ｔ）に応じた排熱ＨｈＥＧ（ｘＥＧＥ（ｔ））を排熱投入型吸収冷温水機（ＥＸＨ）１１２に供給する。電動冷凍機（ＥＲＵ）１０８は、温度ＴＥＲＵ（ｔ）の冷却水で設備を冷却しつつ、指定された負荷率ｘＥＲＵ（ｔ）に応じた冷熱ＨｃＥＲＵ（ｔ）を生成する。排熱投入型吸収冷温水機（ＥＸＨ）１１２は、温度ＴＥＸＨ（ｔ）の冷却水で設備を冷却しつつ、指定された負荷率ｘＥＸＨ（ｔ）に応じた冷熱ＨｃＥＸＨ（ｔ）を生成する。排熱投入型吸収冷温水機（ＥＸＨ）１１２の排熱回収特性は、負荷率ｘＥＸＨ（ｔ）に依存して決まる。ガス吸収式冷温水機（ＧＡＢ）１１４は、温度ＴＧＡＢ（ｔ）の冷却水で設備を冷却しつつ、指定された負荷率ｘＧＡＢ（ｔ）に応じた冷熱ＨｃＧＡＢ（ｔ）を生成する。また、各設備には冷却ポンプ等の補機が装備されており、これらの動力である電気は、受電設備及び発電設備で供給される。
【００２０】
単位時間ｔにおける電力需要をＰｄ（ｔ）、熱需要をＨｄ（ｔ）とすると、熱源機器ネットワーク１００は、
Ｐｒｅｃ（ｔ）＝Ｐｒｅｃｌ（ｔ）＋Ｐｒｅｃｒｅ（ｔ）                  （１）
ＰＥＧ（ｔ）＝ＰＥＧｌ（ｔ）＋ＰＥＧＥＲ（ｔ）                          （２）
ＦＳ（ｔ）＝ＦＳＥＧＥ（ｔ）＋ＦＳＥＸＨ（ｔ）＋ＦＳＧＡＢ（ｔ）      （３）
Ｐｒｅｃｌ（ｔ）＋Ｐｒｅｃｒｅ（ｔ）＋ＰＥＧｌ（ｔ）＋ＰＥＧＥＲ（ｔ）＝Ｐｄ（ｔ）  （４）
ＨｃＥＲＵ（ｔ）＋ＨｃＥＸＨ（ｔ）＋ＨｃＧＡＢ（ｔ）＝Ｈｄ（ｔ）      （５）
を満たす必要がある。
【００２１】
（２）  コジェネレーション・システム
図２は、コジェネレーション・システムの機能ブロック図である。図２において、それぞれ、図２（ａ）は、発電設備（ＥＧＥ）の機能ブロック図、図２（ｂ）は、排熱投入型吸収冷温水機（ＥＸＨ）の機能ブロック図である。
【００２２】
コジェネレーション・システムは、エンジン、タービン、燃料電池といった発電装置を需要地の近くに設置し、熱交換器や吸収式冷温水機といった発電時に発生する排熱を利用する排熱利用機器で発電装置の排熱を有効利用するシステムである。
【００２３】
（i）  発電設備
発電設備（ＥＧＥ）１１０は、機器諸元入力部１３０２、機器諸元格納部１３０４、負荷率計算部１３０６、発電量設定部１３０８、機器特性データ入力部１３１０、機器特性データ格納部１３１２、機器特性モデル生成部１３１４、排熱量計算部１３１６、所要燃料計算部１２１８、補機類稼働判定部１３２０、コスト等計算部１３２２、計算結果格納部１３２４を含む。
【００２４】
発電設備（ＥＧＥ）１１０を動作させる際に、最初に、機器諸元入力部１３０２に発電設備（ＥＧＥ）１１０の定格Ｐ（ｒａｔ）ＥＧを含む諸元を入力する。入力された諸元は、機器諸元格納部１３０４に格納される。同時に、発電量設定部１３０８に所望の発電量を入力すると、発電設備（ＥＧＥ）１１０の負荷率ｘＥＧＥ（ｔ）が負荷率計算部１３０６で計算される。負荷率ｘＥＧＥ（ｔ）を使用することによって、排熱量計算部１３１６で、総排熱量ＨｈＥＧ（ｇｒｏｓｓ）（ｘＥＧＥ（ｔ））及び回収可能な排熱量（すなわち、排熱投入型吸収冷温水機（ＥＸＨ）１１２に供給される排熱量）ＨｈＥＧ（ｘＥＧＥ（ｔ））が計算される。排熱量という絶対値を計算する際には機器の定格が必要となるため、機器諸元格納部１３０４から機器諸元を取得して、それぞれの排熱量を計算する。
【００２５】
発電機は、一般的に、負荷率ｘＥＧＥ（ｔ）に応じて効率（すなわち、損失）ηＥＧＥ（ｘＥＧＥ（ｔ））が定まるようになっており、効率ηＥＧＥ（ｘＥＧＥ（ｔ））が分かればそれぞれの排熱量がわかる仕組みになっている。総排熱量ＨｈＥＧ（ｇｒｏｓｓ）（ｘＥＧＥ（ｔ））が分かれば燃料消費量（所要燃料）ＦＳＥＧＥ（ｔ）がわかるので、所要燃料ＦＳＥＧＥ（ｔ）を、所要燃料計算部１３１８で計算しておく。計算された所要燃料に基づいてコストをコスト等計算部１３２２で計算する。コストを計算する際に、外部から電力料金単価、ガス料金単価、及びＣＯ２排出原単価を取得する。排熱量の計算結果とコストの計算結果は、計算結果格納部１３２４に格納される。回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））の計算結果は、排熱投入型吸収冷温水機（ＥＸＨ）１１２の排熱量入力部１３５６に送られる。
【００２６】
発電設備（ＥＧＥ）１１０は、補機類稼働判定部１３２０を含むが、補機類とは、例えば、発電設備（ＥＧＥ）１１０の冷却ポンプである。冷却ポンプの数は発電設備（ＥＧＥ）１１０の負荷率ｘＥＧＥ（ｔ）によって決まってくるので、負荷率計算部１３０６からの負荷率ｘＥＧＥ（ｔ）を使用することによって、補機類稼働判定部１３２０において補機類が稼働中か否かを判定する。
【００２７】
（ii）  排熱投入型吸収冷温水機
排熱投入型吸収冷温水機（ＥＸＨ）１１２は、機器諸元入力部１３５２、機器諸元格納部１３５４、排熱量入力部１３５６、負荷率計算部１３５８、熱量設定部１３６０、機器特性データ入力部１３６２、機器特性データ格納部１３６４、機器特性モデル生成部１３６６、排熱回収量計算部１３６８、補機類稼働判定部１３７０、ＣＯＰ計算部１３７２、所要燃料計算部１３７４、コスト計算部１３７６、計算結果格納部１３７８を含む。
【００２８】
排熱投入型吸収冷温水機（ＥＸＨ）１１２は、発電設備（ＥＧＥ）１１０から供給された発生した回収可能な排熱ＨｈＥＧ（ｘＥＧＥ（ｔ））を使用して冷房用の冷熱を生成する熱源機器である。排熱投入型吸収冷温水機（ＥＸＨ）１１２を動作させる際に、最初に、機器諸元入力部１３５２に排熱投入型吸収冷温水機（ＥＸＨ）１１２の定格値Ｈ（ｒａｔ）ｃＥＸＨを含む諸元を入力する。入力された諸元は、機器諸元格納部１３５４に格納される。同時に、熱需要Ｈｄ（ｔ）のうち排熱投入型吸収冷温水機（ＥＸＨ）１１２に割り当てられた熱需要ＨｃＥＸＨ（ｒｅｃ）（ｔ）を熱量設定部１３６０に入力すると、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）が負荷率計算部１３５８で計算される。排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収量ＨｈＥＸＨ（ｒｅｃ）（ｘＥＸＨ（ｔ））は、図１５（ａ）に示されているように、負荷率ｘＥＸＨ（ｔ）の関数である。従って、排熱回収量ＨｈＥＸＨ（ｒｅｃ）（ｘＥＸＨ（ｔ））は、排熱量入力部１３５６の回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））と機器特性モデル生成部１３６６の出力と機器諸元格納部１３５４からの機器諸元とを使用して、排熱回収量計算部１３６８で計算される。この際に、排熱回収量ＨｈＥＸＨ（ｒｅｃ）（ｘＥＸＨ（ｔ））という絶対値を求めるのに機器の定格が必要となるため、機器諸元格納部１３５４から機器定格を取得して計算を行う。排熱回収量ＨｈＥＸＨ（ｒｅｃ）（ｘＥＸＨ（ｔ））がわかると消費燃料（所要燃料）ＦＳＥＸＨ（ｔ）がわかるので、燃料供給設備（ＦＳＦ）１０４から供給を受けることが必要となる所要燃料が所要燃料計算部１３７４において計算される。
【００２９】
さらに、計算された所要燃料に基づいてコストをコスト計算部１３７６で計算する。コストを計算する際に、外部から電力料金単価、ガス料金単価、及びＣＯ２排出原単価を取得する。コストの計算結果は、計算結果格納部１３７８に格納される。
【００３０】
（３）  電動冷凍機及びガス吸収式冷温水機
図３は、電動冷凍機（ＥＲＵ）１０８の機能ブロック図である。  
電動冷凍機（ＥＲＵ）１０８は、機器諸元入力部４０２、機器諸元格納部１４０４、負荷率計算部１４０６、熱量設定部１４０８、機器特性データ入力部１４１０、機器特性データ格納部１４１２、機器特性モデル生成部１４１４、冷却水温入力部１４１６、冷却水温格納部１４１８、ＣＯＰ計算部１４２０、所要燃料計算部１４２２、コスト計算部１４２４、計算結果格納部１４２６、補機類稼働判定部１４２８を含む。
【００３１】
電動冷凍機（ＥＲＵ）１０８の動作は冷却水温ＴＥＲＵ（ｔ）に大きく影響を受ける。まず、機器特性データ入力部１４１０にＴ＝１７℃、２０℃、２４℃、２８℃、３２℃についてのＣＯＰ値を入力し、機器特性データ格納部１４１２に入力されたＣＯＰ値を格納しておく。これと並行して、冷却水温入力部１４１６に冷却水温ＴＥＲＵ（ｔ）を入力し、冷却水温格納部１４１８に冷却水温ＴＥＲＵ（ｔ）を格納しておく。機器特性モデル生成部１４１４では、上記のＣＯＰ値を補間することにより、機器特性ＣＯＰ（ｉ，ＴＥＲＵ（ｔ），ｘＥＲＵ（ｔ））を生成する（ｉ＝１，…，ｍは、熱源機器のインデックスを表すものとする）。生成された機器特性ＣＯＰ（ｉ，ＴＥＲＵ（ｔ），ｘＥＲＵ（ｔ））は、ＣＯＰ計算部１４２０に送られ、機器特性ＣＯＰ（ｉ，ＴＥＲＵ（ｔ），ｘＥＲＵ（ｔ））からある特定の負荷率ｘＥＲＵ（ｔ０）と冷却水温ＴＥＲＵ（ｔ０）についてのＣＯＰ値が計算される。ＣＯＰ値がわかれば、燃料消費特性がわかるため、ある特定の熱需要に対する所要燃料を所要燃料計算部１４２２で計算することができる。所要燃料がわかれば、コスト計算部１４２４でコストを計算することができる。
【００３２】
図４は、ガス吸収式冷温水機（ＧＡＢ）１１４の機能ブロック図である。
【００３３】
ガス吸収式冷温水機（ＧＡＢ）１１４は、機器諸元入力部１５０２、機器諸元格納部１５０４、負荷率計算部１５０６、熱量設定部１５０８、機器特性データ入力部１５１０、機器特性データ格納部１５１２、機器特性モデル生成部１５１４、冷却水温入力部１５１６、冷却水温格納部１５１８、ＣＯＰ計算部１５２０、所要燃料計算部１５２２、コスト計算部１５２４、計算結果格納部１５２６、補機類稼働判定部１５２８を含む。
【００３４】
ガス吸収式冷温水機（ＧＡＢ）１１４は、ガスを燃焼させて内部にある冷媒を濃縮することにより、冷房用の冷熱を生成する熱源機器である。ガス吸収式冷温水機（ＧＡＢ）１１４の動作は冷却水温ＴＧＡＢ（ｔ）に大きく影響を受ける。まず、機器特性データ入力部１５１０にＴ＝１７℃、２０℃、２４℃、２８℃、３２℃についてのＣＯＰ値を入力し、機器特性データ格納部１５１２に入力されたＣＯＰ値を格納しておく。これと並行して、冷却水温入力部１５１６に冷却水温ＴＧＡＢ（ｔ）を入力し、冷却水温格納部１５１８に冷却水温ＴＧＡＢ（ｔ）を格納しておく。機器特性モデル生成部１５１４では、入力された上記のＣＯＰ値を補間することにより、機器特性ＣＯＰ（ｉ，ＴＧＡＢ（ｔ），ｘＧＡＢ（ｔ））を生成する。生成された機器特性ＣＯＰ（ｉ，ＴＧＡＢ（ｔ），ｘＧＡＢ（ｔ））は、ＣＯＰ計算部１５２０に送られ、機器特性ＣＯＰ（ｉ，ＴＧＡＢ（ｔ），ｘＧＡＢ（ｔ））からある特定の負荷率ｘＧＡＢ（ｔ０）と冷却水温ＴＧＡＢ（ｔ０）についてのＣＯＰ値が計算される。ＣＯＰ値がわかれば、燃料消費特性がわかるため、ある特定の熱需要に対する所要燃料を所要燃料計算部１５２２で計算することができる。所要燃料がわかれば、コスト計算部１５２４でコストを計算することができる。
【００３５】
（４）  運転計画最適化システムの全体構成
図５は、本願発明に従った熱源機器ネットワークの運転計画を最適化するためのシステムの概念図である。
【００３６】
本発明の最適化システム２００は、電力需要予測ブロック２０２、熱需要予測ブロック２０４、外気冷房（ナイトパージ）ブロック２０６、比エンタルピ予測ブロック２０８、冷却水温予測ブロック２１０、機器特性モデル化ブロック２１２、需要計算ブロック２１４、熱源機器運転計画最適化ブロック２１６、実績データ入力ブロック２１８、計画・予測データ入力ブロック２２０、熱源機器・補機類諸元入力ブロック２２２、熱源機器特性データ入力ブロック２２４、運転計画作成諸元入力ブロック２２６、及び運転計画出力ブロック２２８を含む。
【００３７】
実績データ入力ブロック２１８及び計画・予測データ入力ブロック２２０の入力データは、各種予測ブロック２０２乃至２０８において使用される。実績データ入力ブロック２１８では、カレンダー情報、建物内外の気温θｉｎ及びθｏｕｔの実績値（ｔ−１以前のデータ）、建物内外の湿度ｈｕｍｉｎ及びｈｕｍｏｕｔの実績値、各熱源機器の冷却水温Ｔの実績値、建物内に入る人の数、イベント情報、熱源機器別の需要の実績値、及び換気熱量ＶＱの実績値が入力される。計画・予測データ入力ブロック２２０では、カレンダー情報、建物内外の気温θｉｎ（ｔ）及びθｏｕｔ（ｔ）の予測値、建物内外の湿度ｈｕｍｉｎ（ｔ）及びｈｕｍｏｕｔ（ｔ）の予測値、建物に入る人の数ＮＰ（ｔ）の予測値、イベント情報、及び換気熱量ＶＱ（ｔ）の予測値が入力される。実績データ入力ブロック２１８の入力データは、電力需要予測ブロック２０２、熱需要予測ブロック２０４、及び冷却水温予測ブロック２１０に供給される。計画・予測データ入力ブロック２２０の入力データは、電力需要予測ブロック２０２、熱需要予測ブロック２０４、外気冷房（ナイトパージ）ブロック２０６、比エンタルピ予測ブロック２０８、及び冷却水温予測ブロック２１０に供給される。
【００３８】
機器特性データとしての成績係数（ＣＯＰ）は、機器特性モデル化ブロック２１２に供給される。機器特性モデル化ブロック２１２は、発電設備（ＥＧＥ）１１０の負荷率ｘＥＧＥ（ｔ）の関数として、総排熱量ＨｈＥＧ（ｇｒｏｓｓ）（ｘＥＧＥ（ｔ））を決定し、さらに、負荷率ｘＥＧＥ（ｔ）の関数として、回収可能な排熱量（すなわち、排熱投入型吸収冷温水機（ＥＸＨ）１１２に供給される排熱量）ＨｈＥＧ（ｘＥＧＥ（ｔ））を決定する。機器特性モデル化ブロック２１２は、決定された回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））を使用して、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）の関数として、排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収量ＨｈＥＸＨ（ｒｅｃ）（ｘＥＸＨ（ｔ））を表す。さらに、機器特性モデル化ブロック２１２は、排熱投入型吸収冷温水機（ＥＸＨ）１１２に対して割り当てられた熱需要ＨｃＥＸＨ（ｒｅｃ）（ｔ）と排熱回収量ＨｈＥＸＨ（ｒｅｃ）（ｘＥＸＨ（ｔ））とから、排熱投入型吸収冷温水機（ＥＸＨ）１１２のガス焚き分負荷率をｘ（ｇｄ）ＥＸＨ（ｔ）を求め、ガス焚き分負荷率をｘ（ｇｄ）ＥＸＨ（ｔ）の関数として、排熱投入型吸収冷温水機（ＥＸＨ）１１２のガス焚きＣＯＰ特性をＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））として決定する。機器特性モデル化ブロック２１２は、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）の関数として、排熱投入型吸収冷温水機（ＥＸＨ）１１２の排熱ＣＯＰ特性をＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））として決定する。また、器特性モデル化ブロック２１２は、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の負荷率ｘＥＲＵ（ｔ）及びのｘＧＡＢ（ｔ）、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の冷却水温ＴＥＲＵ（ｔ）及びＴＧＡＢ（ｔ）の関数として、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の熱源機器モデルＣＯＰ（ｘ（ｉ，ｔ），Ｔ（ｉ，ｔ））を計算する（ｉ＝１，…，ｍは、熱源機器のインデックスを表すものとする）。
【００３９】
（５）  電力需要予測ブロック及び熱需要予測ブロック
電力需要予測ブロック２０２での電力需要予測は、周知の方法で行われる。
【００４０】
図６は、需要予測部のハードウェア構成の機能ブロック図を示している。
【００４１】
需要予測部３００は、ＣＰＵ３０２、メモリ３０４、記憶装置３０６、インタフェース３０８、入力装置３１０、及び出力装置３１２を含む。ＣＰＵ３０２は、記憶装置３０６に記憶されているプログラムをメモリ３０４に読み出して実行することにより各種の機能を実現する。メモリ３０４は、例えば、ランダム・アクセス・メモリ（ＲＡＭ：ｒａｎｄｏｍ  ａｃｃｅｓｓ  ｍｅｍｏｒｙ）及び／又は読取り専用メモリ（ＲＯＭ：ｒｅａｄ  ｏｎｌｙ  ｍｅｍｏｒｙ）を含む。記憶装置３０６は、テープ・ドライブ、フロッピー（登録商標）・ドライブ、ハード・ディスク・ドライブ又はコンパクト・ディスク・ドライブを含む。インタフェース３０８は、センサ３１４と接続する、例えば、シリアルケーブルや通信ケーブル等を接続するためのインタフェースである。センサ３１４は、外気温θｏｕｔを測定するセンサであってもよい。需要予測部３００は、インタフェース３０８を介してセンサ３１４から外気温θｏｕｔの測定値を取得することができる。
【００４２】
入力装置３１０は、データの入力を受け付ける、例えば、キーボードやマウス、トラックボール、タッチパネル、マイクロフォン等を含む。出力装置１０６は、データを出力する、例えば、ディスプレイやプリンタ、スピーカ等を含む。
【００４３】
図７は、需要予測部のソフトウェア構成を示す図である。
【００４４】
需要予測部４００は、データベース４０２、実績データ（実測値）期間入力部４０４、実績データ取得部４０６、及びパラメータ推計部４０８を含む。
【００４５】
データベース４０２は、外気温θｏｕｔの実績値を含む、熱需要予測に必要なデータを記憶する。データベース４０２には、年月日、曜日、時刻、休日等のカレンダー情報が記憶されており、そのカレンダー情報の中の日付及びある特定の単位時間ｔに関連させた形で、その単位時間ｔ以前の単位時間、例えば、単位時間ｔ−１及びｔ−２における、外気温θｏｕｔの実績値、建物内空気の温度（内気温）θｉｎの実績値、建物外の湿度（外湿度）ｈｕｍｏｕｔの実績値、建物内の湿度（内湿度）ｈｕｍｉｎの実績値、熱需要Ｈｄの実績値、換気熱量ＶＱの実績値、建物内の目標設定温度、建物に入る人の数ＮＰ、及びイベントの有無ＩＶが記憶される。
【００４６】
換気熱量ＶＱ（ｔ）は、単位時間ｔにおいて建物を換気したことにより変化した熱量として算出される。比エンタルピは温度と湿度とから算出することができる。また、空気密度は温度と大気圧とから算出することができる。
【００４７】
単位時間ｔにおいて建物外から取り入れられた熱量Ｄｉｎ（ｔ）は、建物外の空気の比エンタルピをｈｏｕｔ（ｔ）［ｋＪ／ｋｇ］、建物外の空気の空気密度をρｏｕｔ［ｋｇ／ｍ３］、換気量をＶＶ（ｔ）［ｍ３］とすると、
Ｄｉｎ（ｔ）＝ｈｏｕｔ（ｔ）×ρｏｕｔ×ＶＶ（ｔ）            （６）
により求めることができる。
【００４８】
また、単位時間ｔにおいて建物外に排出された熱量Ｄｏｕｔ（ｔ）は、建物内の空気の比エンタルピをｈｉｎ（ｔ）［ｋＪ／ｋｇ］、建物内の空気の空気密度をρｉｎ［ｋｇ／ｍ３］、換気量をＶＶ（ｔ）［ｍ３］とすると、
Ｄｏｕｔ（ｔ）＝ｈｉｎ（ｔ）×ρｉｎ×ＶＶ（ｔ）                  （７）
により求めることができる。
【００４９】
従って、換気熱量ＶＱ（ｔ）は、
ＶＱ（ｔ）＝Ｄｉｎ（ｔ）−Ｄｏｕｔ（ｔ）                         （８）
により求めることができる。
【００５０】
本発明では、ＣＰＵ３０２は、計画・予測データ入力ブロック２２０から取得した単位時間ｔにおける建物外の気温θｏｕｔ（ｔ）の予測値及び建物外の湿度ｈｕｍｏｕｔ（ｔ）に基づいて建物外の空気の比エンタルピｈｏｕｔ（ｔ）の予測値を算出し、計画・予測データ入力ブロック２２０から取得した単位時間ｔにおける建物内の気温θｉｎ（ｔ）の予測値及び建物内の湿度ｈｕｍｉｎ（ｔ）に基づいて建物内の空気の比エンタルピｈｉｎ（ｔ）の予測値を算出する。ＣＰＵ３０２は、また、単位時間ｔにおいて換気が行われた場合にはＶＶ（ｔ）の値を取得し、換気が行われていない場合にはＶＶ（ｔ）を０とする。ＣＰＵ３０２は、所定の大気圧と、建物内外の気温θｉｎ（ｔ）及びθｏｕｔ（ｔ）に基づいて建物内外の空気の空気密度をρｉｎ及びρｏｕｔを算出する。ＣＰＵ３０２は、上記式に建物外の空気の比エンタルピをｈｏｕｔ（ｔ）、建物外の空気の空気密度をρｏｕｔ及び換気量ＶＶ（ｔ）を適用して、単位時間ｔにおいて建物外から取り入れられた熱量Ｄｉｎ（ｔ）を算出し、上記式に建物内の空気の比エンタルピをｈｉｎ（ｔ）、建物内の空気の空気密度をρｉｎ及び換気量ＶＶ（ｔ）を適用して、単位時間ｔにおいて建物外に排出された熱量Ｄｏｕｔ（ｔ）を算出し、上記式により換気熱量ＶＱ（ｔ）を算出する。
【００５１】
実績データ取得部４０６は、予測を行う単位時間を実績データ期間入力部４０４から取得し、パラメータ推計に用いる実績データのセットをデータベース４０２から取得する。パラメータ推計部４０８は、熱需要の予測モデルに用いるパラメータのセットを推計する。熱需要の予測値Ｈｄ＊（ｔ）の予測モデルは、式
Ｈｄ＊（ｔ）＝ｋ０＋ｋ１・θｏｕｔ（ｔ）＋ｋ２・ｈｕｍｏｕｔ（ｔ）＋ｋ３・ＮＰ（ｔ）＋ＶＱ（ｔ）      （９）
によって表される。ここで、ｋ０は定数であり、ｋ１、ｋ２、及びｋ３は回帰係数である。パラメータ推計部４０８は、変量のセットΩ１＝［θｏｕｔ（ｔ），ｈｕｍｏｕｔ（ｔ），ＮＰ（ｔ），ＶＱ（ｔ）］の多数のサンプルを使用して、定数ｋ０及び回帰係数ｋ１、ｋ２、及びｋ３の推計を行い、推計処理は、一般的な回帰分析の手法により行われる。代替的に、推計処理は、ＬＭＳアルゴリズムを用いて適応的に行われてもよい。
【００５２】
パラメータ推計部４０８は、推計したパラメータのセットＫ＝［ｋ０，ｋ１，ｋ２，ｋ３］をデータベース４０２に格納する。ＣＰＵ３０２は、単位時間ｔにおける変量のセットΩ１＝［θｏｕｔ（ｔ），ｈｕｍｏｕｔ（ｔ），ＮＰ（ｔ），ＶＱ（ｔ）］に格納したパラメータのセットＫ＝［ｋ０，ｋ１，ｋ２，ｋ３］を適用して、熱需要の予測値Ｈｄ＊（ｔ）を計算する。
【００５３】
代替的に、熱需要予測は、ニューラルネットワーク（ＮＮ）、カルマンフィルタ等のモデル、時系列解析、ファジー理論等の手法を用いて行われても良い。
【００５４】
建物に入る人の数ＮＰが増えれば熱需要も増加するが、小売店や飲食店等の店舗においては、経験的に、曜日や時間帯によって来客数が変動する。また、建物のエリア内においてイベントが開催された場合には通常とは異なる人の数ＮＰとなるが、イベントの有無を表すＩＶも考慮される。上記の予測モデルによれば、曜日や時間帯、イベントの有無を考慮することで、熱需要の予測値の精度が向上する。
【００５５】
（６）  外気冷房（ナイトパージ）ブロック
図８は、外気冷房（ナイトパージ）ブロック２０６の動作を示す概念図である。
【００５６】
図９は、比エンタルピ予測ブロック２０８及び外気冷房（ナイトパージ）ブロック２０６によって行われる排熱量予測の概念図である。
【００５７】
図５及び図９を参照すると、比エンタルピ予測ブロック２０８は、実績データ入力ブロック２１８からカレンダー情報、建物内外の気温θｉｎ及びθｏｕｔの実測値（ｔ−１以前のデータ）、及び建物内外の湿度ｈｕｍｉｎ及びｈｕｍｏｕｔの実測値を取得し、計画・予測データ入力ブロック２２０からカレンダー情報、建物内外の気温θｉｎ（ｔ）及びθｏｕｔ（ｔ）の予測値、及び建物内外の湿度ｈｕｍｉｎ（ｔ）及びｈｕｍｏｕｔ（ｔ）の予測値を取得する。上記のように、比エンタルピ予測ブロック２０８は、単位時間ｔにおける建物外の気温θｏｕｔ（ｔ）の予測値及び建物外の絶対湿度ｈｕｍｏｕｔ（ｔ）に基づいて建物外の空気の比エンタルピをｈｏｕｔ（ｔ）の予測値を算出し、単位時間ｔにおける建物内の気温θｉｎ（ｔ）の予測値及び建物内の絶対湿度ｈｕｍｉｎ（ｔ）に基づいて建物内の空気の比エンタルピをｈｉｎ（ｔ）の予測値を算出する。比エンタルピ予測ブロック２０８は、予測した建物内外の空気の比エンタルピｈｉｎ（ｔ）及びｈｏｕｔ（ｔ）を外気冷房（ナイトパージ）ブロック２０６に供給する。
【００５８】
図１０は、外気冷房（ナイトパージ）制御装置を示す機能ブロック図である。
【００５９】
外気冷房（ナイトパージ）ブロック２０６は、実績データ入力ブロック２１８から取得した建物外の気温θｏｕｔの実績値及び建物外の湿度ｈｕｍｏｕｔに基づいて建物外の空気のエンタルピの実績値を算出し、実績データ入力ブロック２１８から取得した建物内の気温θｉｎの実績値及び建物内の湿度ｈｕｍｉｎに基づいて建物内の空気のエンタルピの実績値を算出する。外気冷房（ナイトパージ）ブロック２０６は、比エンタルピ予測ブロック２０８から建物外の空気のエンタルピの予測値及び建物外の空気のエンタルピの予測値を取得する。
【００６０】
ナイトパージ時間帯判定装置７０４は、比エンタルピ予測ブロック２０８から取得した建物外の空気の比エンタルピをｈｏｕｔ（ｔ）の予測値及び建物内の空気の比エンタルピをｈｉｎ（ｔ）の予測値に基づいて各単位時間ｔにおけるエンタルピ差予測値を算出し、ナイトパージ設定装置７０２で設定された基準エンタルピ差以上のエンタルピ差があり、最もエンタルピ差の大きい単位時間を選択する。ナイトパージ実行判定装置７０６は、ナイトパージ設定装置７０２で設定されたナイトパージの基準エンタルピ差に基づいて、ナイトパージ時間帯判定装置７０４によって求められたナイトパージ実行時間帯の予測エンタルピ差が最大となる単位時間から、又は、基準エンタルピ差を越える予測単位時間から、建物内の空気のエンタルピの実績値から建物外の空気のエンタルピの実績値を減じた差が基準エンタルピ差を越えているか否かを判定する。タイマ装置７０８は、ナイトパージ実行判定装置７０６が判定した、建物内外のエンタルピ差の実績値の間の差が基準エンタルピ差を越えている単位時間を検出し、検出した単位時間をナイトパージ実行判定装置７０６に送る。ナイトパージ実行判定装置７０６は、タイマ装置７０８が検出した単位時間からナイトパージ設定装置７０２で設定されたナイトパージ継続時間が経過した場合にナイトパージの終了の判定を行う。
【００６１】
外気冷房（ナイトパージ）ブロック２０６の排熱量予測ブロック６０８は、建物外の空気の比エンタルピｈｏｕｔ（ｔ）の予測値、建物外の空気の空気密度ρｏｕｔ、及び予定換気量ＶＶ（ｔ）に基づいて、単位時間ｔにおいて建物外から取り入れられた熱量Ｄｉｎ（ｔ）の予測値を算出し、建物内の空気の比エンタルピｈｉｎ（ｔ）の予測値、建物内の空気の空気密度ρｉｎ、及び予定換気量ＶＶ（ｔ）に基づいて、単位時間ｔにおいて建物外に排出された熱量Ｄｏｕｔ（ｔ）の予測値を算出する。排熱量予測ブロック６０８は、建物外から取り入れられた熱量Ｄｉｎ（ｔ）の予測値から建物外に排出された熱量Ｄｏｕｔ（ｔ）の予測値を減じて、換気熱量ＶＱ（ｔ）の予測値を算出し、これに基づいて、排熱量の予測値を求める。
【００６２】
図１１は、建物外のθｏｕｔ（ｔ）の予測値、湿度ｈｕｍｏｕｔ（ｔ）の予測値、空気密度ρｏｕｔ、及び比エンタルピｈｏｕｔ（ｔ）、建物内の気温θｉｎ（ｔ）の予測値、湿度ｈｕｍｉｎ（ｔ）の予測値、空気密度ρｉｎ、及び比エンタルピｈｉｎ（ｔ）、換気量ＶＶ（ｔ）、及び排熱量の関係を示す図である。
【００６３】
排熱量の予測値を求めると、熱需要の予測値Ｈｄ＊（ｔ）のうちどれだけを外気冷房で賄えるかを予測することができる。上記のように、電力需要予測ブロック２０２、熱需要予測ブロック２０４、外気冷房（ナイトパージ）ブロック２０６、及び比エンタルピ予測ブロック２０８を取り込むことにより、より正確な電力需要Ｐｄ（ｔ）及び熱需要Ｈｄ（ｔ）を計算することができる。
【００６４】
（７）  冷却水温予測ブロック
図１２は、各熱源機器の冷却水温予測の概念図である。
【００６５】
本発明では、各熱源機器の運転効率を上げるために、各熱源機器の冷却水温Ｔ（ｉ，ｔ）の予測精度の向上を図る。冷熱負荷が大きく、負荷変動の大きいビルやデパート等の空調には複数の熱源機器が設置されており、コストが小さくなるなどの目的に沿って運転計画を作成して、これに基づいて熱源機器の運転を行っている。熱源機器の冷凍能力に関しては、過去の実績を参考にして運転計画を作成しているのが現状であるが、実際の冷凍能力は冷却水温に大きく依存している。運転計画作成のためには冷却水温が重要であるが、運転計画作成の際に冷却水温を予測する先行技術は、現段階では存在しない。
【００６６】
熱源を通過した冷却水は、一般的に、冷却塔へ送水され、そこで冷却される。従って、冷却水温は、熱源の放熱量（すなわち、空調需要）と外気温の影響を大きく受けることが予想される。そこで、本発明では、熱需要予測結果と気象庁等から公表される予測気温値を説明変数として使用し、冷却水温を予測する。
【００６７】
冷却水温予測ブロック２１０は、実績データ入力ブロック２１８からカレンダー情報、建物外の気温θｏｕｔの実績値、各熱源機器（ｉ＝１，…，ｍ）の冷却水温Ｔ（ｉ）の実績値、及び熱需要Ｈｄの実績値を取得し、計画・予測データ入力ブロック２２０からカレンダー情報及び建物外の気温θｏｕｔの予測値θ＊ｏｕｔ（ｔ）を取得する。冷却水温予測ブロック２１０は、また、熱需要予測ブロック２０４から熱需要Ｈｄの予測値Ｈ＊ｄ（ｔ）を取得する。ｌ０を定数、ｌ１、ｌ２、ｌ３、及びｌ４を回帰係数として、冷却水温Ｔ（ｉ，ｔ）の予測モデルは、
Ｔ＊（ｉ，ｔ）＝ｌ０＋ｌ１・θ＊ｏｕｔ（ｔ）＋ｌ２・θｏｕｔ（ｔ−１）＋ｌ３・Ｔ（ｉ，ｔ−１）＋ｌ４・（Ｈｄ（ｉ，ｔ−１）−Ｈｄ（ｉ，ｔ−２））  （ｉ＝１，…、ｍ）    （１０）
で表される。
【００６８】
冷却水温予測ブロック２１０は、変量のセットΩ２＝［θ＊ｏｕｔ（ｔ），θｏｕｔ（ｔ−１），Ｔ（ｉ，ｔ−１），Ｈｄ（ｉ，ｔ−１），Ｈｄ（ｉ，ｔ−２）］の多数のサンプルを使用して、定数ｌ０、及び回帰係数ｌ１、ｌ２、ｌ３、及びｌ４の推計を行い、推計処理は、一般的な回帰分析の手法により行われる。代替的に、推計処理は、ＬＭＳアルゴリズムを用いて適応的に行われてもよい。
【００６９】
冷却水温予測ブロック２１０は、推計したパラメータのセットＬ＝［ｌ０，ｌ１，ｌ２，ｌ３，ｌ４］を図示しないデータベースに格納する。冷却水温予測ブロック２１０は、変量のセットΩ２＝［θ＊ｏｕｔ（ｔ），θｏｕｔ（ｔ−１），Ｔ（ｉ，ｔ−１），Ｈｄ（ｉ，ｔ−１），Ｈｄ（ｉ，ｔ−２）］に格納したパラメータのセットＬ＝［ｌ０，ｌ１，ｌ２，ｌ３，ｌ４］を適用して、冷却水温Ｔ（ｉ）の予測値Ｔ＊（ｉ，ｔ）を計算する。
【００７０】
図１３は、式（１０）で示される予測モデルを用いて予測した冷却水温Ｔ（ｉ）の予測値Ｔ＊（ｉ，ｔ）と実績値との比較を示す図である。
【００７１】
冷却水温予測は、ニューラルネットワーク（ＮＮ）、カルマンフィルタ等のモデル、時系列解析、ファジー理論等の手法を用いて行われても良い。
【００７２】
冷却水の水温を予測するのは、各熱源機器において冷却水は設備冷却後に元の場所に戻されるため、各熱源機器の負荷率ｘ（ｉ，ｔ）及び効率（すなわち、燃料及び電力の消費特性）と相関を持っているためであり、冷却水の水温を予測することは、各熱源機器の機器モデルを決定する上で必要だからである。予測モデル（１０）に熱需要Ｈｄを取り込むのは、熱需要Ｈｄの値が冷却水の水温予測の過程でフィードバックとして作用するからである。式（１０）の最後の項、すなわち、（Ｈｄ（ｉ，ｔ−１）−Ｈｄ（ｉ，ｔ−２））は熱需要Ｈｄの増加率（勾配）であるが、この項は、温度変化の値では補えない部分を補足する補正項として作用する。
【００７３】
上記のように、熱源機器ネットワークの運転計画を最適化するためのシステムに、電力需要予測ブロック２０２、熱需要予測ブロック２０４、外気冷房（ナイトパージ）ブロック２０６、及び比エンタルピ予測ブロック２０８を取り込むことにより、より正確な電力需要及び熱需要の計算が可能となる。需要計算ブロック２１４は、上記のブロックからの情報に基づいて、単位時間ｔにおける電力需要Ｐｄ（ｔ）及び熱需要Ｈｄ（ｔ）を計算する。
【００７４】
（８）  機器特性モデル化ブロック（コジェネレーション・システムの場合）
一般的に、空調負荷が大きく、負荷変動の大きいビルやデパート等の空調には複数の熱源機器が設置されており、コストが小さくなるなどの目的に沿って運転計画を作成して、これに基づいて多種多様な熱源機器からなる熱源機器ネットワークの運転を行っている。
【００７５】
ここで、発電設備と排熱利用吸収冷温水機とからなるコジェネレーション・システムの運転計画作成の際には、実測もしくはシュミュレーションによって得られた特性カーブを機器メーカーから入手し、これを運転計画システムに入力することにより運転計画を作成しているのが現状である。
【００７６】
本発明は、メーカーから提供される特性データを用いて、特性関数を推計してモデル化し、任意の負荷率における特性値を容易に取得できるようにすることを目的の１つとしている。
【００７７】
ここで、コジェネレーション・システムの場合には、発電を行なうためのガスエンジン、ガスタービン、燃料電池等の原動機と、発電時に発生する熱を利用する熱交換器、熱で駆動する吸収式冷温水器等の排熱利用機器から構成され、発電出力と熱の両方を有効に利用しなければ、本来の効率を発揮できない。本発明においても、排熱投入型吸収冷温水機は発電機からの排熱を利用しているため、発電設備（ＥＧＥ）１１０の機器特性と排熱投入型吸収冷温水機（ＥＸＨ）１１２の機器特性をモデル化する必要がある。
【００７８】
本発明においては、メーカーから提供される特性カーブを、負荷率を変数として、発電設備（ＥＧＥ）１１０における排熱量と排熱投入型吸収冷温水機（ＥＸＨ）１１２における排熱回収量及び成績係数（Ｃｏｅｆｆｉｃｉｅｎｔ  ｏｆ  Ｐｅｒｆｏｒｍａｎｃｅ（ＣＯＰ））を回帰して、係数を推計する。すなわち、負荷率がわかれば、精度の高い発電機排熱量を得た後、排熱回収量や成績係数を容易に得ることが可能となる。
【００７９】
（i）  発電機の機器特性モデルの作成及びコスト計算方法
図１４は、発電設備（ＥＧＥ）１１０の機器特性モデルの作成及びコスト計算方法を示すブロック図である。
【００８０】
発電設備（ＥＧＥ）１１０の機器特性モデル化ブロック２１２は、発電設備（ＥＧＥ）１１０の機器特性データ格納部１３１２から発電設備（ＥＧＥ）１１０の機器特性を取得する。機器特性データとして、発電設備（ＥＧＥ）１１０の負荷率ｘＥＧＥ（ｔ）が、ｘＥＧＥ（ｔ）＝５０［％］，７５［％］，１００［％］の各値をとる場合について、それぞれ、入熱量［ＭＪ／ｈ］、軸出力［ＭＪ／ｈ］、冷却水熱量［ＭＪ／ｈ］、排気ガス熱量［ＭＪ／ｈ］、インタークーラー熱量［ＭＪ／ｈ］、放熱量［ＭＪ／ｈ］の各値が与えられている。ここで、入熱量は、燃料供給設備（ＦＳＦ）１０４から発電設備（ＥＧＥ）１１０に供給される燃料ＦＳＥＧＥ（ｔ）を燃焼させた際に発生する熱量を示しており、軸出力は、発電設備（ＥＧＥ）１１０に供給される燃料ＦＳＥＧＥ（ｔ）を燃焼させた際に発生する熱量のうち、発電機のタービンの軸を回すために当該熱量がタービンに対してする仕事量を示しており、冷却水熱量は、冷却水に放熱される熱量を示しており、排気ガス熱量は、排気ガス中に含まれる熱量を示しており、インタークーラー熱量は、冷却フィンから放熱される熱量を示しており、放熱量は、発電設備の壁面等から放熱される熱量を示している。
【００８１】
ここで、発電設備（ＥＧＥ）１１０の総排熱量Ｈ（ｇｒｏｓｓ）ｈＥＧ（ｘＥＧＥ（ｔ））は、冷却水熱量［ＭＪ／ｈ］、排気ガス熱量［ＭＪ／ｈ］、インタークーラー熱量［ＭＪ／ｈ］、及び放熱量［ＭＪ／ｈ］の合計をとることで求められ、
0006641854.tif 000002
となる。ここで、ｋ＝１，…，４は、冷却水熱量［ＭＪ／ｈ］（ｋ＝１）、排気ガス熱量［ＭＪ／ｈ］（ｋ＝２）、インタークーラー熱量［ＭＪ／ｈ］（ｋ＝３）、及び放熱量［ＭＪ／ｈ］（ｋ＝４）のそれぞれに対応しており、ｃｋは定数であり、ｂｋは負荷率ｘＥＧＥ（ｔ）に乗じる係数である。ｘ（ｔ｜ＥＧＥ）はｘＥＧＥ（ｔ）を表している。このうち、回収可能な排熱は、冷却水熱量［ＭＪ／ｈ］（ｋ＝１）及び排気ガス熱量［ＭＪ／ｈ］（ｋ＝２）であるので、回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））は、
0006641854.tif 000003
となる。
【００８２】
発電設備（ＥＧＥ）１１０の機器特性モデル生成部１３１４は、負荷率ｘＥＧＥ（ｔ）＝５０［％］，７５［％］，１００［％］についての冷却水熱量［ＭＪ／ｈ］（ｋ＝１）、排気ガス熱量［ＭＪ／ｈ］（ｋ＝２）、インタークーラー熱量［ＭＪ／ｈ］（ｋ＝３）、及び放熱量［ＭＪ／ｈ］（ｋ＝４）の値を内挿補間して、任意の負荷率ｘＥＧＥ（ｔ）についてパラメータのセット［ｃｋ，ｂｋ］（ｋ＝１，…，４）を推計する。代替的に、発電設備（ＥＧＥ）１１０の機器特性モデル生成部１３１４は、離散的な負荷率ｘＥＧＥ（ｔ）の値に対する冷却水熱量［ＭＪ／ｈ］（ｋ＝１）、排気ガス熱量［ＭＪ／ｈ］（ｋ＝２）、インタークーラー熱量［ＭＪ／ｈ］（ｋ＝３）、及び放熱量［ＭＪ／ｈ］（ｋ＝４）の値の多数のサンプルを使用して重回帰分析を行い、パラメータのセット［ｃｋ，ｂｋ］（ｋ＝１，…，４）を推計しても良い。機器特性モデル生成部１３１４は、推計したパラメータのセット［ｃｋ，ｂｋ］を熱源機器運転計画最適化ブロック２１６に送り、熱源機器運転計画最適化ブロック２１６は、熱源機器モデル記憶ブロック２１６ａに受け取ったパラメータのセット［ｃｋ，ｂｋ］を格納する。
【００８３】
発電設備（ＥＧＥ）１１０の負荷率計算部１３０６は、機器諸元格納部１３０４に格納されている機器定格を取得し、電力需要をＰｄ（ｔ）のうち発電設備（ＥＧＥ）１１０に割り当てられている発電電力ＰＥＧ（ｔ）の値を発電量設定部１３０８から取得し、負荷率ｘＥＧＥ（ｔ）を計算する。計算された負荷率ｘＥＧＥ（ｔ）は、排熱量計算部１３１６と所要燃料計算部１３１８に送られる。所要燃料計算部１３１８では、機器特性モデル生成部１３１４からのパラメータのセット［ｃｋ，ｂｋ］と、負荷率ｘＥＧＥ（ｔ）を使用して、式（１１）に従って発電設備（ＥＧＥ）１１０の総排熱量Ｈ（ｇｒｏｓｓ）ｈＥＧ（ｘＥＧＥ（ｔ））を計算する。所要燃料計算部１３１８は、さらに、発電量ＰＥＧ（ｔ）と発電機の効率（すなわち、損失）ηＥＧＥ（ｘＥＧＥ（ｔ））と総排熱量Ｈ（ｇｒｏｓｓ）ｈＥＧ（ｘＥＧＥ（ｔ））とを用いて、所要燃料ＦＳＥＧＥ（ｔ）を、
ＦＳＥＧＥ（ｔ）＝ＰＥＧ（ｔ）／ηＥＧＥ（ｘＥＧＥ（ｔ））＋Ｈ（ｇｒｏｓｓ）ｈＥＧ（ｘＥＧＥ（ｔ））  （１３）
として求める。求めた所要燃料ＦＳＥＧＥ（ｔ）は、コスト等計算部１３２２に送られ、コスト等計算部１３２２は、発電設備（ＥＧＥ）１１０の総コストをＥＸＥＧＥ（ｔ）、燃料単価をｐｆｕｅ、電力単価をｐｐｏｗ、補機類駆動用電力をＰＥＧＰ（ｔ）とすると、総コストＥＸＥＧＥ（ｔ）を、
ＥＸＥＧＥ（ｔ）＝ＦＳＥＧＥ（ｔ）×ｐｆｕｅ＋ＰＥＧＰ（ｔ）×ｐｐｏｗ                （１４）
として求める。
【００８４】
（II）  排熱投入型吸収冷温水機の機器特性モデルの作成及びコスト計算方法
排熱回収型吸収冷温水機は、いわゆる、ジェネリンクと呼ばれるものであるが、発電設備のガスエンジンの排温水を補助熱源として有効利用することで、主燃料であるガスの使用量を定格時の冷凍出力比率で２０乃至４０パーセント削減可能な冷温水機である。従って、排熱のみを熱源として冷房負荷の冷熱を供給できるときは、排熱のみで動作し、排熱のみで熱需要を賄いきれなくなると、燃料供給設備から供給を受けたガスを使用して動作する。また、排熱を利用できない場合、例えば、発電設備のガスエンジンが停止している等の場合には、ガスのみで動作する。従って、排熱投入型吸収冷温水機のＣＯＰ特性は、（１）  排熱のみで動作する場合（排熱ＣＯＰ特性の場合）、（２）  排熱＋ガスで動作する場合、（３）  ガスのみで動作する場合（ガス焚きＣＯＰ特性の場合）で異なる特性となる。排熱＋ガスで動作する場合のＣＯＰ特性は、排熱ＣＯＰ特性及びガス焚きＣＯＰ特性から導き出すことができる。
【００８５】
図１５は、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率―排熱回収量特性及び負荷率―ＣＯＰ特性を示す図である。図１５（ａ）は、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率―排熱回収容量特性を示している。
【００８６】
排熱投入型吸収冷温水機（ＥＸＨ）１１２の排熱量入力部１３５６は、発電設備（ＥＧＥ）１１０の排熱量計算部１３１６で計算された回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））の計算結果を受け取る。機器諸元入力部１３５２に排熱投入型吸収冷温水機（ＥＸＨ）１１２の冷熱出力ＨｃＥＸＨ（ｔ）の定格値Ｈ（ｒａｔ）ｃＥＸＨ等の定格を表す諸元が入力され、熱需要Ｈｄ（ｔ）のうち排熱投入型吸収冷温水機（ＥＸＨ）１１２に割り当てられた熱需要である熱量設定値ＨｃＥＸＨ（ｒｅｑ）（ｔ）を熱量設定部１３６０に入力すると、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）が負荷率計算部１３５８で計算される。
【００８７】
図１５（ａ）に示されているように、排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））は、負荷率ｘＥＸＨ（ｔ）の関数である。図１５（ａ）において、０≦ｘＥＸＨ（ｔ）≦ｘ０（ｘ０≒５８％）の範囲内では、排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））は、負荷率ｘＥＸＨ（ｔ）に比例して増加するが、負荷率ｘＥＸＨ（ｔ）がｘ０を超えると、負荷率ｘＥＸＨ（ｔ）の３次関数で近似されるプロファイルに従って減少する。排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））の最大値は、おおよそ５８０［ｋＷ］である。従って、ｄ０を定数項、ｄ１を係数、ｅ０を定数項、ｅ１、ｅ２、ｅ３を係数とすると、排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））は、
0006641854.tif 000004
で表される。パラメータのセット［ｄ０，ｄ１］及び［ｅ０，ｅ１，ｅ２，ｅ３］は、離散的な負荷率ｘＥＸＨ（ｔ）の値に対するＨ（ｒｅｃ）ｈＥＸＨ（ｘＥＸＨ（ｔ））の値を内挿補間して、任意の負荷率ｘＥＸＨ（ｔ）についてのＨ（ｒｅｃ）ｈＥＸＨ（ｘＥＸＨ（ｔ））の値を求めることができるように推計することができる。代替的に、排熱投入型吸収冷温水機（ＥＸＨ）１１２の機器特性モデル生成部１３６６は、離散的な負荷率ｘＥＸＨ（ｔ）の値に対するＨ（ｒｅｃ）ｈＥＸＨ（ｘＥＸＨ（ｔ））の値の多数のサンプルを使用して重回帰分析を行い、パラメータのセット［ｄ０，ｄ１］及び［ｅ０，ｅ１，ｅ２，ｅ３］を推計しても良い。
【００８８】
発電設備（ＥＧＥ）１１０からの回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））及び排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））がわかると、排熱回収量計算部１３６８は、排熱投入型吸収冷温水機（ＥＸＨ）１１２での実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））を、
Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））＝ｍｉｎ｛ＨｈＥＧ（ｘＥＧＥ（ｔ）），Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））｝    （１６）
として決定する。すなわち、回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））と排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））のうちのいずれか小さなほうを排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））として決定する。これは、冷房需要Ｈ（ｒｅｃ）ｃＥＸＨ（ｔ）は大きい（この場合には、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）が大きくなる）が、発電設備（ＥＧＥ）１１０が発電する発電電力ＰＥＧ（ｔ）があまり大きくなく（電力需要をＰｄ（ｔ）があまり大きくなく）、排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））よりも回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））自体が小さい場合には、少ない排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））をそのまま使い、残りは、燃料供給設備（ＦＳＦ）１０４から供給された燃料ＦＳＥＸＨ（ｔ）を使用して、冷熱ＨｃＥＸＨ（ｔ）を生成せざるを得ない。逆に、発電設備（ＥＧＥ）１１０からの回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））が、排熱投入型吸収冷温水機（ＥＸＨ）１１２の排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））を超えている場合には、排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））の分だけ目いっぱい回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））を利用するということである。
【００８９】
図１５（ｂ）は、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率―ＣＯＰ特性を示している。排熱投入型吸収冷温水機（ＥＸＨ）１１２のＣＯＰ特性は、大きく分けて、（１）  排熱のみで動作する場合（排熱ＣＯＰ特性の場合）、（２）  排熱＋ガスで動作する場合、（３）  ガスのみで動作する場合（ガス焚きＣＯＰ特性の場合）で異なる特性となる。排熱＋ガスで動作する場合のＣＯＰ特性は、排熱ＣＯＰ特性及びガス焚きＣＯＰ特性から導き出すことができる。排熱ＣＯＰ特性は、投入排熱１単位当たりでどの程度熱需要を充足することができるかを示す尺度である。
【００９０】
上記のように、排熱投入型吸収冷温水機（ＥＸＨ）１１２では、熱量設定部１３６０に、熱需要Ｈｄ（ｔ）のうち排熱投入型吸収冷温水機（ＥＸＨ）１１２に割り当てられた熱需要である熱量設定値Ｈ（ｒｅｃ）ｃＥＸＨ（ｔ）を入力すると、負荷率計算部１３５８が、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）を計算する。ここで、排熱のみで動作する場合（上記の場合（１））の排熱回収ＣＯＰ特性をＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））、ｆ０を定数項、ｆ１、ｆ２、ｆ３を係数とすると、排熱回収ＣＯＰ特性は、
0006641854.tif 000005
で表すことができる。上記と同様に、パラメータのセット［ｆ０，ｆ１，ｆ２，ｆ３］は、離散的な負荷率ｘＥＸＨ（ｔ）の値に対するＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））の値を内挿補間して、任意の負荷率ｘＥＸＨ（ｔ）についてのＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））の値の値を求めることができるように推計することができる。代替的に、排熱投入型吸収冷温水機（ＥＸＨ）１１２の機器特性モデル生成部１３６６は、離散的な負荷率ｘＥＸＨ（ｔ）の値に対するＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））の値の多数のサンプルを使用して重回帰分析を行い、パラメータのセット［ｆ０，ｆ１，ｆ２，ｆ３］を推計しても良い。
【００９１】
排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”が決定されると、排熱回収量計算部１３６８は、排熱投入型吸収冷温水機（ＥＸＨ）１１２の発電機排熱回収分Ｈ（ｒｅｃ）ｃＥＧＥＸＨ（ｘＥＸＨ（ｔ））を、
Ｈ（ｒｅｃ）ｃＥＧＥＸＨ（ｘＥＸＨ（ｔ））＝ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））・Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＸＨ（ｔ））      （１８）
として求めることができる。
【００９２】
図１６は、排熱投入型吸収冷温水機（ＥＸＨ）１１２のコスト計算方法を示すブロック図である。
【００９３】
排熱投入型吸収冷温水機（ＥＸＨ）１１２は、排熱量入力部１３５６で、発電設備（ＥＧＥ）１１０の排熱量計算部１３１６で計算された回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））の計算結果を受け取る。一方、負荷率計算部１３５８は、熱量設定部１３６０から、排熱投入型吸収冷温水機（ＥＸＨ）１１２に割り当てられた熱需要である熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）を受け取り、機器諸元格納部１３５４から冷熱出力ＨｃＥＸＨ（ｔ）の定格値Ｈ（ｒａｔ）ｃＥＸＨを受け取り、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）を計算する。この負荷率ｘＥＸＨ（ｔ）から、式（１５）に従って、排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））を求めることができ、排熱回収量計算部１３６８は、式（１６）に従って、実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））を求めることができる。さらに、この負荷率ｘＥＸＨ（ｔ）から、図１５（ｂ）の（１）に示されている排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”を、式（１７）に従って求めることができ、実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））と排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”から、発電機排熱回収分Ｈ（ｒｅｃ）ｃＥＧＥＸＨ（ｘＥＸＨ（ｔ））を求めることができる。
【００９４】
熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）、実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））、排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”、定格値Ｈ（ｒａｔ）ｃＥＸＨがわかると、負荷率計算部１３５８は、ガス焚き分負荷率をｘ（ｇｄ）ＥＸＨ（ｔ）を、
0006641854.tif 000006
として求めることができる。式（１９）では、熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）から発電機排熱回収分Ｈ（ｒｅｃ）ｃＥＧＥＸＨ（ｘＥＸＨ（ｔ））、すなわち、ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））×Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））を減じ、その減じたものと０とを比較し、大きいものを選択する。そして、その選択したものを定格値Ｈ（ｒａｔ）ｃＥＸＨで割り、その割った値と１とを比較し、小さいほうを選択し、選択したものに１００を乗じてガス焚き分負荷率をｘ（ｇｄ）ＥＸＨ（ｔ）を求める。つまり、満たすべき需要である熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）のうち実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））で賄えない部分を求め、その部分と定格値Ｈ（ｒａｔ）ｃＥＸＨとからガス焚き分負荷率ｘ（ｇｄ）ＥＸＨ（ｔ）を求めるものである。従って、式（１９）から明らかであるように、発電機排熱回収分Ｈ（ｒｅｃ）ｃＥＧＥＸＨ（ｘＥＸＨ（ｔ））が熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）を上回ってしまう場合、すなわち、Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）−ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））×Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））≦０となってしまう場合には、燃料供給設備（ＦＳＦ）１０４から供給された燃料ＦＳＥＸＨ（ｔ）を使う必要はなく、ガス焚き分負荷率をｘ（ｇｄ）ＥＸＨ（ｔ）は０となる。
【００９５】
図１５（ｂ）に示されているように、ガス焚きＣＯＰ特性“ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））”は、ガス焚き分負荷率をｘ（ｇｄ）ＥＸＨ（ｔ）の関数となる。ここで、ｇ０を定数項、ｇ１、ｇ２、ｇ３を係数とすると、ＣＯＰ計算部１３７２は、ガス焚きＣＯＰ特性“ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））”を、
ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））＝ｇ０＋ｇ１・ｘ（ｇｄ）ＥＸＨ（ｔ）＋ｇ２・｛ｘ（ｇｄ）ＥＸＨ（ｔ）｝２＋ｇ３・｛ｘ（ｇｄ）ＥＸＨ（ｔ）｝３    （２０）
として求めることができる。上記と同様に、パラメータのセット［ｇ０，ｇ１，ｇ２，ｇ３］は、離散的な負荷率ｘ（ｇｄ）ＥＸＨ（ｔ）の値に対するＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））の値を内挿補間して、任意の負荷率ｘ（ｇｄ）ＥＸＨ（ｔ）についてのＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））の値を求めることができるように推計することができる。代替的に、排熱投入型吸収冷温水機（ＥＸＨ）１１２の機器特性モデル生成部１３６６は、離散的な負荷率ｘ（ｇｄ）ＥＸＨ（ｔ）の値に対するＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））の値の多数のサンプルを使用して重回帰分析を行い、パラメータのセット［ｇ０，ｇ１，ｇ２，ｇ３］を推計しても良い。
【００９６】
熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）、実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））、排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”、ガス焚きＣＯＰ特性“ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））”がわかると、所要燃料計算部１３７４は、所要燃料ＦＳＥＸＨ（ｔ）を、
ＦＳＥＸＨ（ｔ）＝ｍａｘ［Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）−｛ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））×Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））｝，０］／ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））    （２１）
として求めることができる。式（２１）は、満たすべき需要である熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）のうち実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））で賄えない部分を求め、その部分とガス焚きＣＯＰ特性“ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））”とから所要燃料ＦＳＥＸＨ（ｔ）を求めるものである。従って、式（２１）から明らかであるように、発電機排熱回収分Ｈ（ｒｅｃ）ｃＥＧＥＸＨ（ｘＥＸＨ（ｔ））が熱量設定値Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）を上回ってしまう場合、すなわち、Ｈ（ｒｅｑ）ｃＥＸＨ（ｔ）−ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））×Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ））≦０となってしまう場合には、燃料供給設備（ＦＳＦ）１０４から供給された燃料ＦＳＥＸＨ（ｔ）を使う必要はない。
【００９７】
求めた所要燃料ＦＳＥＸＨ（ｔ）は、コスト等計算部１３７６に送られる。コスト等計算部１３７６は、燃料単価をｐｆｕｅ、電力単価をｐｐｏｗ、補機類駆動用電力をＰＥＸＨＰ（ｔ）とすると、排熱投入型吸収冷温水機（ＥＸＨ）１１２の総コストをＥＸＥＸＨ（ｔ）を、
ＥＸＥＸＨ（ｔ）＝ＦＳＥＸＨ（ｔ）×ｐｆｕｅ＋ＰＥＸＨＰ（ｔ）×ｐｐｏｗ            （２２）
として求める。
【００９８】
（９）  機器特性モデル化ブロック（電動冷凍機及びガス吸収式冷温水機の場合）
図１７は、機器特性モデル化ブロック２１２を用いて電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器モデルを求める際の各機能のブロック図である。
【００９９】
電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の特性は、通常、実測もしくはシミュレーションによって得られた特性カーブを機器メーカーから入手し、これを運転計画最適化システムへ入力することによって作成される。これらの熱源機器の運転効率に影響を与える主な要因は冷却水温と負荷率であるが、メーカーから提供される機器特性は代表的な冷却水温に対する不連続的な機器特性データしかないため、他の冷却水温における特性値は、何らかの方法で推測せざるを得ないのが現状である。この点に関して、従来技術は、その推測方法を提示していない。
【０１００】
本発明においては、不連続的な機器特性データを使用して、統計的な手法に従って代表的な冷却水温以外の水温における機器特性値を推計して、パラメトリックな方法で冷却水温及び負荷率の関数として表すことでモデル化し、任意の冷却水温における機器特性値を容易に取得できるようにする。一例として、負荷率を変数として不連続的な特性カーブについて冷却水温毎に回帰分析を行い、負荷率の回帰係数を推計する。
【０１０１】
電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４等の熱源機器の成績係数（Ｃｏｅｆｆｉｃｉｅｎｔ  ｏｆ  Ｐｅｒｆｏｒｍａｎｃｅ（ＣＯＰ））は、図１８に示されるように、通常、複数の離散的な冷却水温値Ｔについて計測したＣＯＰ値をプロットしたもので与えられる。熱源機器特性データ入力ブロック２２４は、上記の離散的なＴについての不連続的なＣＯＰ値データを取得する。熱源機器特性データ入力ブロック２２４は、機器特性データとして不連続的なＣＯＰ値を機器特性モデル化ブロック２１２に供給する。
【０１０２】
機器特性モデル化ブロック２１２は、熱源機器特性データ入力ブロック２２４から各熱源機器の冷却水温Ｔ及び負荷率ｘ（ｉ，ｔ）、そして離散的な冷却水温値Ｔについての不連続的なＣＯＰ値データを取得する。機器特性モデル化ブロック２１２は、さらに、冷却水温予測ブロック２１０から冷却水温の予測値Ｔ＊（ｉ，ｔ）を取得する。
【０１０３】
例えば、図１８を参照すると、ＣＯＰ値データは、Ｔ＝１７℃、２０℃、２４℃、２８℃、３２℃について、それぞれの冷却水温ごとの曲線として得られる。機器特性モデル化ブロック２１２は、例えば、スプライン補間（ｓｐｌｉｎｅ  ｉｎｔｅｒｐｏｌａｔｉｏｎ）法を使用して、それぞれの冷却水温ごとの各曲線を、負荷率ｘ（ｉ，ｔ）の３次関数として記述する。ＣＯＰ値データの各曲線を表す関数は、例えば、
ＣＯＰ（ｉ，１７，ｘ（ｉ，ｔ））＝ａ０（ｉ，１７）＋ａ１（ｉ，１７）・ｘ（ｉ，ｔ）＋ａ２（ｉ，１７）・｛ｘ（ｉ，ｔ）｝２＋ａ３（ｉ，１７）・｛ｘ（ｉ，ｔ）｝３
ＣＯＰ（ｉ，２０，ｘ（ｉ，ｔ））＝ａ０（ｉ，２０）＋ａ１（ｉ，２０）・ｘ（ｉ，ｔ）＋ａ２（ｉ，２０）・｛ｘ（ｉ，ｔ）｝２＋ａ３（ｉ，２０）・｛ｘ（ｉ，ｔ）｝３
ＣＯＰ（ｉ，２４，ｘ（ｉ，ｔ））＝ａ０（ｉ，２４）＋ａ１（ｉ，２４）・ｘ（ｉ，ｔ）＋ａ２（ｉ，２４）・｛ｘ（ｉ，ｔ）｝２＋ａ３（ｉ，２４）・｛ｘ（ｉ，ｔ）｝３
ＣＯＰ（ｉ，２８，ｘ（ｉ，ｔ））＝ａ０（ｉ，２８）＋ａ１（ｉ，２８）・ｘ（ｉ，ｔ）＋ａ２（ｉ，２８）・｛ｘ（ｉ，ｔ）｝２＋ａ３（ｉ，２８）・｛ｘ（ｉ，ｔ）｝３
ＣＯＰ（ｉ，３２，ｘ（ｉ，ｔ））＝ａ０（ｉ，３２）＋ａ１（ｉ，３２）・ｘ（ｉ，ｔ）＋ａ２（ｉ，３２）・｛ｘ（ｉ，ｔ）｝２＋ａ３（ｉ，３２）・｛ｘ（ｉ，ｔ）｝３
として求めることができる。冷却水温Ｔ＝１７℃、２０℃、２４℃、２８℃、３２℃についての関数を決定する際に、１組のパラメータＡ＝［ａ０，ａ１，ａ２，ａ３］が得られるが、このパラメータＡは、冷却水温Ｔ＝１７℃、２０℃、２４℃、２８℃、３２℃ごとに得られる。機器特性モデル化ブロック２１２は、例えば、スプライン補間法を使用して、このパラメータＡを、Ｔ（ｉ，ｔ）の関数として、
ａ０（ｉ，Ｔ（ｉ，ｔ））＝α０＋α１・Ｔ（ｉ，ｔ）＋α２・｛Ｔ（ｉ，ｔ）｝２＋α３・｛Ｔ（ｉ，ｔ）｝３      （２３ａ）
ａ１（ｉ，Ｔ（ｉ，ｔ））＝β０＋β１・Ｔ（ｉ，ｔ）＋β２・｛Ｔ（ｉ，ｔ）｝２＋β３・｛Ｔ（ｉ，ｔ）｝３      （２３ｂ）
ａ２（ｉ，Ｔ（ｉ，ｔ））＝γ０＋γ１・Ｔ（ｉ，ｔ）＋γ２・｛Ｔ（ｉ，ｔ）｝２＋γ３・｛Ｔ（ｉ，ｔ）｝３      （２３ｃ）
ａ３（ｉ，Ｔ（ｉ，ｔ））＝δ０＋δ１・Ｔ（ｉ，ｔ）＋δ２・｛Ｔ（ｉ，ｔ）｝２＋δ３・｛Ｔ（ｉ，ｔ）｝３        （２３ｄ）
のように求める。
【０１０４】
次いで、機器特性モデル化ブロック２１２は、上記の（１１ａ）乃至（１１ｄ）によって定まるパラメータＡ＝［ａ０，ａ１，ａ２，ａ３］を冷却水温Ｔ＝１７℃、２０℃、２４℃、２８℃、３２℃について決定された関数に適用して、各熱源機器（ｉ＝１，…，ｍ）Ｔ＝１７℃、２０℃、２４℃、２８℃、３２℃についてのＣＯＰ値を、
ＣＯＰ（ｉ，Ｔ（ｉ，ｔ），ｘ（ｉ，ｔ））＝ａ０（ｉ，Ｔ（ｉ，ｔ））＋ａ１（ｉ，Ｔ（ｉ，ｔ））・ｘ（ｉ，ｔ）＋ａ２（ｉ，Ｔ（ｉ，ｔ））・｛ｘ（ｉ，ｔ）｝２＋ａ３（ｉ，Ｔ（ｉ，ｔ））・｛ｘ（ｉ，ｔ）｝３    ｆｏｒ  （ｉ，ｔ）  （ｉ＝１，…，ｍ），（ｔ＝１，…，ｎ）  （２４）
として求める。
【０１０５】
次いで、機器特性モデル化ブロック２１２は、各冷却水温Ｔ＝１７℃、２０℃、２４℃、２８℃、３２℃におけるＣＯＰ値を補間することにより、そのほかの冷却水温Ｔ（ｉ，ｔ）における値を求める。補間された関数は、各熱源機器（ｉ＝１，…，ｍ）ごとに求められ、すべての冷却水温Ｔ（ｉ，ｔ）及び負荷率ｘ（ｉ，ｔ）についての機器特性値を表す関数を求めるように補間が行われるので、機器特性モデルＣＯＰ（ｉ，Ｔ（ｉ，ｔ），ｘ（ｉ，ｔ））におけるｔの依存性は除かれる。
【０１０６】
本実施形態においては、各熱源機器（ｉ＝１，…，ｍ）ごとの熱源機器モデルＣＯＰ（ｉ，Ｔ（ｉ，ｔ），ｘ（ｉ，ｔ））を求めることは、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４等の各熱源機器（ｉ＝１，…，ｍ）ごとに行列
0006641854.tif 000007
を求めることと同じである。
【０１０７】
機器特性モデル化ブロック２１２は、求めた機器特性ＣＯＰ（ｉ，Ｔ（ｉ，ｔ），ｘ（ｉ，ｔ））を表す行列Ｖ（ｉ，ｔ）を熱源機器運転計画最適化ブロック２１６に送り、熱源機器運転計画最適化ブロック２１６は、熱源機器モデル記憶ブロック２１６ａに受け取った行列Ｖ（ｉ，ｔ）を格納する。
【０１０８】
得られた負荷率の回帰係数［ａ０，ａ１，ａ２，ａ３］は冷却水温毎に異なっているので、この負荷率を冷却水温で回帰分析することにより、負荷率の回帰係数は、冷却水温で説明することが可能となる。従って、負荷率と冷却水温がわかれば、精度の高い機器の特性値を容易に得ることが可能となる。このことにより、熱源冷凍能力をより正確に見極めることができるため、熱源機器ネットワーク全体の消費エネルギーの低減を図ることができる。
【０１０９】
図１９乃至図２１は、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器特性を示す図である。
【０１１０】
図２２は、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４のコストの計算の手順を示す機能ブロック図である。
【０１１１】
熱需要Ｈｄ（ｔ）のうち電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４に割り当てられた熱需要である熱量設定値Ｈ（ｒｅｑ）ｃＥＲＵ（ｔ）及びＨ（ｒｅｑ）ｃＧＡＢ（ｔ）が、熱量設定部１４０８及び１５０８に入力され、熱量設定値Ｈ（ｒｅｑ）ｃＥＲＵ（ｔ）及びＨ（ｒｅｑ）ｃＧＡＢ（ｔ）は、負荷率計算部１４０６及び１５０６に送られる。一方、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の定格Ｈ（ｒａｔ）ｃＥＲＵ及びＨ（ｒａｔ）ｃＧＡＢが、機器諸元格納部１４０４及び１５０４から読み出され、負荷率計算部１４０６及び１５０６に入力される。負荷率計算部１４０６及び１５０６は、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の負荷率ｘＥＲＵ（ｔ）及びｘＧＡＢ（ｔ）を計算し、負荷率ｘＥＲＵ（ｔ）及びｘＧＡＢ（ｔ）は、ＣＯＰ計算部１４２０及び１５２０に送られる。一方、ＣＯＰ計算部１４２０及び１５２０は、冷却水温格納部１４１８及び１５１８から電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の冷却水温ＴＥＲＵ（ｔ）及びＴＧＡＢ（ｔ）を取得する。ＣＯＰ計算部１４２０及び１５２０は、電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器特性を、
ＣＯＰＥＲＵ（ＴＥＲＵ（ｔ），ｘＥＲＵ（ｔ））＝ａ０（ＴＥＲＵ（ｔ））＋ａ１（ＴＥＲＵ（ｔ））・ｘＥＲＵ（ｔ）＋ａ２（ＴＥＲＵ（ｔ））・｛ｘＥＲＵ（ｔ）｝２＋ａ３（ＴＥＲＵ（ｔ））・｛ｘＥＲＵ（ｔ）｝３      （２５）
ＣＯＰＧＡＢ（ＴＧＡＢ（ｔ），ｘＧＡＢ（ｔ））＝ａ０（ＴＧＡＢ（ｔ））＋ａ１（ＴＧＡＢ（ｔ））・ｘＧＡＢ（ｔ）＋ａ２（ＴＧＡＢ（ｔ））・｛ｘＧＡＢ（ｔ）｝２＋ａ３（ＴＧＡＢ（ｔ））・｛ｘＧＡＢ（ｔ）｝３      （２６）
に従って計算する。
【０１１２】
電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器特性ＣＯＰＥＲＵ（ＴＥＲＵ（ｔ），ｘＥＲＵ（ｔ））及びＣＯＰＧＡＢ（ＴＧＡＢ（ｔ），ｘＧＡＢ（ｔ））は、所要燃料計算部１４２２及び１５２２に送られる。所要燃料計算部１４２２及び１５２２は、熱量設定部１４０８及び１５０８から受け取った熱量設定値Ｈ（ｒｅｑ）ｃＥＲＵ（ｔ）及びＨ（ｒｅｑ）ｃＧＡＢ（ｔ）と、ＣＯＰ計算部１４２０及び１５２０から受け取った機器特性ＣＯＰＥＲＵ（ＴＥＲＵ（ｔ），ｘＥＲＵ（ｔ））及びＣＯＰＧＡＢ（ＴＧＡＢ（ｔ），ｘＧＡＢ（ｔ））とを使用して、電動冷凍機（ＥＲＵ）１０８の所要電力ＰＥＲＵ（ｔ）及びガス吸収式冷温水機（ＧＡＢ）１１４の所要燃料ＦＳＧＡＢ（ｔ）を、それぞれ、
ＰＥＲＵ（ｔ）＝Ｈ（ｒｅｑ）ｃＥＲＵ（ｔ）／ＣＯＰＥＲＵ（ＴＥＲＵ（ｔ），ｘＥＲＵ（ｔ））      （２７）
ＦＳＧＡＢ（ｔ）＝Ｈ（ｒｅｑ）ｃＧＡＢ（ｔ）／ＣＯＰＧＡＢ（ＴＧＡＢ（ｔ），ｘＧＡＢ（ｔ））    （２８）
に従って計算する。ここで、電動冷凍機（ＥＲＵ）１０８の所要電力ＰＥＲＵ（ｔ）は、受電設備１０２からの電力Ｐｒｅｃｒｅ（ｔ）及び発電設備（ＥＧＥ）１１０からの電力ＰＥＧＥＲ（ｔ）との和となり、
ＰＥＲＵ（ｔ）＝Ｐｒｅｃｒｅ（ｔ）＋ＰＥＧＥＲ（ｔ）            （２９）
の関係が成り立つ。
【０１１３】
求めた所要電力ＰＥＲＵ（ｔ）及び所要燃料ＦＳＧＡＢ（ｔ）は、コスト計算部１４２４及び１５２４に送られる。電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の補機類駆動用電力をそれぞれＰＥＲＵＰ（ｔ）及びＰＧＡＢＰ（ｔ）、燃料単価をｐｆｕｅ、電力単価をｐｐｏｗとすると、電動冷凍機（ＥＲＵ）１０８の総コストＥＸＥＲＵ（ｔ）及びガス吸収式冷温水機（ＧＡＢ）１１４の総コストＥＸＧＡＢ（ｔ）は、それぞれ、
ＥＸＥＲＵ（ｔ）＝｛ＰＥＲＵ（ｔ）＋ＰＥＲＵＰ（ｔ）｝×ｐｐｏｗ      （３０）
ＥＸＧＡＢ（ｔ）＝ＦＳＧＡＢ（ｔ）×ｐｆｕｅ＋ＰＧＡＢＰ（ｔ）×ｐｐｏｗ  （３１）
に従って計算される。
【０１１４】
（１０）  熱源機器運転計画最適化ブロック
図２３及び図２４は、本発明に従った熱源機器ネットワークの運転計画の最適化方法のフローチャートである。
【０１１５】
ステップＳ１６０２において、処理が開始される。
【０１１６】
ステップＳ１６０４において、熱源機器運転計画最適化ブロック２１６は、単位時間を、ｔ＝１に設定する。
【０１１７】
ステップＳ１６０６において、熱源機器運転計画最適化ブロック２１６は、熱源機器・補機類諸元入力ブロック２２２から、機器諸元を取得する。機器諸元は、例えば、発電設備（ＥＧＥ）１１０の定格Ｐ（ｒａｔ）ＥＧ、排熱投入型吸収冷温水機（ＥＸＨ）１１２の定格値Ｈ（ｒａｔ）ｃＥＸＨ、電動冷凍機（ＥＲＵ）１０８の定格Ｈ（ｒａｔ）ｃＥＲＵ、及びガス吸収式冷温水機（ＧＡＢ）１１４の定格値Ｈ（ｒａｔ）ｃＧＡＢを含む。
【０１１８】
ステップＳ１６０８及びステップＳ１６１０において、熱源機器運転計画最適化ブロック２１６は、各熱源機器（ｉ＝１，…，ｍ）の冷却水温Ｔ（ｉ，ｔ）の実績値、冷却水温Ｔ（ｉ，ｔ）の予測値Ｔ＊（ｉ，ｔ）、外気温θｏｕｔの実績値、外気温θｏｕｔの予測値θ＊ｏｕｔ（ｔ）を取得する。冷却水温予測ブロック２１０は、式（１０）によって表される予測モデルを用いて、単位時間ｔにおける冷却水温の予測値Ｔ＊（ｉ，ｔ）を計算する。
【０１１９】
ステップＳ１６１２において、電力需要予測ブロック２０２が、周知の方法で、単位時間ｔにおける電力需要の予測値Ｐ＊ｄ（ｔ）を計算する。
【０１２０】
ステップＳ１６１４において、熱需要予測ブロック２０４が、式（９）によって表される予測モデルを用いて、単位時間ｔにおける熱需要の予測値Ｈｄ＊（ｔ）を計算する。
【０１２１】
ステップＳ１６１６において、熱源機器運転計画最適化ブロック２１６は、制約条件を取得する。制約条件の詳細については、ステップＳ１６２２のところで説明する。また、熱源機器運転計画最適化ブロック２１６は、目的関数Ｆｏｂ（ｊ，ｔ）の値を初期値Ｆ０に設定する。
【０１２２】
ステップＳ１６１８において、熱源機器運転計画最適化ブロック２１６は、運転パターンｐ（ｊ，ｔ）（ｊ＝１，…，ＮＰ）について、ｊ＝１に設定する。
【０１２３】
ステップＳ１６２０において、熱源機器運転計画最適化ブロック２１６は、運転パターンｐ（ｊ，ｔ）を設定する。運転パターンｐ（ｊ，ｔ）は、単位時間ｔにおける電力需要Ｐｄ（ｔ）及び熱需要Ｈｄ（ｔ）を満たすために、各熱源機器をどのように運転すればよいかを示す変数のセット（ベクトル）である。運転パターンｐ（ｊ，ｔ）は、単位時間ｔにおいてＮｐ個存在し、ｊ＝１，…，Ｎｐである。運転パターンｐ（ｊ，ｔ）は、受電電力Ｐｒｅｃ（ｊ，ｔ）、発電電力ＰＥＧ（ｊ，ｔ）、発電設備（ＥＧＥ）１１０の総排熱量Ｈ（ｇｒｏｓｓ）ｈＥＧ（ｘＥＧＥ（ｊ，ｔ））及び回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｊ，ｔ））、排熱投入型吸収冷温水機（ＥＸＨ）１１２での実際の排熱回収量Ｈ（ｒｅｃ）ｈＥＸＨ（ｘＥＧＥ（ｊ，ｔ），ｘＥＸＨ（ｊ，ｔ））、排熱回収型吸収冷温水機出力ＨｃＥＸＨ（ｊ，ｔ）、電動冷凍機出力ＨｃＥＲＵ（ｊ，ｔ）、及び吸収式冷温水機出力ＨｃＧＡＢ（ｊ，ｔ）に基づいて決定される。運転パターンｐ（ｊ，ｔ）は、
ｐ（ｊ，ｔ）＝［ｐｒｅｃ（ｊ，ｔ），ｐＧＥ（ｊ，ｔ），ｈＥＸＨ（ｊ，ｔ），ｈＥＲＵ（ｊ，ｔ），ｈＧＡＢ（ｊ，ｔ）］    （３２）
と表される。運転パターンｐ（ｊ，ｔ）の各値が実数かつ連続値となるのは、本発明の運転計画最適化システムが、各熱源機器の負荷率ｘ（ｉ，ｔ）に基づいてコストを評価することと密接な関連性がある。
【０１２４】
燃料供給設備（ＦＳＦ）１０４から供給された燃料をＦＳＥＸＨ（ｊ，ｔ）、発電設備（ＥＧＥ）１１０からの回収可能な排熱量ＨｈＥＧ（ｘＥＧＥ（ｔ））、排熱投入型吸収冷温水機（ＥＸＨ）１１２の負荷率ｘＥＸＨ（ｔ）、排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”、排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））、排熱投入型吸収冷温水機（ＥＸＨ）１１２のガス焚き分負荷率をｘ（ｇｄ）ＥＸＨ（ｊ、ｔ）、ガス焚きＣＯＰ特性“ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））”とから、排熱回収型吸収冷温水機出力ＨｃＥＸＨ（ｊ，ｔ）は、排熱投入型吸収冷温水機（ＥＸＨ）１１２の出力特性を表す関数をｇ（）とすると、
ＨｃＥＸＨ（ｊ，ｔ）＝ｇ（ＦＳＥＸＨ（ｊ，ｔ），ＨｈＥＧ（ｘＥＧＥ（ｊ，ｔ）），ｘＥＸＨ（ｊ，ｔ），  Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｊ，ｔ）），ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ）），ｘ（ｇｄ）ＥＸＨ（ｊ、ｔ），ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ）））    （３３）
で表される。
【０１２５】
ステップＳ１６２２において、熱源機器運転計画最適化ブロック２１６は、可能な運転パターンｐ（ｊ，ｔ）の集合Λから取得した運転パターンｐ（ｊ，ｔ）が制約条件を満たしているか否かを判定する。制約条件としては、例えば、以下に示す（１）´乃至（５）´、（１１）´乃至（３１）´、及び（３２）が考えられるが、
0006641854.tif 000008
0006641854.tif 000009
他にも、熱源機器の各々の受電電力及び機器出力が０又は上下限値の範囲内にあること、ＣＯ２の排出量が上限値以下であること、及びガス消費量が下限値以上であることを含んでもよい。
【０１２６】
集合Λから取得した運転パターンｐ（ｊ，ｔ）が制約条件を満たしている場合には、処理は、ステップＳ１６２４に進み、満たしていない場合には、処理は、ステップＳ１６２８に進む。
【０１２７】
ステップＳ１６２４において、熱源機器運転計画最適化ブロック２１６は、負荷率ｘ（ｉ，ｔ）及び冷却水温Ｔ（ｉ，ｔ）を用いて、単位時間ｔにおける電動冷凍機（ＥＲＵ）１０８及びガス吸収式冷温水機（ＧＡＢ）１１４の機器特性を表すＣＯＰ（ｉ，Ｔ（ｉ，ｔ），ｘ（ｉ，ｔ））を計算する。また、熱源機器運転計画最適化ブロック２１６は、機器特性モデル生成部１３１４から、発電設備（ＥＧＥ）１１０の排熱量を求めるモデルのパラメータのセット［ｃｋ，ｂｋ］（ｋ＝１，…，４）、発電設備（ＥＧＥ）１１０の発電機の効率ηＥＧＥ（ｘＥＧＥ（ｔ））、排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”、排熱投入型吸収冷温水機（ＥＸＨ）１１２での排熱回収容量Ｈ（ｃａｐ）ｈＥＸＨ（ｘＥＸＨ（ｔ））を求めるパラメータのセット［ｄ０，ｄ１］及び［ｅ０，ｅ１，ｅ２，ｅ３］、排熱投入型吸収冷温水機（ＥＸＨ）１１２の排熱回収ＣＯＰ特性“ＣＯＰ（ｅｘｈ）（ｘＥＸＨ（ｔ））”を求めるパラメータのセット［ｆ０，ｆ１，ｆ２，ｆ３］、排熱投入型吸収冷温水機（ＥＸＨ）１１２のガス焚きＣＯＰ特性“ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ（ｔ））”を求めるパラメータのセット［ｇ０，ｇ１，ｇ２，ｇ３］、排熱投入型吸収冷温水機（ＥＸＨ）１１２のガス焚き分負荷率ｘ（ｇｄ）ＥＸＨ（ｊ、ｔ）を取得する。これらのパラメータは、制約条件を満たしているか否かを判定するステップＳ１６２２において使用されている。
【０１２８】
ステップＳ１６２６において、熱源機器運転計画最適化ブロック２１６は、目的関数Ｆｏｂ（ｊ、ｔ）の最適化を行う。この最適化は、集合Λに含まれるすべての運転パターンｐ（ｊ，ｔ）（但し、制約条件を満たさないものを除く）について目的関数Ｆｏｂ（ｊ、ｔ）を最適化することを含み、最適化は、目的関数Ｆｏｂ（ｊ、ｔ）の最大値を求めるものであってもよく、最小値を求めるものであってもよい。
【０１２９】
例えば、目的関数Ｆｏｂ（ｊ、ｔ）は、
Ｆｏｂ（ｊ、ｔ）＝ＥＸＥＧＥ（ｊ，ｔ）＋ＥＸＥＸＨ（ｊ，ｔ）＋ＥＸＥＲＵ（ｊ，ｔ）＋ＥＸＧＡＢ（ｊ，ｔ）    （３３）
であっても良い。
【０１３０】
ここで、最適化の際の条件を表す変数の組をΓｃｏｎとし、各部でのエネルギーを表す変数の組をΓｅｎｅとすると、
Γｃｏｎ＝［θｉｎ（ｔ），θｏｕｔ（ｔ），ｈｕｍｉｎ（ｔ），ｈｕｍｏｕｔ（ｔ），ＮＰ（ｔ），ＩＶ（ｔ），ＶＱ（ｔ），θ＊ｉｎ（ｔ），θ＊ｏｕｔ（ｔ），ｈｕｍ＊ｉｎ（ｔ），ｈｕｍ＊ｏｕｔ（ｔ），ＮＰ＊（ｔ）ＶＱ＊（ｔ）］；
Γｅｎｅ＝［Ｐｒｅｃ（ｔ），Ｐｒｅｃｌ（ｔ），Ｐｒｅｃｒｅ（ｔ），ＰＰ（ｔ），ＦＳ（ｔ），ＦＳＥＧＥ（ｔ），ＦＳＥＸＨ（ｔ），ＦＳＧＡＢ（ｔ），ＰＥＧ（ｔ），ＰＥＧｌ（ｔ），ＰＥＧＥＲ（ｔ），ＨｈＥＧ（ｔ），ＨｈＥＧ（ｔ），ＨｃＥＲＵ（ｔ），ＨｃＥＸＨ（ｔ），ＨｃＧＡＢ（ｔ），ＴＥＧＥ（ｔ），ｘ（ｉ，ｔ）（ｘＥＧＥ（ｔ），ｘＥＸＨ（ｔ），ｘ（ｇｄ）ＥＸＨ（ｔ），ｘＥＲＵ（ｔ），ｘＧＡＢ（ｔ）），Ｔ（ｉ，ｔ）（ＴＥＧＥ（ｔ），ＴＥＲＵ（ｔ），ＴＥＸＨ（ｔ），ＴＧＡＢ（ｔ）），Ｐｄ（ｔ），Ｈｄ（ｔ），Ｐ＊ｄ（ｔ），Ｈ＊ｄ（ｔ）］
となる。これを用いて目的関数Ｆｏｂ（ｊ，ｔ）を表すと、
Ｆｏｂ（ｊ，ｔ）＝ｆ（ｊ，ｔ，Γｃｏｎ，Γｅｎｅ，ＣＯＰ（ｅｘｈ）（ｘＥＸＨ），ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ），ＣＯＰＥＲＵ（ＴＥＲＵ，ｘＥＲＵ），ＣＯＰＧＡＢ（ＴＧＡＢ，ｘＧＡＢ））  （３４）
となる。従って、単位時間ｔにおけるＦｏｂ（ｊ，ｔ）の最適値をｆｏｐｔ（ｔ）とすると、
0006641854.tif 000010
となる。従って、熱源機器運転計画最適化ブロック２１６は、ｆｏｐｔ（ｔ）を求めるべく、制約条件を満たす各ｊについてｆ（ｊ，ｔ，Γｃｏｎ，Γｅｎｅ，ＣＯＰ（ｅｘｈ）（ｘＥＸＨ），ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ），ＣＯＰＥＲＵ（ＴＥＲＵ，ｘＥＲＵ），ＣＯＰＧＡＢ（ＴＧＡＢ，ｘＧＡＢ））の値を計算する。
【０１３１】
目的関数は、熱源機器の電力及び熱エネルギーの提供のためのコストの合計、契約電力以下という制約条件の下での熱源機器の電力及び熱エネルギーの提供のためのコストの合計、熱源機器の電力及び熱エネルギーの提供のためのＣＯ２排出量の合計、熱源機器の電力及び熱エネルギーの提供のためのＣＯ２排出量の合計が所定の値以下という制約条件の下での熱源機器の電力及び熱エネルギーの提供のためのコストの合計、熱源機器の電力及び熱エネルギーの提供のためのガス消費量の合計が所定の値以下であるという制約条件の下での熱源機器の電力及び熱エネルギーの提供のためのコストの合計、のうちの少なくとも１つから選択されてもよい。
【０１３２】
ステップＳ１６２６における処理が終了すると、処理は、ステップＳ１６２８に進み、ｊがＮｐよりも小さいか否かが判定される。
【０１３３】
ｊがＮｐよりも小さい場合には、ステップＳ１６３０において、ｊ←ｊ＋１の処理が行われる。ｊがＮｐ以上である場合には、その単位時間ｔについてすべての運転パターンｐ（ｊ，ｔ）の判定が終了したということであるので、処理は、ステップＳ１６３２に進む。Ｆｏｂ（ｊ，ｔ）の最適値を与えるｊをｊｏｐｔ（ｔ）とすると、
ｐｏｐｔ（ｔ）＝ｐ（ｊｏｐｔ（ｔ），ｔ）          （３５）
となるので、ステップＳ１６３２において、このｊｏｐｔ（ｔ）とｐｏｐｔ（ｔ）（単位時間ｔにおける最適な運転パターン）を格納する。なお、ここで、
0006641854.tif 000011
である。
【０１３４】
ステップＳ１６３４において、熱源機器運転計画最適化ブロック２１６は、ｔがｎよりも小さいか否かを判定する。
【０１３５】
ｔがｎよりも小さい場合には、ステップＳ１６３６において、ｔ←ｔ＋１の処理が行われる。
【０１３６】
ｔがｎ以上である場合には、熱源機器運転計画最適化ブロック２１６は、ステップＳ１６３８において、最適運転パターンｐのシーケンスを、
Ｐｓ＝［ｐｏｐｔ（１），…，ｐｏｐｔ（ｎ）］      （３６）
として求める。
【０１３７】
本発明の熱源機器ネットワークのコスト関数ｆ（ｊ，ｔ，Γｃｏｎ，Γｅｎｅ，ＣＯＰ（ｅｘｈ）（ｘＥＸＨ），ＣＯＰ（ｇｄ）（ｘ（ｇｄ）ＥＸＨ），ＣＯＰＥＲＵ（ＴＥＲＵ，ｘＥＲＵ），ＣＯＰＧＡＢ（ＴＧＡＢ，ｘＧＡＢ））のモデル化及び最適化、及び運転パターンｐの最適化ｐｏｐｔ（ｔ）＝ｐ（ｊｏｐｔ（ｔ），ｔ）は、構造化ニューラルネットワーク（Ｎｅｕｒａｌ  Ｎｅｔｗｏｒｋ（ＮＮ））によって行われてもよい。
【０１３８】
ステップＳ１６４０で、最適化の処理は終了する。
【０１３９】
本発明に従った上記の方法の各ステップは、ＣＰＵとメモリとを備えた、図示しないオペレータの端末によって行われてもよい。本発明の実施形態においては、上記の方法を実行するのに用いられるソフトウェアプログラムが、デジタル・データ・ストレージ媒体等のプログラム・ストレージ・デバイスに格納されていてもよく、これらのプログラム・ストレージ・デバイスは、機械読取り可能、又はコンピュータ読取り可能であり、また、これらのコンピュータ又は機械は、プログラム命令を機械実行可能プログラム、又はコンピュータ実行可能プログラムとしてエンコードし、エンコードされたプログラム命令は、本発明の方法のステップのうちの一部又は全部を実行する。プログラム・ストレージ・デバイスは、例えば、デジタル・メモリ、磁気ディスク及び磁気テープなどの磁気ストレージ媒体、ハード・ドライブ、又は光学的読取り可能デジタル・データ・ストレージ媒体とすることができる。
【０１４０】
本発明は、ソフトウェア、及び／又はソフトウェアとハードウェアの組合せで、例えば、特定用途向け集積回路（ＡＳＩＣ：ａｐｐｌｉｃａｔｉｏｎ  ｓｐｅｃｉｆｉｃ  ｉｎｔｅｇｒａｔｅｄ  ｃｉｒｃｕｉｔ）、フィールド・プログラマブル・ゲート・アレイ（ＦＰＧＡ：ｆｉｅｌｄ  ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙ）と、汎用コンピュータ、又は他の任意のハードウェア等価物を使用して実装されてもよい。
【０１４１】
上記の説明は、単に、本発明の特定の実施形態の開示を提供しているにすぎず、本発明を上記の実施形態のみに限定するように意図されてはいない。従って、本発明は、上記で説明された実施形態だけに限定されるものではなく、むしろ、当業者が本発明の範囲内に含まれる代替の実施形態を考案し得ることが認識される。
【符号の説明】
【０１４２】
１００    熱源機器ネットワーク
１０２    受電設備
１０４    燃料供給設備
１０６    電気負荷（ＥＬ）
１０８    電動冷凍機（ＥＲＵ）
１１０    発電設備（ＥＧＥ）
１１２    排熱投入型吸収冷温水機（ＥＸＨ）
１１４    ガス吸収式冷温水機（ＧＡＢ）
１１６    空調負荷
２００    最適化システム
２０２    電力需要予測ブロック
２０４    熱需要予測ブロック
２０６    外気冷房（ナイトパージ）ブロック
２０８    比エンタルピ予測ブロック
２１０    冷却水温予測ブロック
２１２    機器特性モデル化ブロック
２１４    需要計算ブロック
２１６    熱源機器運転計画最適化ブロック
２１８    実績データ入力ブロック
２２０    計画・予測データ入力ブロック
２２２    熱源機器・補機類諸元入力ブロック
２２４    熱源機器特性データ入力ブロック
２２６    運転計画作成諸元入力ブロック
２２８    運転計画出力ブロック
【図１】
0006641854.tif 000012
【図２】
0006641854.tif 000013
【図３】
0006641854.tif 000014
【図４】
0006641854.tif 000015
【図５】
0006641854.tif 000016
【図６】
0006641854.tif 000017
【図７】
0006641854.tif 000018
【図８】
0006641854.tif 000019
【図９】
0006641854.tif 000020
【図１０】
0006641854.tif 000021
【図１１】
0006641854.tif 000022
【図１２】
0006641854.tif 000023
【図１３】
0006641854.tif 000024
【図１４】
0006641854.tif 000025
【図１５】
0006641854.tif 000026
【図１６】
0006641854.tif 000027
【図１７】
0006641854.tif 000028
【図１８】
0006641854.tif 000029
【図１９】
0006641854.tif 000030
【図２０】
0006641854.tif 000031
【図２１】
0006641854.tif 000032
【図２２】
0006641854.tif 000033
【図２３】
0006641854.tif 000034
【図２４】
0006641854.tif 000035