(57)【要約】【課題】建物の空調負荷を推定演算する。【解決手段】空調負荷演算部(20)と割合演算部(30)と月毎の空調負荷演算部(30)の演算により演算した最大空調負荷（Ｑimax）及び最小空調負荷（Ｑimin）と、月毎の使用電力量演算部(50)で求めた最大使用電力量（Ｗｐ）及び最小使用電力量（Ｗｂ）とを基に、ベース電力量演算部(60)で、ベース電力量（Ｗｏ）を演算する。空調負荷用電力量演算部(70)で空調負荷のために使用する月毎の電力量（Ｗqi）を求め、最大空調負荷演算部(80)で最大空調負荷（Ｑmax）を求め、空調負荷演算部(90)で各負荷推定時刻での空調負荷（Ｑij）を演算する。【選択図】図１
2011202840.tif 000002
【特許請求の範囲】
【請求項１】
空調機が設置される建物の空調負荷を推定するため、１月から１２月までの月毎に代表日を設定し、各代表日には予め決めた時間毎に負荷推定時刻を設定し、各負荷推定時刻における前記建物の空調負荷を推定する空調負荷推定システムであって、
各負荷推定時刻の気象状況と前記建物に関する物理特性に基づき、前記建物の窓での各負荷推定時刻の空調負荷（Ｑｇ）と、前記建物の壁での各負荷推定時刻の空調負荷（Ｑｗ）と、前記建物の屋根での各負荷推定時刻の空調負荷（Ｑｒ）を演算する空調負荷演算部と、
前記空調機が使用されない時間帯における各負荷推定時刻の前記空調負荷（Ｑｇ）と前記空調負荷（Ｑｗ）と前記空調負荷（Ｑｒ）を強制的に零にしてから前記空調負荷（Ｑｇ）と前記空調負荷（Ｑｗ）と前記空調負荷（Ｑｒ）を加算して各負荷推定時刻の空調負荷(Ｑij)を演算し、前記空調負荷(Ｑij)のうち最大のものを最大空調負荷（Ｑmax）とし、各負荷推定時刻の空調負荷(Ｑij)を最大空調負荷（Ｑmax）で除算することにより各負荷推定時刻の割合（ｒij）を演算する割合演算部と、
前記空調機を稼動する月毎の稼動日数（ｎi）と各負荷推定時刻の空調負荷(Ｑij)を基に、月毎の空調負荷（Ｑi）を演算し、月毎の空調負荷（Ｑi）のうち最大のものを最大空調負荷（Ｑimax）とし最小のものを最小空調負荷（Ｑimin）として選択する月毎の空調負荷演算部と、
月毎の電気料金と料金体系、または、電気料金の明細に示された使用電力量から、月毎の使用電力量（Ｗi）を求め、月毎の使用電力量（Ｗi）のうち最大のものを最大使用電力量（Ｗｐ）とし最小のものを最小使用電力量（Ｗｂ）とする月毎の使用電力量演算部と、
最大空調負荷（Ｑimax）及び最小空調負荷（Ｑimin）と最大使用電力量（Ｗｐ）及び最小使用電力量（Ｗｂ）と空調負荷に関係のないベース電力量（Ｗｏ）が、
（Ｗｂ−Ｗｏ）／（Ｗｐ−Ｗｏ）＝Ｑimin／Ｑimaxの関係にあることを利用してベース電力量（Ｗｏ）を求めるベース電力量演算部と、
月毎の使用電力量（Ｗi）からベース電力量（Ｗｏ）を減算することにより、空調負荷のために使用する月毎の電力量（Ｗqi）を演算する空調負荷用電力量演算部と、
月毎の電力量（Ｗqi）のうち使用電力量が最大になる月の電力量である最大使用電力量（Ｗｐ）と、ベース電力量（Ｗｏ）と、前記空調機の成績係数（ｐ）と、使用電力量が最大になる月の空調機の稼動日数（ｎi）と、使用電力量が最大になる月の各負荷推定時刻の割合（ｒij）を用いて最大空調負荷（Ｑmax）を演算する最大空調負荷演算部と、
前記最大空調負荷演算部で演算した最大空調負荷（Ｑmax）に、各負荷推定時刻の割合（ｒij）を乗算することにより、各負荷推定時刻での空調負荷（Ｑij）を演算する各負荷推定時刻の空調負荷演算部と、
で構成されていることを特徴とする空調負荷推定システム。
【請求項２】
前記空調負荷演算部は、
各負荷推定時刻の日射量ｑと、室外空気の熱伝達率αａと、室外壁表面と室内空気との間の熱抵抗Ｒｗと、室外屋根表面と室内空気との間の熱抵抗Ｒｒと、窓の日射透過率Ｓｇと、壁の表面塗料の熱吸収率εｗと、屋根の表面塗料の熱吸収率εｒと、各負荷推定時刻の外気温度Ｔａと、各負荷推定時刻の室内設定温度Ｔsetと、外気と室内空気との間の熱抵抗Ｒｇと、窓の面積Ａｇと、壁の面積Ａｗと、屋根の面積Ａｒを、下式（１０１），（１０２），（１０３）に適用して、前記空調負荷（Ｑｇ）と、前記空調負荷（Ｑｗ）と、前記空調負荷（Ｑｒ）を演算することを特徴とする請求項１の空調負荷推定システム。
【数１】
2011202840.tif 000009
【請求項３】
前記月毎の空調負荷演算部は、下式（１０４）を用いて、月毎の空調負荷（Ｑi）を演算することを特徴とする空調負荷推定システム。
【数２】
2011202840.tif 000010
【請求項４】
前記最大空調負荷演算部は、下式（１０７）を用いて、最大空調負荷（Ｑimax）を演算することを特徴とする空調負荷推定システム。
【数３】
2011202840.tif 000011
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は空調負荷推定システムに関し、簡易かつ適切に建物の空調負荷を推定することができるように工夫したものである。
【背景技術】
【０００２】
ターボ冷凍機などの大形の冷凍機は、ビルや工場など広い面積を有する建物に設置されている。
冷凍機の製造・販売者（以下「製造者」と称する）が、冷凍機を購入して建物に設置する者（以下「購入者」と称する）に対して、冷凍機を納入・販売をする営業活動をする場合には、最適な機種を選定する必要がある。
このような冷凍機の営業にあたっては、従来では、製造者は、購入者から提示される仕様（主として定格点の仕様（定格容量、定格冷水設定温度、定格冷水流量、定格冷却水温度など））に基づいて、冷凍機の機種を選定する場合が多かった。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００５−２３３４３８号
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、実際に冷凍機が使用される状況は定格点のみの運転ではなく、幅広い負荷に対応した運転を行うこととなる。つまり、定格点（定格負荷）で運転されたり、定格負荷よりも小さい負荷（部分負荷）で運転されたりする。
したがって、定格点を基に選定した機種が、最適な機種であるとは言えない場合もあった。
【０００５】
例えば負荷の定格が５００ｋＷであった場合に、定格容量が５００ｋＷの冷凍機を設置する場合について具体的に説明する。
冷凍機の設置場所（例えば工場）において、昼間でも夜間でも負荷がほぼ定格である場合には、定格容量が５００ｋＷの１台の冷凍機を設置するのが適切である。
一方、冷凍機の設置場所（例えばオフィスビル）において、昼間においては負荷が定格であるが、夜間においては負荷が半減する場合には、定格容量が２５０ｋＷの２台の冷凍機を設置し、昼間には２台の冷凍機を使用し、夜間では１台を休止し残りの１台のみを使用することが、適切である。
更に冷凍機によっては、定格負荷で運転した場合に効率が良い特性を持った機種（例えば回転速度が固定速度の機種）がある一方で、定格負荷のみならず部分負荷で運転した場合であっても効率が良い特性を持った機種（例えばインバータを使用した可変回転速度の機種）がある。
【０００６】
このように空調負荷の変動特性等によって、最適な機種（定格容量、台数、機種特性）が異なるにもかかわらず、従来では最適な機種の冷凍機を選定することができない恐れもあった。
機種選定が最適でない場合には、要求された空調能力を満たすことはできたとしても、年間に使用する総電力量が不要に増加してしまうことがあった。
【０００７】
購入者における年間の空調負荷を推定することができれば、上述したような不具合を防止することができるため、空調負荷を推定（算定）するシステムの発明がすでに幾つか提案されている（例えば特許文献１参照）。
しかし、簡易かつ適切に、空調負荷を推定するシステムが無いのが現状であった。
【０００８】
本発明は、上記従来技術に鑑み、空調機（冷凍機）が設置される建物における使用電力量を演算要件として用いて、簡易かつ適切に、空調負荷を推定することができる、空調負荷推定システムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決する本発明の構成は、
空調機が設置される建物の空調負荷を推定するため、１月から１２月までの月毎に代表日を設定し、各代表日には予め決めた時間毎に負荷推定時刻を設定し、各負荷推定時刻における前記建物の空調負荷を推定する空調負荷推定システムであって、
各負荷推定時刻の気象状況と前記建物に関する物理特性に基づき、前記建物の窓での各負荷推定時刻の空調負荷（Ｑｇ）と、前記建物の壁での各負荷推定時刻の空調負荷（Ｑｗ）と、前記建物の屋根での各負荷推定時刻の空調負荷（Ｑｒ）を演算する空調負荷演算部と、
前記空調機が使用されない時間帯における各負荷推定時刻の前記空調負荷（Ｑｇ）と前記空調負荷（Ｑｗ）と前記空調負荷（Ｑｒ）を強制的に零にしてから前記空調負荷（Ｑｇ）と前記空調負荷（Ｑｗ）と前記空調負荷（Ｑｒ）を加算して各負荷推定時刻の空調負荷(Ｑij)を演算し、前記空調負荷(Ｑij)のうち最大のものを最大空調負荷（Ｑmax）とし、各負荷推定時刻の空調負荷(Ｑij)を最大空調負荷（Ｑmax）で除算することにより各負荷推定時刻の割合（ｒij）を演算する割合演算部と、
前記空調機を稼動する月毎の稼動日数（ｎi）と各負荷推定時刻の空調負荷(Ｑij)を基に、月毎の空調負荷（Ｑi）を演算し、月毎の空調負荷（Ｑi）のうち最大のものを最大空調負荷（Ｑimax）とし最小のものを最小空調負荷（Ｑimin）として選択する月毎の空調負荷演算部と、
月毎の電気料金と料金体系、または、電気料金の明細に示された使用電力量から、月毎の使用電力量（Ｗi）を求め、月毎の使用電力量（Ｗi）のうち最大のものを最大使用電力量（Ｗｐ）とし最小のものを最小使用電力量（Ｗｂ）とする月毎の使用電力量演算部と、
最大空調負荷（Ｑimax）及び最小空調負荷（Ｑimin）と最大使用電力量（Ｗｐ）及び最小使用電力量（Ｗｂ）と空調負荷に関係のないベース電力量（Ｗｏ）が、
（Ｗｂ−Ｗｏ）／（Ｗｐ−Ｗｏ）＝Ｑimin／Ｑimaxの関係にあることを利用してベース電力量（Ｗｏ）を求めるベース電力量演算部と、
月毎の使用電力量（Ｗi）からベース電力量（Ｗｏ）を減算することにより、空調負荷のために使用する月毎の電力量（Ｗqi）を演算する空調負荷用電力量演算部と、
月毎の電力量（Ｗqi）のうち使用電力量が最大になる月の電力量である最大使用電力量（Ｗｐ）と、ベース電力量（Ｗｏ）と、前記空調機の成績係数（ｐ）と、使用電力量が最大になる月の空調機の稼動日数（ｎi）と、使用電力量が最大になる月の各負荷推定時刻の割合（ｒij）を用いて最大空調負荷（Ｑmax）を演算する最大空調負荷演算部と、
前記最大空調負荷演算部で演算した最大空調負荷（Ｑmax）に、各負荷推定時刻の割合（ｒij）を乗算することにより、各負荷推定時刻での空調負荷（Ｑij）を演算する各負荷推定時刻の空調負荷演算部と、
で構成されていることを特徴とする。
【００１０】
また本発明の構成は、
前記空調負荷演算部は、
各負荷推定時刻の日射量ｑと、室外空気の熱伝達率αａと、室外壁表面と室内空気との間の熱抵抗Ｒｗと、室外屋根表面と室内空気との間の熱抵抗Ｒｒと、窓の日射透過率Ｓｇと、壁の表面塗料の熱吸収率εｗと、屋根の表面塗料の熱吸収率εｒと、各負荷推定時刻の外気温度Ｔａと、各負荷推定時刻の室内設定温度Ｔsetと、外気と室内空気との間の熱抵抗Ｒｇと、窓の面積Ａｇと、壁の面積Ａｗと、屋根の面積Ａｒを、下式（１０１），（１０２），（１０３）に適用して、前記空調負荷（Ｑｇ）と、前記空調負荷（Ｑｗ）と、前記空調負荷（Ｑｒ）を演算することを特徴とする。
【数１】
2011202840.tif 000003
【００１１】
また本発明の構成は、
前記月毎の空調負荷演算部は、下式（１０４）を用いて、月毎の空調負荷（Ｑi）を演算することを特徴とする。
【数２】
2011202840.tif 000004
【００１２】
また本発明の構成は、
前記最大空調負荷演算部は、下式（１０７）を用いて、最大空調負荷（Ｑimax）を演算することを特徴とする。
【数３】
2011202840.tif 000005
【発明の効果】
【００１３】
本発明によれば、空調機（冷凍機）が設置される建物における使用電力量を演算要件として用いて、簡易かつ適切に、空調負荷を推定することができる。
これにより最適機種の選択をして空調機（冷凍機）を営業する際の、営業支援をすることができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施例に係る空調負荷推定システムを示す機能ブロック図。
【図２】実施例の演算処理動作を示すフローチャート。
【発明を実施するための形態】
【００１５】
以下、本発明の実施の形態について、実施例に基づき詳細に説明する。
【実施例１】
【００１６】
図１は本発明の実施例１に係る空調負荷推定システム１を示す。この空調負荷推定システム１は、空調機（冷凍機）が設置される建物の空調負荷を、コンピュータシステム演算処理によって推定する空調負荷推定システムであって、コンピュータシステムによって実現される。
具体的には、空調負荷推定システム１は、入力装置１０と、空調負荷演算部２０と、割合演算部３０と、月毎の空調負荷演算部４０と、月毎の使用電力量演算部５０と、ベース電力量演算部６０と、空調負荷用電力量演算部７０と、最大空調負荷演算部８０と、各負荷推定時刻の空調負荷演算部９０と、出力装置１００を備えて構成されている。
なお、空調負荷推定システム１の上記各演算部２０〜９０は、空調負荷推定システム１を構成するコンピュータのハードウェアとそのハードウェア上で実行されるアプリケーションソフトウェアで構成されている。
【００１７】
ここで、空調負荷推定システム１で実行する演算処理において用いる、概念や前提を先に説明しておく。
【００１８】
（１）各月の「代表日」について
本例では、１月から１２月までの各月において、代表日を設定している。
例えば１月の代表日では、代表日の一日の気象変化（温度変化、日照変化等）が、１月の全ての日における気象変化（温度変化、日照変化等）を平均したものとなるように設定している。
同様に、他の２月〜１２月の代表日では、２月〜１２月の各月の代表日の気象変化が、該当する月の全ての日における気象変化を平均したものとなるように設定している。
各月（１月〜１２月）の各日における気象変化（温度変化、日照変化等）は、公的気象データとして公的機関（気象庁等）から発表されているので、この公的気象データを用いる。
【００１９】
（２）「負荷推定時刻」について
１月から１２月までの各月に設定した代表日には、それぞれ、０時から２３時まで１時間毎（予め決めた時間毎）に負荷推定時刻を設定している。
したがって、１年間では、１２カ月分×２４負荷推定時刻（＝２８８負荷推定時刻）を、負荷推定時刻として設定している。
【００２０】
（３）「各負荷推定時刻の日射量」について
各負荷推定時刻の日射量としては、上記の各負荷推定時刻（＝２８８負荷推定時刻）において、空調機（冷凍機）が設置される建物の大きさ及び方位を考慮したこの建物に照射される日射量を設定している。
各月（１月〜１２月）の各日の１時間毎の日射量は、公的気象データとして公的機関（気象庁等）から発表されているので、この公的気象データを基に、上記の各負荷推定時刻の日射量を予め算出している。
例えば、１月の代表日の負荷推定時刻１２時の日射量は、１月の全ての日の１２時の日射量（空調機（冷凍機）が設置される建物の大きさ及び方位を考慮した日射量）を平均した値に設定している。
【００２１】
（４）「各負荷推定時刻の外気温度」について
各負荷推定時刻の外気温度としては、上記の各負荷推定時刻（＝２８８負荷推定時刻）において、空調機が設置される建物の周囲の温度を設定している。
各月（１月〜１２月）の各日の１時間毎の外気温度は、公的気象データとして公的機関（気象庁等）から発表されているので、この公的気象データを基に、上記の各負荷推定時刻の外気温度を予め算出している。
例えば、１月の代表日の負荷推定時刻１２時の外気温度は、１月の全ての日の１２時の外気温度（空調機が設置される建物の周囲の温度）を平均した値に設定している。
【００２２】
（５）空調機が設置される建物において、使用電力量が最大になる月と月毎の空調負荷が最大になる月は合致し、使用電力量が最小になる月と月毎の空調負荷が最小になる月は合致することを前提とする。
本発明は、空調機が建物に既に設置されその空調機が運用されている場合に、その建物の空調負荷を正確に推定するものである。そして、既存の空調機を新規の空調機に取り替える際に、その建物について正確に推定した空調負荷をもとに、最適な機種の空調機（冷凍機）を選択することができるようになる。
これにより最適機種の選択をして空調機（冷凍機）を営業する際の、営業支援をすることができる。
【００２３】
［空調負荷推定の演算処理］
次に、かかる空調負荷推定システム１により、空調機を設置しようとする建物の空調負荷を推定する演算処理手順を、図１の機能ブロック図及び図２のフローチャートに基づき説明する。
【００２４】
［各負荷推定時刻における、日射による単位面積当たりの空調負荷の演算］
入力装置１０により、各負荷推定時刻の日射量ｑと、室外空気の熱伝達率αａと、室外壁表面と室内空気との間の熱抵抗Ｒｗと、室外屋根表面と室内空気との間の熱抵抗Ｒｒと、窓の日射透過率Ｓｇと、壁の表面塗料の熱吸収率εｗと、屋根の表面塗料の熱吸収率εｒを入力すると、空調負荷演算部２０は、各負荷推定時刻における、日射による単位面積当たりの空調負荷の演算をする（ステップ１）。
【００２５】
具体的には、建物の窓、壁、屋根に分けて、各負荷推定時刻における、日射による単位面積当たりの空調負荷を次のように演算する。
（１）窓での、各負荷推定時刻における、日射による単位面積当たりの空調負荷は、各負荷推定時刻の日射量がｑ、窓の日射透過率がＳｇであるので、
Ｓｇ×ｑ
の演算をして求める。
（２）壁での、各負荷推定時刻における、日射による単位面積当たりの空調負荷は、日射量がｑ、壁の表面塗料の熱吸収率がεｗであるので、
（εｗ×ｑ）／（１＋αａ・Ｒｗ）
の演算をして求める。
（３）屋根での、各負荷推定時刻における、日射による単位面積当たりの空調負荷は、日射量がｑ、屋根の表面塗料の熱吸収率がεｒであるので、
（εｒ×ｑ）／（１＋αａ・Ｒｒ）
の演算をして求める。
【００２６】
［各負荷推定時刻における、温度差による単位面積当たりの空調負荷の演算］
入力装置１０により、各負荷推定時刻の外気温度Ｔａと、各負荷推定時刻の室内設定温度Ｔsetと、外気と室内空気との間の熱抵抗Ｒｇを更に入力すると、空調負荷演算部２０は、各負荷推定時刻における、温度差による単位面積当たりの空調負荷の演算をする（ステップ２）。
【００２７】
具体的には、建物の窓、壁、屋根に分けて、各負荷推定時刻における、温度差による単位面積当たりの空調負荷を次のように演算する。
（１）窓での、各負荷推定時刻における、温度差による単位面積当たりの空調負荷は、
（Ｔａ−Ｔset）／Ｒｇ
の演算をして求める。
（２）壁での、各負荷推定時刻における、温度差による単位面積当たりの空調負荷は、
αａ・（Ｔａ−Ｔset）／（１＋αａ・Ｒｗ）
の演算をして求める。
（３）屋根での、各負荷推定時刻における、温度差による単位面積当たりの空調負荷は、
αａ・（Ｔａ−Ｔset）／（１＋αａ・Ｒｒ）
の演算をして求める。
【００２８】
［各負荷推定時刻における、空調負荷の演算］
入力装置１０により、窓の面積Ａｇと、壁の面積Ａｗと、屋根の面積Ａｒを更に入力すると、空調負荷演算部２０は、各負荷推定時刻における、空調負荷の演算をする（ステップ３）。
【００２９】
具体的には、建物の窓、壁、屋根に分けて、各負荷推定時刻における、空調負荷を次のように演算する。
（１）窓での、各負荷推定時刻における、空調負荷Ｑｇは、下式（１０１）の演算をして求める。
式（１０１）に示すように、窓での、各負荷推定時刻における、空調負荷Ｑｇは、各負荷推定時刻における、「日射による単位当たりの空調負荷」と「温度差による単位面積当たりの空調負荷」とを加えたものである。
（２）壁での、各負荷推定時刻における、空調負荷Ｑｗは、下式（１０２）の演算をして求める。
式（１０２）に示すように、壁での、各負荷推定時刻における、空調負荷Ｑｗは、各負荷推定時刻における、「日射による単位当たりの空調負荷」と「温度差による単位面積当たりの空調負荷」とを加えたものである。
（３）屋根での、各負荷推定時刻における、空調負荷Ｑｒは、下式（１０３）の演算をして求める。
式（１０３）に示すように、屋根での、各負荷推定時刻における、空調負荷Ｑｒは、各負荷推定時刻における、「日射による単位当たりの空調負荷」と「温度差による単位面積当たりの空調負荷」とを加えたものである。
【数４】
2011202840.tif 000006
【００３０】
［各負荷推定時刻における、空調負荷の割合の演算］
割合演算部３０は、後に示す演算処理をして、各負荷推定時刻における、各負荷推定時刻の空調負荷と、各空調負荷のうち最大の空調負荷との割合を演算する（ステップ４）。
【００３１】
割合演算部３０では、1月から１２月までの各月を一般的にｉで表し、代表日の０時から２３時までの各負荷推定時刻をｊで表す。
そうすると式（１０１）により求めた窓での各負荷推定時刻における空調負荷Ｑｇは、Ｑgijと表される。式（１０２）により求めた壁での各負荷推定時刻における空調負荷Ｑｗは、Ｑwjiと表される。式（１０３）により求めた屋根での各負荷推定時刻における空調負荷Ｑｒは、Ｑrijと表される。
【００３２】
割合演算部３０は、空調負荷演算部２０から各負荷推定時刻における各部（窓、壁、屋根）の空調負荷Ｑgij、Ｑwji、Ｑrijを取り込み、空調機が運転していない負荷推定時刻では、その負荷推定時刻における各部（窓、壁、屋根）の空調負荷Ｑgij、Ｑwji、Ｑrijの値を強制的に零にする演算処理をする。
つまり、昼夜を問わず空調が必要な工場等の建物では、昼夜を問わず空調機が運転されているため、一日じゅう各負荷推定時刻の各部の空調負荷gij、Ｑwji、Ｑrijは零にはしないが、オフィスのように昼間は空調機が運転されているが夜間では空調機が停止する建物では、夜間における各負荷推定時刻の各部の空調負荷gij、Ｑwji、Ｑrijを強制的に零にする。
空調機が運転されていない負荷推定時刻は、入力装置１０により入力される。空調機が運転されていない負荷推定時間は、個別に入力してもよいし、工場型のパターンとオフィスビル型のパターンに分けて、工場型のパターンでは空調負荷gij、Ｑwji、Ｑrijは零にはしないが、オフィスビル型のパターンでは予め決めた夜間の時間帯において空調負荷gij、Ｑwji、Ｑrijを強制的に零にするようにしてもよい。
【００３３】
その後、割合演算部３０は、各負荷推定時刻の空調負荷Ｑijを、次式により演算する。つまり、各負荷推定時刻の空調負荷Ｑijは、各負荷推定時刻の各部（窓、壁、屋根）の空調負荷Ｑgij、Ｑwji、Ｑrijを加えたものである。なお前記処理により、空調負荷Ｑgij、Ｑwji、Ｑrijの中には、その値が強制的に零になっているものがある。
Ｑij＝Ｑgij＋Ｑwji＋Ｑrij
【００３４】
割合演算部３０は、上述した各空調負荷Ｑijのうち、空調負荷が最大のものを最大空調負荷Ｑmaxとする。例えば、８月の代表日の負荷推定時刻１４時（ｉ＝８月、ｊ＝１４時）における空調負荷が、最大空調負荷Ｑmaxとなる。
そして、各負荷推定時刻（＝２８８負荷推定時刻）における、各負荷推定時刻の空調負荷Ｑijと各空調負荷のうち最大の空調負荷Ｑmaxとの割合ｒijを次式により求める。この割合ｒijは各負荷推定時刻（＝２８８負荷推定時刻）についてそれぞれ求められる。
ｒij＝Ｑij／Ｑmax
【００３５】
［月毎の空調負荷のうち最大空調負荷と最小空調負荷の選択］
入力装置１０により、各月（ｉ）ごとの稼動日数ｎiを入力すると、月毎の空調負荷演算部４０は、月毎の空調負荷のうち、最大の空調負荷と最小の空調負荷を選択する（ステップ５）。
具体的には、月毎の空調負荷演算部４０は、入力された稼動日数ｎiと、割合演算部３０により演算した各負荷推定時刻の空調負荷Ｑijを次式（１０４）に適用して、月毎の空調負荷Ｑiを求める。つまり、１月から１２月（ｉ＝１月〜１２月）までの各月の空調負荷Ｑiを求める。
【数５】
2011202840.tif 000007
更に、月毎の空調負荷演算部４０は、上述した月毎の空調負荷Ｑiのうち、最大のものを最大空調負荷Ｑimax、最小のものを最小空調負荷Ｑiminとして選択する。
【００３６】
［月毎の使用電力量のうち最大使用電力量と最小使用電力量の選択］
入力装置１０により、１年間の月毎（１月から１２月まで）の電気料金と料金体系を入力すると、月毎の使用電力量演算部５０は、月毎の使用電力量を演算し、最大使用電力量と最小使用電力量を選択する（ステップ６）。
具体的には、月毎の使用電力量演算部５０は、月毎の電気料金と料金体系を基に、月ごとの使用電力量Ｗiを求める（なお電気料金の明細で使用電力量が分かる場合はそれを用いる）。
このとき、使用電力量Ｗiのうち、最大のものを最大使用電力量Ｗｐ、最小のものを最小使用電力量Ｗｂとする(このときの最大使用電力量、最小使用電力量の発生月は、空調負荷の最大空調負荷、最小空調負荷の発生月と合致しているものとする)。
【００３７】
［空調負荷に関係のないベース電力量の演算］
ベース電力量演算部６０は、空調負荷に関係のない１カ月当たりのベース電力量Ｗｏの演算をする（ステップ７）。
具体的には、ベース電力量演算部６０は、月毎の空調負荷演算部４０により選択した最大空調負荷Ｑimaxと最小空調負荷をＱiminを取り込むと共に、月毎の使用電力量演算部５０により選択した最大使用電力量Ｗｐと最小電力量Ｗｂを取り込む。
そして、ベース電力量Ｗｏと、取り込んだ上記データとが、次式の関係にあることを利用して、空調負荷に関係のないベース電力量Ｗｏを演算する。
（Ｗｂ−Ｗｏ）／（Ｗｐ−Ｗｏ）＝Ｑimin／Ｑimax
【００３８】
［月毎の使用電力量のうち、空調負荷のために使用する電力量の演算］
空調負荷用電力量演算部７０は、月毎の使用電力量Ｗiのうち、空調負荷のために使用する月毎の電力量Ｗqiを演算する（ステップ８）
具体的には、空調負荷用電力量演算部７０は、ベース電力量演算部６０により演算したベース電力量Ｗｏと、月毎の使用電力量演算部５０により演算した月毎の使用電力量Ｗｉを取り込んで、次式に適用する。これにより、月毎の使用電力量Ｗiのうち、空調負荷のために使用する月毎の電力量Ｗqiを演算することができる。
Ｗqi＝Ｗi−Ｗo
【００３９】
［最大空調負荷を、電力量と成績係数と割合と稼動日数を演算要素として求める演算］
最大空調負荷演算部８０は、空調機（熱源機、例えばターボ冷凍機）のＣＯＰ（成績係数）をｐとして、最大使用電力量Ｗｐとなる月（例えばｉ＝８月）で考えて、最大空調負荷Ｑmaxを、電力量Ｗｐ，Ｗｏと成績係数ｐと割合ｒijと稼動日数ｎｉを演算要素として求める（ステップ９）。
【００４０】
具体的には、最大空調負荷演算部８０は、入力装置１０から成績係数ｐと稼動日数ｎｉが入力され、割合演算部３０により演算した割合ｒijを取り込み、月毎の使用電力量演算部５０により演算した最大使用電力量Ｗｐを取り込み、ベース電力演算部６０により演算したベース電力Ｗｏを取り込み、空調負荷用電力量演算部７０により演算した空調負荷のために使用する月毎の電力量Ｗqiを取り込む。
【００４１】
空調負荷のために使用する月毎の電力量Ｗqiが最大となる月例えばｉ＝８月）、つまり最大使用電力量Ｗｐとなる月（例えばｉ＝８月）で考えると、下式（１０５）が成立する。
式（１０５）の左辺は、最大使用電力量Ｗｐとなる月（例えばｉ＝８月）における空調負荷Ｑｉ（例えばｉ＝８月）であり、右辺は、最大使用電力量Ｗｐとなる月（例えばｉ＝８月）における電力量Ｗqi（例えばｉ＝８月）に成績係数ｐを掛けたものである。
式（１０５）を変形すると下式（１０６）となる。式（１０６）において、ｒijは最大使用電力量Ｗｐとなる月（例えばｉ＝８月）の割合であり、ｎｉは最大使用電力量Ｗｐとなる月（例えばｉ＝８月）の稼動日数である。
更に、式（１０６）を整理すると、式（１０７）となる。
【数６】
2011202840.tif 000008
【００４２】
最大空調負荷演算部８０は、式（１０７）の演算をして、最大空調負荷Ｑmaxを、電力量Ｗｐ，Ｗｏと、成績係数ｐと、使用電力量が最大となる月（例えばｉ＝８月）の割合ｒijと、使用電力量が最大となる月（例えばｉ＝８月）の稼動日数ｎｉを演算要素として求める。
【００４３】
［各負荷推定時刻での空調負荷を求める演算］
各負荷推定時刻の空調負荷演算部９０は、最大空調負荷Ｑmaxと各負荷推定時刻（＝２８８負荷推定時刻）での空調負荷の割合ｒijを用いて、各負荷推定時刻（＝２８８負荷推定時刻）の空調負荷を求める(ステップ１０)。
具体的には各負荷推定時刻の空調負荷演算部９０は、最大空調負荷演算部８０から最大空調負荷Ｑmaxを取り込み、割合演算部３０から各負荷推定時刻の割合ｒijを取り込む。そして、最大空調負荷Ｑmaxに各負荷推定時刻の割合ｒijを乗算した値を、各負荷推定時刻（＝２８８負荷推定時刻）の空調負荷Ｑjiの値として求める。
例えば１月の代表日の負荷推定時刻１２時の空調負荷Ｑi(=1) j(=12)は、最大空調負Ｑmaxに、ｉ＝１月、ｊ＝１２時の割合ｒi(=1)  j(=12)を乗算して求めることができる。
【００４４】
このようにして求めた各負荷推定時刻（＝２８８負荷推定時刻）の空調負荷Ｑjiは、出力装置１００を介して出力されて、表示装置等に表示される。
【００４５】
なお、割合演算部３０においても、各負荷推定時刻の空調負荷Ｑjiを求めているが、割合演算部３０で求めた各負荷推定時刻の空調負荷Ｑjiは、式（１０１），（１０２），（１０３）を基に演算したものであり、空調負荷Ｑjiの増減の傾向は現実の空調負荷の傾向に沿ったものであるが、その値は、現実の空調負荷からシフトしたものとなっている可能性がある。
これは、式（１０１），（１０２），（１０３）が、現実の空調負荷を求めるために設定したものであるが、式と現実とが完全に一致しないことがあるからである。
【００４６】
これに対して、各負荷推定時刻の空調負荷演算部９０は、使用電力と成績係数を基に最大空調負荷演算部８０で求めた最大空調負荷Ｑmaxに対して、各負荷推定時刻の割合ｒijを乗算して、各負荷推定時刻の空調負荷Ｑjiの値を求めている。このため、各負荷推定時刻の空調負荷演算部９０で求めた各負荷推定時刻の空調負荷Ｑjiは、現実の空調負荷に近い正確なものとなっている。
【実施例２】
【００４７】
次に空調負荷推定システム１の変形例を、実施例２として説明する。実施例２は、基本的には実施例１と同様であるが、演算処理が一部異なる。
このため、以下の実施例２の説明では、実施例１と異なる部分を説明し、重複する部分の説明は省略する。
【００４８】
実施例２では、ステップ２において、外気と室内空気との間の熱抵抗Ｒｇとして、ガラスの材質を指定せずに、負荷として厳しい側（最小値）と楽な側（最大値）を与える。このため、窓における各負荷推定時刻の温度差による空調負荷を、負荷として厳しい側と楽な側の２つを演算する。
【００４９】
このためステップ４では、各負荷推定時刻の空調負荷Ｑijは、負荷として厳しい側と楽な側の２つを演算する。
【００５０】
この結果、ステップ５では、最大空調負荷Ｑimaxとして、負荷として厳しい側と楽な側の２つを演算し、最小空調負荷Ｑiminとして、負荷として厳しい側と楽な側の２つを演算する。
【００５１】
したがって、ステップ７では、空調負荷に関係のないベース電力量Ｗｏとして、負荷として厳しい側と楽な側の２つを演算する。
【００５２】
またステップ８では、空調負荷のために使用する月毎の電力量Ｗqiとして、負荷として厳しい側と楽な側の２つを演算する。
【００５３】
更にステップ９では、式（１０５）に用いる成績係数ｐを、一定値ではなく、各負荷推定時刻の負荷見合いのＣＯＰの値の平均値とする。
そして、最大空調負荷Ｑmaxとして、各月ごとに負荷として厳しい側と楽な側の２つを演算する。そして全ての月の、負荷として厳しい側と楽な側の１４の最大空調負荷Ｑmaxの偏差が最も小さくなる値を、最終的な最大空調負荷Ｑmaxとする。
【００５４】
ステップ１０では、ステップ９で求めた最終的な最大空調負荷Ｑmaxを用いて、各負荷推定時刻（＝２８８負荷推定時刻）の空調負荷を求める。
【符号の説明】
【００５５】
１  空調負荷推定システム
１０  入力装置
２０  空調負荷演算部
３０  割合演算部
４０  月毎の空調負荷演算部
５０  月毎の使用電力量演算部
６０  ベース電力量演算部
７０  空調負荷用電力量演算部
８０  最大空調負荷演算部
９０  各負荷推定時刻の空調負荷演算部
１００  出力装置
【図１】
2011202840.tif 000012
【図２】
2011202840.tif 000013