(57)【要約】【課題】インターバルメータやリアルタイムメータ等の高価な測定器を各需要家に設置することなく、需要家毎の電力消費量を容易に推定する。【解決手段】  本発明では、各需要家の時間毎の消費電力量を推定する。本発明は、複数の需要家を電力量の消費のパターンが類似するグループに区分しておき、資源の消費パターンを表した標準負荷カーブを、区分したグループ毎に生成する。本発明では、推定対象となる需要家が属するグループを特定し、特定したグループの標準負荷カーブをその推定対象となる需要家の電力消費規模により補正し、推定対象の需要家に対する任意の期間内の単位時間毎の資源消費量を推定する。【選択図】図４
2006011715.tif 000002
【特許請求の範囲】
【請求項１】
供給者から需要家へネットワークを介して供給される資源の消費量を推定する資源消費量の推定方法において、
複数の需要家を資源量の消費のパターンが類似するグループに区分しておき、
そのグループの資源の消費パターンを表した標準負荷カーブを、区分したグループ毎に生成しておき、
推定対象の需要家が属するグループを特定し、特定したグループの標準負荷カーブを用いて、推定対象の需要家に対する任意の期間内の単位時間毎の資源消費量を推定すること
を特徴とする資源消費量の推定方法。
【請求項２】
グループ毎の標準負荷カーブを生成する際に、そのグループに属するサンプルとなる需要家を選択し、選択した需要家の消費のパターンを実測し、実測した各需要家の負荷のパターンに基づき、上記標準負荷カーブを生成すること
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項３】
上記標準負荷カーブを、そのグループに含まれている需要家の属性を表すインデックスとともに格納しておくこと
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項４】
上記ネットワークの系統単位での資源の消費量を実測し、実測した値に基づき生成したパターンを、その系統の下に含まれている需要家のグループの標準負荷パターンとすること
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項５】
任意のグループに属するサンプルとなる１以上の需要家の単位時間毎の資源の消費量を実測し、実測した値に基づき生成されたパターンを、そのグループの標準負荷カーブとすること
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項６】
運転パターンが推定可能な需要家に対しては、標準負荷カーブをその運転パターンに応じた形状とすること
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項７】
上記標準負荷カーブは、資源の消費量と関連性を有する外的要因に対応付けられており、
推定対象の需要家又は推定対象となる期間の外的要因を特定し、特定した外的要因が対応付けられている標準負荷カーブを用いて、資源消費量を推定すること
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項８】
上記標準負荷カーブは、予測式にモデル化されていること
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項９】
特定したグループの標準負荷カーブを、需要家の消費規模又は積算測定値に応じて補正し、補正した値を需要家が消費した単位時間毎の資源量とすること
を特徴とする請求項１記載の資源消費量の推定方法。
【請求項１０】
供給者から需要家へネットワークを介して供給される資源の消費量を推定する資源消費量の推定装置において、
複数の需要家を資源の消費のパターンが類似するグループに区分し、資源の消費パターンを表した標準負荷カーブをグループ毎に格納した格納部と、
推定対象の需要家が属するグループを特定し、特定したグループの標準負荷カーブを上記格納部から取得し、取得した標準負荷カーブに基づき推定対象の需要家に対する任意の期間内の単位時間毎の資源消費量を推定する推定処理手段と
を備える資源消費量の推定装置。
【請求項１１】
供給者から需要家へネットワークを介して供給される資源の消費量を推定するコンピュータプログラムである資源消費量の推定プログラムにおいて、
複数の需要家を資源の消費のパターンが類似するグループに区分し、資源の消費パターンを表した標準負荷カーブをグループ毎に格納し、
推定対象の需要家が属するグループを特定し、特定したグループの標準負荷カーブを上記格納部から取得し、取得した標準負荷カーブに基づき推定対象の需要家に対する任意の期間内の単位時間毎の資源消費量を推定すること
を特徴とする資源消費量の推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば電気、ガス、水道等の資源の消費量を推定する資源消費量の推定方法及び装置に関するものである。
【背景技術】
【０００２】
電力会社は、家庭、工場、商店、その他ビル等の建物の需要家に対して、送電線のネットワークを介して電力の供給を行っている。電力料金の算出のため、通常、需要家には、電力計が設けられている。
【０００３】
ところで、電力予測や同時同量逸脱量の算出（例えば、特許文献１参照。）等のために、需要家毎に電力消費量の時間変化を計測したい場合がある。このような場合には、単位時間毎に電力量の測定及び記録が可能な計測器（インターバルメータ）や、リアルタイムで電力量の測定及び送信が可能な計測器（リアルタイムメータ）を、個々の需要家に設置すればよい。
【０００４】
しかしながら、インターバルメータ及びリアルタイムメータは通常の積算電力計と比べて非常に高価であり、全ての需要家に対してインターバルメータ又はリアルタイムメータを設置することは、コスト的に非常に困難である。
【０００５】
【特許文献１】特開２００３−１８００３２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
本発明は、インターバルメータやリアルタイムメータ等の高価な測定器を各需要家に設置することなく、需要家毎の例えば電気、ガス、水道等の資源の消費量を、正確且つ簡易に推定することができる資源消費量の推定方法及び装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明に係る資源消費量の推定方法は、供給者から需要家へネットワークを介して供給される資源の消費量を推定する資源消費量の推定方法であって、複数の需要家を資源量の消費のパターンが類似するグループに区分しておき、そのグループの資源の消費パターンを表した標準負荷カーブを、区分したグループ毎に生成しておき、推定対象の需要家が属するグループを特定し、特定したグループの標準負荷カーブを用いて、推定対象の需要家に対する任意の期間内の単位時間毎の資源消費量を推定することを特徴とする。
【０００８】
本発明に係る資源消費量の推定装置は、供給者から需要家へネットワークを介して供給される資源の消費量を推定する資源消費量の推定装置であって、複数の需要家を資源の消費のパターンが類似するグループに区分し、資源の消費パターンを表した標準負荷カーブをグループ毎に格納した格納部と、推定対象の需要家が属するグループを特定し、特定したグループの標準負荷カーブを上記格納部から取得し、取得した標準負荷カーブに基づき推定対象の需要家に対する任意の期間内の単位時間毎の資源消費量を推定する推定処理手段とを備えることを特徴とする。
【０００９】
本発明において、「資源」とは、例えば、電気、ガス、水道等の家庭、工場、商店等の需要家が消費するものであり、且つ、供給源（発電所又は変電所、ガスタンク、貯水設備）から需要家までの間を物理的に接続したネットワーク（送電線、ガス管、水道管）を介して供給されるものである。
【発明の効果】
【００１０】
以上の本発明に係る資源消費量の推定方法及び装置では、需要家のグループ毎にそのグループの資源消費量を表した標準負荷カーブを生成しておき、推定対象の需要家が属するグループを特定し、特定したグループの標準負荷カーブを用いてその需要家の消費資源量を推定する。
【００１１】
このことにより本発明に係る資源消費量の推定方法及び装置では、インターバルメータやリアルタイムメータ等の高価な測定器を各需要家に設置することなく、需要家毎の例えば電気、ガス、水道等の資源の消費量を、正確且つ簡易に推定することができる。
【発明を実施するための最良の形態】
【００１２】
以下、本発明が適用された負荷カーブ推定システムについて説明をする。
【００１３】
１  負荷カーブ推定システムの全体説明
（全体構成）
本発明が適用された負荷カーブ推定システム１０は、図１に示すような電力系統下に存在する需要家２の単位時間毎の消費電力量（負荷カーブ）を、推定するシステムである。
【００１４】
電力系統は、図１に示すように、少なくとも、電力の供給元となる発電所又は変電所（以下、発電所等という。）１と、電力を消費する１以上の需要家２（需要家２-１，２-２，２-３，…）と、発電所等１から各需要家２へ電力を送電する送電ネットワーク３とから構成されている。需要家２とは、家庭、工場、商店、その他ビル等の建物、街路灯等の電力設備、等々の電力を消費する機関の１単位である。
【００１５】
需要家の負荷カーブ（時刻毎の消費電力量）は、３０分毎や１時間毎といった単位時間の消費電力量を、１日や１週間といった一定期間に亘り表したものである。そのため、例えば、負荷カーブは、図２に示すように、横軸に時刻（単位時間刻み）、縦軸に電力消費量（ｋＷｈ）で表したグラフにより表現することができる。
【００１６】
負荷カーブ推定システム１０では、図１に示すような電力系統に存在する任意の一つの需要家２の負荷カーブを、直接測定するのではなく、推定処理を利用して算出するものである。算出する負荷カーブは、対象となる需要家の過去の実績値であってもよいし、将来の予測値であってもよい。
【００１７】
（処理手順）
図３に、負荷カーブ推定システムによる負荷カーブの算出手順を示す。
【００１８】
負荷カーブ推定システムでは、図３に示すように、負荷パターン区分設計処理（Ｓ１）、需要サンプリング調査設計処理（Ｓ２）、需要サンプリング結果集計処理（Ｓ３）、標準負荷カーブ作成処理（Ｓ４）、個別負荷カーブ推定処理（Ｓ５）を、順番に行う。なお、ステップＳ１からステップＳ４までの処理は、負荷パターン区分（詳細後述）毎に行われ、ステップＳ５の処理は、推定対象となる需要家毎に行われる。
【００１９】
--負荷パターン区分設計処理（Ｓ１）--
負荷パターン区分とは、負荷カーブの推定対象となる需要家をグループ分けするための区分である。例えば、地域（需要家の所在地）や負荷率（年平均負荷に対するピーク月の負荷の割合など）、業種等の、その区分に属する需要家同士で負荷カーブの形状が同一又は類似するグループである。負荷カーブ推定システムでは、負荷パターン区分毎に共通の負荷カーブ（標準負荷カーブと呼ぶ。）を作成し、この標準負荷カーブを用いて個々の需要家の負荷カーブを作成する。負荷パターン区分設計処理（Ｓ１）では、この負荷パターン区分を設定する。
【００２０】
--需要サンプリング調査設計処理（Ｓ２）--
上記標準負荷カーブは、例えば予めサンプルとして抽出した需要家（以下、サンプル需要家という。）の負荷カーブに基づき作成する。この場合、需要サンプリング調査設計処理（Ｓ２）では、サンプル需要家を抽出するための調査、必要となるサンプル需要家の数の設定、層別化処理といったサンプル需要家を抽出するための具体的な支援処理を行う。この需要サンプリング調査設計処理（Ｓ２）は、先のステップＳ１で設定された負荷パターン区分毎に実行する。
【００２１】
--需要サンプリング結果集計処理（Ｓ３）--
需要サンプリング結果集計処理（Ｓ３）では、需要サンプリング調査設計処理（Ｓ２）で得たサンプル需要家の各時刻の電力消費量を集計し、サーバ等に格納する。集計方法には、リアルタイム集計及びオンライン集計（後日一括して集計する集計方法。）がある。
【００２２】
--標準負荷カーブ作成処理（Ｓ４）--
標準負荷カーブ作成処理（Ｓ４）では、標準負荷カーブを作成する。標準負荷カーブと
は、負荷カーブを求めたい各対象日の標準的な負荷パターンのことである。標準負荷カーブは、負荷パターン区分設計処理Ｓ１で定めた負荷パターン区分毎に作成する。また、標準負荷カーブを作成するために、本処理（Ｓ４）では、標準負荷カーブモデルが用意され、これに対象日の気象条件、曜日、季節情報等の需要影響要因等を入力し、同負荷カーブを算出する。なお、標準負荷カーブは、個々の需要家の規模（電力消費量の大小）は考慮せず、正規化した（例えば、契約電力の大きさを１とするなど）値である。このため、次の個別負荷カーブ推定処理（Ｓ５）では、個々の需要家の規模に対する補正処理を実施する。
【００２３】
--個別負荷カーブ推定処理（Ｓ５）--
個別負荷カーブ推定処理（Ｓ５）では、先に作成した標準負荷カーブを各需要家の需要規模に応じて補正して、需要家毎の負荷カーブを作成する。需要規模の補正には対象需要家の電力会社との契約電力、月単位の検針電力量などを利用した増加減処理が用いられる。
【００２４】
２  負荷カーブ推定システムの装置構成
つぎに、本発明が適用された負荷カーブ推定システムの具体的な装置構成例について説明をする。
【００２５】
図４に、本発明が適用された負荷カーブ推定システム１０のシステム構成図を示す。
【００２６】
負荷カーブ推定システム１０は、図４に示すように、上記Ｓ１〜Ｓ５の計算処理を行う負荷カーブ作成用主計算機１１と、各種のデータや情報を格納するサーバ１２と、インターネットや電話回線等の外部ネットワーク１３を介して需要家に設置されている電力計量メータ１４の計測値を受信するデータ受信装置１５と、需要家に設置されている電力計量メータ１４の計測値が記録された記録媒体（フレキシブルディスクや光ディスク）から情報を読み出す光学・磁気記録媒体読取装置１６と、操作端末となる端末装置１７とを備えている。負荷カーブ作成用主計算機１１、サーバ１２、データ受信装置１５、光学・磁気記録媒体読取装置１６及び端末装置１７は、イントラネットワーク１８を介して互いにデータの送受信が行われる。
【００２７】
負荷カーブ作成用主計算機１１は、内部に、算出した標準負荷カーブを一時的に格納する標準負荷カーブ一時格納部２０と、負荷パターン区分を記憶する負荷パターン区分格納部２１と、サンプル設計結果を記憶するサンプル設計結果格納部２２と、標準負荷カーブモデルを記憶する標準負荷カーブモデル格納部２３とを備えている。また、サーバ１２は、内部に、需要調査データを記憶する需要調査データ格納部２４と、気象情報等を記憶する気象情報等格納部２５と、顧客情報を記憶する顧客情報格納部２６と、系統需要を記憶する系統需要格納部２７と、検針記録を記憶する検針記録格納部２８と、負荷カーブの推定結果を記憶する負荷カーブ推定結果格納部２９とを備えている。
【００２８】
つぎに、負荷カーブ推定システム１０の負荷カーブ作成用主計算機１１で行われる処理の手順について、図３に示した負荷カーブの算出手順と対応させて説明をする。
【００２９】
なお、図３に示した負荷カーブの算出手順を、ステップＳ１からステップＳ４までの負荷カーブ推定のための準備である負荷カーブ推定準備処理（図５参照）と、個別の対象需要家の負荷カーブを推定する個別負荷カーブ推定処理（図６参照）との二つの処理に分けて説明をする。
【００３０】
（負荷カーブ推定準備処理）
図５に示す負荷カーブ推定準備処理（Ｓ１〜Ｓ４）は、個々の需要家の負荷カーブ推定の基となる標準負荷カーブを作成するために必要となる標準負荷カーブモデルを作成するための一連の処理である。
【００３１】
まず、負荷パターン区分設計処理Ｓ１では、負荷カーブ推定対象である顧客（需要家）のグループ分けを行うための負荷パターン区分を設定する。負荷パターン区分は、利用者が直接指定する他、この負荷パターン区分設計処理Ｓ１により実行することが可能となっている。この処理計算では、サーバ１２内の各格納部２４〜２８に格納されたデータや情報が利用され、結果は、負荷カーブ作成用主計算機１１内の負荷パターン格納部２１に保存される。なお、本処理内容の詳細は後述する。
【００３２】
続いて、需要サンプリング調査設計処理Ｓ２では、上記負荷パターン区分設計処理Ｓ１で決定した負荷パターン区分毎に負荷カーブのサンプリング調査に対する設計（サンプリングする需要家の数、層別化などのサンプリング調査法の詳細設定）の支援計算を実施する。これらの支援計算では、サーバ１２内の各格納部２４〜２８に格納されたデータや情報が利用され、結果は、負荷カーブ作成用主計算機１１内のサンプル設計結果格納部２２に格納される。なお、本処理の詳細は後述する。
【００３３】
続いて、需要サンプリング結果集計処理Ｓ３では、需要サンプリング調査設計処理Ｓ２で得たサンプル需要家の各時刻の電力消費量（サンプル負荷カーブ）を集計する。サンプル負荷カーブは、サンプル需要家毎に設置した電力量計１４で計測される。計測された各サンプル負荷カーブは、外部ネットワーク１３やフレキシブルディスク等の可搬型記録媒体からデータ受信装置１５及び光学・磁気媒体読取装置１６を利用して集計され、サーバ１２の需要調査データ格納部２４に格納される。
【００３４】
続いて、標準負荷カーブモデル作成処理Ｓ４では、収集したサンプル負荷カーブをもとに負荷パターン区分設計処理Ｓ１で設定した各負荷パターン区分毎に標準負荷カーブモデルを設定する。標準負荷カーブモデルは、気象や曜日、各日の系統需要などの諸要因をパラメータとした負荷カーブ推定式（またはロジック）である。これの作成業務を支援するための諸計算が本標準負荷カーブモデル作成処理Ｓ４に含まれる。本処理Ｓ４により作成した標準負荷カーブモデルは、標準負荷カーブモデル格納部２３に格納される。また、標準負荷カーブモデルは、複数のアルゴリズムによるものが用意されており、これらの詳細は後述する。
【００３５】
（負荷カーブ推定準備処理）
図６に示す個別負荷カーブ推定処理（Ｓ５）では、個々の需要家の負荷カーブが作成される。
【００３６】
個別負荷カーブ推定処理（Ｓ５）では、まず始めに、負荷カーブ推定を行う対象顧客（対象需要家）並びに負荷カーブの推定対象期間（過去、現在又は将来の１週間や１ヶ月といった期間を指定）を設定する（Ｓ１５）。次に、各対象需要家に負荷パターン区分を割り当てる（Ｓ１６）。次に、割り当てられた負荷パターン区分に対して標準負荷カーブモデルに基づいて設定した推定対象期間に対する標準負荷カーブが算出される（Ｓ１７）。最後に算出した標準負荷カーブに対して各需要家の需要規模に応じた補正処理が行われ、負荷カーブが算出される（Ｓ１８）。
【００３７】
３  負荷パターン区分設計処理（Ｓ１）
つぎに、負荷パターン区分設計処理Ｓ１の詳細な処理内容について説明をする。
【００３８】
負荷カーブ推定システム１０では、個々の需要家の負荷カーブの形状は、業種や契約種別などのグループ内で共通しており同一とみなすことが可能であることを利用した方法である。すなわち、負荷カーブ推定システム１０では、需要家をいくつかのグループに分類し、グループ単位に共通の負荷カーブを作成（標準負荷カーブ）し、標準負荷カーブを基に個別の需要家の負荷カーブを算出する。負荷パターン区分設計処理Ｓ１では、この需要家グループを設定するためのグルーピング指標（負荷パターン区分）を設計するための計算処理を実施する。
【００３９】
（全体の処理内容）
図７に、負荷パターン区分設計処理Ｓ１の全体の流れを示す。負荷パターン区分設計処理Ｓ１は、負荷カーブ作成用主計算機１１により実行がされる。
【００４０】
まず、負荷パターン区分設計処理Ｓ１では、グループ分けを行う需要家の範囲を指定する（Ｓ４０）。負荷パターン区分設計処理Ｓ１では、指定した需要家範囲をさらにグループ分けするための区分を設計する。
【００４１】
図８に、需要家範囲の指定の一例を示す。需要家範囲の指定は、図８に示すように、所在地、電気料金の契約種別、業種、負荷率などの需要家の属性を表す各種のインデックスに対して、具体的な名称（関東、高圧）やその数値範囲等を指定することにより対象需要の範囲を設定する。指定は、利用者が負荷カーブ作成用主計算機１１の入力デバイスを通じて行う。
【００４２】
続いて、区分設計に用いる指標（区分設計指標）を選択する（Ｓ４１）。区分設計指標とは、需要家をグループ化するためのインデックスであり、顧客情報格納部２６に各需要家毎にこれらのインデックスとその属性情報を格納する。これらのインデックスの中から利用者が区分設計の尺度として用いる指標を一つもしくは複数指定する。
【００４３】
図９に、区分設計指標の一例を示す。区分設計指標には、所在地、契約種別、契約電力、業種などのカテゴリー変数のほか、負荷率などの数値変数も含まれる。
【００４４】
続いて、区分設計アルゴリズムを利用者が指定する（Ｓ４２）。負荷パターン区分設計処理Ｓ１では、凝集アルゴリズムと分割アルゴリズムとの２種類のアルゴリズムを用意している。利用者により指定がされると、いずれか一方のアルゴリズムが選択（Ｓ４３，Ｓ４４）される。
【００４５】
続いて、選択されたアルゴリズムにより、負荷パターン区分の設計が実行される（Ｓ４５，Ｓ４６）。なお、各アルゴリズムでの負荷パターン区分の設計処理手順については、後述する。
【００４６】
続いて、負荷パターン区分の設計が完了すると、設計された負荷パターン区分が指標に対応付けされて一時格納される（Ｓ４７）。
【００４７】
続いて、以上の処理により作成された各負荷パターン区分について、選択された区分指標でさらに細分化を行うか否かを利用者が判断する（Ｓ４８）。
【００４８】
続いて、ひきつづき負荷パターン区分の細分化を行う場合には、分割対象母集団の指定をし（Ｓ４９）、ステップＳ４１に戻り利用者が現在の負荷パターン区分の中から細分化を行う区分設計指標を一つ選択し（Ｓ４１）、以上のステップＳ４８までの処理を行う。なお、母集団や指標の指定は、自動で行うようにしてもよい。
【００４９】
最後に、負荷パターン区分の細分化を終了する場合には、上記全処理を通じて最終的に得られる負荷パターン区分（例えば、図１０参照。）を、負荷パターン区分格納部２１に格納（Ｓ５０）し、処理を終了する。
【００５０】
（負荷パターン区分設計のアルゴリズムの説明）
以下、負荷パターン区分の設計アルゴリズムについて説明する。まず、先の処理で区分設計指標とし指標を一つだけ選んだ場合（例えば、図８の例で所在地のみを選択した場合）のアルゴリズムについて説明する。
【００５１】
凝集アルゴリズムの処理（Ｓ４５）の流れを、図１１のフローチャートに示す。
【００５２】
凝集アルゴリズムでは、まず、選択した区分設計指標のカテゴリー変数化を行う（Ｓ４５−１）。凝集アルゴリズムは、ステップＳ４１で選択する区分設計指標がカテゴリー変数（業種や契約種別など）であることを前提としている。例えば、図９の負荷率などの数値変数が指標として選ばれた場合には、図１２に示すようにカテゴリー変数化を実施する。
【００５３】
続いて、凝集アルゴリズムでは、サーバ１２の需要調査データ格納部２４に予め集約されているサンプル需要家の負荷カーブと同サンプル需要家の顧客情報を取得する（Ｓ４５−２）。顧客情報には、図９に例示した区分設計指標に対する各需要家の属性情報が含まれている。
【００５４】
続いて、凝集アルゴリズムでは、初期クラスターを生成する（Ｓ４５−３）。ステップＳ４５−２で取得したサンプル需要家について、属性情報をもとにステップＳ４１で選択した区分設計指標を用いてクラス分けする。この場合、区分設計指標が示す各カテゴリーがそれぞれがクラスとなり、サンプル需要家はいずれかのクラスに割り当てられる。あるクラスに割り当てたサンプル需要家群のことをクラスターといい、ステップＳ４５−３で最初に生成されたクラスターのことを初期クラスターという。
【００５５】
続いて、凝集アルゴリズムでは、再クラスター化を行う（Ｓ４５−４）。
【００５６】
再クラスター化とは、任意の２つのクラスターを融合し、新たな一つのクラスターに融合する処理である。再クラスター化は、負荷カーブの形状が類似したクラスター同士が融合されるように行う。
具体的には、下記の（I）から（V）の処理を繰り返すことにより、再クラスター化が行われる。
【００５７】
（I）  クラスター（毎に、クラスターに属するサンプル需要家の負荷カーブの類似性を下記式（１）の評価値Ｅｇに基づき判断する。下記式（１）のＥｇは、クラスターに属するサンプル需要家の負荷カーブの偏差二乗和を尺度とした類似性を表す値である。
【００５８】
【数１】
2006011715.tif 000003
【００５９】
（II）  式（１）で得られるクラスターの類似性の評価値Ｅｇを全クラスターについて算出し、算出した評価値Ｅｇに対して、下式（２）に示すような、クラスター毎の需要家母集団数、サンプル需要家数を考慮した重みを加えた加重計算を行い、クラスター全体の類似性を表す尺度とする。
【００６０】
【数２】
2006011715.tif 000004
【００６１】
（III）  全クラスターのうち、任意の２つのクラスターのみを再クラスター化（融合）し、再度下記式（２）で得られるクラスター全体の類似性を表す尺度を計算する。同計算で得た尺度の再クラスター前と後の差（類似性変化値）を求める。この任意の２つのクラスターの組み合わせを再度換え次々に類似性の変化値を求め、変化値がもっとも小さいクラスターの組み合わせを抽出する。
【００６２】
（IV）  （III）の処理で抽出した変化値が最も小さいクラスターの組み合わせについて、クラスターとクラスターとの融合をする。この融合処理によりクラスターの数が一つ減る。
【００６３】
（V）  以上の（III）及び（IV）の処理を繰り返し実施し、クラスターが一つとなるまで処理を行う。
【００６４】
以上の（I）から（V）の処理が、再クラスター化の処理である（Ｓ４５−４）。
【００６５】
続いて、凝集アルゴリズムでは、区分の決定を行う（Ｓ４５−５）。
【００６６】
ここでは、各再クラスター化の各段階のクラスター状況並びに上記式（２）のクラスター全体の評価値との関係を利用者に提示する。例えば、図１３のようなクラスター数に対する融合過程を示す図、及び、図１４に示すような、クラスター数と評価値との関係を示す図を利用者に提示する。この提示を基に、利用者は、負荷パターン区分として利用するクラスターの細かさ（クラスター数）を決定し、計算機１１に入力する。
【００６７】
凝集アルゴリズムでは、以上のように利用者に決定された最終的にクラスターが、負荷パターン区分となる。
【００６８】
次に、分割アルゴリズムの処理（Ｓ４６）の流れを、図１５のフローチャートを参照して説明する。
【００６９】
分割アルゴリズムでは、まず、サーバ１２の需要調査データ格納部２４に予め集約されているサンプル需要家の負荷カーブと同サンプル需要家の顧客情報を取得する（Ｓ４６-１）。顧客情報には、図９に例示した区分設計指標に対する各需要家の属性情報が含まれている。なお、本方法で用いる区分設計指標は負荷率や契約電力などの数値データもしくは順序データのみを対象とする。
【００７０】
続いて、分割アルゴリズムでは、取得したサンプル需要家の負荷カーブをもとに、取得したサンプル全体の類似性を下記式（３）の評価値Ｅｏに基づき評価する（Ｓ４６-２）。
【００７１】
【数３】
2006011715.tif 000005
【００７２】
続いて、分割アルゴリズムでは、サンプル需要家を二つのグループに分割する（Ｓ４６-３）。具体的には、下記の（I）から（IV）の処理を行うことにより分割が行われる。
【００７３】
（I）  選択した区分設計指標について任意の値を選び、選択した値を分割点としてサンプル需要家を２つのグループに分割する。
【００７４】
（II）  分割した各グループに対して、上記（３）式により需要の類似性を評価する。その後、下記式（４）により、分割後のグループ全体の類似性Ｅｏ´を求める。
【００７５】
【数４】
2006011715.tif 000006
【００７６】
（III）  ステップＳ４６−２で算出した分割前のグループ全体の類似性Ｅｏと、上記処理により算出した分割後のグループ全体の類似性Ｅｏ´とを、下記式（５）のように比較し、グループ分割による類似性の向上度Δｏを求める。
【００７７】
【数５】
2006011715.tif 000007
【００７８】
（IV）  分割点を代えて、（I）から（III）までの処理を行い、上記式（５）で得られる向上度Δｏの最も大きな分割点を求め、この点を分割時の最適な分割点として、２つのグループに２分割する。
【００７９】
以上の（I）から（IV）の処理が、グループの分割化の処理である（Ｓ４６−３）。
【００８０】
続いて、分割アルゴリズムでは、分割した各グループについて上記ステップＳ４６−２からステップＳ４６−３の処理を行って再分割化をし、分割数を増加させていく（ステップＳ４６−４）。
【００８１】
続いて、分割アルゴリズムでは、区分の決定を行う（Ｓ４６−５）。
【００８２】
ここでは、各再分割化の各段階のグループ状況並びに上記式（４）のグループ全体の評価値との関係を利用者に提示する。例えば、図１６のようなグループ数に対する分割過程を示す図、及び、図１７に示すようなグループ数と評価値との関係を示す図を利用者に提示する。この提示を基に、利用者は、負荷パターン区分として利用する需要グループの細かさ（分割数）を決定し、計算機１１に入力する。
【００８３】
分割アルゴリズムでは、以上のように利用者に最終的に決定されたグループが、負荷パターン区分となる。
【００８４】
つぎに、分割アルゴリズムを選択した場合の処理で、区分設計指標を複数選んだ場合（例えば、図８の例で所在地と契約種別の２つを選択した場合）を例に負荷パターン区分設計処理での動作について説明する。なお、処理手順は、図１５のフローチャートと同様となる。
【００８５】
分割アルゴリズムでは、まず、サーバ１２の需要調査データ格納部２４に予め集約されているサンプル需要家の負荷カーブと同サンプル需要家の顧客情報を取得する（Ｓ４６-１）。顧客情報には、図９に例示した区分設計指標に対する各需要家の属性情報が含まれている。
【００８６】
続いて、分割アルゴリズムでは、取得したサンプル需要家の負荷カーブをもとに、取得したサンプル全体の類似性を上記式（３）の評価値Ｅｏに基づき評価する（Ｓ４６-２）。
【００８７】
続いて、分割アルゴリズムでは、サンプル需要家を二つのグループに分割する（Ｓ４６-３）。具体的には、下記の（I）から（IV）の処理を行うことにより分割が行われる。
【００８８】
（I）  選択した区分設計指標の中から任意の指標を選び、その任意の指標について任意の値を選び、選択した値を分割点としてサンプル需要家を２つのグループに分割する。
【００８９】
（II）  上記式（４）及び式（５）により最適な分割点を求める。すなわち、分割前の
類似性Ｅｏと分割後の類似性Ｅｏ´の差が最も大きくなる点を分割点とする。
【００９０】
（III）  （I）で選択した区分設計指標とは別の指標を選び、（I）及び（II）の処理を行って最適な分割点を求め、分割を行う。
【００９１】
（IV）  全ての区分設計指標について最適な分割点を求め、式（５）の比較を行い、グループ分割による類似性の向上度Δｏを求め、向上度Δｏの最も高い分割点でグループを２つに分割する。
【００９２】
以上の（I）から（IV）の処理が、グループの分割化の処理である（Ｓ４６−３）。
【００９３】
続いて、分割アルゴリズムでは、分割した各グループについてさらに上記ステップＳ４６−２からステップＳ４６−３の処理を行って再分割化をし、分割数を増加させていく（ステップＳ４６−４）。
【００９４】
続いて、分割アルゴリズムでは、区分の決定を行う（Ｓ４６−５）。
【００９５】
ここでは、各再分割化の各段階のグループ状況並びに上記式（４）のグループ全体の評価値との関係を利用者に提示する。例えば、図１８のようなグループ数に対する分割過程を示す図、及び、図１７に示したグループ数と評価値との関係を示す図を利用者に提示する。この提示を基に、利用者は、負荷パターン区分として利用する需要グループの細かさ（分割数）を決定し、計算機１１に入力する。
【００９６】
分割アルゴリズムでは、以上のように利用者に最終的に決定されたグループが、負荷パターン区分となる。
【００９７】
４  需要サンプリング調査設計処理（Ｓ２）
つぎに、需要サンプリング調査設計処理Ｓ２の詳細な処理内容について説明をする。需要サンプリング調査設計処理Ｓ２は、負荷カーブ作成用主計算機１１により実行がされる。
【００９８】
標準負荷カーブモデルは、サンプルとなる需要家から別途収集したサンプル負荷カーブをもとに作成される。サンプル負荷カーブの収集に際しては、より経済的に（すわなち、できるだけ少数のサンプル需要家を用いて）、また、偏りが少なく（ある特徴的な負荷カーブ形状をもつ需要家のみにサンプルが固まってしまうことなく、負荷カーブ区分として定めた需要家クラスの平均的な負荷カーブの特徴が捉えられるように）、サンプル需要家を選定することが必要である。
【００９９】
需要サンプリング調査設計処理Ｓ２では、負荷カーブのサンプリングのためのサンプル設計（サンプルとして取得する需要家の数やどういう種別の需要家を選択するのか、といった標本調査）を支援する。
【０１００】
図１９に、需要サンプリング調査設計処理Ｓ２の全体の流れを示す。
【０１０１】
需要サンプリング調査設計処理Ｓ２では、まず始めに、サンプル抽出枠を設定する（Ｓ５１）。サンプル抽出枠とは、負荷カーブのサンプリングを目的とする需要家クラスのことである。通常は負荷パターン区分設計処理（Ｓ１）で求めた各負荷パターン区分をサンプル抽出枠として設定する。図２０に、サンプル抽出枠の設定例を示す。図２０に示すように、負荷パターン区分格納部２１に蓄積されている各負荷パターン区分の一つまたは複数の区分を選び、これをサンプル抽出枠として指定する。
【０１０２】
続いて、需要サンプリング調査設計処理Ｓ２では、目標精度を指定する（Ｓ５２）。目標精度は、サンプリングにより得られる負荷カーブのサンプル枠内の全需要家の負荷カーブに対する差を評価するための指標である。本処理では、サンプル枠内の全需要家の負荷カーブの平均値に対するサンプル負荷カーブの差（サンプル需要家を複数選ぶ場合には全サンプル負荷カーブの平均値をサンプル負荷カーブとする）の大きさを確率的に評価する。具体的には、下式のような言語ルールの形で計算機１１に対して目標精度としてα，β，γを指定する。
【０１０３】
精度の指定  ＝
α％の確率水準で発生する誤差を、β×１００％（又はγｋＷｈ）以下とする。
【０１０４】
続いて、需要サンプリング調査設計処理Ｓ２では、サンプルの抽出法を指定する（Ｓ５３）。
【０１０５】
サンプルの抽出法としては、例えば、無作為抽出、系統抽出、集落抽出といった方法がある。無作為抽出は、サンプル抽出枠内の全需要家の中からランダムにサンプル需要家を選ぶ方法である。系統抽出は、サンプル抽出枠内の全需要家の顧客番号をもとに一定ルールで確定的に抽出する方法である。集落抽出は、サンプル抽出枠内の全需要家を電力系統の変圧器などの単位で一括抽出する方法である。
【０１０６】
また、ここでは、併せて層別化の有無の指定をする。なお、層別化については後に説明する。また、上記抽出法により抽出されるサンプル需要家については、図２１に示すように、顧客情報格納部２６内の顧客情報に、フラグが付けられている。
【０１０７】
続いて、需要サンプリング調査設計処理Ｓ２では、上記のサンプル抽出法において層別化を併用する場合には、補助変数の選択を行う（Ｓ５４）。補助変数は、顧客情報格納部に格納されている各需要家の顧客情報の中から一つ、もしくは複数選択される。補助変数の候補は、例えば各需要家の“月間電力量”、“負荷率（月間電力量から算出）”、“業種”等である。
【０１０８】
続いて、需要サンプリング調査設計処理Ｓ２では、サンプル設計を実施する（Ｓ５５）。
【０１０９】
サンプル設計（Ｓ５５）では、サンプル抽出する需要家の数を求めるための計算処理を実施する。このサンプル需要家の抽出件数は、Ｓ５２で指定した目標精度を満たすように自動的に計算が行われる。
【０１１０】
最後に、これらのサンプル設計結果をサンプル設計結果格納部２２に保存し、本処理を終了する。
【０１１１】
次に、ステップＳ５５のサンプル抽出する需要家の件数を求める演算処理についてさらに説明をする。
【０１１２】
まず、層別化を併用しないでサンプル設計を行う場合を例に説明する。
【０１１３】
無作為抽出を指定した場合、目標精度として確率をα、誤差をβ、ｔαを自由度Ｎ−１のｔ分布の確率点αの値、サンプル抽出枠内の全需要家の数をＮとすると、下記式（６）によりサンプル抽出する需要家数ｎが算出される。
【０１１４】
【数６】
2006011715.tif 000008
【０１１５】
ただし、ばらつき（標準偏差）及び平均値を得ることは不可能であるため、下記式（７−１）、（７−２）のように予め事前サンプリングなどにより得た負荷カーブからの値、もしくは利用者の主観に基づいて設定する。
【０１１６】
【数７】
2006011715.tif 000009
【０１１７】
なお、サンプル需要はサンプル需要家の電力消費規模に大きな開きがある場合がある。その場合には必要に応じて次のように各サンプル需要家の負荷カーブの最大値等をもとに正規化する。
【０１１８】
【数８】
2006011715.tif 000010
【０１１９】
次に、層別化をあわせて行う場合について説明する。
【０１２０】
層別化は、サンプル抽出枠として選んだ負荷パターン区分をさらに細分化することで、同じサンプル数でより高いサンプル抽出精度を達成するための方法である。具体的な方法としてサンプル抽出法として無作為抽出を選んだ場合を例に説明する。この場合には、まず層別化を実施した場合のサンプル設計の精度を求めるために、下式（９）に基づき分散を計算する。
【０１２１】
【数９】
2006011715.tif 000011
【０１２２】
次に、下記式（１０）によりサンプル誤差を算出する。
【０１２３】
【数１０】
2006011715.tif 000012
【０１２４】
なお、部分抽出枠内の補助変数のばらつきや平均値を得ることは不可能であるため、式（７−１）、式（７−２）と同様に、予め事前サンプリングなどにより得た補助変数からの値、もしくは利用者の主観に基づいて設定する。
【０１２５】
層別化処理では、上記式（１０）の値と、目標誤差Ｅとが下記式（１１）の条件を満たすようなサンプル抽出枠の分割を実施する。
【０１２６】
【数１１】
2006011715.tif 000013
【０１２７】
また、抽出枠の再分割の方法としては、（Ａ）利用者が任意に分割する方法、（Ｂ）先の負荷パターン区分設計処理を利用する方法。（Ｃ）Dalenius-Hodges法を利用する方法。という３種類の方法を用いることができる。
【０１２８】
（Ｂ）の方法では、負荷パターン区分格納部２１に格納されている需要家を対象母集団として、当該サンプル抽出枠に設定する。その後、区分設計指標の選択を先に指定した補助変数を指定する。その後、同処理フローにそってサンプル抽出枠の分割を行う。
【０１２９】
なお、分割後の部分抽出枠のサンプル需要家の数については、下記式（１２）のネイマン配分等を利用することができる。
【０１３０】
【数１２】
2006011715.tif 000014
【０１３１】
５  標準負荷カーブ作成処理（Ｓ４）
つぎに、標準負荷カーブ作成処理Ｓ４の詳細な処理内容について説明をする。標準負荷カーブ作成処理Ｓ４は、負荷カーブ作成用主計算機１１により実行がされる。
【０１３２】
負荷カーブ推定システム１０では、個別の需要家の負荷カーブを標準負荷カーブモデルにより作成する。標準負荷カーブモデルは、負荷パターン区分毎に計算機１１内に用意された負荷カーブ推定式又は負荷カーブの推定ロジックのことである。
【０１３３】
標準負荷カーブモデルから出力される負荷カーブを、標準負荷カーブと呼ぶ。標準負荷カーブは、個々の需要家の規模（電力消費量の大小）の違いの補正が容易となるように正規化（例えば一年間の最大需要を１とするなど）がされた負荷カーブである。
【０１３４】
（標準負荷カーブモデルの種類）
電力の負荷カーブは、需要家の業種などの他、季節や曜日、時刻、さらには各日の気象条件などの様々な要因の影響を受ける。標準負荷カーブモデルとは、負荷パターン区分毎に、こうした各日の需要影響要因を加味したうえで、対象となる任意の日（または一定期間）の負荷カーブを作成するための負荷カーブ推定式又は負荷カーブの推定ロジックのことである。
【０１３５】
標準負荷カーブモデルの具体例としては、「共通形状モデル」、「即時サンプルモデル」、「運転パターンモデル」、「類似日モデル」、「需要モデリングモデル」といった種類を用意する。利用者は、これらのいずれか１つを選択若しくは複数を選択（この場合には出力される複数の標準負荷カーブモデルの平均化等を行い、一つのモデルに合成する。）して利用する。なお、これらのモデルは一例であり、これら以外のモデルを用いてもよい。また、後述する個別毎の需要家の負荷カーブ作成（Ｓ５）の処理に際しては、これらの中から任意に一つのモデルを選択する他、計算可能なモデルについて全て計算を行うことも可能である。
【０１３６】
以下、「共通形状モデル」、「即時サンプルモデル」、「運転パターンモデル」、「類似日モデル」、「需要モデリングモデル」について、具体的に説明を行う。
【０１３７】
--共通形状モデル--
共通形状モデルは、電力系統全体（発電所の発電出力パターンの合計が相当）やその一部（変電所内の変圧器の二次側電力量の毎時計測値等が相当）など、単位時間毎の電力量の実測パターンから標準負荷カーブを作成する標準負荷カーブモデルである。この標準負荷カーブモデルは、負荷パターン区分に地域単位など大まかな区分が適用された場合に利用することができる。また、系統内にインターバルメータなどの単位時間毎の電力量が測定できる需要家が存在する場合には、その需要分を差し引くなどの処理を利用者の選択によって実施するものとしている。
【０１３８】
図２２に、共通形状モデルによる標準負荷カーブの作成方法の具体例を示す。
【０１３９】
共通形状モデルでは、まず、発電所等１（又は変電所）での対象日（ｘ年ｙ月ｚ日）の電力系統の全体の需要（カーブ１）を測定するとともに、対象日（ｘ年ｙ月ｚ日）のインターバルメータが設定された需要家２−１の負荷の実測値（カーブ２）を直接計測する。
【０１４０】
続いて、インターバルメータが設定された需要家４−１負荷の実測値を、電力系統の全体の需要から差し引く。これにより対象となる需要家群（例えば、地区Ａの需要家２−２〜２−Ｎ）の対象日（ｘ年ｙ月ｚ日）の合計の負荷カーブを計算する。
【０１４１】
続いて、負荷パターンを正規化し、標準負荷カーブを求める。
【０１４２】
正規化方法として、例えば、下記式（１３）に示すような契約電力を用いた方法を用いた方法を適用することができる。なお、この方法に限らず、他の方法を用いて正規化してもよい。
【０１４３】
【数１３】
2006011715.tif 000015
【０１４４】
--即時サンプルモデル--
共通形状モデルでは、系統単位での負荷カーブの計測値を利用したが、即時サンプルモデルでは、負荷パターン区分毎に１または複数のサンプル需要家を設定しておき、設定されたサンプル需要家について対象日（ｘ年ｙ月ｚ日）の単位時間毎の需要電力量を計測し、その計測値に基づき標準負荷カーブを作成する。
【０１４５】
図２３に、即時サンプルモデルによる標準負荷カーブの作成処理の流れを示す。
【０１４６】
即時サンプルモデルでは、まず、サンプル需要家に設置したリアルタイムメータ又はインターバルメータにより、対象日（ｘ年ｙ月ｚ日）の単位時間毎の電力量を計測する（Ｓ１０１）。なお、サンプル需要家は、需要サンプリング調査設計処理（Ｓ２）により設定された需要家である。
【０１４７】
続いて、即時サンプルモデルでは、計測された各サンプル需要家の電力量を単位時間毎に合計し、合計の負荷カーブを作成する（Ｓ１０２）。なお、必要に応じてサンプル需要家の需要規模（契約電力）等を考慮して加重平均してもよい。
【０１４８】
続いて、即時サンプルモデルでは、得られたサンプル需要家の負荷カーブをサンプル需要家の契約電力等を基にして正規化し、これを対象日（ｘ年ｙ月ｚ日）の標準負荷カーブとする（Ｓ１０３）。
【０１４９】
即時サンプルモデルでは、以上のような処理を毎日行うことにより、それぞれの日の標準負荷カーブを生成する。
【０１５０】
--運転パターンモデル--
運転パターンモデルは、街灯等の動作パターンがある程度確定的な需要家について、予め負荷パターンをモデル化し、それを標準負荷カーブとしてそのまま利用するものである。
【０１５１】
運転パターンモデルとして利用できる需要家には、例えば、街路灯、信号機等がある。街路灯は、日没に点灯し、日の出に消灯する。そのため、街路灯の負荷パターンは日没後には、電力量がほぼ所定の値で一定となり、日の出後にはほぼ０で一定となる。従って、図２４（ａ）に示すように地域ごとの日の出、日の入り時刻を変数とした動作パターンとしてモデル化が可能となる。また、信号機は、一日を通じて動作する。従って、図２４（ｂ）に示すように、常にほぼ一定の需要電力量となる。すなわち、機器の定格容量や負荷率等を考慮した動作パターンとしてモデル化できる。
【０１５２】
--類似日モデル--
類似日モデルでは、図２５に示すように、即時サンプルモデルと同様に、負荷パターン区分毎に、各日（ｘ年ｙ月ｚ日）の標準負荷カーブを求める（Ｓ１０１〜Ｓ１０３）。
【０１５３】
さらに、類似日モデルでは、測定日の曜日、気温、天気等の外的要因をインデックスとして作成し（Ｓ１０４）、作成したインデックスと標準負荷カーブとを対応付けてデータベース化する（Ｓ１０５）。
【０１５４】
図２６に、インデックスと対応付けられた標準負荷カーブのデータベース例を示す。このデータベースは、インデックスとして、日付、曜日、気温や天気、さらには系統需要などさまざまなものを用いることが可能であり、これらの選択は利用者任意に選択可能とする。また、曜日や天気といったカテゴリ変数については、電力負荷への影響を考慮したうえで、例えば、晴＝０、曇＝１、雨＝２、雪＝３、月〜金＝０、土日祝日＝１といったように数値化する。
【０１５５】
類似日モデルでは、その後、具体的に個々の需要家の負荷カーブを求めたい日（対象日）の条件に近い過去の日の負荷カーブを、データベースから抽出して標準負荷カーブとして利用する。
【０１５６】
データベースからの標準負荷カーブの抽出は、図２７に示すように、データベース内の過去の各日のインデックス値（Ｓ１１１）と、負荷カーブを求めたい対象日の諸条件（Ｓ１１２）との類似性を計算し（Ｓ１１３）、類似性が高い日の負荷カーブを、標準負荷カーブとして抽出する（Ｓ１１４）。
【０１５７】
類似性は、例えば、下記式（１４）等により計算することができる。
【０１５８】
【数１４】
2006011715.tif 000016
【０１５９】
次に、標準負荷カーブの抽出条件の設定例（ケース１〜ケース３）を示す。
【０１６０】
設定ケース１では、負荷カーブ推定対象日と月、曜日が一致する過去の日に対する標準負荷カーブをデータベースから抽出し利用する。抽出結果が複数の場合には平均化等する。
【０１６１】
設定ケース２では、負荷カーブの算出の対象日の全系統の需要パターン（共通形状モデルでの説明参照）と類似した過去の日に対応する標準負荷カーブをデータベースから抽出して利用する。この場合、図２６に示す系統需要のインデックスの毎時の値を、上記式（１４）の各条件変数に指定することで抽出が可能となる。
【０１６２】
設定ケース３では、負荷カーブの算出の対象日の気象条件と類似した過去の日に対する標準負荷カーブをデータベースから検索して利用する。この場合、過去の日の気象値（気温、天気など）と、対象日の気象値とを、上記式（１４）の各条件変数に指定することで抽出が可能となる。なお、負荷カーブの算出の対象日が未来（例えば翌日など）の場合には、当日の気象値に予報気象の値を利用するなどすればよい。また、気象条件が天候などの場合には、図２８に示すように、予め天気の各条件を数値化した値を利用する。
【０１６３】
--需要モデリングモデル--
需要モデリングモデルは、負荷パターン区分毎に標準負荷カーブの各時刻の値と他の要因との関係を過去の標準負荷カーブ値と同要因との統計データに基づいて数式化した予測式を作成し、これに負荷カーブの算出の対象日の要因の値を代入し、対象日の標準負荷カーブを得る方法である。
【０１６４】
需要モデリングモデルに用いる数式化手法としては、「回帰モデル」、「時系列モデル」、「ニューラルネットワーク」等がある。ここでは、「回帰モデル」を例にとって説明をする。
【０１６５】
回帰モデルでは、まず、下記式（１５）のような、標準負荷カーブの各時刻の値を目的変数、気温などの影響を及ぼす要因を説明変数とした需要モデリングモデルを作成する。
【０１６６】
【数１５】
2006011715.tif 000017
【０１６７】
例えば、比較的単純なモデルとして、下記式（１６）に示すような気温と曜日フラグを説明変数とした回帰モデルを用いることもできる。
【０１６８】
【数１６】
2006011715.tif 000018
【０１６９】
なお、需要モデリングモデルによる推定方法（回帰係数値など）では、予め収集したサンプル需要家の負荷カーブを利用する。この場合、作成した需要モデリングモデルの出力が標準負荷カーブ（正規化された負荷カーブ）となるようにするためには、得られるサンプル需要家の負荷カーブを先に述べた即時サンプルモデルでの処理と同様に、日毎にサンプル需要家の負荷カーブを合計したうえで契約電力などをもとに正規化し、この正規化した負荷カーブを過去データとして、これと気象要因等の説明変数との関係を上記式（１５）によりモデル化すればよい。
【０１７０】
６  個別負荷カーブ推定処理（Ｓ５）
つぎに、個別負荷カーブ推定処理Ｓ５の詳細な処理内容について説明をする。個別負荷カーブ推定処理Ｓ５は、負荷カーブ作成用主計算機１１により実行がされる。
【０１７１】
個別負荷カーブ推定処理（Ｓ５）では、図６に示したように、負荷カーブ推定を行う対象顧客（対象需要家）並びに負荷カーブの推定対象期間を設定し（Ｓ１５）、次に、各対象需要家に負荷パターン区分を割り当て（Ｓ１６）、次に、負荷パターン区分に対して標準負荷カーブモデルに基づいて設定した推定対象期間に対する標準負荷カーブが算出し（Ｓ１７）、最後に、算出した標準負荷カーブに対して各需要家の需要規模に応じた補正処理を行う（Ｓ１８）。このことにより、各需要家に対する個別の負荷カーブを作成することができる。
【０１７２】
--標準負荷カーブの算出処理（Ｓ１７）--
まず、標準負荷カーブの算出処理（Ｓ１７）の処理内容について、標準負荷カーブモデル毎（「共通形状モデル」、「即時サンプルモデル」、「運転パターンモデル」、「類似日モデル」、「需要モデリングモデル」）に説明する。
【０１７３】
（共通形状モデル）
共通形状モデルでの個別負荷カーブ推定処理におけるステップＳ１７の流れを、図２９を参照して説明をする。
【０１７４】
共通形状モデルでは、まず、標準負荷カーブの計測期間、負荷カーブ作成対象となる需要家の負荷パターン区分が設定される（Ｓ２１１）。
【０１７５】
続いて、共通形状モデルとして参照する系統範囲が設定される（Ｓ２１２）。ここで、系統範囲として、一部の系統の負荷カーブを用いる場合には、続くステップＳ２１３に進み、全ての系統の負荷カーブを利用する場合に、続くステップＳ２１４に進む。
【０１７６】
続いて、各系統の負荷カーブの実績値が収集される（Ｓ２１５，Ｓ２１６）。実績値は、発電所又は変電所に設置されている電力計測メータ又は変電所に設定されている電力計測メータにより計測され、系統需要格納部２７に格納されている。
【０１７７】
続いて、指定する各系統内にリアルタイムメータやインターバルメータが設置されている需要家（メータ設置需要家）が存在する場合には、その負荷カーブの実測値が収集される。収集された系統の負荷カーブからメータ設置需要家の負荷カーブを差し引いて、メータ設置需要家以外の需要家の合計の負荷カーブが計算される（Ｓ２１７）。
【０１７８】
続いて、メータ設置需要家以外の需要家の契約電力や積算メータによる検針結果等をもとに、上記式（１３）を用いて合計の負荷カーブが正規化され、標準負荷カーブが作成される（Ｓ２１８）。
【０１７９】
最後に、得られた標準負荷カーブが、標準負荷カーブ一時格納部２０に格納される（Ｓ２１９）。
【０１８０】
（即時サンプルモデル）
即時サンプルモデルでの個別負荷カーブ推定処理におけるステップＳ１７の流れを、図３０を参照して説明をする。
【０１８１】
即時サンプルモデルでは、まず、標準負荷カーブの計測期間、負荷カーブ作成対象となる需要家の負荷パターン区分が設定される（Ｓ２２１）。
【０１８２】
続いて、負荷パターン区分に対するサンプル需要家の計測期間の負荷カーブの実測値を収集し、合計される（Ｓ２２２）。
【０１８３】
続いて、サンプル需要家の契約電力等をもとに、上記式（１３）等を用いて正規化処理が行われ、標準負荷カーブが作成される（Ｓ２２３）。
【０１８４】
最後に、得られた標準負荷カーブが、標準負荷カーブ一時格納部２０に格納される（Ｓ２２４）。
【０１８５】
（運転パターンモデル）
運転パターンモデルでの個別負荷カーブ推定処理におけるステップＳ１７の流れを、図３１を参照して説明をする。
【０１８６】
運転パターンモデルでは、まず、標準負荷カーブの計測期間及び負荷カーブ作成対象となる需要家の種類（電力需要の種類：電気機器など）が設定される（Ｓ２３１）。
【０１８７】
続いて、その需要家に対応した運転パターンモデルが呼び込まれる（Ｓ２３２）。
【０１８８】
続いて、呼び込んだ運転パターンモデルについて対象期間の運転パターンモデルを計算するために必要となる変数の値（日の出、日の入り時刻など）が入力される（Ｓ２３３）。なお、これらの入力値は利用者が入力インタフェースを介して直接入力する他、予め標準負荷カーブモデルとして、前もってカレンダー情報とともに計算機内にデータ化し、これを自動的に与えることも可能である。
【０１８９】
続いて、入力結果をもとに運転パターンモデルが計算され、これを標準負荷カーブとする（Ｓ２３４）。
【０１９０】
最後に、得られた標準負荷カーブが、標準負荷カーブ一時格納部２０に格納される（Ｓ２３５）。
【０１９１】
（類似日モデル）
類似日モデルでの個別負荷カーブ推定処理におけるステップＳ１７の流れを、図３２を参照して説明をする。
【０１９２】
類似日モデルでは、まず、標準負荷カーブの計測期間、負荷カーブ作成対象となる需要家の負荷パターン区分が設定される（Ｓ２４１）。
【０１９３】
続いて、負荷パターン区分に対する、上記式（１４）の類似距離計算に用いる変数及び重み係数が設定される（Ｓ２４２）。なお、変数及び重み係数は、利用者が直接システムの入力インタフェースを介して設定可能であるほか、予め標準負荷カーブモデルとして保存した値を利用することも可能である。
【０１９４】
続いて、対象となる負荷パターン区分について、過去のサンプル需要家から得た各日の標準負荷カーブが抽出対象リストに載せられる（Ｓ２４３）。
【０１９５】
続いて、上記抽出対象リスト内の過去各日と先に指定した計測対象期間内の各日との類似性（類似距離）が計算され、類似距離が小さい過去の日を抽出される（Ｓ２４４）。
【０１９６】
続いて、抽出した過去日について標準負荷カーブをデータベースより呼び込まれる（Ｓ２４５）。
【０１９７】
続いて、抽出した過去の日が複数の場合には標準負荷カーブを平均化し、一つの標準負荷カーブが求められる（Ｓ２４６）。
【０１９８】
最後に、得られた標準負荷カーブが、標準負荷カーブ一時格納部２０に格納される（Ｓ２４７）。
【０１９９】
（需要モデリングモデル）
需要モデリングモデルでの個別負荷カーブ推定処理におけるステップＳ１７の流れを、図３３を参照して説明をする。
【０２００】
需要モデリングモデルでは、まず、標準負荷カーブの計測期間、負荷カーブ作成対象となる需要家の負荷パターン区分が設定される（Ｓ２５１）。
【０２０１】
続いて、該当する負荷パターン区分に対応した需要モデリングモデルが呼び込まれる（Ｓ２５２）。
【０２０２】
続いて、需要モデリングモデルの説明変数に用いている各変数のステップＳ２５１で設定された計測期間内の各時刻毎の値が設定される（Ｓ２５３）。設定は、利用者が入力インタフェースを介して直接入力する他、関連する気象や需要等のデータベース２４，２５に格納されている当該計測期間の変数値により自動設定することも可能である。
【０２０３】
続いて、需要モデリングモデルに基づき計測期間の各時刻の需要値が計算され、標準負荷カーブが求められる（Ｓ２５４）。
【０２０４】
最後に、得られた標準負荷カーブが、標準負荷カーブ一時格納部２０に格納される（Ｓ２４７）。
【０２０５】
（オンライン処理等について）
以上のような標準負荷カーブの作成処理では、「オフライン処理」、「オンライン処理」及び「予測処理」の３つの処理を適用することが可能である。
【０２０６】
オフライン処理とは、過去に遡ってある一定期間若しくはある特定日の標準負荷カーブを作成し、これに基づいて、個々の需要家の負荷カーブを作成する処理のことである。
【０２０７】
オンライン処理とは、現在若しくはごく至近の過去（数時間以内など）に対応した標準負荷カーブを作成し、これに基づいて、個々の需要家の負荷カーブを作成する処理である。
【０２０８】
予測処理とは、将来のある一定期間若しくはある特定日の標準負荷カーブを作成し、これに基づいて、個々の需要家の負荷カーブを作成する処理である。
【０２０９】
「共通形状モデル」、「即時サンプルモデルモデル」、「運転パターンモデル」、「類似日モデル」、「需要モデリングモデル」の各モデルでは、それぞれ下記表１のように対応することが可能である。
【０２１０】
【表１】
2006011715.tif 000019
【０２１１】
--計算対象設定処理（Ｓ１５）--
つぎに、個別負荷カーブ推定処理Ｓ５の計算対象設定処理（図６のＳ１５）について説明をする。
【０２１２】
図３４に、計算対象設定処理（Ｓ１５）の処理の流れを示す。
【０２１３】
計算対象設定処理では、まず、予め登録されている若しくは利用者が直接インタフェイス等を介して指定する顧客名が入力される（Ｓ１５−１）。
【０２１４】
続いて、計算対象設定処理では、顧客名を基に、顧客情報格納部２６から、該当する負荷カーブ作成対象となる顧客の契約種別、契約電力負荷率、電気使用量、最大電力、需要家業種、地域などの負荷パターン区分の設定に必要なデータが読み込まれる（Ｓ１５-２）。
【０２１５】
続いて、計算対象設定処理では、負荷カーブ作成の対象需要家について割り当てられている負荷パターン区分及び負荷カーブ作成期間、単位時間などを設定し、これらの計算条件等を一時的に保存する（１５-３）。
【０２１６】
--負荷パターン区分の割当処理（Ｓ１６）--
つぎに、個別負荷カーブ推定処理Ｓ５の負荷パターン区分の割当処理（図６のＳ１６）について説明をする。
【０２１７】
図３５に、負荷パターン区分の割当処理（Ｓ１６）の処理の流れを示す。
【０２１８】
負荷パターン区分の割当処理では、まず、計算対象設定処理（Ｓ１５）で格納された計算条件を読み出す（Ｓ１６−１）。
【０２１９】
続いて、負荷パターン区分の割当処理では、読み出した計算条件と、負荷パターン区分格納部２１に保管されている負荷パターンの設定条件とを比較して、対象需要家の負荷パターン区分を確定する（Ｓ１６−２）。
【０２２０】
続いて、負荷パターン区分の割当処理では、確定された負荷パターン区分を、顧客名等の情報に対応させて、一時的に格納されている計算条件に追加する（Ｓ１６−３）。この一時格納された諸条件は、上述した標準負荷カーブの算出（Ｓ１７）の条件として利用される。
【０２２１】
--需要家の需要規模に応じた補正処理（Ｓ１８）--
つぎに、個別負荷カーブ推定処理Ｓ５の需要家の需要規模に応じた補正処理（図６のＳ１８）について説明をする。
【０２２２】
図３６に、需要家の需要規模に応じた補正処理（Ｓ１８）の処理の流れを示す。
【０２２３】
需要家の需要規模に応じた補正処理では、まず、需要家の需要規模に応じた補正処理では、標準負荷カーブの算出（Ｓ１７）で算出された標準負荷カーブを呼び込む（Ｓ１８-１）。
【０２２４】
続いて、需要家の需要規模に応じた補正処理では、標準負荷カーブが正規化された値であることから、個々の需要家の規模に応じた補正を行い、対象となる需要家の負荷カーブを算出する（Ｓ１８−２）。
【０２２５】
例えば、式（１３）の契約電力を用いた場合には、下記式（１７）に示すように、各需要家の契約電力を標準負荷カーブに乗ずることで最終的な負荷カーブを算出する。
【０２２６】
【数１７】
2006011715.tif 000020
【０２２７】
最後に、需要家の需要規模に応じた補正処理では、補正して得られた負荷カーブが、負荷カーブ推定結果格納部２９に格納される（Ｓ１８−３）。
【０２２８】
なお、需要家の需要規模に応じた補正処理（Ｓ１８）の説明及び標準負荷カーブモデルに関する説明において、正規化方法として契約電力を用いた方法を述べたが、正規化方法としては、各需要家（サンプル需要家及び負荷カーブ作成対象需要家）の検針値（月使用電力量や年間使用電力量、年間最大月間電力量などを契約電力の替わりに使用など）を用いることも可能である。ただし、契約電力を用いた方法が、個々の需要家の負荷カーブを電力量計によって測定する方法にくらべ、より簡便にまたより低費用での負荷カーブ作成が可能となり、また、予測などのニーズに対応した負荷カーブ測定も可能となる。
【０２２９】
７  本発明が適用された負荷カーブ推定システムの効果等
以上のように、上記負荷カーブ推定システムでは、負荷パターン区分毎に標準負荷カーブを生成しておき、推定対象の需要家が属する負荷パターン区分を特定し、特定した負荷パターン区分の標準負荷カーブを用いてその需要家の負荷カーブを算出する。
【０２３０】
このことにより上記負荷カーブ推定システムでは、インターバルメータやリアルタイムメータ等の高価な測定器を各需要家に設置することなく、需要家毎の電力消費量を、正確且つ簡易に推定することができる。
【０２３１】
なお、本発明の実施の形態として、電力の需要量を測定するシステムについて説明をしたが、本発明は電力に限らず、ガス、水道等の資源の消費量を算出するシステムにも適用することができる。
【０２３２】
つまり、電力に限らず、ガス、水道等の家庭、工場、商店等の需要家が消費するものであり、且つ、供給源（発電所又は変電所、ガスタンク、貯水設備）から需要家までを物理的に接続したネットワーク（送電線、ガス管、水道管）を介して供給されるものであれば、需要量が推定可能であるため、以上説明をした負荷カーブ推定システムを適用することができる。
【０２３３】
なお、本発明は、ハードウェア資源により構成するのみならず、コンピュータにより実行が可能なソフトウェアプログラムにより構成することも可能である。
【図面の簡単な説明】
【０２３４】
【図１】電力系統の構成を示した図である。
【図２】負荷カーブの一例を示した図である。
【図３】本発明が適用された負荷カーブ算出方法の処理手順を示す図である。
【図４】負荷カーブ推定システムのシステム構成図である。
【図５】負荷カーブ推定準備処理を説明するための図である。
【図６】個別負荷カーブ推定処理を説明するための図である。
【図７】負荷パターン区分設計処理のフローチャートである。
【図８】需要家範囲の指定の一例を示す図である。
【図９】区分設計指標の一例を示す図である。
【図１０】負荷パターン区分の一例を示す図である。
【図１１】凝集アルゴリズムを選択した際の区分設計のための処理手順を示すフローチャートである。
【図１２】数値変数をカテゴリー変数化した例を示す図である。
【図１３】クラスター数に対する融合過程を示す図である。
【図１４】クラスター数と評価値との関係を示す図である。
【図１５】分割アルゴリズムを選択した際の区分設計のための処理手順を示すフローチャートである。
【図１６】グループ数に対する分割過程を示す図である。
【図１７】グループ数と評価値との関係を示す図である。
【図１８】複数の区分設計指標を用いた場合のグループ数に対する分割過程を示す図である。
【図１９】需要サンプリング調査設計処理の手順を示すフローチャートである。
【図２０】サンプル抽出枠の設定例を示す図である。
【図２１】サンプル需要家を特定するフラグが付けられた状態の顧客情報を示す図である。
【図２２】標準負荷カーブモデルのうちの共通形状モデルのアルゴリズムを説明するための図である。
【図２３】標準負荷カーブモデルのうちの即時サンプルモデルのアルゴリズムを説明するための図である。
【図２４】標準負荷カーブモデルのうちの運転パターンモデルにモデル化できる需要の動作パターンを示した図である。
【図２５】標準負荷カーブモデルのうちの類似日モデルのアルゴリズムを説明するための図である。
【図２６】インデックスと対応付けられた標準負荷カーブのデータベース例を示す図である。
【図２７】データベースから標準負荷カーブを抽出する処理を示した図である。
【図２８】天気の各条件を数値化した例を示した図である。
【図２９】共通形状モデルを用いた場合の標準負荷カーブの算出処理の流れを示したフローチャートである。
【図３０】即時サンプルモデルを用いた場合の標準負荷カーブの算出処理の流れを示したフローチャートである。
【図３１】運転パターンモデルを用いた場合の標準負荷カーブの算出処理の流れを示したフローチャートである。
【図３２】類似日モデルを用いた場合の標準負荷カーブの算出処理の流れを示したフローチャートである。
【図３３】需要モデリングモデルを用いた場合の標準負荷カーブの算出処理の流れを示したフローチャートである。
【図３４】計算対象設定処理の流れを示したフローチャートである。
【図３５】負荷パターン区分の割当処理の流れを示したフローチャートである。
【図３６】需要家の需要規模に応じた補正処理の流れを示したフローチャートである。
【符号の説明】
【０２３５】
１０  負荷カーブ推定システム、１２  サーバ、２０  標準負荷カーブ一時格納部、２１  負荷パターン区分格納部、２２  サンプル設計結果格納部、２３  標準負荷カーブモデル格納部、２４  需要調査データ格納部、２５  気象等情報格納部、２６  顧客情報格納部、２７  系統需要格納部、２８  検針記録格納部、２９  負荷カーブ推定結果格納部
【図１】
2006011715.tif 000021
【図２】
2006011715.tif 000022
【図３】
2006011715.tif 000023
【図４】
2006011715.tif 000024
【図５】
2006011715.tif 000025
【図６】
2006011715.tif 000026
【図７】
2006011715.tif 000027
【図８】
2006011715.tif 000028
【図９】
2006011715.tif 000029
【図１０】
2006011715.tif 000030
【図１１】
2006011715.tif 000031
【図１２】
2006011715.tif 000032
【図１３】
2006011715.tif 000033
【図１４】
2006011715.tif 000034
【図１５】
2006011715.tif 000035
【図１６】
2006011715.tif 000036
【図１７】
2006011715.tif 000037
【図１８】
2006011715.tif 000038
【図１９】
2006011715.tif 000039
【図２０】
2006011715.tif 000040
【図２１】
2006011715.tif 000041
【図２２】
2006011715.tif 000042
【図２３】
2006011715.tif 000043
【図２４】
2006011715.tif 000044
【図２５】
2006011715.tif 000045
【図２６】
2006011715.tif 000046
【図２７】
2006011715.tif 000047
【図２８】
2006011715.tif 000048
【図２９】
2006011715.tif 000049
【図３０】
2006011715.tif 000050
【図３１】
2006011715.tif 000051
【図３２】
2006011715.tif 000052
【図３３】
2006011715.tif 000053
【図３４】
2006011715.tif 000054
【図３５】
2006011715.tif 000055
【図３６】
2006011715.tif 000056