(57)【要約】【課題】  予測値に基づき一旦作成されたエネルギーシステムの運用計画を、予測値と実測値とのずれが明らかとなった場合に、迅速に且つ運用計画策定時の目標を損なうことなく修正する。【解決手段】  入力パラメータの一部または全部の変化に対応して電力・温熱供給手段２および冷熱供給手段３の一部または全部の構成機器の出力量を調整する修正案を予め用意すると共に、修正案について入力パラメータの変化量と調整対象となる機器の出力調整量との関係式を予め求めておき、予測したパラメータ値から真のパラメータ値への実際の入力パラメータの変化量が明らかとなった場合に、関係式に当該実際の入力パラメータの変化量を代入して調整対象機器の出力調整量を求めて、運用計画における電力・温熱供給手段２および冷熱供給手段３の一部または全部の構成機器の出力量を修正する。【選択図】        図１
2006178626.tif 000002
【特許請求の範囲】
【請求項１】
燃料を動力源として電力と温熱を出力する電力・温熱供給手段と、電力または温熱の一方または双方を動力源として冷熱を出力する冷熱供給手段とを有するエネルギーシステムについて、予め設定された前記電力・温熱供給手段および前記冷熱供給手段の特性関数と、未来の運転時間帯について予測された電力需要および温熱需要および冷熱需要を少なくとも含む入力パラメータとに基づいて、電力需要および温熱需要および冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする前記電力・温熱供給手段および前記冷熱供給手段の当該運転時間帯における出力量が定められた運用計画を修正する方法であり、入力パラメータの一部または全部の変化に対応して前記電力・温熱供給手段および前記冷熱供給手段の一部または全部の構成機器の出力量を調整する修正案候補を予め用意すると共に、前記修正案候補について入力パラメータの変化量と調整対象となる機器の出力調整量との関係式を予め求めておき、予測値から真値への実際の入力パラメータの変化量が明らかとなった場合に、前記関係式に当該実際の入力パラメータの変化量を代入して前記調整対象機器の出力調整量を求めて、前記運用計画を修正することを特徴とするエネルギーシステムの運用計画修正方法。
【請求項２】
前記電力・温熱供給手段は、燃料を動力源として電力と温熱を同時に出力する１台または複数台のコジェネレーション装置と、燃料を動力源として温熱を出力する１台または複数台のボイラー装置とを有することを特徴とする請求項１記載のエネルギーシステムの運用計画修正方法。
【請求項３】
前記冷熱供給手段は、電力を動力源として冷熱を出力する１台または複数台の電動冷凍機と、温熱を動力源として冷熱を出力する１台または複数台の吸収式冷凍機とを有することを特徴とする請求項１または２に記載のエネルギーシステムの運用計画修正方法。
【請求項４】
前記関係式の一部または全部が１次関数であり、入力パラメータの単位変化あたりの調整対象機器の出力調整量を予め求めておき、予測値から真値への実際の入力パラメータの変化量と、予め求めてある入力パラメータの単位変化あたりの調整対象機器の出力調整量との比例関係から、実際の入力パラメータの変化量に対応する前記調整対象機器の出力調整量を求めることを特徴とする請求項１から３のいずれか１つに記載のエネルギーシステムの運用計画修正方法。
【請求項５】
前記特性関数の一部または全部が２次以上の関数であることを特徴とする請求項１から４のいずれか１つに記載のエネルギーシステムの運用計画修正方法。
【請求項６】
前記修正案候補は、入力パラメータ変化の同一ケースについて複数が用意されると共に、前記電力・温熱供給手段および前記冷熱供給手段の一部の構成機器の出力量を前記関係式に基づいて調整し、当該一部以外の構成機器の出力量は前記運用計画のまま一定とするものであり、且つ当該複数の各修正案候補を選択した場合のそれぞれの前記目的関数の値を求めて、複数の修正案候補の中から当該目的関数の値が最小となる又は目的関数の悪化を目的の範囲に抑えられる修正案候補をパラメータの変化に対応するための修正案として採用することを特徴とする請求項１から５のいずれか１つに記載のエネルギーシステムの運用計画修正方法。
【請求項７】
前記入力パラメータは、電気料金および燃料料金の予測値を含み、前記目的関数は、電力需要および温熱需要および冷熱需要を満足し尚且つ前記電力・温熱供給手段および前記冷熱供給手段の稼動に必要な電気料金および燃料料金の合計を表すことを特徴とする請求項１から６のいずれか１つに記載のエネルギーシステムの運用計画修正方法。
【請求項８】
燃料を動力源として電力と温熱を出力する電力・温熱供給手段と、電力または温熱の一方または双方を動力源として冷熱を出力する冷熱供給手段とを有するエネルギーシステムについて、予め設定された前記電力・温熱供給手段および前記冷熱供給手段の特性関数と、未来の運転時間帯について予測された電力需要および温熱需要および冷熱需要を少なくとも含む入力パラメータとに基づいて、電力需要および温熱需要および冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする前記電力・温熱供給手段および前記冷熱供給手段の当該運転時間帯における出力量が定められた運用計画を修正する装置であり、入力パラメータの一部または全部の変化に対応して前記電力・温熱供給手段および前記冷熱供給手段の一部または全部の構成機器の出力量を調整する予め用意された修正案候補について、入力パラメータの変化量と調整対象となる機器の出力調整量との関係式が予め記憶されており、予測値から真値への実際の入力パラメータの変化量が入力されると、該当する前記関係式に当該実際の入力パラメータの変化量を代入して前記調整対象機器の出力調整量を求めて、前記運用計画の修正案を出力することを特徴とするエネルギーシステムの運用計画修正装置。
【請求項９】
燃料を動力源として電力と温熱を出力する電力・温熱供給手段と、電力または温熱の一方または双方を動力源として冷熱を出力する冷熱供給手段とを有するエネルギーシステムについて、予め設定された前記電力・温熱供給手段および前記冷熱供給手段の特性関数と、未来の運転時間帯について予測された電力需要および温熱需要および冷熱需要を少なくとも含む入力パラメータとに基づいて、電力需要および温熱需要および冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする前記電力・温熱供給手段および前記冷熱供給手段の当該運転時間帯における出力量が定められた運用計画を修正する装置としてコンピュータを機能させるプログラムであり、入力パラメータの一部または全部の変化に対応して前記電力・温熱供給手段および前記冷熱供給手段の一部または全部の構成機器の出力量を調整する予め用意された修正案候補について、入力パラメータの変化量と調整対象となる機器の出力調整量との関係式を予め記憶しておく手段と、予測値から真値への実際の入力パラメータの変化量の入力を受け付ける手段と、該当する前記関係式に当該実際の入力パラメータの変化量を代入して前記調整対象機器の出力調整量を求める手段と、前記運用計画の修正案を出力する手段として、コンピュータを機能させることを特徴とするエネルギーシステムの運用計画修正プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一旦作成されたエネルギーシステムの運用計画を修正する方法および装置およびプログラムに関する。さらに詳述すると、本発明は、タービン発電機や燃料電池などのコジェネレーション装置、ボイラー装置、電動冷凍機、吸収式冷凍機などの複数のエネルギー機器で構成されるエネルギーシステムに関し、予測したパラメータ値に基づき一旦作成された当該エネルギーシステムの運用計画を、予測パラメータ値と実測パラメータ値とのずれが明らかとなった場合に、迅速に且つ運用計画策定時の目標を損なうことなく修正する方法および装置およびプログラムに関する。
【背景技術】
【０００２】
図１は、一般的なエネルギーシステムの構成の一例およびエネルギーのフローの一例を示す。このエネルギーシステムは、ガスや灯油またはディーゼルなどの燃料を動力源として、電力と温熱を同時に出力する例えばタービン発電機やエンジン発電機または燃料電池などのコジェネレーション装置と、ガスや灯油またはディーゼルなどの燃料を動力源として温熱を出力するボイラー装置と、電力を動力源として冷熱を出力する電動冷凍機と、温熱を動力源として冷熱を出力する吸収式冷凍機とを有している。エネルギーシステムの構成機器をエネルギー機器とも呼ぶ。このエネルギーシステムでは、コジェネレーション装置が出力する電力で、電動冷凍機を作動するのに必要な電力および電力需要をまかない、コジェネレーション装置の出力電力の不足分を、例えば電力会社から購買する商用電力（購買電力とも呼ぶ。）によって補うようにしている。コジェネレーション装置およびボイラー装置からは温熱が出力され、この温熱は収式冷凍機と、温熱需要（例えば給湯も含む）をまかなう熱交換器に送られる。吸収式冷凍機では与えられた温熱を利用して冷熱を発生する。吸収式冷凍機と電動冷凍機とで冷熱需要をまかなっている。なお、余った温熱は、例えば冷却塔を通して捨てられる。
【０００３】
このエネルギーシステムの運用計画は、設定された目標に従って、例えば電力・冷熱・温熱の各需要を満足し尚且つ最も低コストとなるように、実際の運転時以前に予め策定される。各種システムの最適な運用計画を事前に策定する従来技術として、例えば特許文献１〜６に開示されたものがある。エネルギーシステムの運用計画は、予測値を少なからず含むパラメータを用いて作成される。運用計画を実際の運転時以前に策定するために、電力・冷熱・温熱の各需要は予測値を用いざるを得ないからである。また、最も低コストとなる運用計画を策定するためには、電気料金やガスや灯油またはディーゼルなどの燃料の料金もパラメータとして用いられるが、これらの料金も運用計画策定時と実際の運転時とで一定であるとの仮定に基づいている。予測値に基づくこれらのパラメータを、数理計画法に基づく最適化計算やエネルギーシステムの運転を模擬するシミュレーションを実行するプログラムが実装されたコンピュータに入力することで、電力・冷熱・温熱の各予測需要を満足し尚且つ最も低コストとなる各エネルギー機器の出力値が、当該コンピュータにて算定される。このように運用計画で事前に算定された各エネルギー機器の出力値に基づいて、各エネルギー機器の運転を行うようにしている。
【０００４】
【特許文献１】特開２００１−２１１６９６号
【特許文献２】特開平８−２４９００５号
【特許文献３】特開２０００−１０５０１９号
【特許文献４】特開２００３−１６３７４号
【特許文献５】特開２００３−１１１２７５号
【特許文献６】特開２００３−２１６２０５号
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、エネルギーシステムの運用計画に用いるパラメータ値は予測に基づくものであり、特に電力・冷熱・温熱の各需要は時々刻々と変化し得るものであるため、実際値とのずれは避けられない。また、電気料金などの料金パラメータも時々刻々と変化することが将来的には考えられ、計画策定時とずれが生じる場合があると予想される。特許文献１〜６の従来技術では、このずれに対して如何に各エネルギー機器の出力値を調整するかについては考慮されていない。
【０００６】
上記予測値と実際値とのずれを解消するべく、運転当日など、実際の入力パラメータの値が明らかとなった時点で、当該実際の入力パラメータの値を用いて、数理計画法に基づく最適化計算やエネルギーシステムの運転を模擬するシミュレーションを再度実行し、エネルギーシステムの運用計画を策定し直すことが考えられる。しかし、実際の入力パラメータの値が明らかとなる時点は、運転当日、例えば現実の運転時中もしくは運転時直前であり、時間的制約から、最適化計算やシミュレーションを再度実行することは極めて難しい。加えて、最適化計算やシミュレーションの再実行は、エネルギーシステムを構成する全てのエネルギー機器が対象となっており、最適な各エネルギー機器の出力値を算定するまでの計算負荷が高く、多大な計算時間を要してしまう。
【０００７】
このため、予測値と実際値とのずれが生じた場合、運用者の経験と勘に頼ってエネルギー機器の出力の調整量が決められているのが実状である。この場合、時々刻々変化するパラメータに合わせて、多くのエネルギー機器の出力を同時に調整するのは困難であるため、パラメータ変化に対応するために運転出力を同時に調整するエネルギー機器はいくつかに絞られることとなる。これにより出力調整対象のエネルギー機器をどれにするかで組み合わせができ、その組み合わせによって計画修正の方法には代替案がいくつもあることになる。しかし、エネルギー機器の運転出力の変更は相互に影響を及ぼし合い、さらにエネルギー機器同士が及ぼし合う影響の度合いも運転状態で様々なので、パラメータ変化に見合うだけの運転出力の調整量を簡単に求めることができず、コストの変化を予め見通すことが極めて難しい。
【０００８】
例えば電力需要が増えた場合、計画修正方法の代替案として、電力会社からの買電を増やすか、複数台あるコジェネレーション装置の中から１台もしくは複数台を選んで運転出力を増やす場合が考えられる。買電を増やす場合は、簡単にコストの増加分を算出できるが、コジェネレーション装置の出力を増やす場合、簡単には算出できない。即ち、電力需要の増加をまかなおうとコジェネレーション装置の出力を上げると、同時に発生する温熱も増えるので、増えた温熱分だけ吸収式冷凍機で発生する冷熱も増える。すると、冷熱の供給量が増えるので、電動冷凍機で発生させる冷熱を減少させることができ、その分だけ電力需要も減る。したがって、例えば１００ｋＷ電力需要が増えたとしても、コジェネレーション装置が本当に必要な出力増分は１００ｋＷではなく、コジェネレーション装置で電力需要の増加をまかなおうとすると、一筋縄では行かなくなる。また、複数台あるコジェネレーションは、発生する電力と温熱との関係のみならず、必要とする燃料量との関係もそれぞれ異なり、さらにはそれらの関係の特性式も非線形なので、どのコジェネレーションを使って電力需要の増加に対応すればコストが最も安いのかを直ちに求めるのは不可能である。同様に、冷熱需要が増加した場合、電動冷凍機の出力を上げることで対応したとすると、電力需要が増え、上記と同じ対応をしなくてはならない。一方、吸収式冷凍機の出力を上げて対応する場合、ボイラー装置かコジェネレーション装置の出力を上げて発生する温熱を増やす必要がある。ボイラーの出力を上げる場合は、直ちにコストの増加量が算出できるが、コジェネレーション装置の出力を上げると、発生する電力も増えて買電量を減らすか、電動冷凍機の出力を上げるか、どちらか一方もしくは両方可能となり代替案がいくつもあることになる。また、この代替案では、各エネルギー機器の変化量を算出するのが簡単ではなく、したがって、代替案のコストを算出するのは容易ではない。一方、電気料金が増加した場合、買電量を減らす方がいいが、その場合はコジェネレーション装置で発生する電力を増やす必要があり、同時に発生する温熱も増えるので、吸収式冷凍機で発生する冷熱も増える、という具合に一筋縄では行かなくなる。また、コジェネレーション装置の燃料費が増えた場合も、それに対応するようコジェネレーション装置で発生する電力と温熱を減らすと、吸収式冷凍機を始めとする各エネルギー機器の運転出力に影響を及ぼす。また、燃料費増加に対応する方法にはいくつもの代替案が存在し、各エネルギー機器の入力−出力特性が非線形であるために、簡単に代替案およびそのコストを求めることはできない。
【０００９】
よって、いくつもの代替案の中から、運用計画策定時の目標に最も沿うもの、例えばトータルコストが最小となるもの、を直ちに選択することは難しい。このため、エネルギー機器の出力の調整はパラメータ変化に対応することに専念したものとなり、必ずしも運用計画策定時の目標に沿ったものとはならず、例えばコスト低減を意識した対応とはならず、トータルコストの増加を招く虞がある。
【００１０】
そこで本発明は、予測したパラメータ値に基づき一旦作成されたエネルギーシステムの運用計画を、予測パラメータ値と実測パラメータ値とのずれが明らかとなった場合に、迅速に且つ運用計画策定時の目標を損なうことなく、修正する方法および装置およびプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１１】
かかる目的を達成するため、請求項１記載の発明は、燃料を動力源として電力と温熱を出力する電力・温熱供給手段と、電力または温熱の一方または双方を動力源として冷熱を出力する冷熱供給手段とを有するエネルギーシステムについて、予め設定された前記電力・温熱供給手段および前記冷熱供給手段の特性関数と、未来の運転時間帯について予測された電力需要および温熱需要および冷熱需要を少なくとも含む入力パラメータとに基づいて、電力需要および温熱需要および冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする前記電力・温熱供給手段および前記冷熱供給手段の当該運転時間帯における出力量が定められた運用計画を修正する方法であり、入力パラメータの一部または全部の変化に対応して前記電力・温熱供給手段および前記冷熱供給手段の一部または全部の構成機器の出力量を調整する修正案候補を予め用意すると共に、前記修正案候補について入力パラメータの変化量と調整対象となる機器の出力調整量との関係式を予め求めておき、予測値から真値への実際の入力パラメータの変化量が明らかとなった場合に、前記関係式に当該実際の入力パラメータの変化量を代入して前記調整対象機器の出力調整量を求めて、前記運用計画を修正するようにしている。
【００１２】
また、請求項８記載の発明は、燃料を動力源として電力と温熱を出力する電力・温熱供給手段と、電力または温熱の一方または双方を動力源として冷熱を出力する冷熱供給手段とを有するエネルギーシステムについて、予め設定された前記電力・温熱供給手段および前記冷熱供給手段の特性関数と、未来の運転時間帯について予測された電力需要および温熱需要および冷熱需要を少なくとも含む入力パラメータとに基づいて、電力需要および温熱需要および冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする前記電力・温熱供給手段および前記冷熱供給手段の当該運転時間帯における出力量が定められた運用計画を修正する装置であり、入力パラメータの一部または全部の変化に対応して前記電力・温熱供給手段および前記冷熱供給手段の一部または全部の構成機器の出力量を調整する予め用意された修正案候補について、入力パラメータの変化量と調整対象となる機器の出力調整量との関係式が予め記憶されており、予測値から真値への実際の入力パラメータの変化量が入力されると、該当する前記関係式に当該実際の入力パラメータの変化量を代入して前記調整対象機器の出力調整量を求めて、前記運用計画の修正案を出力するようにしている。
【００１３】
また、請求項９記載の発明は、燃料を動力源として電力と温熱を出力する電力・温熱供給手段と、電力または温熱の一方または双方を動力源として冷熱を出力する冷熱供給手段とを有するエネルギーシステムについて、予め設定された前記電力・温熱供給手段および前記冷熱供給手段の特性関数と、未来の運転時間帯について予測された電力需要および温熱需要および冷熱需要を少なくとも含む入力パラメータとに基づいて、電力需要および温熱需要および冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする前記電力・温熱供給手段および前記冷熱供給手段の当該運転時間帯における出力量が定められた運用計画を修正する装置としてコンピュータを機能させるプログラムであり、入力パラメータの一部または全部の変化に対応して前記電力・温熱供給手段および前記冷熱供給手段の一部または全部の構成機器の出力量を調整する予め用意された修正案候補について、入力パラメータの変化量と調整対象となる機器の出力調整量との関係式を予め記憶しておく手段と、予測値から真値への実際の入力パラメータの変化量の入力を受け付ける手段と、該当する前記関係式に当該実際の入力パラメータの変化量を代入して前記調整対象機器の出力調整量を求める手段と、前記運用計画の修正案を出力する手段として、コンピュータを機能させるようにしている。
【００１４】
したがって、予測値から真値への実際の入力パラメータの値が明らかとなった場合には、当該実際の入力パラメータの値を用いて数理計画法に基づく最適化計算やエネルギーシステムの運転を模擬するシミュレーションを再度実行することなく、入力パラメータの変化量すなわちパラメータの予測値と真値との実際のずれ量を入力すれば、調整対象機器の出力調整量が単純な計算により直ちに求められる。従って、運用者の経験や勘に頼ることなく、目的関数を最小とする又は目的関数の悪化を目的の範囲に抑えられる、例えばトータルコストを最小とする又はコストの増加を許容の範囲に抑えられるといった運用計画策定時の目標を損なうことなく、適切な調整対象機器の出力調整量を瞬時に得ることができる。
【００１５】
また、請求項２記載の発明は、請求項１記載のエネルギーシステムの運用計画修正方法において、前記電力・温熱供給手段は、燃料を動力源として電力と温熱を同時に出力する１台または複数台のコジェネレーション装置と、燃料を動力源として温熱を出力する１台または複数台のボイラー装置とを有するものとしている。従って、コジェネレーション装置とボイラー装置を有する一般的なエネルギーシステムについて、本発明を利用できる。
【００１６】
また、請求項３記載の発明は、請求項１または２に記載のエネルギーシステムの運用計画修正方法において、前記冷熱供給手段は、電力を動力源として冷熱を出力する１台または複数台の電動冷凍機と、温熱を動力源として冷熱を出力する１台または複数台の吸収式冷凍機とを有するものとしている。従って、電動冷凍機と吸収式冷凍機を有するエネルギーシステムについて、本発明を利用できる。
【００１７】
また、請求項４記載の発明は、請求項１から３のいずれか１つに記載のエネルギーシステムの運用計画修正方法において、前記関係式の一部または全部が１次関数であり、入力パラメータの単位変化あたりの調整対象機器の出力調整量を予め求めておき、実際の入力パラメータの変化量と、予め求めてある入力パラメータの単位変化あたりの調整対象機器の出力調整量との比例関係から、実際の入力パラメータの変化量に対応する前記調整対象機器の出力調整量を求めるようにしている。
【００１８】
したがって、予め求めてある「入力パラメータの単位変化あたりの調整対象機器の出力調整量」と、入力パラメータの実際値が明らかとなった場合に入力される「実際の入力パラメータの変化量／入力パラメータの単位変化量」とを単純に掛算することで、「実際の入力パラメータの変化量に対応する調整対象機器の出力調整量」を簡単にかつ迅速に求めることができる。
【００１９】
また、請求項５記載の発明は、請求項１から４のいずれか１つに記載のエネルギーシステムの運用計画修正方法において、前記特性関数の一部または全部が２次以上の関数であるものとしている。したがって、機器の特性を非線形関数を用いてより正確に表現することで、部分負荷のような非線形特性も考慮した最適な運用計画を作成することが可能となる。
【００２０】
また、請求項６記載の発明は、請求項１から５のいずれか１つに記載のエネルギーシステムの運用計画修正方法において、前記修正案候補は、入力パラメータ変化の同一ケースについて複数が用意されると共に、前記電力・温熱供給手段および前記冷熱供給手段の一部の構成機器の出力量を調整し、当該一部以外の構成機器の出力量は前記運用計画のまま一定とするものであり、且つ当該複数の各修正案候補を選択した場合のそれぞれの前記目的関数の値を求めて、当該目的関数の値が最小となる又は目的関数の悪化を目的の範囲に抑えられる修正案候補をパラメータの変化に対応するための修正案として採用するようにしている。
【００２１】
したがって、入力パラメータの変化に一部の構成機器の出力量を調整することで対応するので、全部の構成機器で対応する場合と比べて、出力調整量を算出するまでの計算負荷が軽減され、出力調整量の算定および運用計画の修正が迅速に行える。また、一部の構成機器の出力量を調整する運用計画の修正であるため、機器出力の調整も全部の構成機器で入力パラメータの変化に対応する場合と比べて迅速に行える。また、目的関数の値が最小となる又は目的関数の悪化を目的の範囲に抑えられる修正案が採用されるため、運用計画策定時の目標に沿った計画修正、例えばコスト最小化又はコスト増の抑制を意識した計画修正が行える。
【００２２】
また、請求項７記載の発明は、請求項１から６のいずれか１つに記載のエネルギーシステムの運用計画修正方法において、前記入力パラメータは、電気料金および燃料料金の予測値を含み、前記目的関数は、電力需要および温熱需要および冷熱需要を満足し尚且つ前記電力・温熱供給手段および前記冷熱供給手段の稼動に必要な電気料金および燃料料金の合計を表すものとしている。この場合、運用計画はトータルコストが最も低くなるように作成され、計画修正もトータルコストが最も低くなるように行われる。
【発明の効果】
【００２３】
しかして請求項１記載のエネルギーシステム運用計画の修正方法および請求項８記載のエネルギーシステム運用計画修正装置および請求項９記載のエネルギーシステム運用計画修正プログラムによれば、運用計画の基となった予測に基づく各種パラメータ、例えば電力・冷熱・温熱の各需要や電気料金や燃料の料金が、運転当日等に判明した実際値とずれている場合でも、当該運転当日に、運用計画を再度策定し直すことなく、迅速に適切な調整対象機器の出力調整量を求めることができる。運用者の経験や勘に頼ることなく、例えばトータルコストを最小とするといった運用計画策定時の目標を損なうことなく、パラメータ変化に直ちに対応した適切な運用計画の修正が行える。
【００２４】
さらに請求項２記載の発明によれば、コジェネレーション装置とボイラー装置を有する一般的なエネルギーシステムについて、本発明を利用できる。さらに請求項３記載の発明によれば、電動冷凍機と吸収式冷凍機を有するエネルギーシステムについて、本発明を利用できる。
【００２５】
さらに請求項４記載の発明によれば、予め求めてある「入力パラメータの単位変化あたりの調整対象機器の出力調整量」と、入力パラメータの実際値が明らかとなった場合に入力される「実際の入力パラメータの変化量／入力パラメータの単位変化量」とを単純に掛算することで、「実際の入力パラメータの変化量に対応する調整対象機器の出力調整量」を簡単にかつ迅速に求めることができる。
【００２６】
さらに請求項５記載の発明によれば、機器の特性を２次以上の関数を用いてより正確に表現することで、部分負荷のような非線形特性も考慮した最適な運用計画を作成することが可能となる。
【００２７】
さらに請求項６記載の発明によれば、入力パラメータの変化に一部の構成機器の出力量を調整することで対応するので、全部の構成機器で対応する場合と比べて、出力調整量を算出するまでの計算負荷が軽減され、出力調整量の算定および運用計画の修正が迅速に行える。また、一部の構成機器の出力量を調整する運用計画の修正であるため、機器出力の調整も全部の構成機器で入力パラメータの変化に対応する場合と比べて迅速に行える。また、目的関数の値が最小となる又は目的関数の悪化を目的の範囲に抑えられる修正案が採用されるため、運用計画策定時の目標に沿った計画修正、例えばコスト低減を意識した計画修正が行える。
【００２８】
さらに請求項７記載の発明によれば、運用計画はトータルコストが最も低くなるように作成され、計画修正もトータルコストが最も低くなるように行われる。
【発明を実施するための最良の形態】
【００２９】
以下、本発明の構成を図面に示す実施形態に基づいて詳細に説明する。
【００３０】
図１から図４に本発明のエネルギーシステムの運用計画修正方法の実施の一形態を示す。このエネルギーシステム１は、燃料を動力源として電力と温熱を出力する電力・温熱供給手段２と、電力または温熱の一方または双方を動力源として冷熱を出力する冷熱供給手段３とを有している。
【００３１】
このエネルギーシステム１の運用計画、即ち、ある将来の運転時間帯における電力・温熱供給手段２および冷熱供給手段３の出力量は、次のようにして当該運転時間帯以前に予め決定される。即ち、予め設定された電力・温熱供給手段２および冷熱供給手段３の特性関数と、未来の運転時間帯について予測された電力需要および温熱需要および冷熱需要を少なくとも含む入力パラメータとに基づいて、電力需要および温熱需要および冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする電力・温熱供給手段２および冷熱供給手段３の当該運転時間帯における出力量が定められる。ここで、運用計画の策定対象となる運転時間帯を計画対象時間帯と呼ぶ。計画対象時間帯は例えば１時間、３０分などの一定時間幅を１単位とし、当該１単位における入力パラメータ値は一定とする。エネルギーシステム１の運用計画は、計画対象時間帯の到来前、例えば該当する運転日の前日に、計画対象時間帯の単位毎に、予め策定される。計画対象時間帯を１単位ずつ順次進めていくことで、エネルギーシステム１の連続した運用計画、例えば翌日の運用計画が作成される。
【００３２】
エネルギーシステム１の運用計画の策定に用いる電力需要などの入力パラメータの値は予測に基づくものであり、運転当日、例えば現実の運転時中もしくは運転時直前に明らかになる現実のパラメータ値との間のずれは避けられない。そこで本実施形態では、入力パラメータの予測値と実際値とのずれを埋めるために、既に策定されたエネルギーシステム１の運用計画を次のように修正する。
【００３３】
即ち、入力パラメータの一部または全部の変化に対応して電力・温熱供給手段２および冷熱供給手段３の一部または全部の構成機器の出力量を調整する複数の修正案候補を予め用意すると共に、修正案候補について入力パラメータの変化量と調整対象となる機器の出力調整量との関係式を予め求めておき、予測値から真値への実際の入力パラメータの変化量が明らかとなった場合に、関係式に当該実際の入力パラメータの変化量を代入して調整対象機器の出力調整量を求めて、複数の修正案候補の中から目的関数の値が最小となる又は目的関数の悪化を目的の範囲に抑えられる修正案候補をパラメータの変化に対応するための修正案として採用し、運用計画における電力・温熱供給手段２および冷熱供給手段３の一部または全部の構成機器の出力量を修正するようにしている。
【００３４】
図１は、一般的なエネルギーシステム１の構成の一例およびエネルギーのフローの一例を示す。エネルギーシステム１の構成機器をエネルギー機器４とも呼ぶ。このエネルギーシステム１の電力・温熱供給手段２は、エネルギー機器４として、１台または複数台のコジェネレーション装置４Ａと、１台または複数台のボイラー装置４Ｂとを有している。コジェネレーション装置４Ａは、ガスや灯油またはディーゼルなどの燃料を動力源として、電力と温熱を同時に出力する例えばタービン発電機やエンジン発電機または燃料電池などの発電機である。ボイラー装置４Ｂは、ガスや灯油またはディーゼルなどの燃料を動力源として温熱を出力する装置である。また、エネルギーシステム１の冷熱供給手段３は、電力を動力源として冷熱を出力する１台または複数台の電動冷凍機４Ｃと、温熱を動力源として冷熱を出力する１台または複数台の吸収式冷凍機４Ｄとを有している。
【００３５】
このエネルギーシステム１では、コジェネレーション装置４Ａが出力する電力と、例えば電力会社から購買する商用電力（購買電力とも呼ぶ。）とによって、電動冷凍機４Ｃを作動するのに必要な電力および電力需要をまかなっている。尚、コジェネレーション装置４Ａの出力電力の不足分を購買電力で補うようにしている。コジェネレーション装置４Ａおよびボイラー装置４Ｂからは温熱が出力され、この温熱は吸収式冷凍機４Ｄと、温熱需要（例えば給湯も含む）をまかなう熱交換器５に送られる。吸収式冷凍機４Ｄでは与えられた温熱を利用して冷熱を発生する。吸収式冷凍機４Ｄと電動冷凍機４Ｃとで冷熱需要をまかなっている。なお、余った温熱は、例えば冷却塔６を通して捨てられる。
【００３６】
尚、図１は一般的なシステム構成例を示したものであり、本発明に適用可能なエネルギーシステム１が図１の構成に限定されるものではない。例えば、必ずしもボイラー装置４Ｂを有するものに限らず、１台または複数台のコジェネレーション装置４Ａのみで電力・温熱供給手段２を構成しても構わない。また、電動冷凍機４Ｃと吸収式冷凍機４Ｄの双方を有するものに限らず、電動冷凍機４Ｃと吸収式冷凍機４Ｄの一方のみで冷熱供給手段３を構成しても構わない。
【００３７】
エネルギーシステム１の運用計画は、例えば図４に示すように既存または新規のコンピュータ資源を利用して、作成される。エネルギーシステム１の運用計画の作成に用いるコンピュータを運用計画策定手段７とも呼ぶ。コンピュータに運用計画策定プログラムがインストールされることにより、当該コンピュータが運用計画策定手段７として機能する。
【００３８】
この運用計画策定手段７が備える記憶装置７Ａには、エネルギーシステム１を構成する各エネルギー機器４の特性関数がキーボードなどの入力装置７Ｂを介して予め利用者等により入力され、記憶されている。また、運用計画策定手段７には、例えば入力装置７Ｂを介して、計画対象時間帯に対応する可変パラメータが利用者等により入力される。運用計画策定手段７が備える演算装置７Ｃにより、記憶装置７Ａに記憶されている各エネルギー機器４の特性関数と、入力されたパラメータとに基づいて、予測電力需要および予測温熱需要および予測冷熱需要を満足し尚且つ目的関数を最小とする又は目的関数を目的の範囲とする計画対象時間帯における各エネルギー機器４の出力量が算出される。計画対象時間帯を１単位ずつ順次進めていくと共に当該計画対象時間帯に対応するパラメータを入力していくことで、例えばエネルギーシステム１の翌日の運用計画が作成される。運用計画のデータは、例えばディスプレイやプリンタまたは外部の情報処理装置と通信するためのネットワークインタフェースなどの運用計画策定手段７が備える出力装置７Ｄより出力される。
【００３９】
ここで、本実施形態の入力パラメータは、例えば計画対象時間帯における予測電力需要および予測温熱需要および予測冷熱需要を示す需要パラメータと、電気料金および燃料料金を表す料金パラメータとしている。そして、本実施形態の制約条件は、構成機器の出力上下限と電力需要および温熱需要および冷熱需要を満足することであり、目的関数は電力・温熱供給手段２および冷熱供給手段３の稼動に必要な電気料金および燃料料金の合計を表すものとしている。この場合、運用計画策定手段７に設定される各エネルギー機器４の特性関数は、例えば、コジェネレーション装置４Ａの電力出力と消費燃料費との関係式、コジェネレーション装置４Ａの電力出力と温熱出力との関係式、ボイラー装置４Ｂの温熱出力と消費燃料費との関係式、吸収式冷凍機４Ｄの冷熱出力と消費温熱との関係式、電動冷凍機４Ｃの冷熱出力と消費電力との関係式となる。また、目的関数は、購買電力の料金と、全てのコジェネレーション装置４Ａの総消費燃料費と、全てのボイラー装置４Ｂの総消費燃料費との合計で表される。運用計画策定手段７は、特性関数と入力パラメータとに基づいて、予測電力需要および予測温熱需要および予測冷熱需要を満足し尚且つ目的関数を最小とする即ち最も低コストとなる計画対象時間帯における各エネルギー機器４の出力量を算出する。
【００４０】
但し、目的関数はコストの最小化を図るものに必ずしも限定されず、例えば環境面の最適化を実現するべく二酸化炭素排出量の最小化を図るものであっても良く、さらには低コストと低二酸化炭素排出量の両立を図る目的関数を設定しても良い。コスト最小化ではなく二酸化炭素排出量の最小化を図る場合、特性関数は、例えば、コジェネレーション装置４Ａの電力出力と二酸化炭素排出量との関係式、コジェネレーション装置４Ａの電力出力と温熱出力との関係式、ボイラー装置４Ｂの温熱出力と二酸化炭素排出量との関係式、吸収式冷凍機４Ｄの冷熱出力と消費温熱との関係式、電動冷凍機４Ｃの冷熱出力と消費電力との関係式となり、目的関数は二酸化炭素の総排出量で表される。二酸化炭素排出量の最小化を図る場合、入力パラメータとして電気料金および燃料料金を表す料金パラメータを入力する必要はない。
【００４１】
運用計画策定手段７としてのコンピュータは、目的関数を最小とすること及び需要と供給をバランスさせること（すなわち電力需要および温熱需要および冷熱需要を満足すること）と構成機器の出力上下限を制約条件として、例えば数理計画法により、計画対象時間帯における各エネルギー機器４の出力量を算出する。
【００４２】
需要と供給のバランス、例えば電力需給バランスは、購買電力と全てのコジェネレーション装置４Ａの電力出力との合計が、全ての電動冷凍機４Ｃの消費電力と電力需要との合計に等しい、との制約条件で表される。また、冷熱需給バランスは、全ての電動冷凍機４Ｃの冷熱出力と全ての吸収式冷凍機４Ｄの冷熱出力との合計が冷熱需要に等しい、との制約条件で表される。また、温熱需給バランスは、全てのコジェネレーション装置４Ａの温熱出力と全てのボイラー装置４Ｂの温熱出力との合計が、全ての吸収式冷凍機４Ｄの消費温熱と熱交換器５への供給温熱との合計以上となり、熱交換器５からの出力温熱は温熱需要に等しい、との制約条件で表される。
【００４３】
エネルギー機器４の特性関数が１次関数（線形関数）の場合は線形計画問題となり、特性関数が２次以上の関数（非線形関数）の場合は非線形計画問題となる。運用計画策定手段７としてのコンピュータは、既知又は新規の線形計画法または非線形計画法により、計画対象時間帯における各エネルギー機器４の最適な出力量を算出する。
【００４４】
但し、線形計画法や非線形計画法といった数理計画法により、計画対象時間帯における各エネルギー機器４の最適な出力量を算出するものには限られない。例えば、運用計画策定手段７としてのコンピュータ上で、特性関数と入力パラメータとに基づいてエネルギーシステム１の運転を模擬したシミュレーションを行い、目的関数が目的の範囲となり、且つ電力需要および温熱需要および冷熱需要を満足するように、各エネルギー機器４の出力量を調整しながら定めるようにしても良い。
【００４５】
例えば本実施形態では、特性関数を線形関数のみならず必要に応じて非線形関数で表現し、非線形計画法により、エネルギーシステム１の運用計画を作成するようにしている。エネルギー機器４の特性を線形関数で表現する場合、例えば部分負荷によりエネルギー機器４の効率が定格出力よりも低下する特性を考慮することが難しい。この結果、エネルギー機器４に低効率での運転を強いるような運用計画が作成されかねず、かえってコスト増加を招く虞がある。本実施形態のように、エネルギー機器４の特性を非線形関数を用いてより正確に表現することで、部分負荷のような非線形特性も考慮した最適な運用計画を作成することが可能となる。
【００４６】
エネルギーシステム１の運用計画の修正も、例えば図４に示すように既存または新規のコンピュータ資源を利用して、実行する事が可能である。エネルギーシステム１の運用計画の修正に用いるコンピュータを運用計画修正装置８とも呼ぶ。コンピュータに運用計画修正プログラムがインストールされることにより、当該コンピュータが運用計画修正装置８として機能する。尚、運用計画修正装置８と運用計画策定手段７とを１一台のコンピュータで実現しても良い。
【００４７】
入力パラメータの変化の仕方には、電力需要、温熱需要、冷熱需要、電気料金、燃料料金のどれが変化するのかによって、さらにはこれらのパラメータの２以上について同じ時間帯に同時に変化する場合もあるため、複数のケースが想定される。本実施形態では、入力パラメータ変化の想定し得る各ケース毎に、複数の修正案候補を用意すると共に、各修正案候補は、エネルギーシステム１内の全てのエネルギー機器４で入力パラメータの変化に対応するのではなく、一部のエネルギー機器４を使って入力パラメータの変化に対応するものとしている。即ち、一部のエネルギー機器４の出力を調整して入力パラメータの変化に対応し、他のエネルギー機器４の出力は運用計画策定時のまま一定とする。これは、時々刻々と変化するパラメータに瞬時に対応するために、全てのエネルギー機器４の出力を同時に調整するのは非現実的であるからである。したがって、どの種のエネルギー機器４を何台使うかによって、同じケースの入力パラメータ変化であっても、修正案候補は複数あり得る。
【００４８】
入力パラメータ変化の各ケースの各修正案候補について、入力パラメータの変化量と調整対象となる機器の出力調整量との関係式を予め求めておく。ある入力パラメータがＡ→Ａ＋ΔＡと変化したときに、あるエネルギー機器４の出力がＸ→Ｘ＋ΔＸと変化したとすると、そのエネルギー機器４の出力調整量は、ΔＸである。また、パラメータ１単位あたりの増分出力はΔＸ／ΔＡとなる。ΔＡ→０としてＡについて偏微分を求めると、増分出力は∂Ｘ／∂Ａとなる。∂Ｘ／∂Ａは一般に感度係数と呼ばれる。
【００４９】
一方、テーラー展開より、次式の関係が成り立つ。
【数１】
2006178626.tif 000003
【００５０】
従って、ΔＸは次式で表される。
【数２】
2006178626.tif 000004
【００５１】
結局、入力パラメータがＡ→Ａ＋ΔＡと変化したときのエネルギー機器４の出力調整量ΔＸは、テーラー展開の１次近似であれば数式３で表され、テーラー展開の２次近似であれば数式４で表される。
【数３】
2006178626.tif 000005
【数４】
2006178626.tif 000006
【００５２】
従って、入力パラメータ変化の各ケースの各修正案候補について、数式３または数式４を予め求めておけば、予測値から真値への実際の入力パラメータの変化量ΔＡが明らかとなった場合に、上記予め求めておいた関係式に当該実際の入力パラメータの変化量ΔＡを代入することで、調整対象機器の出力調整量ΔＸを直ちに求めることができる。数式４は、数式３よりも出力調整量ΔＸの精度が高いが、算定式が複雑になるのでより多くの計算時間やより高度な計算機資源を要する。このため、例えば用途に応じて計算速度と精度のどちらを優先させるかを検討し、数式３か数式４のいずれかを適宜選択して良い。本明細書では、便宜上、数式３に基づく方法を簡易計算法と呼び、数式４に基づく方法を詳細計算法と呼ぶ。
【００５３】
簡易計算法を用いる場合、入力パラメータの変化量ΔＡと調整対象機器の出力調整量ΔＸとの関係式は１次関数（線形関数）となり、入力パラメータの単位変化あたりの調整対象機器の出力調整量を予め求めておけば、実際の入力パラメータの変化量が明らかとなった時点で、当該実際の入力パラメータの変化量と、予め求めてある入力パラメータの単位変化あたりの調整対象機器の出力調整量との比例関係から、実際の入力パラメータの変化量に対応する調整対象機器の出力調整量を直ちに求めることができる。即ち、「入力パラメータの単位変化あたりの調整対象機器の出力調整量」と、「実際の入力パラメータの変化量／入力パラメータの単位変化量」とを単純に掛算することで、「実際の入力パラメータの変化量に対応する調整対象機器の出力調整量」を簡単かつ迅速に求めることができる。
【００５４】
例えば入力パラメータの単位変化あたりの調整対象機器の出力調整量を示す感度係数∂Ｘ／∂Ａを予め求めておけば、実際の入力パラメータの変化量ΔＡと、予め求めてある感度係数∂Ｘ／∂Ａとから、数式３に基づき、実際の入力パラメータの変化量ΔＡに対応する調整対象機器の出力調整量ΔＸを直ちに求めることができる。
【００５５】
なお、「入力パラメータの単位変化あたりの調整対象機器の出力調整量」を求める方法には、例えば数理計画法により求められる各エネルギー機器４の出力を表す関数を、対象となる入力パラメータで偏微分して求められる感度係数∂Ｘ／∂Ａを利用する方法の他に、対象となる入力パラメータを１単位だけ増やして、数理計画法による最適化計算やエネルギーシステム１の運転を模擬したシミュレーションを再度実行して求める方法がある。最適化計算やシミュレーションを再度実行する場合は、各修正案候補において出力調整の対象となっていないエネルギー機器４は、その出力を運用計画策定時のままで固定しておく。出力調整の対象となっているエネルギー機器４は、出力を固定せずに最適化を実行する。上記最適化計算やシミュレーションの再実行後に得られた各エネルギー機器４の出力と運用計画策定時の出力とを比較して、その差分が「入力パラメータの単位変化あたりの調整対象機器の出力調整量」となる。尚、入力パラメータの単位変化とは、例えば電力需要パラメータの変化であれば１ｋＷの電力需要増加、例えば冷熱需要パラメータの変化であれば１冷凍トンの冷熱需要増加、のように表される。
【００５６】
入力パラメータ変化の各ケースに対応する修正案の候補は、例えば解析者等の人的能力または解析プログラム等を利用して、計画対象時間帯の到来以前、例えば計画対象時間帯の運用計画作成に併せて、例えば図２のフローチャートに示す手順に従って、予め作成される。先ず、予測値に基づく入力パラメータを用いて運用計画を作成する（ステップ１）。この運用計画をＡとする。次に、どの入力パラメータを変化させるか選択する（ステップ２）。例えばｎ個ある入力パラメータの中から、ある一つの入力パラメータｋ（ｋ＝１，２，３，・・・，ｎ）が１単位変化したとする。他の入力パラメータは固定する。次に、ステップ２で決定された入力パラメータの変化に対応する複数の修正案候補を作成する（ステップ３）。例えば入力パラメータｋの変化に対応する修正案候補をＢ（ｋ，ｉ）と表記する。ｉは１から始まる連続番号（ｉ＝１，２，３，・・・，ｍ）で、入力パラメータｋの変化に対応する想定し得る修正案候補の数ｍまで続く。ｉの値が異なれば、出力調整の対象となるエネルギー機器４の機種または台数が異なるものとなる。次に、ステップ３で作成されたすべての修正案候補について、入力パラメータの１単位の変化量と調整対象となる機器の出力調整量との関係式を求める（ステップ４）。例えば、数式３または数式４に基づき、あるいは運用計画Ａから修正案候補Ｂ（ｋ，ｉ）に計画を変更した場合の各エネルギー機器４の出力の変化量を算定し、これにより「入力パラメータｋの単位変化あたりの調整対象機器の出力調整量」を求める。
【００５７】
上記のステップ２からステップ４の処理を、入力パラメータの変化の全ての組み合わせについて行う（ステップ５）。例えば入力パラメータｋ＝１，２，３，・・・，ｎが個々に変化した場合の全て、さらに複数の入力パラメータが同時に変化した場合の全ての組み合わせについて、ステップ２からステップ４の処理を同様に行う。
【００５８】
尚、図２の処理は、解析者等の人的能力または解析プログラム等を利用して実行される。特にステップ４において、対象となる入力パラメータを１単位だけ増やして、数理計画法による最適化計算やエネルギーシステム１の運転を模擬したシミュレーションを再度実行して、「入力パラメータの単位変化あたりの調整対象機器の出力調整量」を求める場合には、例えば図２の処理を完全にコンピュータで実行することも可能である。この場合、入力パラメータの種類およびエネルギーシステム１を構成するエネルギー機器４の機種や台数は既知であるから、想定される入力パラメータの組み合わせ及びエネルギー機器４の組み合わせの全てについて（ステップ２，５，３）、ステップ４の処理を繰り返せば良い。
【００５９】
運用計画修正装置８の記憶装置８Ａには、上記のように求められた入力パラメータ変化の各ケースの各修正案候補についての入力パラメータの変化量ΔＡと調整対象機器の出力調整量ΔＸとの関係式が、例えばキーボードなどの入力装置８Ｂを介して利用者等により入力され、記憶される（図３のステップ６）。尚、図２の処理を実行したコンピュータより運用計画修正装置８の記憶装置８Ａに関係式のデータを書き込むようにしても良く、運用計画修正装置８により図２の処理を実行して、得られた関係式のデータを記憶装置８Ａに記憶するようにしても良い。
【００６０】
そして、例えば入力パラメータの真値が明らかとなった時点で、当該入力パラメータの１単位の変化量、すなわち入力パラメータの予測値と真値の差分量が、例えば入力装置８Ｂを介して、利用者等により運用計画修正装置８に入力される（ステップ７）。運用計画修正装置８が備える演算装置８Ｃにより、変化が起きた入力パラメータの種別に対応する複数の修正案候補が選択され、選択された各修正案候補について、記憶装置８Ａに記憶されている入力パラメータの１単位の変化量と調整対象機器の出力調整量との関係式に、入力されたパラメータ変化量が代入されて、調整対象機器の出力調整量が算出される（ステップ８）。
【００６１】
ここで本実施形態では、複数の修正案候補の中から次のようにして一つを採用するようにする。即ち、当該複数の各修正案候補を選択した場合のそれぞれの目的関数の値を求めて、当該目的関数の値が最小となる又は目的関数の悪化を目的の範囲に抑えられる修正案候補を採用する。本実施形態の目的関数はコストを表すため、最も低コストとなる修正案候補が採用される。具体的には、各修正案候補に対応する調整対象機器の出力調整量を、既定の運用計画における該当するエネルギー機器４の出力値に足し込み、目的関数を再計算する。修正案候補同士で目的関数値を比較し、最もコストが低くなる修正案候補を選択する（ステップ９）。
【００６２】
選択された修正案に対応する調整対象機器の出力調整量は、例えば計画修正データとして、ディスプレイやプリンタまたは外部の情報処理装置と通信するためのネットワークインタフェースなどの運用計画修正装置８が備える出力装置８Ｄより出力される（ステップ１０）。尚、運用計画策定手段７が出力した運用計画データを、運用計画修正装置８が読み込んで、選択された修正案に対応する調整対象機器の出力調整量を、既定の運用計画における該当するエネルギー機器４の出力値に反映させ、運用計画データを修正するようにしても良い。なお、修正案において調整対象機器として選択されていないエネルギー機器４については、運用計画策定時の出力値のままである。そして、修正後の運用計画に基づいて、エネルギーシステム１の運転が行われる。尚、修正後の運用計画データに係る各エネルギー機器４の出力値を、エネルギーシステム１の制御部に入力して、エネルギーシステム１を自動運転するようにしても良い。さらに、運用計画策定手段７が出力する運用計画データと、運用計画修正装置８が出力する修正データとをエネルギーシステム１の制御部に入力して、エネルギーシステム１を自動運転するようにしても良い。
【００６３】
ここで、複数の修正案候補の中から選択された修正案が、エネルギーシステム１の運用上の問題等（例えば修正後のエネルギー機器４の出力値が上限値を超えてしまっている又は下限値を下回ってしまっている等）により、実行できないものに該当する場合は、当該選択された修正案候補を微調整する（例えばエネルギー機器４の出力値の制限を満足するように調整対象機器の出力調整量を微調整する）ようにしても良く、あるいは次にコストが低くなる修正案候補を選択し直すようにしても良い。さらに、複数の修正案候補および各修正案候補を選択した場合の目的関数値（本実施形態ではトータルコスト）を、運用計画修正装置８が備えるディスプレイなどの出力装置８Ｄに表示して、運用者の判断で修正案を選択できるようにしても良い。
【００６４】
また、複数の修正案候補について、入力パラメータの単位変化あたりの調整対象機器の出力調整量に基づくコストの大きさの順位と、実際の入力パラメータの変化量に対する調整対象機器の出力調整量に基づくコストの大きさの順位が変わらないのであれば、運用計画を修正する時点でコスト比較により複数の修正案候補の中から１つに絞る必要は必ずしも無く、該当運転時の到来前に予め最もコストの低い修正案を決定しておいても良い。この場合、入力パラメータ変化の各ケース毎に想定される複数の修正案候補のすべてについて、入力パラメータの変化量ΔＡと調整対象機器の出力調整量ΔＸとの関係式を、運用計画修正装置８に予め入力しておく必要はなく、入力パラメータ変化の各ケース毎に予め決定された最も低コストとなる１つの修正案について、入力パラメータの変化量ΔＡと調整対象機器の出力調整量ΔＸとの関係式を、運用計画修正装置８に予め入力しておけばよい。
【００６５】
以下に、エネルギーシステム１の運用計画作成方法および作成された運用計画の修正方法について、具体例を挙げて、さらに詳細に説明する。
【００６６】
本実施形態では目的関数を例えば数式５で表現し、電力需給バランスによる制約条件を数式６で表現し、冷熱需給バランスによる制約条件を数式７で表現し、温熱需給バランスを数式８および数式９で表現する。
【００６７】
【数５】
2006178626.tif 000007
【数６】
2006178626.tif 000008
【数７】
2006178626.tif 000009
【数８】
2006178626.tif 000010
【数９】
2006178626.tif 000011
【００６８】
また、購買電力料金を数式１０に示す１次式で表現する。
【数１０】
2006178626.tif 000012
【００６９】
ｉ番目コジェネレーション装置４Ａｉの電力出力と消費燃料費との関係を例えば数式１１に示す２次式で表現し、コジェネレーション装置４Ａｉの電力出力と温熱出力との関係を例えば数式１２に示す２次式で表現する。
【数１１】
2006178626.tif 000013
【数１２】
2006178626.tif 000014
【００７０】
また、ｍ番目吸収式冷凍機４Ｄｍの冷熱出力と消費温熱との関係を例えば数式１３に示す２次式で表現し、ｎ番目電動冷凍機４Ｃｎの冷熱出力と消費電力との関係を例えば数式１４に示す２次式で表現し、ｊ番目ボイラー装置４Ｂｊの温熱出力と消費燃料費との関係を例えば数式１５に示す２次式で表現する。
【数１３】
2006178626.tif 000015
【数１４】
2006178626.tif 000016
【数１５】
2006178626.tif 000017
【００７１】
計画対象時間帯における予測に基づく入力パラメータと、数式１０から数式１５に示す特性関数に基づいて、数式６〜数式９に示す制約条件を満足し、且つ数式５に示す目的関数を最小とする計画対象時間帯における各エネルギー機器４の最適な出力量を算出する。特性関数に非線形関数を用いているため、この最適化問題は非線形計画問題である。非線形計画問題の最適解の求解法は、ラグランジュ関数を使ったラグランジュの未定定数法（参考文献；相吉、志水：「数理計画法演習」、朝倉書店、1985年、福島：「数理計画入門」、朝倉書店、1996年など）が良く用いられる。この方法は、非線形計画問題において最適もしくは局所最適解で成立する条件である、クーンタッカー条件から最適解の算定式を導くものである。
【００７２】
本実施形態の非線形計画問題におけるラグランジュ関数は、次式で表現される。尚、本実施形態では、「熱交換器５に供給される温熱」および「熱交換器５から供給される温熱」および「温熱需要」は等しい、即ちＱＩＨＥ＝ＱｈＨＥ＝ＱｈＤとする。
【数１６】
2006178626.tif 000018
【００７３】
ここに、λ１,λ２,μはラグランジュ乗数と呼ばれる係数である。また、常にμ≧０で、数式９の不等号が成立している場合にμ＝０となり、等号が成立している場合はμ＞０である。本実施形態の非線形計画問題におけるクーンタッカー条件は、数式６〜数式９および数式１６のラグランジュ関数を使って表現される以下の数式１７〜数式２３で記述される。
【数１７】
2006178626.tif 000019
【数１８】
2006178626.tif 000020
【数１９】
2006178626.tif 000021
【数２０】
2006178626.tif 000022
【数２１】
2006178626.tif 000023
【数２２】
2006178626.tif 000024
【数２３】
2006178626.tif 000025
【００７４】
数式１７より次式が成立する。
【数２４】
2006178626.tif 000026
【００７５】
数式１９より次式が成立する。
【数２５】
2006178626.tif 000027
【００７６】
また数式２０と数式２１を整理すると次式が成立する。
【数２６】
2006178626.tif 000028
【００７７】
最適解におけるラグランジュ乗数をλＯ１，λＯ２，μＯとして、数式１０〜数式１５を、数式１７〜数式２１に代入すると、クーンタッカー条件は以下の数式２７〜数式３３となる。
【数２７】
2006178626.tif 000029
【数２８】
2006178626.tif 000030
【数２９】
2006178626.tif 000031
【数３０】
2006178626.tif 000032
【数３１】
2006178626.tif 000033
【数３２】
2006178626.tif 000034
【数３３】
2006178626.tif 000035
【００７８】
これにより、コジェネレーション装置４Ａｉの電力出力は数式３４で表現され、購買電力は数式３５で表現され、ボイラー装置４Ｂｊの温熱出力は数式３６で表現され、吸収式冷凍機４Ｄｍの冷熱出力は数式３７で表現され、電動冷凍機４Ｃｎの冷熱出力は数式３８で表現される。
【数３４】
2006178626.tif 000036
【数３５】
2006178626.tif 000037
【数３６】
2006178626.tif 000038
【数３７】
2006178626.tif 000039
【数３８】
2006178626.tif 000040
【００７９】
ここに、λＯ１は数式３９で表現される。
【数３９】
2006178626.tif 000041
【００８０】
λＯ２は、数式３０と数式３１を数式７に代入して、次式に示すようにμＯの関数となる。
【数４０】
2006178626.tif 000042
【００８１】
結局、最適運用状態における各エネルギー機器４の出力はμＯを決定することで一意に定められる。μＯは例えば二分法などの方法で求めることができる。従って、エネルギーシステム１の１日の運用計画を策定するには、計画対象時間帯毎に、予測電力需要、予測温熱需要、予測冷熱需要を数式３４〜数式４０に入力し、μＯを例えば二分法などで求める、という操作を繰り返せばよい。
【００８２】
なお、商用電力からの電力を購入しない場合は、λＯ１は定数にならない。その場合、数式２８を、Ｅ＝０とした数式６に代入することで、次式が得られる。この場合、λＯ１は定数ではなく、μＯの関数となる。
【数４１】
2006178626.tif 000043
【００８３】
エネルギーシステム１における各種入力パラメータ（電力需要、冷熱需要、温熱需要、電気料金、燃料料金）が変化した事例およびその場合に想定される修正案候補の例を以下に示し、修正案候補毎の各エネルギー機器４の出力調整量を定式化する。但し、以下に事例および修正案候補は一例であって、修正案候補の全てを示したものではない。ここで、一部のケースについての計算式は、簡易計算法としてテーラー展開１次近似、詳細計算法としてテーラー展開２次近似の２種類を示す。また、本事例では、ボイラー装置４Ｂを稼動するとコストメリットが出ないため、ボイラー装置４Ｂの増分はないものとした。
【００８４】
＜事例１：冷熱需要が変化するケース＞
事例１は、策定した運用計画案から冷熱需要ＱｃＤが変化した場合、すなわちＱｃＤ→ＱｃＤ＋ΔＱｃＤとなった場合である。この場合、修正案候補は大きく分けて、以下に示す４つがある。但し、複数台ある同種のエネルギー機器４の中からどの機器を選択するかで、４つの修正案候補がさらに細分化される場合もある。
【００８５】
＜事例１の修正案候補１＞
修正案候補１では、購買電力と電動冷凍機４Ｃｎで冷熱需要の変化に対応する。より具体的には、冷熱需要の変化分はNRE台ある電動冷凍機４Ｃのある１台の出力調整で対応し、それにより変化した消費電力は商用電力からの購買電力を調整して対応する。他のエネルギー機器４の出力は運用計画案策定時のままで一定とする。この条件は、以下の式で表現される。
【数４２】
2006178626.tif 000044
【数４３】
2006178626.tif 000045
【数４４】
2006178626.tif 000046
【００８６】
＜事例１の修正案候補１の簡易計算法＞
数式４２と数式４４から次式が成立する。
【数４５】
2006178626.tif 000047
【００８７】
また、数式４５、数式１４、数式３１から次式が成立する。
【数４６】
2006178626.tif 000048
【００８８】
したがって、購買電力の増分は次式で表される。尚、電動冷凍機４Ｃｎの冷熱出力の増分は、数式４２で表される。
【数４７】
2006178626.tif 000049
【００８９】
＜事例１の修正案候補１の詳細計算法＞
テーラー展開の２次近似より次式が導かれる。
【数４８】
2006178626.tif 000050
【００９０】
数式４５と数式４８より、購買電力の増分は次式で表される。
【数４９】
2006178626.tif 000051
【００９１】
＜事例１の修正案候補２＞
修正案候補２では、コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄｍで冷熱需要の変化に対応する。より具体的には、冷熱需要の変化分はNRS台ある吸収式冷凍機４Ｄのある１台の出力調整で対応し、それにより変化した消費温熱はNCU台あるコジェネレーション装置４Ａのある１台を調整して対応する。また、コジェネレーション装置４Ａｉの出力を上げた分、商用電力からの購買電力を減らす。他のエネルギー機器４の出力は運用計画案策定時のままで一定とする。この条件は、以下の式で表現される。
【数５０】
2006178626.tif 000052
【数５１】
2006178626.tif 000053
【数５２】
2006178626.tif 000054
【数５３】
2006178626.tif 000055
【００９２】
＜事例１の修正案候補２の簡易計算法＞
数式５３より次式が成立する。
【数５４】
2006178626.tif 000056
【００９３】
また、数式１２、数式２８より次式が成立する。
【数５５】
2006178626.tif 000057
【００９４】
また、数式１３、数式３０、数式５０より次式が成立する。
【数５６】
2006178626.tif 000058
【００９５】
数式５５、数式５６を数式５４に代入すると、コジェネレーション装置４Ａｉの電力出力の増分は次式で表される。尚、吸収式冷凍機４Ｄｍの冷熱出力の増分は数式５０で表され、購買電力の増分は数式５２で表される。
【数５７】
2006178626.tif 000059
【００９６】
＜事例１の修正案候補２の詳細計算法＞
テーラー展開の２次近似より次式が導かれる。
【数５８】
2006178626.tif 000060
【数５９】
2006178626.tif 000061
【００９７】
数式５４、数式５８、数式５９より、コジェネレーション装置４Ａｉの電力出力の増分は次式で表される。
【数６０】
2006178626.tif 000062
【数６１】
2006178626.tif 000063
【００９８】
尚、数式１２より次式が導かれる。
【数６２】
2006178626.tif 000064
【００９９】
また、数式６２と数式５５から次式が導かれる。従って、数式５８が導出される。
【数６３】
2006178626.tif 000065
【数６４】
2006178626.tif 000066
【０１００】
＜事例１の修正案候補３＞
修正案候補３では、コジェネレーション装置４Ａｉと電動冷凍機４Ｃｎで冷熱需要の変化に対応する。より具体的には、冷熱需要の変化分はNRE台ある電動冷凍機４Ｃのある１台の出力調整で対応し、それにより変化した消費電力はNCU台あるコジェネレーション装置４Ａのある１台を調整して対応する。コジェネレーション装置４Ａｉの出力調整により、温熱が余った場合は捨てる。尚、温熱が足らない場合はボイラー装置４Ｂの温熱出力を増せばよいが、本事例では、ボイラー装置４Ｂの増分はないものとした。他のエネルギー機器４の出力は運用計画案策定時のままで一定とする。この条件は、以下の式で表現される。
【数６５】
2006178626.tif 000067
【数６６】
2006178626.tif 000068
【数６７】
2006178626.tif 000069
【０１０１】
＜事例１の修正案候補３の簡易計算法＞
数式６５、数式６７より次式が成立する。
【数６８】
2006178626.tif 000070
【０１０２】
また、数式１４、数式３１から次式が成立する。
【数６９】
2006178626.tif 000071
【０１０３】
したがって、コジェネレーション装置４Ａｉの電力出力の増分は数式６８と数式６９より次式で表される。尚、電動冷凍機４Ｃｎの冷熱出力の増分は、数式６５で表される。
【数７０】
2006178626.tif 000072
【０１０４】
＜事例１の修正案候補３の詳細計算法＞
テーラー展開の２次近似と数式６５より次式が導かれる。
【数７１】
2006178626.tif 000073
【０１０５】
数式６８と数式７１より、コジェネレーション装置４Ａｉの電力出力の増分は次式で表される。
【数７２】
2006178626.tif 000074
【０１０６】
＜事例１の修正案候補４＞
修正案候補４では、コジェネレーション装置４Ａｉと複数台の吸収式冷凍機４Ｄと複数台の電動冷凍機４Ｃで冷熱需要の変化に対応する。より具体的には、冷熱需要の変化分は、全て若しくは一部の吸収式冷凍機４Ｄと電動冷凍機４Ｃの出力調整で対応し、それにより変化した消費温熱および消費電力はNCU台あるコジェネレーション装置４Ａのある１台を調整して対応する。他のエネルギー機器４の出力は運用計画案策定時のままで一定とする。この条件は、以下の式で表現される。但し、NRS,NREは出力調整の対象となる吸収式冷凍機４Ｄと電動冷凍機４Ｃの台数とする。
【数７３】
2006178626.tif 000075
【数７４】
2006178626.tif 000076
【数７５】
2006178626.tif 000077
【数７６】
2006178626.tif 000078
【数７７】
2006178626.tif 000079
【数７８】
2006178626.tif 000080
【数７９】
2006178626.tif 000081
【数８０】
2006178626.tif 000082
【０１０７】
＜事例１の修正案候補４の簡易計算法＞
数式７３より次式が成立する。
【数８１】
2006178626.tif 000083
【０１０８】
また、数式７９、数式１３、数式３０、数式１２、数式２８から次式が成立する。
【数８２】
2006178626.tif 000084
【０１０９】
また、数式７７、数式１４、数式３１から次式が成立する。
【数８３】
2006178626.tif 000085
【０１１０】
数式８２、数式８３を数式８１に代入すると、コジェネレーション装置４Ａｉの電力出力の増分は次式で表される。
【数８４】
2006178626.tif 000086
【０１１１】
また、数式８３から複数台の電動冷凍機４Ｃの冷熱出力の増分は次式で表される。
【数８５】
2006178626.tif 000087
【０１１２】
また、数式８２から複数台の吸収式冷凍機４Ｄの冷熱出力の増分は次式で表される。
【数８６】
2006178626.tif 000088
【０１１３】
＜事例１の修正案候補４の詳細計算法＞
電動冷凍機が複数台集まったときの冷熱出力と消費電力の関係を数式８７とする。
【数８７】
2006178626.tif 000089
数式７７、数式７３より次式が成立する。
【数８８】
2006178626.tif 000090
【０１１４】
また、数式７７、数式７９、数式１３、数式３０、数式１４、数式３１より次式が成立する。
【数８９】
2006178626.tif 000091
【０１１５】
また、テーラー展開の２次近似より次式が導かれる。
【数９０】
2006178626.tif 000092
【０１１６】
数式９０を数式８９に代入して次式が得られる。
【数９１】
2006178626.tif 000093
【０１１７】
また、テーラー展開の２次近似より次式が導かれる。
【数９２】
2006178626.tif 000094
【０１１８】
ここで、数式９２の右辺の第２項は第１項よりも充分に小さいことから、第２項のΔＱＲＥの値として、次式に示すテーラー展開の１次近似を用いることにする。
【数９３】
2006178626.tif 000095
【０１１９】
数式９３を数式９２に代入して次式が得られる。
【数９４】
2006178626.tif 000096
【０１２０】
数式９１、数式９４を数式８８に代入して、次式が得られる。
【数９５】
2006178626.tif 000097
【数９６】
2006178626.tif 000098
【０１２１】
数式９５にはΔｐＣＵｉが含まれるが、ΔｐＣＵｉの項の影響は小さいことから、この項のΔｐＣＵｉの値として、次式の１次近似を用いることにする。
【数９７】
2006178626.tif 000099
【０１２２】
数式９７、数式９６を数式９５に代入すると、次式が得られる。
【数９８】
2006178626.tif 000100
【０１２３】
数式９８よりコジェネレーション装置４Ａｉの電力出力の増分は次式で表される。
【数９９】
2006178626.tif 000101
【０１２４】
また、数式７９、数式１３、数式３０より、複数台の吸収式冷凍機４Ｄの冷熱出力の増分は次式で表される。
【数１００】
2006178626.tif 000102
【０１２５】
また、数式７３より複数台の電動冷凍機４Ｃの冷熱出力の増分は次式で表される。
【数１０１】
2006178626.tif 000103
【０１２６】
コジェネレーション装置４Ａｉの温熱出力の増分は次式で表される。
【数１０２】
2006178626.tif 000104
【０１２７】
＜事例２：電力需要が変化するケース＞
事例２は、策定した運用計画案から電力需要ＰＤが変化した場合、すなわちＰＤ→ＰＤ＋ΔＰＤとなった場合である。この場合、修正案候補は大きく分けて、以下に示す２つがある。但し、複数台ある同種のエネルギー機器４の中からどの機器を選択するかで、２つの修正案候補がさらに細分化される場合もある。
【０１２８】
＜事例２の修正案候補１＞
修正案候補１では、コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄと電動冷凍機４Ｃとで電力需要の変化に対応する。より具体的には、コジェネレーション装置４Ａｉの電力出力を電力需要の変化に合わせて変える。それに伴ってコジェネレーション装置４Ａｉからの温熱出力も変わるので、その分だけ吸収式冷凍機４Ｄの冷熱出力も変わる。さらに、吸収式冷凍機４Ｄの冷熱出力が変わった分だけ電動冷凍機４Ｃの冷熱出力を変える。他のエネルギー機器４の出力は運用計画案策定時のままで一定とする。この条件は、以下の式で表現される。尚、冷凍機の入出力はそれぞれ等価機器１台に総合して記述する。
【数１０３】
2006178626.tif 000105
【数１０４】
2006178626.tif 000106
【数１０５】
2006178626.tif 000107
【数１０６】
2006178626.tif 000108
【数１０７】
2006178626.tif 000109
【数１０８】
2006178626.tif 000110
【数１０９】
2006178626.tif 000111
【数１１０】
2006178626.tif 000112
【０１２９】
＜事例２の修正案候補１の簡易計算法＞
数式１０２より次式が導かれる。
【数１１１】
2006178626.tif 000113
【０１３０】
また、数式１０５、数式１１０より次式が導かれる。
【数１１２】
2006178626.tif 000114
【０１３１】
従って、数式１２、数式２８、数式２１、数式２２より、次式が成立する。
【数１１３】
2006178626.tif 000115
【０１３２】
数式１１３を数式１１１に代入すると、次式が得られる。
【数１１４】
2006178626.tif 000116
【０１３３】
従って、コジェネレーション装置４Ａｉの電力出力の増分は次式で表される。
【数１１５】
2006178626.tif 000117
【０１３４】
また、吸収式冷凍機４Ｄ、電動冷凍機４Ｃの冷熱出力の増分は次式より求められる。
【数１１６】
2006178626.tif 000118
【数１１７】
2006178626.tif 000119
【０１３５】
＜事例２の修正案候補１の詳細計算法＞
２次偏微分を求めるのは複雑な計算が必要となるため、詳細計算法は省略する。
詳細計算法は行わない。
【０１３６】
＜事例２の修正案候補２＞
修正案候補２では、電力需要の変化分は購買電力で賄う。全てのエネルギー機器４の出力は運用計画案策定時のままで一定とする。購買電力の増分は次式で表される。尚、この場合は、１次式であるため詳細計算法は行えない。
【数１１８】
2006178626.tif 000120
【０１３７】
＜事例３：電気料金が変化するケース＞
事例３は、策定した運用計画案から電気料金αが変化した場合、すなわちα→α＋Δαとなった場合である。この場合、修正案候補として例えば以下に示すものがある。但し、複数台ある同種のエネルギー機器４の中からどの機器を選択するかで、修正案候補がさらに細分化される場合もある。
【０１３８】
＜事例３の修正案候補＞
コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄｍと電動冷凍機４Ｃｎとで電気料金の変化に対応する。他のエネルギー機器４の出力は運用計画案策定時のままで一定とする。この条件は、以下の式で表現される。
【数１１９】
2006178626.tif 000121
【０１３９】
コジェネレーション装置４Ａｉの温熱出力が変わり、吸収式冷凍機４Ｄｍの冷熱出力もそれにより変わる。また、冷熱需要は変わらないので、冷熱供給量が変わらないよう吸収式冷凍機４Ｄｍの冷熱出力変動分に合わせて電動冷凍機４Ｃｎの出力が変わる。コジェネレーションｉの電力も同時に変わるが、電力需要は変わらないので、コジェネレーションｉの出力変動分だけ購買電力が変わる。
【０１４０】
また、電気料金を１単位だけ増やした時のμＯの増分をΔＭαとする。ΔＭαは、例えばコジェネレーションｉおよび吸収式冷凍機４Ｄｍおよび電動冷凍機４Ｃｎ以外のエネルギー機器４の出力を運用計画策定時のままに固定し、且つΔα＝１として、運用計画策定手段７による最適化計算を再実行することで求める。
【０１４１】
ここで、以下の数式を定義する。
【数１２０】
2006178626.tif 000122
【数１２１】
2006178626.tif 000123
【数１２２】
2006178626.tif 000124
【０１４２】
また、数式３４より次式が導かれる。
【数１２３】
2006178626.tif 000125
【０１４３】
したがって、コジェネレーション装置４Ａｉの電力出力の増分は次式より求められる。
【数１２４】
2006178626.tif 000126
【０１４４】
また、コジェネレーション装置４Ａｉの温熱出力の増分は次式より求められる。
【数１２５】
2006178626.tif 000127
【０１４５】
また、吸収式冷凍機４Ｄｍの冷熱出力の増分は次式より求められる。
【数１２６】
2006178626.tif 000128
【０１４６】
また、電動冷凍機４Ｃｎの冷熱出力の増分は次式で表される。
【数１２７】
2006178626.tif 000129
【０１４７】
電動冷凍機４Ｃｎの消費電力の増分は次式より求められる。
【数１２８】
2006178626.tif 000130
【０１４８】
従って商用電力からの購買電力の増分は次式で表される。
【数１２９】
2006178626.tif 000131
【０１４９】
＜事例４：燃料料金が変化するケース＞
事例４は、燃料料金β（例えばガス料金）が変化した場合、すなわちβ→β＋Δβとなった場合である。この場合、修正案候補として例えば以下に示すものがある。但し、複数台ある同種のエネルギー機器４の中からどの機器を選択するかで、修正案候補がさらに細分化される場合もある。
【０１５０】
＜事例４の修正案候補＞
コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄｍと電動冷凍機４Ｃｎとで燃料料金の変化に対応する。他のエネルギー機器４の出力は運用計画案策定時のままで一定とする。この条件は、以下の式で表現される。
【数１３０】
2006178626.tif 000132
【０１５１】
コジェネレーション装置４Ａｉの温熱出力が変わり、吸収式冷凍機４Ｄｍの冷熱出力もそれにより変わる。また、冷熱需要は変わらないので、冷熱供給量が変わらないよう吸収式冷凍機４Ｄｍの冷熱出力変動分に合わせて電動冷凍機４Ｃｎの出力が変わる。コジェネレーションｉの電力も同時に変わるが、電力需要は変わらないので、コジェネレーションｉの出力変動分だけ購買電力が変わる。
【０１５２】
また、燃料料金を１単位だけ増やした時のμＯの増分をΔＭβとする。ΔＭβは、例えばコジェネレーションｉおよび吸収式冷凍機４Ｄｍおよび電動冷凍機４Ｃｎ以外のエネルギー機器４の出力を運用計画策定時のままに固定し、且つΔβ＝１として、運用計画策定手段７による最適化計算を再実行することで求める。
【０１５３】
燃料料金を考慮するため、数式１１に示したコジェネレーション装置４Ａｉの燃料費関数を次式に書き換える。
【数１３１】
2006178626.tif 000133
【０１５４】
ここで、以下の数式を定義する。
【数１３２】
2006178626.tif 000134
【数１３３】
2006178626.tif 000135
【０１５５】
また、数式３４より次式が導かれる。
【数１３４】
2006178626.tif 000136
【０１５６】
したがって、数式１３４より次式が導かれる。
【数１３５】
2006178626.tif 000137
【０１５７】
数式１３５に基づき、コジェネレーション装置４Ａｉの電力出力の増分は次式より求められる。
【数１３６】
2006178626.tif 000138
【０１５８】
また、コジェネレーション装置４Ａｉの温熱出力の増分は次式より求められる。
【数１３７】
2006178626.tif 000139
【０１５９】
また、吸収式冷凍機４Ｄｍの冷熱出力の増分は次式より求められる。
【数１３８】
2006178626.tif 000140
【０１６０】
また、電動冷凍機４Ｃｎの冷熱出力の増分は次式で表される。
【数１３９】
2006178626.tif 000141
【０１６１】
電動冷凍機４Ｃｎの消費電力の増分は次式より求められる。
【数１４０】
2006178626.tif 000142
【０１６２】
従って商用電力からの購買電力の増分は次式で表される。
【数１４１】
2006178626.tif 000143
【実施例】
【０１６３】
上記に例示した各事例の如く実際に入力パラメータを運用計画策定時から変えて、上記のように定式化した出力調整量と、最適化計算を再実行して求めた値（真の増分とも呼ぶ。）とを比較して、定式化した出力調整量の精度を確認した。また、実際にコストの変化分も計算して、比較した。
【０１６４】
本実施例における図１に示すエネルギーシステム１の構成機器の具体例を表１に示す。
【表１】
2006178626.tif 000144
【０１６５】
また、Ｙ＝ａＸ２＋ｂＸ＋ｃのように２次関数で表現される各エネルギー機器４の特性関数の係数の具体例を表２に示す。
【表２】
2006178626.tif 000145
【０１６６】
また、ある計画対象時間帯について設定した電気料金および燃料（ガス）料金の具体例を表３に示し、当該計画対象時間帯における予測需要の具体例を表４に示す。
【表３】
2006178626.tif 000146
【表４】
2006178626.tif 000147
【０１６７】
先に説明した非線形計画法に基づく最適化計算を行った結果の各エネルギー機器４の出力を表５に示す。
【表５】
2006178626.tif 000148
【０１６８】
＜事例１：冷熱需要が変化するケース＞
冷熱需要が運用計画策定時の予測値から１００USRT増えて２３６０USRTとなった場合の、修正案候補１〜４における真の増分と定式化した出力調整量とを比較する。
【０１６９】
＜事例１の修正案候補１＞
修正案候補１では、購買電力と電動冷凍機４Ｃｎで冷熱需要の変化に対応する。電動冷凍機４Ｃは１台であるので、電動冷凍機４ＣではＲＥ１の出力を調整して、冷熱需要の増加に対応する。ＲＥ１の出力調整量は１００USRTである。表６に計算結果を示す。詳細計算法では真の増分と一致した。
【表６】
2006178626.tif 000149
【０１７０】
＜事例１の修正案候補２＞
修正案候補２では、コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄｍで冷熱需要の変化に対応する。吸収式冷凍機４ＤはＲＳ１を用いた。ＲＳ１の出力を調整して、冷熱需要の増加に対応する。ＲＳ１の出力調整量は１００USRTである。コジェネレーション装置４ＡについてはＣＵ１とＣＵ３を用いた場合に分けて計算した。表７にＣＵ１を用いた場合の計算結果を示し、表８にＣＵ３を用いた場合の計算結果を示す。尚、ＣＵ１を用いた場合を事例１の修正案候補２−１とし、ＣＵ３を用いた場合を事例１の修正案候補２−２とする。詳細計算法ではほぼ真の増分と一致した。簡易計算法でも真の増分に近い値となった。
【表７】
2006178626.tif 000150
【表８】
2006178626.tif 000151
【０１７１】
＜事例１の修正案候補３＞
修正案候補３では、コジェネレーション装置４Ａｉと電動冷凍機４Ｃｎで冷熱需要の変化に対応する。電動冷凍機４Ｃは１台であるので、電動冷凍機４ＣではＲＥ１の出力を調整して、冷熱需要の増加に対応する。ＲＥ１の出力調整量は１００USRTである。コジェネレーション装置４ＡについてはＣＵ１とＣＵ３を用いた場合に分けて計算した。表９にＣＵ１を用いた場合の計算結果を示し、表１０にＣＵ３を用いた場合の計算結果を示す。尚、ＣＵ１を用いた場合を事例１の修正案候補３−１とし、ＣＵ３を用いた場合を事例１の修正案候補３−２とする。詳細計算法では真の増分と一致した。簡易計算法でも出力の調整量はほぼ一致した。
【表９】
2006178626.tif 000152
【表１０】
2006178626.tif 000153
【０１７２】
＜事例１の修正案候補４＞
修正案候補４では、コジェネレーション装置４Ａｉと複数台の吸収式冷凍機４Ｄと複数台の電動冷凍機４Ｃで冷熱需要の変化に対応する。吸収式冷凍機４Ｄは２台とも用い、電動冷凍機４Ｃは１台だけなので１台用いた。コジェネレーション装置４ＡについてはＣＵ１とＣＵ３を用いた場合に分けて計算した。表１１にＣＵ１を用いた場合の計算結果を示し、表１２にＣＵ３を用いた場合の計算結果を示す。尚、ＣＵ１を用いた場合を事例１の修正案候補４−１とし、ＣＵ３を用いた場合を事例１の修正案候補４−２とする。詳細計算法ではコストの変化に若干の違いがあるがほぼ真値に一致している。簡易計算法でも調整量はほぼ真値に一致している。
【表１１】
2006178626.tif 000154
【表１２】
2006178626.tif 000155
【０１７３】
＜事例２：電力需要が変化するケース＞
電力需要が運用計画策定時の予測値から１００ｋＷ増えて３２６０ｋＷとなった場合の、修正案候補１，２における真の増分と定式化した出力調整量とを比較する。
【０１７４】
＜事例２の修正案候補１＞
修正案候補１では、コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄと電動冷凍機４Ｃとで電力需要の変化に対応する。吸収式冷凍機４Ｄは２台とも用い、電動冷凍機４Ｃは１台だけなので１台用いた。コジェネレーション装置４ＡについてはＣＵ１とＣＵ３を用いた場合に分けて計算した。表１３にＣＵ１を用いた場合の計算結果を示し、表１４にＣＵ３を用いた場合の計算結果を示す。コストの変化量に若干の違いがあるがほぼ真値に一致している。
【表１３】
2006178626.tif 000156
【表１４】
2006178626.tif 000157
【０１７５】
＜事例２の修正案候補２＞
修正案候補２では、電力需要の変化分は購買電力で賄う。この場合は、購買電力が１００ｋＷ増えるだけなので、コストの増分は１００ｋＷ×１６円／ｋＷ＝１６００円となる。
【０１７６】
＜事例３：電気料金が変化するケース＞
電気料金が運用計画策定時の設定値から２円／ｋＷ増えて１８円／ｋＷとなった場合の、修正案候補における真の増分と定式化した出力調整量とを比較する。
【０１７７】
＜事例３の修正案候補＞
コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄｍと電動冷凍機４Ｃｎとで電気料金の変化に対応する。吸収式冷凍機４ＤはＲＳ１を用いた。コジェネレーション装置４ＡについてはＣＵ１とＣＵ３を用いた場合に分けて計算した。表１５にＣＵ１を用いた場合の計算結果を示し、表１６にＣＵ３を用いた場合の計算結果を示す。若干の違いがあるものの、どちらもほぼ真値に一致した。
【表１５】
2006178626.tif 000158
【表１６】
2006178626.tif 000159
【０１７８】
＜事例４：燃料料金が変化するケース＞
ガス料金が運用計画策定時の設定値から２円／Ｎｍ３増えて７２円／Ｎｍ３となった場合の、修正案候補における真の増分と定式化した出力調整量とを比較する。
【０１７９】
＜事例４の修正案候補＞
コジェネレーション装置４Ａｉと吸収式冷凍機４Ｄｍと電動冷凍機４Ｃｎとで電気料金の変化に対応する。吸収式冷凍機４ＤはＲＳ１を用い、電動冷凍機４Ｃは１台だけなのでＲＥ１を用いた。コジェネレーション装置４ＡについてはＣＵ１とＣＵ３を用いた場合に分けて計算した。表１７にＣＵ１を用いた場合の計算結果を示し、表１８にＣＵ３を用いた場合の計算結果を示す。若干の違いがあるものの、どちらもほぼ真値に一致した。
【表１７】
2006178626.tif 000160
【表１８】
2006178626.tif 000161
【０１８０】
電力需要、冷熱需、電気料金、ガス料金といった各種パラメータが運用計画策定時の値から変化した場合の各修正案候補における出力調整量とコストの変化量は、表６から表１８に示したとおりである。このように各修正案候補で出力調整量をまず算定し、それを使ってコストの変化量も算定し、該当事例に対応するそれぞれの修正案候補の中で最もトータルコストを安く出来るものを採用すればよい。例えば事例１に対応する各修正案候補のコストの変化量は、以下の表１９に示すものとなる。
【表１９】
2006178626.tif 000162
【０１８１】
この表１９より、簡易計算法で計算した場合には修正案候補１が選択され、詳細計算法で計算した場合には修正案候補４−１が選択されて、前日以前に策定したエネルギーシステム１の運用計画の修正案となる。なお、実際の現場でこの修正案を検討した結果、運用上の都合などの事情で実行できないとなれば、修正案を微調整して実際のエネルギーシステム１の仕様等に合わせるか、２番目、３番目にコストの安い修正案候補を検討する。
【０１８２】
以上のように本発明のエネルギーシステム運用計画の修正方法および装置およびプログラムによれば、運用計画の基となった予測に基づく各種パラメータ、例えば電力・冷熱・温熱の各需要や電気料金や燃料の料金が、運転当日等に判明した実際値とずれている場合でも、当該運転当日に、運用計画を再度策定し直すことなく、迅速に適切な調整対象機器の出力調整量を求めることができる。運用者の経験や勘に頼ることなく、例えばトータルコストを最小とするといった運用計画策定時の目標を損なうことなく、パラメータ変化に直ちに対応した適切な運用計画の修正が行える。
【０１８３】
なお、上述の実施形態は本発明の好適な実施の一例ではあるがこれに限定されるものではなく、本発明の要旨を逸脱しない範囲において種々変形実施可能である。例えば、本実施形態では、エネルギー機器４の出力調整量ΔＸを表す式としてテーラー展開の１次近似及び２次近似を用いているが、これに限られず、テーラー展開の３次以上の近似を用いても良い。
【図面の簡単な説明】
【０１８４】
【図１】エネルギーシステムの一例を示す構成図である。
【図２】本発明のエネルギーシステム運用計画の修正方法の手順の一例を示すフローチャートである。
【図３】本発明のエネルギーシステム運用計画の修正方法の手順の一例を示すフローチャートである。
【図４】本発明のエネルギーシステム運用計画の修正装置の一例を示す構成図である。
【符号の説明】
【０１８５】
１  エネルギーシステム
２  電力・温熱供給手段
３  冷熱供給手段
４  エネルギー機器
４Ａ  コジェネレーション装置
４Ｂ  ボイラー装置
４Ｃ  電動冷凍機
４Ｄ  吸収式冷凍機
５  熱交換器
６  冷却塔
７  運用計画策定手段
８  運用計画修正装置
【図１】
2006178626.tif 000163
【図２】
2006178626.tif 000164
【図３】
2006178626.tif 000165
【図４】
2006178626.tif 000166